Stellar dynamics in a galactic centre surrounded by a massive accretion disc — I. Newtonian description

The long-term evolution of stellar orbits bound to a massive centre is studied in order to understand the cores of star clusters in central regions of galaxies. Stellar trajectories undergo tiny perturbations, the origins of which are twofold: (i) the gravitational field of a thin gaseous disc surrounding the galactic centre, and (ii) cumulative drag arising from successive interactions of the stars with the material of the disc. Both effects are closely related because they depend on the total mass of the disc, assumed to be a small fraction of the central mass. It is shown that, in contrast to previous works, most of the retrograde (with respect to the disc) orbits are captured by the central object, presumably a massive black hole. Initially prograde orbits are also affected, so that statistical properties of the central star cluster in quasi-equilibrium may differ significantly from those deduced in previous analyses.     

Can gas dynamics in centres of galaxies reveal orbiting massive black holes?

If supermassive black holes in centres of galaxies form by merging of black hole remnants of massive Population III stars, then there should be a few black holes of mass one or two orders of magnitude smaller than that of the central ones, orbiting around the centre of a typical galaxy. These black holes constitute a weak perturbation in the gravitational potential, which can generate wave phenomena in gas within a disc close to the centre of the galaxy. Here, we show that a single orbiting black hole generates a three-arm spiral pattern in the central gaseous disc. The density excess in the spiral arms in the disc reaches values of 3–12 per cent when the orbiting black hole is about 10 times less massive than the central black hole. Therefore, the observed density pattern in gas can be used as a signature in detecting the most massive orbiting black holes.     

The light curves of soft X-ray transients

We show that the light curves of soft X-ray transients (SXTs) follow naturally from the disc instability picture, adapted to take account of irradiation by the central X-ray source during the outburst. Irradiation prevents the disc from returning to the cool state until central accretion is greatly reduced. This happens only after most of the disc mass has been accreted by the central object, on a viscous time-scale, accounting naturally for the exponential decay of the outburst on a far longer time-scale (τ20–40 d) than seen in dwarf novae, without any need to manipulate the viscosity parameter α. The accretion of most of the disc mass in outburst explains the much longer recurrence time of SXTs compared with dwarf novae. This picture also suggests an explanation of the secondary maximum seen in SXT light curves about 50–75 d after the start of each outburst, since central irradiation triggers the thermal instability of the outer disc, adding to the central accretion rate one viscous time later. The X-ray outburst decay constant τ should on average increase with orbital period, but saturate at a roughly constant value ∼40 d for orbital periods longer than about a day. The bolometric light curve should show a linear rather than an exponential decay at late times (a few times τ). Outbursts of long-period systems should be entirely in the linear decay regime, as is observed in GRO J1744−28. UV and optical light curves should resemble the X-rays but have decay time-scales up to 2–4 times longer.     

A massive thick disc around a massive black hole and its runaway instability

We investigate the runaway instability of configurations consisting of a massive dense but non-self-gravitating thick disc/torus surrounding a massive black hole (MBH). We limit our model parameters to values that result in a self-consistent thick disc around an MBH. We identify, analytically, the index of the angular momentum distribution that will form a thick disc as an initial equilibrium state, and obtain the mass ratio of the disc to the central black hole for which the disc is dominated by the radiation pressure. In our theoretical framework we find that a self-consistent thick disc with constant angular momentum leads to a runaway instability on a dynamical time-scale. However, even a slight increase of the specific angular momentum outwards has a strong stabilizing effect on the accretion process. Finally, we discuss our results and present possible applications to high-energy emission.     

The impact of a supernova explosion in a very massive binary

We consider the effect of a supernova (SN) explosion in a very massive binary that is expected to form in a portion of Population III stars with the mass higher than  100 M_⊙  . In a Population III binary system, a more massive star can result in the formation of a black hole (BH) and a surrounding accretion disc. Such BH accretion could be a significant source of the cosmic reionization in the early Universe. However, a less massive companion star evolves belatedly and eventually undergoes a SN explosion, so that the accretion disc around a BH might be blown off in a lifetime of companion star. In this paper, we explore the dynamical impact of a SN explosion on an accretion disc around a massive BH, and elucidate whether the BH accretion disc is totally demolished or not. For the purpose, we perform three-dimensional hydrodynamic simulations of a very massive binary system, where we assume a BH of  10³ M_⊙  that results from a direct collapse of a very massive star and a companion star of  100 M_⊙  that undergoes a SN explosion. We calculate the remaining mass of a BH accretion disc as a function of time. As a result, it is found that a significant portion of gas disc can survive through three-dimensional geometrical effects even after the SN explosion of a companion star. Even if the SN explosion energy is higher by two orders of magnitude than the binding energy of gas disc, about a half of disc can be left over. The results imply that the Population III BH accretion disc can be a long-lived luminous source, and therefore could be an important ionizing source in the early Universe.     

