Reprocessing of primary X-ray emission in the accretion disk of black hole X-ray binaries (BHXRBs) produces the reflection spectrum, with the characteristic Fe K$\alpha$ fluorescence line. Strong low frequency quasi-periodic oscillations (QPOs) are observed from BHXRBs, and the dependence of QPO properties (e.g., phase-lag) on the inclination angle suggests that the observed QPO may be associated with a geometrical effect, e.g., the precession of the X-ray source due to frame-dragging near the spining black hole. Here, in the scenario of Lense-Thirring precession of the X-ray source, we use a Monte-Carlo simulation of radiative transfer to study the irradiation/reflection and the resultant spectral properties including the Fe K$\alpha$ line, as a function of precession phase (time). We found that the reflection fraction, i.e., the ratio of the incident flux towards the disk and the direct flux towards to the observer at infinity, is modulated with the precession phase, which depends on the truncation radius (i.e., the spectral state in the truncated disk model) and the inclination angle. The Fe K$\alpha$ line profile also changes as the primary X-ray source precesses, with the line luminosity and the flux-weighted centroid energy varying with the precession phase. The periodically modulated 2-10 keV continuum flux could apparently lag the line luminosity in phase, if the truncation radius is small enough for Doppler effects due to disk orbital motion to significantly affect the observed radiation.

Black hole accretion is one of nature's most efficient energy extraction processes. When gas falls in, a significant fraction of its gravitational binding energy is either converted into radiation or flows outwards in the form of black hole-driven jets and disk-driven winds. Recently, the Event Horizon Telescope (EHT), an Earth-size sub-millimetre radio interferometer, captured the first images of M87's black hole (or M87*). These images were analysed and interpreted using general-relativistic magnetohydrodynamics (GRMHD) models of accretion disks with rotation axes aligned with the black hole spin axis. However, since infalling gas is often insensitive to the black hole spin direction, misalignment between accretion disk and black hole spin may be a common occurrence in nature. In this work, we use the general-relativistic radiative transfer (GRRT) code \texttt{BHOSS} to calculate the first synthetic radio images of (highly) tilted disk/jet models generated by our GPU-accelerated GRMHD code \texttt{HAMR}. While the tilt does not have a noticeable effect on the system dynamics beyond a few tens of gravitational radii from the black hole, the warping of the disk and jet can imprint observable signatures in EHT images on smaller scales. Comparing the images from our GRMHD models to the 43 GHz and 230 GHz EHT images of M87, we find that M87 may feature a tilted disk/jet system. Further, tilted disks and jets display significant time variability in the 230 GHz flux that can be further tested by longer-duration EHT observations of M87.

Black hole and neutron star X-ray binary systems routinely show quasi-periodic oscillations (QPOs) in their X-ray flux. Despite being strong, easily measurable signals, their physical origin has long remained elusive. However, recent observational and theoretical work has greatly improved our understanding. Here, we briefly review the basic phenomenology of the different varieties of QPO in both black hole and neutron star systems before focusing mainly on low frequency QPOs in black hole systems, for which much of the recent progress has been made. We describe the detailed statistical properties of these QPOs and review the physical models proposed in the literature, with particular attention to those based on Lense-Thirring precession. This is a relativistic effect whereby a spinning massive object twists up the surrounding spacetime, inducing nodal precession in inclined orbits. We review the theory describing how an accretion flow reacts to the Lense-Thirring effect, including analytic theory and recent numerical simulations. We then describe recent observational tests that provide very strong evidence that at least a certain type of low frequency QPOs are a geometric effect, and good evidence that they are the result of precession. We discuss the possibility of the spin axis of the compact object being misaligned with the binary rotation axis for a large fraction of X-ray binaries, as is required for QPOs to be driven specifically by Lense-Thirring precession, as well as some outstanding gaps in our understanding and future opportunities provided by X-ray polarimeters and/or high throughput X-ray detectors.

General-relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black-hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multi-faceted optimizations that, collectively, accelerate computation by 2-5 orders of magnitude for a wide range of applications. Firstly, it involves a novel implementation of a spherical-polar grid with 3D adaptive mesh refinement that operates in each of the 3 dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-res computational performance. Secondly, we demonstrate that local adaptive time-stepping (LAT) on a logarithmic spherical-polar grid accelerates computation by a factor of $\lesssim10$ compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of $\sim10^9$ zone-cycles-per-second-per-node on 5,400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We demonstrate its computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of $13$,$440\times4$,$608\times8$,$092$ cells, and a total of $\lesssim22$ billion cells and $\sim0.65\times10^8$ timesteps, it is among the largest astrophysical simulations ever performed. We find that frame-dragging by the black hole tears up the disk into two independently precessing sub-disks. The innermost sub-disk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

Accreting black holes show characteristic 'reflection' features in their X-ray spectra, including the iron K$\alpha$ fluorescence line, which result from X-rays radiated by a compact central corona being reprocessed in the accretion disc atmosphere. The observed line profile is distorted by relativistic effects, providing a diagnostic for disc geometry. Nearly all previous X-ray reflection spectroscopy studies have made the simplifying assumption that the disc ionization state is independent of radius in order to calculate the restframe reflection spectrum. However, this is unlikely to be the case in reality, since the irradiating flux should drop off steeply with radius. Here we analyse a Nuclear Spectroscopic Telescope ARray observation of GRS 1915+105 that exhibits strong reflection features. We find that using a self-consistently calculated radial ionization profile returns a better fit than assuming constant ionization. Our results are consistent with the inner disc being radiation pressure dominated, as is expected from the high inferred accretion rate for this observation. We also find that the assumed ionization profile impacts on the best fitting disc inner radius. This implies that the black hole spin values previously inferred for active galactic nuclei and X-ray binaries by equating the disc inner radius with the innermost stable circular orbit may be changed significantly by the inclusion of a self-consistent ionization profile.

We present the results of X-ray spectral-timing analysis of a 90ks XMM-Newton observation of the nearby, broad line, early-type AGN NGC5273. The X-ray spectrum revealed the clear presence of a reflection component at high energies, with a clear signature of a narrow iron line at 6.4 keV, consistent with distant reflection. Applying a relativistic reflection model, we found only marginal evidence for a broader relativistic line component. However, cross-spectral analysis revealed that, between 4 and 6 $\times 10^{-4}$ Hz, the 5-8 keV band lagged the 2-3 keV band, implying reflection of the iron line from material close to the black hole. From the analysis of the lag-energy spectrum, we found a broad, but skewed line with a peak of $\approx$ 1000s at 7.5 keV relative to the continuum, which we interpret as the iron line in the reverberation spectrum from an illuminated accretion disc. From the asymmetry in the shape of lag-energy spectrum, we also found that the source is consistent with having an inclination $\geq 45^\circ$.

