3D simulations of rotationally-induced line variability from a classical T Tauri star with a misaligned magnetic dipole

We present 3D simulations of rotationally induced line variability arising from complex circumstellar environment of classical T Tauri stars (CTTS) using the results of the 3D magnetohydrodynamics (MHD) simulations of Romanova et al., who considered accretion on to a CTTS with a misaligned dipole magnetic axis with respect to the rotational axis. The density, velocity and temperature structures of the MHD simulations are mapped on to the radiative transfer grid, and corresponding line source function and the observed profiles of neutral hydrogen lines (Hβ, Paβ and Brγ) are computed using the Sobolev escape probability method. We study the dependency of line variability on inclination angles ( i ) and magnetic axis misalignment angles (Θ). We find the line profiles are relatively insensitive to the details of the temperature structure of accretion funnels, but are influenced more by the mean temperature of the flow and its geometry. By comparing our models with the Paβ profiles of 42 CTTS observed by Folha & Emerson, we find that models with a smaller misaligngment angle  (Θ < ∼15°)  are more consistent with the observations which show that majority of Paβ are rather symmetric around the line centre. For a high inclination system with a small dipole misalignment angle  (Θ≈ 15°)  , only one accretion funnel (on the upper hemisphere) is visible to an observer at any given rotational phase. This can cause an anticorrelation of the line equivalent to the width in the blue wing  ( v < 0)  and that in the red wing  ( v > 0)  over half of a rotational period, and a positive correlation over the other half. We find a good overall agreement of the line variability behaviour predicted by our model and those from observations.     

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.     

The eccentricities of the barium stars

We investigate the eccentricities of barium (Ba  ii ) stars formed via a stellar wind accretion model. We carry out a series of Monte Carlo simulations using a rapid binary evolution algorithm, which incorporates full tidal evolution, mass loss and accretion, and nucleosynthesis and dredge-up on the thermally pulsing asymptotic giant branch. We follow the enhancement of barium in the envelope of the accreting main-sequence companion and dilution into its convective envelope once the star ascends the giant branch.
The observed eccentricities of Ba  ii stars are significantly smaller than those of an equivalent set of normal red giants but are nevertheless non-zero. We show that such a distribution of eccentricities is consistent with a wind accretion model for Ba  ii star production with weak viscous tidal dissipation in the convective envelopes of giant stars. We successfully model the distribution of orbital periods and the number of observed Ba  ii stars. The actual distribution of eccentricities is quite sensitive to the strength of the tides, so that we are able to confirm that this strength is close to, but less than, what is expected theoretically and found with alternative observational tests. Two systems – one very short-period but eccentric, and one long-period and highly eccentric – still lie outside the envelope of our models, and so require a more exotic formation mechanism. All our models, even those which were a good fit to the observed distributions, overproduced the number of high-period barium stars, a problem that could not be solved by some combination of the three parameters: tidal strength, tidal enhancement and wind accretion efficiency.     

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.     

Variable accretion and emission from the stellar winds in the Galactic Centre

We present numerical simulations of stellar wind dynamics in the central parsec of the Galactic Centre, studying in particular the accretion of gas on to Sgr A*, the supermassive black hole. Unlike our previous work, here we use state-of-the-art observational data on orbits and wind properties of individual wind-producing stars. Since wind velocities were revised upwards and non-zero eccentricities were considered, our new simulations show fewer clumps of cold gas and no conspicuous disc-like structure. The accretion rate is dominated by a few close 'slow-wind stars' ( v _w≤ 750 km s⁻¹), and is consistent with the Bondi estimate, but variable on time-scales of tens to hundreds of years. This variability is due to the stochastic infall of cold clumps of gas, as in earlier simulations, and to the eccentric orbits of stars. The present models fail to explain the high luminosity of Sgr A* a few hundred years ago implied by Integral observations, but we argue that the accretion of a cold clump with a small impact parameter could have caused it. Finally, we show the possibility of constraining the total mass-loss rate of the 'slow-wind stars' using near infrared observations of gas in the central few arcseconds.     

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

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

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.     

Accretion to magnetized stars through the Rayleigh–Taylor instability: global 3D simulations

