A&A 373, 181-189 (2001)
DOI: 10.1051/0004-6361:20010566

The radio corona of AR Lacertae

C. Trigilio1 - C. S. Buemi1 - G. Umana1 - M. Rodonò2,3 - P. Leto1 - A. J. Beasley4 - I. Pagano2


1 - Istituto di Radioastronomia del C.N.R., Stazione VLBI di Noto, CP 141 Noto, Italy
2 - Osservatorio Astrofisico di Catania, Via S. Sofia 78, 95123 Catania, Italy
3 - Dipartimento di Fisica e Astronomia, Università di Catania, Via S. Sofia 78, 95123 Catania, Italy
4 - Owens Valley Radio Observatory, California Institute of Technology, Big Pine, CA 93513, USA

Received 12 December 2000 / Accepted 12 April 2001

Abstract
We present multifrequency VLA and VLBA observations at 8.4 GHz of the RS CVn system AR Lac, that were performed in autumn 1997 simultaneously with X-ray observations obtained from Rodonò et al. (1999). Our VLBA data indicate a resolved source with dimension close to the system separation, while the study of the flux density curve reveals a small amplitude outside of the eclipse variability. The derived five-frequencies spectra, combined with the size information from VLBA data, are compared with gyrosynchrotron emission from a two component structured source. A comparison with the results of the X-ray observations allow us to exclude the possibility that a thermal gyrosynchrotron is responsible for the radio emission, but it is compatible with the hypothesis of co-spatial X-ray and radio emitting sources.

Key words: stars: coronae - radio continuum: stars - stars: binaries: close - stars: individual: AR Lac - stars: binaries: general - techniques: interferometric


1 Introduction

The RS CVn stars constitute one of the major classes of stellar system radio sources. The intense and highly variable radio emission is strictly related to magnetic activity and originates from the interaction of mildly relativistic electrons and magnetic fields (gyrosynchrotron emission) on one or both components. Radio observations of this kind of system can, therefore, allow us to directly investigate non-thermal phenomena occurring in stellar coronae.

AR Lacertae (G2IV+K0IV) is a totally eclipsing RS CVn binary with a short, almost 2 day, orbital period (see Table 1 for the parameters of the system). Its orbital characteristics and strong coronal emission ( $L_{\rm X}\approx 10^{31}~{\rm erg~s}^{-1}$) make this source an ideal target for the studies of coronal structure and the geometry of RS CVn close binaries.

  
Table 1: AR Lac parameters.
\begin{table}
\par
\begin {tabular}{cll}
\hline
\hline
$R_{\rm K}$\space & 2.81...
...188$\space & Marino et~al. (\cite{Marino}) \\
\hline
\end{tabular}
\end{table}

The AR Lac coronal emission has been widely studied using most of the X-ray satellites launched in the last twenty years. Early observations, carried out by Swank et al. (1981) using the Einstein Solid State Spectrometer, showed two distinct temperature distributions centred at $7\times 10^6$ and 1.5 to $4\times 10^7$ K. On the basis of Einstein IPC observations, Walter et al. (1983) found that low lying coronae ($\sim$ $0.02\,R_\ast$) seem to be associated with both component stars, while a more extended corona ($\sim$$1\,R_\ast$) only with the K subgiant.
 

 
Table 2: Journal of observations.
  VLA VLBA+VLA
  C Band K Band Q Band X Band U Band
           
Nov. 2-3 (UT) 21:00-08:25 20:30-09:25 20:05-09:00 20:05-09:25 21:30-08:37
           
Nov. 3-4 (UT) 21:00-08:25 20:30-09:25 20:05-09:00 20:05-09:25 21:30-08:37
           
Number of scans 10 11 11 7 6
Duration of a scan [m] 19 25 27 44 44
Interval between scans [m] 80 80 80 60 60
Duty cycle (cal-sour-cal) [m] 2-15-2 3-8-3 3-5-3 1.5-3-1.5 1.5-3-1.5
Total Time on source [h] 2.42 2.54 2.4 2.75 2.5


EXOSAT observations by White et al. (1990) confirmed the existence of two different coronal regions associated with different temperature distributions. The observations showed evidence of modulations at low energy (<1 keV), which are not present at high energy (>1 keV). On the other hand, Ottmann et al. (1993) observed evidence of the primary minimum in all energy bandpasses of the ROSAT PSPC.

