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.     

The anticorrelation between the hard X-ray photon index and the Eddington ratio in low-luminosity active galactic nuclei

We find a significant anticorrelation between the hard X-ray photon index Γ and the Eddington ratio   L _bol/ L _Edd  for a sample of low-ionization nuclear emission-line regions and local Seyfert galaxies, compiled from literature with Chandra or XMM–Newton observations. This result is in contrast with the positive correlation found in luminous active galactic nuclei (AGN), while it is similar to that of X-ray binaries (XRBs) in the low/hard state. Our result is qualitatively consistent with the spectra produced from advection-dominated accretion flows (ADAFs). It implies that the X-ray emission of low-luminosity active galactic nuclei (LLAGN) may originate from the Comptonization process in ADAF, and the accretion process in LLAGN may be similar to that of XRBs in the low/hard state, which is different from that in luminous AGN.     

Revisiting the infrared spectra of active galactic nuclei with a new torus emission model

We describe improved modelling of the emission by dust in a toroidal-like structure heated by a central illuminating source within active galactic nuclei (AGNs). We have chosen a simple but realistic torus geometry, a flared disc, and a dust grain distribution function including a full range of grain sizes. The optical depth within the torus is computed in detail taking into account the different sublimation temperatures of the silicate and graphite grains, which solves previously reported inconsistencies in the silicate emission feature in type 1 AGNs. We exploit this model to study the spectral energy distributions (SEDs) of 58 extragalactic (both type 1 and type 2) sources using archival optical and infrared data. We find that both AGN and starburst contributions are often required to reproduce the observed SEDs, although in a few cases they are very well fitted by a pure AGN component. The AGN contribution to the far-infrared luminosity is found to be higher in type 1 sources, with all the type 2 requiring a substantial contribution from a circumnuclear starburst. Our results appear in agreement with the AGN unified scheme, because the distributions of key parameters of the torus models turn out to be compatible for type 1 and type 2 AGNs. Further support to the unification concept comes from comparison with medium-resolution infrared spectra of type 1 AGNs by the Spitzer observatory, showing evidence for a moderate silicate emission around 10 μm, which our code reproduces. From our analysis we infer accretion flows in the inner nucleus of local AGNs characterized by high equatorial optical depths  ( A_V ≃ 100)  , moderate sizes  ( R _max < 100 pc)  and very high covering factors (   f ≃ 80  per cent) on average.     

Variability associated with α accretion disc theory for standard and advection-dominated discs

The α turbulent viscosity formalism for accretion discs must be interpreted as a mean field theory, modelling a steady state only on spatial or time-scales greater than those of the turbulence. The extent of the scale separation determines the relative precision error (RPE) of the predicted luminosity L _ν. Turbulence and the use of α implies that (1) field line stretching gives a magnetic pressure  α²/6 of the total pressure generally, and a one-to-one relation between α and the pressure ratio for thin discs, and (2) large turbulent scales in advection-dominated accretion flows (ADAFs) predict a lower L _ν precision than thin discs for a given observation duration and central mass. The allowed variability (or RPE) at frequency ν increases with the size of the contributing region. For X-ray binary ADAFs, the RPE ∼ 5 per cent at R  ≤ 1000 Schwarzchild radii ( R _s) for averages over  1000 s. However, current data for galaxies like NGC 4258 and M87 give RPEs in L _ν of 50–100 per cent even at R  ≤ 100  R _S. More data are required, but systematic deviations from ADAF predictions are more significant than random deviations, and may constrain properties of the turbulence, the accretion mode, the assumption of a steady state or the accretion rate.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.     

3D simulations of rotationally-induced line variability from a classical T Tauri star with a misaligned magnetic dipole

