ترغب بنشر مسار تعليمي؟ اضغط هنا

Astrophysically robust systematics removal using variational inference: application to the first month of Kepler data

59   0   0.0 ( 0 )
 نشر من قبل Suzanne Aigrain
 تاريخ النشر 2013
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Space-based transit search missions such as Kepler are collecting large numbers of stellar light curves of unprecedented photometric precision and time coverage. However, before this scientific goldmine can be exploited fully, the data must be cleaned of instrumental artefacts. We present a new method to correct common-mode systematics in large ensembles of very high precision light curves. It is based on a Bayesian linear basis model and uses shrinkage priors for robustness, variational inference for speed, and a de-noising step based on empirical mode decomposition to prevent the introduction of spurious noise into the corrected light curves. After demonstrating the performance of our method on a synthetic dataset, we apply it to the first month of Kepler data. We compare the results, which are publicly available, to the output of the Kepler pipelines pre-search data conditioning, and show that the two generally give similar results, but the light curves corrected using our approach have lower scatter, on average, on both long and short timescales. We finish by discussing some limitations of our method and outlining some avenues for further development. The trend-corrected data produced by our approach are publicly available.

قيم البحث

اقرأ أيضاً

118 - S. Aigrain 2017
We present ARC2 (Astrophysically Robust Correction 2), an open-source Python-based systematics-correction pipeline to correct for the Kepler prime mission long cadence light curves. The ARC2 pipeline identifies and corrects any isolated discontinuiti es in the light curves, then removes trends common to many light curves. These trends are modelled using the publicly available co-trending basis vectors, within an (approximate) Bayesian framework with `shrinkage priors to minimise the risk of over-fitting and the injection of any additional noise into the corrected light curves, while keeping any astrophysical signals intact. We show that the ARC2 pipelines performance matches that of the standard Kepler PDC-MAP data products using standard noise metrics, and demonstrate its ability to preserve astrophysical signals using injection tests with simulated stellar rotation and planetary transit signals. Although it is not identical, the ARC2 pipeline can thus be used as an open source alternative to PDC-MAP, whenever the ability to model the impact of the systematics removal process on other kinds of signal is important.
A survey of known threshold excitations of mirror systems suggests a means to estimate masses of nuclear systems that are uncertain or not known, as does a trend in the relative energies of isobaric ground states. Using both studies and known mirror- pair energy differences, we estimate the mass of the nucleus 17-Na and its energy relative to the p+16-Ne threshold. This model-free estimate of the latter is larger than that suggested by recent structure models.
Images taken with modern detectors require calibration via flat fielding to obtain the same flux scale across the whole image. One method for obtaining the best possible flat fielding accuracy is to derive a photometric model from dithered stellar ob servations. A large variety of effects have been taken into account in such modelling. Recently, Moehler et al. (2010) discovered systematic variations in available flat frames for the European Southern Observatorys FORS instrument that change with the orientation of the projected image on the sky. The effect on photometry is large compared to other systematic effects that have already been taken into account. In this paper, we present a correction method for this effect: a generalization of the fitting procedure of Bramich & Freudling (2012) to include a polynomial representation of rotating flat fields. We then applied the method to the specific case of FORS2 photometric observations of a series of standard star fields, and provide parametrised solutions that can be applied by the users. We found polynomial coefficients to describe the static and rotating large-scale systematic flat-field variations across the FORS2 field of view. Applying these coefficients to FORS2 data, the systematic changes in the flux scale across FORS2 images can be improved by ~1% to ~2% of the total flux. This represents a significant improvement in the era of large-scale surveys, which require homogeneous photometry at the 1% level or better.
We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signa ls present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsars sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $sim 60%$ compared to a model that includes DM variations and spin-noise only.
Missions such as WMAP or Planck measure full-sky fluctuations of the cosmic microwave background and foregrounds, among which bright compact source emissions cover a significant fraction of the sky. To accurately estimate the diffuse components, the point-source emissions need to be separated from the data, which requires a dedicated processing. We propose a new technique to estimate the flux of the brightest point sources using a morphological separation approach: point sources with known support and shape are separated from diffuse emissions that are assumed to be sparse in the spherical harmonic domain. This approach is compared on both WMAP simulations and data with the standard local chi2 minimization, modelling the background as a low-order polynomial. The proposed approach generally leads to 1) lower biases in flux recovery, 2) an improved root mean-square error of up to 35% and 3) more robustness to background fluctuations at the scale of the source. The WMAP 9-year point-source-subtracted maps are available online.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا