An analytic model for disc disruption by strong stellar magnetic fields

An analytic model is presented for the inner structure of an accretion disc in the presence of a strong stellar magnetic field. The model is valid inside the radius at which the electron scattering opacity starts to exceed the Kramers opacity. It illustrates how the increasing stellar poloidal field leads to an elevated disc temperature, ultimately causing a breakdown in the vertical equilibrium owing to rapidly increasing radiation pressure which cannot be balanced by the vertical stellar gravity. Viscous instability also occurs. The solution gives an accurate representation of numerical results, and enables useful expressions to be derived for the radius at which the disc is marginally thin and the radius at which viscous instability occurs. The disruption mechanism appears to have general validity for accretion discs around strongly magnetic stars.     

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.     

Magnetic torques between accretion discs and stars

I show in this paper that two types of magnetic torques can appear in the interaction between an accretion disc and a magnetic accretor. There is the well-known torque resulting from the difference in angular velocity between the accretion disc and the star, but in addition there is a torque coming from the interaction between the stellar magnetic field and the disc's own magnetic field. The latter form of magnetic torque decreases in strength more slowly with increasing radius, and will therefore dominate at large radii. The direction of the disc field is not determined by the difference in angular velocity between the star and the disc as in the Ghosh &38; Lamb model, but rather is a free parameter. The magnetic torque may therefore either spin up or spin down the star, and the torque changes sign if the magnetic field in the disc reverses. I suggest that this mechanism can explain the torque reversals that have been observed in some disc-fed X-ray pulsars.     

Accretion disc–stellar magnetosphere interaction: field line inflation and the effect on the spin-down torque

We calculate the structure of a force-free magnetosphere which is assumed to corotate with a central star and which interacts with an embedded differentially rotating accretion disc. The magnetic and rotation axes are aligned, and the stellar field is assumed to be a dipole. We concentrate on the case when the amount of field line twisting through the disc–magnetosphere interaction is large , and consider different outer boundary conditions. In general the field line twisting produces field line inflation (e.g. Bardou & Heyvaerts), and in some cases with large twisting many field lines can become open. We calculate the spin-down torque acting between the star and the disc, and we find that it decreases significantly for cases with large field line twisting. This suggests that the oscillating torques observed for some accreting neutron stars could be caused by the magnetosphere varying between states with low and high field line inflation. Calculations of the spin evolution of T Tauri stars may also have to be revised in the light of the significant effect that field line twisting has on the magnetic torque resulting from star–disc interactions.     

Global 3‐D solar‐type star‐disc dynamo systems: I. MHD modeling

We have carried out global three‐dimensional magnetohydrodynamic simulations of the star‐disc interaction region around a young solar‐type star. The magnetic field is generated and maintained by dynamos in the star as well as in the disc. The developing mass flows possess non‐periodic time‐variable azimuthal structure and are controlled by the nonaxisymmetric magnetic fields. Since the stellar field drives a strong stellar wind, accretion is anti‐correlated with the stellar field strength and disc matter is spiraling onto the star at low latitudes, both contrary to the generally assumed accretion picture. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Reconnection X-winds: spin-down of low-mass protostars

We investigate the interaction of a protostellar magnetosphere with a large-scale magnetic field threading the surrounding accretion disc. It is assumed that a stellar dynamo generates a dipolar-type field with its magnetic moment aligned with the disc magnetic field. This leads to a magnetic neutral line at the disc mid-plane and gives rise to magnetic reconnection, converting closed protostellar magnetic flux into open field lines. These are simultaneously loaded with disc material, which is then ejected in a powerful wind. This process efficiently brakes down the protostar to 10–20 per cent of the break-up velocity during the embedded phase.     

Discovery of drifting high-frequency quasi-periodic oscillations in global simulations of magnetic boundary layers

