Free Access
Issue
A&A
Volume 580, August 2015
Article Number L2
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201425545
Published online 21 July 2015

© ESO, 2015

1. Introduction

The periodical radio emitting stellar system LS I +61°303 consists of a fast rotating Be star and a compact object of unclear nature. One hypothesis is that the compact object is a radio pulsar and that the relativistic electrons responsible for the radio outburst are produced through an interaction between the relativistic wind of the pulsar and the wind of the Be star (Maraschi & Treves 1981; Dhawan et al. 2006; Dubus 2006), i.e. similar to the system PSR B125963. This system is formed by a O9.5Ve star and a compact object with well-detected radio pulses, i.e. a radio pulsar. The radio outburst of PSR B125963 occurring around periastron passage is explained by electrons accelerated at the shock front between the two winds (Chernyakova et al. 2014, and references therein). For the compact object in LS I +61°303 with the radio outburst towards apastron, however, an alternative hypothesis suggests an accreting object, i.e. a black hole or a neutron star (e.g. Taylor et al. 1992; Romero et al. 2007; Massi & Kaufman Bernadó 2009). In this case, the radio emitting electrons belong to a jet typical of microquasars. The Bondi & Hoyle (1944) accretion theory for the eccentric orbit of LS I +61°303 with an orbital period P = 26.496 days (phase Φ) then predicts two accretion maxima and therefore two ejections of relativistic electrons (Taylor et al. 1992; Marti & Paredes 1995; Bosch-Ramon et al. 2006; Romero et al. 2007). During the first ejection near periastron (Φperi = 0.23 ± 0.02; Casares et al. 2005) relativistic electrons suffer severe inverse-Compton (IC) losses due to the proximity to the Be star. Consequently, a high energy outburst is expected, as indeed confirmed by Fermi/LAT observations (Abdo et al. 2009), but no or only negligible radio emission (Bosch-Ramon et al. 2006). For the second accretion peak near apastron it has been demonstrated by Bosch-Ramon et al. (2006) that the larger distance between the compact object and the Be star results in lower IC losses (i.e. a smaller high energy outburst) and the relativistic electrons can propagate out of the orbital plane into a radio jet. Recent Fermi/LAT observations also confirmed this second, smaller high energy outburst (Jaron & Massi 2014).

thumbnail Fig. 1

Top: radio outburst of LS I +61°303: flux density (mJy) vs. time (MJD) and orbital phase Φ as observed with the Effelsberg 100 m telescope at seven frequencies between 2.6 GHz and 32 GHz. Bottom: broad-band spectral index α (all frequencies).