We present results of 3D simulations of magnetohydrodynamics (MHD) instabilities at the accretion disc–magnetosphere boundary. The instability is Rayleigh–Taylor, and develops for a fairly broad range of accretion rates and stellar rotation rates and magnetic fields. It manifests itself in the form of tall, thin tongues of plasma that penetrate the magnetosphere in the equatorial plane. The shape and number of the tongues changes with time on the inner disc dynamical time-scale. In contrast with funnel flows, which deposit matter mainly in the polar region, the tongues deposit matter much closer to the stellar equator. The instability appears for relatively small misalignment angles, Θ≲ 30°, between the star's rotation and magnetic axes, and is associated with higher accretion rates. The hotspots and light curves during accretion through instability are generally much more chaotic than during stable accretion. The unstable state of accretion has possible implications for quasi-periodic oscillations and intermittent pulsations from accreting systems, as well as planet migration.     

Indirect imaging of the accretion stream in eclipsing polars – III. HU Aquarii low state

We apply our technique for indirect imaging of the accretion stream to the polar HU Aqr, using eclipse profiles observed when the system was in a low-accretion state. The eclipse profile is different from that in the high state, and more variable from cycle to cycle. We find that the stream maps are brightest near the white dwarf and there is no significant brightening in the threading region. In the low state the stream threads on to the magnetic field closer to the L₁ point than in the high state, with a footpoint of the accreting field line at high latitude. We then produce maps of the accretion region from polarimetry using Stokes imaging. These show that the majority of the accretion occurs near the equator. The difference between the maps may be explained if most of the stream material is not emitting significantly in the low state. If so, neither the stream eclipse mapping nor Doppler tomography techniques will trace the bulk of the accretion flow between the two stars.     

Indirect imaging of the accretion stream in eclipsing polars – II. HU Aquarii

We apply our technique for indirect imaging of the accretion stream to the polar HU Aquarii, using eclipse profiles observed when the system was in a high accretion state. The accretion stream is relatively luminous, contributing as much as the accretion region on the white dwarf, or more, to the overall system brightness. We model the eclipse profiles using a model stream consisting of a ballistic trajectory from the L1 point followed by a magnetically channelled trajectory that follows a dipole field line out of the orbital plane. We perform model fits using two geometries: a stream that accretes on to both footpoints of the field line, and a stream that accretes only on to the footpoint of the field line above the orbital plane. The stream images indicate that the distribution of emission along the stream is not a simple function of the radial distance from the white dwarf. The stream is redirected by the magnetic field of the white dwarf at a distance 1.0–1.3×10¹⁰ cm from the white dwarf; this implies a mass transfer rate in the range 8–76×10¹⁶ g s⁻¹. The absorption dips in the light curve indicate that the magnetically entrained part of the stream moves from 42° to 48° from the line of centres over the three nights of observation. This is in close agreement with the results of the one-footpoint models, suggesting that this is the more appropriate geometry for these data. The stream images show that, in almost all sections of the stream, the flux peaks in B and is successively fainter in U , V and R .     

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.     

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)     

Hα flares from V404 Cyg in quiescence

We present a spectrophotometric study of short-term optical variability in the quiescent black hole X-ray transient V404 Cyg. This includes two nights of high-time-resolution H α spectroscopy with which we resolve much of the time variability, and a further six nights of archival spectroscopy with lower time resolution but higher spectral resolution. We find significant variability in most of the data considered, with both the H α line and the continuum often varying in a correlated way. This includes both dramatic flares lasting a few hours in which the line flux nearly doubles and lower-level flickering. The strongest flares involve development of asymmetry in the line profile, with the red wing usually strongest independent of orbital phase. It is unclear why this is the case, but we discuss several possible explanations. We consider the energetics of the flares and compare with plausible models including chromospheric activity on the companion star, local magnetic reconnection events within the disc and varying irradiation from close to the black hole. Based on the line profile changes during the flares, we conclude that the most likely origin for the variability is variable photoionization by the central source, although local flares within the disc cannot be ruled out.     

Two-temperature accretion flows in magnetic cataclysmic variables: structures of post-shock emission regions and X-ray spectroscopy

