Issue |
A&A
Volume 517, July 2010
|
|
---|---|---|
Article Number | A72 | |
Number of page(s) | 15 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201014382 | |
Published online | 10 August 2010 |
The double-peaked 2008 outburst of the accreting milli-second X-ray pulsar, IGR J00291+5934![[*]](/icons/foot_motif.png)
F. Lewis1,2,3 - D. M. Russell4 - P. G. Jonker5,6,7 - M. Linares4,8 - V. Tudose9,10,11 - P. Roche1,2,3 - J. S. Clark2 - M. A. P. Torres6 - D. Maitra4,12 - C. G. Bassa13 - D. Steeghs6,14 - A. Patruno4 - S. Migliari15 - R. Wijnands4 - G. Nelemans7 - L. J. Kewley16 - V. E. Stroud1,2,3 - M. Modjaz17,18 - J. S. Bloom17 - C. H. Blake6 - D. Starr17,19
1 - Faulkes Telescope Project, School of Physics and Astronomy, Cardiff University, 5 The Parade, Cardiff, CF24 3AA, Wales, UK
2 - Department of Physics and Astronomy, The Open University, Walton Hall, Milton Keynes, MK7 6AA, UK
3 - Division of Earth, Space and Environment, University of Glamorgan, Pontypridd, CF37 1DL, Wales, UK
4 - Astronomical Institute ``Anton Pannekoek'', University of Amsterdam, PO Box 94249, 1090 GE Amsterdam, The Netherlands
5 - SRON, Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA, Utrecht, The Netherlands
6 - Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
7 - Department of Astrophysics, IMAPP, Radboud University Nijmegen, Toernooiveld 1, 6525 ED, Nijmegen, The Netherlands
8 - MIT Kavli Institute for Astrophysics and Space Research, 70 Vassar Street, Cambridge, MA 02139, USA
9 - Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands
10 - Astronomical Institute of the Romanian Academy, Cutitul de Argint 5, RO-040557 Bucharest, Romania
11 - Research Center for Atomic Physics and Astrophysics, Atomistilor 405, RO-077125 Bucharest, Romania
12 - Department of Astronomy, University of Michigan, 500 Church Street, Ann Arbor, MI 48109, USA
13 - Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK
14 - Department of Physics, University of Warwick, Coventry, CV4 7AL, UK
15 - ESAC, Urb. Villafranca del Castillo, PO Box 50727, 28080 Madrid, Spain
16 - University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA
17 - Department of Astronomy, University of California, Berkeley, 601 Campbell Hall, Berkeley, CA 94720, USA
18 - Miller Fellow
19 - Las Cumbres Observatory Global Telescope Project, 6740 Cortona Dr., Santa Barbara, CA 93117, USA
Received 8 March 2010 / Accepted 3 May 2010
Abstract
Context. In August 2008, the accreting milli-second X-ray pulsar (AMXP), IGR J00291+5934, underwent an outburst lasting 100 days, the first since its discovery in 2004.
Aims. We present data from the 2008 double-peaked outburst of
IGR J00291+5934 from Faulkes Telescope North, the Isaac Newton
Telescope, the Keck Telescope, PAIRITEL, the Westerbork Synthesis Radio
Telescope and the Swift, XMM-Newton and RXTE X-ray missions. We study
the outburst's evolution at various wavelengths, allowing us to probe
accretion physics in this AMXP.
Methods. We study the light curve morphology, presenting the
first radio-X-ray spectral energy Distributions (SEDs) for this source
and the most detailed UV-IR SEDs for any outbursting AMXP. We show
simple models that attempt to identify the emission mechanisms
responsible for the SEDs. We analyse short-timescale optical
variability, and compare a medium resolution optical spectrum with
those from 2004.
Results. The outburst morphology is unusual for an AMXP, comprising two peaks, the second containing a ``plateau'' of 10 days
at maximum brightness within 30 days of the initial activity. This
has implications on duty cycles of short-period X-ray transients. The
X-ray spectrum can be fitted by a single, hard power-law. We detect
optical variability of
0.05 mag,
on timescales of minutes, but find no periodic modulation. In the
optical, the SEDs contain a blue component, indicative of an irradiated
disc, and a transient near-infrared (NIR) excess. This excess is
consistent with a simple model of an optically thick synchrotron jet
(as seen in other outbursting AMXPs), however we discuss other
potential origins. The optical spectrum shows a double-peaked H
profile, a diagnostic of an accretion disc, but we do not clearly see
other lines (e.g. He I, II) that were reported in 2004.
Conclusions. Optical/IR observations of AMXPs appear to be
excellent for studying the evolution of both the outer accretion disc
and the inner jet, and may eventually provide us with tight constraints
to model disc-jet coupling in accreting neutron stars.
Key words: pulsars: general - X-rays: binaries - stars: neutron - accretion, accretion disks
1 Introduction
1.1 Accreting milli-second X-ray pulsars
Low-mass X-ray Binaries (LMXBs) comprise a compact primary (neutron star or black hole) and a low-mass (usually K or M class) secondary with the compact object accreting material from its companion by means of Roche lobe overflow. It was postulated in the early 1980s that accreting milli-second X-ray pulsars would form a sub group of neutron star LMXBs (e.g. Alpar et al. 1982; Radhakrishnan & Srinivasan 1982) and thus prove to be the link between accreting LMXBs and isolated milli-second pulsars. The mechanism believed to be responsible is the spin-up of an old, weakly magnetized, neutron star caused by the accretion of matter and angular momentum from a donor star. This would result in a compact transient X-ray binary (its transient nature demonstrating a period of increased mass transfer rate onto the neutron star) with the neutron star spinning so rapidly that, once the accretion has ceased, it can once again become a rapid radio pulsar (see Wijnands 2006 and references therein). This was subsequently confirmed observationally with the discovery of pulsations from the transient source SAX J1808.4-3658 (Wijnands & van der Klis 1998) in 1998 April and the first LMXB to be detected as a milli-second radio pulsar, FIRST J102347.6+003841 (Archibald et al. 2009).
Properties of AMXPs include weaker peak luminosities in outburst than
``classical'' (i.e. not showing pulsations during outburst) neutron
star LMXBs, extremely low-mass companions and faint quiescent X-ray
luminosities (1032 erg s-1), dominated by a hard power-law component (see Wijnands et al. 2008; Wijnands 2006
and references therein). For most AMXPs, the duty cycle is a few weeks
of outburst every few years. Taking into account the recent (2009
August and September respectively) discoveries of
NGC 6440 X-2 (Altamirano et al. 2009; Heinke et al. 2009) and IGR J17511-305 (Markwardt et al. 2009; Bozzo et al. 2010), there are currently twelve known AMXPs (Wijnands et al. 2008). Observationally, these systems have orbital periods of
min-4 h.
Little is known about the optical counterparts of these systems; many
have not been detected in the optical/infrared regime, even in outburst
(e.g. XTE J1751-305; Jonker et al. 2003).
Some LMXBs suffer from a high level of absorption from interstellar
dust and/or lie at large distances, and those systems with detections,
particularly AMXPs, are often very faint (due to their relatively
low-mass donors), especially when in quiescence (see D'Avanzo et al. 2009).
1.2 The accreting milli-second X-ray pulsar IGR J00291+5934
The accreting millisecond X-ray pulsar IGR J00291+5934 (hereafter
``00291+59'') was first discovered on 2004 December 2 (MJD 53341) in
outburst during a routine Galactic Plane Scan of the Cassiopeia region
by INTEGRAL (Eckert et al. 2004). It was given the co-ordinates RA 00:29.1, Dec +59:34 (J2000.0) with an error radius of 1
5 (l = 120
1,
). In three successive pointings, the IBIS/ISGRI instrument recorded an average flux of
mCrab (20-60 keV). During its only pointing, JEM-X detected the source with a flux of
mCrab (3-10 keV) (Eckert et al. 2004).
On the following day, using the RXTE PCA (Rossi X-ray Timing Explorer
Proportional Counter Array), 00291+59 was confirmed as an X-ray pulsar
with a spin of 598.88 Hz (Markwardt et al. 2004a), corresponding to a spin period of
1.7 ms.
This made it the fastest known AMXP (and the 6th one to be discovered).
A sinusoidal frequency modulation of 147.4 min (2.46 h) found
in X-ray (by RXTE) was interpreted as the orbital period and no
evidence of X-ray eclipses was present (Markwardt et al. 2004b).
The fast X-ray variability was stronger and had lower characteristic
frequencies than any other neutron star LMXB, being more similar to
that of black hole systems (Linares et al. 2007).
Within days, optical and infrared counterparts were identified at RA
00:29:03.06, dec +59:34:19.0 (J2000.0) to an uncertainty of 0
5 (Steeghs et al. 2004; Fox & Kulkarni 2004). The source was bright with magnitudes of
;
;
;
.
Later observations during quiescence gave values of
in 2005 October and
in 2005 January (Torres et al. 2008a). During its decay, variability of
0.3 mag (not associated with the system's orbital period) was noted on timescales of seconds to hours in R-band (Bikmaev et al. 2005).
An optical spectrum of 00291+59 was taken on December 5 by the
William Herschel Telescope (WHT) showing ``weak evidence'' for broad
emission features at H
and He II (4686 Å) (Roelofs et al. 2004; Torres et al. 2008a),
both of which are considered diagnostics for outbursting LMXBs. A
further spectrum was taken on December 12 (still during outburst) using
the 10-m Keck 1 telescope (Reynolds et al. 2006; Fillipenko et al. 2004) which detected broad (
FWHM = 1200 km s-1) emission at H
,
H
and He I (6678 Å) as well as narrow (
FWHM = 300 km s-1) ``very weak'' features at He II
(4686 Å). The optical-IR SED during outburst was dominated by
thermal emission with a NIR excess thought to be due to synchrotron
emission (Torres et al. 2008a). Until the 2008 outburst, no further spectra of 00291+59 had been published.
On 2004 December 4, radio observations with the Ryle Telescope,
Cambridge yielded a ``probable detection'' of 1.1 mJy at
15 GHz (Pooley 2004). 12 h later, at the same frequency, this signal had disappeared below the detection threshold of 0.6 mJy (Fender et al. 2004a),
suggesting that this was indeed the transient source. A day later, the
Westerbork Synthesis Radio Telescope (WSRT) detected a radio source
with mean flux density of
mJy (at 5 GHz), with evidence of this fading during the 10 hour observation (Fender et al. 2004a). Very Large Array (VLA) observations on December 9 gave a detection of
mJy at 4.86 GHz (Rupen et al. 2004).
A reliable distance to the system is still to be determined. An upper
limit of 3.3 kpc based on X-ray absorption has been suggested (Shaw et al. 2005), however Burderi et al. (2007)
suggest a value of between 7.4 and 10.7 kpc, based upon the pulsar's
spin-down rate. Most recently, an estimate of 2-4 kpc has been
made based upon the X-ray luminosity during outburst (Torres et al. 2008a). An upper limit on the system's donor has been estimated at 0.16
,
implying it is most likely a hot brown dwarf (Galloway et al. 2005). This was calculated from its mass function and also its similarities with the best studied AMXP, SAX J1808.4-3658.
Observations of 00291+59 were taken in quiescence with the 3.6 m
Italian TNG telescope. The data were taken over 3 nights in 2005
August (V and R-bands), September (J, H, K-bands) and November (I-band) (D'Avanzo et al. 2007). Folding of I-band images on the known orbital period yielded a semi-amplitude of
mag
with minimum luminosity at phase 0, consistent with the neutron
star being behind the donor. Quiescent observations were also taken
using the WHT in Harris-I band (2006 September), which showed short-term flaring of up to 1 mag above the average magnitude of
(Jonker et al. 2008). This data was folded on the known period and, after removing the flares, a semi-amplitude of
0.06 mag was found, with the system at its brightest at phase
.
Although the 2004 outburst denoted the discovery of 00291+59, it is worth noting that Remillard (2004) retrospectively detected ``marginal'' (5)
similar outbursts from the source as detected by the RXTE All-Sky Monitor (ASM) in 1998 November (
mCrab) and 2001 September (
mCrab) suggesting an outburst from this source every
3 years.
1.3 The 2008 outburst
In 2008 August, 00291+59 again became active. This was first detected by the RXTE PCA on August 13 with a flux of 16 mCrab (2-10 keV). The PCA had been monitoring the source without detection every three days since 2008 May (Chakrabarty et al. 2008).
A previous observation on August 10 had not detected the source.
The August 13 detection was confirmed by Swift X-Ray Telescope (XRT)
observations on August 15 (Markwardt & Swank 2008). Optical photometry was obtained using the Wide Field Camera (WFC) on the Isaac Newton Telescope (INT) (Torres et al. 2008b) and Faulkes Telescope North (FTN) (Russell et al. 2008), noting an increase of 4 mag from quiescence to values comparable with the 2004 outburst (see Fig. 1). A detection in the UV was made by the Swift UltraViolet/Optical Telescope (UVOT) instrument on August 15 (Marshall et al. 2008) in the UVW2 filter (central wavelength, 1928 Å). No radio detections were reported, although the WSRT set a 3
upper limit of 0.16 mJy (Linares et al. 2008) at 4.9 GHz on August 16. After a rapid fade, a second brightening occurred
1 month
later and was detected on September 18, initially in the optical
with FTN, subsequently confirmed by Swift X-ray observations (Lewis et al. 2008b). This second outburst was just
0.2 mag fainter in
-band than the outburst peak on August 15.
![]() |
Figure 1:
FT North and INT |
Open with DEXTER |
Here, we present well-sampled optical light curves of both 2008 outburst peaks, plus a quasi-simultaneous optical spectrum, radio, NIR, UV and X-ray data. We describe the data collection and reduction in Sect. 2. We analyse the light curve morphology, the source's optical variability, spectral properties and broadband behaviour in Sects. 3-5. We present the evolution of the NIR-optical-UV SEDs (which are the most complete of any AMXP to date) and the first broadband (radio-X-ray) SEDs of 00291+59, and compare this outburst to the 2004 outburst and to other AMXPs in Sect. 6. Our conclusions are presented in Sect. 7.
2 Observations and data reduction
2.1 Optical, UV and near infrared data
Data are presented from Faulkes Telescope North (optical), Swift UVOT (UV/optical), INT (optical), an optical spectrum from Keck and PAIRITEL (NIR). The photometric observations are listed in Table 1.
To convert the magnitudes to intrinsic de-reddened flux densities, we use
for the interstellar extinction towards the source (Torres et al. 2008a) and use the extinction curve of Cardelli et al. (1989).
The resulting flux densities will have absolute uncertainties due to
the errors in the value of the extinction, corresponding (since
extinction increases at shorter wavelengths) to 3%, 28% and 109% of the
flux at K, V and UVW2 wavelengths respectively (adopting values of
A K = 0.114 and
AUVW2 = 2.664 from Table 3 and Eq. (4) in Cardelli et al. 1989). These errors are taken into account for the SEDs presented in Sect. 6.
Table 1: Photometric observations (NIR, optical and UV detections) of 00291+59 used in this work.
2.1.1 Faulkes Telescope North
00291+59 has been included as part of a monitoring campaign of 30 LMXBs including 5 AMXPs (Lewis et al. 2008a)
using the two Faulkes Telescopes. Data have been collected using the
2-metre robotic Faulkes Telescope North located at Haleakala on Maui.
It currently uses a Merope camera (EM01). Prior to 2008 August, the
camera used was the Apogee ``Hawkcam'' (EA01). Both cameras were
coupled with an E2V CCD42-40DD CCD giving a 4
7
4
7 field of view and producing images of
pixels binned
to give
pixels at 0.278 arcsecond pixel -1.
Science images are produced using the Faulkes automatic pipeline which
de-biases and flat-fields the raw images. Filters used are B, V, R (Bessell),
(Sloan-SDSS) and y (Pan-STARRS). We monitored the source in
-band, typically acquiring a 200 second exposure every
2 weeks
when the source was visible from FTN. When the outburst was detected on
2008 August 1, we increased the sampling to one observation every
2 days and used B, V, R,
,
y filters to investigate colour changes. This shorter cadence
continued until November 10, near the end of the fade of the
second outburst peak.
Seeing values range from 0
4 to
3
1.
Images were discarded if the signal-to-noise ratio was low (often due
to thin cloud), or if the tracking or focus were poor. During the
outburst, we detected the source in 133 images taken between 2008
August and November. We performed aperture photometry of 00291+59 and
two nearby field stars (see Table 2 and Fig. 2) using the aperture photometry package APPHOT in IRAF.
Point-spread-function (PSF) fitting was not used since minor errors in
tracking on some observations resulted in non-circular or
non-elliptical stars for which aperture photometry was more suitable. A
fixed aperture radius of 6 pixels along with a background annulus
of 10-20 pixels was adopted for all three stars in all filters.
Flux calibration in B, V and R-bands was achieved using photometry of the standard star field PG 0231+051 from the list of Landolt photometric stars (Landolt 1992),
which are observed regularly by FTN. Both the standard and the 00291+59
fields were observed on 2008 September 4, 6 and 14. Four
stars in the same field as PG 0231+051 with known B, V, R magnitudes were also used. Accounting for differences in airmass between target field and standard, we calculated the B, V, R magnitudes
of the two chosen field stars in the 00291+59 field. From the three
dates used, the measured magnitudes varied by only
0.01 mag,
indicating the conditions were photometric on all three dates.
Magnitude errors were estimated from the (night-to-night) range of
measured magnitudes of the two field stars and measured differences in
the relative magnitudes of the stars in the PG 0231+051 field.
Table 2: Optical and NIR magnitudes of the two comparison stars in the 00291+59 field used for flux calibration.
For the -band
flux calibration, the field centred on RA 23:48:20, Dec +00:57:18
(J2000.0) was observed 9 times on 2008 September 27, at airmasses
from 1.58-2.04. Instrumental magnitudes were measured for the star SDSS
J234817.22+005557.2 and its values compared with that given in the
Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6). Comparisons were
made between this star and others in the field to ensure that it was
non-variable and that conditions remained photometric. Having
calculated the change in offset (between instrumental and known
magnitudes) with airmass, and accounting for the difference in exposure
times between this and the target field, we were able to derive
magnitudes for stars 1 and 2 as shown in Table 2. The uncertainties include the instrumental uncertainties given by IRAF
for stars 1, 2 and SDSS J234817.22+005557.2, summed in
quadrature with the value quoted for SDSS J234817.22+005557.2 in SDSS
DR6.
Spectroscopic standards were required for the y-band calibration, since magnitudes in this Pan-STARRS y filter (which has an effective wavelength of 1.004 m) are not typically known for photometric standard stars. The spectroscopic standards BD+28 4211 and Feige 110 (Oke 1990)
were both used on October 27 (4 and 5 exposures,
respectively) and October 31 (one exposure of each). Both
standards have smooth spectra with known AB magnitudes in steps of 2 Å. Using the same method as for the
-band, we estimated the y-band
magnitudes of the two 00291+59 field stars. The uncertainties
incorporate the small differences in measurements between the two
dates.
![]() |
Figure 2:
FT North finder chart (SDSS- |
Open with DEXTER |
Table 3:
Swift UVOT optical/UV magnitudes and 3
upper limits of 00291+59.
2.1.2 Swift UVOT
The Swift UVOT observed 00291+59 on nine dates in 2008. Four detections (at 5)
of the source were made during the 2008 outburst (see Tables 1 and 3); two in B-band and two in UVW2,
which has an effective wavelength of 1928 Å. Individual images
were combined and magnitudes and upper limits were derived using the
standard Swift UVOT routines provided by NASA's High Energy
Astrophysics Science Archive Research Center (HEASARC). Apertures of 3
were used centred on the known coordinates of 00291+59. As a confirmation, we measured the B magnitude of field star 1 to be consistent with that measured by FTN (see Table 2)
to an accuracy of <0.1 mag. On some dates, >10 images
were acquired in the same filter, but the counts of 00291+59 were never
high enough to produce a light curve. In Table 3 we list the UVOT magnitudes when detected and upper limits for non-detections.
2.1.3 INT
Observations of 00291+59 were taken using the WFC at the 2.5 m
Isaac Newton Telescope (INT) at Roque de los Muchachos, La Palma,
Spain, between 2008 August 15 and 19. The WFC has an image
scale of 0.333 arcsec pixel -1. Observations were made using the Sloan
and
filters. A total of 188 observations were taken in
-band over the 5 nights, with a further 92 taken in
-band
over the first 3 nights. A master bias and flat-field were created
and used to reduce the raw images. Photometry was performed in the same
way as for the FTN data (Sect. 2.1.1) with the same comparison
stars, apertures and annuli (in pixels). Flux calibration for
-band was performed as for the FTN
-band data; for
-band, the calibration values of Torres et al. (2008a) were used (see Table 2).
2.1.4 Keck spectrum
A single spectrum was acquired on 2008 Aug. 28 using the Low
Resolution Imaging Spectrometer (LRIS) attached to the Keck 1
telescope. 00291+59 was observed using a 1
0
wide slit with the CCDs on each spectrograph arm unbinned. The exposure
time for both blue and red arms was 900 s. The red arm was used
with the 400/8500 grating. The useful wavelength range covered the
5470-9240 Å interval with a dispersion of 1.86 Å per pixel
and a spectral resolution of about 6 Å FWHM. The blue arm
was used with the 600/4000 grism to yield an useful wavelength
coverage of 3890-5600 Å with a dispersion of 0.62 Å per
pixel. The blue arm data were too noisy to obtain a meaningful
spectrum. Conditions were not considered to be photometric; seeing was
variable from 0
7-1
0. The spectrum was extracted and reduced using KPNOSLIT in IRAF with 1
slit and using the standard star, Feige 110. We use the Starlink page DIPSO to measure equivalent widths.
2.1.5 PAIRITEL
NIR photometry was obtained with the 1.3m Peters Automated Infrared
Imaging Telescope (PAIRITEL) at Fred Lawrence Whipple Observatory
(FLWO) (Bloom et al. 2006) which images simultaneously a 8
5
8
5 field of view in the J, H and Ks
photometric bands. The observations consisted of a large number of
dithered 7.8 s exposures on source. The dithered exposures were
first bias and flat-field corrected to be mosaiced together for each
individual visit (see e.g. Blake et al. 2005).
The observations were made on 8 nights between 2008 September 21 and
October 13, although observations after October 6 yielded
non-detections. Photometry was performed using the same comparison
stars and method as previously described, however with apertures of
radius 3 pixels to compensate for the greater image scale. The J, H, K 2MASS magnitudes of the two field stars (Table 2) were used to calibrate the flux of 00291+59. The magnitudes and upper limits are given in Table 4.
2.2 X-ray data
We use pointed observations from the Swift XRT, XMM-Newton and RXTE PCA and augment these with publicly available data from Swift Burst Alert Telescope (BAT) and RXTE ASM. We obtain a light curve from the RXTE PCA data (14 detections in August and 13 in September). We list the X-ray observations in Table 5.
Table 4: PAIRITEL NIR magnitudes and upper limits of 00291+59.
2.2.1 Swift XRT
We analysed all Swift XRT observations of the field taken between 2008
August and October. We performed the standard screening on all
observations using the latest ``xrtpipeline'' (v. 0.12.3). In each of
the two observations where 00291+59 was detected we extracted a source
spectrum. In both cases, the XRT was in photon counting (PC) mode
and the observations suffered from pile up. To correct for this,
following the Swift Science Data Centre (SSDC) recommendations, we used an annular extraction region centred on 00291+59's position with inner radii of
20
and
10
for observations 00031253001 and 00031253005, respectively.
We thereby excluded the core of the PSF from our analysis, resulting in encircled energy fractions at the inner radius of 50% and
40% for 00031253001 and 00031253005, respectively. The annuli's outer radius was fixed at
120
.
Background spectra were accumulated far from the source using regions
of the same size. We created exposure maps and ancilliary files using
``xrtexpomap'' (v.0.2.5) and ``xrtmkarf'' (v.0.5.6), respectively,
thereby accounting for the energy lost in the PSF core. We grouped the
spectra to have a minimum of 20 counts per energy bin and fitted them
within Xspec (v. 11.3.2ag) using the latest available response
matrices (swxpc0to12s6_20010101v011.rmf) and the 0.5-10 keV energy
range. The spectra were satisfactorily fitted (reduced
2
of 0.9 for 77 and 33 degrees of freedom) with a simple absorbed
power law model. The measured fluxes, absorption and power law
parameters are given in Table 6, along with 3
flux upper limits when the source was not detected.
2.2.2 XMM-Newton
We analysed the X-ray Multi-Mirror satellite (XMM-Newton) observation
of 00291+59 taken on MJD 54703 (2008 August 25), which lasted
approximately 33 kiloseconds. 00291+59 was clearly detected in both the
MOS1 and MOS2 detectors, which operated in large window (medium filter)
mode, at an average net rate of 0.01 counts s-1.
We excluded contaminating flares from the observation, defined as
segments with a count rate above 10 keV higher than 1 count
second-1. We then defined source and background circular regions, both with a radius of 30
.
We took into account the extraction area using the tool ``backscale''
and created response matrices and ancilliary files with ``rmfgen'' and
``arfgen'', respectively. We grouped the channels in the resulting
spectra in order to have a minimum of 15 counts per energy bin,
and fitted them within Xspec using the 0.5-10 keV energy range,
using an absorbed power law model. We note that, as with the majority
of the Swift XRT observations (see Sect. 2.1.1), the low count
rates from the source preclude us from studying variability in detail
within the observation window, although we note that the level of
variability is consistent with the source being at a constant flux.
2.2.3 RXTE PCA
The PCA was used to monitor the X-ray flux during both outburst peaks (Chakrabarty et al. 2008; Galloway et al. 2008).
We use the 16 s time resolution Standard 2 mode data to calculate
the X-ray flux in the 2-16 keV energy band. The energy-channel
conversion is done by using the pca_e2c_e05v02 table provided by the
RXTE team. The deadtime was corrected and the number of background
events within our energy range and time interval was estimated using
the FTOOL ``pcabackest'', following the standard procedure suggested by
the RXTE manual. We reject all data with a measured flux (2-16 keV) of
erg cm-2 s-1
because flux from background contamination was likely significant. This
is most evident from an apparent PCA detection on MJD 54700 of
erg cm-2 s-1 whereas a Swift XRT pointing fours hours earlier did not detect the source, with an upper limit of
erg cm-2 s-1. The nearby intermediate polar V709 Cas is in the PCA's field of view, and contributes some contaminating flux (see Linares et al. 2007) however the flux level of this source is lower than the above flux limit, and therefore its contribution is minimal.
2.2.4 Publicly available data from X-ray monitors
Five 3
detections were also made by the Swift BAT instrument (15-50 keV).
All five detections were during the second peak, in 2008 September; BAT
did not detect the source during the first outburst peak. We include
these public data in our analysis. We also include four 3
detections from the public archive of the RXTE ASM (1.5-12 keV);
three from the first outburst peak and one from the second. All RXTE
ASM fluxes were converted to unabsorbed (adopting
cm-2 from Torres et al. 2008a)
2-10 keV fluxes using the HEASARC tool WebPIMMS to compare with
the Swift XRT and XMM 2-10 keV data. A power law index of 1.6 (as
measured by Swift; see Table 6) was adopted.
2.3 WSRT radio data
00291+59 was observed with the Westerbork Synthesis Radio Telescope on 2008 August 15-16 (Linares et al. 2008). The observation was made between 20:30-06:30 UT
at the median frequency of 4.9 GHz, with a total bandwidth of
160 MHz. The primary calibrator used was 3C 286. The
calibration and analysis of the data were done using MIRIAD (Sault et al. 1995).
No object was detected at the position of the target reported by Rupen et al. (2004). The 3
upper limit to the flux density (measured in the image plane) was 0.16 mJy.
![]() |
Figure 3:
Upper: optical, UV and NIR data from the 2008 outburst. Lower: X-ray and radio data from the same period. We also plot earlier FTN |
Open with DEXTER |
3 Analysis of the 2008 outburst
The multi-wavelength light curve of the 2008 outburst is shown in Fig. 3.
The source displays two separate periods of activity, which we describe
as ``peaks'' within the overall ``outburst'' (Fig. 1).
We display the initial X-ray detections (MJD 54691) of the outburst and
its subsequent detection in the optical/infrared (OIR) at
MJD 54693. FTN observations from long-term monitoring (Lewis et al. 2008a) had provided non-detections, most recently on August 3 and 9 (
and 21.39 respectively), showing that the OIR brightened by at least 0.5 mag day-1.
We show a well-sampled UV-IR fade over the following 30 days
accompanied by RXTE PCA and Swift XRT monitoring. For the first three
days, the X-ray had an approximately constant flux (the optical flux
also seemed approximately constant) before the source faded, more
rapidly in X-ray than in OIR. We present data from Swift XRT
observations in Table 6. The
XRT X-ray spectra are hard, and can be fit with an absorbed power law
with a photon index typical of AMXPs in outburst. Just 11 days
after the first X-ray detection with Swift XRT, an XMM-Newton detection
showed it had drastically faded. From XMM simultaneous fitting of MOS 1
and 2 spectra using the ``wabs*powerlaw'' model, we derive an
unabsorbed 2-10 keV flux of
erg cm-2 s-1,
i.e. more than four orders of magnitude lower than the outburst peak
fluxes. Further details of the XMM spectrum will be reported in Linares
et al. (in preparation). We note that long after the source was no
longer detectable with Swift, the source remained visible in the OIR (
-band), eventually approaching its quiescent value at
MJD 54722 (31 days after the outburst was first detected).
![]() |
Figure 4:
Upper: comparison of the three fades using R-band data from 2004 (Torres et al. 2008a) and |
Open with DEXTER |
A second ``peak'' was detected; this time, first in optical (Lewis et al. 2008b), with the source re-brightening by 4-5 mag to 18.3 (MJD 54727) in
-band,
very similar to the magnitude of the first peak in August. It is clear
that the 2008 August and September peaks differ, with the August peak
fading rapidly and the September peak remaining in a bright `plateau'
at near maximum luminosity for
10 days. This second peak lasted
70% longer than its predecessor (i.e.
50 days as opposed to
30 days). Whereas the
-band
light curve appeared fairly constant during the September plateau, the
X-ray and NIR light curves brightened then faded slightly, before
eventually all bands faded towards quiescence. In the NIR, the 2004
outburst was estimated to reach a maximum of
(Steeghs et al. 2004) from a quiescent
;
in 2008, we measure (during the second peak)
.
This is slightly at odds with the evidence in the R-band which suggests the 2008 outburst was slightly fainter than that of 2004, however we note that the 2004 K-band magnitude was from observations 4 days after the initial detection.
From Fig. 3 and Table 5,
we can see that the RXTE ASM (2-10 keV) detected the first peak on
three occasions, and the second peak only once. In contrast, Swift BAT
(15-50 keV) did not detect the first peak, and yet detected the
second peak on four occasions. However, we note that the respective
sensitivities of the ASM and BAT are variable in time, and indeed we
see that the photon indices as measured by Swift XRT were consistent
with being the same at the start of both peaks (see Table 6). In Fig. 4,
we plot data for the 2004 outburst alongside the two separate peaks of
the 2008 outburst, for the optical (upper panel) and X-ray (lower
panel) fades. The 2004 outburst was measured at its maximum at
(Fox & Kulkarni 2004), as compared with the 2008 outburst, which reached values of
.
In the optical, there is a rapid fade between days
12-20
after all three peaks and the same is evident in X-ray up to at least
day 16. It is remarkable how identical the three X-ray light curve
decays are, given that the initial phases (before day 10 in
Fig. 4) are so different. The
light curves are overlaid such that these fades are aligned (the same
shifts are applied in X-ray and optical). Within the fade from the
first peak, we see possible evidence of optical reflares of
1 mag (
day 20-30, Fig. 4).
This post-peak flaring activity has only been observed from one other
AMXP; SAX J1808.4-3658 in its X-ray light curve, where the system
flares by several orders of magnitude weeks after the initial outburst (Patruno et al. 2009).
The 2008 August fade from peak to quiescence (in -band) is measured as
4.6 mag over a period of 29 days giving an average fade of
mag day-1. The initial fade (from day 2 to day 5 of the outburst) is steeper with a fitted value of
mag day-1. The rise to the second peak over
5 days was measured as
0.72 mag day-1. This second peak remained in a ``plateau'' for
10 days, before fading to quiescence at
mag day-1 over the following 16-28 days (the range being due to the under-sampling of this part of the light curve).
The recurrence of the second peak to a similar maximum flux within 35-40 days
is unique behaviour for an AMXP, although there are some similarities
to the persistent source, HETE J1900.1-2455, which has been in outburst
since its discovery in 2005 (Elebert et al. 2008). Secondary maxima have been noted previously for other LMXBs (e.g. Chen et al. 1997)
who suggested an increase in accretion rate due to re-processing of
X-rays from the initial outburst heating the disc and/or donor star.
Alternatively, Truss et al. (2002) and Truss (2005)
show that tidal instabilities in the disc could cause multiple outburst
peaks, but that this phenomenon is not as relevant for systems with
large mass ratios such as 00291+59, for which the neutron star is at
least a factor of
10 more
massive than its companion. A further scenario is that of propagation
of cooling/heating waves in a disc that may cause additional maxima in
light curves (e.g. Lasota 2001). We note that the
smaller scale of AMXPs (shorter orbital periods and smaller discs) and
their shorter duty cycles (as opposed to ``classical'' LMXBs) means
that we might expect their discs to both empty and re-fill more
quickly, explaining their more rapid duty cycle. However, the previous
outburst of 00291+59 was in fact brighter (in both X-rays and optical)
but did not display this double-peaked behaviour.
Table 6: Swift XRT X-ray fluxes.
![]() |
Figure 5:
Upper: 2008 outburst |
Open with DEXTER |
Table 7: Variability of 00291+59 and a faint field star.
4 Optical variability - orbital, long-term and short-term
We have monitored 00291+59 since 2007 August and prior to 2008 mid-August, had not detected the source with FTN in -band, most notably immediately preceding the outburst on August 3 and 9. In 2008 August and September, 00291+59 reached
on two separate occasions. For comparison, its magnitude near the time of discovery in 2004 was
.
Between the two 2008 peaks, the data are consistent with the source fading to quiescence in
-band (within error bars); the quiescent level in 2004 being
(D'Avanzo et al. 2007; Jonker et al. 2008).
During the August peak, we collected data from the INT in both
and
-bands (see Fig. 5), over timespans of
3-4 h (August 18-19) and runs of
30-60 min per filter (August 15-17 inclusive). We also observed 00291+59 intensively in
-band around the September peak for
4 h and
3 h respectively with FTN. Whilst we note fluctuations of
0.05 mag from one observation to the next (with a cadence of
240 s), we saw no evidence on any occasion for any modulation associated with the orbital, or any other period.
This lack of evidence of orbital modulation is likely the result of the light from the accretion disc swamping much of the lower amplitude orbital variability, and had been noted in the previous outburst (Bikmaev et al. 2005; Torres et al. 2008a; Reynolds et al. 2006). In quiescence, optical observations have twice detected sinusoidal orbital modulation (D'Avanzo et al. 2007; Jonker et al. 2008) consistent with emission from an irradiated companion star; this is similar to the example of the AMXP, XTE J1814-338 (D'Avanzo et al. 2009), but opposed to SAX J1808.4-3658 which displays orbital modulation in both outburst and quiescence (Elebert et al. 2009). From the lower panel of Fig. 5, we see that the amplitude of this modulation seems to be greater as the system fades, even though these three datasets are taken on consecutive nights.
In Table 7, we show the level of variability of 00291+59 and a nearby faint field star. Table 7
demonstrates that the variability in each of the INT runs increases as
the system fades. We note that we see significant variability in -band (3
)
on the first night only (in 2008 August) and that we see less significant variability (2
)
on 3 of the other 4
-band observing runs. For the three
-band runs in 2008 August and both
band runs in 2008 September, we detect variability at less than 2
.
Since the 2008 outburst, we have continued to monitor 00291+59, but to
date (March 2010), have seen no further activity in the system.
5 Spectral analysis
We present a new Keck 1 spectrum (Fig. 6, taken in 2008 August) taken when the source was at R = 21.1 (i.e. when the system was 3 mag
fainter than at its maximum, but still 2 mag above quiescence).
Since the blue arm suffered from poor signal-to-noise, we display the
red arm spectrum only, meaning that we are unable to compare our
spectrum with previous detections of H
,
H
or H
or the He sc ii feature at 4686 Å (Torres et al. 2008a). We note a prominent double-peaked H
feature at 6563 Å with a measured equivalent width (EW) of
Å and
Å (1600 km s-1), indicative of a rotating accretion disc (Frank et al. 2002). The peak-to-peak separation is 18 Å (
820 km s-1) and we see the blue-shifted peak being
48% stronger than that of the redshifted peak. Double-peaked H
emission has previously been seen in this source Torres et al. (2008a), and also in the outbursting AMXP, XTE J1814-338 (Steeghs 2003). Using our peak-to-peak value, and following the method in Sect. 10.2 of Torres et al. (2008a), we derive an inclination of
35
,
which is consistent with their measurement of 27
.
We do not detect a He I feature at 5875 Å which Torres et al. (2008a) detected weakly, nor do we see a He I feature at 6678 Å which
had been previously noted in one single 300 s spectrum (Fillipenko et al. 2004). At 7065 Å, we detect a line with low signal-to-noise, which is likely the He I line as reported in Torres et al. (2008a) from the spectrum in Reynolds et al. (2006).
We are also unable to detect the interstellar Na D doublet
(5889 Å, 5895 Å) but do detect a Diffuse Interstellar Band
(DIB) at 6284 Å (
Å) as in Torres et al. (2008a) as well as telluric features at 6864 Å and 7600 Å.
Table 8 summarises the lines measured in all the published optical spectra of this source. From it, we can see that the H
EW is larger at lower optical and X-ray luminosities. We note that this relationship between H
EW and X-ray luminosity is consistent with the previously found trend for neutron star and black hole XRBs (Fender et al. 2009).
![]() |
Figure 6:
Red arm spectrum from Keck 1, DIB denotes Diffuse Interstellar Bands. For clarity, we mark the wavelengths of He I
lines at 6678 Å and 7065 Å (the latter being a marginal
detection in this spectrum). Insert: zoom in on double-peaked H |
Open with DEXTER |
Table 8: Spectra equivalent widths (Å).
6 Spectral energy distributions (SEDs)
In Fig. 7 we present the
de-reddened (see Sect. 2.1) SEDs compiled from the
quasi-simultaneous NIR, optical and UV data from 2008 August (top left
panel) and September (top right). These are some of the most complete
SEDs ever compiled in this frequency range for neutron star LMXBs (up
to eight flux densities measured, spanning one order of magnitude in
frequency). We have created SEDs whenever the data were acquired within
the same day. The errors in frequency (x-axis) represent the frequency ranges of the filters. For each data point the error in the flux density (y-axis) is derived from the error in the magnitude. We adopt
;
if the extinction is incorrect by 0.3 mag (see Sect. 2.1),
this propagates into additional errors in intrinsic flux density which
are indicated by the separate error bars (crosses) shown across the top
of the top right panel, noting that bluer wavelengths are affected more
than redder ones. The effect is systematic; if the extinction is
then the intrinsic flux densities were underestimated by the amount indicated by these error bars.
Taking into account these uncertainties, we see that there is generally a blue optical-UV SED (
,
where
), with a separate NIR excess which dominates the
y, J, H, K filters, moreso at higher flux
densities. The blue component can be successfully modelled by a
single-temperature blackbody of temperature
K in both the August and September peaks (Fig. 7, bottom panels), which is what we might expect from an illuminated or viscously heated accretion disc (e.g. Hynes 2005).
We are likely sampling a region of the blackbody just redward of the
peak, blueward of the Rayleigh-Jeans tail. However, the errors due to
uncertainties in the extinction could significantly change the slope of
the SED, and hence the derived blackbody temperature. If the extinction
was higher than the measured range of
in Torres et al. (2008a), the intrinsic disc spectrum would be bluer than expected from a Rayleigh-Jeans tail, hence confirming the upper limit of
.
On the initial fade in August, the SED appears less smooth than for the September decline, in particular a bright y-band flux density on August 28. This may be due to rapid variability, which we do detect to be stronger in the August peak compared to the September peak (see Sect. 4).
![]() |
Figure 7: Top row: UV, optical and NIR fluxes for the 2008 peaks. Bottom row: single black-body plus jet contribution model plotted against fluxes for data at each peak. |
Open with DEXTER |
![]() |
Figure 8: Left panel: model for two black-bodies at second outburst peak. Right panel: model showing the dependence on flux for different temperature disc. |
Open with DEXTER |
Our October 4 (MJD 54743) observations show that the blackbody is fading, and can either be modelled (Fig. 8,
right panel) by a cooler disc of the same area or a hot, smaller disc
(implying a change in outer disc size or changes in a disc warp,
resulting in a smaller total disc area but at the same temperature; see
e.g. Hynes et al. 2002). At lower flux
densities later on in the decline of the outburst (October 4, 6),
the SEDs can be approximated by a blackbody of possibly lower
temperature (
K). On October 20 (MJD 54759), the source was about one magnitude above its quiescent level in R-band, and at this point, we cannot easily constrain the spectrum, although in Fig. 8 (right), we overplot a cooler blackbody of
K.
A NIR excess has been reported for almost all AMXPs in outburst. This excess is most apparent near the outburst peak and is always absent at lower luminosities (see Russell et al. 2007 and references therein). It has typically been attributed to synchrotron emission, likely from the jets in the system (e.g. Giles et al. 2005; Wang et al. 2001; Krauss et al. 2005). The NIR excess has sometimes been seen quasi-simultaneously with a radio detection; this was the case for 00291+59 in the 2004 outburst (Torres et al. 2008a).
Here, for the 2008 outburst we also identify a NIR excess above the
disc blackbody which seems to disappear later in the outburst when the
source fades towards quiescence. However, the excess is not consistent
with optically thin synchrotron emission, for which we would expect a
red SED with
.
Instead, we see an excess which is also blue (
), most apparent by a fainter K-band flux density compared to that in H or J. The uncertainties in the extinction cannot explain this fainter K-band
flux density. We therefore conclude that the NIR excess in 2008
September is not due to optically thin synchrotron emission.
![]() |
Figure 9: Evolution of the SED of 00291+59. Left panel: broadband SEDs from the 2004 outburst and both 2008 outburst peaks. Right panel: blackbody plus jet contribution shown to approximate data from left panel. |
Open with DEXTER |
In the lower panels of Fig. 7
we attempt to model the NIR excess from the outburst peaks in August
and September by a simple jet model. Jets in AMXPs usually have a
slightly inverted (
)
radio-to-optical SED (Russell et al. 2007) from the optically thick (self-absorbed) jet. We take this value of
for the optically thick jet, and
for the optically thin jet. If the break between optically thick and optically thin lies around the H-band
then the NIR-to-UV SED from September can be explained by this jet
(producing the NIR excess) and the blackbody (producing the blue
optical-UV emission). Were the break at a lower frequency than H-band, we would expect the K-band
to be brighter than observed since it would now be optically thin.
Equally, were the break at a higher frequency, then we would expect a
smooth increase in flux density (from
0 to
)
as the optically thick jet spectrum joins the disc spectrum. Instead, the low K-band flux density shows that this cannot be the case (although, see below). Only a jet break around H-band
can produce the NIR excess seen in the SED (on several dates), if we
assume the jet is responsible for this NIR excess. The measurement of
the break in the jet spectrum is independent of the spectral indices of
the jet spectrum's optically thick and thin regimes. Although the data
are consistent with this simple model (it is not a fit, but an
approximation), this by no means shows that the origin of the NIR
excess is the jet. More sophisticated modelling of the broadband data
(e.g. Maitra et al. 2009) are required to
confirm whether the NIR excess is consistent or not with a jet origin.
Some models predict an ``optical bump'' in the SED of the synchrotron
emission from the jet (Markoff et al. 2005; Pe'er & Casella 2009) which could be consistent with these data.
In the left panel of Fig. 8 we demonstrate that the NIR excess can also be modelled with a second blackbody component. The fainter K-band flux density is more satisfactorily explained in this case, but the blackbody must have a temperature
K
and have an area around 60 times that of the accretion disc. The
temperature is consistent with (the photosphere of) a low-mass star but
the size is not. If the higher temperature blackbody, which dominates
the optical, can be explained by the underlying viscously heated disc,
the NIR excess may be explained by irradiation of the disc. However,
the inferred temperature of this blackbody (
K)
is far too cool for this to be the case. Alternatively, part of the
disc may produce optically thin emission, where the free-free emission
has a flatter spectra than the blackbody from the disc. We may expect a
circumbinary disc (e.g. Blundell et al. 2008)
surrounding the system to have a temperature and size of this order,
but we would not expect a circumbinary disc to be transient; the
quiescent NIR flux is
two orders of magnitude fainter than this excess. It could be
illuminated by the X-ray photons from the central source, but the NIR
excess varies more rapidly than the disc, and so seems unlikely.
Quasi-simultaneous, broadband radio-to-X-ray SEDs are presented in Fig. 9. We include data from the 2004 outburst when there were radio detections of the source, because although the data have been published, the broadband SEDs have not. Radio data are taken from Pooley (2004); Fender et al. (2004a), Rupen et al. (2004) and NIR, optical and X-ray data are from Torres et al. (2008a) and the RXTE ASM.
During the 2004 outburst, a radio counterpart was detected on three dates (peaking at 1 mJy) and a blue optical/IR SED with an IR-excess was observed at the same time (Fig. 9, top left panel; see also Torres et al. 2008a). The radio-to-optical SED can be reproduced by the simple blackbody + jet model, adopting
for the optically thick jet spectrum and a jet break in the mid-IR. The
temperature of the blackbody is similar to that found near the peaks of
the 2008 outburst, but no UV data were acquired in 2004. Again, this is
simply a demonstration that a jet + disc model can satisfactorily
reproduce the radio-to-optical SED. We are not claiming this is the
broadband spectrum of the jet, or the precise temperature of the disc.
For example, it could be that the optically thick jet has a steeper
spectral index of
and the jet break occurs at lower frequencies. On the other hand, the optically thick jet cannot have a spectral index of
if it is to explain the radio data and the IR-excess because the latter
component is too bright. It is unlikely that the jet model could
explain the X-ray data - using
for the optically thin jet, this component is at least one order of magnitude fainter than the observed X-ray flux. Only if
for the optically thin jet (which is shallower than typically expected
for optically thin synchrotron emission) could this component dominate
the X-ray flux. This is similar to the result obtained for the neutron
star LMXB 4U 0614+09, for which the synchrotron-emitting jet
likely contributes no more than
1% of the X-ray 2-10 keV flux (Migliari et al. 2006,2010).
In the lower six panels of Fig. 9
we present the radio-to-X-ray SEDs of both declines of the 2008
outburst. The unabsorbed X-ray spectra are provided by the Swift XRT
(see Sect. 2.2.1). The simple blackbody + jet models from close to
the outburst peaks are shown in the right hand panels. We again find
that if the IR-excess originates in the jet, the synchrotron jet is
unlikely to produce the majority of the X-ray emission. If the
optically thin jet is extrapolated from the IR-excess to X-ray with
,
the jet could contribute up to
10% of the X-ray flux. In the lower right panel of Fig. 9, we show an alternative jet model, jet (ii), in which the optically thick-thin jet break is at a higher frequency (1015 Hz) and the optically thin spectrum has an index of
.
This jet is able to reproduce the X-ray power law (the index being the same as that measured by Swift XRT; see Table 6),
but struggles to reproduce the shape of the IR-excess and disc
blackbody as well as the first jet model (i) does. We stress that
we are unable to make a direct measurement of the synchrotron jet
contribution to the X-ray flux; we simply note that with these data, it
is possible to reproduce the IR-excess and X-ray powerlaw with a simple
jet model, and that more complex modelling (e.g. Markoff et al. 2005)
and probably more complete SEDs are required to make any solid
measurement. The radio upper limit near the Aug. 2008 outburst
peak puts a solid constraint on the radio-to-optical jet spectral index
of
,
if the IR-excess originates in the jet, which is consistent with the average value for AMXPs of
(see discussion above). In the third right panel the second, cool
blackbody approximating the IR-excess in 2008 September is also shown.
At first glance it appears the radio jet is intrinsically
fainter in 2008 August compared with 2004 December. We investigate this
possibility in Fig. 10,
in which we plot the radio flux density against the quasi-simultaneous
X-ray flux. We find that on the date of the radio observation in 2008
August, the X-ray flux (2-10 keV) is fainter by a factor 2
than on the three dates in 2004 on which a radio detection was made. In
addition, of the three dates when radio emission was detected, the
brightest radio flux was observed on the date of the brightest X-ray
flux.
Figure 10 suggests there may be a relation between radio and X-ray fluxes for 00291+59. The slope of a correlation (,
where
)
is expected to be steeper for neutron star XRBs compared to black hole XRBs (see Gallo et al. 2003 and references therein). This has been observed with some exceptions (Migliari & Fender 2006; Tudose et al. 2009). Migliari & Fender (2006)
include data of 00291+59 in their analysis but they do not include all
the radio detections made in 2004 (or the 2008 upper limit). From the
three radio detections, a very steep correlation slope is inferred;
.
If we neglect the brightest radio detection at 1.1 mJy (see
below), the slope is consistent (within errors) with what may be
expected (e.g. Migliari & Fender 2006);
(for NS XRBs) or
(for BH XRBs). Equally, these slopes are consistent with our
non-detection at a lower X-ray flux than those observed in 2004.
However, from these data, the correlation slope is poorly constrained.
In black hole XRBs, a bright radio flare is often seen when the source makes a transition to the soft X-ray state, possibly due to a fast, discrete jet ejection slamming into slower jet material launched previously (Fender et al. 2004b). Although 00291+59 remained in a hard state throughout its 2004 outburst (Fender et al. 2004b), the 1.1 mJy radio measurement is a factor of three times brighter than the other radio detections (at similar X-ray fluxes). This suggests a discrete jet flare could have been responsible for the brightest detection (as opposed to the supposedly steady, hard state jet). Evidence for ejections in hard states has been found in other XRB transients, e.g. GS 1354-64 (Brocksopp et al. 2001), XTE J1118+480 (Brocksopp et al. 2010), H1743-322 (Jonker et al. 2010).
The upper limit of 0.16 mJy in 2008 August is consistent with a positive radio-X-ray correlation in 00291+59.
![]() |
Figure 10:
Radio flux density against quasi-simultaneous (within one day) X-ray
flux for 2004 and 2008 outbursts. The filled circles represent the
2004 detections (Pooley 2004; Rupen et al. 2004; Fender et al. 2004a);
the cross represents the 2008 non-detection (this paper). The solid
line (without ticks) shows the power-law fit to all three
2004 detections. For comparison, slopes of
|
Open with DEXTER |
7 Conclusions
We present a unique multi-wavelength dataset of the 2008 double-peaked
outburst of the AMXP, IGR J00291+5934, from radio to X-ray, with
coverage of both outburst peaks (2008 August and September). We show
that, optically, the September fade is similar in shape to the previous
2004 outburst, however we note a difference between the 2008 August and
September peaks in that the former exhibits a more rapid fade, while
the latter displays a distinct plateau phase. This plateau, at near
maximum luminosity, lasted 10 days
before the system faded over a longer timescale than in August. The
plateau was also evident in the X-ray light curves; for comparison, the
X-ray flux faded by >4 orders of magnitude in <10 days
during the first outburst peak in August. From Fig. 4,
we see that the morphology of the X-ray fade in both peaks, and in
2004, are very similar. The re-brightening of the system within a few
weeks in 2008 is unique within AMXPs, and at odds to this source's
previous outburst.
We see short-term variability in the optical of 0.05 mag within the outburst down to a timescale possible of
240 s,
but no periodicity, on the previously published orbital period. This
suggests changes in the accretion flow (likely caused by the
irradiation of the disc) dominate over ellipsoidal modulations due to
the companion star when the source is in outburst. In addition, the
optical variability during the September plateau was of a smaller
amplitude than the August peak, at similar flux levels.
We study the light curve morphology and evolution, presenting the first
radio-X-ray spectral energy distributions for this source and the most
detailed UV-IR SEDs for any outbursting AMXP. In the optical, the SEDs
contain a blue component, which can be fitted by a blackbody, likely
from the disc (at maximum luminosity,
K)
and a transient near-infrared excess. This excess is consistent with a
simple model of a synchrotron jet (as seen in other outbursting AMXPs),
however we cannot exclude other potential origins such as the presence
of a second, cooler (
K)
transient blackbody of uncertain origin. We display UV-IR SEDs for
12 dates during the outburst, and find that the IR excess fades
more rapidly than the optical disc flux. Our optical spectrum shows the
double-peaked H
profile of an accretion disc but we do not clearly see other lines (e.g. He I,II) that were observed in 2004.
The lack of a radio detection prevents us drawing many conclusions
about the long wavelength part of the SEDs. However, our upper limit is
consistent with an optically thick jet (with
)
between radio and NIR, as seen in other AMXPs. Our radio-X-ray SEDs suggest that there is a break in the H-band
between an optically thick and optically thin jet. Our simple modelling
is also able to account for the radio-optical SED from 2004 which also
consists of a blue disc and a NIR-excess (Torres et al. 2008a). Extending the optically thin jet from infrared-X-ray suggests that the synchrotron jet could account for up to
10%
of the observed X-ray emission, although we note that one synchrotron
jet model is able to reproduce the observed hard X-ray power law if the
break in the jet spectrum is at a higher frequency than the OIR SED
data suggest.
The unusual nature of this double-peaked outburst has implications for the study of AMXPs in particular, and neutron star X-ray transients in general. Not only is this source's double-peaked outburst intriguing, it is difficult to understand what physical processes were different to cause the previous outburst in 2004 to follow a more ``traditional'' single-peaked morphology. Observing such systems can place tighter constraints on the range of duty cycles of neutron star XRB outbursts. For short period systems such as AMXPs, it allows us to study the evolution of both the outer accretion disc and the inner synchrotron jet, providing data for more accurate modelling of disc-jet coupling in neutron star XRBs. This study illustrates the importance of rapid, regular multi-wavelength monitoring of AMXP outbursts, and of regular optical monitoring of X-ray binaries in quiescence. As was the case with the second outburst peak of IGR J00291+5934, new outbursts can often be detected by such monitoring with 2-metre class telescopes such as the Faulkes Telescopes and YALO/SMARTS (e.g. Jain et al. 2005; Maitra & Bailyn 2007) before X-ray detection.
AcknowledgementsF.L. would like to acknowledge support from the Dill Faulkes Educational Trust. D.M.R. and A.P. acknowledge support from a Netherlands Organization for Scientific Research (NWO) VENI Fellowship. The Faulkes Telescope Project is an educational and research arm of the Las Cumbres Observatory Global Telescope (LCOGT). RXTE/ASM results are provided by the RXTE/ASM teams at MIT and at the RXTE SOF and GOF at NASA's GSFC. The Westerbork Synthesis Radio Telescope is operated by the ASTRON (Netherlands Institute for Radio Astronomy) with support from the Netherlands Foundation for Scientific Research (NWO). P.G.J. acknowledges support from a VIDI grant from the Netherlands Organisation for Scientific Research. D.S. acknowledges a STFC Advanced Fellowship. The Peters Automated Infrared Imaging Telescope (PAIRITEL) is operated by the Smithsonian Astrophysical Observatory (SAO) and was made possible by a grant from the Harvard University Milton Fund, the camera loan from the University of Virginia, and the continued support of the SAO and UC Berkeley. We thank the staff and students of Paulet High School (Burton-on-Trent, England), The Kingsley School (Leamington Spa, England), Czacki High School (Warsaw, Poland), St. Brigid's School (Denbigh, Wales) and St David's Catholic College (Cardiff, Wales) for contributing to the Faulkes LMXB Observing Program and Schools' Initiative (FLOPSI) and Alison Tripp for scheduling these observations. We thank the anonymous referee for their comments and swift reply.
References
- Alpar, M. A., Cheng, A. F., Ruderman, M. A., & Shaham, J. 1982, Nature, 300, 728 [NASA ADS] [CrossRef] [Google Scholar]
- Archibald, A. M., Stairs, I. H., Ransom, S. M., et al. 2009, Science, 324, 1411 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Altamirano, D., Patruno, A., Heinke, C. O., et al. 2010, ApJ, 712, 58 [Google Scholar]
- Bikmaev, I., Suleimanov, V., Galeev, A., et al. 2005, ATel, 395 [Google Scholar]
- Blake, C. H., Bloom, J. S., & Starr, D. L. 2005, Nature, 435, 181 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Bloom, J. S., Starr, D. L., Blake, C. H., Skrutskie, M. F., & Falco, E. E. 2006, in Astronomical Data Analysis Software and Systems XV, ed. Gabriel, C., Arviset, C., Ponz, D., & Solano, E., ASP Conf. Ser., 351, 751 [Google Scholar]
- Blundell, K. M., Bowler M. G., & Schmidtobreick, L. 2008, ApJ, 678, L47 [Google Scholar]
- Bozzo, E., Ferrigno, C., Falanga, M., et al. 2010, A&A, 509, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Brocksopp, C., Jonker, P. G., Fender, R. P., et al. 2001, MNRAS, 323, 517 [NASA ADS] [CrossRef] [Google Scholar]
- Brocksopp, C., Jonker, P. G., Maitra, D., et al. 2010, MNRAS, 404, 908 [NASA ADS] [CrossRef] [Google Scholar]
- Burderi, L., Di Salvo, T., Lavagetto, G., et al. 2007, ApJ, 657, 961 [NASA ADS] [CrossRef] [Google Scholar]
- Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Chakrabarty, D., Swank, J. H., Markwardt, C. B., & Smith, E. 2008, ATel, 1660 [Google Scholar]
- Chen, W., Shrader, C. R., & Livio, M. 1997, ApJ, 491, 312 [NASA ADS] [CrossRef] [Google Scholar]
- D'Avanzo, P. D., Campana, S., Covino, S., et al. 2007, A&A, 472, 881 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- D'Avanzo, P. D., Campana, S., Casares, J., et al. 2009, A&A, 508, 297 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Eckert, D., Walter, R., Kretschmar, P., et al. 2004, ATel, 352 [Google Scholar]
- Elebert, P., Callanan, P. J., Fillippenko, A. V., et al. 2008, MNRAS, 383, 1581 [NASA ADS] [CrossRef] [Google Scholar]
- Elebert, P., Reynolds, M. T., Callanan, P. J., et al. 2009, MNRAS, 395, 884 [NASA ADS] [CrossRef] [Google Scholar]
- Fender, R. P., De Bruyn, G., Pooley, G., & Stappers, B. 2004a, ATel, 361 [Google Scholar]
- Fender, R. P., Belloni, T. M., & Gallo, E. 2004b, MNRAS, 355, 1105 [NASA ADS] [CrossRef] [Google Scholar]
- Fender, R. P., Russell, D. M., Knigge, C., et al. 2009, MNRAS, 393, 1608 [NASA ADS] [CrossRef] [Google Scholar]
- Filippenko, A. V., Foley, R. J., & Callanan, P. J. 2004, ATel, 366 [Google Scholar]
- Fox, D. B., & Kulkarni, S. R. 2004, ATel, 354 [Google Scholar]
- Frank, J., King, A., & Raine, D. 2002, Accretion Power in Astrophysics: Third Edition (Cambridge, England: Cambridge University Press) [Google Scholar]
- Gallo, E., Fender, R. P., & Pooley, G. G. 2003, MNRAS, 344, 60 [NASA ADS] [CrossRef] [Google Scholar]
- Galloway, D. K., Markwardt, C. B., Morgan, E. H., Chakrabarty, D., & Strohmayer, T. E. 2005, ApJ, 622, L45 [NASA ADS] [CrossRef] [Google Scholar]
- Galloway, D. K., Hartman, J. M., Chakrabarty, D. M., Morgan, E. H., Swank, J. H. 2008, ATel, 1786 [Google Scholar]
- Gehrels, N. 1986, ApJ, 303, 336 [NASA ADS] [CrossRef] [Google Scholar]
- Giles, A. B., Greenhill, J. G., Hill, K. M., & Sanders, E. 2005, MNRAS, 361, 1180 [NASA ADS] [CrossRef] [Google Scholar]
- Heinke, C. O., Altamirano, D., Cohn, H. N., et al. 2010, ApJ, 714, 894 [Google Scholar]
- Hynes, R. I., Haswell, C. A., Chaty, S., Shrader, C. R., & Cui, W. 2002, MNRAS, 331, 169 [Google Scholar]
- Hynes, R. I. 2005, ApJ, 623, 1026 [NASA ADS] [CrossRef] [Google Scholar]
- Jain, R. K., Bailyn, C. D., Orosz, J. A., McClintock, J. E., & Remillard, R. A. 2005, ApJ, 554, 181 [Google Scholar]
- Jonker, P. G., Nelemans, G., Wang, Z., et al. 2003, MNRAS, 344, 201 [NASA ADS] [CrossRef] [Google Scholar]
- Jonker, P. G., Torres, M. A. P., & Steeghs, D. 2008, ApJ, 680, 615 [NASA ADS] [CrossRef] [Google Scholar]
- Jonker, P. G., Miller-Jones, J. C. A., & Homan, J. 2010, MNRAS, 401, 1255 [NASA ADS] [CrossRef] [Google Scholar]
- Krauss, M. I., Wang, Z., Dullighan, A., et al. 2005, ApJ, 627, 910 [NASA ADS] [CrossRef] [Google Scholar]
- Landolt, A. U. 1992, AJ, 104, 372 [NASA ADS] [CrossRef] [Google Scholar]
- Lasota, J.-P. 2001, New AR, 45, 449 [Google Scholar]
- Lewis, F., Russell, D. M., Fender., R. P., et al. 2008a, in Proceedings of the VII Microquasar Workshop: Microquasars and Beyond, 1-5 September, Foca, Izmir, Turkey, Proceedings of Science [arXiv:0811:2336] [Google Scholar]
- Lewis, F., Linares, M., Russell, D. M., Wijnands, R., & Roche, P. 2008b, ATel, 1726 [Google Scholar]
- Linares, M., van der Klis, M., & Wijnands, R. 2007, ApJ, 660, 595 [NASA ADS] [CrossRef] [Google Scholar]
- Linares, M., Tudose, V., & Migliari, S. 2008, ATel, 1667 [Google Scholar]
- Maitra D., & Bailyn, C. B. 2007, in Bursts, Pulses and Flickering: wide-field monitoring of the dynamic radio sky, 12-15 June, Kerastari, Tripolis, Greece, Proceedings of Science [Google Scholar]
- Maitra D., Markoff, S., Brocksopp, C., et al. 2009, MNRAS, 398, 1638 [NASA ADS] [CrossRef] [Google Scholar]
- Markoff, S., Nowak, M. A., & Wilms, J. 2005, ApJ, 635, 1203 [NASA ADS] [CrossRef] [Google Scholar]
- Markwardt, C. B., & Swank, J. H. 2008, ATel, 1664 [Google Scholar]
- Markwardt, C. B., Swank, J. H., & Strohmayer, T. E. 2004a, ATel, 353 [Google Scholar]
- Markwardt, C. B., Galloway, D. K., Chakrabarty, D., Morgan, E. H., & Strohmayer, T. E. 2004b, ATel, 360 [Google Scholar]
- Markwardt, C. B., Altamirano, D., Swank, J. H., et al. 2009, ATel, 2197 [Google Scholar]
- Marshall, F. E., Markwardt, C. B., & Chakrabarty, D. 2008, ATel, 1668 [Google Scholar]
- Migliari, S., & Fender, R. P. 2006, MNRAS, 366, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Migliari, S., Tomsick, J. A., Maccarone, T. J., et al. 2006, ApJ, 643, L41 [NASA ADS] [CrossRef] [Google Scholar]
- Migliari, S., Tomsick, J. A., Miller-Jones, J. C. A., et al. 2010, ApJ, 710, 117 [NASA ADS] [CrossRef] [Google Scholar]
- Oke, J. B. 1990, AJ, 99, 1621 [NASA ADS] [CrossRef] [Google Scholar]
- Patruno, A., Watts, A., Klein-Wolt, M., Wijnands, R., & van der Klis, M. 2009, ApJ, 707, 1296 [NASA ADS] [CrossRef] [Google Scholar]
- Pe'er, A., & Casella, P. 2009, ApJ, 699, 1919 [NASA ADS] [CrossRef] [Google Scholar]
- Pooley, G. 2004, ATel, 355 [Google Scholar]
- Radhakrishnan, V., & Srinivasan, G. 1982, Curr. Sci., 51, 1096 [NASA ADS] [Google Scholar]
- Remillard, R. 2004, ATel, 357 [Google Scholar]
- Reynolds, M. T., Elebert, P., Callanan, P. J., et al. 2006, IAUS, 230, 80 [NASA ADS] [Google Scholar]
- Roelofs, G., Jonker, P. G., Steeghs, D., Torres, M., & Nelemans, G. 2004, ATel, 356 [Google Scholar]
- Rupen, M. P., Dhawan, V., & Mioduszewski, A. J. 2004, ATel, 364 [Google Scholar]
- Russell, D. M., Fender, R. P., & Jonker, P. G. 2007, MNRAS, 379, 1108 [NASA ADS] [CrossRef] [Google Scholar]
- Russell, D. M., Lewis, F., Linares, M., Roche, P., & Maitra, D. 2008, ATel, 1666 [Google Scholar]
- Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Data Analysis Software and Systems IV, ed. R. A. Shaw, H. E. Payne, J. J. E. Hayes, ASP Conf. Ser., 77, 433 [Google Scholar]
- Shaw, S. E., Mowlavi, N., Rodriguez, J., et al. 2005, A&A, 432, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Steeghs, D. 2003, IAUC, 8155, 2 [Google Scholar]
- Steeghs, D., Blake, C., Bloom, J. S., et al. 2004, ATel, 363 [Google Scholar]
- Torres, M. A. P., Jonker, P. G., Steeghs, D., et al. 2008a, ApJ, 672, 1079 [NASA ADS] [CrossRef] [Google Scholar]
- Torres, M. A. P., Jonker, P. G., Bassa, C., Nelemans, G., & Steeghs, D. 2008b, ATel, 1665 [Google Scholar]
- Truss, M. R. 2005, MNRAS, 356, 1471 [NASA ADS] [CrossRef] [Google Scholar]
- Truss, M. R., Wynn, G. A., Murray, J. R., & King, A. R. 2002, MNRAS, 337, 1329 [NASA ADS] [CrossRef] [Google Scholar]
- Tudose, V., Fender, R. P., Linares, M., et al. 2009, MNRAS, 400, 2111 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, Z., Chakrabarty, D., Roche, P., et al. 2001, ApJ, 563, L61 [NASA ADS] [CrossRef] [Google Scholar]
- Wijnands, R. 2006, in Trends in Pulsar Research, ed. J. A. Lowry, (New York: Nova Science Publishers) [Google Scholar]
- Wijnands, R., & van der Klis, M. 1998, ATel, 17 [Google Scholar]
- Wijnands, R., Altamirano, D., Soleri, P., et al. 2008, in A Decade of Accreting Millisecond X-Ray Pulsars, ed. R. Wijnands, D. Altamirano, P. Soleri, N. Degenaar, N. Rea, P. Casella, A. Patruno, & M. Linares, AIP Conf. Proc. (Berlin: Springer), 1068 [Google Scholar]
Online Material
Table 5: X-ray observations of 00291+59 used in this work.
Footnotes
- ... J00291+5934
- Table 5 is only available in electronic form at http://www.aanda.org
- ... magnitudes
- http://eso.org/sci/facilities/paranal/sciops/Bessell/PG0231+051.ps
- ... recommendations
- http://www.swift.ac.uk/pileup.shtml
- ... manual
- http://heasarc.gsfc.nasa.gov/docs/xte/recipes/pcabackest.html
All Tables
Table 1: Photometric observations (NIR, optical and UV detections) of 00291+59 used in this work.
Table 2: Optical and NIR magnitudes of the two comparison stars in the 00291+59 field used for flux calibration.
Table 3:
Swift UVOT optical/UV magnitudes and 3
upper limits of 00291+59.
Table 4: PAIRITEL NIR magnitudes and upper limits of 00291+59.
Table 6: Swift XRT X-ray fluxes.
Table 7: Variability of 00291+59 and a faint field star.
Table 8: Spectra equivalent widths (Å).
Table 5: X-ray observations of 00291+59 used in this work.
All Figures
![]() |
Figure 1:
FT North and INT |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
FT North finder chart (SDSS- |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Upper: optical, UV and NIR data from the 2008 outburst. Lower: X-ray and radio data from the same period. We also plot earlier FTN |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Upper: comparison of the three fades using R-band data from 2004 (Torres et al. 2008a) and |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Upper: 2008 outburst |
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Red arm spectrum from Keck 1, DIB denotes Diffuse Interstellar Bands. For clarity, we mark the wavelengths of He I
lines at 6678 Å and 7065 Å (the latter being a marginal
detection in this spectrum). Insert: zoom in on double-peaked H |
Open with DEXTER | |
In the text |
![]() |
Figure 7: Top row: UV, optical and NIR fluxes for the 2008 peaks. Bottom row: single black-body plus jet contribution model plotted against fluxes for data at each peak. |
Open with DEXTER | |
In the text |
![]() |
Figure 8: Left panel: model for two black-bodies at second outburst peak. Right panel: model showing the dependence on flux for different temperature disc. |
Open with DEXTER | |
In the text |
![]() |
Figure 9: Evolution of the SED of 00291+59. Left panel: broadband SEDs from the 2004 outburst and both 2008 outburst peaks. Right panel: blackbody plus jet contribution shown to approximate data from left panel. |
Open with DEXTER | |
In the text |
![]() |
Figure 10:
Radio flux density against quasi-simultaneous (within one day) X-ray
flux for 2004 and 2008 outbursts. The filled circles represent the
2004 detections (Pooley 2004; Rupen et al. 2004; Fender et al. 2004a);
the cross represents the 2008 non-detection (this paper). The solid
line (without ticks) shows the power-law fit to all three
2004 detections. For comparison, slopes of
|
Open with DEXTER | |
In the text |
Copyright ESO 2010
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.