We report on the numerical discovery of quasi-periodic oscillations (QPOs) associated with accretion through a non-axisymmetric magnetic boundary layer in the unstable regime, when two ordered equatorial streams form and rotate synchronously at approximately the angular velocity of the inner disc. The streams hit the star's surface producing hotspots. Rotation of the spots leads to high-frequency QPOs. We performed a number of simulation runs for different magnetospheric sizes from small to tiny, and observed a definite correlation between the inner disc radius and the QPO frequency: the frequency is higher when the magnetosphere is smaller. In the stable regime, a small magnetosphere forms and accretion through the usual funnel streams is observed, and the frequency of the star is expected to dominate the light curve. We performed exploratory investigations of the case in which the magnetosphere becomes negligibly small and the disc interacts with the star through an equatorial belt. We also performed investigation of somewhat larger magnetospheres where one or two ordered tongues may dominate over other chaotic tongues. In application to millisecond pulsars, we obtain QPO frequencies in the range of 350–990 Hz for one spot. The frequency associated with rotation of one spot may dominate if spots are not identical or antipodal. If the spots are similar and antipodal, then the frequencies are twice as high. We show that variation of the accretion rate leads to drift of the QPO peak.     

X-ray and ultraviolet observations of the dwarf nova VW Hyi in quiescence

We present an analysis of X-ray and ultraviolet (UV) data of the dwarf nova VW Hyi that were obtained with XMM–Newton during the quiescent state. The X-ray spectrum indicates the presence of an optically thin plasma in the boundary layer that cools as it settles on to the white dwarf. The plasma has a continuous temperature distribution that is well described by a power law or a cooling flow model with a maximum temperature of 6–8 keV. We estimate from the X-ray spectrum a boundary layer luminosity of  8 × 10³⁰ erg s^-1  , which is only 20 per cent of the disc luminosity. The rate of accretion on to the white dwarf is  5 × 10⁻¹² M_⊙ yr⁻¹  , about half of the rate in the disc. From the high-resolution X-ray spectra, we estimate that the X-ray emitting part of the boundary layer is rotating with a velocity of 540 km s⁻¹, which is close to the rotation velocity of the white dwarf but is significantly smaller than the Keplerian velocity. We detect a 60-s quasi-periodic oscillation of the X-ray flux, which is likely to be due to the rotation of the boundary layer. The X-ray and the UV flux show strong variability on a time-scale of ∼1500 s. We find that the variability in the two bands is correlated and that the X-ray fluctuations are delayed by ∼100 s. The correlation indicates that the variable UV flux is emitted near the transition region between the disc and the boundary layer and that accretion rate fluctuations in this region are propagated to the X-ray emitting part of the boundary layer within ∼100 s. An orbital modulation of the X-ray flux suggests that the inner accretion disc is tilted with respect to the orbital plane. The elemental abundances in the boundary layer are close to their solar values.     

An accretion disc model with a magnetic wind and turbulent viscosity

A model is presented for an accretion disc with turbulent viscosity and a magnetically influenced wind. The magnetic field is generated by a dynamo in the disc, involving the turbulence and radial shear. Disc-wind solutions are found for which the wind mass flux is sufficient to play a major part in driving an imposed steady inflow, but small enough for most material to be accreted on to the central object. Constraints arise for the magnetic Reynolds and Prandtl numbers in terms of the turbulent Mach number and vertical length-scale of the disc's horizontal magnetic field. It is shown that the imposition of a stellar boundary condition enhances the wind mass flux in the very inner region of the disc and may result in jet formation.     

On the effects of the stellar magnetic field on the structure of T Tauri accretion discs

The structure of accretion discs around magnetic T Tauri stars is calculated numerically using a particle hydrodynamical code, in which magnetic interaction is included in the framework of King's diamagnetic blob accretion model. Setting up the calculation so as to simulate the density structure of a quasi-steady disc in the equatorial plane of a T Tauri star, we find that the central star's magnetic field typically produces a central hole in the disc and spreads out the surface density distribution. We argue that this result suggets a promising mechanism for explaining the unusual flatness (IR excess) of T Tauri accretion disc spectra.     

The radial disc structure around a magnetic neutron star: analytic and semi-analytic solutions

The radial structure of a thin accretion disc is calculated in the presence of a central dipole magnetic field aligned with the rotation axis. The problem is treated using a modified expression for the turbulent magnetic diffusion, which allows the angular momentum equation to be integrated analytically. The governing algebraic equations are solved iteratively between 1 and 10⁴ stellar radii. An analytic approximation is provided that is valid near the disruption radius at about 100 stellar radii. At that point, which is approximately 60 per cent of the Alfvén radius and typically about 30 per cent of the corotation radius, the disc becomes viscously unstable. This instability results from the fact that both radiation pressure and opacity caused by electron scattering become important. This in turn is a consequence of the magnetic field which leads to an enhanced temperature in the inner parts. This is because the magnetic field gives rise to a strongly enhanced vertically integrated viscosity, so that the viscous torque can balance the magnetic torque.     

