We present results of temporal and spectral analyses on four XMM-Newton EPIC pn observations of IRAS 17020+4544, a narrow-line Seyfert 1 galaxy with evidence of a radio jet. Analysis of the light curves reveals that this radio-loud source does not behave like the bulk population of its radio-quiet counterparts. A trend of spectral hardening with increased flux is found. Variability is found to increase with energy, though it decreases as the spectrum hardens. The first 40 ks of the most recent observation behave uniquely among the epochs, exhibiting a softer spectral state than at any other time. Possible non-stationarity at low energies is found, with no such effect present at higher energies, suggesting at least two distinct spectral components. A reverberation signature is confirmed, with the lag-frequency, lag-energy, and covariance spectra changing significantly during the soft-state epoch. The temporal analysis suggests a variable power-law in the presence of a reflection component, thus motivating such a fit for the 0.3-10 keV EPIC pn spectra from all epochs. We find an acceptable spectral fit using the timing-motivated parameters and report the detection of a broad Fe K emission line, requiring an additional model component beyond the reflection spectrum. We discuss links between this source and other narrow-line Seyfert 1 sources that show evidence of jet activity, finding similarities among this currently very limited sample of interesting objects.

We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X-1. We model the multi-epoch data with two-component accretion disc models designed to approximate super-Eddington accretion, allowing for both a black hole and a neutron star accretor. With regards to the hotter disc component, the data trace out two distinct tracks in the luminosity-temperature plane, with larger emitting radii and lower temperatures seen at higher observed fluxes. Despite this apparent anti-correlation, each of these tracks individually shows a positive luminosity-temperature relation. Both are broadly consistent with $L\propto{T}^{4}$, as expected for blackbody emission with a constant area, and also with $L\propto{T}^{2}$, as may be expected for an advection-dominated disc around a black hole. We consider a variety of possibilities for this unusual behaviour. Scenarios in which the innermost flow is suddenly blocked from view by outer regions of the super-Eddington disc/wind can explain the luminosity-temperature behaviour, but are difficult to reconcile with the lack of strong variability at higher energies, assuming this emission arises from the most compact regions. Instead, we may be seeing evidence for further radial stratification of the accretion flow than is included in the simple models considered, with a combination of winds and advection resulting in the suppressed high-energy variability.

We explore a disc origin for the highly-blueshifted, variable absorption lines seen in the X-ray spectrum of the Narrow Line Seyfert 1 galaxy IRAS13224-3809. The blueshift corresponds to a velocity of about 0.25c. Such features in other Active Galactic Nuclei are often interpreted as UltraFast Outflows (UFOs). The velocity is of course present in the orbital motions of the inner disk. The absorption lines in IRAS13224-3809 are best seen when the flux is low and the reflection component of the disk is strong relative to the power-law continuum. The spectra are consistent with a model in which the reflection component passes through a thin, highly-ionized absorbing layer at the surface of the inner disc, the blue-shifted side of which dominates the flux due to relativistic aberration (the disc inclination is about 70 deg). No fast outflow need occur beyond the disc.

We present a qualitative search for ultra-fast outflows (UFOs) in excess variance spectra of radio-quiet active galactic nuclei (AGN). We analyse 42 sources from the Tombesi et al. (2010) spectroscopic UFO detection sample, and an additional 22 different sources from the Kara et al. (2016) variability sample. A total of 58 sources have sufficient observational data from XMM-Newton EPIC-pn and variability for an excess variance spectrum to be calculated. We examine these spectra for peaks corresponding to variable blue-shifted H- and He-like ion absorption lines from UFOs. We find good evidence for such outflows in 28% of the AGN sample and weak evidence in a further 31%, meaning that $\sim$ 30-60% of the AGN sample hosts such UFOs. The mean and median blue-shifted velocity is found to be $\sim$ 0.14c and 0.12c, respectively. Current variability methods allow for a fast, model-independent determination of UFOs, however, further work needs to be undertaken to better characterize the statistical significance of the peaks in these spectra by more rigorous modelling. Detecting good evidence for variable UFO lines in a large number of sources also lays the groundwork for detailed analysis of the variability timescales of the absorbers. This will allow us to probe their densities and hence distances from the central super-massive black hole.

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.

