Weakly non-linear stability of a hydrodynamic accretion disc

When the cold accretion disc coupling between neutral gas and a magnetic field is so weak that the magnetorotational instability is less effective or even stops working, it is of prime interest to investigate the pure hydrodynamic origin of turbulence and transport phenomena. As the Reynolds number increases, the relative importance of the non-linear term in the hydrodynamic equation increases. In an accretion disc where the molecular viscosity is too small, the Reynolds number is large enough for the non-linear term to have new effects. We investigate the scenario of the `weakly non-linear' evolution of the amplitude of the linear mode when the flow is bounded by two parallel walls. The unperturbed flow is similar to the plane Couette flow, but with the Coriolis force included in the hydrodynamic equation. Although there is no exponentially growing eigenmode, because of the self-interaction, the least stable eigenmode will grow in an intermediate phase. Later, this will lead to higher-order non-linearity and plausible turbulence. Although the non-linear term in the hydrodynamic equation is energy-conserving, within the weakly non-linear analysis it is possible to define a lower bound of the energy (alpha A(c)(2), where A(c) is the threshold amplitude) needed for the flow to transform to the turbulent phase. Such an unstable phase is possible only if the Reynolds number >= 10(3-4). The numerical difficulties in obtaining such a large Reynolds number might be the reason for the negative result of numerical simulations on a pure hydrodynamic Keplerian accretion disc.

Multidimensional modelling of X-ray spectra for AGN accretion disc outflows - III. Application to a hydrodynamical simulation

We perform multidimensional radiative transfer simulations to compute spectra for a hydrodynamical simulation of a line-driven accretion disc wind from an active galactic nucleus. The synthetic spectra confirm expectations from parametrized models that a disc wind can imprint a wide variety of spectroscopic signatures including narrow absorption lines, broad emission lines and a Compton hump. The formation of these features is complex with contributions originating from many of the different structures present in the hydrodynamical simulation. In particular, spectral features are shaped both by gas in a successfully launched outflow and in complex flows where material is lifted out of the disc plane but ultimately falls back. We also confirm that the strong Fe Ka line can develop a weak, red-skewed line wing as a result of Compton scattering in the outflow. In addition, we demonstrate that X-ray radiation scattered and reprocessed in the flow has a pivotal part in both the spectrum formation and determining the ionization conditions in the wind. We find that scattered radiation is rather effective in ionizing gas which is shielded from direct irradiation from the central source. This effect likely makes the successful launching of a massive disc wind somewhat more challenging and should be considered in future wind simulations. © 2010 The Authors. Journal compilation © 2010 RAS.

Multidimensional modelling of X-ray spectra for AGN accretion disc outflows

We use a multidimensional Monte Carlo code to compute X-ray spectra for a variety of active galactic nucleus (AGN) disc-wind outflow geometries. We focus on the formation of blueshifted absorption features in the Fe K band and show that line features similar to those which have been reported in observations are often produced for lines of sight through disc-wind geometries. We also discuss the formation of other spectral features in highly ionized outflows. In particular, we show that, for sufficiently high wind densities, moderately strong Fe K emission lines can form and that electron scattering in the flow may cause these lines to develop extended red wings. We illustrate the potential relevance of such models to the interpretation of real X-ray data by comparison with observations of a well-known AGN, Mrk 766. Journal compilation © 2008 RAS.

Multidimensional modelling of X-ray spectra for AGN accretion disc outflows - II

Highly ionized fast accretion disc winds have been suggested as an explanation for a variety of observed absorption and emission features in the X-ray spectra of active galactic nuclei. Simple estimates have suggested that these flows may be massive enough to carry away a significant fraction of the accretion energy and could be involved in creating the link between supermassive black holes and their host galaxies. However, testing these hypotheses, and quantifying the outflow signatures, requires high-quality theoretical spectra for comparison with observations. Here, we describe extensions of our Monte Carlo radiative transfer code that allow us to generate realistic theoretical spectra for a much wider variety of disc wind models than that was possible in our previous work. In particular, we have expanded the range of atomic physics simulated by the code so that L- and M-shell ions can now be included. We have also substantially improved our treatment of both ionization and radiative heating such that we are now able to compute spectra for outflows containing far more diverse plasma conditions. We present example calculations that illustrate the variety of spectral features predicted by parametrized outflow models and demonstrate their applicability to the interpretation of data by comparison with observations of the bright quasar PG1211+143. We find that the major features in the observed 2-10 keV spectrum of this object can be well reproduced by our spectra, confirming that it likely hosts a massive outflow. © 2010 The Authors. Journal compilation © 2010 RAS.

The impact of accretion disc winds on the optical spectra of cataclysmic variables

Many high-state non-magnetic cataclysmic variables (CVs) exhibit blueshifted absorption or P-Cygni profiles associated with ultraviolet (UV) resonance lines. These features imply the existence of powerful accretion disc winds in CVs. Here, we use our Monte Carlo ionization and radiative transfer code to investigate whether disc wind models that produce realistic UV line profiles are also likely to generate observationally significant recombination line and continuum emission in the optical waveband. We also test whether outflows may be responsible for the single-peaked emission line profiles often seen in high-state CVs and for the weakness of the Balmer absorption edge (relative to simple models of optically thick accretion discs). We find that a standard disc wind model that is successful in reproducing the UV spectra of CVs also leaves a noticeable imprint on the optical spectrum, particularly for systems viewed at high inclination. The strongest optical wind-formed recombination lines are H alpha and He ii lambda 4686. We demonstrate that a higher density outflow model produces all the expected H and He lines and produces a recombination continuum that can fill in the Balmer jump at high inclinations. This model displays reasonable verisimilitude with the optical spectrum of RW Trianguli. No single-peaked emission is seen, although we observe a narrowing of the double-peaked emission lines from the base of the wind. Finally, we show that even denser models can produce a single-peaked H alpha line. On the basis of our results, we suggest that winds can modify, and perhaps even dominate, the line and continuum emission from CVs.

Two-temperature accretion around rotating black holes: a description of the general advective flow paradigm in the presence of various cooling processes to explain low to high luminous sources

We investigate viscous two-temperature accretion disc flows around rotating black holes. We describe the global solution of accretion flows with a sub-Keplerian angular momentum profile, by solving the underlying conservation equations including explicit cooling processes self-consistently. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. We focus on the set of solutions for sub-Eddington, Eddington and super-Eddington mass accretion rates around Schwarzschild and Kerr black holes with a Kerr parameter of 0.998. It is found that the flow, during its infall from the Keplerian to sub-Kepleria transition region to the black hole event horizon, passes through various phases of advection: the general advective paradigm to the radiatively inefficient phase, and vice versa. Hence, the flow governs a much lower electron temperature similar to 10(8)-10(9.5) K, in the range of accretion rate in Eddington units 0.01 less than or similar to (M) over dot less than or similar to 100, compared to the hot protons of temperature similar to 10(10.2)-10(11.8) K. Therefore, the solution may potentially explain the hard X-rays and gamma-rays emitted from active galactic nuclei (AGNs) and X-ray binaries. We then compare the solutions for two different regimes of viscosity. We conclude that a weakly viscous flow is expected to be cooling dominated, particularly at the inner region of the disc, compared to its highly viscous counterpart, which is radiatively inefficient. With all the solutions in hand, we finally reproduce the observed luminosities of the underfed AGNs and quasars (e. g. Sgr A*) to ultraluminous X-ray sources (e. g. SS433), at different combinations of input parameters, such as the mass accretion rate and the ratio of specific heats. The set of solutions also predicts appropriately the luminosity observed in highly luminous AGNs and ultraluminous quasars (e. g. PKS 0743-67).

Temperature profiles of accretion discs around rapidly rotating strange stars in general relativity: A comparison with neutron stars

We compute the temperature profiles of accretion discs around rapidly rotating strange stars, using constant gravitational mass equilibrium sequences of these objects, considering the full effect of general relativity. Beyond a certain critical value of stellar angular momentum (J), we observe the radius ( $r_{\rm orb}$) of the innermost stable circular orbit (ISCO) to increase with J (a property seen neither in rotating black holes nor in rotating neutron stars). The reason for this is traced to the crucial dependence of ${\rm d}r_{\rm orb}/{\rm d}J$ on the rate of change of the radial gradient of the Keplerian angular velocity at $r_{\rm orb}$ with respect to J. The structure parameters and temperature profiles obtained are compared with those of neutron stars, as an attempt to provide signatures for distinguishing between the two. We show that when the full gamut of strange star equation of state models, with varying degrees of stiffness are considered, there exists a substantial overlap in properties of both neutron stars and strange stars. However, applying accretion disc model constraints to rule out stiff strange star equation of state models, we notice that neutron stars and strange stars exclusively occupy certain parameter spaces. This result implies the possibility of distinguishing these objects from each other by sensitive observations through future X-ray detectors.

A simple disc wind model for broad absorption line quasars

Approximately 20 per cent of quasi-stellar objects (QSOs) exhibit broad, blue-shifted absorption lines in their ultraviolet spectra. Such features provide clear evidence for significant outflows from these systems, most likely in the form of accretion disc winds. These winds may represent the ‘quasar’ mode of feedback that is often invoked in galaxy formation/evolution models, and they are also key to unification scenarios for active galactic nuclei (AGN) and QSOs. To test these ideas, we construct a simple benchmark model of an equatorial, biconical accretion disc wind in a QSO and use a Monte Carlo ionization/radiative transfer code to calculate the ultraviolet spectra as a function of viewing angle. We find that for plausible outflow parameters, sightlines looking directly into the wind cone do produce broad, blue-shifted absorption features in the transitions typically seen in broad absorption line (BAL) QSOs. However, our benchmark model is intrinsically X-ray weak in order to prevent overionization of the outflow, and the wind does not yet produce collisionally excited line emission at the level observed in non-BAL QSOs. As a first step towards addressing these shortcomings, we discuss the sensitivity of our results to changes in the assumed X-ray luminosity and mass-loss rate, Ṁwind. In the context of our adopted geometry, Ṁwind ∼ Ṁacc is required in order to produce significant BAL features. The kinetic luminosity and momentum carried by such outflows would be sufficient to provide significant feedback.

Growing pseudo-eigenmodes and positive logarithmic norms in rotating shear flows

Rotating shear flows, when angular momentum increases and angular velocity decreases as functions of radiation coordinate, are hydrodynamically stable under linear perturbation. The Keplerian flow is an example of such a system, which appears in an astrophysical context. Although decaying eigenmodes exhibit large transient energy growth of perturbation which could govern nonlinearity in the system, the feedback of inherent instability to generate turbulence seems questionable. We show that such systems exhibiting growing pseudo-eigenmodes easily reach an upper bound of growth rate in terms of the logarithmic norm of the involved non-normal operators, thus exhibiting feedback of inherent instability. This supports the existence of turbulence of hydrodynamic origin in the Keplerian accretion disc in astrophysics. Hence, this answers the question of the mismatch between the linear theory and experimental/observed data and helps in resolving the outstanding question of the origin of turbulence therein.

Stochastically driven instability in rotating shear flows

The origin of hydrodynamic turbulence in rotating shear flows is investigated, with particular emphasis on the flows whose angular velocity decreases but whose specific angular momentum increases with the increasing radial coordinate. Such flows are Rayleigh stable, but must be turbulent in order to explain the observed data. Such a mismatch between the linear theory and the observations/experiments is more severe when any hydromagnetic/magnetohydrodynamic instability and then the corresponding turbulence therein is ruled out. This work explores the effect of stochastic noise on such hydrodynamic flows. We essentially concentrate on a small section of such a flow, which is nothing but a plane shear flow supplemented by the Coriolis effect. This also mimics a small section of an astrophysical accretion disc. It is found that such stochastically driven flows exhibit large temporal and spatial correlations of perturbation velocities and hence large energy dissipations of perturbation, which presumably generate the instability. A range of angular velocity (Omega) profiles of the background flow, starting from that of a constant specific angular momentum (lambda = Omega r(2); r being the radial coordinate) to a constant circular velocity (v(phi) = Omega r), is explored. However, all the background angular velocities exhibit identical growth and roughness exponents of their perturbations, revealing a unique universality class for the stochastically forced hydrodynamics of rotating shear flows. This work, to the best of our knowledge, is the first attempt to understand the origin of instability and turbulence in three-dimensional Rayleigh stable rotating shear flows by introducing additive noise to the underlying linearized governing equations. This has important implications to resolve the turbulence problem in astrophysical hydrodynamic flows such as accretion discs.

ULTRACAM photometry of the eclipsing cataclysmic variables GY Cnc, IR Com and HT Cas

We present high-speed, three-colour photometry of the eclipsing cataclysmic variables GY Cnc, IR Com and HT Cas. We find that the sharp eclipses in GY Cnc and IR Com are due to eclipses of the white dwarf. There is some evidence for a bright-spot on the edge of the accretion disc in GY Cnc, but not in IR Com. Eclipse mapping of HT Cas is presented which shows changes in the structure of the quiescent accretion disc. Observations in 2002 show the accretion disc to be invisible except for the presence of a bright-spot at the disc edge. 2003 observations, however, clearly show a bright inner disc and the bright-spot to be much fainter than in 2002. Although no outburst was associated with either set of quiescent observations, the system was similar to 0.6 mJy brighter in 2003, mainly due to the enhanced emission from the inner disc. We propose that these changes are due to variations in the mass-transfer rate from the secondary star and through the disc. The disc colours indicate that it is optically thin in both its inner and outer regions. We estimate the white dwarf temperature of HT Cas to be 15 000 +/- 1000 K in 2002 and 14 000 +/- 1000 K in 2003.

The component masses of the cataclysmic variable V347 Puppis

We present time-resolved spectroscopy and photometry of the double-lined eclipsing cataclysmic variable V347 Pup ( = LB 1800). There is evidence of irradiation on the inner hemisphere of the secondary star, which we correct for using a model to give a secondary-star radial velocity of K R = 198 +/- 5 km s(-1). The rotational velocity of the secondary star in V347 Pup is found to be v sin i = 131 +/- 5 km s(-1) and the system inclination is i = 84degrees.0 +/- 2degrees.3. From these parameters we obtain masses of M-1 = 0.63 +/- 0.04 M for the white dwarf primary and M-2 = 0.52 +/- 0.06 M for the M0.5V secondary star, giving a mass ratio of q = 0.83 +/- 0.05. On the basis of the component masses, and the spectral type and radius of the secondary star in V347 Pup, we find tentative evidence for an evolved companion. V347 Pup shows many of the characteristics of the SW Sex stars, exhibiting single-peaked emission lines, high-velocity S-wave components and phase-offsets in the radial velocity curve. We find spiral arms in the accretion disc of V347 Pup and measure the disc radius to be close to the maximum allowed in a pressureless disc.

Evidence for relativistic features in the X-ray spectrum of Mrk 335

We present an analysis of hard X-ray features in the spectrum of the bright Sy 1 galaxy Mrk 335 observed by the XMM-Newton satellite. Our analysis confirms the presence of a broad, ionized Fe Ka emission line in the spectrum, first found by Gondoin et al. The broad line can be modelled successfully by relativistic accretion disc reflection models. This interpretation is unusually robust in the case of Mrk 335 because of the lack of any ionized ('warm') absorber and the absence a clear narrow core to the line. Partial covering by neutral gas cannot, however, be ruled out statistically as the origin of the broad residuals. Regardless of the underlying continuum we report, for the first time in this source, the detection of a narrow absorption feature at the rest frame energy of ~5.9 keV. If the feature is identified with a resonance absorption line of iron in a highly ionized medium, the redshift of the line corresponds to an inflow velocity of ~0.11-0.15c. We present a simple model for the inflow, accounting approximately for relativistic and radiation pressure effects, and use Monte Carlo methods to compute synthetic spectra for qualitative comparison with the data. This modelling shows that the absorption feature can plausibly be reproduced by infalling gas providing that the feature is identified with Fe xxvi. We require the inflowing gas to extend over a limited range of radii at a few tens of r to match the observed feature. The mass accretion rate in the flow corresponds to 60 per cent of the Eddington limit, in remarkable agreement with the observed rate. The narrowness of the absorption line tends to argue against a purely gravitational origin for the redshift of the line, but given the current data quality we stress that such an interpretation cannot be ruled out. © 2006 The Authors.

On the relativistic iron line and soft excess in the Seyfert 1 galaxy Markarian 335

We report a 133-ks XMM-Newton observation of the Seyfert 1 galaxy Markarian 335. The 0.4-12 keV spectrum contains an underlying power-law continuum, a soft excess below 2 keV, and a double-peaked iron emission feature in the 6-7 keV range. We investigate the possibility that the double-peaked emission might represent the characteristic signature of the accretion disc. Detailed investigations show that a moderately broad accretion disc line is most likely present, but that the peaks may be due to narrower components from more distant material. The peaks at 6.4 and 7 keV can be identified, respectively, with the molecular torus in active galactic nucleus unification schemes, and very highly ionized, optically thin gas filling the torus. The X-ray variability spectra on both long (~100 ks) and short (~1 ks) time-scales do not support the recent suggestion that the soft excess is an artefact of variable, moderately ionized absorption. © 2007 The Authors. Journal compilation © 2007 RAS.

The spotty donor star in the X-ray transient Cen X-4

We accurately determine the fundamental system parameters of the neutron star X-ray transient Cen X-4 solely using phase-resolved high-resolution UV-Visual Echelle Spectrograph spectroscopy. We first determine the radial-velocity curve of the secondary star and then model the shape of the phase-resolved absorption line profiles using an X-ray binary model. The model computes the exact rotationally broadened, phase-resolved spectrum and does not depend on assumptions about the rotation profile, limb-darkening coefficients and the effects of contamination from an accretion disc. We determine the secondary star-to-neutron star binary mass ratio to be 0.1755 ± 0.0025, which is an order of magnitude more accurate than previous estimates. We also constrain the inclination angle to be 32^{+8}_{-2} degrees. Combining these values with the results of the radial-velocity study gives a neutron star mass of 1.94^{+0.37}_{-0.85}M⊙ consistent with previous estimates. Finally, we perform the first Roche tomography reconstruction of the secondary star in an X-ray binary. The tomogram reveals surface inhomogeneities that are due to the presence of cool starspots. A large cool polar spot, similar to that seen in Doppler images of rapidly rotating isolated stars, is present on the Northern hemisphere of the K7 secondary star and we estimate that ~4 percent of the total surface area of the donor star is covered with spots.This evidence for starspots supports the idea that magnetic braking plays an important role in the evolution of low-mass X-ray binaries.