Exo-Transmit: An Open-Source Code for Calculating Transmission Spectra for Exoplanet Atmospheres of Varied Composition

Exo-Transmit: An Open-Source Code for Calculating Transmission Spectra for Exoplanet Atmospheres of Varied Composition

Authors:

Kempton et al

Abstract:

We present Exo-Transmit, a software package to calculate exoplanet transmission spectra for planets of varied composition. The code is designed to generate spectra of planets with a wide range of atmospheric composition, temperature, surface gravity, and size, and is therefore applicable to exoplanets ranging in mass and size from hot Jupiters down to rocky super-Earths. Spectra can be generated with or without clouds or hazes with options to (1) include an optically thick cloud deck at a user-specified atmospheric pressure or (2) to augment the nominal Rayleigh scattering by a user-specified factor. The Exo-Transmit code is written in C and is extremely easy to use. Typically the user will only need to edit parameters in a single user input file in order to run the code for a planet of their choosing. Exo-Transmit is available publicly on Github with open-source licensing at this https URL .

VULCAN: an Open-Source, Validated Chemical Kinetics Python Code for Exoplanetary Atmospheres

VULCAN: an Open-Source, Validated Chemical Kinetics Python Code for Exoplanetary Atmospheres

Authors:

Tsai et al

Abstract:

We present an open-source and validated chemical kinetics code for studying hot exoplanetary atmospheres, which we name VULCAN. It is constructed for gaseous chemistry from 500 to 2500 K using a reduced C- H-O chemical network with about 300 reactions. It uses eddy diffusion to mimic atmospheric dynamics and excludes photochemistry. We have provided a full description of the rate coefficients and thermodynamic data used. We validate VULCAN by reproducing chemical equilibrium and by comparing its output versus the disequilibrium-chemistry calculations of Moses et al. and Rimmer & Helling. It reproduces the models of HD 189733b and HD 209458b by Moses et al., which employ a network with nearly 1600 reactions. Further validation of VULCAN is made by examining the theoretical trends produced when the temperature-pressure profile and carbon-to-oxygen ratio are varied. Assisted by a sensitivity test designed to identify the key reactions responsible for producing a specific molecule, we revisit the quenching approximation and find that it is accurate for methane but breaks down for acetylene, because the disequilibrium abundance of acetylene is not directly determined by transport-induced quenching, but is rather indirectly controlled by the disequilibrium abundance of methane. Therefore, we suggest that the quenching approximation should be used with caution and must always be checked against a chemical kinetics calculation. A one-dimensional model atmosphere with 100 layers, computed using VULCAN, typically takes several minutes to complete. VULCAN is written in Python and is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org).

Modeling the Atmospheres of hot Jupiters with CLOUDY

Investigation of the environment around close-in transiting exoplanets using CLOUDY

Authors:

Turner et al

Abstract:

It has been suggested that hot stellar wind gas in a bow shock around an exoplanet is sufficiently opaque to absorb stellar photons and give rise to an observable transit depth at optical and UV wavelengths. In the first part of this paper, we use the CLOUDY plasma simulation code to model the absorption from X-ray to radio wavelengths by 1-D slabs of gas in coronal equilibrium with varying densities (104−108cm−3) and temperatures (2000−106 K) illuminated by a solar spectrum. For slabs at coronal temperatures (106 K) and densities even orders of magnitude larger than expected for the compressed stellar wind (104−105cm−3), we find optical depths orders of magnitude too small (greater than 3×10−7) to explain the ∼3% UV transit depths seen with Hubble. Using this result and our modeling of slabs with lower temperatures (2000−104K), the conclusion is that the UV transits of WASP-12b and HD 189733b are likely due to atoms originating in the planet, as the stellar wind is too highly ionized. A corollary of this result is that transport of neutral atoms from the denser planetary atmosphere outward must be a primary consideration when constructing physical models. In the second part of this paper, additional calculations using CLOUDY are carried out to model a slab of planetary gas in radiative and thermal equilibrium with the stellar radiation field. Promising sources of opacity from the X-ray to radio wavelengths are discussed, some of which are not yet observed.

Efficient Geometric Probabilities of Multi-Transiting Exoplanetary Systems from CORBITS

Efficient Geometric Probabilities of Multi-Transiting Exoplanetary Systems from CORBITS

Authors:

Brakensiek et al

Abstract:

