Light curve morphology study of UW CrB – evidence for a 5 d superorbital period

Since its discovery in 1990, UW CrB (also known as MS1603+2600) has remained a peculiar source without firm classification. Our current understanding is that it is an accretion disc corona (ADC) low-mass X-ray binary. In this paper, we present results from our photometric campaign dedicated to studying the changing morphology of the optical light curves. We find that the optical light curves show remarkable evidence for strongly evolving light curve shapes. In addition, we find that these changes show a modulation at a period of ∼5 d. We interpret these changes as either due to strong periodic accretion disc warping or due to other geometrical changes because of disc precession at a period of 5 d. Finally, we have detected 11 new optical bursts, the phase distribution of which supports the idea of a vertically extended asymmetric accretion disc.     

Investigating the structure of the accretion disc in WZ Sge from multiwaveband time-resolved spectroscopic observations – I

We present the first of two papers describing an in-depth study of multiwaveband phase-resolved spectroscopy of the unusual dwarf nova WZ Sge. In this paper we present an extensive set of Doppler maps of WZ Sge covering optical and infrared emission lines, and describe a new technique for studying the accretion discs of cataclysmic variables using ratioed Doppler maps. Applying the ratioed Doppler map technique to our WZ Sge data shows that the radial temperature profile of the disc is unlike that predicted for a steady state α disc. Time-averaged spectra of the accretion disc line flux (with the bright spot contribution removed) show evidence in the shapes of the line profiles for the presence of shear broadening in a quiescent non-turbulent accretion disc. From the positions of the bright spots in the Doppler maps of different lines, we conclude that the bright spot region is elongated along the ballistic stream, and that the density of the outer disc is low. The velocity of the outer edge of the accretion disc measured from the H α line is found to be 723±23 km s⁻¹. Assuming that the accretion disc reaches to the 3:1 tidal resonance radius, we derive a value for the primary star mass of 0.82 M_⊙. We discuss the implications of our results on the present theories of WZ Sge type dwarf nova outbursts.     

Stellar parameters of young brown dwarfs

