X-ray irradiation in XTE J1817−330 and the inner radius of the truncated disc in the hard state

The key aspect of the very successful truncated disc model for the low/hard X-ray spectral state in black hole binaries is that the geometrically thin disc recedes from the last stable orbit at the transition to this state. This has recently been challenged by direct observations of the low/hard state disc from CCD data. We reanalyse the Swift and RXTE campaign covering the 2006 outburst of XTE J1817−330, and show that these data actually strongly support the truncated disc model as the transition spectra unambiguously show that the disc begins to recede as the source leaves the disc-dominated soft state. The disc radius inferred for the proper low/hard state is less clear-cut, but we show that the effect of irradiation from the energetically dominant hot plasma leads to an underestimate of the disc radius by a factor of 2–3 in this state. This may also produce the soft excess reported in some hard-state spectra. The inferred radius becomes still larger when the potential difference in stress at the inner boundary, increased colour temperature correction from incomplete thermalization of the irradiation, and loss of observable disc photons from Comptonization in the hot plasma is taken into account. We conclude that the inner disc radius in XTE J1817−330 in the low/hard spectral state is at least six to eight times that seen in the disc-dominated high/soft state, and that recession of the inner disc is the trigger for the soft-hard-state transition, as predicted by the truncated disc models.     

Alignment time-scale of the microquasar GRO J1655−40

The microquasar GRO J1655−40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady-state disc warped by Lense–Thirring precession and find the time-scale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically 50 times higher in the latter case but we find that, in both the cases, the lifetime of the mass-transfer state is at most a few times the alignment time-scale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter aligned after its supernova kick but that it is most likely to be close to alignment rather than counter alignment now.     

Reprocessing of X-rays in the outer accretion disc of the black hole binary XTE J1817−330

We build a simple model of the optical/ultraviolet (UV) emission from irradiation of the outer disc by the inner disc and coronal emission in black hole binaries. We apply this to the broad-band Swift data from the outburst of the black hole binary XTE J1817−330 to confirm previous results that the optical/UV emission in the soft state is consistent with a reprocessing a constant fraction of the bolometric X-ray luminosity. However, this is very surprising as the disc temperature drops by more than a factor of 3 in the soft state, which should produce a marked change in the reprocessing efficiency. The easiest way to match the observed constant reprocessed fraction is for the disc skin to be highly ionized (as suggested 30 yr ago by van Paradijs), so that the bulk of the disc flux is reflected and only the hardest X-rays heat the disc. The constant reprocessed fraction also favours direct illumination of the disc over a scattering origin as the optical depth/solid angle of any scattering material (wind/corona) over the disc should decrease as the source luminosity declines. By contrast, the reprocessed fraction increases very significantly (by a factor of ∼6) as the source enters the hard state. This dramatic change is not evident from X-ray/UV flux correlations as it is masked by bandpass effects. However, it does not necessarily signal a change in emission, for example, the emergence of the jet dominating the optical/UV flux as the reflection albedo must change with the dramatic change in spectral shape.     

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.     

Investigating a fluctuating-accretion model for the spectral-timing properties of accreting black hole systems

The fluctuating-accretion model of Lyubarskii and its extension by Kotov, Churazov & Gilfanov seek to explain the spectral-timing properties of the X-ray variability of accreting black holes in terms of inward-propagating mass accretion fluctuations produced at a broad range of radii. The fluctuations modulate the X-ray emitting region as they move inwards and can produce temporal-frequency-dependent lags between energy bands, and energy-dependent power spectral densities (PSDs) as a result of the different emissivity profiles, which may be expected at different X-ray energies. Here, we use a simple numerical implementation to investigate in detail the X-ray spectral-timing properties of the model and their relation to several physically interesting parameters, namely the emissivity profile in different energy bands, the geometrical thickness and viscosity parameter of the accretion flow, the strength of damping on the fluctuations and the temporal coherence (measured by the 'quality factor', Q ) of the fluctuations introduced at each radius. We find that a geometrically thick flow with large viscosity parameter is favoured, and we confirm that the predicted lags are quite robust to changes in the emissivity profile and physical parameters of the accretion flow, which may help to explain the similarity of the lag spectra in the low/hard and high/soft states of Cyg X-1. We also demonstrate the model regime where the light curves in different energy bands are highly spectrally coherent. We compare model predictions directly to X-ray data from the narrow line Seyfert 1 galaxy NGC 4051 and the black hole X-ray binary (BHXRB) Cyg X-1 in its high/soft state, and we show that this general scheme can reproduce simultaneously the time lags and energy-dependence of the PSD.     