Since the birth of X-ray Astronomy, spectacular advances have been seen in the imaging, spectroscopic and timing studies of the hot and violent X-ray Universe, and further leaps forward are expected in the future. On the other hand, polarimetry is very much lagging behind: after the measurements of the Crab Nebula and Scorpius X-1, obtained by OSO-8 in the 70s, no more observations have been performed in the classical X-ray band, even if some interesting results have been obtained in hard X-rays and in soft gamma-rays. The NASA/ASI mission IXPE, scheduled for the launch in 2021, is going to provide for the first time imaging X-ray polarimetry in the 2-8 keV band thanks to its photoelectric polarimeter, coupled with ~25'' angular resolution X-ray mirrors. Its orders of magnitude improvement in sensitivity with respect to the OSO-8 Bragg polarimeter implies scientifically meaningful polarimetric measurements for at least the brightest specimens of most classes of X-ray sources. In 2027, the Chinese-led mission eXTP should also be launched. In addition to timing and spectroscopic instruments, eXTP will have on board photoelectric polarimeters very similar to those of IXPE, but with a total effective area 2-3 times larger. Building on IXPE results, eXTP will increase the number of sources for which significant polarimetric measurements could be obtained. However, further progresses, such as exploring a broader energy range, considering a larger effective area, improving the angular resolution, and performing wide-field polarization measurements, are needed to reach a mature phase for X-ray polarimetry. In the first part of this White Paper we will discuss a few scientific cases in which a next generation X-ray Polarimetry mission can provide significant advances. In the second part, a possible concept for a medium-class Next Generation X-ray Polarimetry (NGXP) mission will be sketched.

Identifying the compact object in ultraluminous X-ray sources (ULXs) has to-date required detection of pulsations or a cyclotron resonance scattering feature (CRSF), indicating a magnetised neutron star. However, pulsations are observed to be transient and it is plausible that accretion onto the neutron star may have suppressed the surface magnetic field such that pulsations and CRSFs will be entirely absent. We may therefore lack direct means to identify neutron star systems whilst we presently lack an effective means by which to identify black hole ULXs. Here we present a possible method for separating the ULX population by assuming the X-ray, mHz quasi-periodic oscillations (QPOs) and day timescale periods/QPOs are associated with Lense-Thirring precession of the inflow and outflowing wind respectively. The precession timescales combined with the temperature of the soft X-ray component produce planes where the accretor mass enters as a free parameter. Depending on the properties of the wind, use of these planes may be robust to a range in the angular momentum (spin) and, for high accretion rates, essentially independent of the neutron star's surface dipole field strength. Our model also predicts the mHz QPO frequency and magnitude of the phase-lag imprinted due to propagation through the optically thick wind; in the case of NGC 5408 X-1 we subsequently infer a black hole mass and high spin. Finally, we note that observing secular QPO evolution over sufficient baselines may indicate a neutron star, as the precession responds to spin-up which is not readily observable for black hole primaries.

I present analytic error formulae for the energy-dependent cross-spectrum and rms spectrum, which are Fourier statistics widely used to probe the rapid X-ray variability observed from accreting compact objects. The new formulae cover the modulus, phase, real and imaginary parts of the cross-spectrum, and are valid for any value of intrinsic coherence between variability in different energy bands. I show that existing error formulae (including that for the phase lag), which are valid for a single cross-spectrum or power spectrum, lead to over-fitting when applied to the energy-dependent cross-spectrum - which consists of cross-spectra between individual energy channels and a common reference band. I also introduce an optimal, unbiased way to define the reference band and an accurate way to calculate the intrinsic coherence between energy bands. I find that the traditional use of the old formulae has likely had a rather benign impact on the literature, but recommend the use of the new formulae in future wherever appropriate, since they are more accurate and are no harder to implement than existing error estimates. A code to implement the new error formulae on observational data is available online.

We propose the development of X-ray interferometry (XRI), to reveal the universe at high energies with ultra-high spatial resolution. With baselines which can be accommodated on a single spacecraft, XRI can reach 100 $\mu$as resolution at 10 \AA (1.2 keV) and 20 $\mu$as at 2 \AA (6 keV), enabling imaging and imaging-spectroscopy of (for example) X-ray coronae of nearby accreting supermassive black holes (SMBH) and the SMBH `shadow'; SMBH accretion flows and outflows; X-ray binary winds and orbits; stellar coronae within ~100 pc and many exoplanets which transit across them. For sufficiently luminous sources XRI will resolve sub-pc scales across the entire observable universe, revealing accreting binary SMBHs and enabling trigonometric measurements of the Hubble constant with X-ray light echoes from quasars or explosive transients. A multi-spacecraft `constellation' interferometer would resolve well below 1 $\mu$as, enabling SMBH event horizons to be resolved in many active galaxies and the detailed study of the effects of strong field gravity on the dynamics and emission from accreting gas close to the black hole.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

We present the first X-ray reverberation mass measurement of a stellar-mass black hole. Accreting stellar-mass and supermassive black holes display characteristic spectral features resulting from reprocessing of hard X-rays by the accretion disc, such as an Fe K$\alpha$ line and a Compton hump. This emission probes of the innermost region of the accretion disc through general relativistic distortions to the line profile. However, these spectral distortions are insensitive to black hole mass, since they depend on disc geometry in units of gravitational radii. Measuring the reverberation lag resulting from the difference in path length between direct and reflected emission calibrates the absolute length of the gravitational radius. We use a relativistic model able to reproduce the behaviour of the lags as a function of energy for a wide range of variability timescales, addressing both the reverberation lags on short timescales and the intrinsic hard lags on longer timescales. We jointly fit the time-averaged spectrum and the real and imaginary parts of the cross-spectrum as a function of energy for a range of Fourier frequencies to Rossi X-ray Timing Exporer data from the X-ray binary Cygnus X-1. We also show that introducing a self-consistently calculated radial ionisation profile in the disc improves the fit, but requires us to impose an upper limit on ionisation profile peak to allow a plausible value of the accretion disc density. This limit leads to a mass value more consistent with the existing dynamical measurement.

We present the publicly available model \textsc{reltrans} that calculates the light-crossing delays and energy shifts experienced by X-ray photons originally emitted close to the black hole when they reflect from the accretion disk and are scattered into our line-of-sight, accounting for all general relativistic effects. Our model is fast and flexible enough to be simultaneously fit to the observed energy-dependent cross-spectrum for a large range of Fourier frequencies, as well as to the time-averaged spectrum. This not only enables better geometric constraints than only modelling the relativistically broadened reflection features in the time-averaged spectrum, but additionally enables constraints on the mass of supermassive black holes in active galactic nuclei and stellar-mass black holes in X-ray binaries. We include a self-consistently calculated radial profile of the disk ionization parameter and properly account for the effect that the telescope response has on the predicted time lags. We find that a number of previous spectral analyses have measured artificially low source heights due to not accounting for the former effect and that timing analyses have been affected by the latter. In particular, the magnitude of the soft lags in active galactic nuclei may have been under-estimated, and the magnitude of lags attributed to thermal reverberation in X-ray binaries may have been over-estimated. We fit \textsc{reltrans} to the lag-energy spectrum of the Seyfert galaxy Mrk 335, resulting in a best fitting black hole mass that is smaller than previous optical reverberation measurements ($\sim 7$ million compared with $\sim14-26$ million $M_\odot$).

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

