Accretion disc coronae as magnetic reservoirs

Most astrophysical sources powered by accretion on to a black hole, either of stellar mass or supermassive, when observed with hard X-rays show signs of a hot Comptonizing component in the flow, the so-called corona , with observed temperatures and optical depths lying in a narrow range (0.1≲ τ ≲1 and 1×10⁹ K≲ T ≲3×10⁹ K). Here we argue that these facts constitute strong supporting evidence for a magnetically dominated corona. We show that the inferred thermal energy content of the corona, in all black hole systems, is far too low to explain their observed hard X-ray luminosities, unless either the size of the corona is at least of the order of 10³ Schwarzschild radii, or the corona itself is in fact a reservoir , where the energy is mainly stored in the form of a magnetic field generated by a sheared rotator (probably the accretion disc). We briefly outline the main reasons why the former possibility is to be discarded, and the latter preferred.     

Global structure of self-excited magnetic fields arising from the magnetic shear instability

The global structure of a self-excited magnetic field arising from the magnetic shear instability has been simulated in spherical geometry by a 3D fully non-linear approach. In order to model the structure of an accretion disc we prescribe a rotation profile of the Brandt type which is Keplerian in the outer regions but yields rigid rotation at the inner core. We performed a whole series of runs at different dynamo numbers with an increasing number of modes in spectral space, thereby checking the influence of the numerical resolution in our simulations. Starting from arbitrary small perturbations, the magnetic and kinetic energies grow by several orders of magnitude as soon as a certain azimuthal resolution of at least m =15 was used at a dynamo number of order C =10⁵. Several phases of the transition to turbulence are realized and interpretations are given for the respective effects occurring at each stage. The resulting magnetic field is highly non-axisymmetric and possesses a pronounced inhomogeneous vortex structure of twisted flux tubes. The flow is almost axisymmetric but shows a Kolmogorov-like behaviour for small scales. The outer surface of the shell is penetrated by magnetic field lines in spot-like regions, which are located mainly in the equatorial plane. For very high dynamo numbers we find a cyclic behaviour of the averaged magnetic field amplitude. The problem of angular momentum transport is discussed in terms of the ShakuraSunyaev viscosity alpha , which depends exponentially on the radial distance and adopts values in the range 10³10⁵.     

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.     

Structures in a class of magnetized scale-free discs

We construct analytically stationary global configurations for both aligned and logarithmic spiral coplanar magnetohydrodynamics (MHD) perturbations in an axisymmetric background MHD disc with a power-law surface mass density  Σ₀∝ r ^−α  , a coplanar azimuthal magnetic field   B ₀∝ r ^−γ  , a consistent self-gravity and a power-law rotation curve   v ₀∝ r ^−β  , where v ₀ is the linear azimuthal gas rotation speed. The barotropic equation of state  Π∝Σ ⁿ   is adopted for both MHD background equilibrium and coplanar MHD perturbations where Π is the vertically integrated pressure and n is the barotropic index. For a scale-free background MHD equilibrium, a relation exists among  α, β, γ  and n such that only one parameter (e.g. β) is independent. For a linear axisymmetric stability analysis, we provide global criteria in various parameter regimes. For non-axisymmetric aligned and logarithmic spiral cases, two branches of perturbation modes (i.e. fast and slow MHD density waves) can be derived once β is specified. To complement the magnetized singular isothermal disc analysis of Lou, we extend the analysis to a wider range of  −1/4 < β < 1/2  . As an illustrative example, we discuss specifically the  β= 1/4  case when the background magnetic field is force-free. Angular momentum conservation for coplanar MHD perturbations and other relevant aspects of our approach are discussed.     

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.     

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.     

Spectral variability in transonic discs around black holes

Transonic discs with accretion rates relevant to intrinsically bright Galactic X-ray sources ( L ≈10³⁸–10³⁹ erg s⁻¹) exhibit a time-dependent cyclic behaviour due to the onset of a thermal instability driven by radiation pressure. In this paper we calculate radiation spectra emitted from thermally unstable discs to provide detailed theoretical predictions for observationally relevant quantities. The emergent spectrum has been obtained by solving self-consistently the vertical structure and radiative transfer in the disc atmosphere. We focus on four particular stages of the disc evolution, the maximal evacuation stage and three intermediate stages during the replenishment phase. The disc is found to undergo rather dramatic spectral changes during the evolution, emitting mainly in the 1–10 keV band during outburst and in the 0.1–1 keV band off-outburst. Local spectra, although different in shape from a blackbody at the disc effective temperature, may be characterized in terms of a hardening factor f . We have found that f is more or less constant, both in radius and in time, with a typical value ∼ 1.65.     

A possible energy mechanism for cosmological γ-ray bursts

