For the data analysis we use the *modulation index m*, which is
related to the mean flux density <*S*> and the rms flux density
variations
by

(1) |

It provides a measure of the strength of the variations observed, although it does not take into account the error of the individual measurement. As a criterion whether a source is variable or not, we perform a test with

(2) |

where the

To compare the modulation indices at different wavelengths and epochs, and to
obtain a quantity that is independent of the measurement noise, we define the
*variability amplitude* for the variable sources as

(3) |

where

For the polarized intensity the corresponding quantities are

(4) |

(5) |

and for the polarization angle we use

(6) |

(7) |

We note that this description of the polarization variability is somewhat problematic because the polarization measurements are more strongly affected by systematic effects than the total intensity data, and because the statistical errors have a Gaussian distribution in the Stokes parameters

Structure functions can be used for the analysis of the characteristics of the
variability, and to search for typical timescales and periodicities. The first
order structure function
is defined as (see Simonetti et al. 1985):

(8) |

with being the mean over time again. Since our data are unevenly sampled, we derived the structure functions using an interpolation algorithm. For any given time lag , the value of is calculated by linear interpolation of the two adjacent data points. Each structure function was derived twice, firstly starting from the beginning of the time series, and secondly starting at the end and proceeding backwards. This provides a rough assessment of the errors caused by the interpolation process (but not of the systematic uncertainty due to the finite length of the observations).

A typical timescale in the lightcurve (i.e., the time between a minimum and a maximum or vice versa) is indicated by a maximum of the structure function, while a period in the lightcurve causes a minimum of the structure function (e.g. Heidt & Wagner 1996). Following the nomenclature of Heeschen et al. (1987), structure functions can be used to define variability types, depending on the observed timescale. A source whose structure function reaches a maximum within the observing period is called type II (i.e., a fast-varying source). In the case of a monotonically increasing structure function (i.e., when the source shows no typical timescale less than the total observing period), it is assigned type I. Since this definition is problematic for longer observing campaigns as the ones in February 1990 and October 1992, we define a variable source whose structure function reaches a maximum or plateau within days as being of type II (cf. Wagner & Witzel 1995), otherwise type I. Non-variable sources are assigned type 0.

Copyright The European Southern Observatory (ESO)