X-rays from cusps of compact remnants near galactic centres

Compact remnants – stellar mass black holes and neutron stars formed in the inner few parsec of galactic centres are predicted to sink into the central parsec due to dynamical friction on low-mass stars, forming a high concentration cusp. Same physical region may also contain very high-density molecular clouds and accretion discs that are needed to fuel supermassive black hole (SMBH) activity. Here we estimate gas capture rates on to the cusp of stellar remnants, and the resulting X-ray luminosity, as a function of the accretion disc mass. At low disc masses, most compact objects are too dim to be observable, whereas in the high disc case most of them are accreting at their Eddington rates. We find that for low accretion disc masses, compact remnant cusps may be more luminous than the central SMBHs. This 'diffuse' emission may be of importance for local moderately bright active galactic nuclei (AGNs), especially low-luminosity AGNs. We also briefly discuss how this expected emission can be used to put constraints on the black hole cusp near our Galactic Centre.     

Interface modes and their instabilities in accretion disc boundary layers

We study global non-axisymmetric oscillation modes trapped near the inner boundary of an accretion disc. Observations indicate that some of the quasi-periodic oscillations (QPOs) observed in the luminosities of accreting compact objects (neutron stars, black holes and white dwarfs) are produced in the innermost regions of accretion discs or boundary layers. Two simple models are considered in this paper. The magnetosphere–disc model consists of a thin Keplerian disc in contact with a uniformly rotating magnetosphere with and low plasma density, while the star–disc model involves a Keplerian disc terminated at the stellar atmosphere with high density and small density scaleheight. We find that the interface modes at the magnetosphere–disc boundary are generally unstable due to Rayleigh–Taylor and/or Kelvin–Helmholtz instabilities. However, differential rotation of the disc tends to suppress Rayleigh–Taylor instability, and a sufficiently high disc sound speed (or temperature) is needed to overcome this suppression and to attain net mode growth. On the other hand, Kelvin–Helmholtz instability may be active at low disc sound speeds. We also find that the interface modes trapped at the boundary between a thin disc and an unmagnetized star do not suffer Rayleigh–Taylor or Kelvin–Helmholtz instability, but can become unstable due to wave leakage to large disc radii and, for sufficiently steep disc density distributions, due to wave absorption at the corotation resonance in the disc. The non-axisymmetric interface modes studied in this paper may be relevant to the high-frequency QPOs observed in some X-ray binaries and in cataclysmic variables.     

The origin of the density distribution of disc galaxies: a new problem for the standard model of disc formation

We present new models for the formation of disc galaxies that improve upon previous models by following the detailed accretion and cooling of the baryonic mass, and by using realistic distributions of specific angular momentum. Under the assumption of detailed angular momentum conservation, the discs that form have density distributions that are more centrally concentrated than an exponential. We examine the influence of star formation, bulge formation, and feedback on the outcome of the surface brightness distributions of the stars. Low angular momentum haloes yield disc galaxies with a significant bulge component and with a stellar disc that is close to exponential, in good agreement with observations. High angular momentum haloes, on the other hand, produce stellar discs that are much more concentrated than an exponential, in clear conflict with observations. At large radii, the models reveal distinct truncation radii in both the stars and the cold gas. The stellar truncation radii result from our implementation of star formation threshold densities, and are in excellent agreement with observations. The truncation radii in the density distribution of the cold gas reflect the maximum specific angular momentum of the gas that has cooled. We find that these truncation radii occur at H  i surface densities of roughly 1 M_⊙ pc⁻², in conflict with observations. We examine various modifications to our models, including feedback, viscosity, and dark matter haloes with constant-density cores, but show that the models consistently fail to produce bulge less discs with exponential surface brightness profiles. This signals a new problem for the standard model of disc formation: if the baryonic component of the protogalaxies out of which disc galaxies form has the same angular momentum distribution as the dark matter, discs are too compact.     

From discs to planetesimals: Evolution of gas and dust discs

