Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.mrao.cam.ac.uk/wp-content/uploads/2014/04/DMahtani.pdf Дата изменения: Wed Apr 16 18:57:15 2014 Дата индексирования: Sun Apr 10 13:29:25 2016 Кодировка: Поисковые слова: arp 220

Semi-synthetic lightcurves: The effect of systematic noise in Spitzer Space Telescope observations of hot Jupiters
Deepak P. Mahtani and Pierre F. L. Maxted
Keele University, Staffordshire, ST5 5BG
Abstract
Secondary eclipse observations of hot Jupiters is one of the methods most often used to probe their atmospheres. It also allows the orbital eccentricity to be constrained. The depth of the eclipse gives information on the temperature and composition of the atmosphere and the time of mid-eclipse gives the eccentricity information. Spitzer has been the main instrument used to make these measurements to-date. These observations are strongly affected by systematic noise that must be corrected for to enable useful measurements of the eclipse depth and eclipse time. We have injected synthetic eclipses into real Spitzer IRAC data at 3.6 (channel 1) and 4.5 (channel 2) microns to test three methods of recovering the depth and phase of mid-eclipse. We find that the level of systematic noise on the eclipse depth is 0.013% and 0.011% at 3.6 and 4.5 microns respectively. For the time of mid eclipse we find that the error distribution is strongly non-Gaussian so a full analysis of the simulated light curve is required to determine if the orbit is significantly non-circular.

Method
To conduct this analysis we take a flat region of the full orbit phase curve of HAT-P-2b at 3.6 and 4.5 microns and inject into it an eclipse with the same predicted eclipse depth of WASP-35b. We use three methods of recovering the eclipse depth and time of mid eclipse: 1) polynomial fitting, 2) wavelet fitting, 3) a standard MCMC fitting code [5]. For the polynomial and wavelet fitting we inject the eclipse into different portions of the data (4096 data points long) 10,000 times and for each simulation, calculate the difference between the input and output values of the depth of the eclipse and time of mid eclipse. This was conducted for many different decorrelation polynomials. We then ran 10 light curves from each channel through the MCMC code using the standard polynomials no-time, quadratic-position for channel 1 and no-time and linear position for channel 2. We then compare the results of the polynomial wavelet and MCMC fitting.

Results
Figure 1 shows an example of real secondary eclipse lightcurves taken with Spitzer [4] and an example of one of the semi-synthetic lightcurves used in the analysis. The table below shows the input values which determine the shape and mid point of the eclipse. Figure 2 shows the difference between the input values of the depth and time of mid-eclipse and the polynomial, wavelet and MCMC fits of the the depth and time of mid eclipse. It can be seen that there are several values for each parameter that are out by as much as 3 sigma. The reason for this large discrepancy is that the MCMC error bars are too small. The same eclipse is being fit each time into different regions of the data set so we would expect that the results should have a Gaussian distribution. To test this we calculated the p-value and rank for the results of fitting each of the 10 lightcurves for each channel for both parameters of interest. If the results are Gaussian distributed then when one plots p vs rank then the results should be on a 1:1 relation. As can be seen from Figure 3, this is not the case for the depths. We then add in quadrature estimates of the systematic error until we get a value that (using the Kolmogorov-Smirnov test) give a result that is close to a 1:1 relation. For the eclipse depth we find that the systematic error required is about 0.013% and 0.011% at 3.6 and 4.5 microns, respectively . However, for the times of mid-eclipse this does not work because the real distribution of the error in this parameter is not Gaussian (high kurtosis, see table below), whereas the MCMC distribution is. Although adding a systematic error to the MCMC results can widen the distribution, it is not possible to match the shape. To show this figure 4 shows the distributions from the MCMC, polynomial and wavelet fits that were most discrepant. These are run 6 for the channel 1 data (3.2 sigma) and run 7 for channel 2 (2.6 sigma). Input Values R*/a Rp / R b depth tmid
*

Introduction
Observations taken with Spitzer suffer from systematic errors. Some of these systematic errors can be corrected for e.g. Intrapixel Sensitivity Variations (IPSVs), but the origin of some are unknown. Markov Chain Monte Carlo (MCMC) methods are most often used to explore the parameter space when fitting models to these data. However, the error bars that this method gives are usually too small because they do not take into account the remaining systematic error. Many methods have been used to produce more accurate error bars e.g. prayer bead residual permutation methods [1]. Here we present a new method of estimating the amount of extra systematic error in the data. Previous investigations have shown that the level of systematic error is about 0.01% [2] but it has been reported to be as much as 0.05% [3]. The goal of this analysis is to determine to what extent the systematic errors affect the measured eclipse depth and time of mid-eclipse.

ch1 0.120 0.120 0.300

ch2 0.120 0.120 0.300 0.00213 0 .0 0 0

Tmid Kurtosis Poly Wave MCMC

ch1

ch2 9 13 0.06 0.92 1.5 - 0 .2 3

0.00094 0.000

Figure 3: p vs. rank plots showing that with the extra systematic error added the results for the depth in channel 1 (left) and channel 2 (right) become close to Gaussian. Figure 1: (Top) example of real Spitzer observations of the occultation of WASP-26b [4], (bottom) example of the semi-synthetic lightcurves used in this analysis. The red data points are the channel 1 data and the green points are the channel 2 data. The red and green lines are the model fit to the data. (left) raw data, (middle) binned raw data, (right) binned detrended data.

Figure 4: Distributions of the time of mid-eclipse for the channel 1 (left) and channel 2 (right) for each of the decorrelation methods.

Conclusions
The polynomial and wavelet results are consistent with with the MCMC results and hence provide two independent verification of the MCMC method. The wavelet fitting method does not improve on the polynomial method. The MCMC error bars are too small and hence there are large discrepancies between the input and fitted depth and time of mid eclipse. A systematic error of 0.013% and 0.011% at 3.6 and 4.5 microns respectively is required to be added to the error bars on the eclipse depths to give more realistic error bar values. This cannot be applied to the time of mid eclipse as the error distribution is strongly non-Gaussian. The recent work that shows the systematic error is 0.05% has not be reproduced by this analysis- real variability? The time of mid-eclipse gives information on the eccentricity of the planets orbit so a full analysis of the simulated light curve is required to determine if the orbit is significantly non-circular.


Figure 2: comparison between the three methods used to recover the eclipse depth (left) and time of mid eclipse (right) for the channel 1 (top) and channel 2 (bottom).
References [1] Gillon et al. 2007, A&A, 471, L51 [2] Maxted et al. 2013,MNRAS, 428,2645 [3] Hansen et al 2014 eprint arXiv:1402.6699 [4] Mahtani D. P. et al., 2013, MNRAS, 432, 693 [5] Collier Cameron A. et al., 2007, MNRAS, 380, 1230