Search

We present results of time-series analysis of the first year of the Fairall 9 intensive disc-reverberation campaign. We used Swift and the Las Cumbres Observatory global telescope network to continuously monitor Fairall 9 from X-rays to near-infrared at a daily to subdaily cadence. The cross-correlation function between bands provides evidence for a lag spectrum consistent with the τ ∝ λ4/3 scaling expected for an optically thick, geometrically thin blackbody accretion disc. Decomposing the flux into constant and variable components, the variable component’s spectral energy distribution is slightly steeper than the standard accretion disc prediction. We find evidence at the Balmer edge in both the lag and flux spectra for an additional bound-free continuum contribution that may arise from reprocessing in the broad-line region. The inferred driving light curve suggests two distinct components, a rapidly variable (100 d) component with an opposite lag to the reverberation signal.

We present a detailed X-ray timing analysis of the highly variable NLS1 galaxy, IRAS 13224-3809. The source was recently monitored for 1.5 Ms with XMM-Newton which, combined with 500 ks archival data, makes this the best studied NLS1 galaxy in X-rays to date. We apply standard time- and Fourier-domain in order to understand the underlying variability process. The source flux is not distributed lognormally, as would be expected for accreting sources. The first non-linear rms-flux relation for any accreting source in any waveband is found, with $\mathrm{rms} \propto \mathrm{flux}^{2/3}$. The light curves exhibit significant strong non-stationarity, in addition to that caused by the rms-flux relation, and are fractionally more variable at lower source flux. The power spectrum is estimated down to $\sim 10^{-7}$ Hz and consists of multiple peaked components: a low-frequency break at $\sim 10^{-5}$ Hz, with slope $\alpha < 1$ down to low frequencies; an additional component breaking at $\sim 10^{-3}$ Hz. Using the high-frequency break we estimate the black hole mass $M_\mathrm{BH} = [0.5-2] \times 10^{6} M_{\odot}$, and mass accretion rate in Eddington units, $\dot m_{\rm Edd} \gtrsim 1$. The non-stationarity is manifest in the PSD with the normalisation of the peaked components increasing with decreasing source flux, as well as the low-frequency peak moving to higher frequencies. We also detect a narrow coherent feature in the soft band PSD at $0.7$ mHz, modelled with a Lorentzian the feature has $Q \sim 8$ and an $\mathrm{rms} \sim 3$ %. We discuss the implication of these results for accretion of matter onto black holes., Comment: Accepted to MNRAS. 19 pages, 17 figures

Following a tidal disruption event (TDE), the accretion rate can evolve from quiescent to near-Eddington levels and back over months - years timescales. This provides a unique opportunity to study the formation and evolution of the accretion flow around supermassive black holes (SMBHs). We present two years of multi-wavelength monitoring observations of the TDE AT2018fyk at X-ray, UV, optical and radio wavelengths. We identify three distinct accretion states and two state transitions between them. These appear remarkably similar to the behaviour of stellar-mass black holes in outburst. The X-ray spectral properties show a transition from a soft (thermal-dominated) to a hard (power-law dominated) spectral state around L$_{\rm bol} \sim $few $ \times 10^{-2}$ L$_{\rm Edd}$, and the strengthening of the corona over time $\sim$100--200 days after the UV/optical peak. Contemporaneously, the spectral energy distribution (in particular, the UV-to-X-ray spectral slope $\alpha_{ox}$) shows a pronounced softening as the outburst progresses. The X-ray timing properties also show a marked change, initially dominated by variability at long ($>$day) timescales while a high frequency ($\sim$10$^{-3}$ Hz) component emerges after the transition into the hard state. At late times ($\sim$500 days after peak), a second accretion state transition occurs, from the hard into the quiescent state, as identified by the sudden collapse of the bolometric (X-ray+UV) emission to levels below 10$^{-3.4}$ L$_{\rm Edd}$. Our findings illustrate that TDEs can be used to study the scale (in)variance of accretion processes in individual SMBHs. Consequently, they provide a new avenue to study accretion states over seven orders of magnitude in black hole mass, removing limitations inherent to commonly used ensemble studies., Comment: Accepted version following referee comments. 2 new figures compared to previous arxiv version (Figs 9 and 10). Data will be available from the journal webpages, or upon request to the authors

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon., 19 pages, 9 figures, ApJ in press

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., Comment: White Paper submitted in response to ESA's Voyage 2050 call