NASA's Kepler Space Telescope has successfully discovered thousands of exoplanet candidates using the transit method, including hundreds of stars with multiple transiting planets. In order to estimate the frequency of these valuable systems, it is essential to account for the unique geometric probabilities of detecting multiple transiting extrasolar planets around the same parent star. In order to improve on previous studies that used numerical methods, we have constructed an efficient, semi-analytical algorithm called CORBITS which, given a collection of conjectured exoplanets orbiting a star, computes the probability that any particular group of exoplanets can be observed to transit. The algorithm applies theorems of elementary differential geometry to compute the areas bounded by circular curves on the surface of a sphere (see Ragozzine & Holman 2010). The implemented algorithm is more accurate and orders of magnitude faster than previous algorithms, based on comparisons with Monte Carlo simulations. We use CORBITS to show that the present solar system would only show a maximum of 3 transiting planets, but that this varies over time due to dynamical evolution. We also used CORBITS to geometrically debias the period ratio and mutual Hill sphere distributions of Kepler's multi-transiting planet candidates, which results in shifting these distributions towards slightly larger values. In an Appendix, we present additional semi-analytical methods for determining the frequency of exoplanet mutual events, i.e., the geometric probability that two planets will transit each other (Planet-Planet Occultation, relevant to transiting circumbinary planets) and the probability that this transit occurs simultaneously as they transit their star. The CORBITS algorithms and several worked examples are publicly available at this https URL.

SYSTEMATICS-INSENSITIVE PERIODIC SIGNAL SEARCH WITH K2

SYSTEMATICS-INSENSITIVE PERIODIC SIGNAL SEARCH WITH K2

Authors:

Angus et al

Abstract:

From pulsating stars to transiting exoplanets, the search for periodic signals in K2 data, Kepler's two-wheeled extension, is relevant to a long list of scientific goals. Systematics affecting K2 light curves due to the decreased spacecraft pointing precision inhibit the easy extraction of periodic signals from the data. We here develop a method for producing periodograms of K2 light curves that are insensitive to pointing-induced systematics; the Systematics-insensitive Periodogram (SIP). Traditional sine-fitting periodograms use a generative model to find the frequency of a sinusoid that best describes the data. We extend this principle by including systematic trends, based on a set of "eigen light curves," following Foreman-Mackey et al., in our generative model as well as a sum of sine and cosine functions over a grid of frequencies. Using this method we are able to produce periodograms with vastly reduced systematic features. The quality of the resulting periodograms are such that we can recover acoustic oscillations in giant stars and measure stellar rotation periods without the need for any detrending. The algorithm is also applicable to the detection of other periodic phenomena such as variable stars, eclipsing binaries and short-period exoplanet candidates. The SIP code is available at https://github.com/RuthAngus/SIPK2.

Wayne - A Simulator for HST WFC3 IR Grism Spectroscopy

Wayne - A Simulator for HST WFC3 IR Grism Spectroscopy

Authors:

Varley et al

Abstract:

Wayne is an algorithm that simulates Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) grism spectroscopic frames including sources of noise and systematics. It can simulate both staring and spatial scan modes, and observations such as the transit and eclipse of exoplanets. Unlike many other instrument simulators, the focus of Wayne is on creating frames with realistic systematics in order to test the effectiveness of different data analysis methods in a variety of different scenarios. This approach is critical for method validation and optimising observing strategies. In this paper we describe Wayne's implementation for WFC3 in the near-infrared channel with the G141 and G102 grisms. We compare the simulations to real data, obtained for the exoplanet HD 209458b to verify the accuracy of the simulation. The simulated data described in this paper is available now at www.ucl.ac.uk/exoplanets/wayne/. We plan to release this tool to the community as open source software in the near future.

GRAPHIC: Geneva Reduction and Analysis Pipeline for High-contrast Imaging of Planetary Companions

The Geneva Reduction and Analysis Pipeline for High-contrast Imaging of planetary Companions

Authors:

Hagelberg et al

Abstract:

We present graphic, a new angular differential imaging reduction pipeline where all geometric image operations are based on Fourier transforms. To achieve this goal the entire pipeline is parallelized making it possible to reduce large amounts of observation data without the need to bin the data. The specific rotation and shift algorithms based on Fourier transforms are described and performance comparison with conventional interpolation algorithm is given. Tests using fake companions injected in real science frames demonstrate the significant gain obtained by using geometric operations based on Fourier transforms compared to conventional interpolation. This also translates in a better point spread function and speckle subtraction with respect to conventional reduction pipelines, achieving detection limits comparable to current best performing pipelines. Flux conservation of the companions is also demonstrated. This pipeline is currently able to reduce science data produced by Very Large Telescope (VLT)/NACO, Gemini/NICI, VLT/SPHERE, and Subaru/SCExAO.

