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.     

Disc precession in the M31 dipping X-ray binary Bo 158?

We present results from three XMM–Newton observations of the M31 low mass X-ray binary (LMXB) XMMU J004314.4+410726.3 (Bo 158), spaced over 3 d in 2004 July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-h period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anticorrelated with source intensity. In light of the 2004 XMM–Newton observations of Bo 158, we suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio ≲0.3 (period ≲4 h), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of  ∼11 P _orb  , and prograde disc precession on a period of  29 ± 1 P _orb  . This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on an  81 ± 3 h  cycle.     

Discovery of drifting high-frequency quasi-periodic oscillations in global simulations of magnetic boundary layers

We report on the numerical discovery of quasi-periodic oscillations (QPOs) associated with accretion through a non-axisymmetric magnetic boundary layer in the unstable regime, when two ordered equatorial streams form and rotate synchronously at approximately the angular velocity of the inner disc. The streams hit the star's surface producing hotspots. Rotation of the spots leads to high-frequency QPOs. We performed a number of simulation runs for different magnetospheric sizes from small to tiny, and observed a definite correlation between the inner disc radius and the QPO frequency: the frequency is higher when the magnetosphere is smaller. In the stable regime, a small magnetosphere forms and accretion through the usual funnel streams is observed, and the frequency of the star is expected to dominate the light curve. We performed exploratory investigations of the case in which the magnetosphere becomes negligibly small and the disc interacts with the star through an equatorial belt. We also performed investigation of somewhat larger magnetospheres where one or two ordered tongues may dominate over other chaotic tongues. In application to millisecond pulsars, we obtain QPO frequencies in the range of 350–990 Hz for one spot. The frequency associated with rotation of one spot may dominate if spots are not identical or antipodal. If the spots are similar and antipodal, then the frequencies are twice as high. We show that variation of the accretion rate leads to drift of the QPO peak.     

Atlas of tilted accretion discs and source to negative superhumps

Using smoothed particle hydrodynamics, we numerically simulate steady-state accretion discs for cataclysmic variable dwarf novae systems that have a secondary-to-primary mass ratio  0.35 ≤ q ≤ 0.55  . After these accretion discs have come to quasi-equilibrium, we rotate each disc out of the orbital plane by  δ= (1, 2, 3, 4, 5 or 20)°  to induce negative superhumps. For accretion discs tilted  5°  , we generate light curves and associated Fourier transforms for an atlas on negative superhumps and retrograde precession. Our simulation results suggest that accretion discs need to be tilted more than 3° for negative superhumps to be statistically significant. We also show that if the disc is tilted enough such that the gas stream strikes a disc face, then a dense cooling ring is generated near the radius of impact.
In addition to the atlas, we study these artificially tilted accretion discs to find the source to negative superhumps. Our results suggest that the source is additional light from innermost disc annuli, and this additional light waxes and wanes with the amount of gas stream overflow received as the secondary orbits. The nodes, where the gas stream transitions from flowing over to under the disc rim (and vice versa), precess in the retrograde direction.     

Timing of the accreting millisecond pulsar XTE J1814−338

We present a precise timing analysis of the accreting millisecond pulsar XTE J1814−338 during its 2003 outburst, observed by RXTE . A full orbital solution is given for the first time; Doppler effects induced by the motion of the source in the binary system were corrected, leading to a refined estimate of the orbital period,   P _orb= 15 388.7229(2)  s, and of the projected semimajor axis,   a sin  i / c = 0.390633(9)  light-second. We could then investigate the spin behaviour of the accreting compact object during the outburst. We report here a refined value of the spin frequency  (ν= 314.356 108 79(1) Hz)  and the first estimate of the spin frequency derivative of this source while accreting     . This spin-down behaviour arises when both the fundamental frequency and the second harmonic are taken into consideration. We discuss this in the context of the interaction between the disc and the quickly rotating magnetosphere, at accretion rates sufficiently low to allow a threading of the accretion disc in regions where the Keplerian velocity is slower than the magnetosphere velocity. We also present indications of a jitter of the pulse phases around the mean trend, which we argue results from movements of the accreting hotspots in response to variations of the accretion rate.     

Magnetic flares and the optical variability of the X-ray transient XTE J1118+480

The simultaneous presence of a strong quasi-periodic oscillation, of period ∼10 s, in the optical and X-ray light curves of the X-ray transient XTE J1118+480 suggests that a significant fraction of the optical flux originates from the inner part of the accretion flow, where most of the X-rays are produced. We present a model of magnetic flares in an accretion disc corona where thermal cyclo-synchrotron emission contributes significantly to the optical emission, while the X-rays are produced by inverse Compton scattering of the soft photons produced by dissipation in the underlying disc and by the synchrotron process itself. Given the observational constraints, we estimate the values for the coronal temperature, optical depth and magnetic field intensity, as well as the accretion rate for the source. Within our model we predict a correlation between optical and hard X-ray variability and an anticorrelation between optical and soft X-rays. We also expect optical variability on flaring time-scales (∼tens of ms), with a power-density spectrum similar to that observed in the X-ray band. Finally, we use both the available optical/extreme-ultraviolet/X-ray spectral energy distribution and the low-frequency time variability to discuss limits on the inner radius of the optically thick disc.     