Despite the numerous X-ray observations, only a few have been carried out at radio wavelengths. Owen & Spangler (1977) first studied the spatial structure of the radio corona of AR Lac, on the basis of flux curve considerations, through observations performed at 4885 MHz with 5 VLA antennas. They did not find evidence of eclipses of either component. The same results were found by Doiron & Mutel (1984) on the basis of VLA observations carried out at 1.48 and 4.9 GHz, using all the 27 antennas of the array in B configuration. These results suggest that the size of the emitting region is larger than the whole system.

During a VLA 3-frequency monitoring of AR Lac, extended over an orbital period, Walter et al. (1987) detected the decay of a relatively, large flare. The event was not observed at 20 cm and the radio spectra seem to turn over between 2 and 6 cm. In the present paper we present the results of a multifrequency observing campaign on the binary system AR Lacertae (HD 210334), simultaneously performed with the Very Large Array (VLA) and the Very Long Baseline Array (VLBA)[*] radio interferometers in the autumn of 1997. The observations cover two orbital periods and took place contemporaneously with X-ray observations of the same system carried out with the SAX satellite by Rodonò et al. (1999).

2 Observations and data reduction

The observations were performed in 1997 in four sessions, starting each day at 20:00 UT and ending at 9:40 UT of the following day, from Oct. 31 to Nov. 4.

The aim of these observations was to obtain simultaneously radio spectra at five frequencies and high resolution maps of AR Lac at 4 cm and 2 cm. The spectra were obtained with the VLA in D configuration at 5 GHz (6 cm, C band), 8.4 GHz (4 cm, X band), 15 GHz (2 cm, U band), 22 GHz (1.3 cm, K band) and 43 GHz (0.7 cm, Q band), while the high resolution VLBA observations were performed at 8.4 and 15 GHz. In order to reach the best compromise between high sensitivity for the VLBA array and the best temporal coverage for the radio spectrum, we split the VLA into two subarrays, one observing in standard interferometric mode and alternating 3 frequencies, the other in phased array mode as elements of the VLBA.

2.1 VLA data

The first VLA subarray, which consisted of the 13 telescopes equipped with 43 GHz receivers, was used to observe alternately at 43, 22 and 5 GHz, with two independent 50 MHz wide contiguous bands, in Right and Left circular polarization.

As a phase calibrator we used BL Lac, while as a primary flux calibrator we used 3C 286, which was observed at the beginning of each session at all the 3 frequencies.

To minimize the effects of the atmosphere on the phase stability at higher frequencies we observed the phase calibrator more frequently at Q and K bands than at the other frequencies. Other details of the observational strategy are given in Table 2.

Calibration and data editing were performed using the standard procedures of the AIPS package.

Since we intend to study only our radio observations taken simultaneously with the X-ray observations, which started on 1997 Nov. 2, we restrict our analysis to the last two sessions, i.e. Nov. 2-3 and Nov. 3-4. On the first two days the source was found in an active period, with flux densities reaching about 40 mJy at 6 cm. This indicates an intrinsic high variability of the radio emission, which will be analyzed in a following paper.

2.2 VLBA data

The second subarray, consisting of all the remaining 14 antennas, was used to observe in phased-array mode at X and U bands as an element of the VLBA array. The presence of half the phased VLA, corresponding to a 90 m telescope, increased the sensitivity of the VLBA.

The data from the VLA phased-array can be also analyzed as a standard interferometer, therefore total fluxes were computed for all the observed sources. The VLBA observations were performed alternating between X and U bands, with a typical scan lasting 45 min. For each observing frequency, we observed in dual polarization mode, with a total bandwidth of 64 MHz. Since AR Lac has generally a low flux density (from few to some tens of mJy), we used the phase-reference technique (Beasley & Conway 1995), consisting of the rapid switching between a strong calibrator, close to the target source, and the source itself. In this way, in the successive analysis, phase calibration of the calibrator can be applied to the target source. We used BL Lac (3.6 degrees apart from AR Lac) as the reference source. The data tapes were correlated at the VLBA correlator at the Array Operation Center in Socorro (New Mexico).