Negative Reynolds stress generation by accretion disc convection

The phenomenon of negative viscosity-alpha in convectively unstable Keplerian accretion discs is discussed. The convection is considered as a random flow with an axisymmetric mesoscale pattern. Its correlation tensor is computed with a time-averaging procedure using Kley's 2D hydrocode. There is a distinct anisotropy between the turbulence intensities in the radial and azimuthal directions, i.e. the radial velocity rms dominates the azimuthal one. As a consequence, an extra term in the expression for the turbulent transport of angular momentum appears which does not vanish for rigid rotation ('Λ-effect'). It is negative ('inwards transport') and even seems to dominate the positive contribution of the eddy viscosity representing outwards transport of angular momentum. For a turbulence model close to that of the mixing-length theory, the rotational influence on the anisotropy of the turbulence intensities,     , and the covariance  〈 u ' _R u ' _φ 〉  – representing the angular momentum transport – is computed and compared with the accretion disc simulations. Indeed, the negative angular momentum transport can be explained with the observed dominance of the radial turbulence intensity. If, on the other hand, in turbulence fields the azimuthal intensity would dominate or the turbulence is even isotropic, then we always find a positive transport of the angular momentum.     

Numerical simulations of type III planetary migration – I. Disc model and convergence tests

We investigate the fast (type III) migration regime of high-mass protoplanets orbiting in protoplanetary discs. This type of migration is dominated by corotational torques. We study the details of flow structure in the planet's vicinity, the dependence of migration rate on the adopted disc model and the numerical convergence of models (independence of certain numerical parameters such as gravitational softening).
We use two-dimensional hydrodynamical simulations with adaptive mesh refinement, based on the flash code with improved time-stepping scheme. We perform global disc simulations with sufficient resolution close to the planet, which is allowed to freely move throughout the grid. We employ a new type of equation of state in which the gas temperature depends on both the distance to the star and planet, and a simplified correction for self-gravity of the circumplanetary gas.
We find that the migration rate in the type III migration regime depends strongly on the gas dynamics inside the Hill sphere (Roche lobe of the planet) which, in turn, is sensitive to the aspect ratio of the circumplanetary disc. Furthermore, corrections due to the gas self-gravity are necessary to reduce numerical artefacts that act against rapid planet migration. Reliable numerical studies of type III migration thus require consideration of both the thermal and the self-gravity corrections, as well as a sufficient spatial resolution and the calculation of disc–planet attraction both inside and outside the Hill sphere. With this proviso, we find type III migration to be a robust mode of migration, astrophysically promising because of a speed much faster than in the previously studied modes of migration.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Non-linear dynamics of the corotation torque

The excitation of spiral waves by an external perturbation in a disc deposits angular momentum in the vicinity of the corotation resonance (the radius where the speed of a rotating pattern matches the local rotation rate). We use matched asymptotic expansions to derive a reduced model that captures non-linear dynamics of the resulting torque and fluid motions. The model is similar to that derived for forced Rossby wave critical layers in geophysical fluid dynamics. Using the model we explore the saturation of the corotation torque, which occurs when the background potential (specific) vorticity is redistributed by the disturbance. We also consider the effects of dissipation. If there is a radial transport of potential vorticity, the corotation torque does not saturate. The main application is to the creation, growth and migration of protoplanets within discs like the primordial solar nebula. The disturbance also nucleates vortices in the vicinity of corotation, which may spark further epochs of planet formation.     

Magnetized accretion flows: effects of gas pressure