ExoData: A python package to handle large exoplanet catalogue data

ExoData: A python package to handle large exoplanet catalogue data

Author:

Varley

Abstract:

Exoplanet science often involves using the system parameters of real exoplanets for tasks such as simulations, fitting routines, and target selection for proposals. Software that bridges the barrier between the catalogues and code enables users to improve the specific repeatability of results by facilitating the retrieval of exact system parameters used in an articles results along with unifying the equations and software used. As exoplanet science moves towards large data, gone are the days where researchers can recall the current population from memory. An interface able to query the population now becomes invaluable for target selection and population analysis.

ExoData is a python interface and exploratory analysis tool for the Open Exoplanet Catalogue. It allows the loading of exoplanet systems into python as objects (Planet, Star, Binary etc) from which common orbital and system equations can be calculated and measured parameters retrieved. This allows researchers to use tested code of the common equations they require (with units) and provides a large science input catalogue of planets for easy plotting and use in research. Advanced querying of targets are possible using the database and Python programming language. ExoData is also able to parse spectral types and fill in missing parameters according to programmable specifications and equations. Examples of use cases are integration of equations into data reduction pipelines, selecting planets for observing proposals and as an input catalogue to large scale simulation and analysis of planets.

ExoData is a python package freely available on GitHub. It is open source and community contributions are encouraged. The package can be easily installed using "pip install exodata", detailed setup information is provided within.

K2P2: A Photometry Pipeline for the K2 mission

K2P2 − A photometry pipeline for the K2 mission

Authors:

Lund et al

Abstract:

With the loss of a second reaction wheel, resulting in the inability to point continuously and stably at the same field of view, the NASA Kepler satellite recently entered a new mode of observation known as the K2 mission. The data from this redesigned mission present a specific challenge; the targets systematically drift in position on a ~6 hour time scale, inducing a significant instrumental signal in the photometric time series --- this greatly impacts the ability to detect planetary signals and perform asteroseismic analysis. Here we detail our version of a reduction pipeline for K2 target pixel data, which automatically: defines masks for all targets in a given frame; extracts the target's flux- and position time series; corrects the time series based on the apparent movement on the CCD (either in 1D or 2D) combined with the correction of instrumental and/or planetary signals via the KASOC filter (Handberg & Lund 2014), thus rendering the time series ready for asteroseismic analysis; computes power spectra for all targets, and identifies potential contaminations between targets. From a test of our pipeline on a sample of targets from the K2 campaign 0, the recovery of data for multiple targets increases the amount of potential light curves by a factor ≥10. Our pipeline could be applied to the upcoming TESS (Ricker et al. 2014) and PLATO 2.0 (Rauer et al. 2013) missions.

A new Algorithm for Direct Imaging of Exoplanets and Circumstellar Disks

Improving Signal to Noise in the Direct Imaging of Exoplanets and Circumstellar Disks

Authors:

Wahhaj et al

Abstract:

We present a new algorithm designed to improve the signal to noise ratio (SNR) of point and extended source detections in direct imaging data. The novel part of our method is that it finds the linear combination of the science images that best match counterpart images with signal removed from suspected source regions. The algorithm, based on the Locally Optimized Combination of Images (LOCI) method, is called Matched LOCI or MLOCI. We show using data obtained with the Gemini Planet Imager (GPI) and Near-Infrared Coronagraphic Imager (NICI) that the new algorithm can improve the SNR of point source detections by 30-400% over past methods. We also find no increase in false detections rates. No prior knowledge of candidate companion locations is required to use MLOCI. While non-blind applications may yield linear combinations of science images which seem to increase the SNR of true sources by a factor greater than 2, they can also yield false detections at high rates. This is a potential pitfall when trying to confirm marginal detections or to re-detect point sources found in previous epochs. Our findings are relevant to any method where the coefficients of the linear combination are considered tunable, e.g. LOCI and Principal Component Analysis (PCA). Thus we recommend that false detection rates be analyzed when using these techniques.