I review the processes that shape the evolution of protoplanetary discs around young, solar-mass stars. I first discuss observations of protoplanetary discs, and note in particular the constraints these observations place on models of disc evolution. The processes that affect the evolution of gas discs are then discussed, with the focus in particular on viscous accretion and photoevaporation, and recent models which combine the two. I then discuss the dynamics and growth of dust grains in discs, considering models of grain growth, the gas–grain interaction and planetesimal formation, and review recent research in this area. Lastly, I consider the so-called “transitional” discs, which are thought to be observed during disc dispersal. Recent observations and models of these systems are reviewed, and prospects for using statistical surveys to distinguish between the various proposed models are discussed.     

X-ray spectra of accretion discs with dynamic coronae

We compute the X-ray spectra produced by non-static coronae atop accretion discs around black holes and neutron stars. The hot corona is radiatively coupled to the underlying disc (the reflector) and generates an X-ray spectrum which is sensitive to the bulk velocity of the coronal plasma, β = v / c . We show that an outflowing corona reproduces the hard-state spectrum of Cyg X-1 and similar objects. The dynamic model predicts a correlation between the observed amplitude of reflection R and the X-ray spectrum slope Γ since both strongly depend on β . A similar correlation was observed and its shape was well fitted by the dynamic model. The scattering of soft radiation in an outflowing corona can also account for the observed optical–UV polarization pattern in active galactic nuclei.     

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.     

General relativistic spectra of accretion discs around rapidly rotating neutron stars: effect of light bending

We present computed spectra, as seen by a distant observer, from the accretion disc around a rapidly rotating neutron star. Our calculations are carried out in a fully general relativistic framework, with an exact treatment of rotation. We take into account the Doppler shift, gravitational redshift and light-bending effects in order to compute the observed spectrum. We find that light bending significantly modifies the high-energy part of the spectrum. Computed spectra for slowly rotating neutron stars are also presented. These results would be important for modelling the observed X-ray spectra of low-mass X-ray binaries containing fast-spinning neutron stars.     

Viscous overstability and eccentricity evolution in three-dimensional gaseous discs

We investigate the growth or decay rate of the fundamental mode of even symmetry in a viscous accretion disc. This mode occurs in eccentric discs and is known to be potentially overstable. We determine the vertical structure of the disc and its modes, treating radiative energy transport in the diffusion approximation. In the limit of very long radial wavelength, an analytical criterion for viscous overstability is obtained, which involves the effective shear and bulk viscosity, the adiabatic exponent, and the opacity law of the disc. This differs from the prediction of a two-dimensional model. On shorter wavelengths (a few times the disc thickness), the criterion for overstability is more difficult to satisfy because of the different vertical structure of the mode. In a low-viscosity disc a third regime of intermediate wavelengths appears, in which the overstability is suppressed as the horizontal velocity perturbations develop significant vertical shear. We suggest that this effect determines the damping rate of eccentricity in protoplanetary discs, for which the long-wavelength analysis is inapplicable and overstability is unlikely to occur on any scale. In thinner accretion discs and in decretion discs around Be stars overstability may occur only on the longest wavelengths, leading to the preferential excitation of global eccentric modes.     

The potential for Earth-mass planet formation around brown dwarfs

Recent observations point to the presence of structured dust grains in the discs surrounding young brown dwarfs, thus implying that the first stages of planet formation take place also in the substellar regime. Here, we investigate the potential for planet formation around brown dwarfs and very low-mass stars according to the sequential core accretion model of planet formation. We find that, for a brown dwarf mass 0.05 M_⊙, our models predict a maximum planetary mass of  ∼5   M_⊕  , orbiting with semimajor axis ∼ 1 au. However, we note that the predictions for the mass–semimajor axis distribution are strongly dependent upon the models chosen for the disc surface density profiles and the assumed distribution of disc masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, Earth-mass planets might be relatively frequent, while if typical disc masses are only a fraction of Jupiter mass, we predict that planet formation would be extremely rare in the substellar regime. As the observational constraints on disc profiles, mass dependencies and their distributions are poor in the brown dwarf regime, we advise caution in validating theoretical models only on stars similar to the Sun and emphasize the need for observational data on planetary systems around a wide range of stellar masses. We also find that, unlike the situation around solar-like stars, Type II migration is totally absent from the planet formation process around brown dwarfs, suggesting that any future observations of planets around brown dwarfs would provide a direct measure of the role of other types of migration.     