Luminous active galactic nuclei (AGN) and X-Ray binaries (XRBs) tend to be surrounded by geometrically thin, radiatively cooled accretion discs. According to both theory and observations, these are -- in many cases -- highly misaligned with the black hole spin axis. In this work we present the first general relativistic magnetohydrodynamic simulations of very thin ($h/r \sim 0.015-0.05$) accretion discs around rapidly spinning ($a \sim 0.9$) black holes and tilted by 45-65 degrees. We show that the inner regions of the discs with $h/r \lesssim 0.03$ align with the black hole equator, though at smaller radii than predicted by theoretical work. The inner aligned and outer misaligned disc regions are separated by a sharp break in tilt angle accompanied by a sharp drop in density. We find that frame-dragging by the spinning black hole overpowers the disc viscosity, which is self-consistently produced by magnetized turbulence, tearing the disc apart and forming a rapidly precessing inner sub-disc surrounded by a slowly precessing outer sub-disc. We find that at all tilt values the system produces a pair of relativistic jets. At small distances the jets precess rapidly together with the inner sub-disc, whereas at large distances they partially align with the outer sub-disc and precess more slowly. If the tearing radius can be modeled accurately in future work, emission model independent measurements of black hole spin based on precession-driven quasi-periodic oscillations may become possible.

We investigate aperiodic X-ray flux variability in accreting highly magnetized neutron stars - X-ray pulsars (XRPs). The X-ray variability is largely determined by mass accretion rate fluctuations at the NS surface, which replicate accretion rate fluctuations at the inner radius of the accretion disc. The variability at the inner radius is due to fluctuations arising all over the disc and propagating inwards under the influence of viscous diffusion. The inner radius varies with mean mass accretion rate and can be estimated from the known magnetic field strength and accretion luminosity of XRPs. Observations of transient XRPs covering several orders of magnitude in luminosity give a unique opportunity to study effects arising due to the changes of the inner disc radius. We investigate the process of viscous diffusion in XRP accretion discs and construct new analytical solutions of the diffusion equation applicable for thin accretion discs truncated both from inside and outside. Our solutions are the most general ones derived in the approximation of Newtonian mechanics. We argue that the break observed at high frequencies in the power density spectra of XRPs corresponds to the minimal time scale of the dynamo process, which is responsible for the initial fluctuations. Comparing data from the bright X-ray transient A 0535+26 with our model, we conclude that the time scale of initial variability in the accretion disc is a few times longer than the local Keplerian time scale.

The X-ray emission of neutron stars and black holes presents a rich phenomenology that can lead us to a better understanding of their nature and to address more general physics questions: Does general relativity apply in the strong gravity regime? Is spacetime around black holes described by the Kerr metric? This white paper considers how we can investigate these questions by studying reverberation mapping and quasi-periodic oscillations in accreting systems with a combination of high-spectral and high-timing resolution. In the near future, we will be able to study compact objects in the X-rays in a new way: advancements in transition-edge sensors (TES) technology will allow for electron-volt-resolution spectroscopy combined with nanoseconds-precision timing.

The advent of moderately high-resolution X-ray spectroscopy with Chandra and XMM promised to usher in a new age in the study of neutron stars: we thought we would study neutron stars like stars, with resolved absorption spectra revealing their surface chemical composition and physical conditions (e.g. surface gravity, pressure, temperature). Nature, however, did not cooperate in this endeavor, as observations of neutron stars have not revealed verified atomic absorption lines yet. In the near future, advancements in transition-edge sensors (TES) technology will allow for electron-volt-resolution spectroscopy combined with nanoseconds-precision timing. Combining these detectors with collector optics will also us to study neutron stars in much greater detail by achieving high-energy resolution with much larger collecting areas to uncover even weak spectral features over a wide range of the photon energies. Perhaps we will finally be able to study neutron stars like stars.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

We perform a systematic study of the evolution of the waveform of black hole X-ray binary low-frequency QPOs, by measuring the phase difference between their fundamental and harmonic features. This phase difference has been studied previously for small number of QPO frequencies in individual sources. Here, we present a sample study spanning fourteen sources and a wide range of QPO frequencies. With an automated pipeline, we systematically fit power spectra and calculate phase differences from archival Rossi X-ray Timing Explorer (RXTE) observations. We measure well-defined phase differences over a large range of QPO frequencies for most sources, demonstrating that a QPO for a given source and frequency has a persistent underlying waveform. This confirms the validity of recently developed spectral-timing methods performing phase resolved spectroscopy of the QPO. Furthermore, we evaluate the phase difference as a function of QPO frequency. For Type-B QPOs, we find that the phase difference stays constant with frequency for most sources. We propose a simple jet precession model to explain these constant Type-B QPO phase differences. The phase difference of the Type-C QPO is not constant but systematically evolves with QPO frequency, with the resulting relation being similar for a number of high inclination sources, but more variable for low-inclination sources. We discuss how the evolving phase difference can naturally arise in the framework of precession models for the Type-C QPO, by considering the contributions of a direct and reflected component to the QPO waveform.

In the course of its evolution, a black hole (BH) accretes gas from a wide range of directions. Given a random accretion event, the typical angular momentum of an accretion disc would be tilted by $\sim$60$^\circ$ relative to the BH spin. Misalignment causes the disc to precess at a rate that increases with BH spin and tilt angle. We present the first general-relativistic magnetohydrodynamic (GRMHD) simulations spanning a full precession period of highly tilted (60$^\circ$), moderately thin ($h/r=0.1$) accretion discs around a rapidly spinning ($a\simeq0.9$) BH. While the disc and jets precess in phase, we find that the corona, sandwiched between the two, lags behind by $\gtrsim 10^{\circ}$. For spectral models of BH accretion, the implication is that hard non-thermal (corona) emission lags behind the softer (disc) emission, thus potentially explaining some properties of the hard energy lags seen in Type-C low frequency quasi-periodic oscillations in X-Ray binaries. While strong jets are unaffected by this disc-corona lag, weak jets stall when encountering the lagging corona at distances $r \sim 100$ black hole radii. This interaction may quench large-scale jet formation.

It has recently been discovered that a fraction of ultra-luminous X-ray sources (ULXs) exhibit X-ray pulsations, and are therefore powered by super-Eddington accretion onto magnetized neutron stars (NSs). For typical ULX mass accretion rates ($\gtrsim 10^{19}\,{\rm g\,s^{-1}}$), the inner parts of the accretion disc are expected to be in the supercritical regime, meaning that some material is lost in a wind launched from the disc surface, while the rest forms an optically thick envelope around the NS as it follows magnetic field lines from the inner disc radius to the magnetic poles of the star. The envelope hides the central object from a distant observer and defines key observational properties of ULX pulsars: their energy spectrum, polarization, and timing features. The optical thickness of the envelope is affected by the mass losses from the disc. We calculate the mass loss rate due to the wind in ULX pulsars, accounting for the NS magnetic field strength and advection processes in the disc. We argue that the detection of strong outflows from ULX pulsars can be considered evidence of a relatively weak dipole component of the NS magnetic field. We estimate the influence of mass losses on the optical thickness of the envelope and analyze how the envelope affects broadband aperiodic variability in ULXs. We show that brightness fluctuations at high Fourier frequencies can be strongly suppressed by multiple scatterings in the envelope and that the strength of suppression is determined by the mass accretion rate and geometrical size of the magnetosphere.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