Most ultraluminous X-ray sources (ULXs) are thought to be powered by neutron stars and black holes accreting beyond the Eddington limit. If the compact object is a black hole or a neutron star with a magnetic field $\lesssim10^{12}$ G, the accretion disc is expected to thicken and launch powerful winds driven by radiation pressure. Evidence of such winds has been found in ULXs through the high-resolution spectrometers onboard XMM-Newton, but several unknowns remain, such as the geometry and launching mechanism of these winds. In order to better understand ULX winds and their link to the accretion regime, we have undertaken a major campaign with XMM-Newton to study the ULX NGC 1313 X-1, which is known to exhibit strong emission and absorption features from a mildly-relativistic wind. The new observations show clear changes in the wind with a significantly weakened fast component (0.2c) and the rise of a new wind phase which is cooler and slower (0.06-0.08c). We also detect for the first time variability in the emission lines which indicates an origin within the accretion disc or in the wind. We describe the variability of the wind in the framework of variable super-Eddington accretion rate and discuss a possible geometry for the accretion disc.

We present simple XSPEC models for fitting excess variance spectra of AGN. Using a simple Monte-Carlo approach, we simulate a range of spectra corresponding to physical parameters varying, then calculate the resulting variance spectra. Starting from a variable power-law, we build up a set of models corresponding to the different physical processes that can affect the final excess variance spectrum. We show that the complex excess variance spectrum of IRAS 13224-3809 can be well described by such an intrinsic variability model, where the power-law variability is damped by relativistic reflection and enhanced by an ultra fast outflow. The reflection flux is correlated with that of the power-law, but not perfectly. We argue that this correlation is stronger at high frequencies, where reverberation lags are detected, while excess variance spectra are typically dominated by low frequency variability.

Ultraluminous X-ray Sources (ULXs) provide a unique opportunities to probe the geometry and energetics of super-Eddington accretion. The radiative processes involved in super-Eddington accretion are not well understood, and so studying correlated variability between different energy bands can provide insights into the causal connection between different emitting regions. We present a spectral-timing analysis of NGC 1313 X-1 from a recent XMM-Newton campaign. The spectra can be decomposed into two thermal-like components, the hotter of which may originate from the inner accretion disc, and the cooler from an optically thick outflow. We find correlated variability between hard (2-10 keV) and soft (0.3-2 keV) bands on kilosecond timescales, and find a soft lag of ~150 seconds. The covariance spectrum suggests that emission contributing to the lags is largely associated with the hotter of the two thermal-like components, likely originating from the inner accretion flow. This is only the third ULX to exhibit soft lags. The lags range over three orders of magnitude in amplitude, but all three are ~5 to ~20 percent of the corresponding characteristic variability timescales. If these soft lags can be understood in the context of a unified picture of ULXs, then lag timescales may provide constraints on the density and extent of radiatively-driven outflows.

We report the detection of weak pulsations from the archetypal ultraluminous X-ray source (ULX) NGC 1313 X-2. Acceleration searches reveal sinusoidal pulsations in segments of two out of six new deep observations of this object, with a period of $\sim$ 1.5 s and a pulsed fraction of $\sim$ 5%. We use Monte Carlo simulations to demonstrate that the individual detections are unlikely to originate in false Poisson noise detections given their very close frequencies; their strong similarity to other pulsations detected from ULXs also argues they are real. The presence of a large bubble nebula surrounding NGC 1313 X-2 implies an age of order 1 Myr for the accreting phase of the ULX, which implies that the neutron star's magnetic field has not been suppressed over time by accreted material, nor has the neutron star collapsed into a black hole, despite an average energy output into the nebula two orders of magnitude above Eddington. This argues that most of the accreted material has been expelled over the lifetime of the ULX, favouring physical models including strong winds and/or jets for neutron star ULXs.