Quasi-periodic oscillations in low-mass X-ray binaries and constraints on the equation of state of neutron star matter

Recently discovered quasi-periodic oscillations in the X-ray brightness of low-mass X-ray binaries are used to derive constraints on the mass of the neutron star component and the equation of state of neutron star matter. The observations are compared with models of rapidly rotating neutron stars which are calculated by means of an exact numerical method in full relativity. For the equations of state we select a broad collection of models representing different assumptions about the many-body structure and the complexity of the composition of superdense matter. The mass constraints differ from their values in the approximate treatment by ∼10 per cent. Under the assumption that the maximum frequency of the quasi-periodic oscillations originates from the innermost stable orbit, the mass of the neutron star is in the range M ∼1.92–2.25 M_⊙. The quasi-periodic oscillation in the Atoll-source 4U 1820−30 in particular is only consistent with equations of state that are rather stiff at high densities, which is explainable, so far, only with pure nucleonic/leptonic composition. This interpretation contradicts the hypothesis that the protoneutron star formed in SN 1987A collapsed to a black hole, since this would demand a maximum neutron star mass below 1.6 M_⊙. The recently suggested identification of quasi-periodic oscillations with frequencies of about 10 Hz with the Lense–Thirring precession of the accretion disc is found to be inconsistent with the models studied in this work, unless it is assumed that the first overtone of the precession is observed.     

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.     

An irradiated accretion disc in the narrow-line Seyfert 1 RE J1034+396?

We model the optical to X-ray continuum spectrum of the narrow-line Seyfert 1 galaxy RE J1034+396. We show that the flat optical spectrum is consistent with emission from an irradiated accretion disc. The X-ray emission can be modelled with a disc blackbody and a Comptonized component. The temperature at the inner edge of the disc     Using this constraint, we show that the flat optical spectrum is consistent with emission from the irradiatively heated outer part of the accretion disc. We constrain the outer radius of the optically thick disc     and the inner radius of the irradiation-dominated region     . Our optical and X-ray spectral fits indicate a mass     , and do not rule out a low (i.e. face-on) inclination angle for the system.     

Three-dimensional smoothed particle hydrodynamics simulations of radiation-driven warped accretion discs