Prevalent around luminous accreting black holes, thin discs are challenging to resolve in numerical simulations. When the disc and black hole angular momentum vectors are misaligned, the challenge becomes extreme, requiring adaptive meshes to follow the disc proper as it moves through the computational grid. With our new high-performance general relativistic magnetohydrodynamic (GRMHD) code H-AMR we have simulated the thinnest accretion disc to date, of aspect ratio $H/R\approx0.03\approx1.7^\circ$ , around a rapidly spinning ($a=0.9375$) black hole, using a cooling function. Initially tilted at $10^\circ$, the disc warps inside $\sim 5r_g$ into alignment with the black hole, where $r_g$ is the gravitational radius. This is the first demonstration of Bardeen-Petterson alignment in MHD with viscosity self-consistently generated by magnetized turbulence. The disc develops a low-density high-viscosity ($\alpha_{\rm eff}\sim1.0$) magnetic-pressure--dominated inner region at $r\lesssim25r_g$ that rapidly empties itself into the black hole. This inner region may in reality, due to thermal decoupling of ions and electrons, evaporate into a radiatively inefficient accretion flow if, as we propose, the cooling time exceeds the accretion time set by the order unity effective viscosity. We furthermore find the unexpected result that even our very thin disc can sustain large-scale vertical magnetic flux on the black hole, which launches powerful relativistic jets that carry $20{-}50\%$ of the accretion power along the angular momentum vector of the outer tilted disc, providing a potential explanation for the origin of jets in radio-loud quasars.

We present the discovery of a low-frequency $\approx 5.7$ Hz quasi-periodic oscillation (QPO) feature in observations of the black hole X-ray binary MAXI J1535-571 in its soft-intermediate state, obtained in September-October 2017 by the Neutron Star Interior Composition Explorer (NICER). The feature is relatively broad (compared to other low-frequency QPOs; quality factor $Q\approx 2$) and weak (1.9% rms in 3-10 keV), and is accompanied by a weak harmonic and low-amplitude broadband noise. These characteristics identify it as a weak Type A/B QPO, similar to ones previously identified in the soft-intermediate state of the transient black hole X-ray binary XTE J1550-564. The lag-energy spectrum of the QPO shows increasing soft lags towards lower energies, approaching 50 ms at 1 keV (with respect to a 3-10 keV continuum). This large phase shift has similar amplitude but opposite sign to that seen in Rossi X-ray Timing Explorer data for a Type B QPO from the transient black hole X-ray binary GX 339-4. Previous phase-resolved spectroscopy analysis of the Type B QPO in GX 339-4 pointed towards a precessing jet-like corona illuminating the accretion disk as the origin of the QPO signal. We suggest that this QPO in MAXI J1535-571 may have the same origin, with the different lag sign depending on the scale height of the emitting region and the observer inclination angle.

Accreting black holes show characteristic reflection features in their X-ray spectrum, including an iron K$\alpha$ line, resulting from hard X-ray continuum photons illuminating the accretion disk. The reverberation lag resulting from the path length difference between direct and reflected emission provides a powerful tool to probe the innermost regions around both stellar-mass and supermassive black holes. Here, we present for the first time a reverberation mapping formalism that enables modeling of energy dependent time lags and variability amplitude for a wide range of variability timescales, taking the complete information of the cross-spectrum into account. We use a pivoting power-law model to account for the spectral variability of the continuum that dominates over the reverberation lags for longer time scale variability. We use an analytic approximation to self-consistently account for the non-linear effects caused by this continuum spectral variability, which have been ignored by all previous reverberation studies. We find that ignoring these non-linear effects can bias measurements of the reverberation lags, particularly at low frequencies. Since our model is analytic, we are able to fit simultaneously for a wide range of Fourier frequencies without prohibitive computational expense. We also introduce a formalism of fitting to real and imaginary parts of our cross-spectrum statistic, which naturally avoids some mistakes/inaccuracies previously common in the literature. We perform proof-of-principle fits to Rossi X-ray Timing Explorer data of Cygnus X-1.

Many statistical properties of X-ray aperiodic variability from accreting compact objects can be explained by the propagating fluctuations model applied to the accretion disc. The mass accretion rate fluctuations originate from variability of viscosity, which arises at every radius and causes local fluctuations of the density. The fluctuations diffuse through the disc and result in local variability of the mass accretion rate, which modulates the X-ray flux from the inner disc in the case of black holes, or from the surface in the case of neutron stars. A key role in the theoretical explanation of fast variability belongs to the description of the diffusion process. The propagation and evolution of the fluctuations is described by the diffusion equation, which can be solved by the method of Green functions. We implement Green functions in order to accurately describe the propagation of fluctuations in the disc. For the first time we consider both forward and backward propagation. We show that (i) viscous diffusion efficiently suppress variability at time scales shorter than the viscous time, (ii) local fluctuations of viscosity affect the mass accretion rate variability both in the inner and the outer parts of accretion disc, (iii) propagating fluctuations give rise not only to hard time lags as previously shown, but also produce soft lags at high frequency similar to those routinely attributed to reprocessing, (iv) deviation from the linear rms-flux relation is predicted for the case of very large initial perturbations. Our model naturally predicts bumpy power spectra.

The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in the disc (implying surface dipole fields that are sub-magnetar in strength), then the large scale-height flow can precess through the Lense-Thirring effect which could provide an explanation for the observed super-orbital periods. By connecting the details of the Lense-Thirring effect with the observed pulsar spin period, we are able to infer the moment-of-inertia and therefore equation-of-state of the neutron star without relying on the inclination of, or distance to the system. We apply our technique to the case of NGC 7793 P13 and demonstrate that stronger magnetic fields imply stiffer equations of state. We discuss the caveats and uncertainties, many of which can be addressed through forthcoming radiative magnetohydrodynamic (RMHD) simulations and their connection to observation.

Gas falling into a black hole (BH) from large distances is unaware of BH spin direction, and misalignment between the accretion disc and BH spin is expected to be common. However, the physics of tilted discs (e.g., angular momentum transport and jet formation) is poorly understood. Using our new GPU-accelerated code H-AMR, we performed 3D general relativistic magnetohydrodynamic simulations of tilted thick accretion discs around rapidly spinning BHs, at the highest resolution to date. We explored the limit where disc thermal pressure dominates magnetic pressure, and showed for the first time that, for different magnetic field strengths on the BH, these flows launch magnetized relativistic jets propagating along the rotation axis of the tilted disc (rather than of the BH). If strong large-scale magnetic flux reaches the BH, it bends the inner few gravitational radii of the disc and jets into partial alignment with the BH spin. On longer time scales, the simulated disc-jet system as a whole undergoes Lense-Thirring precession and approaches alignment, demonstrating for the first time that jets can be used as probes of disc precession. When the disc turbulence is well-resolved, our isolated discs spread out, causing both the alignment and precession to slow down.