SS433 is the only Galactic binary system known to accrete at highly super-critical rates, analogous to tidal disruption events, and needed to explain the mass of some high redshift quasars. Probing the inner regions of SS433 in the X-rays is crucial to understanding this system, and super-critical accretion in general, but has not yet been possible due to obscuration. NuSTAR observed SS433 in the hard X-ray band across multiple phases of its super-orbital precession period. Spectral-timing tools have allowed us to confirm that the hard X-ray emission from the inner regions is scattered towards us by the walls of the wind-cone. By comparing to numerical models, we determine an intrinsic X-ray luminosity of >= 3x10^37 erg/s and that, if viewed face on, the apparent luminosity would be > 1x10^39 erg/s, confirming its long-suspected nature as an ultraluminous X-ray source (ULX). A lag due to absorption by Fe XXV/XXVI in outflowing material travelling at least 0.14-0.29c matches absorption lines seen in ULXs and - in future - will allow us to map a super-critical outflow for the first time.

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.

Accreting compact objects show variations in source flux over a broad range of timescales and in all wavebands. The light curves typically show a lognormal distribution of flux and a linear relation between flux and rms. It has been demonstrated that an exponential transform of an underlying (and unobserved) Gaussian stochastic process provides a very good description of the light curves with these observed properties. Recently, a non-stationary power spectrum was observed on fast timescales (~ days) in the active galaxy, IRAS 13224--3809, as well as a non-lognormal flux distribution and non-linear rms-flux relation. Here, we investigate the affects of piecewise non-stationary power spectra on the resultant flux distribution and rms-flux relation. We demonstrate that the simple "exponentiation" model successfully reproduces the observed quantities, even when the light curves are non-stationary. We also demonstrate how non-lognormal flux distributions and rms-flux relations inconsistent with a linear model can be erroneously produced from poorly sampled PSDs. This is of particular importance for AGN surveys where very long baselines are required to sample the PSD down to low enough frequencies.

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.

The narrow line Seyfert 1 (NLS1) IRAS 13224-3809 is the most X-ray variable active galactic nucleus (AGN), exhibiting 0.3-10 keV flux changes of over an order of magnitude within an hour. We report on the results of the 1.5 Ms 2016 XMM-Newton/NuSTAR observing campaign, which revealed the presence of a 0.24c ultra-fast outflow in addition to the well-known strong relativistic reflection. We also summarise other key results of the campaign, such as the first detection of a non-linear RMS-flux relation in an accreting source, correlations between outflow absorption strength/velocity and source flux, and a disconnect between the X-ray and UV emission. Our results are consistent with a scenario where a disk wind is launched close to the black hole, imprinting absorption features into the spectrum and variability.

