Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userK2SC: Flexible systematics correction and detrending of K2 light curves using Gaussian Process regression
71
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present K2SC (K2 Systematics Correction), a Python pipeline to model instrumental systematics and astrophysical variability in light curves from the K2 mission. K2SC uses Gaussian process regression to model position-dependent systematics and time-dependent variability simultaneously, enabling the user to remove both (e.g., for transit searches) or to remove systematics while preserving variability (for variability studies). For periodic variables, K2SC automatically computes estimates of the period, amplitude and evolution timescale of the variability. We apply K2SC to publicly available K2 data from campaigns 3--5, showing that we obtain photometric precision approaching that of the original Kepler mission. We compare our results to other publicly available K2 pipelines, showing that we obtain similar or better results, on average. We use transit injection and recovery tests to evaluate the impact of K2SC on planetary transit searches in K2 PDC (Pre-search Data Conditioning) data, for planet-to-star radius ratios down Rp/Rstar = 0.01 and periods up to P = 40 d, and show that K2SC significantly improves the ability to distinguish between correct and false detections, particularly for small planets. K2SC can be run automatically on many light curves, or manually tailored for specific objects such as pulsating stars or large amplitude eclipsing binaries. It can be run on ASCII and FITS light curve files, regardless of their origin. Both the code and the processed light curves are publicly available, and we provide instructions for downloading and using them. The methodology used by K2SC will be applicable to future transit search missions such as TESS and PLATO.
rate research
Read More
The use of Gaussian processes (GPs) as models for astronomical time series datasets has recently become almost ubiquitous, given their ease of use and flexibility. GPs excel in particular at marginalization over the stellar signal in cases where the variability due to starspots rotating in and out of view is treated as a nuisance, such as in exoplanet transit modeling. However, these effective models are less useful in cases where the starspot signal is of primary interest since it is not obvious how the parameters of the GP model are related to the physical properties of interest, such as the size, contrast, and latitudinal distribution of the spots. Instead, it is common practice to explicitly model the effect of individual starspots on the light curve and attempt to infer their properties via optimization or posterior inference. Unfortunately, this process is degenerate, ill-posed, and often computationally intractable when applied to stars with more than a few spots and/or to ensembles of many light curves. In this paper, we derive a closed-form expression for the mean and covariance of a Gaussian process model that describes the light curve of a rotating, evolving stellar surface conditioned on a given distribution of starspot sizes, contrasts, and latitudes. We demonstrate that this model is correctly calibrated, allowing one to robustly infer physical parameters of interest from one or more stellar light curves, including the typical radii and the mean and variance of the latitude distribution of starspots. Our GP has far-ranging implications for understanding the variability and magnetic activity of stars from both light curves and radial velocity (RV) measurements, as well as for robustly modeling correlated noise in both transiting and RV exoplanet searches. Our implementation is efficient, user-friendly, and open source, available as the Python package starry-process.
In this note we present the starry_process code, which implements an interpretable Gaussian process (GP) for modeling variability in stellar light curves. As dark starspots rotate in and out of view, the total flux received from a distant star will change over time. Unresolved flux time series therefore encode information about the spatial structure of features on the stellar surface. The starry_process software package allows one to easily model the flux variability due to starspots, whether one is interested in understanding the properties of these spots or marginalizing over the stellar variability when it is treated as a nuisance signal. The main difference between the GP implemented here and typical GPs used to model stellar variability is the explicit dependence of our GP on physical properties of the star, such as its period, inclination, and limb darkening coefficients, and on properties of the spots, such as their radius and latitude distributions. This code is the Python implementation of the interpretable GP algorithm developed in Luger, Foreman-Mackey, and Hedges (2021).
Instrumental data are affected by systematic effects that dominate the errors and can be relevant when searching for small signals. This is the case of the K2 mission, a follow up of the Kepler mission, that, after a failure on two reaction wheels, has lost its stability properties rising strongly the systematics in the light curves and reducing its photometric precision. In this work, we have developed a general method to remove time related systematics from a set of light curves, that has been applied to K2 data. The method uses the Principal Component Analysis to retrieve the correlation between the light curves due to the systematics and to remove its effect without knowing any information other than the data itself. We have applied the method to all the K2 campaigns available at the Mikulski Archive for Space Telescopes, and we have tested the effectiveness of the procedure and its capability in preserving the astrophysical signal on a few transits and on eclipsing binaries. One product of this work is the identification of stable sources along the ecliptic plane that can be used as photometric calibrators for the upcoming Atmospheric Remote-sensing Exoplanet Large-survey mission.
Data from the Transiting Exoplanet Survey Satellite (TESS) has produced of order one million light curves at cadences of 120 s and especially 1800 s for every ~27-day observing sector during its two-year nominal mission. These data constitute a treasure trove for the study of stellar variability and exoplanets. However, to fully utilize the data in such studies a proper removal of systematic noise sources must be performed before any analysis. The TESS Data for Asteroseismology (TDA) group is tasked with providing analysis-ready data for the TESS Asteroseismic Science Consortium, which covers the full spectrum of stellar variability types, including stellar oscillations and pulsations, spanning a wide range of variability timescales and amplitudes. We present here the two current implementations for co-trending of raw photometric light curves from TESS, which cover different regimes of variability to serve the entire seismic community. We find performance in terms of commonly used noise statistics to meet expectations and to be applicable to a wide range of different intrinsic variability types. Further, we find that the correction of light curves from a full sector of data can be completed well within a few days, meaning that when running in steady-state our routines are able to process one sector before data from the next arrives. Our pipeline is open-source and all processed data will be made available on TASOC and MAST.
We present the first data release of the Kepler Smear Campaign, using collateral smear data obtained in the Kepler four-year mission to reconstruct light curves of 102 stars too bright to have been otherwise targeted. We describe the pipeline developed to extract and calibrate these light curves, and show that we attain photometric precision comparable to stars analyzed by the standard pipeline in the nominal Kepler mission. In this paper, aside from publishing the light curves of these stars, we focus on 66 red giants for which we detect solar-like oscillations, characterizing 33 of these in detail with spectroscopic chemical abundances and asteroseismic masses as benchmark stars. We also classify the whole sample, finding nearly all to be variable, with classical pulsations and binary effects. All source code, light curves, TRES spectra, and asteroseismic and stellar parameters are publicly available as a Kepler legacy sample.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