For VLA data, calibration and data editing were performed using the standard procedures of the AIPS package. The flux calibration was performed by using the measured system temperature and gain curves for VLBA telescopes and from the measurements of the ratio $T_{\rm ant}/T_{\rm sys}$ made on the calibrator BL Lac for the phased-array, once the flux density of this source was determined from the VLA data. Delay and delay rate were determined for BL Lac, which was then self-calibrated. The final map of the self-calibration was used as a model for the final global fringe fitting, whose solutions were applied to our target. A preliminary map of AR Lac at 8.4 GHz was made to check its position, and it was found about 14 milliarcsec (mas) away from the phase tracking center. This was due to the fact that the coordinates used at the correlator were the heliocentric ones, and the displacement of the source was due to the annual parallax. The visibilities of AR Lac were then phase-rotated in order to bring it to the phase center. The operation makes it possible to time average the data over a relatively long time (tens of minutes) without loss of signal in the longest baselines, allowing us to analyze the amplitude of the visibility as a function of the baseline length. If the source is not at the phase center, a rapid change of phase would occur at the longest baseline.

  \begin{figure}
\par\includegraphics[width=15.4cm,clip]{ms10553f1a.eps} \includegraphics[width=15.4cm,clip]{ms10553f1b.eps}\end{figure} Figure 1: Top panel: VLA flux curves of AR Lac binned over 20 min. Dot-dashed vertical lines indicate the times of the spectra labelled by the letters on the top and shown in Fig. 2. Lower panel: the contemporaneous SAX light curves (from Rodonò et al. 1999).
Open with DEXTER

The source was successfully detected at 8.4 GHz, but not at 15 GHz, probably because of a poor phase stability at this frequency, which we could not correct with the phase referencing technique.

  \begin{figure}
\par\includegraphics[width=15.8cm,clip]{ms10553f2.eps}\end{figure} Figure 2: AR Lac VLA spectra obtained at times indicated by dashed vertical lines in Fig. 1.
Open with DEXTER

3 Results

3.1 Total flux and Spectra from VLA data

In order to investigate the temporal radio flux behaviour of the source and locate possible rapid flux variations, we analysed each scan with the AIPS program DFTPL. This procedure performs the direct Fourier transform of the visibility function in a limited portion of the UV plane. The results of such an analysis led to the conclusion that, even if the radio flux of AR Lac is variable, it does not change significantly on timescales shorter than 45 min, which is about the time length of each scan.

In Fig. 1 the flux curves, obtained by averaged flux in 20 min bins and using the new ephemeris by Marino et al. (1998), are shown.

The source was detected at all frequencies during both sessions, except during the last session in the Q-band, where only an upper limit of 0.6 mJy for the flux density can be given. This value corresponds to the rms of the cleaned map integrated over the entire 14 hours. A clear flux variation, more evident at the lowest frequencies (up to a factor of 2 in C-band), seems to indicate the presence of inhomogeneous features on one or both of the system's components. Unfortunately, AR Lac was not visible from the VLA during eclipses. We observed a clear flux maximum at phase $\approx$0.88, and a clear decay at the beginning of each of the two sessions, suggesting a flux modulation versus phase. No significant circular polarization at any observed frequency was detected.

Since each band was not observed continuously, it is not possible to build radio spectra with simultaneous flux measurements at all the frequencies. However, the variations of the flux are slow, and it was possible to make reasonable interpolations in order to obtain the radio spectra. Six spectra corresponding to different times are shown in Fig. 2. They all show negative spectral indexes $\alpha$ ( $S \propto\, \nu^{\alpha}$), implying that the source is optically thin in the available range of frequencies ($\alpha$ varies from -0.16 to -0.45). There is no evidence of a turn-over frequency, which is presumable lower than 5 GHz. The spectra show approximately the same shape, and the main difference is variation in flux density.

3.2 Spatial structure from VLBA data


  \begin{figure}
\leavevmode
\resizebox{17cm}{!}{
\includegraphics{ms10553f3a.eps}
\includegraphics{ms10553f3b.eps}} \end{figure} Figure 3: VLBA maps of AR Lac at 8.4 GHz for Nov. 2-3 (left) and Nov. 3-4 (right). The rms of both maps is 0.06 mJy/beam and the peak intensity are 1.26 and 1.07 mJy, respectively. For the first map the contours correspond to 0.25, 0.38, 0.50, 0.63, 0.76, 0.88, 1. 01, 1.13, 1.21 mJy, for the second map correspond to 0, 0.21, 0.32, 0.43, 0.54, 0.64, 0.75, 0.86, 0.96, 1.03 mJy. On the lower left corner the restoring beam is shown.
Open with DEXTER


  \begin{figure}
