Possible quasi-periodic oscillations from unstable accretion: 3D magnetohydrodynamic simulations

We investigate the photometric variability of magnetized stars, particularly neutron stars, accreting through a magnetic Rayleigh–Taylor-type instability at the disc–magnetosphere interface, and compare it with the variability during stable accretion, with the goal of looking for possible quasi-periodic oscillations (QPOs). The light curves during stable accretion show periodicity at the star's frequency and sometimes twice that, due to the presence of two funnel streams that produce antipodal hotspots near the magnetic poles. On the other hand, light curves during unstable accretion through tongues penetrating the magnetosphere are more chaotic due to the stochastic behaviour of the tongues, and produce noisier power spectra. However, the power spectra do show some signs of quasi-periodic variability. Most importantly, the rotation frequency of the tongues and the resulting hotspots are close to the inner-disc orbital frequency, except in the most strongly unstable cases. There is therefore a high probability of observing QPOs at that frequency in longer simulations. In addition, the light curves in the unstable regime show periodicity at the star's rotation frequency in many of the cases investigated here, again except in the most strongly unstable cases which lack funnel flows and the resulting antipodal hotspots. The noisier power spectra result in the fractional rms amplitudes of the Fourier peaks being smaller.
We also study in detail the effect of the misalignment angle between the rotation and magnetic axes of the star on the variability, and find that at misalignment angles  ≳25°  the star's period always appears in the light curves.     

Magnetic activity in accretion disc boundary layers

We use three-dimensional magnetohydrodynamic simulations to study the structure of the boundary layer between an accretion disc and a non-rotating, unmagnetized star. Under the assumption that cooling is efficient, we obtain a narrow but highly variable transition region in which the radial velocity is only a small fraction of the sound speed. A large fraction of the energy dissipation occurs in high-density gas adjacent to the hydrostatic stellar envelope, and may therefore be reprocessed and largely hidden from view of the observer. As suggested by Pringle , the magnetic field energy in the boundary layer is strongly amplified by shear, and exceeds that in the disc by an order of magnitude. These fields may play a role in generating the magnetic activity, X-ray emission and outflows in disc systems where the accretion rate is high enough to overwhelm the stellar magnetosphere.     

Launching of conical winds and axial jets from the disc–magnetosphere boundary: axisymmetric and 3D simulations

We investigate the launching of outflows from the disc–magnetosphere boundary of slowly and rapidly rotating magnetized stars using axisymmetric and exploratory 3D magnetohydrodynamic simulations. We find long-lasting outflows in the following cases. (1) In the case of slowly rotating stars , a new type of outflow, a conical wind , is found and studied in simulations. The conical winds appear in cases where the magnetic flux of the star is bunched up by the disc into an X-type configuration. The winds have the shape of a thin conical shell with a half-opening angle  θ∼ 30°–40°  . About 10–30 per cent of the disc matter flows from the inner disc into the conical winds. The conical winds may be responsible for episodic as well as long-lasting outflows in different types of stars. There is also a low-density, higher velocity component (a jet) in the region inside the conical wind. (2) In the case of rapidly rotating stars (the 'propeller regime'), a two-component outflow is observed. One component is similar to the conical winds. A significant fraction of the disc matter may be ejected into the winds. The second component is a high-velocity, low-density magnetically dominated axial jet where matter flows along the opened polar field lines of the star. The jet has a mass flux of about 10 per cent of that of the conical wind, but its energy flux (dominantly magnetic) can be larger than the energy flux of the conical wind. The jet's angular momentum flux (also dominantly magnetic) causes the star to spin down rapidly. Propeller-driven outflows may be responsible for the jets in protostars and for their rapid spin-down. The jet is collimated by the magnetic force while the conical winds are only weakly collimated in the simulation region. Exploratory 3D simulations show that conical winds are axisymmetric about the rotational axis (of the star and the disc), even when the dipole field of the star is significantly misaligned.     

The magnetohydrodynamics model of twin kilohertz quasi-periodic oscillations in low-mass X-ray binaries

We suggest an explanation for the twin kilohertz quasi-periodic oscillations (kHz QPOs) in low-mass X-ray binaries (LMXBs) based on magnetohydrodynamics (MHD) oscillation modes in neutron star magnetospheres. Including the effect of the neutron star spin, we derive several MHD wave modes by solving the dispersion equations, and propose that the coupling of the two resonant MHD modes may lead to the twin kHz QPOs. This model naturally relates the upper, lower kHz QPO frequencies with the spin frequencies of the neutron stars, and can well account for the measured data of six LMXBs.     

Accretion to magnetized stars through the Rayleigh–Taylor instability: global 3D simulations