An X-ray Interferometer (XRI) has recently been proposed as a theme for ESA's Voyage 2050 planning cycle, with the eventual goal to observe the X-ray sky with an unprecedented angular resolution better than 1 micro arcsec (5 prad) [1]. A scientifically very interesting mission is possible on the basis of a single spacecraft [2], owing to the compact 'telephoto' design proposed earlier by Willingale [3]. Between the practical demonstration of X-ray interferometry at 1 keV by Cash et al. [4] with a 1 mm baseline and 0.1 arcsec effective resolution to a mission flying an interferometer with a baseline of one or more meters, an effective collecting area of square meters and micro arcsec resolution lie many milestones. The first important steps to scale up from a laboratory experiment to a viable mission concept will have to be taken on a scalable and flexible testbed set-up. Such a testbed cannot singularly focus on the optical aspects, but should simultaneously address the thermal and mechanical stability of the interferometer. A particular challenge is the coherent X-ray source, which should provide a wavefront at the entrance of the interferometer that is transversely coherent over a distance at least equal to the baseline, and bright enough. In this paper, we will explore the build-up of a testbed in several stages, with increasing requirements on optical quality and associated thermo-mechanical control and source sophistication, with the intent to guide the technological development of X-ray interferometry from the lab to space in a sequence of achievable milestones., Comment: 12 pages, 4 figures, final draft manuscript for the SPIE Astronomical Telescopes + Instrumentation online conference, 2020

We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$\beta$, C IV, and Ly$\alpha$ broad emission lines. We find an extended disk-like H$\beta$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$\alpha$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$\alpha$ emission arising at smaller radii than the H$\beta$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$\beta$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole., Comment: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

We studied the outburst evolution and timing properties of the recently discovered X-ray transient MAXI J1348-630 as observed with NICER. We produced the fundamental diagrams commonly used to trace the spectral evolution, and power density spectra to study the fast X-ray variability. The main outburst evolution of MAXI J1348-630 is similar to that commonly observed in black hole transients. The source evolved from the hard state, through hard- and soft-intermediate states, into the soft state in the outburst rise, and back to the hard state in reverse during the outburst decay. At the end of the outburst, MAXI J1348-630 underwent two reflares with peak fluxes ~1 and ~2 orders of magnitude fainter than the main outburst, respectively. During the reflares, the source remained in the hard state only, without undergoing any state transitions, which is similar to the so-called "failed outbursts". Different types of quasi-periodic oscillations (QPOs) are observed at different phases of the outburst. Based on our spectral-timing results, we conclude that MAXI J1348-630 is a black hole candidate., 12 pages, 8 figures, 1 table, accepted for publication in MNRAS

We present a detailed X-ray spectral and variability study of the full 2018 outburst of MAXI J1727-203 using NICER observations. The outburst lasted approximately four months. Spectral modelling in the 0.3-10 keV band shows the presence of both a soft thermal and a hard Comptonised component. The analysis of these components shows that MAXI J1727-203 evolved through the soft, intermediate and hard spectral states during the outburst. We find that the soft (disc) component was detected throughout almost the entire outburst, with temperatures ranging from ~0.4 keV, at the moment of maximum luminosity, to ~0.1 keV near the end of the outburst. The power spectrum in the hard and intermediate states shows broadband noise up to 20 Hz, with no evidence of quasi-periodic oscillations. We also study the rms spectra of the broadband noise at 0.3-10 keV of this source. We find that the fractional rms increases with energy in most of the outburst except during the hard state, where the fractional rms remains approximately constant with energy. We also find that, below 3 keV, the fractional rms follows the same trend generally observed at energies >3 keV, a behaviour known from previous studies of black holes and neutron stars. The spectral and timing evolution of MAXI J1727-203, as parametrised by the hardness-intensity, hardness-rms, and rms-intensity diagrams, suggest that the system hosts a black hole, although we could not rule out a neutron star., 15 pages, 9 figures, 1 table, Accepted for publication in MNRAS main journal

We present Neutron Star Interior Composition Explorer X-ray and Arcminute Microkelvin Imager Large Array radio observations of a rapid hard-to-soft state transition in the black hole X-ray transient MAXI J1820+070. During the transition from the hard state to the soft state a switch between two particular types of quasiperiodic oscillations (QPOs) was seen in the X-ray power density spectra, from type-C to type-B, along with a drop in the strength of the broadband X-ray variability and a brief flare in the 7-12 keV band. Soon after this switch (~1.5-2.5 hr) a strong radio flare was observed that corresponded to the launch of superluminal ejecta. Although hints of a connection between QPO transitions and radio flares have been seen in other black hole X-ray transients, our observations constitute the strongest observational evidence to date for a link between the appearance of type-B QPOs and the launch of discrete jet ejections., Comment: Accepted for publication in ApJ Letters

