GX~339$-$4 is a black hole X-ray binary that is a key focus of accretion studies since it goes into outburst roughly every two-to-three years. Tracking of its radio, IR and X-ray flux during multiple outbursts reveals tight broadband correlations. The radio emission originates in a compact, self-absorbed jet, however the origin of the X-ray emission is still debated: jet base or corona? We fit 20 quasi-simultaneous radio, IR, optical and X-ray observations of GX~339$-$4 covering three separate outbursts in 2005, 2007, 2010--2011, with a composite corona + jet model, where inverse Compton emission from both regions contributes to the X-ray emission. Using a recently-proposed identifier of the X-ray variability properties known as power-spectral hue, we attempt to explain both the spectral and evolving timing characteristics, with the model. We find the X-ray spectra are best fit by inverse Compton scattering in a dominant hot corona ($kT_{\rm e}\sim$ hundreds of keV). However, radio and IR-optical constraints imply a non-negligible contribution from inverse Compton scattering off hotter electrons ($kT_{\rm e} \ge 511$~keV) in the base of the jets, ranging from a few up to $\sim50$\% of the integrated 3--100~keV flux. We also find that the physical properties of the jet show interesting correlations with the shape of the broadband X-ray variability of the source, posing intriguing suggestions for the connection between the jet and corona.

Shemar et al. (2016) presented results based on the output of a feasibility study for the European Space Agency (ESA) on the use of X-ray pulsars for deep space navigation, a concept often referred to as 'XNAV'. Here we describe some of the key results as well as providing additional information which includes navigation uncertainties and the potential X-ray technology that could be used. For a conventional deep space mission, an X-ray navigation system must be practical to implement as a spacecraft subsystem and to this end it must meet restrictive mass, volume and power consumption requirements. The implementation of an X-ray observatory sized instrument is unrealistic in this case. The Mercury Imaging X-ray Spectrometer (MIXS) instrument, due to be launched on the ESA/JAXA BepiColombo mission to Mercury in 2018, is an example of an instrument that may be further developed as a practical telescope for XNAV. Simulations involving different pulsar combinations and navigation strategies are used to estimate the navigation uncertainties that may be achievable using such an instrument. Possible options for future developments in terms of simpler, lower-cost Kirkpatrick-Baez optics are discussed, in addition to the principal design and development challenges that must be addressed in order to realise an operational XNAV system.

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.

Among the population of known galactic black hole X-ray binaries, GRS 1915+105 stands out in multiple ways. It has been in continuous outburst since 1992, and has shown a wide range of different states that can be distinguished by their timing and spectral properties. These states, also observed in IGR J17091-3624, have in the past been linked to accretion dynamics. Here, we present the first comprehensive study into the long-term evolution of GRS 1915+105, using the entire data set observed with RXTE over its sixteen-year lifetime. We develop a set of descriptive features allowing for automatic separation of states, and show that supervised machine learning in the form of logistic regression and random forests can be used to efficiently classify the entire data set. For the first time, we explore the duty cycle and time evolution of states over the entire sixteen-year time span, and find that the temporal distribution of states has significantly changed over the span of the observations. We connect the machine classification with physical interpretations of the phenomenology in terms of chaotic and stochastic processes.

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