We present the first results from a detailed spectral-timing analysis of a long ($\sim$130 ks) XMM-Newton observation and quasi-simultaneous NuSTAR and Swift observations of the highly-accreting narrow-line Seyfert 1 galaxy Mrk 1044. The broadband (0.3$-$50 keV) spectrum reveals the presence of a strong soft X-ray excess emission below $\sim$1.5 keV, iron K$_{\alpha}$ emission complex at $\sim$6$-$7 keV and a `Compton hump' at $\sim$15$-$30 keV. We find that the relativistic reflection from a high-density accretion disc with a broken power-law emissivity profile can simultaneously explain the soft X-ray excess, highly ionized broad iron line and the Compton hump. At low frequencies ($[2-6]\times10^{-5}$ Hz), the power-law continuum dominated 1.5$-$5 keV band lags behind the reflection dominated 0.3$-$1 keV band, which is explained with a combination of propagation fluctuation and Comptonization processes, while at higher frequencies ($[1-2]\times10^{-4}$ Hz), we detect a soft lag which is interpreted as a signature of X-ray reverberation from the accretion disc. The fractional root-mean-squared (rms) variability of the source decreases with energy and is well described by two variable components: a less variable relativistic disc reflection and a more variable direct coronal emission. Our combined spectral-timing analyses suggest that the observed broadband X-ray variability of Mrk~1044 is mainly driven by variations in the location or geometry of the optically thin, hot corona.

We present a detailed spectral analysis of the recent 1.5\,Ms XMM-Newton observing campaign on the narrow line Seyfert 1 galaxy IRAS~13224$-$3809, taken simultaneously with 500\,ks of NuSTAR data. The X-ray lightcurve shows three flux peaks, registering at about 100 times the minimum flux seen during the campaign, and rapid variability with a time scale of kiloseconds. The spectra are well fit with a primary powerlaw continuum, two relativistic-blurred reflection components from the inner accretion disk with very high iron abundance, and a simple blackbody-shaped model for the remaining soft excess. The spectral variability is dominated by the power law continuum from a corona region within a few gravitational radii from the black hole. Additionally, blueshifted Ne \textsc{x}, Mg \textsc{xii}, Si \textsc{xiv} and S \textsc{xvi} absorption lines are identified in the stacked low-flux spectrum, confirming the presence of a highly ionized outflow with velocity up to $v= 0.267$ and $0.225$\,c. We fit the absorption features with \texttt{xstar} models and find a relatively constant velocity outflow through the whole observation. Finally, we replace the \texttt{bbody} and supersolar abundance reflection models by fitting the soft excess successfully with the extended reflection model \texttt{relxillD}, which allows for higher densities than the standard \texttt{relxill} model. This returns a disk electron density $n_{\rm e}>10^{18.7}$\,cm$^{-3}$ and lowers the iron abundance from $Z_{\rm Fe}=24^{+3}_{-4}Z_\odot$ with $n_{\rm e}\equiv10^{15}$\,cm$^{-3}$ to $Z_{\rm Fe}=6.6^{+0.8}_{-2.1}Z_\odot$.

Ultrafast outflows (UFOs) are the most extreme winds launched by active galactic nuclei (AGN) due to their mildly-relativistic speeds (~0.1-0.3c) and are thought to significantly contribute to galactic evolution via AGN feedback. Their nature and launching mechanism are however not well understood. Recently, we have discovered the presence of a variable UFO in the narrow-line Seyfert 1 IRAS 13224-3809. The UFO varies in response to the brightness of the source. In this work we perform flux-resolved X-ray spectroscopy to study the variability of the UFO and found that the ionisation parameter is correlated with the luminosity. In the brightest states the gas is almost completely ionised by the powerful radiation field and the UFO is hardly detected. This agrees with our recent results obtained with principal component analysis. We might have found the tip of the iceberg: the high ionisation of the outflowing gas may explain why it is commonly difficult to detect UFOs in AGN and possibly suggest that we may underestimate their actual feedback. We have also found a tentative correlation between the outflow velocity and the luminosity, which is expected from theoretical predictions of radiation-pressure driven winds. This trend is rather marginal due to the Fe XXV-XXVI degeneracy. Further work is needed to break such degeneracy through time-resolved spectroscopy.

We present results from the optical, ultraviolet and X-ray monitoring of the NLS1 galaxy IRAS 13224-3809 taken with Swift and XMM-Newton during 2016. IRAS 13224-3809 is the most variable bright AGN in the X-ray sky and shows strong X-ray reflection, implying that the X-rays strongly illuminate the inner disc. Therefore, it is a good candidate to study the relationship between coronal X-ray and disc UV emission. However, we find no correlation between the X-ray and UV flux over the available ~40 day monitoring, despite the presence of strong X-ray variability and the variable part of the UV spectrum being consistent with irradiation of a standard thin disc. This means either that the X-ray flux which irradiates the UV emitting outer disc does not correlate with the X-ray flux in our line of sight and/or that another process drives the majority of the UV variability. The former case may be due to changes in coronal geometry, absorption or scattering between the corona and the disc.

We present an analysis of the long-term X-ray variability of the extreme narrow-line Seyfert 1 (NLS1) galaxy IRAS 13224-3809 using principal component analysis (PCA) and fractional excess variability (Fvar) spectra to identify model-independent spectral components. We identify a series of variability peaks in both the first PCA component and Fvar spectrum which correspond to the strongest predicted absorption lines from the ultra-fast outflow (UFO) discovered by Parker et al. (2017). We also find higher order PCA components, which correspond to variability of the soft excess and reflection features. The subtle differences between RMS and PCA results argue that the observed flux-dependence of the absorption is due to increased ionization of the gas, rather than changes in column density or covering fraction. This result demonstrates that we can detect outflows from variability alone, and that variability studies of UFOs are an extremely promising avenue for future research.

In recent work with high-resolution grating spectrometers (RGS) aboard XMM-Newton Pinto et al. (2016) have discovered that two bright and archetypal ultraluminous X-ray sources (ULXs) have strong relativistic winds in agreement with theoretical predictions of high accretion rates. It has been proposed that such winds can become optically thick enough to block and reprocess the disk X-ray photons almost entirely, making the source appear as a soft thermal emitter or ultraluminous supersoft X-ray source (ULS). To test this hypothesis we have studied a ULX where the wind is strong enough to cause significant absorption of the hard X-ray continuum: NGC 55 ULX. The RGS spectrum of NGC 55 ULX shows a wealth of emission and absorption lines blueshifted by significant fractions of the light speed (0.01 - 0.20)c indicating the presence of a powerful wind. The wind has a complex dynamical structure with the ionization state increasing with the outflow velocity, which may indicate launching from different regions of the accretion disk. The comparison with other ULXs such as NGC 1313 X-1 and NGC 5408 X-1 suggests that NGC 55 ULX is being observed at higher inclination. The wind partly absorbs the source flux above 1 keV, generating a spectral drop similar to that observed in ULSs. The softening of the spectrum at lower (~ Eddington) luminosities and the detection of a soft lag agree with the scenario of wind clumps crossing the line of sight, partly absorbing and reprocessing the hard X-rays from the innermost region.

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.

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.

We analyse a sample of 21 active galactic nuclei (AGN) using data from the Swift satellite to study the variability properties of the population in the X-ray, UV and optical band. We find that the variable part of the UV-optical emission has a spectrum consistent with a powerlaw, with an average index of $-2.21\pm0.13$, as would be expected from central illumination of a thin disc (index of -7/3). We also calculate the slope of a powerlaw from UV to X-ray variable emission, $\alpha_{\rm OX,Var}$; the average for this sample is $\alpha_{\rm OX,Var} = -1.06 \pm 0.04$. The anticorrelation of $\alpha_{\rm OX}$ with the UV luminosity, $L_{\rm UV}$, previously found in the average emission is also present in the variable part: $\alpha_{\rm OX,Var} = (-0.177 \pm 0.083) {\rm log} (L_{\rm \nu,Var} (2500\,\AA)) + (3.88 \pm 2.33)$. Correlated variability between the emission in X-rays and UV is detected significantly for 9 of the 21 sources. All these cases are consistent with the UV lagging the X-rays, as would be seen if the correlated UV variations were produced by the reprocessing of X-ray emission. The observed UV lags are tentatively longer than expected for a standard thin disc.

There have recently been several reports of apparently periodic variations in the light curves of quasars, e.g. PG 1302-102 by Graham et al. (2015a). Any quasar showing periodic oscillations in brightness would be a strong candidate to be a close binary supermassive black hole and, in turn, a candidate for gravitational wave studies. However, normal quasars -- powered by accretion onto a single, supermassive black hole -- usually show stochastic variability over a wide range of timescales. It is therefore important to carefully assess the methods for identifying periodic candidates from among a population dominated by stochastic variability. Using a Bayesian analysis of the light curve of PG 1302-102, we find that a simple stochastic process is preferred over a sinusoidal variations. We then discuss some of the problems one encounters when searching for rare, strictly periodic signals among a large number of irregularly sampled, stochastic time series, and use simulations of quasar light curves to illustrate these points. From a few thousand simulations of steep spectrum (`red noise') stochastic processes, we find many simulations that display few-cycle periodicity like that seen in PG 1302-102. We emphasise the importance of calibrating the false positive rate when the number of targets in a search is very large.

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 XMM-Newton archive above a nominal rms variability and exposure level (a total of 43 sources). 50 per cent of source 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 timescales, 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 vs. 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.