We suggest that an extreme Kerr black hole with a mass ∼10⁶ M_⊙, a dimensionless angular momentum     and a marginally stable orbital radius     located in a normal galaxy may produce a γ -ray burst (GRB) by capturing and disrupting a star. During the capture period, a transient accretion disc is formed and a strong transient magnetic field ∼     lasting for     may be produced at the inner boundary of the accretion disc. A large amount of rotational energy of the black hole is extracted and released in an ultrarelativistic jet with a bulk Lorentz factor Γ larger than 10³ via the Blandford–Znajek process. The relativistic jet energy can be converted into γ -radiation via an internal shock mechanism. The GRB duration should be the same as the lifetime of the strong transient magnetic field. The maximum number of sub-bursts is estimated to be     because the disc material is likely to break into pieces with a size about the thickness of the disc h at the cusp     The shortest risetime of the burst estimated from this model is ∼     The model GRB density rate is also estimated.     

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.     

Thermal conduction between an accretion disc and a corona in active galactic nuclei: vertical structure of the transition layer

We study the vertical structure of the transition layer between an accretion disc and a corona in the context of the existence of a two-phase medium in thermally unstable regions. The disc is illuminated by hard X-ray radiation, and satisfies the condition of hydrostatic equilibrium. We take into account the energy exchange between the hot, Compton-heated corona (∼10⁸ K) and cool disc (∼10⁴ K) arising from both radiative processes and thermal conduction. In the case including thermal conduction, we perform a local stability analysis, and conclude that thermal conduction does not suppress thermal instability. In spite of the continuous temperature profile T ( τ ) there are regions of strong temperature gradient, in which spontaneous perturbations can lead to cloud condensation in the transition layer. We determine the minimum size λ _TC of such a perturbation.     

Gravitational wave emission from a companion black hole in the presence of an accretion disc around a super-massive Kerr black hole

Gravitational wave signal characteristics from a binary black hole system in which the companion moves through the accretion disc of the primary are studied. We chose the primary to be a super-massive  ( M = 10⁸ M_⊙)  Kerr black hole and the companion to be a massive black hole  ( M = 10⁵ M_⊙)  to clearly demonstrate the effects. We show that the drag exerted on the companion by the disc is sufficient to reduce the coalescence time of the binary. The drag is primarily due to the fact that the accretion disc on a black hole deviates from a Keplerian disc and becomes sub-Keplerian due to inner boundary condition on the black hole horizon. We consider two types of accretion rates on to the companion. The companion is deeply immersed inside the disc and it can accrete at the Bondi rate which depends on the instantaneous density of the disc. However, an accretion disc can also form around the smaller black hole and it can accrete at its Eddington rate. Thus, this case is also studied and the results are compared. We find that the effect of the disc will be significant in reducing the coalescence time and one needs to incorporate this while interpreting gravitational wave signals emitted from such a binary system.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.     

Super-Eddington accretion discs around Kerr black holes

We calculate the structure of the accretion disc around a rapidly rotating black hole with a super-Eddington accretion rate. The luminosity and height of the disc are reduced by the advection effect. In the case of large viscosity parameter, α>0.03, the accretion flow deviates strongly from thermodynamic equilibrium and overheats in the central region. With increasing accretion rate, the flow temperature steeply increases, reaches maximum, and then falls off. The maximum is achieved in the advection-dominated regime of accretion. The maximum temperature in the disc around a massive black hole of M =10⁸ M⊙ with α=0.3 is of order 3×10⁸ K. The discs with large accretion rates can emit X-rays in quasars as well as in galactic black hole candidates.     

Triaxial stellar systems following the r1/n luminosity law: an analytical mass–density expression, gravitational torques and the bulge/disc interplay

We have investigated the structural and dynamical properties of triaxial stellar systems whose surface brightness profiles follow the   r ^{1/ n}   luminosity law – extending the analysis by Ciotti, who explored the properties of spherical   r ^{1/ n}   systems. A new analytical expression that accurately reproduces the spatial (i.e., deprojected) luminosity density profiles (error less than 0.1 per cent) is presented for detailed modelling of the Sérsic family of luminosity profiles. We evaluate both the symmetric and the non-axisymmetric components of the gravitational potential and force, and compute the torques as a function of position. For a given triaxiality, stellar systems with smaller values of n have a greater non-axisymmetric gravitational field component . We also explore the strength of the non-axisymmetric forces produced by bulges with differing n and triaxiality on systems having a range of bulge-to-disc ratios. The increasing disc-to-bulge ratio with increasing galaxy type (decreasing n ) is found to greatly reduce the amplitude of the non-axisymmetric terms, and therefore reduce the possibility that triaxial bulges in late-type systems may be the mechanism or perturbation for non-symmetric structures in the disc.
Using seeing-convolved   r ^{1/ n}   -bulge plus exponential-disc fits to the K -band data from a sample of 80 nearby disc galaxies, we probe the relations between galaxy type, Sérsic index n and the bulge-to-disc luminosity ratio. These relations are shown to be primarily a consequence of the relation between n and the total bulge luminosity. In the K band, the trend of decreasing bulge-to-disc luminosity ratio along the spiral Hubble sequence is predominantly, though not entirely, a consequence of the change in the total bulge luminosity; the trend between the total disc luminosity and Hubble type is much weaker.     