Accretion disc outbursts: a new version of an old model

We have developed 1D time-dependent numerical models of accretion discs, using an adaptive grid technique and an implicit numerical scheme, in which the disc size is allowed to vary with time. The code fully resolves the cooling and heating fronts propagating in the disc. We show that models in which the radius of the outer edge of the disc is fixed produce incorrect results, from which probably incorrect conclusions about the viscosity law have been inferred. In particular we show that outside-in outbursts are possible when a standard bimodal behaviour of the Shakura–Sunyaev viscosity parameter α is used. We also discuss to what extent insufficient grid resolution has limited the predictive power of previous models. We find that the global properties (magnitudes, etc.) of transient discs can be addressed by codes using a high, but reasonable, number of fixed grid points. However, the study of the detailed physical properties of the transition fronts generally requires resolutions which are out of reach of fixed grid codes. It appears that most time-dependent models of accretion discs published in the literature have been limited by resolution effects, improper outer boundary conditions, or both.     

Analysing the atolls: X-ray spectral transitions of accreting neutron stars

We systematically analyse all the available X-ray spectra of disc accreting neutron stars (atolls and millisecond pulsars) from the RXTE data base. We show that while all these have similar spectral evolution as a function of mass accretion rate, there are also subtle differences. There are two different types of hard/soft transition, those where the spectrum softens at all energies, leading to a diagonal track on a colour–colour diagram, and those where only the higher energy spectrum softens, giving a vertical track. The luminosity at which the transition occurs is correlated with this spectral behaviour, with the vertical transition at   L / L _Edd∼ 0.02  while the diagonal one is at ∼0.1. Superimposed on this is the well-known hysteresis effect, but we show that classic, large-scale hysteresis occurs only in the outbursting sources, indicating that its origin is in the dramatic rate of change of mass accretion rate during the disc instability. We show that the long-term mass accretion rate correlates with the transition behaviour, and speculate that this is due to the magnetic field being able to emerge from the neutron star surface for low average mass accretion rates. While this is not strong enough to collimate the flow except in the millisecond pulsars, its presence may affect the inner accretion flow by changing the properties of the jet.     

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.     

The outburst properties of AM CVn stars