We present X‐shooter observations of two brown dwarf candidates. We focus on the determination of stellar parameters and their errors. The targets, an accreting class II and a non‐accreting class III objects, are members of the σ Orionis star‐forming region. We derive the spectroscopic spectral types from the VIS spectrum and the stellar parameters. We find that the uncertainties on the stellar parameters have a minor effect on the determination of the mass accretion rate for the accreting star, thus confirming that the discrepancies between the mass accretion rate estimates found with different (simultaneous) tracers are probably due to different physical conditions where the accretion/wind indicators are produced (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Relation between blue/ultraviolet continuum shape and the ratio of radio-to-optical emission for B3-VLA quasars

The low-frequency radio luminosity is believed to be an indicator of jet power, while the optical/ultraviolet (UV) emission is probably from accretion discs in the nuclei of steep-spectrum radio quasars. We present a correlation between the ratio of radio-to-optical luminosities and the continuum spectral index in blue/UV bands, which might indicate that the continuum shape in blue/UV bands is related to the ratio of jet to accretion power. The results may imply that the spectra and structure of accretion discs are probably affected by the interactions between jets and discs.     

Investigating the structure of the accretion disc in WZ Sge from multiwaveband time-resolved spectroscopic observations – II

We present our second paper describing multiwaveband time-resolved spectroscopy of WZ Sge. We analyse the evolution of both optical and IR emission lines throughout the orbital period and find evidence, in the Balmer lines, for an optically thin accretion disc and an optically thick hotspot. Optical and IR emission lines are used to compute radial velocity curves. Fits to our radial velocity measurements give an internally inconsistent set of values for K ₁, γ and the phase of red-to-blue crossing. We present a probable explanation for these discrepancies, and provide evidence for similar behaviour in other short orbital period dwarf novae. Selected optical and IR spectra are measured to determine the accretion disc radii. Values for the disc radii are found to be strongly dependent on the assumed WD mass and binary orbital inclination. However, the separation of the peaks in the optical emission line (i.e., an indication of the outer disc radius) has been found to be constant during all phases of the supercycle period over the last 40 years.     

Observations of the SW Sextantis star UU Aquarii

We present 14 nights of medium resolution (1–2 Å) spectroscopy of the eclipsing cataclysmic variable UU Aquarii obtained during a high accretion state in 1995 August–October. UU Aqr appears to be an SW Sextantis (SW Sex) star, as noted by Baptista, Steiner & Horne, and we discuss its spectroscopic behaviour in the context of the SW Sex phenomenon. Emission-line equivalent width curves, Doppler tomography, and line profile simulation provide evidence for the presence of a bright spot at the impact site of the accretion stream with the edge of the disc, and a non-axisymmetric, vertically and azimuthally extended absorbing structure in the disc. The absorption has maximum depth in the emission lines around orbital phase 0.8, but is present from φ≈0.4 to φ≈0.95. An origin is explored for this absorbing structure (as well as for the other spectroscopic behaviour of UU Aqr) in terms of the explosive impact of the accretion stream with the disc.     

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

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

Magnetic field and radius of the innermost stable circular orbit near super‐massive black holes in active galactic nuclei

Magnetic fields in an accretion disk around the central black hole can modify the size of the innermost stable circular orbit (ISCO) and can produce a difference to the classical Novikov‐Thorne radius. We estimated the ISCO magnetic field strength from the polarimetric observations of the accretion‐disk radiation. This estimate is obtained taking into account the effect of the Faraday rotation of the polarization plane at the distance of the mean free path of photons between successive electron scattering events. We present the new method for estimating the ISCO radius in the accretion disk, i.e. in the nearest vicinity of a central black hole. Our estimates confirmed the Frolov, Shoom & Tzounis (2014) and Ranea‐Sandoval & Garcia (2015) conclusion that the magnetic field in the accretion disk decreases the size of the innermost stable circular orbit. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Iron line profiles and self-shadowing from relativistic thick accretion discs

We present Fe Kα line profiles from and images of relativistic discs with finite thickness around a rotating black hole using a novel code. The line is thought to be produced by iron fluorescence of a relatively cold X-ray-illuminated material in the innermost parts of the accretion disc and provides an excellent diagnostic of accretion flows in the vicinity of black holes. Previous studies have concentrated on the case of a thin, Keplerian accretion disc. This disc must become thicker and sub-Keplerian with increasing accretion rates. These can affect the line profiles and in turn can influence the estimation of the accretion disc and black hole parameters from the observed line profiles. We here embark on, for the first time, a fully relativistic computation which offers key insights into the effects of geometrical thickness and the sub-Keplerian orbital velocity on the line profiles. We include all relativistic effects such as frame-dragging, Doppler boost, time dilation, gravitational redshift and light bending. We find that the separation and the relative height between the blue and red peaks of the line profile diminish as the thickness of the disc increases. This code is also well suited to produce accretion disc images. We calculate the redshift and flux images of the accretion disc and find that the observed image of the disc strongly depends on the inclination angle. The self-shadowing effect appears remarkable for a high inclination angle, and leads to the black hole shadow being in this case, completely hidden by the disc itself.     

Return mapping of phases and the analysis of the gravitational clustering hierarchy

We present spectroscopic and high-speed photometric data of the eclipsing polar V895 Cen. We find that the eclipsed component is consistent with it being the accretion regions on the white dwarf. This is in contrast to Stobie et al. who concluded that the eclipsed component was not the white dwarf. Further, we find no evidence for an accretion disc in our data. From our Doppler tomography results, we find that the white dwarf has   M ≳0.7 M_⊙  . Our indirect imaging of the accretion stream suggests that the stream is brightest close to the white dwarf. When we observed V895 Cen in its highest accretion state, emission was concentrated along field lines leading to the upper pole. There is no evidence for enhanced emission at the magnetic coupling region.     

On the accretion mode of the intermediate polar V1025 Centauri

The long white-dwarf spin periods in the magnetic cataclysmic variables EX Hya and V1025 Cen imply that if the systems possess accretion discs then they cannot be in equilibrium. It has been suggested that instead they are discless accretors in which the spin-up torques resulting from accretion are balanced by the ejection of part of the accretion flow back towards the secondary. We present phase-resolved spectroscopy of V1025 Cen aimed at deducing the nature of the accretion flow, and compare this with simulations of a discless accretor. We find that both the conventional disc-fed model and the discless-accretor model have strengths and weaknesses, and that further work is needed before we can decide which applies to V1025 Cen.     

Evolution of Class 0 protostars: models versus observations

The rates at which mass accumulates into protostellar cores can now be predicted in numerical simulations. Our purpose here is to develop methods to compare the statistical properties of the predicted protostars with the observable parameters. This requires (1) an evolutionary scheme to convert numerically derived mass accretion rates into evolutionary tracks and (2) a technique to compare the tracks to the observed statistics of protostars. Here, we use a 3D Kolmogorov–Smirnov test to quantitatively compare model evolutionary tracks and observations of Class 0 protostars.
We find that the wide range of accretion functions and time-scales associated with gravoturbulent simulations naturally overcome difficulties associated with schemes that use a fixed accretion pattern. This implies that the location of a protostar on an evolutionary track does not precisely determine the present age or final accrued mass. Rather, we find that predictions of the final mass for protostars from observed   T _bol– L _bol  values are uncertain by a factor of 2 and that the bolometric temperature is not always a reliable measure of the evolutionary stage. Furthermore, we constrain several parameters of the evolutionary scheme and estimate a lifetime of Class 0 sources of  2–6 × 10⁴ yr  , which is related to the local free-fall time and thus to the local density at the onset of the collapse. Models with Mach numbers smaller than six are found to best explain the observational data. Generally, only a probability of 70 per cent was found that our models explain the current observations. This is caused by not well-understood selection effects in the observational sample and the simplified assumptions in the models.     

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‐shooting Herbig Ae/Be stars: Accretion probed by near‐infrared He I emission

The Herbig Ae/Be stars are intermediate mass pre‐main sequence stars that bridge the gap between the low mass T Tauri stars and the Massive Young Stellar Objects. In this mass range, the acting star forming mechanism switches from magnetically controlled accretion to an as yet unknown mechanism, but which is likely to be direct disk accretion onto the star. We observed a large sample of Herbig Ae/Be stars with X‐shooter to address this issue from a multi‐wavelength perspective. It is the largest such study to date, not only because of the number of objects involved, but also because of the large wavelength coverage from the blue to the near‐infrared. This allows many accretion diagnostics to be studied simultaneously. By correlating the various properties with mass, temperature and age, we aim to determine where and whether the magnetically controlled mass accretion mechanism halts and the proposed direct disk accretion takes over. Here, we will give an overview of the background, present some observations and discuss our initial results. We will introduce a new accretion diagnostic for the research of Herbig Ae/Be stars, the HeI 1.083 μm line (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Towards steady-state solutions for supersonic wind accretion on to gravitating objects

We discuss the results of a numerical simulation of the hydrodynamic Bondi–Hoyle accretion obtained on the basis of a high-resolution numerical scheme with the proper entropy correction procedure necessary to avoid the 'weak' non-monotonicity intrinsic in these schemes. Both the axisymmetric and plane cases are considered. The axisymmetric accretion problem turns out to have steady-state solutions for all determining parameters. Steady-state solutions for the plane accretion are found in certain cases even for sinusoidally perturbed and strongly non-uniform winds. Non-stationary phenomena do exist but they are not very violent and can sometimes be attributed as well to numerical as to physical reasons.     

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.     

Thermal synchrotron radiation and its Comptonization in compact X-ray sources

We investigate the process of synchrotron radiation from thermal electrons at semirelativistic and relativistic temperatures. We find an analytic expression for the emission coefficient for random magnetic fields with an accuracy significantly higher than those derived previously. We also present analytic approximations to the synchrotron turnover frequency, treat Comptonization of self-absorbed synchrotron radiation, and give simple expressions for the spectral shape and the emitted power. We also consider modifications of the above results by bremsstrahlung.
We then study the importance of Comptonization of thermal synchrotron radiation in compact X-ray sources. We first consider emission from hot accretion flows and active coronae above optically thick accretion discs in black hole binaries and active galactic nuclei (AGNs). We find that for plausible values of the magnetic field strength, this radiative process is negligible in luminous sources, except for those with hardest X-ray spectra and stellar masses. Increasing the black hole mass results in a further reduction of the maximum Eddington ratio from this process. Then, X-ray spectra of intermediate-luminosity sources, e.g. low-luminosity AGNs, can be explained by synchrotron Comptonization only if they come from hot accretion flows, and X-ray spectra of very weak sources are always dominated by bremsstrahlung. On the other hand, synchrotron Comptonization can account for power-law X-ray spectra observed in the low states of sources around weakly magnetized neutron stars.     

Two- and three-dimensional numerical simulations of accretion discs in a close binary system

We perform 2D and 3D numerical simulations of an accretion disc in a close binary system using the simplified flux vector splitting (SFS) finite volume method. In our calculations, the gas is assumed to be ideal with γ =1.01, 1.05, 1.1 and 1.2 . The mass ratio of the mass-losing star to the mass-accreting star is unity. Our results show that spiral shocks are formed on the accretion disc in all cases. In 2D calculations we find that the smaller γ is, the more tightly the spiral winds. We observe this trend in 3D calculations as well in a somewhat weaker sense. Mach numbers in our discs are less than 10. These values are lower than the values in observed accretion discs in close binary systems.
Recently, Steeghs, Harlaftis & Horne found the first convincing evidence for spiral structure in the accretion disc of the eclipsing dwarf nova binary IP Pegasi, using the technique known as Doppler tomography. Although the Mach numbers in present calculations are rather low, we may claim that the spiral structure that we discovered in earlier numerical simulations is now found observationally.     

SPH simulations of negative (nodal) superhumps: a parametric study

Negative superhumps in cataclysmic variable systems result when the accretion disc is tilted with respect to the orbital plane. The line of nodes of the tilted disc precesses slowly in the retrograde direction, resulting in a photometric signal with a period slightly less than the orbital period. We use the method of smoothed particle hydrodynamics to simulate a series of models of differing mass ratio and effective viscosity to determine the retrograde precession period and superhump period deficit  ɛ₋  as a function of system mass ratio q . We tabulate our results and present fits to both  ɛ₋  and  ɛ₊  versus q , as well as compare the numerical results with those compiled from the literature of negative superhump observations. One surprising result is that while we find negative superhumps most clearly in simulations with an accretion stream present, we also find evidence for negative superhumps in simulations in which we shut off the mass transfer stream completely, indicating that the origin of the photometric signal is more complicated than previously believed.