A&A 451, 187-194 (2006)
DOI: 10.1051/0004-6361:20054235
N. Mowlavi1,2 - I. Kreykenbohm1,3 - S. E. Shaw1,4 - K. Pottschmidt5 - J. Wilms6 - J. Rodriguez7,1 - N. Produit1 - S. Soldi1,2 - S. Larsson8 - P. Dubath1,2
1 - INTEGRAL Science Data Centre, Ch. d'Écogia, 1290 Versoix, Switzerland
2 -
Observatoire de Genève, Ch. des Maillettes, 1290 Versoix, Switzerland
3 -
Institut für Astronomie und Astrophysik - Astronomie, Sand 1, 72076 Tübingen, Germany
4 -
School of Physics and Astronomy, University of Southampton, SO171BJ, UK
5 -
Center for Astrophysics and Space Sciences, University of California at San Diego, La Jolla, CA 92093-0424, USA
6 -
Department of Physics, University of Warwick, Coventry, CV4 7AL, UK
7 -
CEA Saclay, DSM/DAPNIA/SAp(CNRS UMR 7158 AIM), 91191 Gif-Sur-Yvette,
France
8 -
Stockholm Observatory, Alba Nova, 106 91 Stockholm, Sweden
Received 21 September 2005 / Accepted 14 December 2005
Abstract
The decline of the high mass X-ray transient V 0332+53 during the Dec. 2004 to Feb. 2005 outburst
is analysed from the data recorded by INTEGRAL.
The flux is shown to decrease exponentially until 2005 Feb. 10, with a decay time scale
of 30 days above 20 keV and
20 days at lower energies, and to decrease linearly thereafter.
The energy spectrum is well modelled throughout the decay by a power law with a
folding energy of 7.5 keV, and with two cyclotron absorption features.
The folding energy does not vary significantly over the decay, but the spectrum becomes harder
with time.
Most importantly, we show that the parameters describing the fundamental cyclotron line around 27 keV do
vary with time:
its energy and depth increase (by about 17% for the energy in
6 weeks), while its width decreases.
These changes of the cyclotron line parameters are interpreted as resulting from a change in the extent of the
cyclotron scattering region.
Two quasi-periodic oscillations are also observed at various times during the observations,
one at 0.05 Hz and another one near the pulsation frequency around 0.23 Hz.
Key words: X-ray: stars - stars: pulsars: individual: V0332+53 - stars: magnetic fields - X-ray: bursts
The recurrent X-ray transient V 0332+53 experienced an X-ray outburst from December 2004 to February 2005. The outburst was predicted one year earlier from an optical brightening of the massive O8-9Ve companion star, BQ Cam (Goranskij & Barsukova 2004) and was first detected by the Rossi X-ray Timing Explorer (RXTE) in November 2005 (Swank et al. 2004). The INTernational Gamma-Ray Astrophysics Laboratory (INTEGRAL) started to observe the source on 2005 Jan. 6 (Kreykenbohm et al. 2005, hereafter Paper I) and monitored its X-ray activity until the end of February.
The observation of this fourth detected outburst of V 0332+53 in thirty years
revealed the presence of three cyclotron lines
(Coburn et al. 2005; Paper I; Pottschmidt et al. 2005).
The fundamental has an energy of 28 keV, and the first and second harmonics are at 50 keV and 71 keV, respectively.
The magnetic field strength is
G (Paper I, Pottschmidt et al. 2005).
The spectrum is otherwise described by a power law with a high-energy cut-off, typical of the spectra of
X-ray pulsars.
A new Quasi Periodic Oscillation (QPO) around the 0.23 Hz spin frequency of the neutron star has also been
reported in the RXTE data (Qu et al. 2005), in addition to the 0.05 Hz QPO discovered by Takeshima et al. (1994) in the Ginga data of the 1989 outburst.
In Paper I, we presented the INTEGRAL observation of V 0332+53 around the outburst peak. This paper presents the evolution of V 0332+53 during its decline phase. The observations and data analysis method are described in Sects. 2 and 3, respectively. We then present in Sect. 4 the evolution during the decline of the fluxes in several energy bands, of the energy spectrum, and of the variability spectrum of V 0332+53. The results are then discussed in Sect. 5.
Table 1: Summary of observations. From first to last column: (1) Revolution and pointing numbers; (2) Earth date of the start of the observation in UTC / IJD (INTEGRAL Julian Day = MJD - 51544) / MJD (Modified Julian Day); (3) Earth date of the end of the observation in MJD; (4) Individual pointing duration; (5) Total IBIS observation time; (6) Sum of IBIS Good Time Intervals; (7) Total JEM-X1 observation time; (8) Sum of JEM-X1 Good Time Intervals.
The INTEGRAL satellite carries three high energy instruments, the imager IBIS (operating in the range 15 keV to 10 MeV), the spectrometer SPI (20 keV to 8 MeV) and two identical X-ray instruments JEM-X1 and JEM-X2 (3 to 35 keV) operational in turn (see Winkler et al. 2003 for more information on the INTEGRAL mission and instruments). IBIS is further composed of two modules, ISGRI operational down to 15 keV and PICsIT above 200 keV.
INTEGRAL observed V 0332+53 from 2005 Jan. 6 to Feb. 21.
Three main observation sets
were scheduled, in revolutions 272-274 (120 ks), 278 (
80 ks) and 284 (
150 ks),
followed by four
15 ks observations from revolution 285 to 288 to monitor the tail of the outburst decay.
The observations are summarized in Table 1.
Except for the observation in revolution 272, which
was performed in a staring mode, all observations were made following
a hexagonal dithering pattern.
The hexagonal dithering pattern consists of one pointing (stable
spacecraft attitude) with
V 0332+53 positioned at the center of the field of view of the
instruments, and six
pointings with V 0332+53 off-axis in steps of 2 degrees in a
hexagonal pattern.
The integration time for each pointing may vary from one observation
to the next and is given in
Table 1.
The V 0332+53 observations were affected by solar activity which impacted on the quality of the data. In particular, a strong solar flare lasting from Jan. 20 to 23 lead to a split of the IBIS data of revolution 278 into 7760 small Good Time Intervals (GTIs) of only few seconds each. The JEM-X GTIs were greater than 97% of the observation time for that revolution, and in general greater than 90% for all revolutions (see Table 1).
We analyse all the public data collected by the ISGRI and the JEM-X1 instruments on board of INTEGRAL. PICsIT is not useful for our analysis because the spectrum of V0332+53 falls off above 100 keV, and SPI data are not used either since its sensitivity in the 20 to 100 keV range is less favorable than ISGRI (Paper I). Pointings close to perigee passage times are further discarded from the analysis as they are contaminated by the Earth radiation belt.
The ISGRI and JEM-X data are processed with version 5.0 of the Offline Scientific Analysis (OSA) package delivered by the INTEGRAL Science Data Center (http://isdc.unige.ch) using default settings. Images and lightcurves are produced in the 3-5 keV, 5-10 keV and 10-15 keV energy bands for JEM-X, and in the 20-30 keV, 30-40 keV and 40-60 keV bands for ISGRI.
Contamination from other sources is not a problem as
V 0332+53 is the only bright X-ray source in the ISGRI field of view.
Systematic errors are however known to affect the data analysis.
In order to estimate the amplitude of those systematic errors,
and to calibrate the results, an analysis of the Crab was performed
with OSA parameters identical to those used for the analysis of V 0332+53.
Crab public data are taken from revolution 239, pointings 3 to 27
(Modified Julian Day
to 53276.1326),
during which the Crab is located within 3 degrees from the center of
the field of view.
The resulting mean Crab fluxes extracted from the images in the 3-5,
5-10, 10-15 keV bands of JEM-X are 43.5 cps (14% error), 43.2 cps (10%) and 16.4 cps (12%), respectively.
The quoted errors represent the maximum deviations measured from the
mean fluxes over the 25 pointings.
The Crab fluxes in the 20-30, 30-40 and 40-60 keV energy bands of ISGRI are 78.8 cps (13%), 39.4 cps (7%)
and 45.9 cps (4%), respectively.
![]() |
Figure 1:
a) Pointing averaged ISGRI flux evolution of V 0332+53
as a function of time in the different energy bands as
labelled in the figure.
The fluxes have been normalised to Crab units.
The error bars displayed in the figure represent only
the statistical errors; neither do they take into account the
uncertainties in the Crab fluxes in the respective energy bands.
The solid curves are exponential (from
![]() ![]() ![]() ![]() |
Open with DEXTER |
The fluxes of V 0332+53 in the ISGRI and JEM-X energy bands are plotted in Fig. 1 as a function of time. They are normalised to the Crab fluxes in the respective bands.
The decline phase can be described by an exponential decay up to
,
followed by a linear decrease. The exponential decay is less
rapid at the ISGRI energies than at the JEM-X energies.
From a simple estimation, the exponential decay time is about 30 days above 20 keV (solid curves in Fig. 1a)
and about 20 days below 15 keV (solid curves in Fig. 1b).
The exponential decay is also confirmed by the 2-12 keV lightcurve
recorded by the All Sky Monitor (ASM) on board of RXTE,
plotted in Fig. 1b for comparison with the JEM-X
lightcurves
(see Levine et al. 1996 for more information on the
RXTE satellite).
The evolution of the ISGRI and JEM-X hardness ratios are shown in
Fig. 2.
Above 20 keV the spectrum of V 0332+53 softens with time, with the
[30-40] keV / [20-30] keV
ratio decreasing by about 36% from an average value of 0.55 in
the first observation to
0.35 in the last observations.
On the other hand, below 15 keV, the spectrum becomes harder, with
the [10-15] keV / [5-10] keV
ratio increasing by about 28% from an average value of
1.8 to
2.3 during the decline.
As will be seen in Sect. 4.2, the opposite behaviour
of the hardness ratios above 20 keV
and below 15 keV are understood from the evolution of the spectral
shape of V 0332+53 during the decline.
The ISGRI/JEM-X hardness ratio ([20-30] keV/ [5-15] keV) also
evolves with time, increasing from an average of 0.7 to
1.5 during the decline.
This results from the increase of the spectral power law index
presented in the next section.
![]() |
Figure 2: Hardness ratios computed from the pointing-averaged (and Crab normalized) fluxes displayed in Fig. 1. Upper triangles: ISGRI [30-40 keV]/[20-30 keV]; lower triangles: JEM-X [10-15 keV]/ [5-10 keV]; squares: ISGRI [20-30 keV] / JEM-X [5-10 keV]. |
Open with DEXTER |
Table 2: Fitting parameters. From first to last column: (1) revolution number (see Table 1); (2) normalisation factor; (3) powerlaw index; (4) folding energy (cutoffpl in XSPEC); (5) energy of the fundamental cyclotron line; (6) width of the line; (7) optical depth at the line center; (8), (9) and (10) same as 5th, 6th and 7th columns, respectively, but for the first harmonic cyclotron line; (11) flux in the 2-10 keV band; (12) flux in the 20-40 keV band.
The JEM-X and ISGRI spectra of V 0332+53, normalised to the Crab
spectrum,
are shown in Fig. 3 at four different times during the decline of the source.
The maximum flux relative to the Crab is recorded between 10 and 20 keV,
at an energy
which increases with time.
Fitting Gaussian curves to the JEM-X spectra in a 10 keV width range around the maximum leads
to
,
15.4, 16.8 and 18.0 keV for revolutions 273, 278, 284 and 287, respectively.
This shift of
explains the hardening (softening) at low (high) energies noticed
from the JEM-X (ISGRI) hardness ratios presented in Sect. 4.1.
The spectra of accreting X-ray pulsars include contributions from at least
an accretion column and from the hotspots at the surface of the neutron star.
In the absence of a consistent model, various empirical models are commonly used
in the literature to fit such spectra (see for example Kreykenbohm et al. 1999 for a discussion).
In this study, we use the simple "cutoffpl'' model of the form
![]() |
Figure 3: JEM-X and ISGRI spectra integrated, from upper to lower curves, during revolutions 273, 278, 284 and 287, respectively. The spectra have been normalised to the Crab spectrum of revolution 239 (see text for details). |
Open with DEXTER |
To fit the cyclotron lines we use Gaussian optical depth profiles of the form
![]() |
Figure 4: Evolution of the specral fit parameters describing the fundamental cyclotron line. a) Line energy, b) width and c) depth. The horizontal bars on the points indicate the duration of each observation. |
Open with DEXTER |
The main fit parameters are summarised in Table 2.
The 2-10 keV and 20-40 keV fluxes derived from the fitted models confirm
the two phases set forth in Sect. 4.1, i.e. first an
exponential decay followed by a linear decrease of the flux.
The folding energy remains about constant around 7.5 keV, which is
consistent with the results
of Paper I (Table 1) for revolution 273-274.
The power law index ,
on the other hand, decreases from -0.18
in the early revolutions to -0.4 in the later ones, in agreement
with the hardening observed
from the hardness ratios (Fig. 2).
![]() |
Figure 5:
a) Combined spectrum and folded model of data obtained with JEM-X ( left) and ISGRI ( right) for revolution 284;
b) residuals of the model with
![]() ![]() ![]() |
Open with DEXTER |
The cyclotron line parameters for the fundamental are displayed in
Fig. 4 as a function of time, and an example of the
unfolded spectrum with the fit and residuals for revolution 284 is
shown in Fig. 5.
The energy of the fundamental cyclotron line increases by 7% from
the first to the last INTEGRAL observations, while the depth
of the line increases by 16%
and the width decreases by 12%. The changes are significant;
attempting to fit the spectrum
in revolution 284 with a line at 27.5 keV leads to strong residuals
around the line (Fig. 5e). The changes in the line
parameters cannot be attributed to variations of the continuum
either; fixing
in revolution 284 to its value of -0.17 from
revolution 273 leads to strong residuals at energies above 45 keV
(Fig. 5f).
The determination of the parameters describing the 50 keV line
becomes problematic from revolution 284 on, due to poor statistics.
Therefore, for those revolutions, we fix the line energy to the value
keV obtained for Rev. 273 where the statistics
are the best, and fit the width and depth of the line to obtain the
best residals. The results are shown in Table 2. The
width and depth of line are seen to be rather stable throughout the
whole decay of the outburst, around 9 keV and 2, respectively. The
large error bars on the numbers result from the poor statistics.
![]() |
Figure 6: 5-20 keV power spectra, averaged over the indicated revolutions and normalised using the method described by Miyamoto et al. (1991). Poisson noise has not been subtracted (see text for more details). The power spectra of revolutions 278 and 273 have been multiplied by factors of 5 and 25, respectively. Error bars have been omitted for clarity. |
Open with DEXTER |
Several X-ray pulsars and cyclotron line sources are known to exhibit
low frequency quasi-periodic oscillations (QPOs) (see, e.g.,
Takeshima et al. 1991; Takeshima et al. 1994;
Heindl et al. 1999). Takeshima et al. (1994)
first discovered a QPO with a centroid frequency of 0.05 Hz and
a relative root mean square (rms) amplitude of
5% in
2.3-37.2 keV data from the 1989 outburst of V 0332+53. This
feature has recently also been reported in the 2-60 keV
RXTE-PCA power spectral density functions (PSDs) obtained
during the 2004 outburst decay (Qu et al. 2005). Furthermore
Qu et al. found a QPO at
0.22 Hz, i.e. centered on
the fundamental pulsar frequency, in addition to the narrow peak at
that frequency. The rms strength of this QPO also amounts to a few per
cent. In the following we show that the JEM-X and ISGRI observations confirm the detection of both features during the
outburst decay, and we analyse their changes with time.
Note that, to our knowledge, these results represent the first QPOs
reported from INTEGRAL JEM-X and ISGRI data in the literature.
The background subtracted 5-20 keV and 20-40 keV lightcurves of each revolution were rebinned to a time resolution of 0.5 s and power spectra were produced from 512 s long lightcurve segments using Fast Fourier Transforms. Averaged power spectra for selected revolutions are shown in Figs. 6 and 7.
For coded mask instruments the determination of the white noise level
as well as the influence of the background flux on the power spectrum
can be complex. The white noise component in JEM-X lightcurves,
e.g., is known to be non-Poissonian.
This is also true for the ISGRI which consists of 16 000 pixels, each
having different noise charactersistics.
Using the normalization described
by Leahy et al. (1983) we find the noise level to be
about one order of
magnitude higher than expected from Poisson noise. However, source
features are well apparent above the noise level, e.g., the 0.05 Hz
QPO is clearly visible in revolution 284. Therefore, and since the
broad-band PSD (which is not to be discussed here) is dominated by the
noise component up to at least 0.1 Hz, we did not attempt to correct
the power spectra for white noise. The JEM-X power spectra
shown in Fig. 6 are normalised according to
Miyamoto et al. (1991) and a comparison of the
strength of features
above the noise level between the JEM-X PSDs and the
quasi-simultaneous PCA PSDs of Qu et al. (2005) finds
comparable values.
![]() |
Figure 7: 20-40 keV power spectra, averaged over the indicated revolutions. The spectra are arbitrarily shifted in the y-axis for easy comparison. Note that the ISGRI power spectra of revolution 278 are not shown due to solar flaring giving rise to many gaps in the lightcurve after GTI filtering. |
Open with DEXTER |
Strong changes in the shape of the power spectrum can be seen over the outburst decay, where the varying amplitude ratios at the pulsar frequency and its harmonics reflect a strong change of the pulse profile (Zhang et al. 2005). The PSDs shown in Figs. 6 and 7 represent the main stages of the evolution. Although, for the reasons discussed earlier, the white noise level is not accurately estimated in the INTEGRAL data, the JEM-X and ISGRI PSDs for revolutions 0273 and 0284 do compare well and show similar features as a function of frequency.
In revolution 273 (Jan. 8-9) the power spectrum is dominated by
signatures of the pulsation frequency, especially a very strong first
harmonic. In addition weak red noise is present. The peak at the
fundamental pulsation frequency is slightly broadened at the bottom,
especially in the ISGRI PSD. However, the 0.22 Hz QPO is not as
obvious as a few days later in the Jan. 15-19 PCA power spectrum of
Qu et al. (2005). In revolution 278 (Jan. 23-24) the fundamental
pulsation at 0.22 Hz is not seen anymore but a weak broad QPO feature
is present at this frequency. The first harmonic is still clearly
visible. The strength of the red noise component more than doubled at
low frequencies and there is an indication of additional power in the
frequency region of the 0.05 Hz QPO. Principally the overall PSD shape is very similar to the Jan. 15-19 PCA power spectrum mentioned
above. Finally, in revolution 284 (Feb. 9-10) the 0.05 Hz QPO is
clearly visible, especially in the JEM-X power spectrum. The
red noise is still strong while the pulsation fundamental is back as a
narrow feature, stronger than the first harmonic. Fitting the low
frequency QPO in revolution 284 with a Lorentzian gives (
) Hz, with a FWHM of (
) Hz for JEM-X
and (
) Hz with a FWHM of (
) Hz for
ISGRI, where the underlying continuum can be modelled with a flat power
law component and a very low frequency Lorentzian.
V 0332+53 is a highly magnetised pulsar which has a magnetic field
of about
G (Paper I, Pottschmidt et al. 2005) and orbits a massive O8-9Ve companion with a
period of about 34 days.
The X-ray outbursts from this system are observed with a time recurrence of about 10 years.
They are believed to be triggered by a sudden increase of the mass ejection rate from the massive star,
the signature of which has been observed in the optical brightening
of the system in 2004 (Goranskij & Barsukova 2004).
The ejected mass feeds an accretion disk around the pulsar, from
which it is funneled by the magnetic field towards the polar caps and
falls on the surface of the neutron star.
The analysis of the 2005 outburst presented in Sect. 4 provides new hints towards a better understanding of V 0332+53 within this basic picture. The X-ray lightcurve (Sect. 4.1) displays a characteristic exponential decay with a time scale of 20 to 30 days, followed by a linear decrease. This feature may provide some information on the disk properties. The X-ray spectrum (Sect. 4.2) is well fitted by a folding power law modified by two cyclotron lines. The 27 keV energy of the fundamental cyclotron line is shown to increase by 7% during the outburst decay and to become narrower by 12%. At the same time, its depth increases by 16%. This information would provide some insight into the resonance scattering region near the polar caps. Finally, the appearance of the 0.05 Hz QPO only from the middle of the decay on (Sect. 4.3) provides further information on the properties of the matter around the neutron star. Each of those observational facts are successively discussed in the following sections.
An exponential decay of the flux is observed in dwarf novae and in soft X-ray transients on timescales from a few days for the former to several weeks for the latter (e.g. Lasota 1996; Tanaka & Shibazaki 1996). King & Ritter (1998) show how illuminated disks in low-mass X-ray binary stars can produce such an exponential decay on the time scales observed for the soft X-ray transients. In that picture, the central star illuminates the disk and keeps its temperature high. The disk empties proportionally to its mass, leading to an exponential decay of the emitted flux. They further show that the lightcurve would be characterised by a linear decrease at later times. The change to a linear regime is triggered by a change in the conditions in the disk. In low-mass X-ray binaries the change corresponds to the decrease of the disk temperature below the temperature of partial H ionisation, such that the high temperature cannot hold over the entire disk anymore. This scenario of X-ray irradiated accretion disks has subsequently been confirmed by numerical calculations (Dubus et al. 1999, 2001).
It is tempting to make the parallel with the results found for V 0332+53. The source displays an exponentially decreasing flux followed by a linear decrease, similar to the disk-irradiated systems. This is the first time, to our knowledge, that such a behaviour is observed in a high-mass X-ray binary. Since the X-ray flux is proportional to the mass accretion rate, this picture suggests that the accretion rate is proportional to the mass of the disk during the exponential decay phase. The transition to a linear decrease phase would be triggered by an as yet un-identified change of conditions in the disk.
The two different exponential decay times observed in the lightcurves of V 0332+53, of 30 days above 20 keV and 20 days below (Sect. 4.1) suggest the presence of at least two regions contributing to the X-ray continuum of that source. The harder component would originate from the accretion column near the resonant scattering region, while the softer component would originate from a region located, for example, higher in the accretion column. The observed increase of the depth of the cyclotron line with time is also compatible with this scenario; the softer component decays more rapidly and hence its tail would contribute relatively less to the high energy part of the spectrum from which the photons are scattered by the electrons in the magnetic field. Basically, one time scale governs the rate of mass flow onto the neutron star through the disk, leading to the exponential decay followed by a linear decrease (see above); another time scale governs the spectral changes of the emission that may come from different heights within the accretion column.
Time evolution of cyclotron resonance line parameters has already
been reported in the literature for some binary systems.
Mihara et al. (1998) report a relation between the
line energy and the X-ray luminosity for 5 pulsars observed by Ginga.
The resonance energy of 4U 0115+53, for example, increased by 40% as
the luminosity decreased by a factor of 6 between 1990 and 1991. The
1990 outburst of V 0332+53 is also reported in that paper, with a 10% decrease of the line energy from
keV to
keV as the Ginga count rate in the 3-37 keV band increased
by a factor of 1.6.
These authors attribute the change of the line energy to that of the
height of the scattering region in a dipole magnetic field, leading
to a height change of
1.1 km for 4U 0115+53 and
330 m
for V 0332+53.
Kreykenbohm et al. (2004), on the other hand,
report a variation of the cyclotron energy of GX 301-2 with the
pulse phase of the pulsar, and attribute the change to different
viewing angles of the accretion column where the line originates.
In the case of the 2005 outburst decay of V 0332+53 analysed in this
paper, the increase with time of the energy and depth of the
cyclotron fundamental line, together with the decrease of its width
can also be understood as resulting from the variation of the
scattering region characteristics in a dipole magnetic field above
the polar caps of the pulsar.
It is interesting to note that Pottschmidt et al. (2005), analyzing the RXTE data of V 0332+53
taken two weeks before the first INTEGRAL observation,
derive
E1=25.2 keV for 2004 December 24-26, which confirms the
trend of increasing line energy with time found in this paper.
Let us thus assume that the region where the cyclotron resonance
scattering takes place extends over some fraction of the
accretion column. Then the region would cover some range of magnetic strengths,
and a broad line shape would result from the superposition of a succession of narrow lines,
each originating from a different height and magnetic field strength in the accretion column.
As the accretion rate and hence the flux decreases with time, the extension of the
resonance scattering region is expected to decrease along with the total width of the line.
At the end of the decline only the region close to the neutron star surface, where the magnetic
field is stronger, would contribute to the resonance scattering; this
would lead to a narrower line shape at a higher energy.
The observed energy E1 indeed increases by 17% from the RXTE
observation of December 2004 to the last INTEGRAL observation in February 2005.
If we consider a dipole approximation for the magnetic field, we have
in a non-relativistic description, where r is the
radius to the center of the neutron star.
The extent of the resonance scattering region at the peak of the
outburst would then, in this first scenario, be about 6% of the
stellar radius, i.e. about 500 meters for typical radii of neutron
stars. A change in the scattering energy between two different
observations of Her X-1, from 34 keV to 41 keV, has been discussed
by Gruber et al. (2001); they conclude that the range
over which the scattering occurs is of the order of 1 km in height,
if their observed line width of 5 keV is largely due to a mix of
magnetic field strengths in a dipole field. In the case of V 0332+53,
we additionally observe an evolution of the width and depth of
the line with time.
Alternatively, we can argue that the fundamental line originates from a unique narrow region in the column density, but that the line shape is quite complex and variable with time, resulting in variations of the line width and depth. An increase in the line energy would then be related to a displacement with time of the scattering region. In the case of V 0332+53, the observed 17% increase of E1 would imply a move of the location of the scattering region by 500 m in the accretion column down to the surface of the neutron star. However, in this second scenario, the line profile is expected to remain approximately constant with time, which is not supported by the INTEGRAL observations. Also, if the scattering region gets closer to the neutron star, higher order terms in the B-field may increase the line width, which is opposite to the trend observed by INTEGRAL. We therefore would rather support the first scenario.
The 0.22 Hz QPO observed around the spin frequency in both the INTEGRAL (Sect. 4.3) and RXTE (Qu et al. 2005) data must originate from a region around the neutron star co-rotating with the pulsar. This region can either be located at the inner radius of the accretion disk near the Alfvén radius, or in the corona above the hot spots, or somewhere in between along the magnetic field through which the matter is tunneled.
The origin of the lower frequency QPO at 0.05 Hz, on the other hand,
is more debatable.
Low frequency QPO-type variabilities are known in several X-ray
pulsars (see van der Klis 2004 for a review).
Although there is no general consensus on the origin of QPOs, most models do invoke the
presence of some sort of inhomogeneities in the accretion disk.
With a typical neutron star mass
and a spin of
s for V 0332+53,
the co-rotation radius where the angular velocity of the
magnetosphere and the Keplerian velocity of the disk are equal is
![]() |
(4) |