On the complex disc–corona interactions in the soft spectral states of soft X-ray transients

Accreting black holes show a complex and diverse behaviour in their soft spectral states. Although these spectra are dominated by a soft, thermal component which almost certainly arises from an accretion disc, there is also a hard X-ray tail indicating that some fraction of the accretion power is instead dissipated in hot, optically thin coronal material. During such states, best observed in the early outburst of soft X-ray transients, the ratio of power dissipated in the hot corona to that in the disc can vary from ∼ 0 (pure disc accretion) to ∼ 1 (equal power in each). Here we present results of spectral analyses of a number of sources, demonstrating the presence of complex features in their energy spectra. Our main findings are: (1) the soft components are not properly described by a thermal emission from accretion discs: they are appreciably broader than can be described by disc blackbody models even including relativistic effects, and (2) the spectral features near     commonly seen in such spectra can be well described by reprocessing of hard X-rays by optically thick, highly ionized, relativistically moving plasma.     

MS 1603.6+2600: an accretion disc corona source?

We have observed the eclipsing low-mass X-ray binary MS 1603.6+2600 with Chandra for 7 ks. The X-ray spectrum is well fit with a single absorbed power law with an index of ∼2. We find a clear sinusoidal modulation in the X-ray light curve with a period of  1.7 ± 0.2 h  , consistent with the period of 1.85 h found before. However, no (partial) eclipses were found. We argue that if the X-ray flare observed in earlier X-ray observations was a type I X-ray burst, then the source can only be an accretion disc corona source at a distance of ∼11–24 kpc (implying a height above the Galactic disc of ∼8–17 kpc). It has also been proposed in the literature that MS 1603.76+2600 is a dipper at ∼75 kpc. We argue that, in this dipper scenario, the observed optical properties of MS 1603.6+2600 are difficult to reconcile with the optical properties one would expect on the basis of comparisons with other high-inclination, low-mass X-ray binaries, unless the X-ray flare was not a type I X-ray burst. In that case, the source can be a nearby soft X-ray transient accreting at a quiescent rate, as was proposed by Hakala et al., or a high-inclination source at ∼15–20 kpc.     

Simultaneous X-ray/optical observations of GX 9+9 (4U 1728−16)

We report on the results of the first simultaneous X-ray ( RXTE ) and optical [South African Astronomical Observatory (SAAO)] observations of the luminous low-mass X-ray binary (LMXB) GX 9+9 in 1999 August. The high-speed optical photometry revealed an orbital period of 4.1958 h and confirmed previous observations, but with greater precision. No X-ray modulation was found at the orbital period. On shorter time-scales, a possible 1.4-h variability was found in the optical light curves which might be related to the MHz quasi-periodic oscillations seen in other LMXBs. We do not find any significant X-ray/optical correlation in the light curves. In X-rays, the colour–colour and hardness-intensity diagrams indicate that the source shows characteristics of an atoll source in the upper banana state, with a correlation between intensity and spectral hardness. Time-resolved X-ray spectroscopy suggests that two-component spectral models give a reasonable fit to the X-ray emission. Such models consist of a blackbody component which can be interpreted as the emission from an optically thick accretion disc or an optically thick boundary layer, and a hard Comptonized component for an extended corona.     

Keck II spectroscopy of mHz quasi-periodic oscillations in Hercules X-1

We present Keck II spectroscopy of optical mHz quasi-periodic oscillations (QPOs) in the light curve of the X-ray pulsar binary Hercules X-1. In the power spectrum it appears as 'peaked noise', with a coherency ∼2, a central frequency of 35 mHz and a peak-to-peak amplitude of 5 per cent. However, the dynamic power spectrum shows it to be an intermittent QPO, with a lifetime of ∼100 s, as expected if the lifetime of the orbiting material is equal to the thermal time-scale of the inner disc. We have decomposed the spectral time series into constant and variable components and used blackbody fits to the resulting spectra to characterize the spectrum of the QPO variability and constrain possible production sites. We find that the spectrum of the QPO is best fitted by a small hot region, possibly the inner regions of the accretion disc, where the ballistic accretion stream impacts on to the disc. The lack of any excess power around the QPO frequency in the X-ray power spectrum, created using simultaneous light curves from RXTE , implies that the QPO is not simply reprocessed X-ray variability.     

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

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