We present a multi-epoch X-ray spectral analysis of the Seyfert 1 galaxy Fairall 9. Our analysis shows that Fairall 9 displays unique spectral variability in that its ratio residuals to a simple absorbed power law in the 0.5-10 keV band remain constant with time in spite of large variations in flux. This behavior implies an unchanging source geometry and the same emission processes continuously at work at the timescale probed. With the constraints from NuSTAR on the broad-band spectral shape, it is clear that the soft excess in this source is a superposition of two different processes, one being blurred ionized reflection in the innermost parts of the accretion disk, and the other a continuum component such as spatially distinct Comptonizing region. Alternatively, a more complex primary Comptonization component together with blurred ionized reflection could be responsible.

Quasi-periodic oscillations (QPOs) are coherent peaks of variability power observed in the X-ray power spectra (PSDs) of stellar mass X-ray binaries (XRBs). A scale invariance of the accretion process implies they should be present in the active galactic nuclei. The first robust detection was a $\sim 1$ hr periodicity in the Seyfert galaxy RE J1034+396 from a $\sim 90$ ks XMM-Newton observation, however, subsequent observations failed to detect the QPO in the $0.3-10.0$ keV band. In this talk we present the recent detection of the $\sim 1$ hr periodicity in the $1.0-4.0$ keV band of 4 further low-flux/spectrally-harder observations of RE J1034+396 (see Alston et al 2014). We also present recent work on the discovery of a QPO in the Seyfert galaxy, MS 22549-3712, which again is only detected in energy bands associated with the primary power-law continuum emission (Alston et al 2015). We conclude these features are most likely analogous to the high-frequency QPOs observed in XRBs. In both sources, we also see evidence for X-ray reverberation at the QPO frequency, where soft X-ray bands and Iron K$\alpha$ emission lag the primary X-ray continuum. These time delays may provide another diagnostic for understanding the underlying QPO mechanism observed in accreting black holes.

