Peranso: Light Curve and Period Analysis Software

Peranso - Light Curve and Period Analysis Software

Authors:

Paunzen et al

Abstract:

A time series is a sample of observations of well-defined data points obtained through repeated measurements over a certain time range. The analysis of such data samples has become increasingly important not only in natural science but also in many other fields of research. Peranso offers a complete set of powerful light curve and period analysis functions to work with large, astronomical data sets. Substantial attention has been given to ease-of-use and data accuracy, making it one of the most productive time series analysis software available. In this paper, we give an introduction to Peranso and its functionality.

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.

Exposure-based Algorithm for Removing Systematics out of the CoRoT Light Curves

Exposure-based Algorithm for Removing Systematics out of the CoRoT Light Curves

Authors:

Guterman et al

Abstract:

The CoRoT space mission was operating for almost 6 years, producing thousands of continuous photometric light curves. The temporal series of exposures are processed by the production pipeline, correcting the data for known instrumental effects. But even after these model-based corrections, some collective trends are still visible in the light curves. We propose here a simple exposure-based algorithm to remove instrumental effects. The effect of each exposure is a function of only two instrumental stellar parameters, position on the CCD and photometric aperture. The effect is not a function of the stellar flux, and therefore much more robust. As an example, we show that the ∼2% long-term variation of the early run LRc01 is nicely detrended on average. This systematics removal process is part of the CoRoT legacy data pipeline.

Using Machine Learning to ID Exoplanet Transit

A Machine Learning Technique to Identify Transit Shaped Signals

Authors:

Thompson et al

Abstract:

We describe a new metric that uses machine learning to determine if a periodic signal found in a photometric time series appears to be shaped like the signature of a transiting exoplanet. This metric uses dimensionality reduction and k-nearest neighbors to determine whether a given signal is sufficiently similar to known transits in the same data set. This metric is being used by the Kepler Robovetter to determine which signals should be part of the Q1-Q17 DR24 catalog of planetary candidates. The Kepler Mission reports roughly 20,000 potential transiting signals with each run of its pipeline, yet only a few thousand appear sufficiently transit shaped to be part of the catalog. The other signals tend to be variable stars and instrumental noise. With this metric we are able to remove more than 90% of the non-transiting signals while retaining more than 99% of the known planet candidates. When tested with injected transits, less than 1% are lost. This metric will enable the Kepler mission and future missions looking for transiting planets to rapidly and consistently find the best planetary candidates for follow-up and cataloging.

The Transit Detection Efficiency of Kepler With one Year's Worth of Data

Measuring Transit Signal Recovery in the Kepler Pipeline II: Detection Efficiency as Calculated in One Year of Data

Authors:

Christiansen et al

Abstract:

The Kepler planet sample can only be used to reconstruct the underlying planet occurrence rate if the detection efficiency of the Kepler pipeline is known, here we present the results of a second experiment aimed at characterising this detection efficiency. We inject simulated transiting planet signals into the pixel data of ~10,000 targets, spanning one year of observations, and process the pixels as normal. We compare the set of detections made by the pipeline with the expectation from the set of simulated planets, and construct a sensitivity curve of signal recovery as a function of the signal-to-noise of the simulated transit signal train. The sensitivity curve does not meet the hypothetical maximum detection efficiency, however it is not as pessimistic as some of the published estimates of the detection efficiency. For the FGK stars in our sample, the sensitivity curve is well fit by a gamma function with the coefficients a = 4.35 and b = 1.05. We also find that the pipeline algorithms recover the depths and periods of the injected signals with very high fidelity, especially for periods longer than 10 days. We perform a simplified occurrence rate calculation using the measured detection efficiency compared to previous assumptions of the detection efficiency found in the literature to demonstrate the systematic error introduced into the resulting occurrence rates. The discrepancies in the calculated occurrence rates may go some way towards reconciling some of the inconsistencies found in the literature.

Using Ray-tracing to Study Refraction in Exoplanetary Atmospheres

Refraction in planetary atmospheres: improved analytical expressions and comparison with a new ray-tracing algorithm

Authors:

Bétrémieux et al

Abstract:

Atmospheric refraction affects to various degrees exoplanet transit, lunar eclipse, as well as stellar occultation observations. Exoplanet retrieval algorithms often use analytical expressions for the column abundance along a ray traversing the atmosphere as well as for the deflection of that ray, which are first-order approximations valid for low densities in a spherically symmetric homogeneous isothermal atmosphere. We derive new analytical formulae for both of these quantities, which are valid for higher densities, and use them to refine and validate a new ray-tracing algorithm which can be used for arbitrary atmospheric temperature–pressure profiles. We illustrate with simple isothermal atmospheric profiles the consequences of our model for different planets: temperate Earth-like and Jovian-like planets, as well as HD 189733b, and GJ1214b. We find that, for both hot exoplanets, our treatment of refraction does not make much of a difference to pressures as high as 10 atm, but that it is important to consider the variation of gravity with altitude for GJ1214b. However, we find that the temperate atmospheres have an apparent scaleheight significantly smaller than their actual density scaleheight at densities larger than 1 amagat, thus increasing the difficulty of detecting spectral features originating in these regions. These denser atmospheric regions form a refractive boundary layer where column abundances and ray deflection increases dramatically with decreasing impact parameter. This refractive boundary layer mimics a surface, and none of the techniques mentioned above can probe atmospheric regions denser than about 4 amagat on these temperate planets.