Gravitational instability of discs with dissipative coronae around supermassive black holes

We study the dynamical structure of a self-gravitating disc with coronae around a supermassive black hole. Assuming that the magnetorotational instability responsible for generating the turbulent stresses inside the disc is also the source for a magnetically dominated corona, a fraction of the power released when the disc matter accretes is transported to and dissipated in the corona. This has a major effect on the structure of the disc and its gravitational (in)stability according to our analytical and self-consistent solutions. We determine the radius where the disc crosses the inner radius of gravitational instability and forms the first stars. Not only the location of this radius which may extend to very large distances from the central black hole, but also the mass of the first stars highly depends on the input parameters, notably the viscosity coefficient, the mass of the central object and the accretion rate. For accretion discs around quasi-stellar objects (QSOs) and the Galactic Centre, we determine the self-gravitating radius and the mass of the first clumps. Comparing the cases with a corona and without a corona for typical discs around QSOs or the Galactic Centre, when the viscosity coefficient is around 0.3, we show that the self-gravitating radius decreases by a factor of approximately 2, but the mass of the fragments increases with more or less the same factor. The existence of a corona implies a more gravitationally unstable disc according to our results. The effect of a corona on the instability of the disc is more effective when the viscosity coefficient increases.     

Progress on a model for β Lyrae

It is generally accepted that β Lyrae is a massive mass-transfer system in which the more massive gainer is embedded in a thick accretion disc. No detailed, quantitative accretion disc model of the system, which fits both the available photometric and spectroscopic data, has been published.
This paper considers a thick accretion disc model. An empirical two-temperature model for the accretion disc rim produces an accurate fit to the UBV and OAO2 photometry. The two-temperature rim model leads to synthetic spectra with a much larger Balmer jump than is present in spectrophotometric scans. The same model fails to produce IR light curves with the 0.5^p eclipse deeper than the 0.0^p eclipse, as shown in the observations of Jameson & Longmore.
The paper next reports initial results for a model in which an extra source of continuum radiation arises from Thomson scattering of radiation from the gainer, which is hidden from the observer by the thick accretion disc.
The adopted rate of mass transfer cannot supply the energy radiated by the accretion disc rim, nor can any credible larger rate. The most likely energy source is the mass gainer.     

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.     

On the wake generated by a planet in a disc

A planet of low mass orbiting in a two-dimensional gaseous disc generates a one-armed spiral wake. We explain this phenomenon as the result of constructive interference between wave modes in the disc, somewhat similar to the Kelvin wedge produced in the wake of a ship. The same feature is not expected in a three-dimensional disc with thermal stratification.     

Magnetized, mass-loaded, rotating accretion flows

We present a semi-analytical investigation of a simple one-dimensional, steady-state model for a mass-loaded, rotating, magnetized, hydrodynamical flow. Our approach is analogous to one used in early studies of magnetized winds. The model represents the infall towards a central point mass of the gas generated in a cluster of stars surrounding it, as is likely to occur in some active nuclei and starburst galaxies. We describe the properties of the different classes of infall solutions. We find that the flow becomes faster than the fast-mode speed, and hence decoupled from the centre, only for a limited range of parameter values, and when magnetic stresses are ineffective. Such flow is slowed as it approaches a centrifugal barrier, implying the existence of an accretion disc. When the flow does not become super-fast and the magnetic torque is insufficient, no steady solution extending inward to the centre exists. Finally, with a larger magnetic torque, solutions representing steady sub-Alfvénic flows are found, which can resemble spherical hydrodynamical infall. Such solutions, if applicable, would imply that rotation is not important and that any accretion disc formed would be of very limited size.     

The effect of radiation stress on gas motions in AGN accretion discs