Updates for the PlanetPack Software

PlanetPack software tool for exoplanets detection: coming new features

Authors:

Baluev et al

Abstract:

We briefly overview the new features of PlanetPack2, the forthcoming update of PlanetPack, which is a software tool for exoplanets detection and characterization from Doppler radial velocity data. Among other things, this major update brings parallelized computing, new advanced models of the Doppler noise, handling of the so-called Keplerian periodogram, and routines for transits fitting and transit timing variation analysis.

Tau Retrieval for Exoplanets: A Framework for Detecting Exoplanetary Atmospheres

Tau-REx I: A next generation retrieval code for exoplanetary atmospheres

Authors:

Waldmann et al

Abstract:

Spectroscopy of exoplanetary atmospheres has become a well established method for the characterisation of extrasolar planets. We here present a novel inverse retrieval code for exoplanetary atmospheres. Tau-REx (Tau Retrieval for Exoplanets) is a line-by-line radiative transfer fully Bayesian retrieval framework. It differentiates itself from existing codes in the literature through: 1) the optimised use of molecular line-lists from the ExoMol project; 2) an unbiased atmospheric composition prior selection, through custom built pattern recognition software; 3) the use of two independent algorithms to fully sample the Bayesian likelihood space: nested sampling as well as a more classical Markov Chain Monte Carlo approach; 4) iterative Bayesian parameter and model selection using the full Bayesian Evidence as well as the Savage-Dickey Ratio for nested models, and 5) the ability to fully map very large parameter spaces through optimal code parallelisation and scalability to cluster computing. In this publication we outline the Tau-REx framework and demonstrate, using a theoretical hot-Jupiter transmission spectrum, the parameter retrieval and model selection. We investigate the impact of Signal-to-Noise and spectral resolution on the retrievability of individual model parameters, both in terms of error bars on the temperature and molecular mixing ratios as well as its effect on the model's global Bayesian evidence.

BinHab: Tool for the Calculation of S/P-Type Habitable Zones in Binary Systems

BinHab: A Numerical Tool for the Calculation of S/P-Type Habitable Zones in Binary Systems

Authors:

Cuntz et al

Abstract:

The aim of this contribution is to introduce the numerical tool BinHab, a publicly accessible code, available at The University of Texas at Arlington, that allows the calculation of S-type and P-type habitable zones of general binary systems.

ECCSAMPLES: Software for Finding the Proper Exoplanet Eccentricity

Bayesian priors for the eccentricity of transiting planets

Authors:

Kipping et al

Abstract:

Planets on eccentric orbits have a higher geometric probability of transiting their host star. By application of Bayes' theorem, we reverse this logic to show that the eccentricity distribution of transiting planets is positively biased. Adopting the flexible Beta distribution as the underlying prior for eccentricity, we derive the marginalized transit probability as well as the a-priori joint probability distribution of eccentricity and argument of periastron, given that a planet is known to transit. These results allow to demonstrate that most planet occurrence rate calculations using Kepler data have overestimated the prevalence of planets by ~10%. Indeed, the true occurrence of planets from transit surveys is fundamentally intractable without a prior assumption for the eccentricity distribution. Further more, we show that previously extracted eccentricity distributions using Kepler data are positively biased. In cases where one wishes to impose an informative eccentricity prior, we provide a recursive algorithm to apply inverse transform sampling of our joint prior probability distribution. Computer code of this algorithm, ECCSAMPLES, is provided to enable the community to sample directly from the prior.

TRAnsits and Dynamics of Exoplanetary Systems Software for Detecting, Refining Transit Detections

TRADES: a new software to derive orbital parameters from observed Transit Times and Radial Velocities. Revisiting Kepler-11 and Kepler-9

Authors:

Borsato et al

Abstract:

Aims.

With the purpose of determining the orbital parameters of exoplanetary systems from observational data, we have developed a software, named TRADES (TRAnsits and Dynamics of Exoplanetary Systems) to simultaneously fit observed radial velocities and transit times data.

Methods.

We implemented a dynamical simulator for N-body system which also fits the available data during the orbital integration and determines the best combination of the orbital parameters using grid search, χ2 minimization, genetic algorithms, particle swarm optimization, and bootstrap analysis.

Results.

