Formation and evolution of inclined accretion disks in intermediate polars

The results of 3D modeling of the formation of the accretion disks of intermediate polars are presented. A model with misaligned rotation axes of accretor and the orbit is onsidered, in which it is assumed that the white dwarf has a dipolar magnetic field with its symmetry axis inclined to the whitedwarf rotation and orbital axes. The computations show that, in the early stages of formation of the disk, the action of magnetic field is able to create the initial (seed) inclination of the disk. This inclination is then supported mainly by the dynamical pressure of the flow from the inner Lagrangian point L₁. As themass of the disk increases, the inclination disappears. Under certain conditions, the disk inclination does not arise in systems with misaligned white-dwarf rotation and orbital axes. The influence of the magnetic field and asynchronous rotation of the accretor may result in the formation of spiral waves in the disk with amplitudes sufficient to be detected observationally.     

A possible mechanism for the formation of tilted disks in intermediate polars

Using 3D gas dynamics, we numerically simulate accretion-disk formation in typical cataclysmic variable intermediate polars with dipolar magnetic fields (B _a = 10⁵?5 × 10⁵ G) and misaligned white-dwarf magnetic and rotation axes. Our simulations confirm that a significant misalignment of the axes results in a significant misalignment of the disk to the orbital plane. However, over time, this disk tilt disappears: early in the simulation, the initial particle positions in the rarefied tilted disk are governed solely by the magnetic field of the white dwarf. Due to the increasing disk mass and hence increasing disk gas pressure, the tilted disk eventually becomes decoupled from the magnetic field. The tidal action of the donor leads to a retrograde (i.e., nodal) precession of the tilted disk’s streamlines, and the disk becomes twisted. When the disk tilt is greater than 4°, the incoming gas stream no longer strikes the disk rim (i.e., bright shocked region). Matter is now transported over and under the disk rim to the inner regions of the disk. Over time, the increased mass of inner parts of the disk due to the action of the colinear gas stream returns the inner-disk regions to a colinear configuration. Meanwhile, the outer regions of the tilted, twisted disk become warped. Our simulations suggest that the lifetime of an intermediate polar’s tilted disk could be several tens to thousands of orbital periods.     

Three-Dimensional Numerical Simulation of a Flow Structure in the Asynchronous Polar CD Ind in the Approximation of an Offset Dipole Magnetic Field of a White Dwarf

The flow structure in the asynchronous polar CD Ind was investigated in the approximation of an offset dipole. Computation results made it possible to identify such system features as the drift of hot spots over the white dwarf surface, the flow structure dependence on the phase of the beat period, and the magnetic pole switching. To study the flow structure, a three-dimensional numerical MHD model based on the approximation of modified magnetic hydrodynamics was used. Numerical calculations were performed for ten phases of the beat period at a constant rate of mass transfer. In addition, for a more detailed study of the poles switching, additional calculation series were carried out in the corresponding phase ranges of the beat period. The energy release zones during the spin-orbit period shift in longitude on average by 20°, which corresponds to 0.05 of the orbital period phase. The magnetic pole switching occurs in a time that does not exceed 0.1 of the beat period. In the middle of this process, accretion occurs on both poles with the same intensity, and the flow resembles an arch. Based on the calculation results, synthetic light curves were constructed for the optical range of the spectrum.     

Thermal cyclotron radiation by isolated magnetic white dwarfs and constraints on the parameters of their coronas

An efficient method for the detection and estimation of the parameters of the coronas of isolated white dwarfs possessing magnetic fields of about 10⁷ G is tested. This method is based on the detection of thermal radiation of the coronal plasma at harmonics of the electron gyrofrequency, which is manifest as a polarized infrared excess. The Stokes parameters for the thermal cyclotron radiation from the hot corona of a white dwarf with a dipolar magnetic field are calculated. A new upper limit for the electron density, 10¹⁰ cm^?3, in a corona with a temperature of ?10⁶ K is found for the white dwarf G99-47 (WD 0553+053). This limit is a factor of 40 lower than the value derived earlier from ROSAT X-ray observations. Recommendations for subsequent infrared observations of isolated magnetic white dwarfs aimed at detecting their coronas or deriving better constraints on their parameters are presented.     

Spectropolarimetry and IR photometry of magnetic white dwarfs: Vacuum polarization or rydberg states in their magnetic fields?