The evolution of massive stars is influenced by the mass lost to stellar winds over their lifetimes. These winds limit the masses of the stellar remnants (such as black holes) that the stars ultimately produce. We use radio astrometry to refine the distance to the black hole X-ray binary Cygnus X-1, which we find to be $2.22^{+0.18}_{-0.17}$ kiloparsecs. When combined with previous optical data, this implies a black hole mass of $21.2\pm2.2$ solar masses, higher than previous measurements. The formation of such a high-mass black hole in a high-metallicity system constrains wind mass loss from massive stars., Published online in Science on 2021 February 18; Main (3 figures; 1 Table) + Supplementary (11 figures; 3 Tables)

X-ray reverberation echoes are assumed to be produced in the strongly distorted spacetime around accreting supermassive black holes. This signal allows us to spatially map the geometry of the inner accretion flow - a region which cannot yet be spatially resolved by any telescope - and provides a direct measure of the black hole mass and spin. The reverberation timescale is set by the light travel path between the direct emission from a hot X-ray corona and the reprocessed emission from the inner edge of the accretion disc. However, there is an inherent degeneracy in the reverberation signal between black hole mass, inner disc radius and height of the illuminating corona above the disc. Here, we use a long X-ray observation of the highly-variable active galaxy, IRAS 13224-3809, to track the reverberation signal as the system evolves on timescales of a day. With the inclusion of all the relativistic effects, modelling reveals that the height of the X-ray corona increases with increasing luminosity, providing a dynamic view of the inner accretion region. This simultaneous modelling allows us to break the inherent degeneracies and obtain an independent timing-based estimate for the mass and spin of the black hole. The uncertainty on black hole mass is comparable to the leading optical reverberation method, making X-ray reverberation a powerful technique, particularly for sources with low optical variability., Comment: Accepted Nat. Ast. version

We present the first short time-scale ($\sim$hours to days) optical variability study of a large sample of Active Galactic Nuclei (AGN) observed with the Kepler/K2 mission. The sample contains 252 AGN observed over four campaigns with $\sim 30$ minute cadence selected from the Million Quasar Catalogue with R magnitude $< 19$. We performed time series analysis to determine their variability properties by means of the power spectral densities (PSDs) and applied Monte Carlo techniques to find the best model parameters that fit the observed power spectra. A power-law model is sufficient to describe all the PSDs of our sample. A variety of power-law slopes were found indicating that there is not a universal slope for all AGN. We find that the rest-frame amplitude variability in the frequency range of $6\times10^{-6}-10^{-4}$ Hz varies from $1-10$ % with an average of 1.7 %. We explore correlations between the variability amplitude and key parameters of the AGN, finding a significant correlation of rest-frame short-term variability amplitude with redshift. We attribute this effect to the known "bluer when brighter" variability of quasars combined with the fixed bandpass of Kepler data. This study also enables us to distinguish between Seyferts and Blazars and confirm AGN candidates. For our study we have compared results obtained from light curves extracted using different aperture sizes and with and without de-trending. We find that limited de-trending of the optimal photometric precision light curve is the best approach, although some systematic effects still remain present., Published in MNRAS 16 February 2018

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., Comment: 52 pages, 15 figures, accepted, invited review (to appear in New Astronomy Reviews), v3: updated figure and text

Active galactic nucleus (AGN) feedback is the process by which supermassive black holes in the centres of galaxies may moderate the growth of their hosts. Gas outflows from supermassive black holes release huge quantities of energy into the interstellar medium, clearing the surrounding gas. The most extreme of these, the ultra-fast outflows (UFOs), are the subset of X-ray detected outflows with velocities higher than 10,000 km/s, believed to originate in relativistic disc winds, a few hundred gravitational radii from the black hole. The absorption features produced by these outflows are variable, but no clear link has been found between the behaviour of the X-ray continuum and the energy or equivalent width of the outflow features due to the long time-scales of quasar variability. Here, we present the detection of multiple absorption lines from an extreme ultra-fast gas flow in the X-ray spectrum of the active galactic nucleus IRAS 13224-3809, at 0.236+/-0.006 times the speed of light (71,000 km/s), where the absorption is strongly anti-correlated with the emission from the inner regions of the accretion disk. If the gas flow is identified as a genuine outflow then it is in the fastest 5 per cent of such winds, and its variability is hundreds of times faster than in other variable winds, allowing us to observe in hours what would take months in a quasar. We find signatures of the wind simultaneously in both low and high energy detectors, which are consistent with a single ionized outflow, linking the two phenomena. The detection of the wind responding to the emission from the inner disk demonstrates a connection between accretion processes occurring on very different scales, with the X-rays from within a few gravitational radii of the black hole ionizing the relativistically outflowing gas as the flux rises., Comment: 13 pages, 7 figures, published in Nature 02/03/17

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., Comment: Accepted for publication in MNRAS

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state., Comment: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