Whilst astronomy as a science is historically founded on observations at optical wavelengths, studying the Universe in other bands has yielded remarkable discoveries, from pulsars in the radio, signatures of the Big Bang at submm wavelengths, through to high energy emission from accreting, gravitationally-compact objects and the discovery of gamma-ray bursts. Unsurprisingly, the result of combining multiple wavebands leads to an enormous increase in diagnostic power, but powerful insights can be lost when the sources studied vary on timescales shorter than the temporal separation between observations in different bands. In July 2015, the workshop "Paving the way to simultaneous multi-wavelength astronomy" was held as a concerted effort to address this at the Lorentz Center, Leiden. It was attended by 50 astronomers from diverse fields as well as the directors and staff of observatories and spaced-based missions. This community white paper has been written with the goal of disseminating the findings of that workshop by providing a concise review of the field of multi-wavelength astronomy covering a wide range of important source classes, the problems associated with their study and the solutions we believe need to be implemented for the future of observational astronomy. We hope that this paper will both stimulate further discussion and raise overall awareness within the community of the issues faced in a developing, important field.

We present a timing analysis of threeRossi X-ray Timing Explorer observations of the black hole binary Cygnus X-1 with the propagating mass accretion rate fluctuations model PROPFLUC. The model simultaneously predicts power spectra, time lags, and coherence of the variability as a function of energy. The observations cover the soft and the hard state of the source, and the transition between the two. We find good agreement between model pre-dictions and data in the hard and in the soft state. Our analysis suggests that in the soft state the fluctuations propagate in an optically thin hot flow extending up to large radii above and below a stable optically thick disc. In the hard state, our results are consistent with a truncated disc geometry, where the hot flow extends radially inside the inner radius of the disc. In the transition from soft to hard state, the characteristics of the rapid variability are too complex to be successfully described with PROPFLUC. The surface density profile of the hot flow predicted by our model and the lack of QPOs in the soft and hard state, suggest that the spin of the black hole is aligned with the inner accretion disc and therefore probably with the rotational axis of the binary system.

The upcoming launch of the first space based X-ray polarimeter in $\sim 40$ years will provide powerful new diagnostic information to study accreting compact objects. In particular, analysis of rapid variability of the polarisation degree and angle will provide the opportunity to probe the relativistic motions of material in the strong gravitational fields close to the compact objects, and enable new methods to measure black hole and neutron star parameters. However, polarisation properties are measured in a statistical sense, and a statistically significant polarisation detection requires a fairly long exposure, even for the brightest objects. Therefore, the sub-minute timescales of interest are not accessible using a direct time-resolved analysis of polarisation degree and angle. Phase-folding can be used for coherent pulsations, but not for stochastic variability such as quasi-periodic oscillations. Here, we introduce a Fourier method that enables statistically robust detection of stochastic polarisation variability for arbitrarily short variability timescales. Our method is analogous to commonly used spectral-timing techniques. We find that it should be possible in the near future to detect the quasi-periodic swings in polarisation angle predicted by Lense-Thirring precession of the inner accretion flow. This is contingent on the mean polarisation degree of the source being greater than $\sim 4-5\%$, which is consistent with the best current constraints on Cygnus X-1 from the late 1970s.

Timing properties of black hole X-ray binaries in outburst can be modeled with mass accretion rate fluctuations propagating towards the black hole. Such models predict time lags between energy bands due to propagation delays. First application of a propagating fluctuations model to black hole power spectra showed good agreement with the data. Indeed, hard lags observed from these systems appear to be in agreement with this generic prediction. Our PROPFLUC code allows to simultaneously predict power spectra, time lags, and coherence of the variability as a function of energy. This was successfully applied to Swift data on the black hole MAXIJ1659-152, fitting jointly the power spectra in two energy bands and the cross-spectrum between these two bands. In the current work, we attempt to to model two high signal to noise Rossi X-ray Timing Explorer (RXTE) observations of the black hole XTE J1550-564. We find that neither observation can be adequately explained by the model even when considering, additionally to previous PROPFLUC versions, different propagation speeds of the fluctuations. After extensive exploration of model extensions, we tentatively conclude that the quantitative and qualitative discrepancy between model predictions and data is generic to the propagating fluctuations paradigm. This result encourages further investigation of the fundamental hypotheses of the propagating fluctuations model. We discuss some of these hypotheses with an eye to future works.

In June 2015, the black hole X-ray binary (BHXRB) V404 Cygni went into outburst for the first time since 1989. Here, we present a comprehensive search for quasi-periodic oscillations (QPOs) of V404 Cygni during its recent outburst, utilizing data from six instruments on board five different X-ray missions: Swift/XRT, Fermi/GBM, Chandra/ACIS, INTEGRAL's IBIS/ISGRI and JEM-X, and NuSTAR. We report the detection of a QPO at 18 mHz simultaneously with both Fermi/GBM and Swift/XRT, another example of a rare but slowly growing new class of mHz-QPOs in BHXRBs linked to sources with a high orbital inclination. Additionally, we find a duo of QPOs in a Chandra/ACIS observation at 73 mHz and 1.03 Hz, as well as a QPO at 136 mHz in a single Swift/XRT observation that can be interpreted as standard Type-C QPOs. Aside from the detected QPOs, there is significant structure in the broadband power, with a strong feature observable in the Chandra observations between 0.1 and 1 Hz. We discuss our results in the context of current models for QPO formation.

Magnetized neutron stars power at least some ultra-luminous X-ray sources. The accretion flow in these cases is interrupted at the magnetospheric radius and then reaches the surface of a neutron star following magnetic field lines. Accreting matter moving along magnetic field lines forms the accretion envelope around the central object. We show that, in case of high mass accretion rates $\gtrsim 10^{19}\,{\rm g\,s^{-1}}$ the envelope becomes closed and optically thick, which influences the dynamics of the accretion flow and the observational manifestation of the neutron star hidden behind the envelope. Particularly, the optically thick accretion envelope results in a multi-color black-body spectrum originating from the magnetospheric surface. The spectrum and photon energy flux vary with the viewing angle, which gives rise to pulsations characterized by high pulsed fraction and typically smooth pulse profiles. The reprocessing of radiation due to interaction with the envelope leads to the disappearance of cyclotron scattering features from the spectrum. We speculate that the super-orbital variability of ultra-luminous X-ray sources powered by accreting neutron stars can be attributed to precession of the neutron star due to interaction of magnetic dipole with the accretion disc.

We model the energy dependence of a quasi periodic oscillation (QPOs) produced by Lense-Thirring precession of a hot inner flow. We use a fully 3-dimensional Monte-Carlo code to compute the Compton scattered flux produced by the hot inner flow intercepting seed photons from an outer truncated standard disc. The changing orientation of the precessing torus relative to the line of sight produces the observed modulation of the X-ray flux. We consider two scenarios of precession. First, we assume that the precession axis is perpendicular to the plane of the outer disc. In this scenario the relative geometry of the cold disc and the hot torus does not change during precession, so the emitted spectrum does not change, and the modulation is solely due to the changing viewing angle. In the second scenario the precession axis is tilted with respect to the outer disc plane. This leads to changes in the relative geometry of the hot flow and cold plasma, possibly resulting in variations of the plasma temperature and thus generating additional spectral variability, which combines with the variations due to the viewing angle changes.