We suggest two mechanisms to explain IR photometric and spectropolarimetric observations of magnetic white dwarfs: vacuum polarization and the existence of Rydberg atomic states with large dipole moments arising due to atomic collisions in the strong magnetic field of the white dwarf (so-called magnetic collision-induced absorption, or magnetic CIA). Both mechanisms can explain the observed rotations of the polarization ellipses and the depression of the IR spectral energy distribution. We present the results of spectropolarimetric observations of several magnetic white dwarfs with the Special Astrophysical Observatory 6-m telescope, together with photometric observations in the near-IR obtained with the Russian-Italian AZT-24 Telescope at Campo Imperatore.     

Features of the Flow Structure in the Vicinity of the Inner Lagrangian Point in Polars

The structures of plasma flows in close binary systems whose accretors have strong intrinsic magnetic fields are studied. A close binary system with the parameters of a typical polar is considered. The results of three-dimensional numerical simulations of the matter flow from the donor into the accretor Roche lobe are presented. Special attention is given to the flow structure in the vicinity of the inner Lagrangian point, where the accretion flow is formed. The interaction of the accretion-flow material from the donor’s envelope with the magnetic field of the accretor results in the formation of a hierarchical structure of the magnetosphere, because less dense areas of the accretion flow are stopped by the magnetic field of the white dwarf earlier than more dense regions. Taking into account this kind of magnetosphere structure can affect analysis results and interpretation of the observations.     

Impact of the magnetic field on the structure of accretion disks in semi-detached binaries

We discuss characteristic features of the magnetic gas-dynamical structure of the flows in a semi-detached binary system obtained from three-dimensional simulations, assuming that the intrinsic magnetic field of the accreting star is dipolar. The turbulent diffusion of the magnetic field is taken into account. The SS Cyg system is considered as an example. Including the magnetic field can alter the basic parameters of the accretion disk, such as the accretion rate and the characteristic density. The magnetic field in the disk is primarily toroidal.     

Active longitudes and north-south asymmetry of the activity the Sun as manifestations of its relic magnetic field

The distributions of dominant magnetic polarities in synoptic maps of photospheric magnetic fields and their extrapolations to the corona based on Stanford Observatory data are studied. Both dipolar and quadrupolar magnetic patterns are detected in the distributions of dominant polarities in the near-equatorial region of the photosphere for activity cycles 21, 22, and 23. The field in these patterns often has opposite signs on opposite sides of the equator, with this sign changing from cycle to cycle. A longitude-time analysis of variations of the mean solar magnetic field shows that the contribution of the large-scale magnetic patterns to the total field does not exceed 20 µT. The most stable magnetic structures at a quasi-source surface in the solar corona are separated by approximately 180° in heliographic longitude and are close to dipolar. The nature and behavior of these large-scale magnetic patterns are interpreted as a superposition of cyclic dynamo modes and the nonaxially symmetric relic field of the Sun. The contribution of the relic field to the mean solar magnetic field appears as a weak but stable rotational modulation whose amplitude does not exceed 8 µT.     

Features of the accretion in the EX Hydrae system: Results of numerical simulation

A two-dimensional numerical model in the axisymmetric approximation that describes the flow structure in the magnetosphere of the white dwarf in the EX Hya system has been developed. Results of simulations show that the accretion in EX Hya proceeds via accretion columns, which are not closed and have curtain-like shapes. The thickness of the accretion curtains depends only weakly on the thickness of the accretion disk. This thickness developed in the simulations does not agree with observations. It is concluded that the main reason for the formation of thick accretion curtains in the model is the assumption that the magnetic field penetrates fully into the plasma of the disk. An analysis based on simple estimates shows that a diamagnetic disk that fully or partially shields the magnetic field of the star may be a more attractive explanation for the observed features of the accretion in EX Hya.     

Photometry of the intermediate polar 1RXS J062518.2+733433: Characteristics of its periodic brightness variations

We present the results of our photometry of the recently discovered intermediate polar 1RXS J062518.2+733433. The observations were made using the 70 cm telescope of the Astronomical Observatory of the Ural State University with a multichannel photometer, and were carried out during 13 nights in March–May 2004. Our analysis reveals brightness variations with periods of 19.788 ± 0.003 min, corresponding to the rotational period of the white dwarf, 21.273 ± 0.003 min, corresponding to the orbital sideband, and 283.3 ± 0.5 min, which is the orbital period. The variations with the white dwarf’s rotational period show a stable amplitude and a quasi-sinusoidal pulse shape persistent over a long time. In contrast, the orbital-sideband variations have an unstable amplitude and a significantly nonsinusoidal pulse shape that varies with time. Variations of the amplitude and pulse shape of the orbital-sideband variations can be related to structural changes of the accretion disk surrounding the white dwarf.     

Formation of accretion disks in close-binary systems with magnetic fields