The geometry of the accretion flow around stellar-mass black holes can change on timescales of days to months. When a black hole emerges from quiescence (that is, it "turns on" after accreting material from its companion) it has a very hard (high-energy) X-ray spectrum produced by a hot corona positioned above its accretion disk, and then transitions to a soft (lower-energy) spectrum dominated by emission from the geometrically thin accretion disk, which extends to the innermost stable circular orbit. Much debate persists over how this transition occurs and whether it is driven largely by a reduction in the truncation radius of the disk or by a reduction in the spatial extent of the corona. Observations of X-ray reverberation lags in supermassive black-hole systems suggest that the corona is compact and that the disk extends nearly to the central black hole. Observations of stellar-mass black holes, however, reveal equivalent (mass-scaled) reverberation lags that are much larger, leading to the suggestion that the accretion disk in the hard X-ray state of stellar-mass black holes is truncated at a few hundreds of gravitational radii from the black hole. Here we report X-ray observations of the black-hole transient MAXI J1820+070. We find that the reverberation time lags between the continuum-emitting corona and the irradiated accretion disk are 6 to 20 times shorter than previously seen. The timescale of the reverberation lags shortens by an order of magnitude over a period of weeks, whereas the shape of the broadened iron K emission line remains remarkably constant. This suggests a reduction in the spatial extent of the corona, rather than a change in the inner edge of the accretion disk., Published in Nature January 9, 2019

Swift intensive accretion disk reverberation mapping of four AGN yielded light curves sampled $\sim$200-350 times in 0.3-10 keV X-ray and six UV/optical bands. Uniform reduction and cross-correlation analysis of these datasets yields three main results: 1) The X-ray/UV correlations are much weaker than those within the UV/optical, posing severe problems for the lamp-post reprocessing model in which variations in a central X-ray corona drive and power those in the surrounding accretion disk. 2) The UV/optical interband lags are generally consistent with $ \tau \propto \lambda^{4/3} $ as predicted by the centrally illuminated thin accretion disk model. While the average interband lags are somewhat larger than predicted, these results alone are not inconsistent with the thin disk model given the large systematic uncertainties involved. 3) The one exception is the U band lags, which are on average a factor of $\sim$2.2 larger than predicted from the surrounding band data and fits. This excess appears due to diffuse continuum emission from the broad-line region (BLR). The precise mixing of disk and BLR components cannot be determined from these data alone. The lags in different AGN appear to scale with mass or luminosity. We also find that there are systematic differences between the uncertainties derived by javelin vs. more standard lag measurement techniques, with javelin reporting smaller uncertainties by a factor of 2.5 on average. In order to be conservative only standard techniques were used in the analyses reported herein., Comment: 25 pages, 9 tables, 6 figures, published in ApJ. The paper has been revised to conform to the published version. Please note that the data (Table 2) are available at the Digital Repository at the University of Maryland, https://drum.lib.umd.edu/handle/1903/21536

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment., Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

We report on simultaneous sub-second optical and X-ray timing observations of the low mass X-ray binary black hole candidate MAXI J1820+070. The bright 2018 outburst rise allowed simultaneous photometry in five optical bands ($ugriz_s$) with HiPERCAM/GTC (Optical) at frame rates over 100 Hz, together with NICER/ISS observations (X-rays). Intense (factor of two) red flaring activity in the optical is seen over a broad range of timescales down to $\sim$10 ms. Cross-correlating the bands reveals a prominent anti-correlation on timescales of $\sim$seconds, and a narrow sub-second correlation at a lag of $\approx$+165 ms (optical lagging X-rays). This lag increases with optical wavelength, and is approximately constant over Fourier frequencies of $\sim$0.3-10 Hz. These features are consistent with an origin in the inner accretion flow and jet base within $\sim$5000 Gravitational radii. An additional $\sim$+5 s lag feature may be ascribable to disc reprocessing. MAXI J1820+070 is the third black hole transient to display a clear $\sim$0.1s optical lag, which may be common feature in such objects. The sub-second lag $variation$ with wavelength is novel, and may allow constraints on internal shock jet stratification models., MNRAS Letters, In Press

MAXI J1727-203 is a new X-ray transient discovered on 5 June 2018. A hard-to-soft state transition at the beginning of the outburst led to the identification as a black hole candidate. MAXI J1727-203 was monitored with the Neutron Star Interior Composition Explorer (NICER) on an almost daily basis from the beginning of the outburst. We present a spectral and timing analysis of the full outburst of the source, which lasted approximately four months. A preliminary spectral analysis suggest that the accretion disk component can was detected throughout the entire outburst, with temperatures ranging from ~0.4 keV (in the soft state), down to ~0.2 keV near the end of the outburst when the source was in the hard state. The power spectrum in the hard state shows broadband noise up to 10 Hz, with no detection of any quasi-periodic oscillations. We argue that the system's characteristics are not consistent with those expected for a neutron star and that they are particularly reminiscent of the black hole X-ray binaries XTE J1118+480 and Cyg X-1.

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.