Synchrotron radiation emitted from one relativistic electron population with density Nrel gyrating along a magnetic field B produces a power-law radio spectrum (Sνα) with spectral slope αthin< 0 above a critical frequency νbreak that depends on B and Nrel. Below this critical frequency self-absorption effects become important and the emission becomes optically thick with a spectral slope αthick = 2.5 (e.g. Kaiser 2006). For instance, the spectrum of the radio outburst in the pulsar system PSR B125963 has a spectral index α = −0.7, which is consistent with optically thin synchrotron emission (Fig. 3 in Connors et al. 2002). In contrast, the jets of microquasars often show flat/inverted radio spectra as demonstrated by Fender (2001). Plasma and magnetic field variations along the jet in fact may create regions with different B and Nrel values resulting in spectral components with different turn-over frequencies νbreak. The overlap of the optically thin part of one spectral component with the self-absorbed part of the adjacent one will result in a flat spectrum when observed with a spatial resolution insufficient to resolve the jet (Kaiser 2006). In the microquasar Cygnus X-1, radio emission has been measured with a flat spectrum, i.e. α = −0.06 ± 0.05 and α = 0.07 ± 0.04 (Fender et al. 2000). In Cygnus X-1 the black hole is powered by accretion of the stellar wind of its supergiant companion star. Since the companion is close to filling its Roche lobe, the wind is not symmetric but strongly focused towards the black hole (Miškovičová et al. 2011, and references therein). In contrast, as discussed above in the context of an accretion scenario, accretion in LS I +61°303 occurs in two particular orbital phases due to the eccentric orbit around the Be star. Is the accretion phase and subsequent ejection in LS I +61°303 able to develop a microquasar jet with a flat spectrum? As shown by Corbel et al. (2013), in the early phase of jet formation the low density particles in the jet produce an optically thin synchrotron power-law spectrum at GHz frequencies. Only when the density of the jet plasma increases, a transition to higher optical depth occurs resulting in flat/inverted radio spectra (α ≥ 0). In LS I +61°303 strong evidence for a flat radio spectrum does exist. Early VLA observations at four epochs revealed a flat spectrum between 1.5 and 22 GHz (Gregory et al. 1979). Furthermore, Strickman et al. (1998) carried out multi-frequency VLA observations sparsely covering one orbit and found deviations in the radio spectrum from a simple power law during the outburst. Finally, Massi & Kaufman Bernadó (2009) measured even inverted, optically thick spectra between 2.2 GHz and 8.3 GHz. The primary aim of the current work is to extend previous radio observations of LS I +61°303 and systematically study – for the first time – broad-band cm/mm (2.6–32 GHz) radio spectra and their evolution during a complete outburst. The densely sampled observations of LS I +61°303 were performed over a period of about three weeks (a total of 24 days) using the Effelsberg 100 m telescope at a total of seven frequency bands. The paper is structured as follows. In Sect. 2 we briefly present the observations and data reduction procedures. In Sect. 3 we present the results. Section 4 provides a short discussion and our conclusions.

thumbnail Fig. 2

Quasi-simultaneous broad-band radio continuum spectra obtained with the Effelsberg 100 m telescope. The spectral index α resulting from a power-law fit is indicated in each panel together with the time of observation. The measurement errors are comparable to or smaller than the symbol size. We note that up to MJD 56 017.2 the spectral index is flat with an average of −0.10 ± 0.02 (see text).

2. Observations and data reduction

The multi-frequency flux density measurements of the periodical (Jaron & Massi 2013) radio outburst of LS I +61°303 were performed approximately every 12 h over 14 consecutive days between 2012 March 27, orbital phase Φ = 0.31, and 2012 April 8, Φ = 0.78 (Zimmermann 2013). The orbital phase Φ is equal to the fractional part of (t−43 366.275) /Porbital (Porbital = 26.4960 ± 0.0028 d, Gregory 2002). Two additional measurements were obtained on 2012 April 12 (Φ = 0.94) and 19 (Φ = 0.2) bringing the total observing period to 24 days. At each observing session lasting for about 2 h, LS I +61°303 was observed together with primary calibrators consecutively at seven frequencies, namely 2.64, 4.85, 8.35, 10.45, 14.60, 23.02, and 32.00 GHz (110, 60, 36, 28, 20, 13, and 9 mm wavelengths) using the secondary focus receivers of the Effelsberg 100 m telescope. Flux density measurements were performed using the cross-scan technique, i.e. progressively slewing over the source position in azimuthal and elevation direction with the number of sub-scans matching the source brightness at a given frequency. The data reduction from raw telescope data to calibrated flux densities/spectra was done in the standard manner as described in Fuhrmann et al. (2008). The overall, final measurement uncertainties are of the order of 1% and 5% at lower and higher frequencies, respectively.

3. Results