In this paper we look at one of the effects of irradiation on a warped accretion disc in the context of active galactic nuclei (AGN). A warp will catch a substantial amount of the radiation emitted by the central object. We consider the fluid motions that may arise inside a warped disc when the surface is subject to a radiation stress, and also the net mass flows that result. We find that, to first order, we have a balance of the viscous and Coriolis-type forces. The radial radiation stress causes outward motion of the surface layer, but only the azimuthal Poynting–Robertson drag leads to an increase in the net accretion rate. We investigate the distribution of the velocity perturbations and find them to be significant in determining the local structure of the disc.
An unexpected result is that the picture changes significantly when we take into account the periodic illumination of the warped disc. A type of resonance at the local Keplerian rotation frequency causes a flow that penetrates the whole thickness of the disc; these flows are faster than the flows due to unchanging illumination. They totally dominate the induced flows in terms of sheer mass, but significant impact on disc structure still occurs only near the surface, where velocity perturbations typically go up to some kilometres per second.     

Substellar companions and isolated planetary-mass objects from protostellar disc fragmentation

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time-scale that is short compared with the orbital period. We use a combination of hydrodynamic simulations and N -body orbit integrations to study the long-term evolution of a fragmenting disc with an initial mass ratio to the star of   M _disc/ M _*= 0.1  . For a disc that is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar-mass star) up to  ≈0.01 M_⊙  . Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 – in which the companions are close to or beyond the deuterium burning limit – appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ∼1 au should differ from that of stars with lower mass planetary companions.     

The jet power extracted from a magnetized accretion disc

We consider the power of a relativistic jet accelerated by the magnetic field of an accretion disc. It is found that the power extracted from the disc is mainly determined by the field strength and configuration of the field far from the disc. Comparing it with the power extracted from a rotating black hole, we find that the jet power extracted from a disc can dominate over that from the rotating black hole. However, in some cases, the jet power extracted from a rapidly rotating hole can be more important than that from the disc, even if the poloidal field threading the hole is not significantly larger than that threading the inner edge of the disc. The results imply that the radio-loudness of quasars may be governed by its accretion rate, which might be regulated by the central black hole mass. It is proposed that the different disc field generation mechanisms might be tested against observations of radio-loud quasars if their black hole masses are available.     

Pseudo-viscous modelling of self-gravitating discs and the formation of low mass ratio binaries

We present analytic models for the local structure of self-regulated self-gravitating accretion discs that are subject to realistic cooling. Such an approach can be used to predict the secular evolution of self-gravitating discs (which can usefully be compared with future radiation hydrodynamical simulations) and to define various physical regimes as a function of radius and equivalent steady state accretion rate. We show that fragmentation is inevitable, given realistic rates of infall into the disc, once the disc extends to radii >70 au (in the case of a solar mass central object). Owing to the outward redistribution of disc material by gravitational torques, we also predict fragmentation at >70 au even in the case of low angular momentum cores which initially collapse to a much smaller radius. We point out that 70 au is close to the median binary separation and propose that such delayed fragmentation, at the point that the disc expands to >70 au, ensures the creation of low mass ratio companions that can avoid substantial further growth and consequent evolution towards unit mass ratio. We thus propose this as a promising mechanism for producing low mass ratio binaries, which, while abundant observationally, are severely underproduced in hydrodynamical models.     

Stability of accretion discs threaded by a strong magnetic field

We study the stability of poloidal magnetic fields anchored in a thin accretion disc. The two-dimensional hydrodynamics in the disc plane is followed by a grid-based numerical simulation including the vertically integrated magnetic forces. The three-dimensional magnetic field outside the disc is calculated in a potential field approximation from the magnetic flux density distribution in the disc. For uniformly rotating discs we confirm numerically the existence of the interchange instability as predicted by Spruit, Stehle & Papaloizou . In agreement with predictions from the shearing sheet model, discs with Keplerian rotation are found to be stabilized by the shear, as long as the contribution of magnetic forces to support against gravity is small. When this support becomes significant, we find a global instability which transports angular momentum outwardly and allows mass to accrete inwardly. The instability takes the form of a m =1 rotating 'crescent', reminiscent of the purely hydrodynamic non-linear instability previously found in pressure-supported discs. A model where the initial surface mass density Σ( r ) and B _z ( r ) decrease with radius as power laws shows transient mass accretion during about six orbital periods, and settles into a state with surface density and field strength decreasing approximately exponentially with radius. We argue that this instability is likely to be the main angular momentum transport mechanism in discs with a poloidal magnetic field sufficiently strong to suppress magnetic turbulence. It may be especially relevant in jet-producing discs.     