We present simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR ) and Suzaku observations of the X-ray binary Cygnus X-1 in the hard state. This is the first time this state has been observed in Cyg X-1 with NuSTAR, which enables us to study the reflection and broad-band spectra in unprecedented detail. We confirm that the iron line cannot be fit with a combination of narrow lines and absorption features, and instead requires a relativistically blurred profile in combination with a narrow line and absorption from the companion wind. We use the reflection models of Garcia et al. (2014) to simultaneously measure the black hole spin, disk inner radius, and coronal height in a self-consistent manner. Detailed fits to the iron line profile indicate a high level of relativistic blurring, indicative of reflection from the inner accretion disk. We find a high spin, a small inner disk radius, and a low source height, and rule out truncation to greater than three gravitational radii at the 3{\sigma} confidence level. In addition, we find that the line profile has not changed greatly in the switch from soft to hard states, and that the differences are consistent with changes in the underlying reflection spectrum rather than the relativistic blurring. We find that the blurring parameters are consistent when fitting either just the iron line or the entire broad-band spectrum, which is well modelled with a Comptonized continuum plus reflection model.

We report the discovery of a $\sim 1.5 \times 10^{-4}$ Hz ($\sim 2$ hr) X-ray quasi-periodic oscillation (QPO) in the active galaxy MS 2254.9-3712, using a $\sim 70$ ks XMM-Newton observation. The QPO is significantly detected ($\sim 3.3 \sigma$) in the $1.2 - 5.0$ keV band only, connecting its origin with the primary X-ray power-law continuum. We detect a highly coherent soft lag between the $0.3 - 0.7$ keV and $1.2 - 5.0$ keV energy bands at the QPO frequency and at a frequency band in a 3:2 ratio, strongly suggesting the presence of a QPO harmonic. An iron K$\alpha$ reverberation lag is found at the harmonic frequency, indicating the reflecting material subtends some angle to the primary continuum, which is modulated by the QPO mechanism. Frequency resolved spectroscopy reveals the QPO and harmonic to have a hard energy dependence. These properties of the QPO variability, together with the current black hole mass estimate, $M_{\rm bh} \sim 4 \times 10^{6} M_{\rm sun}$, are consistent with the QPO originating from the same process as the high frequency QPO phenomenon observed in black hole X-ray binaries. Principle component analysis reveals the spectral variability in MS 2254.9-3712 is similar to that of the active galaxy RE J1034+396, a source which also displays an X-ray QPO. This suggests a distinct spectral variability pattern for accreting black holes when in a state where QPOs are present.