We present 3D simulations of rotationally induced line variability arising from complex circumstellar environment of classical T Tauri stars (CTTS) using the results of the 3D magnetohydrodynamics (MHD) simulations of Romanova et al., who considered accretion on to a CTTS with a misaligned dipole magnetic axis with respect to the rotational axis. The density, velocity and temperature structures of the MHD simulations are mapped on to the radiative transfer grid, and corresponding line source function and the observed profiles of neutral hydrogen lines (Hβ, Paβ and Brγ) are computed using the Sobolev escape probability method. We study the dependency of line variability on inclination angles ( i ) and magnetic axis misalignment angles (Θ). We find the line profiles are relatively insensitive to the details of the temperature structure of accretion funnels, but are influenced more by the mean temperature of the flow and its geometry. By comparing our models with the Paβ profiles of 42 CTTS observed by Folha & Emerson, we find that models with a smaller misaligngment angle  (Θ < ∼15°)  are more consistent with the observations which show that majority of Paβ are rather symmetric around the line centre. For a high inclination system with a small dipole misalignment angle  (Θ≈ 15°)  , only one accretion funnel (on the upper hemisphere) is visible to an observer at any given rotational phase. This can cause an anticorrelation of the line equivalent to the width in the blue wing  ( v < 0)  and that in the red wing  ( v > 0)  over half of a rotational period, and a positive correlation over the other half. We find a good overall agreement of the line variability behaviour predicted by our model and those from observations.     

Self-similar solution of hot accretion flows with ordered magnetic field and outflow

Observations and numerical magnetohydrodynamic (MHD) simulations indicate the existence of outflows and ordered large-scale magnetic fields in the inner region of hot accretion flows. In this paper, we present the self-similar solutions for advection-dominated accretion flows (ADAFs) with outflows and ordered magnetic fields. Stimulated by numerical simulations, we assume that the magnetic field has a strong toroidal component and a vertical component in addition to a stochastic component. We obtain the self-similar solutions to the equations describing the magnetized ADAFs, taking into account the dynamical effects of the outflow. We compare the results with the canonical ADAFs and find that the dynamical properties of ADAFs such as radial velocity, angular velocity and temperature can be significantly changed in the presence of ordered magnetic fields and outflows. The stronger the magnetic field is, the lower the temperature of the accretion flow will be and the faster the flow rotates. The relevance to observations is briefly discussed.     

Inflow and outflow from the accretion disc of the microquasar SS 433: UKIRT spectroscopy

A succession of near-infrared (near-IR) spectroscopic observations, taken nightly throughout an entire cycle of SS 433's orbit, reveal (i) the persistent signature of SS 433's accretion disc, having a rotation speed of  ∼500 km s⁻¹  , (ii) the presence of circumbinary disc recently discovered at optical wavelengths by Blundell, Bowler & Schmidtobreick (2008) and (iii) a much faster outflow than has previously been measured for the disc wind, with a terminal velocity of  ∼1500 km s⁻¹  . The increased wind terminal velocity results in a mass-loss rate of  ∼10⁻⁴ M_⊙ yr⁻¹  . These, together with the newly (upwardly) determined masses for the components of the SS 433 system, result in an accurate diagnosis of the extent to which SS 433 has super-Eddington flows. Our observations imply that the size of the companion star is comparable with the semiminor axis of the orbit which is given by     , where e is the eccentricity. Our relatively spectral resolution at these near-IR wavelengths has enabled us to deconstruct the different components that comprise the Brackett-γ (Brγ) line in this binary system, and their physical origins. With this line being dominated throughout our series of observations by the disc wind, and the accretion disc itself being only a minority (∼15 per cent) contribution, we caution against use of the unresolved Brγ line intensity as an 'accretion signature' in X-ray binaries or microquasars in any quantitative way.     

Hot accretion discs with thermal Comptonization and advection in luminous black hole sources

We solve for the structure of a hot accretion disc with unsaturated thermal Comptonization of soft photons and with advection, generalizing the classical model of Shapiro et al. The upper limit on the accretion rate due to advection constrains the luminosity to ≲ 0.15 y^3/5 α^7/5 of the Eddington limit, where y and α are the Compton and viscosity parameters, respectively. The characteristic electron temperature and Thomson optical depth of the inner flow at accretion rates within an order of magnitude of that upper limit are ∼ 10⁹ K and ∼ 1, respectively. The resulting spectra are then in close agreement with the X-ray and soft γ-ray spectra from black hole binaries in the hard state and Seyferts. At low accretion rates, bremsstrahlung becomes the dominant radiative process.     

Magnetohydrodynamical non-radiative accretion flows in two dimensions

We present the results of axisymmetric, time-dependent magnetohydrodynamic simulations of accretion flows around black holes. The calculations begin from a rotationally supported thick torus which contains a weak poloidal field. Accretion is produced by growth and saturation of the magnetorotational instability (MRI) provided that the wavelength of the fastest growing mode is less than the thickness of the torus. Using a computational grid that spans more than two decades in radius, we compare the time-averaged properties of the flow with previous hydrodynamical simulations. The net mass accretion rate is small compared with the mass inflow and outflow rates at large radii associated with turbulent eddies. Turbulence is driven by the MRI rather than convection. The two-dimensional structure of the time-averaged flow is significantly different compared with the hydrodynamical case. We discuss the limitations imposed on our results by the assumption of axisymmetry and the relatively small radial domain.     