A comparative analysis of standard accretion discs spectra: an application to ultraluminous X-ray sources

We compare standard models of accretion discs around black holes (BHs) that include the appropriate zero-torque inner boundary condition and relativistic effects on the emission and propagation of radiation. The comparison is performed adopting the multicolour disc blackbody model (MCD) as reference and looking for the parameter space in which it is in statistical agreement with 'more physical' accretion disc models. We find simple 'recipes' that can be used for adjusting the estimates of the physical inner radius of the disc, the BH mass and the accretion rate inferred using the parameters of the MCD fits. We applied these results to four ultraluminous X-ray sources for which MCD spectral fits of their X-ray soft spectral components have been published and find that, in three cases (NGC 1313 X-1, X-2 and M 81 X-9), the BH masses inferred for a standard disc around a Schwarzschild BH are in the interval  ∼100–200 M_⊙  . Only if the BH is maximally rotating are the masses comparable to the much larger values previously derived in the literature.     

X-ray spectral evolution of low-mass X-ray binary GX 349+2

We present the results of a systematic investigation of spectral evolution in the Z source GX 349+2, using data obtained during 1998 with the Proportional Counter Array (PCA) on-board the RXTE satellite. The source traced a extended normal branch (NB) and flaring branch (FB) in the colour–colour diagram (CD) and the hardness-intensity diagram (HID) during these observations. The spectra at different positions of the Z-track were best fitted by a model consisting of a disc blackbody and a Comptonized spectrum. A broad (Gaussian) iron line at ∼6.7 keV is also required to improve the fit. The spectral parameters showed a systematic and significant variation with the position along the Z-track. The evolution in spectral parameters is discussed in view of the increasing mass accretion rate scenario, proposed to explain the motion of Z sources in the CD and the HID.     

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.     

XMM–Newton observations of UW CrB: detection of X-ray bursts and evidence for accretion disc evolution

UW CrB (MS 1603+2600) is a peculiar short-period X-ray binary that exhibits extraordinary optical behaviour. The shape of the optical light curve of the system changes drastically from night to night, without any changes in overall brightness. Here we report X-ray observations of UW CrB obtained with XMM–Newton . We find evidence for several X-ray bursts, confirming a neutron star primary. This considerably strengthens the case that UW CrB is an accretion disc corona system located at a distance of at least 5–7 kpc (3–5 kpc above the Galactic plane). The X-ray and Optical Monitor (ultraviolet–optical) light curves show remarkable shape variation from one observing run to another, which we suggest are due to large-scale variations in the accretion disc shape resulting from a warp that periodically obscures the optical and soft X-ray emission. This is also supported by the changes in phase-resolved X-ray spectra.     

Energy-dependent variability from accretion flows

We develop a formalism to calculate energy-dependent fractional variability (rms) in accretion flows. We consider rms spectra resulting from radial dependences of the level of local variability (as expected from the propagation of disturbances in accretion flows) assuming the constant shape of the spectrum emitted at a given radius. We consider the cases when the variability of the flow is either coherent or incoherent between different radial zones. As an example of local emission, we consider blackbody, Wien and thermal Comptonization spectra. In addition to numerical results, we present a number of analytical formulae for the resulting rms. We also find an analytical formula for the disc Wien spectrum, which we find to be a very good approximation to the disc blackbody. We compare our results to the rms spectrum observed in an ultrasoft state of GRS 1915+105.     

X-ray and ultraviolet observations of the dwarf nova VW Hyi in quiescence