We analyse a sample of 26 active galactic nuclei with deep XMM-Newton observations, using principal component analysis (PCA) to find model independent spectra of the different variable components. In total, we identify at least 12 qualitatively different patterns of spectral variability, involving several different mechanisms, including five sources which show evidence of variable relativistic reflection (MCG-6-30-15, NGC 4051, 1H 0707-495, NGC 3516 and Mrk 766) and three which show evidence of varying partial covering neutral absorption (NGC 4395, NGC 1365, and NGC 4151). In over half of the sources studied, the variability is dominated by changes in a power law continuum, both in terms of changes in flux and power law index, which could be produced by propagating fluctuations within the corona. Simulations are used to find unique predictions for different physical models, and we then attempt to qualitatively match the results from the simulations to the behaviour observed in the real data. We are able to explain a large proportion of the variability in these sources using simple models of spectral variability, but more complex models may be needed for the remainder. We have begun the process of building up a library of different principal components, so that spectral variability in AGN can quickly be matched to physical processes. We show that PCA can be an extremely powerful tool for distinguishing different patterns of variability in AGN, and that it can be used effectively on the large amounts of high-quality archival data available from the current generation of X-ray telescopes.

Using new XMM-Newton observations we detect hard X-ray time lags in the rapid variability of the Compton-thin Seyfert 2 galaxy IRAS 18325-5926. The higher-energy X-ray variations lag behind correlated lower-energy variations by up to ~3ks and the magnitude of the lag increases clearly with energy separation between the energy bands. We find that the lag-energy spectrum has a relatively simple log(E) shape. This is quite different in both shape and magnitude from the lags predicted by simple reflection models, but very similar to the hard X-ray lags often seen in black hole X-ray binaries. We apply several spectral models to the lag-energy spectrum and rule out simple reflection as an origin for the hard lags. We find that both propagating fluctuations embedded in the accretion flow and electron scattering from material embedded in or behind a cold absorbing medium offer equally good fits to the observed low-frequency hard X-ray lags and are both consistent with the time-averaged spectrum. Such models will likely look very different outside of XMM-Newton's observable bandpass, paving the way for future studies with NuSTAR.

MCG-6-30-15 is one of the most observed Narrow Line Seyfert 1 galaxies in the X-ray band. In this paper we examine the X-ray time lags in this source using a total of 600 ks in observations (440 ks exposure) taken with the XMM-Newton telescope (300 ks in 2001 and 300 ks in 2013). Both the old and new observations show the usual hard lag that increases with energy, however, the hard lag turns over to a soft lag at frequencies below ~1e-4 Hz. The highest frequencies (~1e-3 Hz) in this source show a clear soft lag, as previously presented for the first 300 ks observation, but no clear iron K lag is detected in either the old or new observation. The soft lag is more significant in the old observation than the new. The observations are consistent with a reverberation interpretation, where the soft, reflected emission is delayed with respect to the hard powerlaw component. These spectral timing results suggest that two distinct variability mechanisms are important in this source: intrinsic coronal variations (which lead to correlated variability in the reprocessed emission), and geometrical changes in the corona. Variability due to geometrical changes does not result in correlated variability in the reflection, and therefore inhibits the clear detection of an iron K lag.

The X-ray quasi-periodic oscillation (QPO) at $\sim 2.6 \times 10^{-4}$ Hz in RE J1034+396 has been robustly observed in the 0.2-10 keV band of one $\sim 90$ ks XMM-Newton observation, but was not detected in subsequent observations. Here, we investigate the power spectral density (PSD) of 8 archival XMM-Newton observations of RE J1034+396, and search for the presence of QPOs in three energy bands: soft (0.3-0.8 keV); hard (1.0-4.0 keV) and total (0.2-10.0 keV). We find a significant detection of a QPO feature in the hard-band PSD of 5 low flux/spectrally-harder observations. The QPO frequency has remained persistent at $\sim 2.6 \times 10^{-4}$ Hz over $\sim 4$ years of observations, though it is no longer detected in the soft band, except in the original observation. This result increases the duration where the QPO is present by a factor of three (now $\sim 250$ ks), allowing for a better understanding of the QPO phenomenon observed in both Active Galactic Nuclei (AGN) and black hole X-ray binaries (BHBs).