We have developed a three-dimensional numerical model and applied it to simulate plasma flows in semi-detached binary systems whose accretor possesses a strong intrinsic magnetic field. The model is based on the assumption that the plasma dynamics are determined by the slow mean flow, which forms a backdrop for the rapid propagation of MHD waves. The equations describing the slow motion of matter were obtained by averaging over rapidly propagating pulsations. The numerical model includes the diffusion of magnetic field by current dissipation in turbulent vortices, magnetic buoyancy, and wave MHD turbulence. A modified three-dimensional, parallel, numerical code was used to simulate the flow structure in close binary systems with various accretor magnetic fields, from 10⁵ to 10⁸ G. The conditions for the formation of the accretion disk and the criteria distinguishing the two types of flow corresponding to intermediate polars and polars are discussed.     

AE Aquarii represents a new subclass of Cataclysmic Variables

We analyze properties of the unique nova-like star AE Aquarii identified with a close binary system containing a red dwarf and a very fast rotating magnetized white dwarf. It cannot be assigned to any of the three commonly adopted sub-classes of Cataclysmic Variables: Polars, Intermediate Polars, and Accreting non-magnetized White Dwarfs. Our study has shown that the white dwarf in AE Aqr is in the ejector state and its dipole magnetic moment is ???1.5 × 10³⁴ G cm³. It switched into this state due to intensive mass exchange between the system components during a previous epoch. A high rate of disk accretion onto the white dwarf surface resulted in temporary screening of its magnetic field and spin-up of the white dwarf to its present spin period. Transition of the white dwarf to the ejector state had occurred at a final stage of the spin-up epoch as its magnetic field emerged from the accreted plasma due to diffusion. In the frame of this scenario AE Aqr represents a missing link in the chain of Polars evolution and the white dwarf resembles a recycled pulsar.     

Spectropolarimetric Observations of Magnetic White Dwarfs with the SAO 6-m Telescope

The results of spectropolarimetric observations of a number of magnetic white dwarfs obtained on the 6-m optical telescope of the Special Astrophysical Observatory are presented. The observations were carried out using the SCORPIO focal aperture-ratio reducer in a spectropolarimetric regime. Two characteristic dependences of the degree of polarization on the wavelength are observed. For one group of objects, the degree of linear polarization grows with wavelength, suggesting that the alignment of atoms and molecules in Rydberg states in the atmosphere of the white dwarf due to the action of its magnetic field influences scattering processes. The second group of objects displays an increase in the degrees of both linear and circular polarization with wavelength, providing evidence for the presence of protoplanetary disks around these magnetic white dwarfs, in which the alignment of circumstellar grains leads to the observed behavior.     

Magnetic-field structure in the accretion disks of semi-detached binary systems

The results of three-dimensional MHD numerical simulations are used to investigate the characteristic properties of the magnetic-field structures in the accretion disks of semi-detached binary systems. It is assumed that the intrinsic magnetic field of the accretor star is dipolar. Turbulent diffusion of the magnetic field in the disk is taken into account. The SS Cyg system is considered as an example. The results of the numerical simulations show the intense generation of a predominantly toroidal magnetic field in the accretion disk. Magnetic zones with well defined structures for the toroidal magnetic field form in the disk, which are separated by current sheets in which there ismagnetic reconnection and current dissipation. Possible observational manifestations of such structures are discussed. It is shown that the interaction of a spiral precessional wave with the accretor’s magnetosphere could lead to quasi-periodic oscillations of the accretion rate.     

Population synthesis of rapidly rotating main-sequence stars with companions

Usingthe “Scenario Machine” (a specialized numerical code formodeling the evolution of large ensembles of binary systems), we have studied the physical properties of rapidly rotating main-sequence binary stars (Be stars) with white-dwarf companions and their abundance in the Galaxy. The calculations are the first to take into account the cooling of the compact object and the effect of synchronization of the rotation on the evolution of Be stars in close binaries. The synchronization time scale can be shorter than the main-sequence lifetime of a Be star formed during the first mass transfer. This strongly influences the distribution of orbital periods for binary Be stars. In particular, it can explain the observed deficit of short-period Be binaries. According to our computations, the number of binary systems in the Galaxy containing a Be star and white dwarf is large: 70–80% of all Be stars in binaries should have degenerate dwarf companions. Based on our calculations, we conclude that the compact components in these systems have high surface temperatures. Despite their high surface temperatures, the detection of white dwarfs in such systems is hampered by the fact that the entire orbit of the white dwarf is embedded in the dense circumstellar envelope of the primary, and all the extreme-UV and soft X-ray emission of the compact object is absorbed by the Be star’s envelope. It may be possible to detect the white dwarfs via observations of helium emission lines of Be stars of not very early spectral types. The ultraviolet continuum energies of these stars are not sufficient to produce helium line emission. We also discuss numerical results for Be stars with other evolved companions, such as helium stars and neutron stars, and suggest an explanation for the absence of Be-black-hole binaries.     