We present three-dimensional smoothed particle hydrodynamics calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc, resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and, in general, the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios, using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios   q ∼ 0.1  a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (cf. X 1916−053). For less extreme mass ratios, the entire disc tilts out of the orbital plane (cf. Her X–1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction, we obtain both positive and negative superhumps simultaneously in the dissipation light curve (cf. V603 Aql).     

Weighing the black holes in ultraluminous X-ray sources through timing

We describe a new method to estimate the mass of black holes in Ultraluminous X-ray Sources (ULXs). The method is based on the recently discovered 'variability plane', populated by Galactic stellar-mass black-hole candidates (BHCs) and supermassive active galactic nuclei (AGNs), in the parameter space defined by the black-hole mass, accretion rate and characteristic frequency. We apply this method to the two ULXs from which low-frequency quasi-periodic oscillations have been discovered, M82 X-1 and NGC 5408 X-1. For both sources we obtain a black-hole mass in the range  100–1300 M_⊙  , thus providing evidence for these two sources to host an intermediate-mass black hole.     

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.     

Possible quasi-periodic oscillations from unstable accretion: 3D magnetohydrodynamic simulations

We investigate the photometric variability of magnetized stars, particularly neutron stars, accreting through a magnetic Rayleigh–Taylor-type instability at the disc–magnetosphere interface, and compare it with the variability during stable accretion, with the goal of looking for possible quasi-periodic oscillations (QPOs). The light curves during stable accretion show periodicity at the star's frequency and sometimes twice that, due to the presence of two funnel streams that produce antipodal hotspots near the magnetic poles. On the other hand, light curves during unstable accretion through tongues penetrating the magnetosphere are more chaotic due to the stochastic behaviour of the tongues, and produce noisier power spectra. However, the power spectra do show some signs of quasi-periodic variability. Most importantly, the rotation frequency of the tongues and the resulting hotspots are close to the inner-disc orbital frequency, except in the most strongly unstable cases. There is therefore a high probability of observing QPOs at that frequency in longer simulations. In addition, the light curves in the unstable regime show periodicity at the star's rotation frequency in many of the cases investigated here, again except in the most strongly unstable cases which lack funnel flows and the resulting antipodal hotspots. The noisier power spectra result in the fractional rms amplitudes of the Fourier peaks being smaller.
We also study in detail the effect of the misalignment angle between the rotation and magnetic axes of the star on the variability, and find that at misalignment angles  ≳25°  the star's period always appears in the light curves.     

The impact of a supernova explosion in a very massive binary

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

How galaxies lose their angular momentum

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disc galaxies that are too compact with respect to the observations. High-resolution N -body/SPH simulations in a cosmological context are presented including cold gas and dark matter (DM). A halo with quiet merging activity since redshift   z ∼ 3.8  and with a high-spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disc. We show that the gas and the DM have similar specific angular momenta until a merger event occurs at   z ∼ 2  with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and dynamical friction processes falls quickly into the centre. In contrast, gas infall through small subclumps or accretion does not lead to catastrophic angular momentum loss. In fact, a new extended disc begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disc galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered starburst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected, late-type galaxies could form with a dominant extended disc component, resulting from late infall, a small bulge-to-disc ratio and a low baryon fraction, in agreement with observations.     

The minimum gap-opening planet mass in an irradiated circumstellar accretion disc

We consider the minimum mass planet, as a function of radius, that is capable of opening a gap in an α-accretion disc. We estimate that a half-Jupiter mass planet can open a gap in a disc with accretion rate     for viscosity parameter  α= 0.01  , and solar mass and luminosity. The minimum mass is approximately proportional to     . This estimate can be used to rule out the presence of massive planets in gapless accretion discs. We identify two radii at which an inwardly migrating planet may become able to open a gap and so slow its migration; the radius at which the heating from viscous dissipation is similar to that from stellar radiation in a flared disc, and the radius at which the disc becomes optically thin in a self-shadowed disc. In the inner portions of the disc, we find that the minimum planet mass required to open a gap is only weakly dependent on radius. If a migrating planet is unable to open a gap by the time it reaches either of the transition radii, then it is likely to be lost on to the star. If a gap-opening planet cuts off disc accretion allowing the formation of a central hole or clearing in the disc then we would estimate that the clearing radius would approximately be proportional to the stellar mass.     

Doppler signatures of Hα flares in AE Aquarii

We present high-speed spectrophotometric observations of the Hα emission in the exotic cataclysmic variable AE Aqr. Over 10 000 spectra with a sampling rate of ∼ 6 s were obtained. The well-known 16.5- and 33-s oscillations were not detected in either the continuum or the line, nor were any quasi-periodic oscillations detected. The Hα line shows remarkably complex behaviour, with changes in flux of up to a factor of 3 in a few hours. The line does not act as a single entity — the instantaneous line profile is highly asymmetric and consists of a superposition of emission components that are localized in velocity. We find that during flares the Hα emission is delayed and decays more slowly than the continuum. The Hα variations are not linearly correlated with the continuum variations. Examination of the trailed spectrograms and Doppler tomograms reveals no signatures of an accretion disc, but rather supports the idea that mass is being ejected by a propeller mechanism.     

RXTE–PCA observations of XMMU J054134.7−682550

We analysed Rossi X-ray Timing Explorer Proportional Counter Array observations of a recent outburst of the X-ray pulsar XMMU J054134.7−682550. We calculated the pulse frequency history of the source. We found no sign of a binary companion. The source spins up when the X-ray flux is higher, with a correlation between the spin-up rate and X-ray flux, which may be interpreted as a sign of an accretion disc. On the other hand, the source was found to have an almost constant spin frequency when the X-ray flux is lower without any clear sign of a spin-down episode. The decrease in pulsed fraction with decreasing X-ray flux was interpreted as a sign of accretion geometry change, but we did not find any evidence of a transition from accretor to propeller regimes. The source was found to have variable pulse profiles. Two peaks in pulse profiles were usually observed. We studied the X-ray spectral evolution of the source throughout the observation. Pulse-phase-resolved analysis does not provide any further evidence for a cyclotron line, but may suggest a slight variation of intensity and width of the 6.4 keV iron line with phase.     

Feedback and the formation of dwarf galaxy stellar haloes

Stellar population studies show that low-mass galaxies in all environments exhibit stellar haloes that are older and more spherically distributed than the main body of the galaxy. In some cases, there is a significant intermediate age component that extends beyond the young disc. We examine a suite of Smoothed Particle Hydrodynamic simulations and find that elevated early star formation activity combined with supernova feedback can produce an extended stellar distribution that resembles these haloes for model galaxies ranging from   v ₂₀₀= 15  to 35 km s⁻¹, without the need for accretion of subhaloes.     

Shock oscillation model for quasi-periodic oscillations in stellar mass and supermassive black holes

We numerically examine centrifugally supported shock waves in 2D rotating accretion flows around a stellar mass  (10 M_⊙)  and a supermassive  (10⁶ M_⊙)  black holes over a wide range of input accretion rates of     . The resultant 2D shocks are unstable with time and the luminosities show quasi-periodic oscillations (QPOs) with modulations of a factor of 2–3 and with periods of a tenth of a second to several hours, depending on the black hole masses. The shock oscillation model may explain the intermediate frequency QPOs with 1–10 Hz observed in the stellar mass black hole candidates and also suggest the existence of QPOs with the period of hours in active galactic nuclei. When the accretion rate     is low, the luminosity increases in proportion to the accretion rate. However, when     greatly exceeds the Eddington critical rate     , the luminosity is insensitive to the accretion rate and is kept constantly around  ∼3 L _E  . On the other hand, the mass-outflow rate     increases in proportion to     and it amounts to about a few per cent of the input mass-flow rate.