We study how axisymmetric magnetohydrodynamic (MHD) accretion flows depend on γ adiabatic index in the polytropic equation of state. This work is an extension of Mościbrodzka & Proga, where we investigated the γ dependence of two-dimensional Bondi-like accretion flows in the hydrodynamical (HD) limit. Our main goal is to study if simulations for various γ can give us insights into the problem of various modes of accretion observed in several types of accretion systems, such as black hole binaries (BHBs), active galactic nuclei (AGN) and gamma-ray bursts. We find that for  γ≳ 4/3  , the fast-rotating flow forms a thick torus that is supported by rotation and gas pressure. As shown before for  γ= 5/3  , such a torus produces a strong, persistent bipolar outflow that can significantly reduce the polar funnel accretion of a slowly rotating flow. For low γ, close to 1, the torus is thin and is supported by rotation. The thin torus produces an unsteady outflow which is too weak to propagate throughout the polar funnel inflow. Compared to their HD counterparts, the MHD simulations show that the magnetized torus can produce an outflow and does not exhibit regular oscillations. Generally, our simulations demonstrate how the torus thickness affects the outflow production. They also support the notion that the geometrical thickness of the torus correlates with the power of the torus outflow. Our results, applied to observations, suggest that the torus ability to radiatively cool and become thin can correspond to a suppression of a jet as observed in the BHBs during a transition from a hard/low to soft/high spectral state and a transition from a quiescent to hard/low state in AGN.     

On the evolution of giant protoplanets forming in circumbinary discs

We present the results of hydrodynamic simulations of Jovian mass protoplanets that form in circumbinary discs. The simulations follow the orbital evolution of the binary plus protoplanet system acting under their mutual gravitational forces, and forces exerted by the viscous circumbinary disc. The evolution involves the clearing of the inner circumbinary disc initially, so that the binary plus protoplanet system orbits within a low density cavity. Continued interaction between disc and protoplanet causes inward migration of the planet towards the inner binary. Subsequent evolution can take three distinct paths: (i) the protoplanet enters the 4 : 1 mean motion resonance with the binary, but is gravitationally scattered through a close encounter with the secondary star; (ii) the protoplanet enters the 4 : 1 mean motion resonance, the resonance breaks, and the planet remains in a stable orbit just outside the resonance; (iii) when the binary has initial eccentricity   e _bin≥ 0.2  , the disc becomes eccentric, leading to a stalling of the planet migration, and the formation of a stable circumbinary planet.
These results have implications for a number of issues in the study of extrasolar planets. The ejection of protoplanets in close binary systems provides a source of 'free-floating planets', which have been discovered recently. The formation of a large, tidally truncated cavity may provide an observational signature of circumbinary planets during formation. The existence of protoplanets orbiting stably just outside a mean motion resonance (4 : 1) in the simulations indicate that such sites may harbour planets in binary star systems, and these could potentially be observed. Finally, the formation of stable circumbinary planets in eccentric binary systems indicates that circumbinary planets may not be uncommon.     

Eccentric instability of a two-dimensional fluid disc with a central massive object

In this paper we investigate, by linear modal analysis, the one-armed dynamical instability of a two-dimensional fluid disc that has a massive object at its centre. The model of the disc is chosen to avoid the artificial instabilities that originate from the unrealistic disc configurations that have been adopted in previous studies. We find a one-armed instability for which the central massive object is displaced from the centre, which is generally called the 'eccentric instability'. However, to excite the eccentric instability, the mass of the central object should be appreciably smaller than that of the disc, and this mass ratio is far smaller than what was originally proposed. The instability shown in this paper is likely to be excited in a stellar system with a central massive object, e.g. a galactic nucleus harbouring a massive black hole, and further studies are desirable via techniques such as numerical simulations.     

The aperiodic broad-band X-ray variability of Cygnus X-3