Thanks to the large collecting area (3 x ~1500 cm$^2$ at 1.5 keV) and wide field of view (30' across in full field mode) of the X-ray cameras on board the European Space Agency X-ray observatory XMM-Newton, each individual pointing can result in the detection of hundreds of X-ray sources, most of which are newly discovered. Recently, many improvements in the XMM-Newton data reduction algorithms have been made. These include enhanced source characterisation and reduced spurious source detections, refined astrometric precision, greater net sensitivity and the extraction of spectra and time series for fainter sources, with better signal-to-noise. Further, almost 50\% more observations are in the public domain compared to 2XMMi-DR3, allowing the XMM-Newton Survey Science Centre (XMM-SSC) to produce a much larger and better quality X-ray source catalogue. The XMM-SSC has developed a pipeline to reduce the XMM-Newton data automatically and using improved calibration a new catalogue version has been produced from XMM-Newton data made public by 2013 Dec. 31 (13 years of data). Manual screening ensures the highest data quality. This catalogue is known as 3XMM. In the latest release, 3XMM-DR5, there are 565962 X-ray detections comprising 396910 unique X-ray sources. For the 133000 brightest sources, spectra and lightcurves are provided. For all detections, the positions on the sky, a measure of the quality of the detection, and an evaluation of the X-ray variability is provided, along with the fluxes and count rates in 7 X-ray energy bands, the total 0.2-12 keV band counts, and four hardness ratios. To identify the detections, a cross correlation with 228 catalogues is also provided for each X-ray detection. 3XMM-DR5 is the largest X-ray source catalogue ever produced. Thanks to the large array of data products, it is an excellent resource in which to find new and extreme objects.

ULX-7, in the northern spiral arm of M51, demonstrates unusual behaviour for an ultraluminous X-ray source, with a hard X-ray spectrum but very high short-term variability. This suggests that it is not in a typical ultraluminous state. We analyse the source using archival data from XMM-Newton, Chandra and NuSTAR, and by examining optical and radio data from HST and VLA. Our X-ray spectral analysis shows that the source has a hard power-law spectral shape with a photon index Gamma~1.5, which persists despite the source's X-ray luminosity varying by over an order of magnitude. The power spectrum of the source features a break at 6.5^{+0.5}_{-1.1}x10^-3 Hz, from a low-frequency spectral index of alpha_1=-0.1^{+0.5}_{-0.2} to a high-frequency spectral index of alpha_2=0.65^{+0.05}_{-0.14}, making it analogous to the low-frequency break found in the power spectra of low/hard state black holes (BHs). We can take a lower frequency limit for a corresponding high-frequency break to calculate a BH mass upper limit of 1.6x10^3 solar masses. Using the X-ray/radio fundamental plane we calculate another upper limit to the BH mass of 3.5x10^4 solar masses for a BH in the low/hard state. The hard spectrum, high rms variability and mass limits are consistent with ULX-7 being an intermediate-mass BH; however we cannot exclude other interpretations of this source's interesting behaviour, most notably a neutron star with an extreme accretion rate.

Recent evidence - in particular the hard X-ray spectra obtained by NuSTAR, and the large amplitude hard X-ray variability observed when ultraluminous X-ray sources (ULXs) show soft spectra - reveals that common ULX behaviour is inconsistent with known sub-Eddington accretion modes, as would be expected for an intermediate-mass black hole (IMBH). Instead, it appears that the majority of ULXs are powered by super-Eddington accretion onto stellar-mass black holes. Here, we will review work that delves deeper into ULX spectral-timing behaviour, demonstrating it remains consistent with the expectations of super-Eddington accretion. One critical missing piece from this picture is the direct detection of the massive, radiatively-driven winds expected from ULXs as atomic emission/absorption line features in ULX spectra; we will show it is very likely these have already been detected as residuals in the soft X-ray spectra of ULXs. Finally, we will discuss ULXs that do not appear to conform to the emerging ULX behaviour patterns. In particular we discuss the implications of the identification of a good IMBH candidate as a background QSO; and the confirmation of an IMBH/ULX candidate in the galaxy NGC 2276 via the radio/X-ray fundamental plane.

The lack of unambiguous detections of atomic features in the X-ray spectra of ultraluminous X-ray sources (ULXs) has proven a hindrance in diagnosing the nature of the accretion flow. The possible association of spectral residuals at soft energies with atomic features seen in absorption and/or emission and potentially broadened by velocity dispersion could therefore hold the key to understanding much about these enigmatic sources. Here we show for the first time that such residuals are seen in several sources and appear extremely similar in shape, implying a common origin. Via simple arguments we assert that emission from extreme colliding winds, absorption in a shell of material associated with the ULX nebula and thermal plasma emission associated with star formation are all highly unlikely to provide an origin. Whilst CCD spectra lack the energy resolution necessary to directly determine the nature of the features (i.e. formed of a complex of narrow lines or intrinsically broad), studying the evolution of the residuals with underlying spectral shape allows for an important, indirect test for their origin. The ULX NGC 1313 X-1 provides the best opportunity to perform such a test due to the dynamic range in spectral hardness provided by archival observations. We show through highly simplified spectral modelling that the strength of the features (in either absorption or emission) appears to anti-correlate with spectral hardness, which would rule out an origin via reflection of a primary continuum and instead supports a picture of atomic transitions in a wind or nearby material associated with such an outflow.

We present an energy dependent X-ray variability study of the 2010 outburst of the black hole X-ray binary MAXI J1659-152 with the Swift X-ray Telescope (XRT). The broad-band noise components and the quasi periodic oscillations (QPO) observed in the power spectra show a strong and varied energy dependence. Combining Swift XRT data with data from the Rossi X-ray Timing Explorer, we report, for the first time, an rms spectrum (fractional rms amplitude as a function of energy) of these components in the 0.5-30 keV energy range. We find that the strength of the low-frequency component (< 0.1 Hz) decreases with energy, contrary to the higher frequency components (> 0.1 Hz) whose strengths increase with energy. In the context of the propagating fluctuations model for X-ray variability, we suggest that the low-frequency component originates in the accretion disk (which dominates emission below ~ 2 keV) and the higher frequency components are formed in the hot flow (which dominates emission above ~ 2 keV). As the properties of the QPO suggest that it may have a different driving mechanism, we investigate the Lense-Thirring precession of the hot flow as a candidate model. We also report on the QPO coherence evolution for the first time in the energy band below 2 keV. While there are strong indications that the QPO is less coherent at energies below 2 keV than above 2 keV, the coherence increases with intensity similar to what is observed at energies above 2 keV in other black-hole X-ray binaries.

We demonstrate a new method of variability classification using observations of black hole X-ray binaries. Using `power colours' -- ratios of integrated power in different Fourier frequency bands -- we can clearly differentiate different canonical black hole states as the objects evolve during outburst. We analyse (~ 2400) Rossi X-ray Timing Explorer observations of 12 transient low mass black hole X-ray binaries and find that the path taken around the power colour-colour diagram as the sources evolve is highly consistent from object to object. We discuss how the consistency observed in the power colour-colour diagram between different objects allows for easy state classification based on only a few observations, and show how the power-spectral shapes can be simply classified using a single parameter, the power-spectral `hue'. To illustrate the benefits of our simple model-independent approach, we show that the persistent high mass X-ray binary Cyg X-1 shows very similar power-spectral evolution to the transient black hole sources, with the main difference being caused by a combination of a lack of quasi-periodic oscillations and an excess of low-frequency power-law noise in the Cyg X-1 power spectra during the transitional state. We also compare the transient objects to the neutron star atoll source Aquila X-1, demonstrating that it traces a different path in the power colour-colour plot. Thus, power-colours could be an effective method to classify newly discovered X-ray binaries.

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.

Ultraluminous X-ray sources (ULXs) with luminosities lying between ~3x10^{39} - 2x10^{40} erg/s represent a contentious sample of objects as their brightness, together with a lack of unambiguous mass estimates for the vast majority of the central objects, leads to a degenerate scenario where the accretor could be a stellar remnant (black hole or neutron star) or intermediate mass black hole (IMBH). Recent, high-quality observations imply that the presence of IMBHs in the majority of these objects is unlikely unless the accretion flow somehow deviates strongly from expectation based on objects with known masses. On the other-hand, physically motivated models for super-critical inflows can re-create the observed X-ray spectra and their evolution, although have been lacking a robust explanation for their variability properties. In this paper we include the effect of a partially inhomogeneous wind that imprints variability onto the X-ray emission via two distinct methods. The model is heavily dependent on both inclination to the line-of-sight and mass accretion rate, resulting in a series of qualitative and semi-quantitative predictions. We study the time-averaged spectra and variability of a sample of well-observed ULXs, finding that the source behaviours can be explained by our model in both individual cases as well as across the entire sample, specifically in the trend of hardness-variability power. We present the covariance spectra for these sources for the first time, which shed light on the correlated variability and issues associated with modelling broad ULX spectra.

Quasi-periodic oscillations (QPOs) observed in the giant flares of magnetars are of particular interest due to their potential to open up a window into the neutron star interior via neutron star asteroseismology. However, only three giant flares have been observed. We therefore make use of the much larger data set of shorter, less energetic recurrent bursts. Here, we report on a search for QPOs in a large data set of bursts from the two most burst-active magnetars, SGR 1806-20 and SGR 1900+14, observed with the Rossi X-ray Timing Explorer (RXTE). We find a single detection in an averaged periodogram comprising 30 bursts from SGR 1806-20, with a frequency of 57 Hz and a width of 5 Hz, remarkably similar to a giant flare QPO observed from SGR 1900+14. This QPO fits naturally within the framework of global magneto-elastic torsional oscillations employed to explain the giant flare QPOs. Additionally, we uncover a limit on the applicability of Fourier analysis for light curves with low background count rates and strong variability on short timescales. In this regime, standard Fourier methodology and more sophisticated Fourier analyses fail in equal parts by yielding an unacceptably large number of false positive detections. This problem is not straightforward to solve in the Fourier domain. Instead, we show how simulations of light curves can offer a viable solution for QPO searches in these light curves.

The discovery of quasi-periodic oscillations (QPOs) in magnetar giant flares has opened up prospects for neutron star asteroseismology. The scarcity of giant flares makes a search for QPOs in the shorter, far more numerous bursts from Soft Gamma Repeaters (SGRs) desirable. In Huppenkothen et al (2013), we developed a Bayesian method for searching for QPOs in short magnetar bursts, taking into account the effects of the complicated burst structure, and have shown its feasibility on a small sample of bursts. Here, we apply the same method to a much larger sample from a burst storm of 286 bursts from SGR J1550-5418. We report a candidate signal at 260 Hz in a search of the individual bursts, which is fairly broad. We also find two QPOs at 93 Hz and one at 127 Hz, when averaging periodograms from a number of bursts in individual triggers, at frequencies close to QPOs previously observed in magnetar giant flares. Finally, for the first time, we explore the overall burst variability in the sample, and report a weak anti-correlation between the power-law index of the broadband model characterising aperiodic burst variability, and the burst duration: shorter bursts have steeper power law indices than longer bursts. This indicates that longer bursts vary over a broader range of time scales, and are not simply longer versions of the short bursts.

The luminosities of ultraluminous X-ray sources (ULXs) require an exotic solution with either super-critical accretion modes onto stellar mass black holes or sub-critical accretion onto intermediate mass black holes (IMBHs) being invoked. Discriminating between the two is non-trivial due to the present lack of a direct mass measurement. A key expectation of the super-critical mode of accretion is the presence of powerful radiatively-driven winds. Here we analyse XMM-Newton data from NGC 5408 X-1 and NGC 6946 X-1 and find that strong soft residuals present in the X-ray spectra can be reconciled with broadened, blue-shifted absorption by a partially ionised, optically thin phase of this wind. We derive initial values for the physical parameters of the wind; we also discuss other possible origins for the observed features.

We have investigated the short term linear relation between the rms variability and the flux in 1,961 observations of 9 black hole X-ray binaries. The rms-flux relation for the 1-10 Hz range is ubiquitously observed in any observation with good variability signal to noise (> 3 % 1-10 Hz fractional rms). This concurs with results from a previous study of Cygnus X-1 (Gleissner et. al. 2004), and extends detection of the rms-flux relation to a wider range of states. We find a strong dependence of the flux intercept of the rms-flux relation on source state; as the source transitions from the hard state into the hard intermediate state the intercept becomes strongly positive. We find little evidence for flux dependence of the broad-band noise within the PSD shape, excepting a small subset of observations from one object in an anomalous soft-state. We speculate that the ubiquitous linear rms-flux relation in the broad band noise of this sample, representing a range of different states and objects, indicates that its formation mechanism is an essential property of the luminous accretion flow around black holes.

We present an analysis of the short timescale variations in the properties of the strong (type "C") quasi-periodic oscillation observed in XTE J1550-564 during its 1998 outburst. In particular, the QPO shows a correlation between absolute rms amplitude and mean source flux over timescales shorter than ~ 3 ksec. A linear rms-flux relation has been observed to be a common property of broad-band noise but here we report the first detection of rms-flux dependence in a QPO. The gradient of the rms-flux relation is correlated with the QPO peak frequency: from a strong positive correlation when the QPO peak frequency is below ~ 4 Hz, through no correlation, to a strong negative correlation when the peak frequency is above 6 Hz. This is the first time a negative short term rms-flux relation has been observed in any component of the power spectrum. Previous work on both the broad-band noise and QPOs in a range of sources have suggested the presence of a filter reducing the amplitude of QPOs with increasing frequency. We attempt to remove the possible effects of this filter and find that the previously negative rms-flux relations above ~ 5 Hz become constant.

We report the first detection of a linear correlation between rms variability amplitude and flux in the Ultraluminous X-ray source NGC 5408 X-1. The rms-flux relation has previously been observed in several Galactic black hole X-ray binaries (BHBs), several Active Galactic Nuclei (AGN) and at least one neutron star X-ray binary. This result supports the hypothesis that a linear rms-flux relation is common to all luminous black hole accretion and perhaps even a fundamental property of accretion flows about compact objects. We also show for the first time the cross-spectral properties of the variability of this ULX, comparing variations below and above 1 keV. The coherence and time delays are poorly constrained but consistent with high coherence between the two bands, over most of the observable frequency range, and a significant time delay (with hard leading soft variations). The magnitude and frequency dependence of the lags are broadly consistent with those commonly observed in BHBs, but the direction of the lag is reversed. These results indicate that ULX variability studies, using long X-ray observations, hold great promise for constraining the processes driving ULXs behaviour, and the position of ULXs in the scheme of black hole accretion from BHBs to AGN.

We present results from a study of short term variability in 19 archival observations by XMM-Newton of 16 Ultraluminous X-ray Sources (ULXs). Eight observations (six sources) showed intrinsic variability with power spectra in the form of either a power law or broken power law-like continuum and in some cases quasi-periodic oscillations (QPOs). The remaining observations were used to place upper limits on the strength of possible variability hidden within. Seven observations (seven sources) yielded upper limits comparable to, or higher than, the values measured from those observations with detectable variations. These represented the seven faintest sources all with f_x < 3x10^-12 erg cm^-2 s^-1. In contrast there are four observations (three sources) that gave upper limits significantly lower than both the values measured from the ULX observations with detectable variations, and the values expected by comparison with luminous Galactic black hole X-ray binaries (BHBs) and Active Galactic Nuclei (AGN) in the observed frequency bandpass (10^-3 - 1 Hz). This is the case irrespective of whether one assumes characteristic frequencies appropriate for a stellar mass (10 M_sun) or an intermediate mass (1000 M_sun) black hole, and means that in some ULXs the variability is significantly suppressed compared to bright BHBs and AGN. We discuss ways to account for this unusual suppression in terms of both observational and intrinsic effects and whether these solutions are supported by our results.