- 2.1 XMM-Newton EPIC
- 2.2 ROSAT PSPC
- 2.3 Timing analysis with a Bayesian change point detection method
- 2.4 HST observations
- 2.5 OPTIMA observations

2 Observations and data reduction

DP Leo was observed using XMM-Newton on 22 of November 2000 for a net exposure time of 19 949 s. DP Leo was detected in all three EPIC detectors (Turner et al. 2001; Strüder et al. 2001). The thin filter was used and the CCDs were read out in full window mode.

Before extracting source photons, the data were processed using the current release of the XMM-Newton Science Analysis System (version 5.1). Standard procedures of data screening (creation of an image and a background light curve) revealed time intervals with enhanced particle background. These intervals were excluded from the subsequent timing analysis using an approach described below (see Sect. 2.3). This reduces the accepted exposure time to 15 827 s with the EPIC-PN detector. The observations were performed without any interruption, i.e. full phase-coverage of the s binary was achieved with an average exposure of 150 s per 0.01 phase unit.

Figure 1:
Phase-averaged X-ray light curves of the ROSAT
and the XMM-Newton observations. Phase bins have a size of
0.01 phase units of the
s binary. |

The field of DP Leo was also observed with the ROSAT PSPC on May 30, 1992 (ROR 300169, PI: Cordova) and on May 30, 1993 (ROR 600263, PI: Petre) for a net exposure time of 8580 s and 23 169 s, correspondingly. Results of the 1992 observations were presented by Robinson & Cordova (1994), the results of the much more extended observations of 1993 are unpublished.

Although the net exposure time of the two ROSAT observations was larger than the binary period, in neither case was complete phase coverage achieved due to the close proximity of the periods of the satellite and of the binary.

There are further X-ray observations reported by Biermann et al. (1985), and Schaaf et al. (1987), respectively, with the EINSTEIN and EXOSAT satellites. We make use of the timing of the X-ray eclipses detected with Einstein, the EXOSAT data have a too low count-rate and are not used further in this paper.

2.3 Timing analysis with a Bayesian change point detection method

In order to study the abrupt changes of the X-ray count rate
particularly at eclipse ingress and egress we performed a
timing analysis of the datasets using a Bayesian change point
detection method^{} developed by Scargle (1998, 2000).
This method is well suited for a
statistical examination when the arrival times of individual
X-ray photons are registered (see Hambaryan et al. 1999).
It is superior to methods which work on binned data, since it requires
no a priori knowledge of the relevant time-scale of the structure
which will be investigated.

The method is applicable to data that are known to originate from a nearly ideal Poisson process, i.e. a class of independent, identically distributed processes, having zero lengths of dead time. The data gathered in XMM-Newton EPIC-PN and ROSAT PSPC observations allow the measurement of arrival times of individual X-ray photons with a resolution of 73.3 ms and 0.1 ms, respectively, a resolution much smaller than the ingress and egress time scale which is of the order of seconds. The EPIC-MOS data cannot be used for the study of the eclipse length, since they provides a resolution of only 2.6 s.

Scargle's (1998, 2000) method
decomposes a given set of photon counting data into
Bayesian blocks with piecewise constant count-rate according to
Poisson statistics.
Bayesian blocks are built by a Cell Coalescence algorithm
(Scargle 2000), which begins with a fine-grained segmentation.
It uses a Voronoi tessellation^{} of data points, where neighboring cells
are merged if allowed by the corresponding marginal likelihoods
(see Scargle 2000).

We repeat here the essential parts of the method, expanding upon particular
modifications of the original method as used in the present application.
Assume that during a continuous
observational interval of length *T*, consisting of *m* discrete
moments in time (spacecraft's "clock tick''), a set of
photon arrival times *D* (
*t*_{i}, *t*_{i+1},...,*t*_{i+n}) is registered.
Suppose now that we want to use these data to
compare two competing hypotheses, The first hypothesis is that the
data are generated from a
constant rate Poisson process (model *M*_{1})
and the second one from two-rate Poisson
process (model *M*_{2}). Evidently, model *M*_{1} is described by only
one parameter
(the count rate)
of the one rate Poisson process while the model *M*_{2} is described by
parameters ,
and .
The parameter is the time when the Poisson process switches from
to
during the total time *T* of observation,
which thus is divided in intervals *T*_{1} and *T*_{2}.

By taking as a background information (*I*) the proposition that one
of the models under consideration is true and by using Bayes' theorem
we can calculate the posterior probability
of each model by (the probability that *M*_{k}(*k*=1,2) is the correct model, see, e.g., Jaynes 1997)

where

The first factor on the right-hand side of Eq. (2) is the ratio of the

where

The second factor on the right-hand side of Eq. (2) is the prior odds ratio, which will often be equal to 1 (see below), representing the absence of an a priori preference for either model.

It follows that the Bayes factor is equal to the posterior odds when the
prior odds is equal to 1. When
,
the data favor *M*_{2}over *M*_{1}, and when
the data favor *M*_{1}.

If we have calculated the odds ratio *O*_{21}, in favor of model *M*_{2}over *M*_{1}, we can find the probability for model *M*_{2} by inverting
Eq. (2), giving

Applying this approach to the observational data set, Scargle's (1998, 2000) method returns an array of rates, , and a set of so called "change points'' , giving the times when an abrupt change in the rate is determined, i.e. a significant variation. This is the most probable partitioning of the observational interval into blocks during which the photon arrival rate displayed no statistically significant variations.

We determined the timing accuracy of these change points through simulations. We generated 1000 data sets (photon arrival times) with one change point each. The data in the two segments obeyed Poisson statistics. Each simulated data set had approximately the same characteristics as the observed data in terms of number of registered counts, spanned time, characteristic time scales of expected variations, and was analyzed exactly in the same way. The standard deviation of the distribution of change points was found to be -3) s, if phase-folded data are used. The uncertainty was larger, when data in original time sequence were used due to the smaller total number of photons involved. We adopted an uncertainty of 2.5 s for the observationally determined change points which were used to derive the eclipse length.