\par\resizebox{16.2cm}{!}{
\includegraphics{ms10553f4a.eps}
\includegraphics{ms10553f4b.eps}} \end{figure} Figure 4: Correlated flux at 3.6 cm on the VLA-VLBA baselines as function of u-v distance, for Nov. 2-3 and Nov. 3-4. The lines represent the visibility of the core-halo model (continuous lines), of the single core (dot lines) and halo (dashed lines) as discussed in Sect. 4.
Open with DEXTER

In Fig. 3 the VLBA cleaned map of AR Lac for the two sessions, which were obtained from the phase reference calibration, are shown. The FWHMs of the beam for both the maps are $1.9 \times 1.2$ mas, with a position angle of $-8\hbox{$.\!\!^\circ$ }78$ and $\sigma\approx\,0.06$ mJy. We did not measure significant circular polarisation from the VLBA data. The source appears resolved and the two maps do not show significant changes in their structure. The evidence of a resolved structure is confirmed by the analysis of the visibility amplitude versus the projected baseline spacing for the VLBA baselines involving VLA, which is shown in Fig. 4. The correlated flux has been first coherently averaged through the scan and then the amplitudes normalized to the zero spacing value, which was assumed to be equal to the total flux obtained from the VLA measurements in the same scan. Then, the visibility amplitudes have been averaged for each baseline, whose length variation during the observation is indicated by the horizontal bars. In order to measure the size of the emitting regions, or to derive the upper limit, we fitted the data with a two dimensional Gaussian model by using the AIPS routines UVFIT and JMFIT. The first routine gives an acceptable fit of UV-data with an elliptical Gaussian with FWHM of $1.4 \times 1.2\pm 0.2$ mas, while the second one fits the source image with a Gaussian component having FWHM of $1.3 \times 1.1 \pm 0.2$ mas. The two models are then consistent with each other and with a circular symmetric source. Using a circular Gaussian to fit the data, we obtain a source size of $1.2 \pm 0.2$ mas, corresponding to $4\,R_{\rm K}$ (see Table 1) or 1.2 times the system separation ($\simeq$1 mas). Our data are also compatible with two component models, where the first component is resolved, with a dimension between 1.3 and 1.6 mas, while the second component is a point source. On the basis of $\chi^2$ results we cannot exclude the presence of such a compact core.
 

 
Table 3: Parameters of the core-halo.
  CORE HALO
B [Gauss] 200 30
$N_{\rm rel}$ [cm-3] $8\times 10^5$ $3\times 10^4$
Spectrum d d
  [mas] [mas]
a 0.16 1.4
b 0.14 1.1
c 0.15 1.3
d 0.14 1.2
e 0.16 1.3
f 0.13 1.1
$<R>/R_\odot$ 0.7 5.4
EM [cm-3]$^{\dag }$ $2.3\times10^{53}$ $3.1\times10^{53}$
T [K]$^{\dag }$ $(0.93-1.02)\times10^{7}$ $(2.47-2.73)\times10^{7}$
$N_{\rm e}$ [cm-3] $2.3\times10^{10}$ $1.2\times10^{9}$

$^{\dag }$
Rodonò et al. (1999).

   
4 Discussion

The availability of both microwave spectra and spatial information at the same time gives us the opportunity to model the coronae of AR Lac. As shown in previous papers (Umana et al. 1993; Umana et al. 1999), the observed flat radio spectra of Algols and RS CVn type binary systems cannot be reproduced by an homogeneous source model. VLBI observations of the close binary systems Algol and UX Arietis pointed out the existence of a two component structure in the coronal layers (Mutel et al. 1985; Lestrade et al. 1998): a compact core, coinciding with the active star, and a larger halo, having approximatively the size of the entire system.

We used the core-halo model developed by Umana et al. (1993) to fit the observed spectra of AR Lac, in order to check if the flux variability can be attributed to the variation of one of the parameters. Although it is not possible to derive a unique solution, we have limited the sizes of core and halo on the basis of the VLBA observation results. In particular, for the halo we assumed the size measured from the VLBA data, and for the core a size smaller than the resolution limit of the interferometer.

  \begin{figure}
\par\includegraphics[width=15.8cm,clip]{ms10553f5.eps} \end{figure} Figure 5: Comparison between the observed radio spectra of AR Lac shown in Fig. 2 and the computed spectra obtained by assuming a core-halo structure for the radio source (thick line). The contribution of the halo (dot-dashed line) and core (dashed line) to the composite spectrum are also shown.
Open with DEXTER

