A&A 372, 557-562 (2001)
DOI: 10.1051/0004-6361:20010522

A model for the optical high state light curve of AMHerculis[*]

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 $\sim $ $5\times10^{16}~{\rm cm^2}$ and sustaining an inflow of $\sim $0.06 $\rm g\;cm^{-2}\,s^{-1}$. 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


1 Introduction

In AMHerculis stars, or polars, a strongly magnetic white dwarf ( $B\mathrel{\mathchoice {\vcenter{\offinterlineskip\halign{\hfil
$\displaystyle ...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 ( $\dot m\mathrel{\mathchoice {\vcenter{\offinterlineskip\halign{\hfil
$\displayst...
...criptscriptstyle ...), 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 ( $\dot m>10\,\mbox{$\rm g\;cm^{-2}\,s^{-1}$ }$), 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 ( $\mbox{$P_{\rm spin}$ }=\mbox{$P_{\rm orb}$ }$). 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.

2 Observations

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 $\phi\,(=\mbox{$\phi_{\rm orb}$ })=0.0$ as the inferior conjunction of the secondary star (Southwell et al. 1995)[*]. In AMHer, the conversion between the phase conventions is given by $\mbox{$\phi_{\rm orb}$ }=\mbox{$\phi_{\rm mag}$ }+0.367$. The high precision of the ephemeris of Heise & Verbunt (1988) ensures that the error in the computed phases is $\delta\mbox{$\phi_{\rm orb}$ }\le2\times10^{-3}$, 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).

3 Modelling the light curve

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 $\mbox{$\phi_{\rm orb}$ }\approx0.0$ and a secondary minimum near $\mbox{$\phi_{\rm orb}$ }\approx0.5$. 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.


 

 
Table 1: Optical photometry of AMHer in 1998. Listed are the start time and the the total exposure time.

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



  \begin{figure}
\par\includegraphics[angle=270,width=8.7cm]{ms10055.f1}
\end{figure} 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

3.1 The system geometry

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, $\beta$, 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, $\psi$. For our present purpose, the exact knowledge of $\psi$ is not crucial, as it defines only a phase offset between our model light curve and the observed light curve, and we will adjust $\psi$ accordingly below.

The inclination and the co-latitude of AMHer have been subject of dispute. Brainerd & Lamb (1985) derived from the observed polarization properties $i=35^{\circ}\pm5^{\circ}$and $\beta=58^{\circ}\pm5^{\circ}$. Another polarimetric study indicates $i\approx50^{\circ}$ and $\beta\approx50$-60$^{\circ}$(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 $i+\beta\mathrel{\mathchoice {\vcenter{\offinterlineskip\halign{\hfil
$\displays...
...erlineskip\halign{\hfil$\scriptscriptstyle ... 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 $i+\beta=105^{\circ}$.

We will use the combinations ( $i=50^{\circ},\ \beta=55^{\circ}$) and ( $i=35^{\circ},\ \beta=70^{\circ}$) in our models of the B and Vlight curves of AMHer in order to reflect the uncertainties in the exact geometry, while satisfying $i+\beta=105^{\circ}$.

3.2 The heated white dwarf

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 $\mbox{$T_{\rm wd}$ }=20\,000$K. They chose a temperature distribution decreasing linearly in angle from the spot centre merging into $\mbox{$T_{\rm wd}$ }$ at the spot radius. The best-fit spot parameters implied a central temperature of 47000K and a spot size of $f\sim0.09$ of the total white dwarf surface. With a distance of 90pc, the inferred radius of the white dwarf was $\mbox{$R_{\rm wd}$ }=1.12\times10^9$ 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).

  
3.3 The accretion stream

The near-ultraviolet emission ( $\lambda\mathrel{\mathchoice {\vcenter{\offinterlineskip\halign{\hfil
$\displays...
...offinterlineskip\halign{\hfil$\scriptscriptstyle ...Å) 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.

3.4 The secondary star

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.

3.5 The cyclotron emission

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).


  \begin{figure}
\par\includegraphics[angle=270,width=8.6cm]{ms10055.f2}
\end{figure} Figure 2: Sample cyclotron model spectra from the accretion column in AMHer for several angles $\vartheta $ between the magnetic field line and the line-of-sight and $\dot
m=0.06$ $\rm g\;cm^{-2}\,s^{-1}$. Plotted are also the B and V response functions (arbitarily scaled). See Fig. 3 for the variation of  $\vartheta (\phi _{\rm orb})$.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=8.5cm]{ms10055.f3}
\end{figure} Figure 3: Angle between the magnetic field line and the line-of-sight as function of the orbital phase. Angles $\vartheta $exceeding $90^{\circ }$ are expressed as $180^{\circ }-\vartheta $. The full line corresponds to $i=50^{\circ }$, $\beta =55^{\circ }$, the dashed to $i=35^{\circ }$, $\beta =70^{\circ }$.
Open with DEXTER


  \begin{figure}
\par {\includegraphics[width=4.1cm]{ms10055.f4} }\end{figure} Figure 4: Viewing geometry for an inclination $i=50^{\circ }$ and a co-latitude $\beta =55^{\circ }$. $\omega $ denotes the rotation axis, B the magnetic field line in the accretion column. The accretion column with a height of $\sim $0.1  $R_{\rm wd}$ is plotted as bold line segment. Indicated are the lines of sight for inferior conjunction ( $\mbox{$\phi_{\rm orb}$ }\approx0.0$, primary minimum) and superior conjunction ( $\mbox{$\phi_{\rm orb}$ }\approx0.5$, secondary minimum) of the secondary star.
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 $\sim $(1-3) $\times10^{-10}$  $\rm erg\;cm^{-2}\,s^{-1}$ (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 $\dot
M=10^{15}$-1016 $\rm g\;s^{-1}$ for the stand-off shock region[*]. The magnetic field strength of AMHer is $B=14\,\rm MG$(Bailey et al. 1991). For the accretion rate estimated above and a dipolar field geometry, a cross section of the accretion column of $A_{\rm acc} \approx 5\times10^{16}\,\rm cm^2$ is estimated near the magnetic pole (Lubow & Shu 1975, 1976; Heerlein et al. 1999). From these numbers, a mean local mass flow density $\dot m
\approx 0.02$-0.2  $\rm g\;cm^{-2}\,s^{-1}$ 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 $\vartheta $ 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 $\dot m$ 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, $\vartheta $. For a given set of inclination i, co-latitutde $\beta$ and orbital phase $\mbox{$\phi_{\rm orb}$ }$

\begin{displaymath}\cos \vartheta = \cos i \cos\beta - \sin i \sin\beta\cos(\mbox{$\phi_{\rm orb}$ }+\pi/2-\psi).
\end{displaymath} (1)

Figure 3 shows the variation of $\vartheta $ with $\mbox{$\phi_{\rm orb}$ }$for ( $i=50^{\circ},\ \beta=55^{\circ}$) and for ( $i=35^{\circ},\ \beta=70^{\circ}$), the viewing geometry is illustrated in Fig. 4. We chose $\psi=-12^{\circ}$ so that the minimum in $\vartheta $ aligns with the minimum of the observed Vlight curve, i.e. the magnetic pole leads the secondary by $\psi$. This value is in good agreement with the azimuth of the hot accretion region derived from X-ray and ultraviolet light curves, $\psi_{\rm HS}\approx-10^{\circ}{\rm ~to~}0^{\circ}$ (Paerels et al. 1996; Gänsicke et al. 1998). $\vartheta $ reaches a minimum near $\mbox{$\phi_{\rm orb}$ }\approx-0.03$, when the observer looks down along the accretion funnel, passes through $90^{\circ }$ at $\mbox{$\phi_{\rm orb}$ }\approx0.3$ and reaches the maximum of $105^{\circ}$ at $\mbox{$\phi_{\rm orb}$ }\approx0.47$. As the cyclotron beaming is sensitive only to $\vert\cos\vartheta\vert$ we plot $180^{\circ }-\vartheta $when $\vartheta>90^{\circ}$. From Fig. 3 it is clear that a reduction of the cyclotron flux is expected in the range $\mbox{$\phi_{\rm orb}$ }\approx0.3$-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 $\mbox{$\phi_{\rm mag}$ }\approx0.0$ and 0.3, which corresponds to the phases when $\vartheta\simeq90^{\circ}$.

  
3.6 Composite model light curves

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 $\dot m$ which determines the (phase-dependent) contribution of the cyclotron emission and the (phase-independent) contribution of the accretion stream, $B_{\rm stream}$ and $V_{\rm stream}$. 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 $\dot
m=0.06$  $\rm g\;cm^{-2}\,s^{-1}$, $B_{\rm stream}=15.0$, and $V_{\rm stream}=14.1$ (Figs. 56). Figures 5 and 6 show the model light curves for the two different choices of ($i, \beta$). The contribution of the heated white dwarf to the observed B and V high state light curves is small (Figs. 56), 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 $\mbox{$\phi_{\rm orb}$ }=0.0$ and a shallow secondary minimum at $\mbox{$\phi_{\rm orb}$ }=0.5$ in V and practically no orbital modulation in B. However, for ( $i=50^{\circ},\ \beta=55^{\circ}$) the angle $\vartheta $ drops below $40^{\circ}$ during $\phi=0.8$-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.


  \begin{figure}
\par\includegraphics[angle=270,width=8.4cm]{ms10055.f5}
\end{figure} Figure 5: The observed light curves from Fig. 1 have been sampled into $\Delta\mbox{$\phi_{\rm orb}$ }=0.03$ orbital phase bins (grey points). The composite model light curves assuming $i=50^{\circ }$, $\beta =55^{\circ }$ are shown as thick black lines. Contributions to the total fluxes come from the heated white dwarf, from the accretion column, and from the accretion stream (the secondary star contributes too little to show up in the graph). The angles $\vartheta $ indicate the angle between the magnetic field and the line-of-sight (cf. Fig. 3).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=8.4cm]{ms10055.f6}
\end{figure} Figure 6: As Fig. 5, but assuming $i=35^{\circ }$, $\beta =70^{\circ }$.
Open with DEXTER

  
4 Discussion

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. $\gg$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 ${\rm H\alpha}$ and ${\rm H\beta}$. 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 $B\approx14$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 $i+\beta=105^{\circ}$ structures extending less than than $\sim $0.035  $R_{\rm wd}$ above the white dwarf surface will be totally eclipsed at $\mbox{$\phi_{\rm orb}$ }\approx0.5$. However, for the low magnetic field and the rather low mass flow density, the shock is expected to stand $\sim $0.1  $R_{\rm wd}$ 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 $\mbox{$\phi_{\rm orb}$ }\approx0.5$ (i.e. when $\vartheta $ 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 $\mbox{$\phi_{\rm orb}$ }\approx0.6$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 $\vartheta(\mbox{$\phi_{\rm orb}$ })$, 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 $\mbox{$\phi_{\rm orb}$ }\simeq0.5$ (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.

References

 


Copyright ESO 2001