The total flux density evolution at all frequencies is shown in Fig. 1 (top). At MJD 56 013.2 (Φ = 0.29), the flux density rises steeply from around 20 mJy to 100 mJy within 24 h (until MJD 56 014.2) and remains approximately at that level. After another 24 h (MJD 56 015.7), it rises yet again to the highest flux density level (around 200 mJy at the lowest frequencies) at MJD 56 016.7, i.e. the outburst reaches its peak within about 3.5 days after onset. We further note that the peak of the outburst at frequencies 10.45 GHz is reached at the same epoch, i.e. no significant delay is obvious within our time sampling of ~12 h. At higher frequencies the situation is less clear given the lower data sampling and larger measurement uncertainties. Subsequently, the slow decay of the outburst takes place (MJD  56 017.2) though intersected by several superimposed events. These sub-flares occur with rise/decay times of ~12 to 24 h and with amplitudes of a few 10 mJy. We note that during the last sub-flare (around MJD 56 024), the peak at higher frequencies (here: 10.45 and 14.60 GHz) appears to occur 24 h earlier. In Fig. 2, the broad-band radio spectra of all epochs are shown. In order to quantify the spectral behaviour of LS I +61°303, broad-band spectral indices α (Sνα) were obtained by performing a power-law fit to all available multi-frequency flux densities of a given epoch and the resulting fits are superimposed in Fig. 2. The spectral index evolution with time is shown in Fig. 1 (bottom: α vs. MJD) and is characterized by an overall long-term decreasing trend with (i) the flattest spectral indices occurring during the early raise/peak of the outburst and (ii) increasingly steeper values during the outburst decay with the steepest values reached towards the end (α ~ −0.5 to 0.7). A linear fit to the spectral indices reveals that on average a spectral steepening of Δα = 0.03 per day takes place over the total observing period of 24 days. Examining the spectral index behaviour in more detail we basically identify three stages (vertical dotted lines in Fig. 1, bottom). During the rise until shortly after the main peak, i.e. from MJD 56 013.2 to MJD 56 017.2, the broad-band spectrum remains almost flat (average spectral index α = −0.10 ± 0.02) and an even inverted spectrum is observed at high frequencies during some epochs (i.e. α14 GHz/23 GHz ~0.6, 0.4, and 0.6 at MJD 56 013.8, MJD 56 015.2, and MJD 56 015.7, respectively). Subsequently, a stage covering the decay phase follows with the superimposed sub-flares lasting until MJD 56 023.1 and with spectral indices dropping to values as low as 0.39 (mean value: −0.29 ± 0.03). Here, however, α appears to oscillate between 0.4 and 0.1 in a quasi-regular fashion, whereas the local peaks in spectral index roughly coincide with the peaks of the sub-flares seen in the total intensity light curves (Fig. 1, top). We further note that also during the decay stage several epochs show an inverted spectrum at high frequencies. During the last stage of the decay (MJD 56 025 to MJD 56 037), the spectral index finally drops to the steepest observed values ranging between 0.5 and 0.7, i.e. clearly optically thin, with an average of −0.5 ± 0.1 during this period.

4. Discussion and conclusions

During our Effelsberg campaign in 2012 March–April, LS I +61°303 was monitored at seven frequencies between 2.6 and 32 GHz over a total period of 24 days with a cadence of about 12 h (over the first 12 days). Our observations reveal the following results:

  • 1.

    The broad-band data show a main peak with nearly flat spectrafollowed by a sequence of small amplitude, optically thin events.This sequence of outbursts with different spectral characteristicsis also found in Cygnus X-3: Fig. 3b shows GBIobservations of Cygnus X-3 (Waltman etal. 1995), where a first outburst with almost flatspectrum(α ~ 0.05) is seen, followed by a minor optically thin flare. Figure 3c displays a first outburst exhibiting an inverted spectrum (α ~ 0.6) followed by a bursting phase of optically thin emission, as also seen for LS I +61°303 in the GBI observations of Fig. 3a. Clearly both sources, LS I +61°303 and Cygnus X-3, show the same sequence of events: a first outburst with flat/inverted spectrum is followed by one or more optically thin flares. A third microquasar with the same complex outburst behaviour is XTE J1752-223 (Brocksopp et al. 2013). In Fig. 3d we show one outburst at 5.5/9 GHz (see also Fig. 1 in Brocksopp et al. 2013), where the first peak with α ~ 0.2 is clearly distinguished from the subsequent optically thin flares. Brocksopp et al. (2013) associate the first flare with the compact XTE J1752-223 jet and the sequence of optically thin flares with discrete ejections or shocks.

    thumbnail Fig. 3

    Structured radio outbursts in LS I +61°303, Cygnus X-3, and XTE 1752-223: an outburst with flat/inverted spectrum is followed by optically thin events. a) LS I +61°303 (GBI data, black squares: 8.3 GHz, red circles: 2.2 GHz; see also Fig. 3 in Massi 2015). The black dashed line at MJD 50 888 indicates the optically thin flare of LS I +61°303 observed at 1.25 mm wavelength by Paredes et al. (2000), b)c) Cygnus X-3 (GBI data, black squares: 8.3 GHz, red circles: 2.2 GHz). Here, quasi-periodic oscillations of the spectral index are shown in the inset, d) XTE J1752-223 (ATCA data, black squares: 9 GHz, red circles: 5.5 GHz, whereas errors are smaller than the symbols size).

  • 2.

    During the sequence of LS I +61°303 sub-flares presented in Fig. 1, α appears to oscillate in a quasi-regular fashion. Such spectral variability has so far only been discussed in the microquasar GRS 1915+105 (see Fig. 2 in Fender et al. 2002). As apparent from the inset of Fig. 3c (2.2/8.3 GHz spectral index vs. time), Cygnus X-3 also shows similar spectral index oscillations during the three, minor sub-flares occurring between MJD 48 470 and MJD 48 480.