We report on radio and X-ray monitoring observations of the BHC Swift J1753.5-0127 taken over a ~10 year period. Presented are daily radio observations at 15 GHz with the AMI-LA and X-ray data from Swift XRT and BAT. Also presented is a deep 2hr JVLA observation taken in an unusually low-luminosity soft-state (with a low disk temperature). We show that although the source has remained relatively radio-quiet compared to XRBs with a similar X-ray luminosity in the hard-state, the power-law relationship scales as $\zeta=0.96\pm0.06$ i.e. slightly closer to what has been considered for radiatively inefficient accretion disks. We also place the most stringent limit to date on the radio-jet quenching in an XRB soft-state, showing the connection of the jet quenching to the X-ray power-law component; the radio flux in the soft-state was found to be $, Comment: 8 pages, 4 figures, Accepted for publication in MNRAS

We present results from the spectral fitting of the candidate black hole X-ray binary Swift J1753.5-0127 in an accretion state previously unseen in this source. We fit the 0.7-78 keV spectrum with a number of models, however the preferred model is one of a multi-temperature disk with an inner disk temperature $\mathrm{k}T_\mathrm{in}=0.252\pm0.003$ keV scattered into a steep power-law with photon index $\Gamma=6.39^{+0.08}_{-0.02}$ and an additional hard power law tail ($\Gamma=1.79\pm0.02$). We report on the emergence of a strong disk-dominated component in the X-ray spectrum and we conclude that the source has entered the soft state for the first time in its ~10 year prolonged outburst. Using reasonable estimates for the distance to the source ($3$ kpc) and black hole mass ($5M_{\odot}$), we find the unabsorbed luminosity (0.1-100 keV) to be $\approx0.60$% of the Eddington luminosity, making this one of the lowest luminosity soft states recorded in X-ray binaries. We also find that the accretion disk extended towards the compact object during its transition from hard to soft, with the inner radius estimated to be $R_{\mathrm{in}}=28.0^{+0.7}_{-0.4} R_g$ or ~$12R_g$, dependent on the boundary condition chosen, assuming the above distance and mass, a spectral hardening factor $f=1.7$ and a binary inclination $i=55^{\circ}$., Comment: 10 pages, 5 figures, accepted for publication in MNRAS

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., Accepted for publication in ApJ Letters

We present a detailed analysis of the X-ray/IR fast variability of the Black-Hole Transient GX 339-4 during its low/hard state in August 2008. Thanks to simultaneous high time-resolution observations made with the VLT and RXTE, we performed the first characterisation of the sub-second variability in the near-infrared band - and of its correlation with the X-rays - for a low-mass X-ray binary, using both time and frequency-domain techniques. We found a power-law correlation between the X-ray and infrared fluxes when measured on timescales of 16 seconds, with a marginally variable slope, steeper than the one found on timescales of days at similar flux levels. We suggest the variable slope - if confirmed - could be due to the infrared flux being a non-constant combination of both optically thin and optically thick synchrotron emission from the jet, as a result of a variable self-absorption break. From cross spectral analysis we found an approximately constant infrared time lag of $\approx$ 0.1s, and a very high coherence of $\approx$ 90 per cent on timescales of tens of seconds, slowly decreasing toward higher frequencies. Finally, we report on the first detection of a linear rms-flux relation in the emission from a low-mass X-ray binary jet, on timescales where little correlation is found between the X-rays and the jet emission itself. This suggests that either the inflow variations and jet IR emission are coupled by a non-linear or time-variable transform, or that the IR rms-flux relation is not transferred from the inflow to the jet, but is an intrinsic property of emission processes in the jet., 11 pages, 8 figures. Accepted for publication in MNRAS

We use Rossi X-ray Timing Explorer observations to conduct a population study of the timing properties of accretion-powered neutron star Low Mass X-ray Binaries (LMXBs), following a model-independent approach that was originally applied to black hole systems. The ratios of integrated power in four equally spaced Fourier frequency bands allow power spectral shapes to be parameterised with two `power colour' values, providing a simple way of tracking the evolution in timing properties across observations. We show that neutron star LMXBs follow a remarkably similar power spectral evolution to black hole LMXBs, confirming that the broadband noise variability seen in both types of system has a common origin in the accretion flow. Both neutron stars and black holes follow a similar clear track in the power colour-colour diagram as they transition from the hard through soft states. Quantifying the position on this oval track using a single parameter, the power-spectral `hue', we find that the transition in X-ray spectral shape occurs at the same hue for both neutron star and black hole systems. The power colours of Z sources map on to those of soft state black holes, while those of atoll sources correspond to all spectral states. There is no clear dependence of power colour on neutron star spin, or on whether the neutron star is clearly magnetised (determined by ms X-ray pulsations)., Comment: 21 pages, 19 figures, accepted for publication in MNRAS