To validate TRADES, we tested the code on a synthetic three-body system and on two real systems discovered by the Kepler mission: Kepler-9 and Kepler-11. These systems are good benchmarks to test multiple exoplanet systems showing transit time variations (TTVs) due to the gravitational interaction among planets. We have found orbital parameters of Kepler-11 planets in good agreement with the values proposed in the discovery paper and with a a recent work from the same authors. We analyzed the first three quarters of Kepler-9 system and found parameters in partial agreement with discovery paper. Analyzing transit times (T0s) covering 12 quarters of Kepler data we have found a new best-fit solution. This solution outputs masses that are about the 55% of the values proposed in the discovery paper; this leads to a reduced semi-amplitude of the radial velocities of about 12.80 m/s.

Automatic Detection of Exoplanet Transit Signals

Automatic Classification of Kepler Threshold Crossing Events

Authors:

McCauliff et al

Abstract:

The Kepler Science Operations Center detects interesting, exoplanet transit-like signals while searching over 211,000 distinct light curves. The mission has produced four catalogs of interesting objects with planet transit-like features known as Kepler Objects of Interest (KOI). The total number of objects with transit-like features identified in the light curves has increased to as many as approximately 18,000, just examining the first three years of data. This number of significant detections has become difficult for human beings to inspect by eye in a thorough and timely fashion. In order to accelerate the process by which new planet candidates are classified and to provide an independent assessment of planet candidates, we propose a machine learning approach to establish a preliminary list of planetary candidates ranked from most credible to least credible. The classifier must distinguish between three classes of detections: non-transiting phenomena, astrophysical false positives, and planet candidates. We use random forests, a supervised classification algorithm, that has an error rate of 1.34 percent with some qualifications.

PynPoint Analysis Software for Angular Differential Imaging Released

PynPoint Code for Exoplanet Imaging

Authors:

Amara et al

Abstract:

We announce the public release of PynPoint, a Python package that we have developed for analysing exoplanet data taken with the angular differential imaging observing technique. In particular, PynPoint is designed to model the point spread function of the central star and to subtract its flux contribution to reveal nearby faint companion planets. The current version of the package does this correction by using a principal component analysis method to build a basis set for modelling the point spread function of the observations. We demonstrate the performance of the package by reanalysing publicly available data on the exoplanet beta Pictoris b, which consists of close to 24,000 individual image frames. We show that PynPoint is able to analyse this typical data in roughly 1.5 minutes on a Mac Pro, when the number of images is reduced by co-adding in sets of 5. The main computational work parallelises well as a result of a reliance on SciPy and NumPy functions. For this calculation the peak memory load is 6Gb, which can be run comfortably on most workstations. A simpler calculation, by co-adding over 50, takes 3 seconds with a peak memory usage of 600 Mb. This can be performed easily on a laptop. In developing the package we have modularised the code so that we will be able to extend functionality in future releases, through the inclusion of more modules, without it affecting the users application programming interface. We distribute the PynPoint package through the central PyPi sever, and the documentation is available online.

Introducing Planet Analysis and Small Transit Investigation Software

PASTIS: Bayesian extrasolar planet validation. I. General framework, models, and performance

Authors:

Díaz et al

Abstract:

A large fraction of the smallest transiting planet candidates discovered by the Kepler and CoRoT space missions cannot be confirmed by a dynamical measurement of the mass using currently available observing facilities. To establish their planetary nature, the concept of planet validation has been advanced. This technique compares the probability of the planetary hypothesis against that of all reasonably conceivable alternative false-positive (FP) hypotheses. The candidate is considered as validated if the posterior probability of the planetary hypothesis is sufficiently larger than the sum of the probabilities of all FP scenarios. In this paper, we present PASTIS, the Planet Analysis and Small Transit Investigation Software, a tool designed to perform a rigorous model comparison of the hypotheses involved in the problem of planet validation, and to fully exploit the information available in the candidate light curves. PASTIS self-consistently models the transit light curves and follow-up observations. Its object-oriented structure offers a large flexibility for defining the scenarios to be compared. The performance is explored using artificial transit light curves of planets and FPs with a realistic error distribution obtained from a Kepler light curve. We find that data support for the correct hypothesis is strong only when the signal is high enough (transit signal-to-noise ratio above 50 for the planet case) and remains inconclusive otherwise. PLATO shall provide transits with high enough signal-to-noise ratio, but to establish the true nature of the vast majority of Kepler and CoRoT transit candidates additional data or strong reliance on hypotheses priors is needed.