Radiation mechanisms and geometry of Cygnus X-1 in the soft state

We present X-ray/ γ -ray spectra of Cyg X-1 observed during the transition from the hard to the soft state and in the soft state by ASCA , RXTE and CGRO /OSSE in 1996 May and June. The spectra consist of a dominant soft component below ∼2 keV and a power-law-like continuum extending to at least ∼800 keV. We interpret them as emission from an optically thick, cold accretion disc and from an optically thin, non-thermal corona above the disc. A fraction f ≳0.5 of total available power is dissipated in the corona.
We model the soft component by multicolour blackbody disc emission taking into account the torque-free inner-boundary condition. If the disc extends down to the minimum stable orbit, the ASCA RXTE data yield the most probable black hole mass of M _X≈10 M_⊙ and an accretion rate,     , locating Cyg X-1 in the soft state in the upper part of the stable, gas-pressure-dominated, accretion-disc solution branch.
The spectrum of the corona is well modelled by repeated Compton scattering of seed photons from the disc off electrons with a hybrid, thermal/non-thermal distribution. The electron distribution can be characterized by a Maxwellian with an equilibrium temperature of kT _e∼30–50 keV, a Thomson optical depth of τ ∼0.3 and a quasi-power-law tail. The compactness of the corona is 2≲ℓ_h≲7, and a presence of a significant population of electron–positron pairs is ruled out.
We find strong signatures of Compton reflection from a cold and ionized medium, presumably an accretion disc, with an apparent reflector solid angle, Ω/2π∼0.5–0.7. The reflected continuum is accompanied by a broad iron K α line.     

X-ray emissions from two-temperature accretion flows within a dipole magnetic funnel

We investigate the hydrodynamics of accretion channelled by a dipolar magnetic field (funnel flows). We consider situations in which the electrons and ions in the flow cannot maintain thermal equilibrium [two-temperature (2T) effects] due to strong radiative loss, and determine the effects on the keV X-ray properties of the systems. We apply this model to investigate the accretion shocks of white dwarfs in magnetic cataclysmic variables (mCVs). We have found that the incorporation of 2T effects could harden the keV X-rays. Also, the dipolar model yields harder X-ray spectra than the standard planar model if white dwarf is sufficiently massive  (≳1 M_⊙)  . When fitting observed keV X-ray spectra of mCVs, the inclusion of 2T hydrodynamics and a dipolar accretion geometry lowers estimates for white dwarf masses when compared with masses inferred from models excluding these effects. We find mass reductions ≲9 per cent in the most massive cases.     

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.     

Accretion in stellar clusters and the initial mass function

We present a simple physical mechanism that can account for the observed stellar mass spectrum for masses M ∗≳0.5 M_⊙ . The model depends solely on the competitive accretion that occurs in stellar clusters where each star's accretion rate depends on the local gas density and the square of the accretion radius. In a stellar cluster, there are two different regimes depending on whether the gas or the stars dominate the gravitational potential. When the cluster is dominated by cold gas, the accretion radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a ρ  ∝  R ⁻² distribution. Accretion in this regime results in a mass spectrum with an asymptotic limit of γ =−3/2 (where Salpeter is γ =−2.35) . Once the stars dominate the potential and are virialized, which occurs first in the cluster core, the accretion radius is the Bondi–Hoyle radius. The resultant mass spectrum has an asymptotic limit of γ =−2 with slightly steeper slopes ( γ ≈−2.5) if the stars are already mass-segregated. Simulations of accretion on to clusters containing 1000 stars show that, as expected, the low-mass stars accumulate the majority of their masses during the gas-dominated phase whereas the high-mass stars accumulate the majority of their masses during the stellar-dominated phase. This results in a mass spectrum with a relatively shallow γ ≈3/2 power law for low-mass stars and a steeper power law for high-mass stars −2.5≲ γ ≤−2 . This competitive accretion model also results in a mass-segregated cluster.     

The terminal bulk Lorentz factor of relativistic electron–positron jets