Magneto-rotational instability in the accreting envelope of a protostar and the formation of the large-scale magnetic field

We investigate the role of the magnetic field in the collapse of a gas-dust cloud into a massive gravitating object. Observations of one such object (G31.41) indicate that the magnetic field has an hourglass shape oriented along the rotation axis of the matter, due to the freezing-in of the magnetic-field lines in the accreting matter. It is believed that accretion in stellar disks is associated with the transport of angular momentum from the center to the periphery, which could be initiated by large-scale vortex structures arising in the presence of unstable rotational flows of matter. The numerical simulations have established that the equilibrium configuration of a gas-dust disk rotating in a spherically symmetrical gravitational potential is subject to the development of strong instability in the presence of a weak magnetic field. It is shown that the development of instability leads to a transport of angular momentum to the disk periphery by large-scale vortex structures, together with the accretion of matter onto the gravitating object. The magnetic-field lines near the equator take on a chaotic character, but an hourglass configuration is observed near the rotation axis, in agreement with observations.     

The influence of the inclination of the accretor’s magnetic axis on the accretion-disk structure in intermediate polars

Typical changes of the accretion-disk structures in intermediate polars are studied as a function of the inclination of the accretor’s magnetic field. Thre-dimensional numerical modeling was performed for seven differentmagnetic-axis inclinations. The results showthat the system forms a magnetosphere region, and that column accretion occurs. The action of the magnetic field tilts the inner parts of the disk along the magnetic axis of the accretor. The magnetic-field inclination appreciably influences matter transfer in the disk and accretion processes. Generation of toroidal magnetic field, magnetic braking, and alignment of the dipole magnetic field result in oscillations of the accretion rate. A direct relationship between the field inclination and the oscillation amplitude is found, as well as an inverse relationship between the field inclination and the oscillation period.     

Periodic instabilities in protostellar disks with azimuthal magnetic fields

The stability of magnetohydrodynamic oscillations in a protostellar disk with a toroidal magnetic field is analyzed. It is shown that, apart from the aperiodic magnetorotational instability, two other types of periodic instabilities of non-axisymmetric perturbations can exist. The simultaneous presence of azimuthal and vertical components of the wave vector are necessary for these to exist. One instability is due to the inductive winding-up of the azimuthal magnetic field of the wave, and the other arises when the field amplitude is increased by a comoving Hall wave, transferring magnetic field into a region of enhanced field intensity. The bandwidths of the unstable wave numbers are analyzed as a function of the Hall current, the β parameter of a plasma, and the angle between the direction of wave propagation and the plane of the disk. Regions in the accretion disks typical of T Tauri stars are indentified where these instabilities could be most active.     

Features of the matter flows in the peculiar cataclysmic variable AE Aquarii

The structure ofmatter flows in close binary systems in which one of the components is a rapidly rotating magnetic white dwarf is studied. Themain example considered is the AEAquarii system; the period of the white dwarf’s rotation is about a factor of 1000 shorter than the orbital period, and the magnetic field on the white-dwarf surface is of order 50MG. The matter flows in this system were analyzed via numerical solution of a systemofmagnetohydrodynamical equatons. These computations show that the white dwarf’s magnetic field does not significantly influence the velocity field of the matter in its Roche lobe in the case of a laminar flow regime, so that the field does not hinder the formation of a transient disk (ring) surrounding the magnetosphere. However, the efficiency of the energy and angular-momentum exchange between the white dwarf and the surrounding matter increases considerably with the development of turbulent motions in the matter, accelerating the matter at the magnetosphere boundary and leading to a high escape rate from the system. The time scales for the system’s transition between the laminar and turbulent modes are close to those for the transition of AE Aquarii between its quiet and active phases.     

Light-curve synthesis for close binaries: An extended shock as an analog of a hot spot in cataclysmic variables

An algorithm is presented for the synthesis of the light curve of a close binary system consisting of a red dwarf that fills its Roche lobe and a spherical white dwarf. The spherical component is surrounded by an elliptical, geometrically thick accretion disk. The code models an extended shock located along the edge of the stream near the outer boundary of the disk. The observational manifestations of the shock show that it can be considered as an analog of a hot spot at the edge of the disk. Synthetic light curves for the SU UMa system OY Car at various phases of its activity indicate that the model can describe both typical and peculiar light curves for this cataclysmic variable reasonably well.