next up previous
Next: The overall time delays Up: Application to MCG-6-30-15 Previous: The RXTE data

Searching for lags and leads

We now apply the procedure outlined in Section 2 to these lightcurves. To begin with, we must estimate the structure function for these data. Figure 2a shows a pair-wise estimate of the continuum band structure function obtained following the method of PRH92. This figure also shows our analytic approximation which is given by eqns (13) and (14) with

eqnarray201

Using this covariance model, the PRH92 reconstruction was applied to the continuum light curve using N=5000 data points. A portion of the resulting reconstructed light curve is shown in Fig. 2b.

The next stage in the procedure is to convolve the reconstructed continuum band light curve with a trial transfer function and compare the result with the line band light curve in a tex2html_wrap_inline893 sense. We can then minimize the tex2html_wrap_inline893 statistic in order to constrain free parameters in the trial transfer function. We also minimize tex2html_wrap_inline893 over multiplicative and additive offsets between the continuum and line band light curves, i.e. we set

equation205

and minimize over B and K as well as the parameters describing the trial transfer function tex2html_wrap_inline1037 .

In this work, we choose two trial transfer functions. The first represents the case where some fraction tex2html_wrap_inline1039 of the line band flux is a delayed copy of the continuum band with a time delay tex2html_wrap_inline1041 :

equation211

The second represents the case where some fraction tex2html_wrap_inline1039 of the line band flux is a delayed and smeared copy of the continuum band flux, where a Gaussian kernel is used:

equation217

No extrapolations were performed during this procedure. In order to avoid extrapolating, the tex2html_wrap_inline893 statistic was calculated using a subset of data points. For the trial transfer function tex2html_wrap_inline895 , only data during times tex2html_wrap_inline1049 were used to compute tex2html_wrap_inline893 , where tex2html_wrap_inline1053 and tex2html_wrap_inline1055 are the times of the start and end of the reconstructed continuum light curve. For tex2html_wrap_inline917 , tex2html_wrap_inline893 is computed based upon data from times tex2html_wrap_inline1061 .

Figure 3 shows the tex2html_wrap_inline893 surfaces and confidence contours once this procedure has been performed. When displaying the tex2html_wrap_inline893 surfaces, we plot tex2html_wrap_inline897 in order to highlight the topography of the surface near the global minimum in the surface. It can be seen that the minimum of the tex2html_wrap_inline1069 2 surface corresponds to the two lines tex2html_wrap_inline1071 and tex2html_wrap_inline1073 , i.e. no time delayed component of the line band light curve is detected. Here we only show the results for tex2html_wrap_inline895 -- the tex2html_wrap_inline917 results are trivial (i.e. tex2html_wrap_inline893 surface is completely flat) since the preferred solution always have tex2html_wrap_inline1071 . The best fit values of the multiplicative and additive constants are B=0.78 and K=0.90.

 

figure243


Figure 3: Results for MCG-6-30-15: tex2html_wrap_inline893 surfaces and confidence contours resulting from applying trial transfer function tex2html_wrap_inline895 to the reconstructed continuum light curves and comparing with the line band light curve. Surfaces are plotted using tex2html_wrap_inline897 as the ordinate in order to display the topography of the region near the minimum. Contours are shown the following levels: tex2html_wrap_inline899 . The first three of these contours correspond to tex2html_wrap_inline901 , 90% and 95% for two interesting parameters and are shown in bold.


next up previous
Next: The overall time delays Up: Application to MCG-6-30-15 Previous: The RXTE data

Chris Reynolds
Tue Jan 11 17:27:37 MST 2000