We present results of 3D simulations of magnetohydrodynamics (MHD) instabilities at the accretion disc–magnetosphere boundary. The instability is Rayleigh–Taylor, and develops for a fairly broad range of accretion rates and stellar rotation rates and magnetic fields. It manifests itself in the form of tall, thin tongues of plasma that penetrate the magnetosphere in the equatorial plane. The shape and number of the tongues changes with time on the inner disc dynamical time-scale. In contrast with funnel flows, which deposit matter mainly in the polar region, the tongues deposit matter much closer to the stellar equator. The instability appears for relatively small misalignment angles, Θ≲ 30°, between the star's rotation and magnetic axes, and is associated with higher accretion rates. The hotspots and light curves during accretion through instability are generally much more chaotic than during stable accretion. The unstable state of accretion has possible implications for quasi-periodic oscillations and intermittent pulsations from accreting systems, as well as planet migration.     

A magnetic scaleheight: the effect of toroidal magnetic fields on the thickness of accretion discs

We consider accreting systems in which the central object interacts, via the agency of its magnetic field, with the disc that surrounds it. The disc is turbulent and, so, has a finite effective conductivity. The field sweeps across the face of the disc, thereby forming a current that is directed radially within the disc. In turn, this disc current creates a toroidal field, where the interaction between the disc current and the toroidal field produces a Lorentz force that compresses the disc. We investigate this compression, which creates a magnetic scaleheight of the disc that can be much smaller than the conventional scaleheight. We derive an analytic expression for the magnetic scaleheight and apply it to fully ionized discs.     

Dwarf nova oscillations and quasi-periodic oscillations in cataclysmic variables – V. Results from an extensive survey

We present observations of dwarf nova oscillations (DNOs), longer-period dwarf nova oscillations (lpDNOs), and quasi-periodic oscillations (QPOs) in 13 cataclysmic variable stars. In the six systems, WW Cet, BP CrA, BR Lup, HP Nor, AG Hya and V1193 Ori, rapid, quasi-coherent oscillations are detected for the first time. For the remainder of the systems discussed, we have observed more classes of oscillations, in addition to the rapid oscillations they were already known to display, or previously unknown aspects of the behaviour of the oscillations. The period of a QPO in RU Peg is seen to change by 84 per cent over the 10 nights of the decline from outburst – the largest evolution of a QPO period observed to date. A period–luminosity relation similar to the relation that has long been known to apply to DNOs is found for lpDNOs in X Leo; this is the first clear case of the lpDNO frequency scaling with accretion luminosity. WX Hyi and V893 Sco are added to the small list of dwarf novae that have shown oscillations in quiescence.     

Dwarf nova oscillations and quasi-periodic oscillations in cataclysmic variables – I. Observations of VW Hyi

From archived and recent high-speed photometry of VW Hyi we find dwarf nova oscillations (DNOs) occasionally present throughout outburst, evolving from a 14.06-s period at maximum to >40 s near the end of outburst. A relatively slow increase of period is followed by a rapid increase and a subsequent decrease.
Quasi-periodic oscillations (QPOs) are seen at periods of hundreds of seconds. For the first time, the evolution of a QPO period is seen, increasing steadily during the final decline of an outburst. The occasional presence of two DNOs, separated in frequency by the QPO frequency, suggests reprocessing of the rotating DNO beam by a 'wall' rotating progradely in the disc at the QPO period.     

A boundary layer origin for dwarf nova oscillations

Dwarf nova oscillations (DNOs) have been observed in a number of cataclysmic variables. I propose that these oscillations could be produced by a non-axisymmetric bulge at the transition between the optically thick disc and the optically thin boundary layer region. This would naturally explain the observed oscillation periods and the dependence of the oscillation amplitude on photon energy. The transition radius moves inward and outward with changing mass accretion rate, which explains the correlation between period and flux, and the time-scale for period variations. The underlying cause of the non-axisymmetry that produces the oscillations is not known, so it is not possible to predict the oscillation amplitude from first principles.     

Dwarf nova oscillations and quasi-periodic oscillations in cataclysmic variables – II. A low-inertia magnetic accretor model

The dwarf nova oscillations observed in cataclysmic variable (CV) stars are interpreted in the context of a low-inertia accretor model, in which accretion on to an equatorial belt of the white dwarf primary causes the belt to vary its angular velocity. The rapid deceleration phase is attributed to propellering. Evidence that temporary expulsion rather than accretion of gas occurs during this phase is obtained from the large drop in extreme ultraviolet flux.
We show that the quasi-periodic oscillations are most probably caused by a vertical thickening of the disc, moving as a travelling wave near the inner edge of the disc. This alternately obscures and 'reflects' radiation from the central source, and is visible even in quite low inclination systems. A possible excitation mechanism, caused by winding up and reconnection of magnetic field lines, is proposed.
We apply the model, deduced largely from VW Hyi observations, to re-interpret observations of SS Cyg, OY Car, UX UMa, V2051 Oph, V436 Cen and WZ Sge. In the last of these we demonstrate the existence of a 742-s period in the light curve, arising from obscuration by the travelling wave, and hence show that the two principal oscillations are a dwarf nova oscillation and its reprocessed companion.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.