Mass loss from the inner regions of accretion discs due to centrifugally driven magnetic wind flows

Wind flows and collimated jets are believed to be a feature of a range of disc accreting systems. These include active galactic nuclei, T Tauri stars, X-ray binaries and cataclysmic variables. The observed collimation implies large-scale magnetic fields and it is known that dipole-symmetry fields of sufficient strength can channel wind flows emanating from the surfaces of a disc. The disc inflow leads to the bending of the poloidal magnetic field lines, and centrifugally driven magnetic winds can be launched when the bending exceeds a critical value. Such winds can result in angular momentum transport at least as effective as turbulent viscosity, and hence they can play a major part in driving the disc inflow.
It is shown here that if the standard boundary condition of vanishing viscous stress close to the stellar surface is applied, together with the standard connection between viscosity and magnetic diffusivity, then poloidal magnetic field bending increases as the star is approached with a corresponding increase in the wind mass loss rate. A significant amount of material can be lost from the system via the enhanced wind from a narrow region close to the stellar surface. This occurs for a Keplerian angular velocity distribution and for a modified form of angular velocity, which allows for matching of the disc and stellar rotation rates through a boundary layer above the stellar surface. The enhanced mass loss is significantly affected by the behaviour of the disc angular velocity as the stellar surface is approached, and hence by the stellar rotation rate. Such a mechanism may be related to the production of jets from the inner regions of disc accreting systems.     

Disc structure around strongly magnetic accretors: a full disc solution with buoyancy diffusivity

A full numerical solution is found for the effect of a strongly magnetic star on its accretion disc, for the case of magnetic buoyancy diffusion. As in the previously considered case of turbulent diffusion, the disc becomes disrupted when magnetic and viscous stresses become comparable. A magnetically induced temperature elevation leads to electron scattering opacity and radiation pressure becoming significant far from the stellar surface, with consequent viscous instability and vertical disruption of the disc. This, together with the previous turbulent case, suggests that such a disruption mechanism owing to strongly magnetic accretors is generally operable.     

The X-ray and extreme-ultraviolet flux evolution of SS Cygni throughout outburst

We present the most complete multiwavelength coverage of any dwarf nova outburst: simultaneous optical, Extreme Ultraviolet Explorer and Rossi X-ray Timing Explorer observations of SS Cygni throughout a narrow asymmetric outburst. Our data show that the high-energy outburst begins in the X-ray waveband 0.9–1.4 d after the beginning of the optical rise and 0.6 d before the extreme-ultraviolet rise. The X-ray flux drops suddenly, immediately before the extreme-ultraviolet flux rise, supporting the view that both components arise in the boundary layer between the accretion disc and white dwarf surface. The early rise of the X-ray flux shows that the propagation time of the outburst heating wave may have been previously overestimated.
The transitions between X-ray and extreme-ultraviolet dominated emission are accompanied by intense variability in the X-ray flux, with time-scales of minutes. As detailed by Mauche & Robinson, dwarf nova oscillations are detected throughout the extreme-ultraviolet outburst, but we find they are absent from the X-ray light curve.
X-ray and extreme-ultraviolet luminosities imply accretion rates of  3 × 10¹⁵ g s⁻¹  in quiescence,  1 × 10¹⁶ g s⁻¹  when the boundary layer becomes optically thick, and  ∼10¹⁸ g s⁻¹  at the peak of the outburst. The quiescent accretion rate is two and a half orders of magnitude higher than predicted by the standard disc instability model, and we suggest this may be because the inner accretion disc in SS Cyg is in a permanent outburst state.     

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.     

Stellar dynamo driven wind braking versus disc coupling

Star‐disc coupling is considered in numerical models where the stellar field is not an imposed perfect dipole, but instead a more irregular self‐adjusting dynamo‐generated field. Using axisymmetric simulations of the hydromagnetic mean‐field equations, it is shown that the resulting stellar field configuration is more complex, but significantly better suited for driving a stellar wind. In agreement with recent findings by a number of people, star‐disc coupling is less efficient in braking the star than previously thought. Moreover, stellar wind braking becomes equally important. In contrast to a perfect stellar dipole field, dynamo‐generated stellar fields favor field‐aligned accretion with considerably higher velocity at low latitudes, where the field is weaker and originating in the disc. Accretion is no longer nearly periodic (as it is in the case of a stellar dipole), but it is more irregular and episodic. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Constraining jet/disc geometry and radiative processes in stellar black holes XTE J1118+480 and GX 339−4