Emission from the active nucleus in the core of the brightest cluster galaxy of the Perseus cluster, NGC 1275, has varied dramatically over the past four decades. Prompted by the Fermi detection of flaring in the gamma-ray band, we present the recent increased activity of this source in the context of its past radio and gamma-ray output. The broad correspondence between the high-frequency radio data and the high-energy (HE) emission is striking. However, on short timescales this correlation breaks down and the 1.3 mm Submillimeter Array flux is apparently unaffected during Fermi-detected flaring activity. The fact that NGC 1275 is also detected at TeV energies during the periods of HE gamma-ray flaring suggests that the short-timescale variation might be primarily related to changes in the inverse Compton scattering of photons by the electron population in the jet. The longer-timescale changes suggest a 30--40 year variation in the fuelling of the black hole, that affects the power of the inner jet. NCG 1275 is a laboratory for the class of brightest cluster galaxies, and its variability on these timescales has implications for our understanding of massive galaxies in cooling-core clusters. The case of NGC 1275 highlights the need for wide coverage across the radio band to correctly account for the contribution to emission from a synchrotron self-absorbed core (for example when considering contamination of Sunyaev-Zel'dovich effect observations), and the danger of variability biases in radio surveys of galaxies.

We analyse the X-ray time-lags in the Narrow Line Seyfert 1 (NLS1) galaxy PG 1244+026 ($M_{\rm BH}\sim 10^7 M_{\rm sun}$, $L/L_{\rm Edd}\sim 1$). The time delay between the soft (0.3--0.7~keV) and harder (1.2--4.0~keV) variations shows shows the well established switch from hard lags at low frequencies to soft lags at high frequencies. The low frequency hard lags are qualitatively consistent with the propagation of fluctuations model, with some long-timescale response of the reflection component. The high frequency soft lag appears to extend over a wide frequency band, that we divide this into two narrow frequency ranges, and examine the lag as a function of energy for each of these. At high frequencies the soft excess emission is delayed with respect to the harder energy bands, without any corresponding strong, hard X-ray reflection signature. At even higher frequencies a soft lag is seen at the softest energies, as well as tentative evidence for an iron K$\alpha$ reverberation signal. These results point to the importance of reprocessing as well as reflection in determining the lags in NLS1s.

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.

We analyse 15 XMM-Newton observations of the Seyfert galaxy NGC 4051 obtained over 45 days to determine the ultraviolet (UV) light curve variability characteristics and search for correlated UV/X-ray emission. The UV light curve shows variability on all time scales, however with lower fractional rms than the 0.2-10 keV X-rays. On days-weeks timescales the fractional variability of the UV is Fvar ~ 8%, and on short (~ hours) timescales Fvar ~ 2%. The within-observation excess variance in 4 of the 15 UV observations was found be much higher than the remaining 11. This was caused by large systematic uncertainties in the count rate masking the intrinsic source variance. For the four "good" observations we fit an unbroken power-law model to the UV power spectra with slope -2.0 +/- 0.5. We compute the UV/X-ray Cross-correlation function for the "good" observations and find a correlation of ~ 0.5 at time lag of ~ 3 ks, where the UV lags the X-rays. We also compute for the first time the UV/X-ray Cross-spectrum in the range 0-28.5 ks, and find a low coherence and an average time lag of ~ 3 ks. Combining the 15 XMM-Newton and the Swift observations we compute the DCF over +/-40 days but are unable to recover a significant correlation. The magnitude and direction of the lag estimate from the 4 "good" observations indicates a scenario where ~ 25 % of the UV variance is caused by thermal reprocessing of the incident X-ray emission.