We present an analysis of X-ray and ultraviolet (UV) data of the dwarf nova VW Hyi that were obtained with XMM–Newton during the quiescent state. The X-ray spectrum indicates the presence of an optically thin plasma in the boundary layer that cools as it settles on to the white dwarf. The plasma has a continuous temperature distribution that is well described by a power law or a cooling flow model with a maximum temperature of 6–8 keV. We estimate from the X-ray spectrum a boundary layer luminosity of  8 × 10³⁰ erg s^-1  , which is only 20 per cent of the disc luminosity. The rate of accretion on to the white dwarf is  5 × 10⁻¹² M_⊙ yr⁻¹  , about half of the rate in the disc. From the high-resolution X-ray spectra, we estimate that the X-ray emitting part of the boundary layer is rotating with a velocity of 540 km s⁻¹, which is close to the rotation velocity of the white dwarf but is significantly smaller than the Keplerian velocity. We detect a 60-s quasi-periodic oscillation of the X-ray flux, which is likely to be due to the rotation of the boundary layer. The X-ray and the UV flux show strong variability on a time-scale of ∼1500 s. We find that the variability in the two bands is correlated and that the X-ray fluctuations are delayed by ∼100 s. The correlation indicates that the variable UV flux is emitted near the transition region between the disc and the boundary layer and that accretion rate fluctuations in this region are propagated to the X-ray emitting part of the boundary layer within ∼100 s. An orbital modulation of the X-ray flux suggests that the inner accretion disc is tilted with respect to the orbital plane. The elemental abundances in the boundary layer are close to their solar values.     

Discovery of a broad iron line in the black hole candidate Swift J1753.5−0127, and the disc emission in the low/hard state revisited

We analysed simultaneous archival XMM–Newton and Rossi X-ray Timing Explorer observations of the X-ray binary and black hole candidate Swift J  1753.5−0127  . In a previous analysis of the same data, a soft thermal component was found in the X-ray spectrum, and the presence of an accretion disc extending close to the innermost stable circular orbit was proposed. This is in contrast with the standard picture in which the accretion disc is truncated at large radii in the low/hard state. We tested a number of spectral models and found that several of them fit the observed spectra without the need of a soft disc-like component. This result implies that the classical paradigm of a truncated accretion disc in the low/hard state cannot be ruled out by these data. We further discovered a broad iron emission line between 6 and 7 keV in these data. From fits to the line profile we found an inner disc radius that ranges between ∼6 and 16 gravitational radii, which can be in fact much larger, up to ∼250 gravitational radii, depending on the model used to fit the continuum and the line. We discuss the implications of these results in the context of a fully or partially truncated accretion disc.     

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.     

On the observability of spiral structures in cataclysmic variable accretion discs

We use the grid of hydrodynamic accretion disc calculations of Stehle to construct orbital phase‐dependent emission‐line profiles of thin discs carrying spiral density waves. The observational signatures of spiral waves are explored to establish the feasibility of detecting spiral waves in cataclysmic variable discs using prominent emission lines in the visible range of the spectrum. For high Mach number accretion discs ( M v _φ c _s≃ 15 – 30), we find that the spiral shock arms are so tightly wound that they leave few obvious fingerprints in the emission lines. Only a minor variation of the double peak separation in the line profile at a level of ∼8 per cent is produced. For accretion discs in outburst ( M ≃ 5 – 20) however, the lines are dominated by the emission from an m =2 spiral pattern in the disc. We show that reliable Doppler tomograms of spiral shock patterns can be reconstructed provided that a signal‐to‐noise ratio of at least 15, a wavelength resolution of ∼80 km s⁻¹ and a time resolution of ∼50 spectra per binary orbit are achieved. We confirm that the observed spiral pattern in the disc of IP Pegasi can be reproduced by tidal density waves in the accretion disc and demands the presence of a large, hot disc, at least in the early outburst stages.     

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.     