We use a simple one-dimensional parameterisation of timing properties to show that hard and hard-intermediate state transient black hole X-ray binaries with the same power-spectral shape have systematically harder X-ray power-law emission in higher-inclination systems. We also show that the power-spectral shape and amplitude of the broadband noise (with low-frequency quasi-periodic oscillations, QPOs, removed) is independent of inclination, confirming that it is well-correlated with the intrinsic structure of the emitting regions and that the "type C" QPO, which is inclination-dependent, has a different origin to the noise, probably geometric. Our findings suggest that the power-law emission originates in a corona which is flattened in the plane of the disc, and not in a jet-like structure which would lead to softer spectra at higher inclinations. However, there is tentative evidence that the inclination-dependence of spectral shape breaks down deeper into the hard state. This suggests either a change in the coronal geometry and possible evidence for contribution from jet emission, or alternatively an even more optically thin flow in these states., Comment: 6 Pages, 4 Figures, accepted as a Letter by MNRAS

The light curves of a diverse range of accreting objects show characteristic linear relationships between the short-term rms amplitude of variability and the flux as measured on longer time-scales. This behaviour is thought to be imprinted on the light curves by accretion rate fluctuations on different time-scales, propagating and coupling together through the accretion flow. Recently, a simple mathematical interpretation has been proposed for the rms-flux relation, where short-term variations are modulated by a single slower process. Here we show that this model was already considered and ruled out by another publication on the grounds that it did not produce the observed broad time-scale dependence of the rms-flux relation and associated lognormal flux distribution. We demonstrate the problems with the model via mathematical arguments and a case-study of Cyg X-1 data compared with numerical simulations. We also highlight another conclusion of our original work, which is that a linear rms-flux relation is easy to produce using a variety of models with positively skewed flux distributions. Observing such a relation in a non-accreting object (e.g. in solar flares) does not necessarily imply a phenomenological connection with the behaviour of accretion flows, unless the relation is seen over a similarly broad range of time-scales.

XMM-Newton is capable of making a transformational advance in our understanding of how luminous accreting black holes work, by dedicating about 10 per cent of future observing time to long observations, of order Megaseconds, to X-ray variable Active Galactic Nuclei (AGN) research. This would enable reverberation studies, already a commonplace feature of AGN, to proceed to the next level and follow the behaviour of the powerful dynamic corona. Such a dedicated legacy programme can only be carried out with XMM-Newton., Comment: 5 pages, 5 figures for publication in Astronomische Nachrichten as part of the proceedings of the workshop "XMM-Newton: The Next Decade"

