arxiv:astro-ph/ v1 6 Oct 2004

Transcription

1 **FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Outbursts of WZ Sagittae arxiv:astro-ph/ v1 6 Oct 2004 O. M. Matthews 1, M. R. Truss 2, G. A. Wynn 1, R Speith 3 1 Department of Physics & Astronomy, University of Leicester, University Road, Leicester, LE1 7RH, United Kingdom 2 School of Physics & Astronomy, University of St Andrews, North Haugh, St Andrews, Fife, KY16 9SS, United Kingdom 3 Institut für Astronomie und Astrophysik, Universität Tübingen, Auf der Morgenstelle 10C, D Tübingen, Germany Abstract. The outbursts of the dwarf nova WZ Sagittae occur on timescales which are much longer than expected from the standard thermal viscous disc instability model, which correctly reproduces the outburst behaviour of many other dwarf novae. In addition short echo outbursts have sometimes been observed after the main outburst. The increased timescales could be explained by a magnetic propeller, anchored to the primary star, which clears the inner disc. We argue that the propeller, combined with the resonances of the binary orbit, may also be the cause of the echo outbursts. We present numerical models which reproduce the echo outbursts. 1. Introduction WZSagittaeisamemberoftheSUUrsaeMajorisclassofdwarfnovae. Theseare short period systems which exhibit long superoutbursts in addition to ordinary outbursts. During these superoutbursts periodic modulations of the optical light curve, known as superhumps are observed. These features have a frequency which is slightly longer than that of the binary orbit, and which is usually attributed to a precessing disc. Superoutbursts in the SU UMa systems typically occur with recurrence times of weeks or months. Stars from the subclass to which WZ Sge belongs are sometimes known as TOADs or Tremendous Outburst Amplitude Dwarf novae. These systems have very long recurrence times and correspondingly bright superoutbursts. WZ Sge itself showed a 33 year recurrence time, until this pattern was broken by the 2001 outburst which occurred after only 23 years. WZ Sge, the brightest of all dwarf novae, is also notable for its extreme mass ratio. According to Skidmore et al. (2002) M 2 = ± M and M 1 = 0.74 ± 0.07 M. It has a binary period of P orb = 81.6 min (Krzeminski 1962) and a primary spin period of P spin = s. In this outburst (e.g Patterson et al. 2002) superhumps were observed, but prior to these a modulation at the frequency of the binary orbit was observed by Patterson et al. (2002) as discussed by e.g. Kunze & Speith (2004). In addition, the main outburst was followed by twelve echo outbursts. Echoes have also been observed in the similar system EG Cancri (Osaki, Meyer & Meyer-Hofmeister 2001). 1

2 2 Matthews et al. Numerical models incorporating the thermal viscous disc instability model (DIM) due to Hoshi (1979), have reproduced the outburst behaviour of many dwarf noave. These models have not however, without using extremely low viscosities, been able to explain the long recurrence time of WZ Sge type systems. Osaki (1989) extended the DIM to binaries with small mass ratios when Whitehurst(1988) showed that discs in such systems could become tidally unstable. Osaki proposed that, during the normal outburst cycle, material is allowed to build up in the outer disc. When the disc becomes sufficiently large that it may access the 3:1 resonance of the binary orbit, it will become eccentric, enhancing dissipation and the loss of angular momentum, leading to a superoutburst. Simulations (e.g. Truss, Murray & Wynn 2001) supported this theory. The extended recurrence times of the WZ Sge type systems may be explained by the use of a magnetic propeller to evacuate the inner regions of the disc. This mechanism, similar to that proposed by Matthews, Speith & Wynn (2004a)fortheFUOrionisstars, isdetailedintheforthcomingpaperbywynn et al. (2004). The magnetic propeller is a mechanism by which the inner regions of an accretion disc can be evacuated. The co-rotation radius R co is the radius at whichakeplerian orbithas thesameangularfrequencyas that of thestar s spin. If the primary magnetic field rotates as a solid body then this radius can be said to divide the magnetic accretion and magnetic propeller regimes. Particles at radii smaller than the co-rotation radius lose angular momentum as a result of their interactions with the stellar magnetic field, and fall rapidly onto the star. This magnetically enhanced accretion occurs because the particles are moving more rapidly than the magnetic field. Elements of plasma which exist outside corotation are, conversely, accelerated by the magnetic field, and gain angular momentum at the star s expense. This is the magnetic propeller regime. An example of a strong magnetic propeller is that of AE Aquarii (Wynn, King & Horne 1997), in which accretion is completely prevented by the magnetic field. What is proposed here is a far weaker propeller, driven by a rapidly spinning, moderately magnetic white dwarf. Simulations using weaker fields have been shown to warp accretion discs (O Sullivan et al. 2004). Lasota, Kuulkers & Charles (1999) also used a magnetically truncated disc in their model of WZ Sge. The cause of the echo outbursts is as yet uncertain. Osaki, Meyer & Meyer-Hofmeister (2001) suggested that the mass which accumulates near the 3:1 resonance in the disc is released gradually after the main outburst, leading to echo outbursts. In this paper we model a full outburst of WZ Sge using smoothed particle hydrodynamics (SPH). The rapidly spinning white dwarf creates a magnetic propeller which, in addition to evacuating the inner disc, aids the creation of a reservoir of material in the outer disc. Echo outbursts arise naturally from this configuration. 2. Numerical Method The simulations in this paper were performed using a fully three dimensional SPH code with viscosity switching and a magnetic acceleration as described below. The specific code is fully described by Truss et al. (2000) whilst Monaghan (1992) provides a general review of SPH. SPH is a Lagrangian technique for the modelling of fluid dynamics. It works by discretising the fluid into a series of