The dashed line in Fig. 3a indicates the epoch where Paredes et al. (2000) obtained a steep LS I +61°303 cm/mm spectrum up to 250 GHz (1998 March 15) during one of the optically thin sub-flares following the main flare. As a follow-up of our work, it is certainly of interest to observe the main outburst of LS I +61°303 also at higher radio frequencies. Indeed, flat/inverted spectra in microquasars have also been established at short millimeter wavelengths (Fender et al. 2000; Fender 2001) and it will be important to clarify if the short accretion phase in the eccentric orbit of LS I +61°303 is able to maintain a flat spectrum up to such high frequencies.

Acknowledgments

This research is based on observations with the 100 m telescope of the MPIfR (Max-Planck-Institut für Radioastronomie) at Effelsberg. We thank Frédéric Jaron and Jürgen Neidhofer for several interesting discussions and the anonymous referee for useful comments. The work is partly supported by the German Excellence Initiative via the Bonn Cologne Graduate School and the International Max Planck Research School for Astronomy and Astrophysics.

References

  1. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, ApJ, 701, L123 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bondi, H., & Hoyle, F. 1944, MNRAS, 104, 273 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bosch-Ramon, V., Paredes, J. M., Romero, G. E., & Ribó, M. 2006, A&A, 459, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Brocksopp, C., Corbel, S., Tzioumis, A., et al. 2013, MNRAS, 432, 931 [NASA ADS] [CrossRef] [Google Scholar]
  5. Casares, J., Ribas, I., Paredes, J. M., Martí, J., & Allen de Prieto, C. 2005, MNRAS, 360, 1105 [NASA ADS] [CrossRef] [Google Scholar]
  6. Chernyakova, M., Abdo, A. A., Neronov, A., et al. 2014, MNRAS, 439, 432 [NASA ADS] [CrossRef] [Google Scholar]
  7. Connors, T. W., Johnston, S., Manchester, R. N., & McConnell, D. 2002, MNRAS, 336, 1201 [NASA ADS] [CrossRef] [Google Scholar]
  8. Corbel, S., Aussel, H., Broderick, J. W., et al. 2013, MNRAS, 431, L107 [NASA ADS] [CrossRef] [Google Scholar]
  9. Dhawan, V., Mioduszewski, A., & Rupen, M. 2006, Proc. VI Microquasar Workshop, 52.1 [Google Scholar]
  10. Dubus, G. 2006, A&A, 456, 801 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Fender, R. P. 2001, MNRAS, 322, 31 [NASA ADS] [CrossRef] [Google Scholar]
  12. Fender, R. P., Pooley, G. G., Durouchoux, P., Tilanus, R. P. J., & Brocksopp, C. 2000, MNRAS, 312, 853 [NASA ADS] [CrossRef] [Google Scholar]
  13. Fender, R. P., Rayner, D., Trushkin, S. A., et al. 2002, MNRAS, 330, 212 [NASA ADS] [CrossRef] [Google Scholar]
  14. Fuhrmann, L., Krichbaum, T. P., Witzel, A., et al. 2008, A&A, 490, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Gregory, P. C. 2002, ApJ, 575, 427 [NASA ADS] [CrossRef] [Google Scholar]
  16. Gregory, P. C., Taylor, A. R., Crampton, D., et al. 1979, AJ, 84, 1030 [NASA ADS] [CrossRef] [Google Scholar]
  17. Jaron, F., & Massi, M. 2013, A&A, 559, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Jaron, F., & Massi, M. 2014, A&A, 572, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Kaiser, C. R. 2006, MNRAS, 367, 1083 [NASA ADS] [CrossRef] [Google Scholar]
  20. Maraschi, L., & Treves, A. 1981, MNRAS, 194, 1 [Google Scholar]
  21. Marti, J., & Paredes, J. M. 1995, A&A, 298, 151 [NASA ADS] [Google Scholar]
  22. Massi, M. 2015, PoS(EVN 2014)062, [arXiv:1502.07543] [Google Scholar]
  23. Massi, M., & Jaron, F. 2013, A&A, 554, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Massi, M., & Kaufman Bernadó, M. 2009, ApJ, 702, 1179 [NASA ADS] [CrossRef] [Google Scholar]
  25. Massi, M., & Torricelli-Ciamponi, G. 2014, A&A, 564, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Miškovičová, I., Hanke, M., Wilms, J., et al. 2011, Acta Polytechnica, 51, 85 [Google Scholar]
  27. Paredes, J. M., Martí, J., Peracaula, M., Pooley, G., & Mirabel, I. F. 2000, A&A, 357, 507 [NASA ADS] [Google Scholar]
  28. Romero, G. E., Okazaki, A. T., Orellana, M., & Owocki, S. P. 2007, A&A, 474, 15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Strickman, M. S., Tavani, M., Coe, M. J., et al. 1998, ApJ, 497, 419 [NASA ADS] [CrossRef] [Google Scholar]
  30. Taylor, A. R., Kenny, H. T., Spencer, R. E., & Tzioumis, A. 1992, ApJ, 395, 268 [NASA ADS] [CrossRef] [Google Scholar]
  31. Waltman, E. B., Ghigo, F. D., Johnston, K. J., et al. 1995, AJ, 110, 290 [NASA ADS] [CrossRef] [Google Scholar]
  32. Zimmermann, L. 2013, Variability of radio and TeV emitting X-ray binary systems The case of LS I +61°303. Bonn, Univ., Diss., URN:urn:nbn:de:hbz:5n-33175 [Google Scholar]