Quasi-periodic oscillations (QPOs) with frequencies from $\sim0.05$-$30$ Hz are a common feature in the X-ray emission of accreting black hole binaries. As the QPOs originate from the innermost accretion flow, they provide the opportunity to probe the behaviour of matter in extreme gravity. In this paper, we present a systematic analysis of the inclination dependence of phase lags associated with both Type-B and Type-C QPOs in a sample of 15 Galactic black hole binaries. We find that the phase lag at the Type-C QPO frequency strongly depends on inclination, both in evolution with QPO frequency and sign. Although we find that the Type-B QPO soft lags are associated with high inclination sources, the source sample is too small to confirm this as a significant inclination dependence. These results are consistent with a geometrical origin of Type-C QPOs and a different origin for Type-B and Type-C QPOs. We discuss the possibility that the phase lags originate from a pivoting spectral power law during each QPO cycle, while the inclination dependence arises from differences in dominant relativistic effects. We also search for energy dependences in the Type-C QPO frequency. We confirm this effect in the three known sources (GRS 1915+105, H1743-322 and XTE J1550-564) and newly detect it in XTE J1859+226. Lastly, our results indicate that the unknown inclination sources XTE J1859+226 and MAXI J1543-564 are most consistent with a high inclination.

Accreting stellar mass black holes (BHs) routinely exhibit Type-C quasi-periodic oscillations (QPOs). These are often interpreted as Lense-Thirring precession of the inner accretion flow, a relativistic effect whereby the spin of the BH distorts the surrounding space-time, inducing nodal precession. The best evidence for the precession model is the recent discovery, using a long joint XMM-Newton and NuSTAR observation of H 1743-322, that the centroid energy of the iron fluorescence line changes systematically with QPO phase. This was interpreted as the inner flow illuminating different azimuths of the accretion disc as it precesses, giving rise to a blue/red shifted iron line when the approaching/receding disc material is illuminated. Here, we develop a physical model for this interpretation, including a self-consistent reflection continuum, and fit this to the same H 1743-322 data. We use an analytic function to parameterise the asymmetric illumination pattern on the disc surface that would result from inner flow precession, and find that the data are well described if two bright patches rotate about the disc surface. This model is preferred to alternatives considering an oscillating disc ionisation parameter, disc inner radius and radial emissivity profile. We find that the reflection fraction varies with QPO phase (3.5 sigma), adding to the now formidable body of evidence that Type-C QPOs are a geometric effect. This is the first example of tomographic QPO modelling, initiating a powerful new technique that utilizes QPOs in order to map the dynamics of accreting material close to the BH.

The power spectrum of the X-ray fluctuations of accreting black holes often consists of two broad humps. We quantitatively investigate the hypothesis that the lower frequency hump orig- inates from variability in a truncated thin accretion disc, propagating into a large scale-height inner hot flow which, in turn, itself is the origin of the higher frequency hump. We extend the propagating mass accretion rate fluctuations model PROPFLUC to accommodate double hump power spectra in this way. Furthermore, we extend the model to predict the cross-spectrum between two energy bands in addition to their power spectra, allowing us to constrain the model using the observed time lags, which in the model result from both propagation of fluc- tuations from the disc to the hot flow, and inside the hot flow. We jointly fit soft and hard power spectrum, and the cross-spectrum between the two bands using this model for 5 Swift X-ray Telescope observations of MAXI J1659-152. The new double hump model provides a better fit to the data than the old single hump model for most of our observations. The data show only a small phase lag associated with the low frequency hump. We demonstrate quan- titatively that this is consistent with the model. We compare the truncation radius measured from our fits with that measured purely by spectral fitting and find agreement within a factor of two. This analysis encompasses the first joint fits of stellar-mass black hole cross-spectra and power spectra with a single self-consistent physical model.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