3 Outbursts of WZ Sagittae 3 particles, normally with equal mass, which represent fluid elements. These particles carry values for fluid properties such as density. The properties of the fluid may then be interpolated for any point in space by taking a weighted average of the properties of nearby particles. This is done using a smoothing kernel, which usually takes the form of a cubic spline function. The viscosity in the disc is an artificial viscosity which is introduced into the velocity and energy equations of the SPH formulation. The extra term is composed of shear and bulk viscosity elements. However, providing that the local velocity divergence is negligible, it can be shown(murray 1996) that the viscosity acts inamannervery similar totheshakura & Sunyaev (1973)αshearviscosity. The SPH viscosity parameter ζ, which is proportional to α, is switched between two values, mimicking the behaviour of the DIM. Switching occurs at two trigger surface densities. Any annulus of the disc which exceeds Σ crit high will begin to have its ζ gradually altered to the high state value. Conversely, an annulus which falls below Σ crit low will have its ζ value slowly reduced to the lower, quiescent value. This is discussed by Truss et al. (2000). The positions of the triggers, which are roughly proportional to radius, are dervied by Cannizzo, Shafter & Wheeler (1988). Dendy (1990) expresses the magnetic tension acceleration as a mag 1 ρr c ( ) Bz B φ 4π, (1) where ρ is density and r c represents the radius of curvature of the magnetic field lines. The azimuthal and vertical components of the magnetic field are representedbyb φ andb z respectively. Wewill maketheapproximation that the radius of curvature of the field lines is of the order of the disc scale-height thus, r c H (e.g. Pearson, Wynn & King 1997). The ratio of vertical and azimuthal field strengths is related to the shear between the disc and the magnetic field. If it is assumed that the field rotates rigidly with the star then this ratio can be expressed in the form (e.g. Livio & Pringle 1992) B φ B z (Ω k Ω ) Ω k, (2) where Ω indicates the angular frequency of the accreting star while Ω k represents the Keplerian angular frequency at a given angle. Equations (1) and (2) are combined. Additionally it is assumed that the field is dipolar so that the magnetic field can be related to the magnetic moment and distance from the star by the expression B z µ R 3 where R represents the radius of the star. It is also assumed that surface density can be approximated by Σ ρh, where H is the scale height of the accretion disc. This allows us to make the following approximation a mag B2 z (Ω k Ω ) = µ2 R 6 (v k v ) k m R 6(v k v ). (3) 4πρH Ω k 4πΣ v k v k In this fomulation the acceleration falls to zero when v k = v at co-rotation as expected. It may be useful to compare the value constant k m with the β used

4 4 Matthews et al. Figure 1. Upper frame: Light curve of WZ Sge in 2001, showing a rapid rise and subsequent decay punctuated by an episode of 12 echo outbursts. A freehand curve has been added to improve visibility. This article originally appeared in Patterson et al. (2002). Copyright 2002, Astronomical Society of the Pacific; reproduced with permission of the Editors. Lower frame: Results of simulations showing the occurrence of echo outbursts. in Matthews et al. (2004b). It can be shown that k m = β GM 2Σ where β = µ 2 2π GM, (4) where G is the universal constant of gravitation and M is the mass of the star. 3. Results Simulations were performed using the code described above. The results presented used a very high viscosity, such that the Shakura-Sunyaev alpha parameters in the hot and cold states are 0.6 and 0.06 respectively. Thetriggers are also placed rather close together and these initial calculations were performed using a low resolution of 30,000 particles, since they were part of a large parameter study. The very high viscosity has the effect of speeding up the viscous processes and hence reducing processing time. The viscosity switching triggers were activated only after the disc had been filled. A weak magnetic field of B G was used. As shown in Figure 1, the main outburst in our simulation is followed by a series of echo outbursts which are of a similar form to those in the observation which is also plotted. The timescale for the echo outbursts is an order of magnitude shorter than that observed. This is as expected since physical disc

