A&A 372, 557-562 (2001)
DOI: 10.1051/0004-6361:20010522
B. T. Gänsicke1 - A. Fischer1 - R. Silvotti2 - D. de Martino2
1 -
Universitäts-Sternwarte Göttingen,
Geismarlandstr. 11, 37083 Göttingen, Germany
2 -
Osservatorio di Capodimonte, Via Moiariello 16, 80131 Napoli, Italy
Received 28 June 2000 / Accepted 4 April 2001
Abstract
We present a simple quantitative model that can describe
the photometric B and V band light curves of AMHerculis obtained
during a high state. The double-humped shape of the V band light
curve is dominated by cyclotron emission from a region at the main
accreting pole with an area of
and
sustaining an inflow of
0.06
.
The almost unmodulated B band is dominated by emission from the accretion stream. The
contribution of the heated white dwarf to the optical emission is
small in the B band, but comparable to that of the accretion stream
in the V band. The emission of the secondary star is negligible
both in B and in V.
Key words:
accretion - novae,
cataclysmic variables -
stars: binaries close -
stars: individual: AM Her -
stars: white dwarfs
In AMHerculis stars, or polars, a strongly magnetic white dwarf
(
MG) accretes from a Roche-lobe filling late type secondary
star. The flow of matter leaving the secondary star is threaded by the
magnetic field once the magnetic pressure exceeds the ram pressure in
the accretion stream. The kinetic energy is converted into heat in a
shock near the footpoints of the flux tubes on the white dwarf
surface. For low and intermediate local mass flow densities
(
), this shock occurs above the white
dwarf surface, and the cooling thermal bremsstrahlung and cyclotron
radiation can escape unhindered into space. For high mass flow densities
(
), the shock is submerged in the
white dwarf atmosphere, and the primary accretion luminosity is
reprocessed into the soft X-ray regime.
The strong magnetic field of the white dwarf synchronises its rotation
with the binary orbital period (
). As an observational
consequence, the emission of polars is strongly modulated at the
orbital period in almost all wavelength bands.
While the shape of the hard X-ray light curve of AMHerculis (and
most other polars) can be easily interpreted in terms of the
changing geometric aspect of the hot plasma below the shock, the
optical through infrared emission may show a more complicated phase
dependence, as different components within the binary contribute to the
observed emission at a given wavelength, and as the cyclotron emission
from the accretion column is subject to wavelength-dependent beaming.
In this paper we develop a simple quantitative model which takes into account the various emission sites in AMHerculis, and which can quantitatively describe the observed B and V band high state light curves with a minimum of free parameters.
We obtained B and V high state photometry of AMHer in August
1998 at the Loiano 1.5m telescope using a 2-channel photoelectric
photometer and a time resolution of 1s (Table 1). The
observations were interrupted every 60-120 min for sky measurements
in both channels. Data reduction included sky subtraction,
differential extinction, absolute calibration using Landolt standards,
and binning into 5s bins.
The timings have been converted to orbital phases using the ephemeris
of Heise & Verbunt (1988), but defining
as
the inferior conjunction of the secondary star
(Southwell et al. 1995)
. In AMHer, the conversion
between the phase conventions is given by
.
The
high precision of the ephemeris of Heise & Verbunt (1988) ensures
that the error in the computed phases is
,
which is of no concern for the
purpose of the present paper.
A phenomenological description of the photometric data used here has
been included in the analysis of quasi-simultaneous BeppoSAX
observations by Matt et al. (2000).
The observed light curves (Fig. 1) are typical of AMHer in
the high state (e.g. Olson 1977; Szkody & Brownlee 1977). The V band shows a broad and round minimum near
and a
secondary minimum near
.
The B band is almost
unmodulated. Both light curves exhibit a significant amount of
flickering on timescales of several minutes, which is probably due to
inhomogeneities in the accretion flow (Matt et al. 2000).
We identify four possible emission sites that will contribute to the observed light curves: the (heated) white dwarf, the secondary star, the accretion stream feeding material from the secondary star towards the white dwarf, and the accretion column just above the white dwarf surface near the magnetic pole. In the following, we will discuss the system geometry of AMHer and the four individual emission components.
Date | Time (UT) | Exp. time [s] | Filter | Mean magnitude |
Aug. 20 | 20:45:21 | 13320 | V | 13.50 |
Aug. 22 | 19:45:20 | 14040 | B | 13.64 |
![]() |
Figure 1: Calibrated B and V light curves obtained August 20 and 22, respectively, at the Loiano 1.5 m telescope. Indicated are the orbital phase (bottom) and magnetic phase (top). Both observations cover the binary orbit more than once. The two consecutive orbits are shown by black and grey points to highlight the cycle-to-cycle variations. |
Open with DEXTER |
The orientation of the accretion column with respect to the observer
can be described by a set of three angles: the inclination of the
system, i, the co-latitude of the accretion column/spot, ,
measured as
the angle between the rotation axis and the field line threading the
accretion column and the azimuth of the magnetic axis from
the line connecting the centres of the two stars,
.
For our present purpose, the exact knowledge of
is not crucial,
as it defines only a phase offset between our model light curve and
the observed light curve, and we will adjust
accordingly below.
The inclination and the co-latitude of AMHer have been
subject of dispute. Brainerd & Lamb (1985) derived
from the observed polarization properties
and
.
Another polarimetric study
indicates
and
-60
(Wickramasinghe et al. 1991). The higher value of the inclination is
favoured by spectroscopy of the secondary and of the accretion stream
(Davey & Smith 1996; Southwell et al. 1995; Gänsicke et al. 1998).
An additional constraint
comes from the
observed self-eclipse of the X-ray/EUV/ultraviolet emission from the
accretion region and from the heated surrounding white dwarf
photosphere. Gänsicke et al. (1998) and Sirk & Howell (1998)
derive independently
.
We will use the combinations (
)
and
(
)
in our models of the B and Vlight curves of AMHer in order to reflect the uncertainties in the
exact geometry, while satisfying
.
IUE and HST observations show that the ultraviolet flux of
AMHer is dominated by the emission from the accretion-heated white
dwarf during both low states and high states
(Gänsicke et al. 1995). Gänsicke et al. (1998) could
quantitatively model the quasi-sinusoidal flux modulation observed with
HST during a high state with a large, moderately hot spot
covering a white dwarf of
K. They chose a temperature
distribution decreasing linearly in angle from the spot centre merging
into
at the spot radius. The best-fit spot parameters implied a
central temperature of 47000K and a spot size of
of
the total white dwarf surface. With a distance of 90pc, the inferred
radius of the white dwarf was
cm.
We use here synthetic phase-resolved spectra computed with the 3D
white dwarf model described by Gänsicke et al. (1998), using their
best-fit parameters.
We note for completeness that during a low state the U and Blight curves of AMHer can very well be dominated by emission from the heated white dwarf (see e.g. Fig. 1 of Bonnet-Bidaud et al. 2000).
The near-ultraviolet emission (
Å) of AMHer
during the high state is dominated by continuum emission from the
accretion stream. Based on phase-resolved IUE spectroscopy,
Gänsicke et al. (1995) showed that this near-ultraviolet emission
of the accretion stream does not vary noticeably with the orbital
phase. In the following, we assume therefore that the accretion
stream emission does not vary over the binary orbit, and treat the
(constant) stream contributions in B and V as a free
parameters. The limitations implied by this assumption will be
discussed in Sect. 4.
The secondary star in AMHer can be detected during low states, and both its spectral type and magnitude are easily measured (Schmidt et al. 1981; Gänsicke et al. 1995). For completeness, we use here B=18.43 and V=16.83, as derived from the low state data, to describe the contribution of the secondary. It is, however, clear that the secondary star contribution to the observed B and V high state light curves is negligible.
The accreted material releases its kinetic energy in a hydrodynamical shock close to the white dwarf surface. The ions are heated to the shock temperature - a few 107 K - while the electrons retain at first their pre-shock temperature because of their low mass. Coulomb interactions between the particles transfer energy from the ions to the electrons, such that their temperature increases below the shock. The electrons then cool radiatively through emission of cyclotron radiation and thermal bremsstrahlung. This interplay between heating and cooling results in a temperature distribution along the post-shock flow which is characterized by its maximum temperature (Woelk & Beuermann 1996).
![]() |
Figure 2:
Sample cyclotron model spectra from the
accretion column in AMHer for several angles ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3:
Angle between the magnetic field line and the
line-of-sight as function of the orbital phase. Angles ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
Viewing geometry for an inclination
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In order to compute the emission from the accretion column we approximated the volume between the shock and the stellar surface with a given cross section and height. The total emission of this column is computed with a ray-tracing algorithm, where we solve the radiative transfer equation for each ray passing through the structured column using an updated version of the code described by van Teeseling et al. (1999). For the computation of the angle-dependent cyclotron emission we assumed that the magnetic field within the accretion column is perpendicular to the stellar surface and does not vary over the considered height of the column. A full description of the code used to solve the radiative and hydrodynamical problem is given elsewhere (Fischer & Beuermann 2001).
High state X-ray observations of AMHer give bolometric hard X-ray
fluxes of (1-3)
(Gänsicke et al. 1995; Ishida et al. 1997; Matt et al. 2000). Assuming that
the cyclotron emission can contribute roughly the same amount of flux
(Gänsicke et al. 1995), we estimate an accretion rate
-1016
for the stand-off shock region
.
The magnetic field strength of AMHer is
(Bailey et al. 1991). For the accretion rate estimated above and a
dipolar field geometry, a cross section of the accretion column of
is estimated near the
magnetic pole (Lubow & Shu 1975, 1976; Heerlein et al. 1999).
From these numbers, a mean local mass flow density
-0.2
is derived. We have computed a grid of
spectra for several mass flow densities in the given range and for a
number of angles between the line of sight and the magnetic field. A
representative sample of spectra for different
is shown in
Fig. 2.
From the computed set of model spectra, we synthesized B and Vlight curves of the cyclotron emission from the accretion column,
where the mass flow density
was taken as a free parameter.
The fundamental parameter in the computation of the observed cyclotron
flux for a given orbital phase is the angle between the line of sight
and the magnetic field axis,
.
For a given set of
inclination i, co-latitutde
and orbital phase
![]() |
(1) |
Figure 3 shows the variation of
with
for (
)
and for (
), the viewing geometry is illustrated in
Fig. 4. We chose
so that the
minimum in
aligns with the minimum of the observed Vlight curve, i.e. the magnetic pole leads the secondary by
.
This value is in good agreement with the azimuth of the hot
accretion region derived from X-ray and ultraviolet light curves,
(Paerels et al. 1996; Gänsicke et al. 1998).
reaches a minimum near
,
when the
observer looks down along the accretion funnel, passes through
at
and reaches the maximum of
at
.
As the cyclotron beaming is
sensitive only to
we plot
when
.
From Fig. 3 it is clear that a
reduction of the cyclotron flux is expected in the range
-0.6 from cyclotron beaming alone, without taking
into account a possible occultation of the accretion column by the
body of the white dwarf (see also Sect. 4).
A very similar qualitative description of the phase-dependent
cyclotron emission has been given by Bailey et al. (1984) in order
to explain the polarimetric properties of AMHer: linear polarization
pulses and a change of sign of the circular polarization are observed
at
and 0.3, which corresponds to the phases when
.
We obtain the model light curves by summing the four
individual emission components, white dwarf, accretion stream,
secondary star, and accretion column.
To fit the flux and the shape of the observed B and V light
curves, we use two free parameters: the mean local mass flow density
which determines the (phase-dependent) contribution of the
cyclotron emission and the (phase-independent) contribution of the
accretion stream,
and
.
As
mentioned above, the contributions of the heated white dwarf and of
the secondary star are fixed. The observed light curves are reasonably
well reproduced for
,
,
and
(Figs. 5, 6).
Figures 5 and 6 show the model light curves for
the two different choices of (
). The contribution of the
heated white dwarf to the observed B and V high state light curves
is small (Figs. 5, 6), the contribution of the
secondary is negligible.
Both fits grossly reproduce the shape of the observed light curves,
i.e. a broad and deep minimum at
and a shallow secondary
minimum at
in V and practically no orbital modulation
in B. However, for (
)
the angle
drops below
during
-1.1, with the
result that cyclotron radiation contributes very little to the
total V band emission. Consequently, this model light curve has a
flat bottom, which contrasts with the observed round shape.
![]() |
Figure 5:
The observed light curves from
Fig. 1 have been sampled into
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
As Fig. 5, but assuming
![]() ![]() |
Open with DEXTER |
We have made a number of simplifications and assumptions in our light curve model that need some further attention.
(1) We have assumed that the accretion stream emission does not vary over the orbit, and have treated its contribution as a free parameter to fit the observed flux level. The computation of realistic accretion stream spectra is a problem of great complexity (e.g. Ferrario & Wehrse 1999), and certainly beyond the scope of our paper. However, a short inspection of the general spectral features of accretion stream emission lends support to our assumption.
The optical high state spectrum of AMHer is characterized by strong
Balmer and Helium emission lines and by a strong Balmer jump in
emission (e.g. Schachter et al. 1991). The origin of such a
spectrum is clearly a plasma in which the optical depth in the continuum
is below unity (optically thin in the continuum), and the
optical depth in the lines is much larger than in the continuum,
i.e. 1 (optically thick in the lines). The details of
the spectrum depend on the temperature, the density, and the extension
of the plasma region, but the general characteristics are (a)
significant continuum emission in the near-ultraviolet part of the
Balmer continuum, (b) a strong contribution in the B band, which
covers the higher lines of the Balmer series, Ca K+H in
emission, as well as the Balmer jump, and (c) a weak contribution in
the V band, which covers practically only the weak continuum between
and
.
As already stated in Sect. 3.3,
Gänsicke et al. (1995) detected strong near-ultraviolet continuum
emission with a weak phase-dependence, and identified the accretion
stream as the source of this emission. In AMHer, with
MG, cyclotron emission cannot contribute significantly
at such short wavelengths, the secondary star can be excluded, and the
white dwarf contribution was already modelled in the analysis of
Gänsicke et al. (1995). We feel, therefore, it is safe to identify the
observed near-ultraviolet continuum flux with the accretion stream
feature (a) defined above.
We find in our present analysis that the accretion stream contributes more in B than in V. In fact, the stream emission dominates the practically flat B band light curve. This result appears entirely plausible, as our above considerations show that the stream is indeed expected to be brigher in B than in V. We note that weak photometric modulation at blue wavelengths is a hallmark of AMHer during the high state (e.g. Olson 1977). A more sophisticated approach would need simultaneous phase-resolved spectroscopy, so that the stream contribution in the B and V bands could be quantitatively estimated e.g. from the orbital variation of the emission line fluxes. Hence, in summary, the observed weak phase dependence of the near-ultraviolet stream emission and the fact that the stream is expected to be brighter in B than in V strongly supports our assumption of a phase-independent stream contribution. This in turn implies, as shown in Sect. 3.6, that cyclotron emission is the dominant phase-dependent contribution to the observed V band high state light curve.
(2) We ignore the geometric occultation of the accretion column by the
body of the white dwarf. For
structures
extending less than than
0.035
above the white dwarf
surface will be totally eclipsed at
.
However, for
the low magnetic field and the rather low mass flow density, the shock
is expected to stand
0.1
above the white dwarf surface
(Fischer & Beuermann 2001). The X-ray data of AMHer
clearly show that at least part of the accretion region is
eclipsed at
(i.e. when
reaches its
maximum) by the white dwarf (e.g. Hearn & Richardson 1977). In
contrast to this, significant polarized emission is still detected
during this phase, indicating either that the cyclotron emission
arises at a greater height above the white dwarf surface than the
X-rays (Bailey et al. 1984), or that additional cyclotron emission
comes from the opposite magnetic pole (Wickramasinghe et al. 1991).
An indication that an eclipse of the accretion column by the white
dwarf might be relevant comes from the UBVRI polarimetry of
Piirola et al. (1985): the secondary minimum at
has practically the same shape in V, R, and I. This is
unexpected, as the cyclotron beaming decreases towards the lower
cyclotron harmonics (i.e. longer wavelengths). In fact, R and I band
light curves calculated with our model give reasonable magnitudes, but
do not show a strong secondary minimum.
(3) We have neglected that the accretion column has a finite lateral
extent, which will cause a variation of the magnetic field strength
and, more importantly, of its orientation within the cyclotron emitting
region. The expected result of such variation is a smoothing of
,
and, consequently, of the cyclotron light curve.
In addition, a significant lateral extent of the cyclotron
emitting region would counteract to some extent the self-eclipse
just described above.
(4) We deliberately restrict our model to a single cyclotron region near the main (upper) pole. The X-ray emission of AMHer is known to switch between (at least, see below) two different patterns. In the normal mode, hard and soft X-rays are emitted in phase from the upper pole. During the reversed mode, the soft X-ray emission originates predominantly at the lower pole, while the hard X-rays are still emitted from the upper pole (Heise et al. 1985). Despite this apparently drastic change in the accretion pattern, the optical light curve of AMHer does not change between normal and reversed mode (Mazeh et al. 1986). These observations strongly suggest that only the high mass flow densities may switch back and forth between the poles, while the upper pole is fed during both states by low mass flow densities - associated with the emission of hard X-ray and cyclotron emission. This interpretation supports our one-pole cyclotron emission model.
An ultimate puzzle, though, remains: while our optical data from 1998,
August 20/22 presented here are completely typical of AMHer in
either the normal mode or the reversed mode, they were obtained only
about a week after our BeppoSAX observations (1998, August 12) which
showed the system in an hitherto unknown state with significant hard
X-ray emission during
(Matt et al. 2000). RXTE
and EUVE observations obtained on 1998, August 4 confirmed the
atypical X-ray light curve (Christian 2000). We are left with the
possibility that AMHer returned into one of the previously known
- normal or reversed - modes by the time of our optical campaign.
Considering that also the optical high state light curve of AMHer
noticeably changed on a (to our knowledge) single occasion
(Szkody 1978), it appears that AMHer occasionally drops into
short-lived "atypical'' behaviour. A multiwavelength campaign
covering one of these "atypical'' states could reveal whether they are
related to the switching between the normal and the reversed accretion
mode.
Acknowledgements
We thank Klaus Beuermann for many discussions on this topic and the referee Steve Howell for a number of helpful comments. This research was supported by the DLR under grant 50OR99036.