We study the soft X-ray variability of Cygnus X-3. By combining data from the All-Sky Monitor and Proportional Counter Array instruments on the RXTE satellite with EXOSAT /Medium Energy (ME) detector observations, we are able to analyse the power density spectrum (PDS) of the source from 10⁻⁹ to 0.1 Hz, thus covering time-scales from seconds to years. As the data on the longer time-scales are unevenly sampled, we combine traditional power spectral techniques with simulations to analyse the variability in this range. The PDS at higher frequencies  (≳10⁻³ Hz)  are for the first time compared for all states of this source. We find that it is for all states well described by a power law, with index  ∼−2  in the soft states and a tendency for a less steep power law in the hard state. At longer time-scales, we study the effect of the state transitions on the PDS, and find that the variability below  ∼10⁻⁷ Hz  is dominated by the transitions. Furthermore, we find no correlation between the length of a high/soft-state episode and the time since the previous high/soft state. On intermediate time-scales, we find evidence for a break in the PDS at time-scales of the order of the orbital period. This may be interpreted as evidence for the existence of a tidal resonance in the accretion disc around the compact object, and constraining the mass ratio to   M ₂/ M ₁≲ 0.3  .     

Launching of conical winds and axial jets from the disc–magnetosphere boundary: axisymmetric and 3D simulations

We investigate the launching of outflows from the disc–magnetosphere boundary of slowly and rapidly rotating magnetized stars using axisymmetric and exploratory 3D magnetohydrodynamic simulations. We find long-lasting outflows in the following cases. (1) In the case of slowly rotating stars , a new type of outflow, a conical wind , is found and studied in simulations. The conical winds appear in cases where the magnetic flux of the star is bunched up by the disc into an X-type configuration. The winds have the shape of a thin conical shell with a half-opening angle  θ∼ 30°–40°  . About 10–30 per cent of the disc matter flows from the inner disc into the conical winds. The conical winds may be responsible for episodic as well as long-lasting outflows in different types of stars. There is also a low-density, higher velocity component (a jet) in the region inside the conical wind. (2) In the case of rapidly rotating stars (the 'propeller regime'), a two-component outflow is observed. One component is similar to the conical winds. A significant fraction of the disc matter may be ejected into the winds. The second component is a high-velocity, low-density magnetically dominated axial jet where matter flows along the opened polar field lines of the star. The jet has a mass flux of about 10 per cent of that of the conical wind, but its energy flux (dominantly magnetic) can be larger than the energy flux of the conical wind. The jet's angular momentum flux (also dominantly magnetic) causes the star to spin down rapidly. Propeller-driven outflows may be responsible for the jets in protostars and for their rapid spin-down. The jet is collimated by the magnetic force while the conical winds are only weakly collimated in the simulation region. Exploratory 3D simulations show that conical winds are axisymmetric about the rotational axis (of the star and the disc), even when the dipole field of the star is significantly misaligned.     

Unipolar outflows and global meridional circulations in rotating accretion flows

Using two-dimensional simulations of non-radiative viscous rotating black hole accretion flows, we show that the flows with α ∼0.1–0.3 self-organize to form stationary unipolar or bipolar outflows accompanied by global meridional circulations. The required energy comes, with efficiency ∼0.001–0.01, from the matter directly accreted on to the black hole. Observational implications are discussed.     

Evolution of black hole intermediate-mass X-ray binaries: the influence of a circumbinary disc

We report the results of a systematic timing analysis of RXTE observations of GRS 1915+105 when the source was in its variability class θ, characterized by alternating soft and hard states on a time-scale of a few hundred seconds. The aim was to examine the high-frequency part of the power spectrum in order to confirm the hectohertz quasi-periodic oscillations (QPO) previously reported from observations from mixed variability behaviours. During the hard intervals (corresponding to state C in the classification of Belloni et al.), we find a significant QPO at a frequency of ∼170 Hz, although much broader (Q∼2) than previously reported. No other significant peak is observed at frequencies >30 Hz. A time-resolved spectral analysis of selected observations shows that the hard intervals from class θ show a stronger and steeper  (Γ= 2.8–3.0)  power-law component than hard intervals from other classes. We discuss these results in the framework of hectohertz QPOs reported from GRS 1915+105 and other black hole binaries.     

Accretion disc viscosity: how big is alpha?

We consider observational and theoretical estimates of the accretion disc viscosity parameter α. We find that in thin, fully ionized discs, the best observational evidence suggests a typical range α∼ 0.1–0.4, whereas the relevant numerical simulations tend to derive estimates for α which are an order of magnitude smaller. We discuss possible reasons for this apparent discrepancy.     

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.     

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.