We report the results of intensive X-ray, UV and optical monitoring of the Seyfert 1 galaxy NGC 4593 with Swift. There is no intrinsic flux-related spectral change in the the variable components in any band with small apparent variations due only to contamination by a second constant component, possibly a (hard) reflection component in the X-rays and the (red) host galaxy in the UV/optical bands. Relative to the shortest wavelength band, UVW2, the lags of the other UV and optical bands are mostly in agreement with the predictions of reprocessing of high energy emission from an accretion disc. The U-band lag is, however, far larger than expected, almost certainly because of reprocessed Balmer continuum emission from the more distant broad line region gas. The UVW2 band is well correlated with the X-rays but lags by ~6x more than expected if the UVW2 results from reprocessing of X-rays on the accretion disc. However, if the lightcurves are filtered to remove variations on timescales >5d, the lag approaches the expectation from disc reprocessing. MEMEcho analysis shows that direct X-rays can be the driver of most of the variations in the UV/optical bands as long as the response functions for those bands all have long tails (up to 10d) in addition to a strong peak (from disc reprocessing) at short lag (, 18 pages, 17 figures. Submitted to Monthly Notices of the Royal Astronomical Society

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548., Comment: 22 pages, 13 figures, accepted to ApJ

High-frequency iron K reverberation lags, where the red wing of the line responds before the line centroid, are a robust signature of relativistic reflection off the inner accretion disc. In this letter, we report the discovery of the Fe K lag in PG 1244+026 from ~120 ks of data (1 orbit of the XMM-Newton telescope). The amplitude of the lag with respect to the continuum is 1000 s at a frequency of ~1e-4 Hz. We also find a possible frequency-dependence of the line: as we probe higher frequencies (i.e. shorter timescales from a smaller emitting region) the Fe K lag peaks at the red wing of the line, while at lower frequencies (from a larger emitting region) we see the dominant reflection lag from the rest frame line centroid. The mean energy spectrum shows a strong soft excess, though interestingly, there is no indication of a soft lag. Given that this source has radio emission and it has little reported correlated variability between the soft excess and the hard band, we explore one possible explanation in which the soft excess in this source is dominated by the steep power-law like emission from a jet, and that a corona (or base of the jet) irradiates the inner accretion disc, creating the blurred reflection features evident in the spectrum and the lag. General Relativistic ray-tracing models fit the Fe K lag well, with the best-fit giving a compact X-ray source at a height of 5 gravitational radii and a black hole mass of 1.3e7 Msun., Comment: 6 pages, 6 figures, resubmitted to MNRAS after moderate revisions. This paper focuses on the discovery of the Fe K reverberation lag in PG 1244+026. We point the interested reader to Alston, Done & Vaughan (See today: arXiv:submit/0851673), which focuses on the soft lags in this source

Radio galaxies are among the largest and most powerful single objects known and are found at variety of redshifts, hence they are believed to have had a significant impact on the evolving Universe. Their relativistic jets inject considerable amounts of energy into the environments in which the sources reside; thus the knowledge of the fundamental properties (such as kinetic luminosities, lifetimes and ambient gas densities) of these sources is crucial for understanding AGN feedback in galaxy clusters. In this work, we explore the intrinsic and extrinsic fundamental properties of Fanaroff-Riley II (FR II) objects through the construction of multidimensional Monte Carlo simulations which use complete, flux limited radio catalogues and semi-analytical models of FR IIs' time evolution to create artificial samples of radio galaxies. This method allows us to set better limits on the confidence intervals of the intrinsic and extrinsic fundamental parameters and to investigate the total energy produced and injected to the clusters' environments by populations of FR IIs at various cosmological epochs (0.0

We use archival XMM-Newton observations of Ark 564 and Mrk 335 to calculate the frequency dependent time-lags for these two well-studied sources. We discover high-frequency Fe K lags in both sources, indicating that the red wing of the line precedes the rest frame energy by roughly 100 s and 150 s for Ark 564 and Mrk 335, respectively. Including these two new sources, Fe K reverberation lags have been observed in seven Seyfert galaxies. We examine the low-frequency lag-energy spectrum, which is smooth, and shows no feature of reverberation, as would be expected if the low-frequency lags were produced by distant reflection off circumnuclear material. The clear differences in the low and high frequency lag-energy spectra indicate that the lags are produced by two distinct physical processes. Finally, we find that the amplitude of the Fe K lag scales with black hole mass for these seven sources, consistent with a relativistic reflection model where the lag is the light travel delay associated with reflection of continuum photons off the inner disc., 10 pages, 12 figures, accepted for publication in MNRAS

We present an analysis of the X-ray time lags for the highly variable Seyfert 1 galaxy NGC 4051, based on a series of XMM-Newton observations taken in 2009. We investigate the Fourier frequency dependent time lags in the light curves between the 0.3--1.0 keV and 2.0--5.0 keV energy bands as a function of source flux, including simultaneous modelling of the resulting lag-frequency spectra. We find the shape of the lag-frequency spectra to vary significantly and systematically with source flux. We model the lag-frequency spectra using simple transfer functions, and find that two time lag components are required, one in each energy band. The simplest acceptable fits have only the relative contribution of the lagged component in the hard band varying with flux level, which can be associated with changes in the energy spectrum. We discuss the interpretation of these results in terms of the currently popular models for X-ray time lags., Comment: 9 pages, 4 figures. Accepted for publication in MNRAS

Reverberation lags in AGN were first discovered in the NLS1 galaxy, 1H0707-495. We present a follow-up analysis using 1.3 Ms of data, which allows for the closest ever look at the reverberation signature of this remarkable source. We confirm previous findings of a hard lag of ~100 seconds at frequencies v ~ [0.5 - 4] e-4 Hz, and a soft lag of ~30 seconds at higher frequencies, v ~ [0.6 - 3] e-3 Hz. These two frequency domains clearly show different energy dependences in their lag spectra. We also find evidence for a signature from the broad Fe K line in the high frequency lag spectrum. We use Monte Carlo simulations to show how the lag and coherence measurements respond to the addition of Poisson noise and to dilution by other components. With our better understanding of these effects on the lag, we show that the lag-energy spectra can be modelled with a scenario in which low frequency hard lags are produced by a compact corona responding to accretion rate fluctuations propagating through an optically thick accretion disc, and the high frequency soft lags are produced by short light-travel delay associated with reflection of coronal power-law photons off the disc., 11 pages, 10 figures. Accepted for publication in MNRAS

We report the first study of Fourier-frequency-dependent coherence and phase/time lags at optical wavelengths of cataclysmic variables (MV Lyr and LU Cam) displaying typical flickering variability in white light. Observations were performed on the William Herschel Telescope using ULTRACAM. Lightcurves for both systems have been obtained with the SDSS filters $u'$, $g'$ and $r'$ simultaneously with cadences between $\approx0.5-2$ seconds, and allow us to probe temporal frequencies between ~10^{-3} Hz and ~1 Hz. We find high levels of coherence between the u', g' and r' lightcurves up to at least ~10^{-2} Hz. Furthermore we detect red/negative lags where the redder bands lag the bluer ones at the lowest observed frequencies. For MV Lyr time lags up to ~3 seconds are observed, whilst LU Cam displays larger time lags of ~10 seconds. Mechanisms which seek to explain red/negative lags observed in X-ray binaries and Active Galactic Nuclei involve reflection of photons generated close to the compact object onto the surface layers of the accretion disk, where the lag delay is simply the light-travel time from the emitting source to the reflecting accretion disk area. Although this could be a viable explanation for the lags observed in MV Lyr, the lags observed in LU Cam are too large to be explained by reflection from the disk and/or the donor star. We suggest reprocessing on the thermal timescale of boundary layer photons onto the accretion disk as a possible mechanism to explain the lags observed in accreting white dwarfs, or reverse (inside-out) shocks within the disk travelling through cooler disk regions as they move outwards., 8 pages, 3 figures, 1 table. Accepted for publication in MNRAS

We investigated, using spectral-timing analysis, the characterization of highly ionized outflows in Seyfert galaxies, the so-called warm absorbers. Here, we present our results on the extensive ~ 600 ks of XMM-Newton archival observations of the bright and highly variable Seyfert 1 galaxy NGC 4051, whose spectrum has revealed a complex multicomponent wind. Making use of both RGS and EPIC-pn data, we performed a detailed analysis through a time-dependent photoionization code in combination with spectral and Fourier spectral-timing techniques. The source light curves and the warm absorber parameters obtained from the data were used to simulate the response of the gas due to variations in the ionizing flux of the central source. The resulting time variable spectra were employed to predict the effects of the warm absorber on the time lags and coherence of the energy dependent light curves. We have found that, in the absence of any other lag mechanisms, a warm absorber with the characteristics of the one observed in NGC 4051, is able to produce soft lags, up to 100 s, on timescales of ~ hours. The time delay is associated with the response of the gas to changes in the ionizing source, either by photoionization or radiative recombination, which is dependent on its density. The range of radial distances that, under our assumptions, yield longer time delays are comparable to the existing estimates of the location of the warm absorber in NGC 4051. For this reason, we suggest that it is likely that the observed X-ray time lags may carry a signature of the warm absorber response time, to changes in the ionizing continuum. These results highlight the importance of understanding the contribution of the warm absorber to the AGN X-ray time lags, since it is also vital information for interpreting the lags associated with propagation and reverberation effects in the inner emitting regions., Comment: 13 pages, 12 figures, Accepted for publication by A&A

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., 16 pages, 10 figures; accepted for publication in MNRAS

X-ray reverberation, where light-travel time delays map out the compact geometry around the inner accretion flow in supermassive black holes, has been discovered in several of the brightest, most variable and well-known Seyfert galaxies. In this work, we expand the study of X-ray reverberation to all Seyfert galaxies in the $\textit{XMM–Newton}$ archive above a nominal rms variability and exposure level (a total of 43 sources). Approximately 50 per cent of sources exhibit iron K reverberation, in that the broad iron K emission line responds to rapid variability in the continuum. We also find that on long time-scales, the hard band emission lags behind the soft band emission in 85 per cent of sources. This ‘low-frequency hard lag’ is likely associated with the coronal emission, and so this result suggests that most sources with X-ray variability show intrinsic variability from the nuclear region. We update the known iron K lag amplitude versus black hole mass relation, and find evidence that the height or extent of the coronal source (as inferred by the reverberation time delay) increases with mass accretion rate., EK thanks Ari Laor for interesting discussions on this work and acknowledges support from the International Space Science Institute. This work is based on observations obtained with XMM–Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. EK thanks the Gates Cambridge Scholarship and the Hubble Fellowship Program. Support for Program number HSTHF2-51360.001-A was provided by NASA through a Hubble Fellowship grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. EK, WNA, and ACF acknowledge support from the European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 312789, StrongGravity. CSR acknowledges support from NASA under grant NNX14AF86G. EMC gratefully acknowledges support from the National Science Foundation through CAREER award number 1351222.