For the magnetic field strength we assumed $B\sim B_{\rm phot}(\frac{R}{R_\ast})^{-3}$, where $B_{\rm phot}\sim 600-1000$ Gauss, as derived for other RS CVn systems (Gondoin et al. 1985; Donati et al. 1992).

According to these constraints, we derived the best core-halo model fit labeled a in Fig. 2. We then used the derived values of the average magnetic field strength Band energetic electron number density ( $N_{\rm rel}$), in fitting the other 5 spectra, under the hypothesis that these physical properties of the coronal emitting regions are stationary. Then we tried to fit the other spectra by varying only the structure size. The results of our analysis are summarised in Table 3 and are shown in Fig. 5. For the spectra from a to d, we can get a good agreement between the observed data and a core-halo structure by assuming that the flux variations are due to structure size changes in the range of 1.1 and 1.4 times the stellar radius, corresponding a variation of 0.3 mas, which is below the errors of our measurements.

The low s/n ratio of the VLBA data does not permit us to investigate the behaviour of the emitting region size for each single scan, and so as a function of the orbital phase. The source size is very close to the beam width, so it was not possible to obtain an estimate of the source size for the time ranges during which the flux density was lower then $\approx$2 mJy. Nevertheless, the results obtained for the scans with sufficient s/n ratio, that are shown in Fig. 6 by diamond symbols, suggest that the source size remains almost constant, as the measured changes fall within errors.

  \begin{figure}
\resizebox{12cm}{!}{\includegraphics{ms10553f6.eps}} %
\end{figure} Figure 6: Lower panel: source size estimates as a function of orbital phase. The diamonds represent the value obtained from UVFIT, the triangles represent the value of source size we used to fit the observed spectra with the core-halo model. Upper panel: the corresponding total flux density at 8.4 GHz.
Open with DEXTER

On the other hand, it should be noted that the slow modulation of the modelled halo size is in good agreement with the results obtained from an independent analysis of the VLBA data and plotted in Fig. 6. To fit the spectra e and h, we had to assume an "ageing'' of the relativistic electron population and a variation of its number density from $8.0\times 10^{5}$ to $7.0\times 10^{5}$ cm-3, that is needed to explain the faster decay at the higher frequency.

If the core-halo model is able to account for the radio corona of AR Lac, the visibilities of the VLBA data should fit by a two Gaussian model corresponding to the core and halo. For the Nov. 2-3 data, the average size of the core from the analysis of the spectra a, b, c, d is 0.15 mas, and the average flux density at 8.4 GHz is 1.21 mJy; for the halo, 1.25 mas and 1.74 mJy. For the Nov. 3-4 data, the same parameters from the analysis of the spectra e and f are 0.145 mas, 1.05 mJy for the core and 1.20 mas, 1.17 mJy for the halo. We then model the normalized visibility function at 8.4 GHz with that corresponding to the core-halo model (Fig. 4). It is evident that the VLBA data are consistent with the core-halo scenario derived from the analysis of the radio spectra.

VLBA data indicate a source size close to the separation of the binary components, suggesting the possibility of an emitting region located between the system components. UV emission from plasma close to the Lagrangian point in between the system components has been suggested for AR Lac (Pagano et al. 2001) and other RS CVn-type systems (Busà et al. 1999).

In partial overlap to our observations, X-ray observations of AR Lac were performed with the Beppo SAX satellite (Rodonò et al. 1999). This gave us the opportunity to determine whether the physical parameters of the radio emitting regions, derived from the comparisons between the observations and the core-halo models, are consistent with a co-spatial model for both the X-ray and radio emitting source.

Spectral analyses performed by several authors (Swank et al. 1981; Singh et al. 1995) showed that the X-ray emission from close binary systems requires at least two plasma components characterised by different temperature and volumetric emission measures (EM) to be modelled. On the basis of the first observation run, that started on Nov. 2 at 06:07 and ended on Nov. 4 at 17:50, Rodonò et al. (1999) derived for the first component $T_1=(0.93-1.02)\times 10^7$ K and $EM_1=2.3\times 10^{53}$ cm-3 and for the second $T_2=(2.47-2.73)\times 10^7$ K and $EM_2=3.1\times 10^{53}$ cm-3.