Bayesian Exoplanet Validation

PASTIS: Bayesian extrasolar planet validation II. Constraining exoplanet blend scenarios using spectroscopic diagnoses

Authors:

Santerne et al

Abstract:

The statistical validation of transiting exoplanets proved to be an efficient technique to secure the nature of small exoplanet signals which cannot be established by purely spectroscopic means. However, the spectroscopic diagnoses are providing us with useful constraints on the presence of blended stellar contaminants. In this paper, we present how a contaminating star affects the measurements of the various spectroscopic diagnoses as function of the parameters of the target and contaminating stars using the model implemented into the PASTIS planet-validation software. We find particular cases for which a blend might produce a large radial velocity signal but no bisector variation. It might also produce a bisector variation anti-correlated with the radial velocity one, as in the case of stellar spots. In those cases, the full width half maximum variation provides complementary constraints. These results can be used to constrain blend scenarios for transiting planet candidates or radial velocity planets. We review all the spectroscopic diagnoses reported in the literature so far, especially the ones to monitor the line asymmetry. We estimate their uncertainty and compare their sensitivity to blends. Based on that, we recommend the use of BiGauss which is the most sensitive diagnosis to monitor line-profile asymmetry. In this paper, we also investigate the sensitivity of the radial velocities to constrain blend scenarios and develop a formalism to estimate the level of dilution of a blended signal. Finally, we apply our blend model to re-analyse the spectroscopic diagnoses of HD16702, an unresolved face-on binary which exhibits bisector variations.

Using CoRoT Style Analysis on Kepler, Plato, etc Data Sets

CoRoT Data Reduction By Example

Author:

Weingrill

Abstract:

Data reduction techniques published so far for the CoRoT N2 data product were targeted primarily on the detection of extrasolar planets. Since the whole dataset has been released, specific algorithms are required to process the lightcurves from CoRoT correctly. Though only unflagged datapoints must be chosen for scientific processing, some flags might be reconsidered. The reduction of data along with improving the signal-to-noise ratio can be achieved by applying a one dimensional drizzle algorithm. Gaps can be filled by linear interpolated data without harming the frequency spectrum. Magnitudes derived from the CoRoT color channels might be used to derive additional information about the targets. Depending on the needs, various filters in the frequency domain remove either the red noise background or high frequency noise. The autocorrelation function or the least squares periodogram are appropriate methods to identify periodic signals.The methods described here are not strictly limited to CoRoT data but may also be applied on Kepler data or the upcoming Plato mission.

Kepler Host Stars Prepped in Automatically for Asteroseismic Analysis

Automated preparation of Kepler time series of planet hosts for asteroseismic analysis

Authors:

Handberg et al

Abstract:

One of the tasks of the Kepler Asteroseismic Science Operations Center (KASOC) is to provide asteroseismic analyses on Kepler Objects of Interest (KOIs). However, asteroseismic analysis of planetary host stars presents some unique complications with respect to data preprocessing, compared to pure asteroseismic targets. If not accounted for, the presence of planetary transits in the photometric time series often greatly complicates or even hinders these asteroseismic analyses. This drives the need for specialised methods of preprocessing data to make them suitable for asteroseismic analysis. In this paper we present the KASOC Filter, which is used to automatically prepare data from the Kepler/K2 mission for asteroseismic analyses of solar-like planet host stars. The methods are very effective at removing unwanted signals of both instrumental and planetary origins and produce significantly cleaner photometric time series than the original data. The methods are automated and can therefore easily be applied to a large number of stars. The application of the filter is not restricted to planetary hosts, but can be applied to any solar-like or red giant stars observed by Kepler/K2.

Analyzing the Dynamics of an Exoplanet in the HD 196885 AB System

Chaotic dynamics of the planet in HD 196885 AB

Authors:

Satyal et al

Abstract:

Depending on the planetary orbit around the host star(s), a planet could orbit either one or both stars in a binary system as S type or P type, respectively. We have analysed the dynamics of the S-type planetary system in HD 196885 AB with an emphasis on a planet with a higher orbital inclination relative to the binary plane. The mean exponential growth factor of nearby orbits (MEGNO) maps are used as an indicator to determine regions of periodicity and chaos for the various choices of the planet's semimajor axis, eccentricity and inclination with respect to the previously determined observational uncertainties. We have quantitatively mapped out the chaotic and quasi-periodic regions of the system's phase space which indicate a likely regime of the planet's inclination. In addition, we inspect the resonant angle to determine whether alternation between libration and circulation occurs as a consequence of Kozai oscillations, a probable mechanism that can drive the planetary orbit to a very large inclination. Also, we demonstrate the possible higher mass limit of the planet and improve upon the current dynamical model based on our analysis.

