Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPrincipal Component Analysis to correct data systematics. Case study: K2 light curves
84
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
With growing data from ongoing and future supernova surveys it is possible to empirically quantify the shapes of SNIa light curves in more detail, and to quantitatively relate the shape parameters with the intrinsic properties of SNIa. Building such relationship is critical in controlling systematic errors associated with supernova cosmology. Based on a collection of well-observed SNIa samples accumulated in the past years, we construct an empirical SNIa light curve model using a statistical method called the functional principal component analysis (FPCA) for sparse and irregularly sampled functional data. Using this method, the entire light curve of an SNIa is represented by a linear combination of principal component functions, and the SNIa is represented by a few numbers called principal component scores. These scores are used to establish relations between light curve shapes and physical quantities such as intrinsic color, interstellar dust reddening, spectral line strength, and spectral classes. These relations allow for descriptions of some critical physical quantities based purely on light curve shape parameters. Our study shows that some important spectral feature information is being encoded in the broad band light curves, for instance, we find that the light curve shapes are correlated with the velocity and velocity gradient of the Si II $lambda$6355 line. This is important for supernova surveys, e.g., LSST and WFIRST. Moreover, the FPCA light curve model is used to construct the entire light curve shape, which in turn is used in a functional linear form to adjust intrinsic luminosity when fitting distance models.
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.
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 describe the results of principal component analysis (PCA) of up-the-ramp sampled IR array data from the HST WFC3 IR, JWST NIRSpec, and prototype WFIRST WFI detectors. These systems use respectively Teledyne H1R, H2RG, and H4RG-10 near-IR detector arrays with a variety of IR array controllers. The PCA shows that the Legendre polynomials approximate the principal components of these systems (i.e. they roughly diagonalize the covariance matrix). In contrast to the monomial basis that is widely used for polynomial fitting and linearization today, the Legendre polynomials are an orthonormal basis. They provide a quantifiable, compact, and (nearly) linearly uncorrelated representation of the information content of the data. By fitting a few Legendre polynomials, nearly all of the meaningful information in representative WFC3 astronomical datacubes can be condensed from 15 up-the-ramp samples down to 6 compressible Legendre coefficients per pixel. The higher order coefficients contain time domain information that is lost when one projects up-the-ramp sampled datacubes onto 2-dimensional images by fitting a straight line, even if the data are linearized before fitting the line. Going forward, we believe that this time domain information is potentially important for disentangling the various non-linearities that can affect IR array observations, i.e. inherent pixel non-linearity, persistence, burn in, brighter-fatter effect, (potentially) non-linear inter-pixel capacitance (IPC), and perhaps others.
High-Resolution Spectroscopy (HRS) has been used to study the composition and dynamics of exoplanetary atmospheres. In particular, the spectrometer CRIRES installed on the ESO-VLT has been used to record high-resolution spectra in the Near-IR of gaseous exoplanets. Here we present a new automatic pipeline to analyze CRIRES data-sets. Said pipeline is based on a novel use of Principal Component Analysis (PCA) and Cross-Correlation Function (CCF). The exoplanetary atmosphere is modeled with the $tau$-REx code using opacities at high temperature from the ExoMol project. In this work, we tested our analysis tools on the detection of CO and H$_2$O in the atmospheres of the hot-Jupiters HD209458b and HD189733b. The results of our pipeline are in agreement with previous results in the literature and other techniques.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