We use a two-temperature hydrodynamical formulation to determine the temperature and density structures of the post-shock accretion flows in magnetic cataclysmic variables (mCVs) and calculate the corresponding X-ray spectra. The effects of two-temperature flows are significant for systems with a massive white dwarf and a strong white-dwarf magnetic field. Our calculations show that two-temperature flows predict harder keV spectra than one-temperature flows for the same white-dwarf mass and magnetic field. This result is insensitive to whether the electrons and ions have equal temperature at the shock, but depends on the electron–ion exchange rate, relative to the rate of radiative loss along the flow. White-dwarf masses obtained by fitting the X-ray spectra of mCVs using hydrodynamic models including the two-temperature effects will be lower than those obtained using single-temperature models. The bias is more severe for systems with a massive white dwarf.     

Comptonization and QPO origins in accreting neutron star systems

We develop a simple, time-dependent Comptonization model to probe the origins of spectral variability in accreting neutron star systems. In the model, soft 'seed photons' are injected into a corona of hot electrons, where they are Compton upscattered before escaping as hard X-rays. The model describes how the hard X-ray spectrum varies when the properties of either the soft photon source or the Comptonizing medium undergo small oscillations. Observations of the resulting spectral modulations can determine whether the variability is due to (i) oscillations in the injection of seed photons, (ii) oscillations in the coronal electron density, or (iii) oscillations in the coronal energy dissipation rate. Identifying the origin of spectral variability should help clarify how the corona operates and its relation to the accretion disc. It will also help in finding the mechanisms underlying the various quasi-periodic oscillations (QPOs) observed in the X-ray outputs of many accreting neutron star and black hole systems. As a sample application of our model, we analyse a kilohertz QPO observed in the atoll source 4U 1608–52. We find that the QPO is driven predominantly by an oscillation in the electron density of the Comptonizing gas.     

The accretion flow in the discless intermediate polar V2400 Ophiuchi

RXTE observations confirm that the X-ray light curve of V2400 Oph is pulsed at the beat cycle, as expected in a discless intermediate polar. There are no X-ray modulations at the orbital or spin cycles, but optical line profiles vary with all three cycles. We construct a model for line-profile variations in a discless accretor, based on the idea that the accretion stream flips from one magnetic pole to the other, and show that this accounts for the observed behaviour over the spin and beat cycles. The minimal variability over the orbital cycle implies that (1) V2400 Oph is at an inclination of only ≈10°, and (2) much of the accretion flow is not in a coherent stream, but is circling the white dwarf, possibly as a ring of denser, diamagnetic blobs. We discuss the light that this sheds on disc formation in intermediate polars.     

Global 3‐D solar‐type star‐disc dynamo systems: I. MHD modeling

We have carried out global three‐dimensional magnetohydrodynamic simulations of the star‐disc interaction region around a young solar‐type star. The magnetic field is generated and maintained by dynamos in the star as well as in the disc. The developing mass flows possess non‐periodic time‐variable azimuthal structure and are controlled by the nonaxisymmetric magnetic fields. Since the stellar field drives a strong stellar wind, accretion is anti‐correlated with the stellar field strength and disc matter is spiraling onto the star at low latitudes, both contrary to the generally assumed accretion picture. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

On the variability and spectral distortion of fluorescent iron lines from black hole accretion discs

We investigate the properties of fluorescent iron lines that arise as a result of the illumination of a black hole accretion disc by an X-ray source located above the disc's surface. We study in detail the light-bending model of the variability of the lines, extending previous work on the subject. We indicate that the bending of photon trajectories to the equatorial plane (a distinct property of the Kerr metric) is the most feasible effect underlying the reduced variability of the lines observed in several objects. A model involving an X-ray source with a varying radial distance, located within a few central gravitational radii around a rapidly rotating black hole, close to the disc's surface, may explain both the elongated red wing of the line profile and the complex variability pattern observed in MCG–6-30-15 by XMM–Newton . We also point out that illumination by radiation that returns to the disc (following the previous reflection) contributes significantly to the formation of the line profile in some cases. As a result of this effect, the line profile always has a pronounced blue peak (which is not observed in the deep minimum state in MCG–6-30-15), unless the reflecting material is absent within the innermost 2–3 gravitational radii.