We present a numerical simulation of the bulk Lorentz factor of a relativistic electron–positron jet driven by the Compton rocket effect from accretion disc radiation. The plasma is assumed to have a power-law distribution n _e(γ) ∝ γ^{− s} with 1 < γ < γ_max and is continuously reheated to compensate for radiation losses. We include the full Klein–Nishina (hereafter KN) cross-section, and study the role of the energy upper cut-off γ_max, spectral index s and source compactness. We determine the terminal bulk Lorentz factor in the cases of supermassive black holes, relevant to AGN, and stellar black holes, relevant to galactic microquasars. In the latter case, Klein–Nishina cross-section effects are more important and induce a terminal bulk Lorentz factor smaller than in the former case. Our result are in good agreement with bulk Lorentz factors observed in Galactic (GRS 1915+105, GRO J1655−40) and extragalactic sources. Differences in scattered radiation and acceleration mechanism efficiency in the AGN environment can be responsible for the variety of relativistic motion in those objects. We also take into account the influence of the size of the accretion disc; if the external radius is small enough, the bulk Lorentz factor can be as high as 60.     

Testing the locality of transport in self-gravitating accretion discs – II. The massive disc case

We present hydrogen emission-line profile models of magnetospheric accretion on to classical T Tauri stars. The models are computed under the Sobolev approximation using the three-dimensional Monte Carlo radiative-transfer code torus . We have calculated four illustrative models in which the accretion flows are confined to azimuthal curtains – a geometry predicted by magnetohydrodynamic simulations. Properties of the line profile variability of our models are discussed, with reference to dynamic spectra and cross-correlation images. We find that some gross characteristics of the observed line profile variability are reproduced by our models, although in general the level of variability predicted is larger than that observed. We conclude that this excessive variability probably excludes dynamical simulations that predict accretion flows with low degrees of axisymmetry.     

An unconventional accretion disc dynamo

In the light of recent results from numerical simulations of accretion disc MHD turbulence, we revisit the problem of the configuration of large-scale magnetic fields resulting from an α Ω dynamo operating in a thin accretion disc. In particular, we analyse the consequences of the peculiar sign of the α -effect suggested by numerical simulations . We determine the symmetry of the fastest-growing modes in the kinematic dynamo approximation and, in the framework of an ' α -quenched' dynamo model, study the evolution of the magnetic field. We find that the resulting field for this negative polarity of the α -effect generally has dipole symmetry with respect to the disc midplane, although the existence of an equilibrium configuration depends on the properties of the turbulence. The role of magnetic field dragging is discussed and, finally, the presence of an external uniform magnetic field is included to address the issue of magneto centrifugal wind launching from accretion discs.     

Resistive relaxation of a magnetically confined mountain on an accreting neutron star

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to  τ_I∼ 10⁵–10⁸ yr  (depending on the conductivity of the neutron star crust) for an accreted mass of   M _a= 1.2 × 10⁻⁴ M_⊙  . The magnetic dipole moment simultaneously re-emerges as the screening currents dissipate over  τ_I  . For non-axisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfvén time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease  τ_I  by one order of magnitude relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.     

Hydrodynamical non-radiative accretion flows in two dimensions

Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows that cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus with constant specific angular momentum. Accretion is induced by the addition of a small anomalous azimuthal shear stress which is characterized by a function ν . We study the flows generated as the amplitude and form of ν are varied. A spherical polar grid which spans more than two orders of magnitude in radius is used to resolve the flow over a wide range of spatial scales. We find that convection in the inner regions produces significant outward mass motions that carry away both the energy liberated by and a large fraction of the mass participating in the accretion flow. Although the instantaneous structure of the flow is complex and dominated by convective eddies, long-time averages of the dynamical variables show remarkable correspondence to certain steady-state solutions. The two-dimensional structure of the time-averaged flow is marginally stable to the Høiland criterion, indicating that convection is efficient. Near the equatorial plane, the radial profiles of the time-averaged variables are power laws with an index that depends on the radial scaling of the shear stress. A stress in which ν ∝ r ^1/2 recovers the widely studied self-similar solution corresponding to an ' α -disc'. We find that, regardless of the adiabatic index of the gas, or the form or magnitude of the shear stress, the mass inflow rate is a strongly increasing function of radius, and is everywhere nearly exactly balanced by mass outflow. The net mass accretion rate through the disc is only a fraction of the rate at which mass is supplied to the inflow at large radii, and is given by the local, viscous accretion rate associated with the flow properties near the central object.