Accreting stellar-mass black holes often show a `Type-C' quasi-periodic oscillation (QPO) in their X-ray flux, and an iron emission line in their X-ray spectrum. The iron line is generated through continuum photons reflecting off the accretion disk, and its shape is distorted by relativistic motion of the orbiting plasma and the gravitational pull of the black hole. The physical origin of the QPO has long been debated, but is often attributed to Lense-Thirring precession, a General Relativistic effect causing the inner flow to precess as the spinning black hole twists up the surrounding space-time. This predicts a characteristic rocking of the iron line between red and blue shift as the receding and approaching sides of the disk are respectively illuminated. Here we report on XMM-Newton and NuSTAR observations of the black hole binary H 1743-322 in which the line energy varies systematically over the ~4 s QPO cycle (3.70 sigma significance), as predicted. This provides strong evidence that the QPO is produced by Lense-Thirring precession, constituting the first detection of this effect in the strong gravitation regime. There are however elements of our results harder to explain, with one section of data behaving differently to all the others. Our result enables the future application of tomographic techniques to map the inner regions of black hole accretion disks.

We present a model-independent analysis of the short-timescale energy dependence of low frequency quasi-periodic oscillations (QPOs) in the X-ray flux of GRS 1915+105. The QPO frequency in this source has previously been observed to depend on photon energy, with the frequency increasing with energy for observations with a high ($\gtrsim 2$ Hz) QPO frequency, and decreasing with energy for observations with a low ($\lesssim 2$ Hz) QPO frequency. As this observed energy dependence is currently unexplained, we investigate if it is intrinsic to the QPO mechanism by tracking phase lags on (sub)second timescales. We find that the phase lag between two broad energy bands systematically increases for $5$ - $10$ QPO cycles, after which the QPO becomes decoherent, the phase lag resets and the pattern repeats. This shows that the band with the higher QPO frequency is running away from the other band on short timescales, providing strong evidence that the energy dependence of the QPO frequency is intrinsic. We also find that the faster the QPO decoheres, the faster the phase lag increases, suggesting that the intrinsic frequency difference contributes to the decoherence of the QPO. We interpret our results within a simple geometric QPO model, where different radii in the inner accretion flow experience Lense-Thirring precession at different frequencies, causing the decoherence of the oscillation. By varying the spectral shape of the inner accretion flow as a function of radius, we are able to qualitatively explain the energy-dependent behaviour of both QPO frequency and phase lag.

Black hole binary systems can emit very bright and rapidly varying X-ray signals when material from the companion accretes onto the black hole, liberating huge amounts of gravitational potential energy. Central to this process of accretion is turbulence. In the propagating mass accretion rate fluctuations model, turbulence is generated throughout the inner accretion flow, causing fluctuations in the accretion rate. Fluctuations from the outer regions propagate towards the black hole, modulating the fluctuations generated in the inner regions. Here, I present the theoretical motivation behind this picture before reviewing the array of statistical variability properties observed in the light curves of black hole binaries that are naturally explained by the model. I also discuss the remaining challenges for the model, both in terms of comparison to data and in terms of including more sophisticated theoretical considerations.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

It has long been recognised that quasi-periodic oscillations (QPOs) in the X-ray light curves of accreting black hole and neutron star binaries have the potential to be powerful diagnostics of strong field gravity. However, this potential cannot be fulfilled without a working theoretical model, which has remained elusive. Perhaps the most promising model associates the QPO with Lense-Thirring precession of the inner accretion flow, with the changes in viewing angle and Doppler boosting modulating the flux over the course of a precession cycle. Here, we consider the polarization signature of a precessing inner accretion flow. We use simple assumptions about the Comptonization process generating the emitted spectrum and take all relativistic effects into account, parallel transporting polarization vectors towards the observer along null geodesics in the Kerr metric. We find that both the degree of linear polarization and the polarization angle should be modulated on the QPO frequency. We calculate the predicted absolute rms variability amplitude of the polarization degree and angle for a specific model geometry. We find that it should be possible to detect these modulations for a reasonable fraction of parameter space with a future X-ray polarimeter such as NASA's Polarization Spectroscopic Telescope Array (PolSTAR: the satellite incarnation of the recent balloon experiment X-Calibur).

Doppler tomography can provide a powerful means of determining black hole spin when our view to the central regions are revealed and obscured by optically thick orbiting material, and can provide an independent estimate that does not suffer as many degeneracies as traditional methods. For low mass AGN, time-dependent obscuration is expected to leave a signature in the changing spectrum of the disc emission which extends into the soft X-ray bandpass. We create a spectral model incorporating Doppler tomography and apply it to the case of the low mass (8 $\times$ 10$^{5}$ M$_{\odot}$) AGN, RX J1301.9+2747 which shows unusual timing properties in the form of short-lived flares that we argue are best explained by the orbit of a window through an optically thick wind. Modelling the phase-resolved spectrum over the course of the highest data quality flare indicates a very low spin even when we relax our constraints. This is the lowest mass AGN for which a spin has been measured and the first via this technique. We note that, as the mass and spin are very low, this appears to favour supermassive black hole (SMBH) growth by chaotic rather than constant accretion.

X-ray radiation from black hole binary (BHB) systems regularly displays quasi-periodic oscillations (QPOs). In principle, a number of suggested physical mechanisms can reproduce their power spectral properties, thus more powerful diagnostics which preserve phase are required to discern between different models. In this paper, we first find for two Rossi X-ray Timing Explorer (RXTE) observations of the BHB GRS 1915+105 that the QPO has a well defined average waveform. That is, the phase difference and amplitude ratios between the first two harmonics vary tightly around a well defined mean. This enables us to reconstruct QPO waveforms in each energy channel, in order to constrain QPO phase-resolved spectra. We fit these phase resolved spectra across 16 phases with a model including Comptonisation and reflection (Gaussian and smeared edge components) to find strong spectral pivoting and a modulation in the iron line equivalent width. The latter indicates the observed reflection fraction is changing throughout the QPO cycle. This points to a geometric QPO origin, although we note that the data presented here do not entirely rule out an alternative interpretation of variable disc ionisation state. We also see tentative hints of modulations in the iron line centroid and width which, although not statistically significant, could result from a non-azimuthally symmetric QPO mechanism.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

The relativistic precession model (RPM) can be used to obtain a precise measurement of the mass and spin of a black hole when the appropriate set of quasi periodic oscillations is detected in the power-density spectrum of an accreting black hole. However, in previous studies the solution of the RPM equations could be obtained only through numerical methods at a price of an intensive computational effort. Here we demonstrate that the RPM system of equations can be solved analytically, drastically reducing the computational load, now limited to the Monte-Carlo simulation necessary to estimate the uncertainties. The analytical method not only provides an easy solution to the RPM system when three oscillations are detected, but in all the cases where the detection of two simultaneous oscillations is coupled with an independent mass measurement. We also present a computationally inexpensive method to place limits on the black hole mass and spin when only two oscillations are observed.

We present a spectral and timing analysis of the black hole candidate MAXI J1543-564 during its 2011 outburst. As shown in previous work, the source follows the standard evolution of a black hole outburst. During the rising phase of the outburst we detect an abrupt change in timing behavior associated with the occurrence of a type-B quasi-periodic oscillation (QPO). This QPO and the simultaneously detected radio emission mark the transition between hard and soft intermediate state. We fit power spectra from the rising phase of the outburst using the recently proposed model propfluc. This assumes a truncated disc / hot inner flow geometry, with mass accretion rate fluctuations propagating through a precessing inner flow. We link the propfluc physical parameters to the phenomenological multi-Lorentzian fit parameters. The physical parameter dominating the QPO frequency is the truncation radius, while broad band noise characteristics are also influenced by the radial surface density and emissivity profiles of the flow. In the outburst rise we found that the truncation radius decreases from $r_o \sim 24$ to $10 R_g$, and the surface density increases faster than the mass accretion rate, as previously reported for XTE J1550-564. Two soft intermediate state observations could not be fitted with propfluc, and we suggest that they are coincident with the ejection of material from the inner regions of the flow in a jet or accretion of these regions into the BH horizon, explaining the drop in QPO frequency and suppression of broad band variability preferentially at high energy bands coincident with a radio flare.

The power spectra of black hole binaries have been well studied for decades, giving a detailed phenomenological picture of the variability properties and their correlation with the energy spectrum (spectral state) of the source. Here we take the truncated disc/hot inner flow picture which can describe the spectral changes, and show that propagating mass accretion rate fluctuations in the hot flow can match the broad band power spectral properties seen in black hole binaries, i.e. give approximately band limited noise between a low and high frequency break. The low frequency break marks the viscous timescale at the outer edge of the hot inner flow, which is the inner edge of the truncated disc. The model also predicts the Lense-Thirring precession timescale of the hot flow, as this is set by the surface density of the flow which is self consistently calculated from the propagating fluctuations. We show that this naturally gives the observed relation between the low frequency break and QPO frequency as the outer radius of the flow moves inwards, and that this model predicts many of the observed QPO properties such as correlation of coherence with frequency, and of the recently discovered correlation of frequency with flux on short timescales. We fit this total model of the variability to a sequence of 5 observed power spectra from the bright black hole binary XTE J1550-564 as the source transitioned from a low/hard to very high state. This is the first time that a power spectrum from a black hole binary has been fit with a physical model for the variability. The data are well fit if the inner radius of the flow remains constant, while the outer radius sweeps inwards from ~60-12 Rg. This range of radii agrees with models of the energy spectral evolution, giving the first self consistent description of the evolution of both the spectrum and variability of BHB.

We show that the low frequency QPO seen in the power density spectra of black hole binaries (and neutron stars) can be explained by Lense-Thirring precession. This has been proposed many times in the past, and simple, single radius models can qualitatively match the observed increase in QPO frequency by decreasing a characteristic radius, as predicted by the truncated disc models. However, this also predicts that the frequency is strongly dependent on spin, and gives a maximum frequency at the last stable orbit which is generally much higher than the remarkably constant maximum frequency at ~10Hz observed in all black hole binaries. The key aspect of our model which makes it match these observations is that the precession is of a radially extended region of the hot inner flow. The outer radius is set by the truncation radius of the disc as above, but the inner radius lies well outside of the last stable orbit at the point where numerical simulations show that the density drops off sharply for a misaligned flow. Physically motivated analytic estimates for this inner radius show that it increases with a_*, decreasing the expected frequency in a way which almost completely cancels the expected increase with spin. This ties the maximum predicted frequency to around 10Hz irrespective of a_*, as observed. This is the first QPO model which explains both frequencies and spectrum in the context of a well established geometry for the accretion flow.

Recent observations of accreting black holes reveal the presence of quasi-periodic oscillations (QPO) in the optical power density spectra. The corresponding oscillation periods match those found in the X-rays, implying a common origin. Among the numerous suggested X-ray QPO mechanisms, some may also work in the optical. However, their relevance to the broadband -- optical through X-ray -- spectral properties have not been investigated. For the first time, we discuss the QPO mechanism in the context of the self-consistent spectral model. We propose that the QPOs are produced by Lense-Thirring precession of the hot accretion flow, whose outer parts radiate in the optical wavelengths. At the same time, its innermost parts are emitting the X-rays, explaining the observed connection of QPO periods. We predict that the X-ray and optical QPOs should be either in phase or shifted by half a period, depending on the observer position. We investigate the QPO harmonic content and find that the variability amplitudes at the fundamental frequency are larger in the optical, while the X-rays are expected to have strong harmonics. We then discuss the QPO spectral dependence and compare the expectations to the existing data.

Many statistical properties of the aperiodic variability observed in X-ray radiation from accreting compact objects can be naturally explained by the propagating fluctuations model. This considers variations in mass accretion rate to be stirred up throughout the accretion flow. Variations from the outer regions of the accretion flow will propagate towards the central object, modulating the variations from the inner regions and eventually modulating the radiation, giving rise to the observed linear RMS-flux relation and also Fourier frequency dependent time lags. Previous treatments of this model have relied on computationally intensive Monte Carlo simulations which can only yield an estimate of statistical properties such as the power spectrum. Here, we find exact and analytic expressions for the power spectrum and lag spectrum predicted by the same model. We use our calculation to fit the model of Ingram & Done (2012) to a power spectrum of XTE J1550-564. The result we present here will apply to any treatment of the propagating fluctuations model and thus provides a very powerful tool for future theoretical modelling.

We report on 6 RXTE observations taken during the 2010 outburst of the 11 Hz accreting pulsar IGR J17480-2446 located in the globular cluster Terzan 5. During these observations we find power spectra which resemble those seen in Z-type high-luminosity neutron star low-mass X-ray binaries, with a quasi-periodic oscillation (QPO) in the 35-50 Hz range simultaneous with a kHz QPO and broad band noise. Using well known frequency-frequency correlations, we identify the 35-50 Hz QPOs as the horizontal branch oscillations (HBO), which were previously suggested to be due to Lense-Thirring precession. As IGR J17480-2446 spins more than an order of magnitude more slowly than any of the other neutron stars where these QPOs were found, this QPO can not be explained by frame dragging. By extension, this casts doubt on the Lense-Thirring precession model for other low-frequency QPOs in neutron-star and perhaps even black-hole systems.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

The clear characteristic timescale picked out by the low frequency quasi-periodic oscillations (QPOs) seen in many black hole and neutron star binaries has the potential to provide a very powerful diagnostic of the inner regions of the accretion flow. However, this potential cannot be realised without a quantitative model for the QPO. We have recently shown that the same truncated disc/hot inner flow geometry which is used to interpret the spectral transitions can also directly produce the QPO from Lense-Thirring (vertical) precession of the hot inner flow. This correctly predicts both the frequency and spectrum of the QPO, and the tight correlation of these properties with the total spectrum of the source via a changing truncation radius between the disc and hot flow. This model predicts a unique iron line signature as a vertically tilted flow illuminates different azimuths of the disc as it precesses. The iron line arising from this rotating illumination is blue shifted when the flow irradiates the approaching region of the spinning disc and red shifted when the flow irradiates the receding region of the disc. This gives rise to a characteristic rocking of the iron line on the QPO frequency which is a necessary (and probably sufficient) test of a Lense-Thirring origin. This is also an independent test of disc truncation models for the low/hard state, as vertical precession cannot occur if there is a disc in the midplane. We show that it may be possible to observe this effect using archival data from the Rossi X-ray timing explorer (RXTE) or XMM Newton. However, a clean test requires a combination of moderate resolution and good statistics, such as would be available from a long XMM-Newton observation or with data from the proposed ESA mission LOFT.

Black hole accretion flows show rapid X-ray variability. The Power Spectral Density (PSD) of this is typically fit by a phenomenological model of multiple Lorentzians for both the broad band noise and Quasi-Periodic Oscillations (QPOs). Our previous paper (Ingram & Done 2011) developed the first physical model for the PSD and fit this to observational data. This was based on the same truncated disc/hot inner flow geometry which can explain the correlated properties of the energy spectra. This assumes that the broad band noise is from propagating fluctuations in mass accretion rate within the hot flow, while the QPO is produced by global Lense-Thirring precession of the same hot flow. Here we develop this model, making some significant improvements. Firstly we specify that the viscous frequency (equivalently, surface density) in the hot flow has the same form as that measured from numerical simulations of precessing, tilted accretion flows. Secondly, we refine the statistical techniques which we use to fit the model to the data. We re-analyse the PSD from the 1998 rise to outburst of XTE J1550-564 with our new model in order to assess the impact of these changes. We find that the derived outer radii of the hot flow (set by the inner radius of the truncated disc) are rather similar, changing from ~68-13 Rg throughout the outburst rise. However, the more physical assumptions of our new model also allow us to constrain the scale height of the flow. This decreases as the outer radius of the flow decreases, as expected from the spectral evolution. The spectrum steepens in response to the increased cooling as the as the truncation radius sweeps in, so gas pressure support for the flow decreases. The new model, propfluc, is publically available within the xspec spectral fitting package.

We propose a physical framework for interpreting the characteristic frequencies seen in the broad band power spectra from black hole and neutron star binaries. We use the truncated disc/hot inner flow geometry, and assume that the hot flow is generically turbulent. Each radius in the hot flow produces fluctuations, and we further assume that these are damped on the viscous frequency. Integrating over radii gives broad band continuum noise power between low and high frequency breaks which are set by the viscous timescale at the outer and inner edge of the hot flow, respectively. Lense-Thirring (vertical) precession of the entire hot flow superimposes the low frequency QPO on this continuum power. We test this model on the power spectra seen in the neutron star systems (atolls) as these have the key advantage that the (upper) kHz QPO most likely independently tracks the truncation radius. These show that this model can give a consistent solution, with the truncation radius decreasing from 20-8 Rg while the inner radius of the flow remains approximately constant at ~4.5 Rg i.e. 9.2 km. We use this very constrained geometry to predict the low frequency QPO from Lense-Thirring precession of the entire hot flow from r_o to r_i. The simplest assumption of a constant surface density in the hot flow matches the observed QPO frequency to within 25 per cent. This match can be made even better by considering that the surface density should become increasingly centrally concentrated as the flow collapses into an optically thick boundary layer during the spectral transition. The success of the model opens up the way to use the broad band power spectra as a diagnostic of accretion flows in strong gravity.