All Figures

thumbnail Fig. 1

Top: radio outburst of LS I +61°303: flux density (mJy) vs. time (MJD) and orbital phase Φ as observed with the Effelsberg 100 m telescope at seven frequencies between 2.6 GHz and 32 GHz. Bottom: broad-band spectral index α (all frequencies).

In the text
thumbnail Fig. 2

Quasi-simultaneous broad-band radio continuum spectra obtained with the Effelsberg 100 m telescope. The spectral index α resulting from a power-law fit is indicated in each panel together with the time of observation. The measurement errors are comparable to or smaller than the symbol size. We note that up to MJD 56 017.2 the spectral index is flat with an average of −0.10 ± 0.02 (see text).

In the text
thumbnail Fig. 3

Structured radio outbursts in LS I +61°303, Cygnus X-3, and XTE 1752-223: an outburst with flat/inverted spectrum is followed by optically thin events. a) LS I +61°303 (GBI data, black squares: 8.3 GHz, red circles: 2.2 GHz; see also Fig. 3 in Massi 2015). The black dashed line at MJD 50 888 indicates the optically thin flare of LS I +61°303 observed at 1.25 mm wavelength by Paredes et al. (2000), b)c) Cygnus X-3 (GBI data, black squares: 8.3 GHz, red circles: 2.2 GHz). Here, quasi-periodic oscillations of the spectral index are shown in the inset, d) XTE J1752-223 (ATCA data, black squares: 9 GHz, red circles: 5.5 GHz, whereas errors are smaller than the symbols size).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.