5 Outbursts of WZ Sagittae 5 Figure 2. Upper frames: Snapshots of the accretion disc from our simulations showing log surface density and the viscous parameter ζ = 10α respectively. The binary separation is represented by a in both plots. Lower frame: A plot of azimuthally averaged surface density as a function of radius also showing the viscosity switching triggers. alphas are thought to be of order 0.1 and The echoes do not at present switch off, but continue up to the end of the run. Figure 2 shows the state of the disc during the echo outburst phase. The upper panels show surface density and the viscous alpha in the disc. The two circles represent the 2:1 and 3:1 orbital resonances. The lower panel is a plot of the azimuthally averaged surface density as a function of radius. The central hole in the accretion disc, due to the magnetic propeller, is visible in all these plots. This latter plot also illustrates the positions of the two viscosity triggers in the simulation. Together the plots show how the echo outbursts occur in the simulations. After a main outburst, some material remains in the outer disc. Spiral waves are excited, as a result of the binary potential. The associated density enhancements, when augmented by continuing mass transfer, aid the transport of matter towards the centre where, as accretion is inhibited by the magnetic propeller, an echo outburst occurs. The echo fades when the propeller is overcome by the viscous forces and accretion onto the star resumes. Echoes do not appear if the simulation is performed without the magnetic field.

6 6 Matthews et al. 4. Discussion In this paper echo outbursts are reproduced in an SPH model of WZ Sge. These echoes are similar to those observed in the light curve of the system However, it is not yet clear whether this mechanism is likely to be the one which occurs in reality. In order to establish this there are a number of refinements which will be made to the model. Firstly the viscous alphas used in the model will be brought closer to their accepted values. It is not desirable to rely on time scaling arguments in this system since neither magnetic nor resonant effects scale with alpha. The 30,000 particles used are not sufficient to fully resolve such a disc in three dimensions, so higher resolution calculations will be made. In addition a mechanism for the termination of the echoes will be sought. Acknowledgements The computations reported here were performed using the UK Astrophysical Fluids Facility (UKAFF). Research in theoretical astrophysics at the University of Leicester is supported by a PPARC rolling grant. OMM gratefully acknowledges support through a PPARC research studentship. He would also like to thank the University of St Andrews School of Physics and Astronomy for their hospitality during a visit, when parts of this work were initiated. MRT acknowledges a PPARC postdoctoral fellowship. References Cannizzo J.K, Shafter R.A., Wheeler J.C., 1988, ApJ, 333, 227 Dendy R.O., 1990, Plasma Dynamics, Oxford University Press, Oxford Hoshi R., 1979, Prog. Theor. Phys., 61, 1307 Krzeminski W., 1962, PASP, 74, 66 Kunze S., Speith R., 2004, PASP, submitted Lasota J.-P., Kuulers E., Charles P., MNRAS, 305, 473 Livio M., Pringle J.E., 1992, MNRAS, 259, 23P Matthews O.M., Speith R., Wynn G.A., 2004, MNRAS, 347, 873 Matthews O.M., Speith R., Truss M.R., Wynn G.A., 2004, MNRAS, submitted Monaghan J.J., 1992, ARA&A, 30, 543 Murray J.R., 1996, MNRAS, 279, 402 Osaki Y., 1989, PASP, 41, 1005 Osaki Y., Meyer F., Meyer-Hofmeister E., 2001, A&A, 370, 488 O Sullivan M., Truss M.R., Walker C., Wood K., Matthews O.M., Whitney B.A., Bjorkman J.E., 2004, MNRAS, submitted Patterson J. et al., 2002, PASP, 114, 721 Pearson K.J., Wynn G.A., King A.R., 1997, MNRAS, 288, 421 Shakura N.I., Sunyaev R.A., 1973, A&A, 24, 337 Skidmore W., Wynn G.A., Leach R., Jameson R.F., 2002, MNRAS, 336, 1223 Truss M.R., Murray J.R., Wynn G.A., Edgar, R., 2000, MNRAS, 319,467 Truss M.R., Murray J.R., Wynn G.A., 2001, MNRAS, 324, L1 Whitehurst R., 1988, MNRAS, 232, 35 Wynn G.A., King A.R., Horne K., 1997, MNRAS, 286, 436 Wynn G.A., Leach R., Truss, M.R., Matthews O.M., 2004, MNRAS, in prep.