Assuming that the higher temperature component is associated with the halo and the the lower temperature component with the core, we can check if the magnetic field, as derived from the radio data, is strong enough to contain the X-ray source. This means that $\beta$, i.e. the ratio between the density of kinetic energy ( $P_{\rm c}\approx 2N_{\rm e}kT$) and the density of magnetic energy ( $P_{\rm M} \approx \frac{B^2}{8 \pi}$) has to be less then unity. If the plasma density $N_{\rm e}$ is constant over the emitting volume V, $EM=V\times N_{\rm e}^2$, and assuming the size (radius) from the analysis of the radio data of $4.7\times 10^{10}$ and $3.8\times 10^{11}$ cm for the core and the halo respectively, we get $N_{\rm e}\approx 2.3\times 10^{10}$ and $\approx$ $1.2\times 10^{9}\ {\rm cm}^{-3}$(see Table 3).

We obtain $\beta = 0.04$ for the core and $\beta =0.23$ for the halo. The physical parameters obtained from our analysis are therefore consistent with the hypothesis of a co-spatial X-ray and radio source.

We will furthermore test the possibility that the radio emission can be attributed to the same thermal electron population responsible for the observed X-ray emission. The brightness temperature of the resolved radio source at 3.6 cm, obtained from the relation

\begin{eqnarray*}T_{\rm B}=1.97\times 10^6
\frac{F_{\rm mJy}\lambda^2_{\rm cm}}{\theta_{\rm mas}^2}
\end{eqnarray*}


is $T_{\rm B}\, \simeq \, 5.39\times 10^7$ K and $T_{\rm B}\, \simeq \, 3.70\times 10^7$ K for the first and the second session, respectively. These values are higher by a factor of two than the temperature derived from the Beppo SAX data, but the order of magnitude is the same. We simulated the spectra of a core-halo structure by adopting thermal gyrosynchrotron emission and using the expression given by Dulk (1985) for the emission and absorption coefficients. By combining the physical parameters derived from Rodonò et al. (1999) for the thermal coronal plasma with our VLBA results, we can fix the temperature, the thermal plasma density and the maximum size of the emitting regions. The only free parameter left is, thus, the magnetic field strength. The spectrum obtained for B=600 Gauss for the core and B=200 Gauss for the halo is plotted in Fig. 7 and compared with the spectra obtained from the VLA maps over the two sessions.
  \begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ms10553f7.eps}} \end{figure} Figure 7: Simulated spectra from a core-halo structure for thermal gyrosynchrotron obtained with a magnetic field of 600 Gauss in the core and 200 Gauss in the halo (thick line) and $10\,000$ Gauss in the core and 300 Gauss (dashed line) for the halo. Overlaid are the average spectra for the two sessions.
Open with DEXTER

It is evident that gyrosynchrotron emission from the same thermal population responsible for the X-ray emission is not able to account for the observed spectra, unless magnetic field strengths higher then 1000 Gauss are considered. Moreover, even assuming such an intense magnetic field, it is not possible to reproduce the quite flat observed spectra. Beasley & Guedel (2000) reached a similar result from simultaneous radio and X-ray observations of the RS CVn-type binary system UX Ari during quiescence.

5 Conclusions

In this paper we presented the results of VLA multiwavelength observations and the first 8.4 GHz VLBA images of AR Lac we obtained in November 1997. These images show a spatially resolved structure with a diameter of $1.2 \pm 0.2$ mas, as derived from model fits. Since the source was very compact, with dimension close to the beam size, and the flux density level was very low, it was not viable to derive the source diameter for each phase. Still, whenever possible, the source size was derived (within 1-2 $\sigma$). No significant variations of the source size are apparent.

The flux density curves show a slight flux modulation, more evident at lower frequencies, which suggests the presence of inhomogeneous structures. The five-frequency spectra show a slightly negative spectral index, which is characteristic of an optically thin source.

By combining the spectral information, the estimate of the source size and the results obtained from simultaneous X-ray observations (Rodonò et al. 1999), we conclude that:

The application of the core-halo model to the observed spectra allows us to reproduce the flux modulation observed between phases 0.7 and 0.9, by keeping fixed the physical parameters which characterise the two structures but changing the size within the range defined from the VLBA resolution. Further observations will be able to verify the stability of the structure that can cause the modulation. In particular, since the numerous X-ray observations of this binary system seem to indicate structural changes in its corona (see Rodonò et al. 1999 for a summary) it will be extremely important to find out if this is also true for the radio corona and to check if our hypothesis of a co-spatial emitting region is confirmed.

In addition, multiwavelength VLA observations to be carried out during eclipses would be a crucial test of the hypothesis of an active region located in between the system components.

References

 


Copyright ESO 2001