The runaway instability of thick discs around black holes I. The constant angular momentum case

We present results from a numerical study of the runaway instability of thick discs around black holes. This instability is an important issue for most models of cosmic gamma-ray bursts, where the central engine responsible for the initial energy release is such a system consisting of a thick disc surrounding a black hole. We have carried out a comprehensive number of time-dependent simulations aimed at exploring the appearance of the instability. Our study has been performed using a fully relativistic hydrodynamics code. The general relativistic hydrodynamic equations are formulated as a hyperbolic flux-conservative system and solved using a suitable Godunov-type scheme. We build a series of constant angular momentum discs around a Schwarzschild black hole. Furthermore, the self-gravity of the disc is neglected and the evolution of the central black hole is assumed to be that of a sequence of exact Schwarzschild black holes of varying mass. The black hole mass increase is thus determined by the mass accretion rate across the event horizon. In agreement with previous studies based on stationary models, we find that by allowing the mass of the black hole to grow the disc becomes unstable. Our hydrodynamical simulations show that for all disc-to-hole mass ratios considered (between 1 and 0.05), the runaway instability appears very fast on a dynamical time-scale of a few orbital periods, typically a few 10 ms and never exceeding 1 s for our particular choice of the mass of the black hole (2.5 M) and a large range of mass fluxes  ( m 10^-3 M s^-1)  . The implications of our results in the context of gamma-ray bursts are briefly discussed.     

The brown dwarf desert as a consequence of orbital migration

We show that the dearth of brown dwarfs in short-period orbits around Solar-mass stars – the brown dwarf desert – can be understood as a consequence of inward migration within an evolving protoplanetary disc. Brown dwarf secondaries forming at the same time as the primary star have masses which are comparable to the initial mass of the protoplanetary disc. Subsequent disc evolution leads to inward migration, and destruction of the brown dwarf, via merger with the star. This is in contrast with massive planets, which avoid this fate by forming at a later epoch when the disc is close to being dispersed. Within this model, a brown dwarf desert arises because the mass at the hydrogen-burning limit is coincidentally comparable to the initial disc mass for a Solar mass star. Brown dwarfs should be found in close binaries around very low mass stars, around other brown dwarfs, and around Solar-type stars during the earliest phases of star formation.     

Orbital motion of a massive object in a dense stellar system

In this paper we show the positional oscillation of a massive object in a dense stellar system by numerical N -body simulations. We found that the central massive object, which at first is placed at rest at the centre of the surrounding spherical stellar system, promptly departs from the centre and rotates in accordance with the rotation of the stellar system, if the stellar system has an appreciable rotation. This oscillatory motion continues for a long time because of the absence of dynamical friction. Such a long-lasting oscillation may explain the asymmetric structure observed in the centres of M31 and NGC 4486B, may cause the secular flow of gaseous elements distributed in the central regions of galaxies on to the massive object, and may ignite activity in the centres of galaxies.     

Planetary migration to large radii

There is evidence for the existence of massive planets at orbital radii of several hundred au from their parent stars where the time-scale for planet formation by core accretion is longer than the disc lifetime. These planets could have formed close to their star and then migrated outwards. We consider how the transfer of angular momentum by viscous disc interactions from a massive inner planet could cause significant outward migration of a smaller outer planet. We find that it is in principle possible for planets to migrate to large radii. We note, however, a number of effects which may render the process somewhat problematic.     

On the migration of a system of protoplanets

The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter-sized planets is studied numerically. The disc is assumed to be flat and non-self-gravitating; this is modelled by the planar (two-dimensional) Navier–Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass M _Jup each, and the radial distances from the star are one and two semimajor axes of Jupiter, respectively.
During the evolution a joint wide annular gap is created by the planets. Both planets increase their mass owing to accretion of gas from the disc: after about 2500 orbital periods of the inner planet it has reached a mass of 2.3  M _Jup, while the outer planet has reached a mass of 3.2  M _Jup. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc. The semimajor axis of the inner planet remains constant as there is very little gas left in its vicinity to induce any migration. When the distance of close approach eventually becomes smaller than the mutual Hill radius, the eccentricities increase strongly and the system may become unstable.
If disc depletion occurs rapidly enough before the planets come too close to each other, a stable system similar to our own Solar system may remain. Otherwise the orbits may become unstable and produce systems like υ And.