Constraints on the ionizing flux emitted by T Tauri stars

We present the results of an analysis of ultraviolet observations of T Tauri stars (TTs). By analysing emission measures taken from the literature, we derive rates of ionizing photons from the chromospheres of five classical TTs in the range  ∼10⁴¹–10⁴⁴  photon s⁻¹, although these values are subject to large uncertainties. We propose that the He  ii /C  iv line ratio can be used as a reddening-independent indicator of the hardness of the ultraviolet spectrum emitted by TTs. By studying this line ratio in a much larger sample of objects, we find evidence for an ionizing flux which does not decrease, and may even increase, as TTs evolve. This implies that a significant fraction of the ionizing flux from TTs is not powered by the accretion of disc material on to the central object, and we discuss the significance of this result and its implications for models of disc evolution. The presence of a significant ionizing flux in the later stages of circumstellar disc evolution provides an important new constraint on disc photoevaporation models.     

Numerical simulations of type I planetary migration in non-turbulent magnetized discs

Using 2D magnetohydrodynamic (MHD) numerical simulations performed with two different finite-difference Eulerian codes, we analyse the effect that a toroidal magnetic field has on low-mass planet migration in non-turbulent protoplanetary discs. The presence of the magnetic field modifies the waves that can propagate in the disc. In agreement with a recent linear analysis, we find that two magnetic resonances develop on both sides of the planet orbit, which contribute to a significant global torque. In order to measure the torque exerted by the disc on the planet, we perform simulations in which the latter is either fixed on a circular orbit or allowed to migrate. For a     planet, when the ratio β between the square of the sound speed and that of the Alfven speed at the location of the planet is equal to 2, we find inward migration when the magnetic field   B _φ  is uniform in the disc, reduced migration when   B _φ  decreases as   r ⁻¹  and outward migration when   B _φ  decreases as   r ⁻²  . These results are in agreement with predictions from the linear analysis. Taken as a whole, our results confirm that even a subthermal stable field can stop inward migration of an earth-like planet.     

Photon bubbles in accretion discs

We show that radiation-dominated accretion discs are likely to suffer from a 'photon bubble' instability similar to that described by Arons in the context of accretion on to neutron star polar caps. The instability requires a magnetic field for its existence. In an asymptotic regime appropriate to accretion discs, we find that the overstable modes obey the remarkably simple dispersion relation
ο²=−i gkF ( B , k ).
Here g is the vertical gravitational acceleration, B is the magnetic field, and F is a geometric factor of order unity that depends on the relative orientation of the magnetic field and the wavevector. In the non-linear outcome it seems likely that the instability will enhance vertical energy transport and thereby change the structure of the innermost parts of relativistic accretion discs.     

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.     

Molecular hydrogen emission in Cygnus A

We present J , H and K -band spectroscopy of Cygnus A, spanning 1.0–2.4 μm in the rest-frame and hence several rovibrational H₂, H recombination and [Fe  ii ] emission lines. The lines are spatially extended by up to 6 kpc from the nucleus, but their distinct kinematics indicate that the three groups (H, H₂ and [Fe  ii ]) are not wholly produced in the same gas. The broadest line, [Fe  ii ] λ 1.644, exhibits a non-Gaussian profile with a broad base (FWHM≃1040 km s⁻¹), perhaps because of the interaction with the radio source. Extinctions to the line-emitting regions substantially exceed earlier measurements based on optical H recombination lines.
Hard X-rays from the quasar nucleus are likely to dominate the excitation of the H₂ emission. The results of Maloney, Hollenbach & Tielens are thus used to infer the total mass of gas in H₂ v=1–0 S(1)-emitting clouds as a function of radius, for gas densities of 10³ and 10⁵ cm⁻³, and stopping column densities N _H=10²²–10²⁴ cm⁻². Assuming azimuthal symmetry, at least 2.3×10⁸ M_⊙ of such material is present within 5 kpc of the nucleus, if the line-emitting clouds see an unobscured quasar spectrum. Alternatively, if the bulk of the X-ray absorption to the nucleus inferred by Ueno et al. actually arises in a circumnuclear torus, the implied gas mass rises to ∼10¹⁰ M_⊙. The latter plausibly accounts for 10⁹ yr of mass deposition from the cluster cooling flow, for which within this radius.