In Fig. 2 we visualize the outcome of the process for the ROSAT observations performed in 1993 and for the XMM-Newton observations (EPIC-PN data only). The top panel shows the distribution of the reciprocal of the time interval between neighbouring photons (1993 data), in the two lower panels the X-ray light curves binned in intervals of 5 s are plotted. The change points determined by our method are indicated by vertical lines.

2.3.1 XMM-Newton EPIC PN

As a first step, we applied the change point detection method to a background region free of any X-ray source. This allowed us to determine time intervals where the background showed no significant variation. These were regarded as good time intervals and further used for the timing analysis of DP Leo.

We extracted 3180 EPIC-PN photon events from the source, whose arrival times were corrected to the solar system barycenter using the "barycen'' task, as implemented in SAS version 5.1.

We used phase-folded data and data in original time sequence in order to determine different quantities. The length of the bright phase and the eclipse length were measured in phase-folded data, the times of individual eclipses (for a period update) were measured in original time sequence.

The mean bright-phase count-rate, the eclipse length, bright phase center and length of bright phase are listed in Table 1, whereas the times of mid-eclipse of the individual eclipses are listed in in Table 2. The times given there are barycentric Julian ephemeris days, i.e. they take into account the 14 leap-seconds introduced between the first and the last data point. Since leap-seconds were omitted by seemingly all authors in the past (all timings in the literature are given in HJD only), we computed for all eclipse times we found in the literature the leap-second correction and include those times in the table for consistency and future work.

Epoch | Mission/Det | CR | |||

[s^{-1}] |
[s] | ||||

1992.4 | ROSAT/PSPC | 0.35 | - | ||

1993.4 | ROSAT/PSPC | 0.50 | 0.57 | ||

2000.9 | XMM-Newton/EPIC | 0.25 | 0.57 |

We also performed a timing analysis of the ROSAT PSPC observations. In both cases we used a circular region with radius to extract source photons. In total 2705 and 8385 source counts were extracted, respectively. The radius chosen encompasses 85% of the events in the ROSAT point spread function. Photon arrival times were corrected to the solar barycenter, as implemented in the Extended Scientific Analysis Software System (Zimmerman et al. 1998).

These data were treated in the same manner as the XMM-Newton data and the corresponding results are also listed in Tables 1 and 2.

2.4 HST observations

Ultraviolet spectra of DP Leo were obtained on 1991 October 31 on three consecutive HST orbits. The RAPID mode was chosen to provide time resolution of 1.6914 s per low-resolution spectrum. In each orbit one complete eclipse was covered. Due to the close proximity of the orbital periods of the spacecraft and of the binary star, almost the same binary phase was covered during the three HST-orbits and full-phase coverage was not achieved. The data were published by Stockman et al. (1994). They modeled the observed spectra at certain phases with black-body spectra. Eclipse times were derived from the maxima of the time derivatives of the light curve. The data are re-analysed in the context of this paper, in order to derive a consistent set of epochs for the times of conjunction of the white dwarf. The eclipse part of the phase-folded light curve used for our analysis is shown in Fig. 3. It was derived from the original data by averaging the spectral range 1340-2400 Å, which was found free of emission lines.

2.5 OPTIMA observations

DP Leo was re-observed with the photon-counting camera OPTIMA (Straubmeier et al. 2001) attached to the 3.5 m telescope at Calar Alto on January 10, 2002. The observations started UT 02:54:57 and lasted 156 min, i.e. two eclipses were completely covered. The observing conditions were good with no significant seeing or transparency variations. The original data were binned in chunks of 1 s effective integration time and were performed in white light, i.e. without any filter. In the context of this paper we are concentrating on the eclipse properties only. The light curve, binned in 1 s time intervals of the first eclipse (cycle 56 307) is shown in Fig. 3, too. Since the first eclipse seemingly has the somewhat better noise properties, we are using in our analyse only those data. The second eclipse is consistent with the first, the results are essentially unchanged when data from the second eclipse are included.

Figure 3:
HST/FOS and OPTIMA light curves of DP Leo obtained in October 31, 1991,
and
January 10, 2002, respectively. The FOS data shown are mean fluxes
in the wavelength range 1300-2400 Å in units of
10^{-14} erg cm^{-2} s^{-1}, OPTIMA data were obtained in
white light (units are arbitrarily scaled count rates).
Phase bins are 2.7 s (FOS) and 1 s (OPTIMA), respectively.
The dotted line in the lower panel indicates the fit to the stream
component (see text for details). The data were phased according to the
final ephemeris of Eq. (5). |

Copyright ESO 2002