Outbursts of irradiated accretion discs

I solve analytically the viscous evolution of an irradiated accretion disc, as seen during outbursts of soft X-ray transients. The solutions predict steep power-law X-ray decays L _X ∼ (1 + t/t_visc)⁻⁴, changing to L _X ∼ (1 − t/t'_visc)⁴ at late times, where t _visc, t '_visc are viscous time-scales. These forms closely resemble the approximate exponential and linear decays inferred by King and Ritter in these two regimes. The decays are much steeper than for unirradiated discs because the viscosity is a function of the central accretion rate rather than of local conditions in the disc.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

PDS 456: an extreme accretion rate quasar?

We present quasi-simultaneous ASCA and RXTE observations of the most luminous known active galactic nucleus in the local ( z <0.3) Universe, the recently discovered quasar PDS 456. Multiwavelength observations have been conducted that show that PDS 456 has a bolometric luminosity of ∼10⁴⁷ erg s⁻¹ peaking in the ultraviolet part of the spectrum. In the X-ray band the 2–10 keV (rest-frame) luminosity is 10⁴⁵ erg s⁻¹. The broad-band X-ray spectrum obtained with ASCA and RXTE contains considerable complexity. The most striking feature observed is a very deep, ionized iron K edge, observed at 8.7 keV in the quasar rest-frame. We find that these features are consistent with reprocessing from highly ionized matter, probably the inner accretion disc. PDS 456 appeared to show a strong (factor of ∼2.1) outburst in just ∼17 ks, although non-intrinsic sources cannot be completely ruled out. If confirmed, this would be an unusual event for such a high-luminosity source, with a light-crossing-time corresponding to ∼2 R _S . The implication would be that flaring occurs within the very central regions, or else that PDS 456 is a 'super-Eddington' or relativistically beamed system. Overall we conclude on the basis of the extreme blue/UV luminosity, the rapid X-ray variability and from the imprint of highly ionized material on the X-ray spectrum, that PDS 456 is a quasar with an unusually high accretion rate.     

Frequency-resolved spectroscopy of Cyg X-1: fast variability of the reflected emission in the soft state

Using RXTE /PCA data, we study the fast variability of the reflected emission in the soft spectral state of Cyg X-1 by means of Fourier frequency-resolved spectroscopy. We find that the rms amplitude of variations of the reflected emission has the same frequency dependence as the primary radiation down to time-scales of ≲30–50 ms. This might indicate that the reflected flux reproduces, with nearly flat response, variations of the primary emission. Such behaviour differs notably from that of the hard spectral state, in which variations of the reflected flux are significantly suppressed in comparison with the primary emission, on time-scales shorter than ∼0.5–1 s.
If related to the finite light-crossing time of the reflector, these results suggest that the characteristic size of the reflector, presumably an optically thick accretion disc, in the hard spectral state is larger by a factor of ≳5–10 than in the soft spectral state. Modelling the transfer function of the disc, we estimate the inner radius of the accretion disc to be R _in∼100 R _g in the hard state and R _in≲10 R _g in the soft state for a 10-M_⊙ black hole.     

Correlated multiwavelength emission from the X-ray-bright Seyfert galaxy III Zw 2

The X-ray-bright Seyfert 1 galaxy III Zw 2 was observed with XMM–Newton in 2000 July. Its X-ray spectrum can be described by a power law of photon index Γ= 1.7 and an extremely broad (FWHM∼ 140 000 km  s⁻¹  ) Fe Kα line at 6.44 keV. The iron line has an equivalent width of ∼800 eV. To study the long-term X-ray behaviour of the source we have analysed 25 yr of data, from 1975 to 2000. There is no evidence of significant intrinsic absorption within the source or of a soft X-ray excess in the XMM or archival data. We do not detect rapid X-ray variability (a few  × 10³ s  ) during any of the individual observations; however, on longer time-scales (a few years) the X-ray light curve shows 10-fold flux variations. We infer a black hole mass of  ∼10⁹ M_⊙  (from Hβ FWHM) for III Zw 2 which is much higher than some previous estimates.
A comparison of X-ray variability with light curves at other wavelengths over a 25-yr period reveals correlated flux variations from radio to X-ray wavelengths. We interpret the variable radio to optical emission as synchrotron radiation, self-absorbed in the radio/millimetre region, and the X-rays mainly as a result of Compton up-scattering of low-energy photons by the population of high-energy electrons that give rise to the synchrotron radiation.     

The X-ray and extreme-ultraviolet flux evolution of SS Cygni throughout outburst

We present the most complete multiwavelength coverage of any dwarf nova outburst: simultaneous optical, Extreme Ultraviolet Explorer and Rossi X-ray Timing Explorer observations of SS Cygni throughout a narrow asymmetric outburst. Our data show that the high-energy outburst begins in the X-ray waveband 0.9–1.4 d after the beginning of the optical rise and 0.6 d before the extreme-ultraviolet rise. The X-ray flux drops suddenly, immediately before the extreme-ultraviolet flux rise, supporting the view that both components arise in the boundary layer between the accretion disc and white dwarf surface. The early rise of the X-ray flux shows that the propagation time of the outburst heating wave may have been previously overestimated.
The transitions between X-ray and extreme-ultraviolet dominated emission are accompanied by intense variability in the X-ray flux, with time-scales of minutes. As detailed by Mauche & Robinson, dwarf nova oscillations are detected throughout the extreme-ultraviolet outburst, but we find they are absent from the X-ray light curve.
X-ray and extreme-ultraviolet luminosities imply accretion rates of  3 × 10¹⁵ g s⁻¹  in quiescence,  1 × 10¹⁶ g s⁻¹  when the boundary layer becomes optically thick, and  ∼10¹⁸ g s⁻¹  at the peak of the outburst. The quiescent accretion rate is two and a half orders of magnitude higher than predicted by the standard disc instability model, and we suggest this may be because the inner accretion disc in SS Cyg is in a permanent outburst state.     

The redshift evolution of clustering in the Hubble Deep Field

We present a correlation function analysis for the catalogue of photometric redshifts obtained from the Hubble Deep Field image by Fernandez-Soto, Lanzetta & Yahil. By dividing the catalogue into redshift bins of width Δ z =0.4 we measured the angular correlation function w ( θ ) as a function of redshift up to z ∼4.8. From these measurements we derive the trend of the correlation length r ₀. We find that r ₀( z ) is roughly constant with look-back time up to z ≃2, and then increases to higher values at z ≳2.4. We estimate the values of r ₀, assuming ξ ( r , z )=[ r r ₀( z )]^{− γ} , γ =1.8 and various geometries. For Ω₀=1 we find r ₀( z =3)≃7.00±4.87  h ⁻¹ Mpc, in good agreement with the values obtained from analysis of the Lyman break galaxies.     

Supernova 1994aj: a probe for pre-supernova evolution and mass-loss from the progenitor

Extensive photometric and spectroscopic observations of SN 1994aj until 540 d after maximum light have been obtained. The photometry around maximum suggests that the SN belongs to the Type II Linear class, with a peak absolute magnitude of M _V∼−17.8 (assuming H ₀=75 km s⁻¹ Mpc⁻¹). The spectra of SN 1994aj were unusual, with the presence of a narrow line with a P Cygni profile on top of the broad Balmer line emission. This narrow feature is attributed to the presence of a dense superwind surrounding the SN. At 100–120 d after maximum light the SN ejecta start to interact with this circumstellar material. The SN luminosity decline rates slowed down [γ _R =0.46 mag (100 d)⁻¹], becoming less steep than the average late luminosity decline of normal SN II [∼1 mag (100 d)⁻¹]. This dense ( ˙M / u _W∼10¹⁵ g cm⁻¹) wind was confined to a short distance from the progenitor ( R _out=∼5×10¹⁶ cm), and results from a very strong mass-loss episode ( ˙M =10⁻³ M⊙ yr⁻¹), which terminated shortly before explosion (∼5–10 yr).     

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.     

Spectral variability of ultraluminous X-ray sources

We study spectral variability of 11 ultraluminous X-ray sources (ULX) using archived XMM–Newton and Chandra observations. We use three models to describe the observed spectra: a power law, a multicolour disc (MCD) and a combination of these two models. We find that seven ULXs show a correlation between the luminosity L _X and the photon index Γ. Furthermore, four out of these seven ULXs also show spectral pivoting in the observed energy band. We also find that two ULXs show an   L _X–Γ  anticorrelation. The spectra of four ULXs in the sample can be adequately fitted with a MCD model. We compare these sources to known black hole binaries (BHB) and find that they follow similar paths in their luminosity–temperature diagrams. Finally, we show that the 'soft excess' reported for many of these ULXs at ∼0.2 keV seems to roughly follow a trend   L _soft∝ T ^−3.5  when modelled with a power law plus a 'cool' MCD model. This is contrary to the   L ∝ T ⁴  relation that is expected from theory and what is seen for many accreting BHBs. The observed trend could instead arise from disc emission beamed by an outflowing wind around a  ∼10 M_⊙  black hole.     

The evolution of galactic discs

We use recent observations of high-redshift galaxies to study the evolution of galactic discs over the redshift range 0 <  z ≲1. The data are inconsistent with models in which discs were already assembled at z  = 1 and have evolved only in luminosity since that time. Assuming that disc properties change with redshift as powers of 1 +   z and analysing the observations assuming an Einstein–de Sitter universe, we find that for given rotation speed, disc scalelength decreases with z as ∼ (1 +  z )⁻¹, total B -band mass-to-light ratio decreases with z as ∼ (1 +  z )⁻¹, and disc luminosity (again in B ) depends only weakly on z . These scalings are consistent with current data on the evolution of disc galaxy abundance as a function of size and luminosity. Both the scalings and the abundance evolution are close to the predictions of hierarchical models for galaxy formation. If different cosmogonies are compared, the observed evolution in disc size and disc abundance favours a flat low-Ω₀ universe over an Einstein–de Sitter universe.     

High-frequency X-ray variability as a mass estimator of stellar and supermassive black holes

There is increasing evidence that supermassive black holes in active galactic nuclei (AGN) are scaled-up versions of Galactic black holes. We show that the amplitude of high-frequency X-ray variability in the hard spectral state is inversely proportional to the black hole mass over eight orders of magnitude. We have analysed all available hard-state data from RXTE of seven Galactic black holes. Their power density spectra change dramatically from observation to observation, except for the high-frequency (≳10 Hz) tail, which seems to have a universal shape, roughly represented by a power law of index −2. The amplitude of the tail,   C _M   (extrapolated to 1 Hz), remains approximately constant for a given source, regardless of the luminosity, unlike the break or quasi-periodic oscillation frequencies, which are usually strongly correlated with luminosity. Comparison with a moderate-luminosity sample of AGN shows that the amplitude of the tail is a simple function of black hole mass,   C _M = C / M   , where   C ≈ 1.25 M_⊙ Hz⁻¹  . This makes   C _M   a robust estimator of the black hole mass which is easy to apply to low- to moderate-luminosity supermassive black holes. The high-frequency tail with its universal shape is an invariant feature of a black hole and, possibly, an imprint of the last stable orbit.     

Large-Scale Coherent Dipole Anisotropy?

We have reanalysed and compared the dipoles of the 1.2-Jy and 0.6-Jy (QDOT) IRAS galaxy samples. We find strong indications from both samples for (a) significant contributions to the gravitational field that shapes the Local Group motion from depths up to ∼ 170 h ⁻¹ Mpc and (b) a large-scale coherence of the dipole anisotropy, indications provided mainly by the fact that the differential dipoles of large equal-volume shells are aligned with the CMB dipole and exhibit significant dipole signals. The two IRAS dipoles are indistinguishable within 50 h ⁻¹ Mpc and beyond ∼ 130 h ⁻¹ Mpc while the QDOT dipole, having a lower flux limit, continues growing with respect to the 1.2-Jy sample up to ∼ 100 h ⁻¹ Mpc in agreement with Rowan-Robinson et al.     

Application of a clumpy core model to NGC 2071

The multi-transitional observations of CS molecules towards the NGC 2071 core have been re-analysed by using a tri-dimensional Monte Carlo radiative transfer code. Better agreement with the observations is made by an introduction of clumpiness to this model than by smoothly varying density to the 1D microturbulent one. The best-fitting model shows that, when a unique density is assumed for clumps, the volume filling factor of the clumps varies as r ⁻² with an average of ∼5 per cent over the entire core, and that the H₂ number density and the CS abundance of the clump relative to H₂ are ∼ 2 × 10⁶ cm⁻³ and ∼ 6 × 10⁻¹⁰, respectively. The radial density gradient ∝ r ⁻² obtained from our clumpy core model is steeper than that (∝ r ^−1.3) obtained from the microturbulent model. Since all clumps are subject to random bulk motions in this 3D clumpy macroturbulent model, synthesized line profiles do not show self-absorption dips even for opaque transitions and the resulting linewidth is in good accordance with the observations.     

Circumplanetary disc properties obtained from radiation hydrodynamical simulations of gas accretion by protoplanets

We investigate the properties of circumplanetary discs formed in three-dimensional, self-gravitating radiation hydrodynamical models of gas accretion by protoplanets. We determine disc sizes, scaleheights, and density and temperature profiles for different protoplanet masses, in solar nebulae of differing grain opacities.
We find that the analytical prediction of circumplanetary disc radii in an evacuated gap  ( R _Hill/3)  from Quillen & Trilling yields a good estimate for discs formed by high-mass protoplanets. The radial density profiles of the circumplanetary discs may be described by power laws between   r ⁻²  and   r ^−3/2  . We find no evidence for the ring-like density enhancements that have been found in some previous models of circumplanetary discs. Temperature profiles follow a  ∼ r ^−7/10  power law regardless of protoplanet mass or nebula grain opacity. The discs invariably have large scaleheights  ( H / r > 0.2)  , making them thick in comparison with their encompassing circumstellar discs, and they show no flaring.     

A good long look at the black hole candidates LMC X-1 and LMC X-3

LMC X-1 and LMC X-3 are the only known persistent stellar-mass black-hole candidates that have almost always shown spectra that are dominated by a soft, thermal component. We present here results from 170-ks-long Rossi X-ray Timing Explorer ( RXTE ) observations of these objects, taken in 1996 December, where their spectra can be described by a disc blackbody plus an additional soft     high-energy power law (detected up to energies of 50 keV in LMC X-3). These observations, as well as archival Advanced Satellite for Cosmology and Astrophysics ( ASCA ) observations, constrain any narrow Fe line present in the spectra to have an equivalent width ≲90 eV. Stronger, broad lines (≈150 eV EW,     are permitted. We also study the variability of LMC X-1. Its X-ray power spectral density (PSD) is approximately proportional to     between 10⁻³ and 0.3 Hz with a root-mean-square (rms) variability of ≈7 per cent. At energies >5 keV, the PSD shows evidence of a break at     possibly indicating an outer disc radius of ≲1000  GM c ² in this likely wind-fed system. Furthermore, the coherence function     a measure of the degree of linear correlation between variability in the >5 keV band and variability in the lower energy bands, is extremely low (≲50 per cent). We discuss the implications of these observations for the mechanisms that might be producing the soft and hard X-rays in these systems.     

Stellar-mass black hole binaries as ultraluminous X-ray sources

Ultraluminous X-ray sources (ULXs) with   L _x > 10³⁹ erg s⁻¹  have been discovered in great numbers in external galaxies with ROSAT , Chandra and XMM-Newton . The central question regarding this important class of sources is whether they represent an extension in the luminosity function of binary X-ray sources containing neutron stars and stellar-mass black holes (BHs), or a new class of objects, e.g. systems containing intermediate-mass BHs  (100–1000 M_⊙)  . We have carried out a theoretical study to test whether a large fraction of the ULXs, especially those in galaxies with recent star formation activity, can be explained with binary systems containing stellar-mass BHs. To this end, we have applied a unique set of binary evolution models for BH X-ray binaries, coupled to a binary population synthesis code, to model the ULXs observed in external galaxies. We find that for donor stars with initial masses  ≳10 M_⊙  the mass transfer driven by the normal nuclear evolution of the donor star is sufficient to potentially power most ULXs. This is the case during core hydrogen burning and, to an even more pronounced degree, while the donor star ascends the giant branch, although the latter phases last only ∼5 per cent of the main-sequence phase. We show that with only a modest violation of the Eddington limit, e.g. a factor of ∼10, both the numbers and properties of the majority of the ULXs can be reproduced. One of our conclusions is that if stellar-mass BH binaries account for a significant fraction of ULXs in star-forming galaxies, then the rate of formation of such systems is  ∼3 × 10⁻⁷ yr⁻¹  normalized to a core-collapse supernova rate of 0.01 yr⁻¹.     

Multiwavelength study of the nuclei of a volume-limited sample of galaxies – II. Optical, infrared and radio observations

We present optical and infrared broad-band images, radio maps, and optical spectroscopy for the nuclear region of a sample of nearby galaxies. The galaxies have been drawn from a complete volume-limited sample for which we have already presented X-ray imaging. We modelled the stellar component of the spectroscopic observations to determine the star formation history of our targets. Diagnostic diagrams were used to classify the emission-line spectra and determine the ionizing mechanism driving the nuclear regions. All those sources classified as active galactic nuclei present small Eddington ratios  (∼10⁻³–10⁻⁶)  , implying a very slow growth rate of their black holes. We finally investigate the relative numbers of active and normal nuclei as a function of host galaxy luminosity and find that the fraction of active galaxies slowly rises as a function of host absolute magnitude in the   M _B ∼−12  to −22 range.     

Sakurai's object: the ionized nebula at radio wavelengths

Sakurai's object (V4334 Sgr) is a planetary nebula nucleus which is undergoing its final helium shell flash. This is the first of these rare and important events to be observable with non-optical instruments. We report the first radio detection, using a short (2-h) observation with the Very Large Array (VLA) at 4.86 GHz. The radio emission structure is coincident with the 34-arcsec diameter planetary nebula seen in optical emission lines. We find a statistical distance ∼ 3.8 ± 0.6 kpc, with a range of 1.9 <  D  < 5.3 kpc, depending on the planetary nebula (PN) mass. While we have no direct evidence for a new (post-flash) stellar wind, we estimate an upper limit to the mass-loss rate due to any such wind of 1.7 × 10⁻⁷ M⊙ yr⁻¹. The number of emitting knots in the radio-visible nebula indicates an electron density of ∼ 2 × 10⁸ m⁻³ in those knots, and a total emitting ionized mass of ∼ 0.15 M⊙, at an assumed distance of 3.8 kpc. The radio flux density indicates an Hβ flux of ∼ 6 × 10⁻¹⁶ W m⁻², suggesting an extinction E ( B  −  V ) ∼ 1.15, comparable with reddening estimates in the direction of V4334 Sgr.     

ASCA observations of two steep soft X-ray quasars

Steep soft X-ray (0.1–2 keV) quasars share several unusual properties: narrow Balmer lines, strong Fe  II emission, large and fast X-ray variability, and a rather steep 2–10 keV spectrum. These intriguing objects have been suggested to be the analogues of Galactic black hole candidates in the high, soft state. We present here results from ASCA observations for two of these quasars: NAB 0205 + 024 and PG 1244 + 026.   Both objects show similar variations (factor of ∼ 2 in 10 ks), despite a factor of ∼ 10 difference in the 0.5–10 keV luminosity (7.3 × 10⁴³ erg s⁻¹ for PG 1244 + 026 and 6.4 × 10⁴⁴ erg s⁻¹ for NAB 0205 + 024, assuming isotropic emission, H ₀ = 50.0 and q ₀ = 0.0).   The X-ray continuum of the two quasars flattens by 0.5–1 going from the 0.1–2 keV band towards higher energies, strengthening recent results on another half-dozen steep soft X-ray active galactic nuclei.   PG 1244 + 026 shows a significant feature in the '1-keV' region, which can be described either as a broad emission line centred at 0.95 keV (quasar frame) or as edge or line absorption at 1.17 (1.22) keV. The line emission could be a result of reflection from a highly ionized accretion disc, in line with the view that steep soft X-ray quasars are emitting close to the Eddington luminosity. Photoelectric edge absorption or resonant line absorption could be produced by gas outflowing at a large velocity (0.3–0.6 c ).