Ultraluminous X-ray source 1E 0953.8+6918 (M81 X-9): an intermediate-mass black hole candidate and its environs

We present a ROSAT and ASCA study of the Einstein source X-9 and its relation to a shock-heated shell-like optical nebula in a tidal arm of the M81 group of interacting galaxies. Our ASCA observation of the source shows a flat and featureless X-ray spectrum well described by a multicolour disc blackbody model. The source most likely represents an optically thick accretion disc around an intermediate-mass black hole  ( M ∼10² M_⊙)  in its high/soft state, similar to other variable ultraluminous X-ray sources observed in nearby disc galaxies. Using constraints derived from both the innermost stable orbit around a black hole and the Eddington luminosity, we find that the black hole is fast-rotating and that its mass is between ∼80 M_⊙–1.5×10² M_⊙. The inferred bolometric luminosity of the accretion disc is ∼(1.1×10⁴⁰ erg s⁻¹)/(cos  i ). Furthermore, we find that the optical nebula is very energetic and may contain large amounts of hot gas, accounting for a soft X-ray component as indicated by archival ROSAT PSPC data. The nebula is apparently associated with X-9; the latter may be powering the former and/or they could be formed in the same event (e.g. a hypernova). Such a connection, if confirmed, could have strong implications for understanding both the birth of intermediate-mass black holes and the formation of energetic interstellar structures.     

Multicolour observations of V404 Cyg with ULTRACAM

We present high-time-resolution multicolour observations of the quiescent soft X-ray transient V404 Cyg obtained with ULTRACAM. Superimposed on the ellipsoidal modulation of the secondary star are large flares on time-scales of a few hours, as well as several distinct rapid flares on time-scales of tens of minutes. The rapid flares, most of which show further variability and unresolved peaks, cover shorter time-scales than those reported in previous observations. The power density spectrum of the 5-s time-resolution data shows a quasi-periodic oscillation (QPO) feature at 0.78 mHz (=21.5 min). Assuming this periodicity represents the Keplerian period at the transition between the thin and advective disc regions, we determine the transition radius. We discuss the possible origins for the QPO feature in the context of the advection-dominated accretion flow model.
We determine the colour of the large flares and find that the i '-band flux per unit frequency interval is larger than that in the g ' band. The colour is consistent with optically thin gas with a temperature of ∼8000 K arising from a region with an equivalent blackbody radius of at least  2 R_⊙  , which covers 3 per cent of the surface of the accretion disc. Our timing and spectral analysis results support the idea that the rapid flares (i.e. the QPO feature) most likely arise from regions near the transition radius.     

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.     

Evidence of a change in the long-term spin-down rate of the X-ray pulsar 4U 1907+09

We analysed RXTE archival observations of 4U 1907+09 between 1996 February 17 and 2002 March 6. The pulse timing analysis showed that the source stayed at almost constant period around 1998 August and then started to spin-down at a rate of  (−1.887 ∓ 0.042) × 10⁻¹⁴ Hz s⁻¹  which is ∼0.60 times lower than the long-term (∼15 yr) spin-down rate. Our pulse-frequency measurements for the first time resolved significant spin-down rate variations since the discovery of the source. We also presented orbital phase resolved X-ray spectra during two stable spin-down episodes during 1996 November–1997 December and 2001 March–2002 March. The source has been known to have two orbitally locked flares. We found that X-ray flux and spectral parameters except hydrogen column density agreed with each other during the flares. We interpreted the similar values of X-ray fluxes as an indication of the fact that the source accretes not only via transient retrograde accretion disc but also via the stellar wind of the companion, so that the variation of the accretion rate from the disc does not cause significant variation in the observed X-ray flux. Lack of significant change in spectral parameters except hydrogen column density was interpreted as a sign of the fact that the change in the spin-down rate of the source was not accompanied by a significant variation in the accretion geometry.     

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.