The Distribution of Exoplanetary Radii Around Cool Stars

THE RADIUS DISTRIBUTION OF PLANETS AROUND COOL STARS

Authors:

Morton et al

Abstract:

We calculate an empirical, non-parametric estimate of the shape of the period-marginalized radius distribution of planets with periods less than 150 days using the small yet well-characterized sample of cool (T eff less than 4000 K) dwarf stars in the Kepler catalog. In particular, we present and validate a new procedure, based on weighted kernel density estimation, to reconstruct the shape of the planet radius function down to radii smaller than the completeness limit of the survey at the longest periods. Under the assumption that the period distribution of planets does not change dramatically with planet radius, we show that the occurrence of planets around these stars continues to increase to below 1 R ⊕, and that there is no strong evidence for a turnover in the planet radius function. In fact, we demonstrate using many iterations of simulated data that a spurious turnover may be inferred from data even when the true distribution continues to rise toward smaller radii. Finally, the sharp rise in the radius distribution below ~3 R ⊕ implies that a large number of planets await discovery around cool dwarfs as the sensitivities of ground-based transit surveys increase.

Using Wavelet Transforms on Kepler & CoRoT Light Curves

Wavelets: a powerful tool for studying rotation, activity, and pulsation in Kepler and CoRoT stellar light curves

Authors:

Bravo et al

Abstract:

Aims.

The wavelet transform has been used as a powerful tool for treating several problems in astrophysics. In this work, we show that the time-frequency analysis of stellar light curves using the wavelet transform is a practical tool for identifying rotation, magnetic activity, and pulsation signatures. We present the wavelet spectral composition and multiscale variations of the time series for four classes of stars: targets dominated by magnetic activity, stars with transiting planets, those with binary transits, and pulsating stars.

Methods.

We applied the Morlet wavelet (6th order), which offers high time and frequency resolution. By applying the wavelet transform to the signal, we obtain the wavelet local and global power spectra. The first is interpreted as energy distribution of the signal in time-frequency space, and the second is obtained by time integration of the local map.

Results.

Since the wavelet transform is a useful mathematical tool for nonstationary signals, this technique applied to Kepler and CoRoT light curves allows us to clearly identify particular signatures for different phenomena. In particular, patterns were identified for the temporal evolution of the rotation period and other periodicity due to active regions affecting these light curves. In addition, a beat-pattern signature in the local wavelet map of pulsating stars over the entire time span was also detected.

The Methods for Analyzing Tidal Dissipation in a Homogeneous Spherical Body

Tidal dissipation in a homogeneous spherical body. I. Methods

Authors:

Efroimsky et al

Abstract:

A formula for the tidal dissipation rate in a spherical body is derived from first principles, to correct some mathematical inaccuracies found in the literature. The development is combined with the Darwin-Kaula formalism for tides. Our intermediate results are compared with those by Zschau (1978) and Platzman (1984). When restricted to the special case of an incompressible spherical planet spinning synchronously without libration, our final formula can be compared with the commonly used expression from Peale & Cassen (1978, Eqn. 31). The two turn out to differ. In our expression, the contributions from all Fourier modes are positive-definite, this not being the case of the formula from Ibid. (The presence of negative terms in their formula was noticed by Makarov 2013.) Examples of application of our expression for the tidal damping rate are provided in the work by Makarov and Efroimsky (2014).

Multiple Approaches for Analyzing Radial-Velocity Data From Gaia Mission

A multi-method approach to radial-velocity measurement for single-object spectra

Authors:

David et al

Abstract:

The derivation of radial velocities from large numbers of spectra that typically result from survey work, requires automation. However, except for the classical cases of slowly rotating late-type spectra, existing methods of measuring Doppler shifts require fine-tuning to avoid a loss of accuracy due to the idiosyncrasies of individual spectra. The radial velocity spectrometer (RVS) on the Gaia mission, which will start operating very soon, prompted a new attempt at creating a measurement pipeline to handle a wide variety of spectral types.

The present paper describes the theoretical background on which this software is based. However, apart from the assumption that only synthetic templates are used, we do not rely on any of the characteristics of this instrument, so our results should be relevant for most telescope-detector combinations.

We propose an approach based on the simultaneous use of several alternative measurement methods, each having its own merits and drawbacks, and conveying the spectral information in a different way, leading to different values for the measurement. A comparison or a combination of the various results either leads to a "best estimate" or indicates to the user that the observed spectrum is problematic and should be analysed manually.

We selected three methods and analysed the relationships and differences between them from a unified point of view; with each method an appropriate estimator for the individual random error is chosen. We also develop a procedure for tackling the problem of template mismatch in a systematic way. Furthermore, we propose several tests for studying and comparing the performance of the various methods as a function of the atmospheric parameters of the observed objects. Finally, we describe a procedure for obtaining a knowledge-based combination of the various Doppler-shift measurements.

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.