We present results from modelling of quasi-simultaneous broad-band (radio through X-ray) observations of the Galactic stellar black hole (BH) transient X-ray binary (XRB) systems XTE J1118+480 and GX 339−4 using an irradiated disc + compact jet model. In addition to quantifying the physical properties of the jet, we have developed a new irradiated disc model which also constrains the geometry and temperature of the outer accretion disc by assuming a disc heated by viscous energy release and X-ray irradiation from the inner regions. For the source XTE J1118+480, which has better spectral coverage of the two in optical and near-infrared (OIR) wavelengths, we show that the entire broad-band continuum can be well described by an outflow-dominated model + an irradiated disc. The best-fitting radius of the outer edge of the disc is consistent with the Roche lobe geometry of the system, and the temperature of the outer edge of the accretion disc is similar to those found for other XRBs. Irradiation of the disc by the jet is found to be negligible for this source. For GX 339−4, the entire continuum is well described by the jet-dominated model only, with no disc component required. For the two XRBs, which have very different physical and orbital parameters and were in different accretion states during the observations, the sizes of the jet base are similar and both seem to prefer a high fraction of non-thermal electrons in the acceleration/shock region and a magnetically dominated plasma in the jet. These results, along with recent similar results from modelling other galactic XRBs and AGNs, may suggest an inherent unity in diversity in the geometric and radiative properties of compact jets from accreting black holes.     

Disc structure around strongly magnetic accretors: a full disc solution with turbulent diffusivity

The problem of the effect of a strongly magnetic star on a surrounding accretion disc is considered. For stellar rotation periods greater than a critical value, a numerical solution is found for a steady disc with turbulent magnetic diffusion, including electron scattering opacity and radiation pressure. Inside the corotation radius, the extraction of disc angular momentum by magnetic coupling to the star becomes strong and this leads to enhanced viscous stress and dissipation. The resulting elevated temperature causes electron scattering opacity and radiation pressure to become significant further from the star than in the absence of its magnetic field. The disc ends as its height increases rapidly due to the large central pressure, its density decreases and magnetically induced viscous instability occurs.     

S 235 B explained: an accreting Herbig Be star surrounded by reflection nebulosity

The intent of this study is to determine the nature of the star and associated nebulosity S 235 B, which are located in a region of active star formation still heavily obscured by the parent molecular cloud. Low-resolution  ( R = 400)  long-slit spectra of the star and nebulosity, and medium-  ( R = 1800)  and high-resolution  ( R = 60 000)  spectra of the central star are presented along with the results of Fabry–Perot interferometric imaging of the entire region. Based on the long-slit and Fabry–Perot observations, the nebulosity appears to be entirely reflective in nature, with the stellar component S 235 B★ providing most of the illuminating flux. The stellar source itself is classified here as a B1V star, with emission-line profiles indicative of an accretion disc. S 235 B★ thus belongs to the relatively rare class of early-type Hebrig Be stars. Based on the intensity of the reflected component, it is concluded that the accretion disc must be viewed nearly edge-on. Estimates of the accretion rate of S 235 B★ from the width of the Hα profile at 10 per cent of maximum intensity, a method which has been used lately for T Tauri stars and Brown Dwarfs, appear to be inconsistent with the mass outflow rate and accretion rate implied from previous infrared observations by Felli et al., suggesting this empirical law does not extend to higher masses.     

Electromagnetic fields in jets

The magnetic fields and energy flows in an astronomical jet described by our earlier model are calculated in detail. Though the field distribution varies with the external pressure function   p ( z )  , it depends only weakly on the other boundary conditions. Individual field lines were plotted; the lines become nearly vertical at the bottom and are twisted at the top. An animation of a field line's motion was made, which shows the line being wound up by the accretion disc's differential rotation and rising as a result of this. The distribution of Poynting flux within the jet indicates that much of the energy flows up the jet from the inside of the accretion disc but a substantial fraction flows back down to the outside.