We briefly summarise the observational properties of ultra‐compact binaries called AM CVn stars. We analyse their outbursts originating from the thermal‐viscous instability in helium accretion disc. We present our preliminary results in applying the model of Dwarf Novae outbursts to helium discs. We can calculate models of outbursts of reasonable amplitude of 2 mag with a constant α parameter throughout the calculation. Setting the mass transfer rate close to its upper critical value produces model lightcurves that resemble short superoutbursts (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The photoevaporation of discs around young stars in massive clusters

We present models in which the photoevaporation of discs around young stars by an external ultraviolet source (as computed by Adams et al.) is coupled with the internal viscous evolution of the discs. These models are applied to the case of the Orion Nebula Cluster (ONC), where the presence of a strong ultraviolet field from the central OB stars, together with a detailed census of circumstellar discs and photoevaporative flows, is well established. In particular we investigate the constraints that are placed on the initial disc properties in the ONC by the twin requirement that most stars possess a disc on a scale of a few astronomical unit (au), but that only a minority (<20 per cent) are resolved by Hubble Space Telescope ( HST ) at a scale of 50 au. We find that these requirements place very weak constraints on the initial radius distribution of circumstellar discs: the resulting size distribution readily forgets the initial radius distribution, owing to the strong positive dependence of the photoevaporation rate on disc radius. Instead, the scarcity of large discs reflects the relative scarcity of initially massive discs (with mass  >0.1 M_⊙  ). The ubiquity of discs on a small scale, on the other hand, mainly constrains the time-span over which the discs have been exposed to the ultraviolet field (<2 Myr). We argue that the discs that are resolved by HST represent a population of discs in which self-gravity was important at the time that the dominant central OB star switched on, but that, according to our models, self-gravity is unlikely to be important in these discs at the present time. We discuss the implications of our results for the so-called proplyd lifetime problem.     

An accretion disc surrounding a component of epsilon aurigae

The low-mass model of Epsilon Aurigae consisting of an F supergiant star a Main-Sequence B star which is surrounded by a cool accretion disc has been investigated. The disc requires a high accretion rate if we assume the disc to be completely heated up by friction alone. On the other hand, discs with low accretion rates are possible when the disc are heated up from the radiation of both stars. A swollen crescent portion produced by the radiation of the F star may form a bulky body sufficient to obscure the radiation of the F star.Astronomical Institute, Tohoku University, Sendai, Japan.     

Modelling the behaviour of accretion flows in X-ray binaries

We review how the recent increase in X-ray and radio data from black hole and neutron star binaries can be merged together with theoretical advances to give a coherent picture of the physics of the accretion flow in strong gravity. Both long term X-ray light curves, X-ray spectra, the rapid X-ray variability and the radio jet behaviour are consistent with a model where a standard outer accretion disc is truncated at low luminosities, being replaced by a hot, inner flow which also acts as the launching site of the jet. Decreasing the disc truncation radius leads to softer spectra, as well as higher frequencies (including quasi periodic oscillations, QPOs) in the power spectra, and a faster jet. The collapse of the hot flow when the disc reaches the last stable orbit triggers the dramatic decrease in radio flux, as well as giving a qualitative (and often quantitative) explanation for the major hard–soft spectral transition seen in black holes. The neutron stars are also consistent with the same models, but with an additional component due to their surface, giving implicit evidence for the event horizon in black holes. We review claims of observational data which conflict with this picture, but show that these can also be consistent with the truncated disc model. We also review suggested alternative models for the accretion flow which do not involve a truncated disc. The most successful of these converge on a similar geometry, where there is a transition at some radius larger than the last stable orbit between a standard disc and an inner, jet dominated region, with the X-ray source associated with a mildly relativistic outflow, beamed away from the disc. However, the observed uniformity of properties between black holes at different inclinations suggests that even weak beaming of the X-ray emission may be constrained by the data. After collapse of the hot inner flow, the spectrum in black hole systems can be dominated by the disc emission. Its behaviour is consistent with the existence of a last stable orbit, and such data can be used to estimate the black hole spin. By contrast, these systems can also show very different spectra at these high luminosities, in which the disc spectrum (and probably structure) is strongly distorted by Comptonization. The structure of the accretion flow becomes increasingly uncertain as the luminosity approaches (and exceeds) the Eddington luminosity, though there is growing evidence that winds may play an important role. We stress that these high Eddington fraction flows are key to understanding many disparate and currently very active fields such as ULX, Narrow Line Seyfert 1’s, and the growth of the first black holes in the Early Universe.     

Numerical simulations of protostellar encounters — II. Coplanar disc–disc encounters

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. Such interactions must have a significant impact upon the evolution of the protostars and their discs.   We have carried out a series of simulations of coplanar encounters between two stars, each possessing a massive circumstellar disc, using an SPH code that models gravitational, hydrodynamic and viscous forces. We find that during a coplanar encounter, disc material is swept up into a shock layer between the two interacting stars, and the layer then fragments to produce new protostellar condensations. The truncated remains of the discs may subsequently fragment; and the outer regions of the discs may be thrown off to form circumbinary disc-like structures around the stars. Thus coplanar disc–disc encounters lead efficiently to the formation of multiple star systems and small- N clusters, including substellar objects.