Source author record

Chris Done

Chris Done appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

80works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2024arXiv

Disc corona radii and QPO frequencies in black hole binaries: testing Lense-Thirring precession origin

Stellar-mass black hole binary systems in the luminous X-ray states show a strong quasi-periodic oscillation (QPO) in their Comptonised emission. The frequency of this feature correlates with the ratio of a disc to Comptonised emission rather than with total luminosity. Hence it changes dramatically during spectral transitions between the hard and soft states. Its amplitude is also strongest in these intermediate states, making them an important test of QPO models. However, these have complex spectra which generally require a disc and two separate Comptonisation components, making it difficult to uniquely derive the spectral parameters. We build a new energy-conserving model of the accretion flow, SSsed model, which assumes a fixed radial emissivity but with a changing emission mechanism. This is similar to the agnsed model in xspec but tuned to be more suitable for stellar mass black holes. It uses a combination of the disc luminosity and temperature to constrain the inner radius of the (colour temperature corrected) blackbody disc, separating this from the more complex Comptonisation spectra emitted inwards of this radius. We show a pilot study of this model fit to hundreds of RXTE spectra of the black hole binary XTE J1550-564. We show that the derived disc radius tightly anti-correlates with the central frequencies of the low-frequency QPO detected in the same observations. The relation is consistent with the quantitative predictions of Lense-Thirring precession of the entire inner Comptonisation regions for the assumed system parameters. This supports the scenario that low-frequency QPOs are caused by Lense-Thirring precession.

preprint2023arXiv

The radio detection and accretion properties of the peculiar nuclear transient AT 2019avd

AT 2019avd is a nuclear transient detected from infrared to soft X-rays, though its nature is yet unclear. The source has shown two consecutive flaring episodes in the optical and the infrared bands and its second flare was covered by X-ray monitoring programs. During this flare, the UVOT/Swift photometries revealed two plateaus: one observed after the peak and the other one appeared ~240 days later. Meanwhile, our NICER and XRT/Swift campaigns show two declines in the X-ray emission, one during the first optical plateau and one 70-90 days after the optical/UV decline. The evidence suggests that the optical/UV could not have been primarily originated from X-ray reprocessing. Furthermore, we detected a timelag of ~16-34 days between the optical and UV emission, which indicates the optical likely comes from UV reprocessing by a gas at a distance of 0.01-0.03 pc. We also report the first VLA and VLBA detection of this source at different frequencies and different stages of the second flare. The information obtained in the radio band - namely a steep and a late-time inverted radio spectrum, a high brightness temperature and a radio-loud state at late times - together with the multiwavelength properties of AT 2019avd suggests the launching and evolution of outflows such as disc winds or jets. In conclusion, we propose that after the ignition of black hole activity in the first flare, a super-Eddington flaring accretion disc formed and settled to a sub-Eddington state by the end of the second flare, associated with a compact radio outflow.

preprint2022arXiv

A full spectral-timing model to map the accretion flow in black hole binaries: the low/hard state of MAXI J1820+070

The nature and geometry of the accretion flow in the low/hard state of black hole binaries is currently controversial. While most properties are generally explained in the truncated disc/hot inner flow model, the detection of a broad residual around the iron line argues for strong relativistic effects from an untruncated disc. Since spectral fitting alone is somewhat degenerate, we combine it with the additional information in the fast X-ray variability and perform a full spectral-timing analysis for NICER and NuSTAR data on a bright low/hard state of MAXI J1820+070. We model the variability with propagating mass accretion rate fluctuations by combining two separate current insights: that the hot flow is spectrally inhomogeneous, and that there is a discontinuous jump in viscous time-scale between the hot flow and variable disc. Our model naturally gives the double-humped shape of the power spectra, and the increasing high-frequency variability with energy in the second hump. Including reflection and reprocessing from a disc truncated at a few tens of gravitational radii quantitatively reproduces the switch in the lag-frequency spectra, from hard lagging soft at low frequencies (propagation through the variable flow) to the soft lagging hard at the high frequencies (reverberation from the hard X-ray continuum illuminating the disc). The viscous time-scale of the hot flow is derived from the model, and we show how this can be used to observationally test ideas about the origin of the jet.

preprint2022arXiv

MAXI J1820+070 X-ray spectral-timing reveals the nature of the accretion flow in black hole binaries

Black hole X-ray binaries display significant stochastic variability on short time-scales (0.01-100 seconds), with a complex pattern of lags in correlated variability seen in different energy bands. This behaviour is generally interpreted in a model where slow fluctuations stirred up at large radii propagate down through the accretion flow, modulating faster fluctuations generated at smaller radii. Coupling this scenario with radially-stratified emission opens the way to measure the propagation time-scale from data, allowing direct tests of the accretion flow structure. We previously developed a model based on this picture and showed that it could fit NICER (0.5-10 keV) data from the brightest recent black hole transient, MAXI J1820+070. However, here we show it fails when extrapolated to higher energy variability data from Insight-HXMT. We extend our model so that the spectrum emitted at each radius changes shape in response to fluctuations (pivoting) rather than just changing normalisation. This gives the strong suppression of fractional variability as a function of energy seen in the data. The derived propagation time-scale is slower than predicted by a magnetically arrested disc (MAD), despite this system showing a strong jet. Our new model jointly fits the spectrum and variability up to 50 keV, though still cannot match all the data above this. Nonetheless, the good fit for 3-40 keV means the QPO can most easily be explained as an extrinsic modulation of the flow, such as produced in Lense-Thirring precession, rather than arising in an additional spectral-timing component such as the jet.

preprint2022arXiv

Multi-wavelength Campaign on the Super-Eddington NLS1 RX J0134.2-4258 -- I. Peculiar X-ray Spectra and Variability

We have conducted a new long-term multi-wavelength campaign on one of the most super- Eddington narrow-line Seyfert 1s (NLS1s) known, namely RX J0134.2-4258. In this first paper, we report deep simultaneous X-ray observations performed by XMM-Newton and NuSTAR in 2019-12-19, during which RX J0134.2-4258 was fortuitously at one of its lowest X-ray flux states. However, there is a clear rise above 4 keV which implies that the intrinsic source flux may be higher. The X-ray spectra observed between 1996 and 2019 show drastic variability, probably due to complex, variable absorption along the line of sight. Unusually, the soft X-ray excess appears extremely weak in all these spectra, even when the hard X-ray spectrum has a steep slope of $Γ\simeq 2.2$. We explore the spectral-timing properties of the new (low X-ray flux) and archival (high X-ray flux) XMM-Newton data, fitting their time-average, rms and lag spectra simultaneously. The variability spectra indicate the presence of a very weak soft X-ray Comptonisation component, whose shape is similar to the soft excess in normal super-Eddington NLS1s, but with flux relative to the power law which is lower by more than one order of magnitude. Above 4 keV the low-flux data are dominated by a different component, which lags with respect to the lower energy emission. This is consistent with an origin of reflection or partial covering absorption from low ionization material located within 100 Rg. We interpret this as a further indication of the presence of a clumpy disc wind.

preprint2020arXiv

Geometry of the X-ray source 1H 0707-495

We investigate constraints for the size and location of the X-ray source in 1H 0707-495 determined from the shape of the relativistically smeared reflection from the accretion disc. We develop a new code to model an extended X-ray source and we apply it to all archival XMM observations of 1H 0707-495. Contrary to Wilkins et al. we find that the relativistic reflection in this source is not consistent with an extended uniform corona. Instead, we find that the X-ray source must be very compact, with the size of at most a gravitational radius, and located at most at a few gravitational radii from the black hole horizon. A uniform extended corona does indeed produce an emissivity which is like a twice broken power law, but the inner emissivity is fixed by the source geometry rather than being a free parameter. In 1H0707-495, reflection from the inner disc is much stronger than expected for a uniformly extended source. Including the effect of ionised absorption from a wind does not change this conclusion, but including scattered emission (and more complex absorption) from the wind can dramatically change the reflection parameters.

preprint2020arXiv

Qwind code release: a non-hydrodynamical approach to modelling line-driven winds in active galactic nuclei

Ultraviolet (UV) line driven winds may be an important part of the active galactic nucleus (AGN) feedback process, but understanding their impact is hindered by the complex nature of the radiation hydrodynamics. Instead, we have taken the approach pioneered by Risaliti & Elvis, calculating only ballistic trajectories from radiation forces and gravity, but neglecting gas pressure. We have completely re-written their QWIND code using more robust algorithms, and can now quickly model the acceleration phase of these winds for any AGN spectral energy distribution spanning UV and X-ray wavebands. We demonstrate the code using an AGN with black hole mass $10^8\, M_\odot$ emitting at half the Eddington rate and show that this can effectively eject a wind with velocities $\simeq (0.1-0.2)\, c$. The mass loss rates can be up to $\simeq 0.3 M_\odot$ per year, consistent with more computationally expensive hydrodynamical simulations, though we highlight the importance of future improvements in radiation transfer along the multiple different lines of sight illuminating the wind. The code is fully public, and can be used to quickly explore the conditions under which AGN feedback can be dominated by accretion disc winds.

preprint2020arXiv

Re-observing the NLS1 Galaxy RE J1034+396. I. the Long-term, Recurrent X-ray QPO with a High Significance

RE J1034+396 is a narrow-line Seyfert 1 galaxy (NLS1) in which the first significant X-ray quasi-periodic oscillation (QPO) in an active galactic nuclei (AGN) was observed in 2007. We report the detection of this QPO in a recent XMM-Newton observation in 2018 with an even higher significance. The quality factor of this QPO is 20, and its period is 3550 $\pm$ 80 s, which is 250 $\pm$ 100 s shorter than in 2007. While the QPO's period has no significant energy dependence, its fractional root-mean-square (rms) variability increases from 4% in 0.3-1 keV to 12% in 1-4 keV bands. An interesting phenomenon is that the QPO in 0.3-1 keV leads that in the 1-4 keV bands by 430 $\pm$ 50 s with a high coherence, opposite to the soft X-ray lag reported for the observation in 2007. We speculate that the QPO has an intrinsic hard lag, while the previous reported soft lag is caused by the interference of stochastic variability. This soft X-ray lead in the new data supports the idea that the QPO of RE J1034+396 is a possible AGN counterpart of the 67 Hz high-frequency QPO seen in the black hole binary (BHB) GRS 1915+105. We also search for QPO harmonics, but do not find any significant signals. Our new data reinforce previous results that the QPO is seen in a specific spectral state, as the only 2 observations showing no significant QPO signal exhibit an even stronger soft X-ray excess than the other 6 observations which display the QPO. Therefore, our results imply that the QPO in RE J1034+396 is physically linked to a soft X-ray component.

preprint2020arXiv

Reverberation Reveals the Truncated Disc in the Hard State of GX 339-4

The nature and geometry of the hard state in black hole binaries is controversial. The broadband continuum spectrum and fast variability properties can be explained in a model where the inner disc evaporates into a geometrically thick, hot flow. However these models are challenged by the persistent detection of an extremely broad iron line, which requires that the disc extends down to the last stable orbit of a high spin black hole. This line width can be considerably reduced if the Comptonisation continuum is multi-component rather than single temperature, but such models are highly degenerate. Here we show a specific model of a radially stratified continuum coupled to a model of propagating fluctuations, fit to some of the best hard state data from GX 339-4. This full spectral-timing model can fit the time averaged spectrum, the power spectra in different energy bands, and the frequency dependent lags between these bands. For the first time we also include disc reverberation and show that this same spectral-timing model successfully \textit{predicts} the lag-energy spectra on all timescales. This gives a more robust method to determine the inner radius of the disc, which is of order $20~R_g$, i.e. significantly truncated. This opens up the way to use the fast variability spectral-timing data to trace the source geometry of black hole binaries in all states.

preprint2020arXiv

State-of-the-art AGN SEDs for Photoionization Models: BLR Predictions Confront the Observations

The great power offered by photoionization models of Active Galactic Nuclei (AGN) emission-line regions has long been mitigated by the fact that very little is known about the spectral energy distribution (SED) between the Lyman limit, where intervening absorption becomes a problem, and 0.3 keV, where soft x-ray observations become possible. The emission lines themselves can, to some degree, be used to probe the SED, but only in the broadest terms. This paper employs a new generation of theoretical SEDs which are internally self-consistent, energy-conserving, and tested against observations, to infer properties of the emission-line regions. The SEDs are given as a function of the Eddington ratio, allowing emission-line correlations to be investigated on a fundamental basis. We apply the simplest possible tests, based on the foundations of photoionization theory, to investigate the implications for the geometry of the emission-line region. The SEDs become more far-ultraviolet bright as the Eddington ratio increases, so the equivalent widths of recombination lines should also become larger, an effect which we quantify. The observed lack of correlation between Eddington ratio and equivalent width shows that the cloud covering factor must decrease as Eddington ratio increases. This would be consistent with recent models proposing that the broad-line region is a failed dusty wind off the accretion disc.

preprint2020arXiv

The thermal-radiative wind in low mass X-ray binary H 1743-322: II. iron line predictions from Monte Carlo radiation transfer

We show the best current simulations of the absorption and emission features predicted from thermal-radiative winds produced from X-ray illumination of the outer accretion disc in binary systems. We use the density and velocity structure derived from a radiation hydrodynamic code as input to a Monte-Carlo radiation transport calculation. The initial conditions are matched to those of the black hole binary system H1743-322 in its soft, disc dominated state, where wind features are seen in Chandra grating data. Our simulation fits well to the observed line profile, showing that these physical wind models can be the origin of the absorption features seen, rather than requiring a magnetically driven wind. We show how the velocity structure is the key observable discriminator between magnetic and thermal winds. Magnetic winds are faster at smaller radii, whereas thermal winds transition to a static atmosphere at smaller radii. New data from XRISM (due for launch Jan 2022) will give an unprecedented view of the physics of the wind launch and acceleration processes, but the existence of static atmospheres in small disc systems already rules out magnetic winds which assume self-similar magnetic fields from the entire disc as the origin of the absorption features seen.

preprint2020arXiv

The thermal-radiative wind in the neutron star low mass X-ray binary GX 13+1

We fit the observed high ionisation X-ray absorption lines in the neutron star binary GX13+1 with a full simulation of a thermal-radiative wind. This uses a radiation hydrodynamic code coupled to Monte Carlo radiation transfer to compute the observed line profiles from Hydrogen and Helium-like iron and Nickel, including all strong Kα and Kβ transitions. The wind is very strong as this object has a very large disc and is very luminous. The absorption lines from Fe Kα are strongly saturated as the ion columns are large, so the line equivalent widths (EWs) depend sensitively on the velocity structure. We additionally simulate the lines including isotropic turbulence at the level of the azimuthal and radial velocities. We fit these models to the Fe xxv and xxvi absorption lines seen in the highest resolution Chandra third order HETGS data. These data already rule out the addition of turbulence at the level of the radial velocity of ~500 km/s. The velocity structure predicted by the thermal-radiative wind alone is a fairly good match to the observed profile, with an upper limit to additional turbulence at the level of the azimuthal velocity of ~100 km/s. This gives stringent constraints on any remaining contribution from magnetic acceleration.

preprint2019arXiv

Discarding the disc in a changing state AGN: the UV/X-ray relation in NGC 4151

Recent monitoring campaigns designed to map the accretion regime in AGN show major discrepancies with models where the optical/ultraviolet (UV) is produced by X-ray-illuminated, optically thick disc material within a few hundred gravitational radii. However, these campaigns only monitored X-rays below $10$ keV, whereas the bolometric luminosity for most of these AGN peaks above $50$ keV. We use data from the recent multiwavelength campaign by \cite{E17} on NGC 4151 - the only AGN bright enough to be monitored at higher energies with \textit{Swift} BAT. We develop a spectral-timing model with a hot corona, warm Comptonisation, and outer standard disc. This fits the time-averaged spectrum well, but completely fails to match the UV variability predicted from the X-ray lightcurve. However, it reveals that NGC 4151 had a bolometric luminosity around $1.4\%$ of the Eddington luminosity during this campaign, close to the luminosity at which AGN show a `changing-state' transition, where the broad optical lines disappear. Stellar mass black holes show a state transition at a similarly low Eddington fraction, which is broadly interpreted as the inner disc being replaced by an optically thin flow. We find that the UV lightcurve can instead be matched by reprocessing of the X-ray flux on size scales of the broad line region (BLR; $1.5-20$ light-days) and rule out there being optically thick material inwards of this, as expected if the thin disc is replaced by the flow below the inner radius of the BLR. These results emphasise the need for even longer-timescale, multiwavelength monitoring campaigns on variable AGN.

preprint2016arXiv

A disk wind interpretation of the strong Fe Kα features in 1H 0707-495

1H 0707-495 is the most convincing example of a supermassive black hole with an X-ray spectrum being dominated by extremely smeared, relativistic reflection, with the additional requirement of strongly supersoler iron abundance. However, here we show that the iron features in its 2--10 keV spectrum are rather similar to the archetypal wind dominated source, PDS 456. We fit all the 2--10 keV spectra from 1H 0707-495 using the same wind model as used for PDS 456, but viewed at higher inclination so that the iron absorption line is broader but not so blueshifted. This gives a good overall fit to the data from 1H 0707-495, and an extrapolation of this model to higher energies also gives a good match to the NuSTAR data. Small remaining residuals indicate that the iron line emission is stronger than in PDS 456. This is consistent with the wider angle wind expected from a continuum driven wind from the super-Eddington mass accretion rate in 1H 0707-495, and/or the presence of residual reflection from the underlying disk though the presence of the absorption line in the model removes the requirement for highly relativistic smearing, and highly supersoler iron abundance. We suggest that the spectrum of 1H 0707-495 is sculpted more by absorption in a wind than by extreme relativistic effects in strong gravity.

preprint2016arXiv

A quasi-periodic modulation of the iron line centroid energy in the black hole binary H 1743-322

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.

preprint2016arXiv

A variable ULX and possible IMBH candidate in M51a

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.

preprint2016arXiv

An optically-thick disk wind in GRO J1655-40?

We revisited the unusual wind in GRO J1655-40 detected with Chandra in 2005 April, using long-term RXTE X-ray data and simultaneous optical/near-infrared photometric data. This wind is the most convincing case for magnetic driving in black hole binaries, as it has an inferred launch radius which is a factor of 10 smaller than the thermal wind prediction. However, the optical and near-infrared fluxes monotonically increase around the Chandra observation, whereas the X-ray flux monotonically decreases from 10 days beforehand. Yet the optical and near-infrared fluxes are from the outer, irradiated disk, so for them to increase implies that the X-rays likewise increased. We applied a new irradiated disk model to the multi-wavelength spectral energy distributions (SEDs). Fitting the optical and near-infrared fluxes, we estimated the intrinsic luminosity at the Chandra epoch was > ~0.7 L_Edd, which is more than one order of magnitude larger than the observed X-ray luminosity. These results could be explained if a Compton-thick, almost completely ionized gas was present in the wind and strong scattering reduced the apparent X-ray luminosity. The effects of scattering in the wind should then be taken into account for discussion of the wind-driving mechanism. Radiation pressure and Compton heating may also contribute to powering the wind at this high luminosity.

preprint2016arXiv

Energy spectra of X-ray quasi-periodic oscillations in the Lense-Thirring precession model

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

preprint2016arXiv

Reaching the Peak of the quasar spectral energy distribution - II. Exploring the accretion disc, dusty torus and host galaxy

We continue our study of the spectral energy distributions (SEDs) of 11 AGN at 1.5 < z < 2.2, with optical-NIR spectra, X-ray data and mid-IR photometry. In a previous paper we presented the observations and models; in this paper we explore the parameter space of these models. We first quantify uncertainties on the black hole masses (M$_{\rm BH}$) and degeneracies between SED parameters. The effect of BH spin is tested, and we find that while low to moderate spin values (a$_*$ $\leq$ 0.9) are compatible with the data in all cases, maximal spin (a$_*$ = 0.998) can only describe the data if the accretion disc is face-on. The outer accretion disc radii are well constrained in 8/11 objects, and are found to be a factor ~5 smaller than the self-gravity radii. We then extend our modelling campaign into the mid-IR regime with WISE photometry, adding components for the host galaxy and dusty torus. Our estimates of the host galaxy luminosities are consistent with the M$_{\rm BH}$-bulge relationship, and the measured torus properties (covering factor and temperatures) are in agreement with earlier work, suggesting a predominantly silicate-based grain composition. Finally, we deconvolve the optical-NIR spectra using our SED continuum model. We claim that this is a more physically motivated approach than using empirical descriptions of the continuum such as broken power-laws. For our small sample, we verify previously noted correlations between emission linewidths and luminosities commonly used for single-epoch M$_{\rm BH}$ estimates, and observe a statistically significant anti-correlation between [O III] equivalent width and AGN luminosity.

preprint2016arXiv

Testing the Completeness of the SDSS Colour Selection for Ultramassive, Slowly Spinning Black Holes

We investigate the sensitivity of the colour-based quasar selection algorithm of the Sloan Digital Sky Survey to several key physical parameters of supermassive black holes (SMBHs), focusing on BH spin ($a_{\star}$) at the high BH-mass regime ($M_{BH} \geqslant10^9\, M_{\odot}$). We use a large grid of model spectral energy distribution, assuming geometrically-thin, optically-thick accretion discs, and spanning a wide range of five physical parameters: BH mass $M_{BH}$, BH spin $a_{\star}$, Eddington ratio $L / L_{Edd}$ , redshift $z$, and inclination angle $inc$. Based on the expected fluxes in the SDSS imaging ugriz bands, we find that $\sim 99.8\%$ of our models with $M_{BH} \leqslant 10^{9.5}\, M_{\odot}$ are selected as quasar candidates and thus would have been targeted for spectroscopic follow-up. However, in the extremely high-mass regime, $\geqslant 10^{10} M_{\odot}$, we identify a bias against slowly/retrograde spinning SMBHs. The fraction of SEDs that would have been selected as quasar candidates drops below $\sim50\%$ for $a_{\star} <0$ across $0.5<z<2$. For particularly massive BHs, with $M_{BH} \simeq 3\times10^{10}\, M_{\odot}$, this rate drops below $\sim20\%$, and can be yet lower for specific redshifts. We further find that the chances of identifying any hypothetical sources with $M_{BH} = 10^{11}\, M_{\odot}$ by colour selection would be extremely low at the level of $\sim 3\%$. Our findings, along with several recent theoretical arguments and empirical findings, demonstrate that the current understanding of the SMBH population at the high-$M_{BH}$, and particularly the low- or retrograde-spinning regime, is highly incomplete.

preprint2016arXiv

The Mass and Spin of The Extreme Narrow Line Seyfert 1 Galaxy 1H 0707-495 and Its Implications for The Trigger for Relativistic Jets

Relativistic reflection models of the X-ray spectrum of the `complex' Narrow Line Seyfert 1 (NLS1) 1H 0707-495 require a high spin, moderate inclination, low mass black hole. With these parameters fixed, the observed optical/UV emission directly determines the mass accretion rate through the outer disc and hence predicts the bolometric luminosity. This is $140-260~\times$ the Eddington limit. Such a disc should power a strong wind, and winds are generically expected to be clumpy. Changing inclination angle with respect to a clumpy wind structure gives a possible explanation for the otherwise puzzling difference between `complex' NLS1 such as 1H 0707-495 and `simple' ones like PG 1244+026. Lines of sight which intercept the wind show deep absorption features at iron from the hot phase of the wind, together with stochastic dips and complex absorption when the clumps occult the X-ray source (complex NLS1), whereas both these features are absent for more face-on inclination (simple NLS1). This geometry is quite different to the clean view of a flat disc which is assumed for the spin measurements in relativistic reflection models, so it is possible that even 1H 0707-495 has low spin. If so, this re-opens the simplest and hence very attractive possibility that high black hole spin is a necessary and sufficient condition to trigger highly relativistic (bulk Lorentz factor $\sim 10-15$) jets.

preprint2016arXiv

The Quiescent Intracluster Medium in the Core of the Perseus Cluster

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

preprint2016arXiv

Tracking the energetics of the non-thermal disc-corona-jet in the very high state GX 339-4

The dramatic hard-soft spectral transition in Black Hole Binaries is important as it is associated with the collapse of the jet and with the strongest low frequency QPOs. These transition spectra (intermediate and very high state: VHS) are complex, with soft but distinctly non-thermal Comptonisation which merges smoothly into the disc emission. Here we develop a physical model for the accretion flow which can accommodate all these features, with an outer standard disc, which can make a transition to an energetically coupled disc-corona region, and make a further transition to a hot inner flow which can be radiatively inefficient if required. The code explicitly uses fully relativistic emissivity (Novikov-Thorne), and all Comptonisation is calculated with a hybrid (thermal and non-thermal) electron distribution. We fit this to a VHS spectrum from GX 339-4. We show that the complex continuum curvature produced by a hybrid electron distribution is enough to remove the strong constraint on black hole spin derived from reflection using simpler Comptonisation models. More fundamentally, we show that the VHS cannot be fit with the same Novikov-Thorne emissivity which can fit the disc dominated spectrum but instead requires that the inner flow is somewhat radiatively inefficient. This is consistent with an accretion powered jet, but simultaneous radio data show that the jet has already collapsed at the time of our data. Instead, it could point to truncation of the inner flow at radii larger than the innermost stable circular orbit, as predicted by the Lense-Thirring QPO models.

preprint2015arXiv

Complex Narrow Line Seyfert 1s: High spin or high inclination?

Complex narrow line Seyfert 1s (NLS1s), such as 1H0707-495, differ from simple NLS1s like PG1244+026 by showing stronger broad spectral features at Fe K and larger amplitude flux variability. These are correlated: the strongest Fe K features are seen during deep dips in the light curves of complex NLS1s. There are two competing explanations for these features, one where a compact X-ray source on the spin axis of a highly spinning black hole approaches the horizon and the consequent strong relativistic effects focus the intrinsic flux onto the inner edge of a thin disc, giving a dim, reflection dominated spectrum. The other is that the deep dips are caused by complex absorption by clumps close to the hard X-ray source. The reflection dominated model is able to reproduce the very short 30s soft lag from reverberation seen in the complex NLS1 1H0707-495. However, it does not explain the characteristic switch to hard lags on longer timescales. Instead, a full model of propagating fluctuations coupled to reverberation can explain the switch in the simple NLS1 PG1244+026 using a low spin black hole. However PG1244+026 has a longer reverberation lag of $\sim 200$s. Here we extend the successful propagation-reverberation model for the simple NLS1 PG1244+026 to include the effect of absorption from clumps in a turbulent region above the disk. The resulting occultations of the inner accretion flow can introduce additional hard lags when relativistic effects are taken into account. This dilutes the soft lag from reverberation and shifts it to higher frequencies, making a smooth transition between the 200s lags seen in simple NLS1s to the 30s lags in complex NLS1s. These two classes of NLS1 could then be determined by inclination angle with respect to a clumpy, probably turbulent, failed wind structure on the disc.

preprint2015arXiv

Minimum X-ray source size for a lamp-post corona in light-bending models for AGN

The lamppost model is often used to describe the X-ray source geometry in AGN, where an infinitesimal point source is located on the black hole spin axis. This is especially invoked for Narrow Line Seyfert 1 (NLS1) galaxies, where an extremely broad iron line seen in episodes of low X-ray flux can both be explained by extremely strong relativistic effects as the source approaches the black hole horizon. The most extreme spectrum seen from the NLS1 1H0707-495 requires that the source is less than 1Rg above the event horizon in this geometry. However, the source must also be large enough to intercept sufficient seed photons from the disk to make the hard X-ray Compton continuum which produces the observed iron line/reflected spectrum. We use a fully relativistic ray tracing code to show that this implies that the source must be substantially larger than 1Rg in 1H0707-495 if the disk is the source of seed photons. Hence the source cannot fit as close as 1Rg to the horizon, so the observed spectrum and variability are not formed purely by effects of strong gravity but probably also by changes in corona and inner accretion flow geometry.

preprint2015arXiv

Reaching the Peak of the Quasar Spectral Energy Distribution I: Observations and Models

We perform a spectral analysis of a sample of 11 medium redshift (1.5 < z < 2.2) quasars. Our sample all have optical spectra from the SDSS, infrared spectra from GNIRS and TSPEC, and X-ray spectra from XMM-Newton. We first analyse the Balmer broad emission line profiles which are shifted into the IR spectra to constrain black hole masses. Then we fit an energy-conserving, three component accretion model of the broadband spectral energy distribution (SED) to our multi-wavelength data. Five out of the 11 quasars show evidence of an SED peak, allowing us to constrain their bolometric luminosity from these models and estimate their mass accretion rates. Based on our limited sample, we suggest that estimating bolometric luminosities from L_5100A and L_2-10keV may be unreliable, as has been also noted for a low-redshift, X-ray selected AGN sample.

preprint2015arXiv

The Accreting Black Hole Swift J1753.5-0127 from Radio to Hard X-Ray

(abridged) We report on multi-wavelength measurements of Swift J1753.5-0127 in the hard state at L=2.7e36 erg/s (assuming d=3 kpc) in 2014. The radio emission is optically thick synchrotron, presumably from a compact jet. We take advantage of the low extinction and model the near-IR to UV emission with a multi-temperature disk model. Assuming a BH mass of M_BH=5 Msun and a system inclination of 40 deg, the fits imply an inner radius for the disk of Rin/Rg>212 d_3 (5Msun/M_BH). The outer radius is R_out/R_g=90,000 d_3 (5Msun/M_BH), which corresponds to 6.6e10 d_3 cm, consistent with the expected size of the disk. The 0.5-240 keV spectrum measured by Swift/XRT, Suzaku, and NuSTAR is relatively well characterized by a power-law with a photon index of Gamma=1.722+/-0.003, but a significant improvement is seen when a second continuum component is added. Reflection is a possibility, but no iron line is detected, implying a low iron abundance. We are able to fit the entire SED with a multi-temperature disk component, a Comptonization component, and a broken power-law, representing the emission from the compact jet. The broken power-law cannot significantly contribute to the soft X-ray emission, and this may be related to why Swift J1753.5-0127 is an outlier in the radio/X-ray correlation. The broken power-law might dominate above 20 keV, which would constrain the break frequency to be between 2.4e10 Hz and 3.6e12 Hz. Although the fits to the full SED do not include significant thermal emission in the X-ray band, previous observations have consistently seen such a component, and we find that there is evidence at the 3.1-sigma level for a disk-blackbody component with a temperature of 150(+30)(-20) eV and an inner radius of 5-14 R_g. If this component is real, it might imply the presence of an inner accretion disk in addition to the strongly truncated (R_in>212 R_g) disk.

preprint2015arXiv

The Best Constraints on A Super-Eddington Accretion Flow: XMM-Newton Observations of An Intermediate-mass Black Hole

RX J1140.1+0307 (hereafter RX1140) is a Narrow Line Seyfert 1 (NLS1) with one of the lowest black hole masses known in an AGN (M ${\le} 10^6$ M$_{\odot}$). We show results from two new {\it XMM-Newton} observations, showing soft 2-10~keV spectra, a strong excess at lower energies, and fast X-ray variability as is typical of this class. The soft excess can be equally well fit by either low temperature Comptonisation or highly smeared, ionised reflection models, but we use a covariance analysis of the fast X-ray variability as well as lag and coherence spectra to show that the low temperature Comptonisation model gives a better description of the break in variability properties between soft and hard X-rays. Both models also require an additional component at the softest energies, as expected from the accretion disc. However, this inner disc spectrum does not join smoothly onto the variable optical and far UV emission (which should be produced in the outer disc) unless the mass is underestimated by an order of magnitude. The variable optical and far UV emission instead suggests that $L/L_{Edd}\sim 10$ through the outer disc, in which case advection and/or wind losses are required to explain the observed broadband spectral energy distribution. However, the similarity of the X-ray properties of RX1140 to other simple NLS1 such as PG 1244+026, RE J1034+396 and RX~J0136 means it is likely that these are also super-Eddington sources. This means their spectral energy distribution cannot be used to determine black hole spin despite appearing to be disc dominated. It also means that the accretion geometry close to the black hole is unlikely to be a flat disc as assumed in the new X-ray reverberation mapping techniques.

preprint2015arXiv

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

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

preprint2014arXiv

A compact, metal-rich, kpc-scale outflow in FBQS J0209-0438: Detailed diagnostics from HST/COS extreme UV observations

We present HST/COS observations of highly ionized absorption lines associated with a radio-loud QSO at $z=1.1319$. The absorption system has multiple velocity components, tracing gas that is largely outflowing from the QSO at velocities of a few 100 km s$^{-1}$. There is an unprecedented range in ionization, with detections of HI, NIII, NIV, NV, OIV, OIV*, OV, OVI, NeVIII, MgX, SV and ArVIII. We estimate the total hydrogen number density from the column density ratio N(OIV*)/N(OIV) to be $\log(n_{\textrm{H}}/\textrm{cm}^3)\sim 3$. Assuming photoionization equilibrium, we derive a distance to the absorbing complex of $2.3<R<6.0$ kpc from the centre of the QSO. A range in ionization parameter, covering $\sim 2$ orders of magnitude, suggest absorption path lengths in the range $10^{-4.5}<l_{\textrm{abs}}<1$ pc. In addition, the absorbing gas only partially covers the background emission from the QSO continuum, which suggests clouds with transverse sizes $l_{\textrm{trans}}<10^{-2.5}$ pc. Widely differing absorption path lengths, combined with covering fractions less than unity across all ions pose a challenge to models involving simple cloud geometries. These issues may be mitigated by the presence of non-equilibrium effects, together with the possibility of multiple gas temperatures. The dynamics and expected lifetimes of the gas clouds suggest that they do not originate from close to the AGN, but are instead formed close to their observed location. Their inferred distance, outflow velocities and gas densities are broadly consistent with scenarios involving gas entrainment or condensations in winds driven by either supernovae, or the supermassive black hole accretion disc. In the case of the latter, the present data most likely does not trace the bulk of the outflow by mass, which could instead manifest itself as an accompanying warm absorber, detectable in X-rays.

preprint2014arXiv

A physical model for the continuum variability and QPO in accreting black holes

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

preprint2014arXiv

A Physical Model for the X-ray Time Lags of Narrow Line Seyfert Type 1 Active Galactic Nuclei

We study the origin of the soft X-ray excess seen in the 'simple' Narrow Line Seyfert 1 galaxy PG1244+026 using all available spectral-timing information. This object shows the standard switch between soft leading the hard band on long timescales, to the opposite behaviour on short timescales. This is interpreted as a combination of intrinsic fluctuations propagating down through the accretion flow giving the soft lead, together with reflection of the hard X-rays giving the soft lag. We build a full model of the spectral and time variability including both propagation and reflection, and compare our model with the observed power spectra, coherence, covariance, lag-frequency and lag-energy spectra. We compare models based on a separate soft excess component with those based on reflection dominated soft emission. Reflection dominated spectra have difficulty reproducing the soft lead at low frequency since reflection will always lag. They also suffer from high coherence and nearly identical hard and soft band power spectra in disagreement with the observations. This is a direct result of the power law and reflection components both contributing to the hard and soft energy bands, and the small radii over which the relativistically smeared reflection is produced transmitting too much high frequency power into the soft band. Conversely, we find the separate soft excess models (where the inner disc radius is $>6R_g$) have difficulty reproducing the soft lag at high frequency, as reflected flux does not contribute enough signal to overwhelm the soft lead. However, reflection should also be accompanied by reprocessing and this should add to the soft excess at low energies. This model can quantitatively reproduce the switch from soft lead to soft lag seen in the data and reproduces well the observed power spectra and other timing features which reflection dominated models cannot.

preprint2014arXiv

Irradiated, colour-temperature-corrected accretion discs in ultraluminous X-ray sources

Although attempts have been made to constrain the stellar types of optical counterparts to ULXs, the detection of optical variability instead suggests that they may be dominated by reprocessed emission from X-rays which irradiate the outer accretion disc. Here, we report results from a combined X-ray and optical spectral study of a sample of ULXs, which were selected for having broadened disc-like X-ray spectra, and known optical counterparts. We simultaneously fit optical and X-ray data from ULXs with a new spectral model of emission from an irradiated, colour-temperature-corrected accretion disc around a black hole, with a central Comptonising corona. We find that the ULXs require reprocessing fractions of $\sim 10^{-3}$, which is similar to sub-Eddington thermal dominant state BHBs, but less than has been reported for ULXs with soft ultraluminous X-ray spectra. We suggest that the reprocessing fraction may be due to the opposing effects of self-shielding in a geometrically thick super-critical accretion disc, and reflection from far above the central black hole by optically thin material ejected in a natal super-Eddington wind. Then, the higher reprocessing fractions reported for ULXs with wind-dominated X-ray spectra may be due to enhanced scattering onto the outer disc via the stronger wind in these objects. Alternatively, the accretion discs in these ULXs may not be particularly geometrically thick, rather they may be similar in this regard to the thermal dominant state BHBs.

preprint2014arXiv

Scaling accretion flow models from BHB to AGN - Why doesn't it work?

Black holes depend only on mass and spin, while what we see from the accretion flow in steady state depends also on mass accretion rate and (weakly) inclination. Hence we should be able to scale the accretion flow properties from the stellar to the supermassive black holes. But the data show significant differences between these two types of systems, suggesting that we are missing some crucial physics in AGN. One of these differences is the soft X-ray excess which is seen ubiquitously in bright AGN, but only occasionally in BHB. Another is the much faster variability seen in the high energy tail of high mass accretion rate AGN compared to that seen in the tail of BHB. We show that while this variability is not understood, it can be used via the new spectral-timing techniques to constrain the nature of the soft X-ray excess. The coherence, lag-frequency and lag-energy results strongly support this being an additional low temperature Comptonisation component rather than extreme relativistically smeared reflection in the 'simple' Narrow Line Seyfert 1 PG1244+026.

preprint2014arXiv

Spectral and Timing Properties of the Black Hole X-ray Binary H 1743-322 in the Low/hard State Studied with Suzaku

We report on the results from Suzaku observations of the Galactic black hole X-ray binary H 1743-322 in the low/hard state during its outburst in 2012 October. We appropriately take into account the effects of dust-scattering to accurately analyze the X-ray spectra. The time-averaged spectra in the 1-200 keV band are dominated by a hard power-law component of a photon index of \approx 1.6 with a high-energy cutoff at \approx 60 keV, which is well described with Comptonization of the disk emission by hot corona. We estimate the inner disk radius from the multi-color disk component, and find that it is 1.3-2.3 times larger than the radius in the high/soft state. This suggests that the standard disk was not extended to the innermost stable circular orbit. A reflection component from the disk is detected with R = Ω/2π\approx 0.6 (Ωis the solid angle). We also successfully estimate the stable disk component in a way independent of the time-averaged spectral modeling, by analyzing short-term spectral variability on the \sim 1-sec timescale. A weak low-frequency quasi-periodic oscillation (LF QPO) at 0.1-0.2 Hz is detected, whose frequency is found to correlate with the X-ray luminosity and photon index. This result may be explained by the evolution of the disk truncation radius.

preprint2014arXiv

Suzaku Observation of the Black Hole Binary 4U 1630--47 in the Very High State

We report the results from an X-ray and near-infrared observation of the Galactic black hole binary 4U 1630--47 in the very high state, performed with {\it Suzaku} and IRSF around the peak of the 2012 September-October outburst. The X-ray spectrum is approximated by a steep power law, with photon index of 3.2, identifying the source as being in the very high state. A more detailed fit shows that the X-ray continuum is well described by a multi-color disc, together with thermal and non-thermal Comptonization. The inner disc appears slightly truncated by comparison with a previous high/soft state of this source, even taking into account energetic coupling between the disc and corona, although there are uncertainties due to the dust scattering correction. The near-infrared fluxes are higher than the extrapolated disc model, showing that there is a contribution from irradiation in the outer disk and/or the companion star at these wavelengths. Our X-ray spectra do not show the Doppler shifted iron emission lines indicating a baryonic jet which were seen four days previously in an XMM-Newton observation, despite the source being in a similar state. There are also no significant absorption lines from highly ionized irons as are seen in the previous high/soft state data. We show that the increased source luminosity is not enough on its own to make the wind so highly ionized as to be undetectable. This shows that the disc wind has changed in terms of its launch radius and/or density compared to the high/soft state.

preprint2014arXiv

The ASTRO-H X-ray Astronomy Satellite

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

preprint2014arXiv

The origin of ultra-fast outflows in AGN: Monte-Carlo simulations of the wind in PDS 456

Ultra-fast outflows (UFOs) are seen in many AGN, giving a possible mode for AGN feedback onto the host galaxy. However, the mechanism(s) for the launch and acceleration of these outflows are currently unknown, with UV line driving apparently strongly disfavoured as the material along the line of sight is so highly ionised that it has no UV transitions. We revisit this issue using the Suzaku X-ray data from PDS 456, an AGN with the most powerful UFO seen in the local Universe. We explore conditions in the wind by developing a new 3-D Monte-Carlo code for radiation transport. The code only handles highly ionised ions, but the data show the ionisation state of the wind is high enough that this is appropriate, and this restriction makes it fast enough to explore parameter space. We reproduce the results of earlier work, confirming that the mass loss rate in the wind is around 30% of the inferred inflow rate through the outer disc. We show for the first time that UV line driving is likely to be a major contribution to the wind acceleration. The mass loss rate in the wind matches that predicted from a purely line driven system, and this UV absorption can take place out of the line of sight. Continuum driving should also play a role as the source is close to Eddington. This predicts that the most extreme outflows will be produced from the highest mass accretion rate flows onto high mass black holes, as observed.

preprint2014arXiv

X-ray time delays in the Narrow Line Seyfert 1 galaxy PG 1244+026

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$α$ reverberation signal. These results point to the importance of reprocessing as well as reflection in determining the lags in NLS1s.

preprint2013arXiv

A Long XMM-Newton Observation of An Extreme Narrow Line Seyfert 1: PG 1244+026

We explore the origin of the strong soft X-ray excess in Narrow Line Seyfert 1 galaxies using spectral-timing information from a 120ks {\it XMM-Newton} observation of PG 1244+026. Spectral fitting alone cannot distinguish between a true additional soft X-ray continuum component and strongly relativistically smeared reflection, but both models also require a separate soft blackbody component. This is most likely intrinsic emission from the disc extending into the lowest energy X-ray bandpass. The {\it RMS} spectra on short timescales (200-5000s) contain both (non-disk) soft excess and power law emission. However, the spectrum of the variability on these timescales correlated with the 4-10 keV lightcurve contains only the power law. Together these show that there is fast variability of the soft excess which is independent of the 4-10 keV variability. This is inconsistent with a single reflection component making the soft X-ray excess as this necessarily produces correlated variability in the 4-10 keV bandpass. Instead, the {\it RMS} and covariance spectra are consistent with an additional cool Comptonisation component which does not contribute to the spectrum above 2 keV.

preprint2013arXiv

A new way to measure supermassive black hole spin in accretion disc dominated Active Galaxies

We show that disc continuum fitting can be used to constrain black hole spin in a subclass of Narrow Line Seyfert 1 (NLS1) AGN as their low mass and high mass accretion rate means that the disc peaks at energies just below the soft X-ray bandpass. We apply the technique to the NLS1 PG1244+026, where the optical/UV/X-ray spectrum is consistent with being dominated by a standard disc component. This gives a best estimate for black hole spin which is low, with a firm upper limit of $a_*<0.86$. This contrasts with the recent X-ray determinations of (close to) maximal black hole spin in other NLS1 based on relativistic smearing of the iron profile. While our data on PG1244+026 does not have sufficient statistics at high energy to give a good measure of black hole spin from the iron line profile, cosmological simulations predict that black holes with similar masses have similar growth histories and so should have similar spins. This suggests that there is a problem either in our understanding of disc spectra, or/and X-ray reflection or/and the evolution of black hole spin.

preprint2013arXiv

Accretion disk and ionized absorber of the 9.7-hour dipping black hole binary MAXI J1305-704

We report the results from X-ray studies of the newly discovered black hole candidate MAXI J1305-704 based on Suzaku and Swift observations in the low/hard and high/soft states, respectively. The long Suzaku observation shows two types of clear absorption dips, both of which recur on a dip interval of 9.74 +- 0.04 hours, which we identify with the orbital period. There is also partially ionized absorption in the non-dip (persistent) emission in both the high/soft state and, very unusually, the low/hard state. However, this absorption (in both states) has substantially lower ionization than that seen in other high inclination systems, where the material forms a homogeneous disk wind. Here instead the absorption is most probably associated with clumpy, compact structures associated with the dipping material, which we see uniquely in this source likely because we view it at a very large inclination angle. A large inclination angle is also favored, together with a low black hole mass, to explain the high disk temperature seen in the fairly low luminosity high/soft state, as Doppler boosting enhances the disk temperature at high inclination. The disk radius inferred from these data is significantly smaller than that of the soft component seen in the low/hard state, supporting models where the disk is truncated at low luminosities. We find, however, that the lack of variability power on time scales of ~50 sec in the Suzaku low/hard state data is difficult to explain, even with a low mass black hole.

preprint2013arXiv

An imperfect double: probing the physical origin of the low-frequency QPO and its harmonic in black hole binaries

We extract the spectra of the strong low-frequency quasi-periodic oscillation (QPO) and its harmonic during the rising phase of an outburst in the black-hole binary XTE J1550-564. We compare these frequency resolved spectra to the time-averaged spectrum and the spectrum of the rapid (<0.1s) variability. The spectrum of the time averaged emission can be described by a disc, a Compton upscattered tail, and its reflection. The QPO spectrum contains no detectable disc, and the Compton spectrum is generally harder than in the time averaged emission, and shows less reflection, making it very similar to the spectrum of the rapid variability. The harmonic likewise contains no detectable disc component, but has a Compton spectrum which is systematically softer than the QPO, softer even than the Compton tail in the time averaged emission. We interpret these results in the context of the Lense-Thirring model for the QPO, where a precessing hot flow replaces the inner disc, and the harmonic is produced by the angular dependence of Compton scattering within the hot flow. We extend these models to include stratification of the hot flow, so that it is softer (lower optical depth) at larger radii closer to the truncated disc, and harder (higher optical depth) in the innermost parts of the flow where the rapid variability is produced. The different optical depth with radius gives rise to different angular dependence of the Comptonised emission, weighting the fundamental to the inner parts of the hot flow, and the harmonic to the outer. This is the first model which can explain both the spectrum of the QPO and its harmonic in a self consistent geometry.

preprint2013arXiv

Evidence for a Cool Disk and Inhomogeneous Coronae from Wide-band Temporal Spectroscopy of Cygnus X-1 with Suzaku

Unified X-ray spectral and timing studies of Cygnus X-1 in the low/hard and hard intermediate state were conducted in a model-independent manner, using broadband Suzaku data acquired on 25 occasions from 2005 to 2009 with a total exposure of ~ 450 ks. The unabsorbed 0.1--500 keV source luminosity changed over 0.8--2.8% of the Eddington limit for 14.8 solar masses. Variations on short (1--2 seconds) and long (days to months) time scales require at least three separate components: a constant component localized below ~2 keV, a broad soft one dominating in the 2--10 keV range, and a hard one mostly seen in 10--300 keV range. In view of the truncated disk/hot inner flow picture, these are respectively interpreted as emission from the truncated cool disk, a soft Compton component, and a hard Compton component. Long-term spectral evolution can be produced by the constant disk increasing in temperature and luminosity as the truncation radius decreases. The soft Compton component likewise increases, but the hard Compton does not, so that the spectrum in the hard intermediate state is dominated by the soft Compton component; on the other hand, the hard Compton component dominates the spectrum in the dim low/hard state, probably associated with a variable soft emission providing seed photons for the Comptonization.

preprint2013arXiv

Jets and the Accretion Flow in Low Luminosity Black Holes

The X-ray spectra of black hole binaries in the low/hard state first harden as the flux decreases, then soften. This change in behaviour has been variously attributed to either the X-rays switching from being produced in the flow to being dominated by the jet, or to the flow switching seed photons from the disc to self generated seed photons from cyclo-synchrotron. Here we build a simple truncated disc, hot inner flow, plus standard conical synchrotron jet model to explore what this predicts for the X-ray emission mechanism as a function of mass accretion rate. We find that the change in X-ray spectral index can be quantitatively (not just qualitatively) explained by the seed photon switch in the hot flow i.e. this supports models where the X-rays are always produced by the hot flow. By contrast, standard conical jet models are as radiatively inefficient as the hot flow so there is no transition in X-ray production mechanism with $\dot{m}$. Including the effects of electron cooling allows the jet X-rays to drop more slowly with accretion rate and hence overtake the X-rays from the hot flow, however this produces a corresponding change in the radio-X-ray correlation, which is not observed. We argue that the unbroken radio-X-ray correlation down to quiescence rules out the jet transition model as an explanation for the trend in X-ray spectral index. Our favoured model is then a truncated disc with an inner, hot, radiatively inefficient flow which always dominates the hard X-rays, coupled to a conical synchrotron jet which produces the radio emission. However, even this has issues at low $\dot{m}$ as the low optical depth and high temperature of the flow means that the Compton spectrum is not well approximated by a power law. This shows the need for a more sophisticated model for the electron distribution in the hot flow.

preprint2013arXiv

The effect of advection at luminosities close to Eddington: The ULX in M31

The transient, ultra-luminous X-ray source CXOM31 J004253.1+411422 in the Andromeda galaxy is most likely a 10 solar mass black hole, with super-Eddington luminosity at its peak. The XMM-Newton spectra taken during the decline then trace luminosities of 0.86-0.27 L_Edd. These spectra are all dominated by a hot disc component, which roughly follows a constant inner radius track in luminosity and temperature as the source declines. At the highest luminosity the disc structure should change due to advection of radiation through the disc. This advected flux can be partly released at lower radii thus modifying the spectral shape. To study the effect of advection at luminosities close to Eddington we employ a fully relativistic slim disc model, SLIMBH, that includes advective cooling and full radiative transfer through the photosphere based on TLUSTY. The model also incorporates relativistic photon ray-tracing from the proper location of the disc photosphere rather than the mid-plane as the slim disc is no longer geometrically thin. We find that these new models differ only slightly from the non-advective (standard) BHSPEC models even at the highest luminosities considered here. While both discs can fit the highest luminosity data, neither is a very good fit to the *lower* luminosities. This could indicate a missing physical process that acts in low luminosity discs and subsides as the disc luminosity approaches the Eddington limit.

preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

preprint2013arXiv

What Powers the Most Relativistic Jets? I: BL Lacs

The dramatic relativistic jets pointing directly at us in BL Lac objects can be well modelled by bulk motion beaming of synchrotron self-Compton emission powered by a low Eddington fraction accretion flow. Nearly 500 of these AGN are seen in the 2nd Fermi Large Area Telescope catalogue of AGN. We combine the jet models which describe individual spectra with the expected jet parameter scalings with mass and mass accretion rate to predict the expected number of Fermi detected sources given the number densities of AGN from cosmological simulations. We select only sources with Eddington scaled mass accretion rate $< 0.01$ (i.e. radiatively inefficient flows), and include cooling, orientation effects, and the effects of absorption from pair production on the extragalactic IR background. These models overpredict the number of Fermi detected BL Lacs by a factor of 1000! This clearly shows that one of the underlying assumptions is incorrect, almost certainly that jets do not scale simply with mass and accretion rate. The most plausible additional parameter which can affect the region producing the Fermi emission is black hole spin. We can reproduce the observed numbers of BL Lacs if such relativistic jets are only produced by the highest spin ($a_*>0.8$) black holes, in agreement with the longstanding spin-jet paradigm. This also requires that high spins are intrinsically rare, as predicted by the cosmological simulations for growing black hole mass via chaotic (randomly aligned) accretion episodes, where only the most massive black holes have high spin due to black hole-black hole mergers.

preprint2012arXiv

A combined Optical and X-ray Spectra Study for Type 1 AGN. III. Broadband SED Properties

In this third paper in a series of three, we present a detailed study of the AGN broadband SED based on a nearby unobscured Type 1 AGN sample. We perform a systematic cross-correlation study of the following key parameters: $Γ_{2-10keV}$, $L_{2-10keV}$, $L_{bol}$, $L_{bol}/L_{Edd}$, $κ_{2-10keV}$, $κ_{5100A}$, FWHM$_{Hβ}$, M$_{BH}$, $α_{ox}$, $α_{X}$ and $α_{UV}$, and identify various strong correlations among these parameters. The principal component analysis (PCA) is performed on the correlation matrix of the above parameters, which shows that the three physical parameters, i.e. black hole mass, mass accretion rate and Eddington ratio, drive the majority of the correlations. This is consistent with PCA results found from previous optical spectral studies. We produce various mean SEDs classified by each of the key parameters. Most parameters, except L$_{bol}$, show similar systematic changes in the mean SEDs such that the temperature at which the disc peaks is correlated with the ratio of power in the disc versus the Comptonised components and the hard X-ray spectral index. This underlying change in SED shape shows that AGN do exhibit intrinsically different spectral states. This is superficially similar to the SED differences in BHB seen as $λ_{Edd}$ increases, but the analogy does not hold in detail. Only objects with the highest $λ_{Edd}$ appear to correspond to a BHB spectral state (the disc dominated high/soft state). The AGN with typical mass accretion rates have spectra which do not match well with any state observed in BHB. We speculate that this could be due to the presence of a powerful UV line driven disc wind, which complicates simple mass scaling between stellar and supermassive black holes.

preprint2012arXiv

Bright radio emission from an ultraluminous stellar-mass microquasar in M31

A subset of ultraluminous X-ray sources (those with luminosities < 10^40 erg/s) are thought to be powered by the accretion of gas onto black holes with masses of ~5-20 M_solar, probably via an accretion disc. The X-ray and radio emission are coupled in such Galactic sources, with the radio emission originating in a relativistic jet thought to be launched from the innermost regions near the black hole, with the most powerful emission occurring when the rate of infalling matter approaches a theoretical maximum (the Eddington limit). Only four such maximal sources are known in the Milky Way, and the absorption of soft X-rays in the interstellar medium precludes determining the causal sequence of events that leads to the ejection of the jet. Here we report radio and X-ray observations of a bright new X-ray source whose peak luminosity can exceed 10^39 erg/s in the nearby galaxy, M31. The radio luminosity is extremely high and shows variability on a timescale of tens of minutes, arguing that the source is highly compact and powered by accretion close to the Eddington limit onto a stellar mass black hole. Continued radio and X-ray monitoring of such sources should reveal the causal relationship between the accretion flow and the powerful jet emission.

preprint2012arXiv

The ASTRO-H X-ray Observatory

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe via a suite of four instruments, covering a very wide energy range, from 0.3 keV to 600 keV. These instruments include a high-resolution, high-throughput spectrometer sensitive over 0.3-2 keV with high spectral resolution of Delta E < 7 eV, enabled by a micro-calorimeter array located in the focal plane of thin-foil X-ray optics; hard X-ray imaging spectrometers covering 5-80 keV, located in the focal plane of multilayer-coated, focusing hard X-ray mirrors; a wide-field imaging spectrometer sensitive over 0.4-12 keV, with an X-ray CCD camera in the focal plane of a soft X-ray telescope; and a non-focusing Compton-camera type soft gamma-ray detector, sensitive in the 40-600 keV band. The simultaneous broad bandpass, coupled with high spectral resolution, will enable the pursuit of a wide variety of important science themes.

preprint2012arXiv

The effect of frame dragging on the iron K alpha line in X-ray binaries

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

preprint2012arXiv

The truncated disk from Suzaku data of GX 339-4 in the extreme very high state

We report on the geometry of accretion disk and high energy coronae in the strong Comptonization state (the very high/steep power law/hard intermediate state) based on a Suzaku observation of the famous Galactic black hole GX 339-4. These data were taken just before the peak of the 2006-2007 outburst, and the average X-ray luminosity in the 0.7-200 keV band is estimated to be 2.9E38 erg/s for a distance of 8 kpc. We fit the spectrum with both simple (independent disk and corona) and sophisticated (energetically coupled disk and corona) models, but all fits imply that the underlying optically thick disk is truncated significantly before the innermost stable circular orbit around the black hole. We show this directly by a comparison with similarly broadband data from a disk dominated spectrum at almost the same luminosity observed by XMM-Newton and RXTE 3 days after the Suzaku observation. During the Suzaku observation, the QPO frequency changes from 4.3 Hz to 5.5 Hz, while the spectrum softens. The energetically coupled model gives a corresponding 5+/- 8 % decrease in derived inner radius of the disk. While this is not significant, it is consistent with the predicted change in QPO frequency from Lense-Thirring precession of the hot flow interior to the disk and/or a deformation mode of this flow, as a higher QPO frequency implies a smaller size scale for the corona. This is consistent with the truncated disk extending further inwards towards the black hole.

preprint2011arXiv

Intrinsic disc emission and the Soft X-ray Excess in AGN

(Abridged) Narrow Line Seyfert 1 (NLS1) galaxies have low mass black holes and mass accretion rates close to (or exceeding) Eddington, so a standard blackbody accretion disc should peak in the EUV. However, the lack of true absorption opacity in the disc means that the emission is better approximated by a colour temperature corrected blackbody, and this colour temperature correction is large enough ($\sim 2.4$) that the bare disc emission from a zero spin black hole can extend into the soft X-ray bandpass. Part of the soft X-ray excess seen in these objects must be intrinsic emission from the disc unless the vertical structure is very different to that predicted. However, the soft excess is much broader than predicted by a bare disc spectrum, indicating some Compton upscattering by cool, optically thick material. We associate this with the disc itself, so it must ultimately be powered by mass accretion. We build an energetically self consistent model assuming that the emission thermalises at large radii, but that at smaller radii the gravitational energy is split between powering optically thick Comptonised disc emission (forming the soft X-ray excess) and an optically thin corona above the disc (forming the tail to higher energies). We show examples of this model fit to the extreme NLS1 REJ1034+396, and to the much lower Eddington fraction Broad Line Seyfert 1 PG1048+231. We use these to guide our fits and interpretations of three template spectra made from co-adding multiple sources to track out a sequence of AGN spectra as a function of $L/L_{Edd}$. The new model is publically available within the {\sc xspec} spectral fitting package.

preprint2011arXiv

Modelling the high mass accretion rate spectra of GX 339-4: Black hole spin from reflection?

We extract all the XMM-Newton EPIC pn burst mode spectra of GX 339-4, together with simultaneous/contemporaneous RXTE data. These include three disc dominated and two soft intermediate spectra, and the combination of broad bandpass/moderate spectral resolution gives some of the best data on these bright soft states in black hole binaries. The disc dominated spectra span a factor three in luminosity, and all show that the disc emission is broader than the simplest multicolour disc model. This is consistent with the expected relativistic smearing and changing colour temperature correction produced by atomic features in the newest disc models. However, these models do not match the data at the 5 per cent level as the predicted atomic features are not present in the data, perhaps indicating that irradiation is important even when the high energy tail is weak. Whatever the reason, this means that the data have smaller errors than the best physical disc models, forcing use of more phenomenological models for the disc emission. We use these for the soft intermediate state data, where previous analysis using a simple disc continuum found an extremely broad residual, identified as the red wing of the iron line from reflection around a highly spinning black hole. However, the iron line energy is close to where the disc and tail have equal fluxes, so using a broader disc continuum changes the residual 'iron line' profile dramatically. With a broader disc continuum model, the inferred line is formed outside of 30 ${\rm{R_g}}$, so cannot constrain black hole spin. We caution that a robust determination of black hole spin from the iron line profile is very difficult where the disc makes a significant contribution at the iron line energy i.e. in most bright black hole states.

preprint2011arXiv

Modelling variability in black hole binaries: linking simulations to observations

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

preprint2011arXiv

Searching for the trigger of the AGN QPO: 8 years of RE J1034+396

RE J1034+396 is one of the most extreme Narrow-line Seyfert 1s detected thus far, showing the only quasi-periodic oscillation (QPO) reliably detected in an Active Galactic Nucleus (AGN). Comparison with similar spectral and timing properties observed in the black hole X-ray binary (BHB) GRS 1915+105 suggests that RE J1034+396 is a super-Eddington accretor. A more complete understanding of the behaviour of RE J1034+396 can therefore lead to a unification of the accretion physics between such extreme AGN and super-Eddington BHBs. Here we report on our latest XMM-Newton observations of RE J1034+396, which no longer show the QPO, indicating that this source shows a non-stationary power-spectrum. We use spectral and temporal analysis across all five XMM-Newton observations of the source to probe the evolution of the object. The combination of the shape of the fractional variability with energy and the inferred velocity of absorbing material in the line-of-sight rules out an absorption-only method of creating the QPO. Instead the periodically changing absorption may be produced by the QPO causing a change in ionization state. We extend our analysis by including the covariance spectra which give much better signal to noise than an rms spectrum. These reveal a new aspect of the QPO, which is that there is also a small contribution from a soft component which is hotter than the soft excess seen in the mean spectrum. Folding the lightcurve on the QPO period shows that this component lags behind the hard X-rays. If this is due to re-processing then the lag corresponding to a light travel time across 30 R_g. Some of the remaining observations have similar energy spectra and covariance spectra, but none of them show a significant QPO, so we conclude that none of these features are the trigger for the appearance of the QPO in this object.

preprint2011arXiv

Spectral and Timing Studies of Cyg X-1 in the Low/Hard State with Suzaku

From 2005 to 2009, 25 observations of Cyg X-1 were performed with Suzaku, achieving a total exposure of 446 ks. In all observations, the source was found in the low/hard state, while the 1.5--12.0 keV count rate of the All-Sky Monitor onboard RXTE varied by a factor of ~ 3. In each observation, the 10--60 keV HXD-PIN spectrum and the 60--400 keV HXD-GSO spectrum were fitted successfully by a thermal Comptonization model plus reflection by a thick neutral material. As the soft X-ray intensity increased, the Compton y-parameter was found to decrease from 1.0 to 0.6, while the solid angle of reflection to increase by ~ 30%. Also conducted was timing analysis over a frequency range of 0.001--10 Hz. As the source became brighter in soft X-rays, the characteristic frequency of hard X-ray variation increased from 0.03 to 0.3 Hz, while the fractional hard X-ray variation integrated over 0.001--0.01 Hz decreased by a factor of ~ 5. The signals in the 60--200 keV band were generally found to vary on shorter time scales than those in the 10--60 keV band. These spectral and timing results can be consistently interpreted by presuming that increases in the mass accretion rate cause the Comptonizing hot corona to shrink, while the optically-thick disk to intrude deeper therein.

preprint2011arXiv

The missing link: a low mass X-ray binary in M31 seen as an ultraluminous X-ray source

A new, transient ultraluminous X-ray source (ULX) was recently discovered by Chandra in M31 with a luminosity at ~ 5 x 10^39 erg/s. Here we analyse a series of five subsequent XMM-Newton observations. The steady decline in X-ray luminosity over 1.5 months gives an observed e-fold timescale of ~40 days; similar to the decay timescales seen in soft X-ray transients in our own Galaxy. This supports the interpretation of this ULX as a stellar mass black hole in a low-mass X-ray binary (LMXB), accreting at super Eddington rates. This is further supported by the lack of detection of an O/B star in quiescence and the XMM-Newton spectral data being dominated by a disc-like component rather than the power-law expected from a sub-Eddington intermediate-mass black hole. These data give the best sequence of high Eddington fraction spectra ever assembled due to the combination of low absorption column to M31 and well calibrated bandpass down to 0.3 keV of XMM-Newton. The spectra can be roughly described by our best current disc model, BHSPEC, assuming a 10 M_solar black hole with best fit spin of ~0.4. However, the data are better described by a two component model, where the disc emission is significantly affected by advection, and with an additional low temperature Comptonisation component at high energies. This could simply indicate the limitations of current disc models, though the energy-dependent variability also weakly supports a two component interpretation of the data. Irrespective of the detailed interpretation of the spectral properties, these data support the presence of accretion onto a stellar mass black hole in a LMXB accreting in the Eddington-regime. This allows an unambiguous connection of this object, and, by extension, similar low luminosity ULXs, to `standard' X-ray binaries.

preprint2011arXiv

Unlocking the nature of ultraluminous X-ray sources using their X-ray spectra

We explore the nature of ultraluminous X-ray sources through detailed investigations of their spectral shape using some of the highest quality data available in the XMM-Newton public archives. Phenomenological models allow us to characterise their spectra, while more 'physically-motivated' models enable us to explore the physical processes underlying these characteristics. These physical models imply the presence of extreme (probably super-Eddington) accretion on to stellar mass black holes.

preprint2010arXiv

A physical interpretation of the variability power spectral components in accreting neutron stars

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

preprint2010arXiv

A re--analysis of the iron line in the XMM-Newton data from the low/hard state in GX339--4

The detection of an extremely broad iron line in XMM-Newton MOS data from the low/hard state of the black hole binary GX339-4 is the only piece of evidence which unambiguously conflicts with the otherwise extremely successful truncated disc interpretation of this state. However, it also conflicts with some aspect of observational data for all other alternative geometries of the low/hard state, including jet models, making it very difficult to understand. We re-analyse these data and show that they are strongly affected by pileup even with extensive centroid removal as the source is ~200x brighter than the recommended maximum countrate. Instead, we extract the simultaneous PN timing mode data which should not be affected by pileup. These show a line which is significantly narrower than in the MOS data. Thus these data are easily consistent with a truncated disc, and indeed, strongly support such an interpretation.

preprint2010arXiv

A strong and broad iron line in the XMM-Newton spectrum of the new X-ray transient and black-hole candidate XTE J1652-453

We observed the new X-ray transient and black-hole candidate XTE J1652-453 simultaneously with XMM-Newton and the Rossi X-ray Timing Explorer (RXTE). The observation was done during the decay of the 2009 outburst, when XTE J1652-453 was in the hard-intermediate state. The spectrum shows a strong and broad iron emission line with an equivalent width of ~ 450 eV. The profile is consistent with that of a line being produced by reflection off the accretion disk, broadened by relativistic effects close to the black hole. The best-fitting inner radius of the accretion disk is ~ 4 gravitational radii. Assuming that the accretion disk is truncated at the radius of the innermost stable circular orbit, the black hole in XTE J1652-453 has a spin parameter of ~ 0.5. The power spectrum of the RXTE observation has an additional variability component above 50 Hz, which is typical for the hard-intermediate state. No coherent quasi-periodic oscillations at low frequency are apparent in the power spectrum, which may imply that we view the system at a rather low inclination angle.

preprint2010arXiv

A study of spectral and timing properties of Cyg X-1 based on a large sample of pointed RXTE observations

We study a large sample of RXTE PCA/HEXTE observations of Cyg X-1. We characterize the spectra by soft and hard X-ray colours (which define the spectral states), and fit them with a physical model of hybrid, thermal/non-thermal Comptonization. We then fit the power spectra by a sum of Lorentzians. We show the resulting correlations between the spectral colour (defining the hard, intermediate and soft spectral states), the Comptonization amplification factor, the fractional Compton reflection strength, and the peak frequency of the Lorentzians. We also calculate the fractional variability (rms) as a function of photon energy. We fit the obtained rms dependencies by physical models of varying accretion rate, a varying soft-photon input from an outer disc to a hot inner plasma, and a varying dissipation rate in a hot corona above a disc.

preprint2010arXiv

Limits on spin determination from disc spectral fitting in GX 339-4

We attempt to constrain the black hole spin in GX 339-4 from spectral fitting of disc dominated data using RXTE spectra from the three most recent outbursts. We use the best current models for the disc emission, including full radiative transfer through the photosphere rather than assuming that the intrinsic emission from each radius has a (colour temperature corrected) blackbody spectrum. The results strongly depend on the poorly known binary system parameters, but we find a strict upper limit of a<0.9 for any distance greater than 6 kpc, assuming that the orbital inclination is the same as that of the inner disc. By contrast, the higher spin of 0.935 +/- 0.01 (statistical) +/-0.01 (systematic) claimed from fitting the iron line profile in this object requires that the inner disc is misaligned by over 20 degrees from the orbital inclination. While some of these datasets are distorted by instrumental pileup, the same spin/inclination constraints are derived from data which are not piled up, so there is a real conflict between the two techniques to measure spin. We argue that the disc spectral fits are more likely to be robust hence that there are still issues to be understood in the iron line profile.

preprint2010arXiv

Observational characteristics of accretion onto black holes

These notes resulted from a series of lectures at the IAC winter school. They are designed to help students, especially those just starting in subject, to get hold of the fundamental tools used to study accretion powered sources. As such, the references give a place to start reading, rather than representing a complete survey of work done in the field. I outline Compton scattering and blackbody radiation as the two predominant radiation mechanisms for accreting black holes, producing the hard X-ray tail and disc spectral components, respectively. The interaction of this radiation with matter can result in photo-electric absorption and/or reflection. While the basic processes can be found in any textbook, here I focus on how these can be used as a toolkit to interpret the spectra and variability of black hole binaries (hereafter BHB) and Active Galactic Nuclei (AGN). I also discuss how to use these to physically interpret real data using the publicly available XSPEC spectral fitting package (Arnaud et al 1996), and how this has led to current models (and controversies) of the accretion flow in both BHB and AGN.

preprint2010arXiv

Quasi-periodic oscillations under wavelet microscope: the application of Matching Pursuit algorithm

We zoom in on the internal structure of the low-frequency quasi-periodic oscillation (LF QPO) often observed in black hole binary systems to investigate the physical nature of the lack of coherence in this feature. We show the limitations of standard Fourier power spectral analysis for following the evolution of the QPO with time and instead use wavelet analysis and a new time-frequency technique -Matching Pursuit algorithm- to maximise the resolution with which we can follow the QPO behaviour. We use the LF QPO seen in a very high state of XTE J1550-564 to illustrate these techniques and show that the best description of the QPO is that it is composed of multiple independent oscillations with a distribution of lifetimes but with constant frequency over this duration. This rules out models where there is continual frequency modulation, such as multiple blobs spiralling inwards. Instead it favours models where the QPO is excited by random turbulence in the flow.

preprint2010arXiv

Testing Accretion Disk Structure with Suzaku data of LMC X-3

The Suzaku observation of LMC X-3 gives the best data to date on the shape of the accretion disk spectrum. This is due to the combination of very low absorbing column density along this line of sight which allows the shape of the disk emisison to be constrained at low energies by the CCD's, while the tail can be simultaneously determined up to 30 keV by the high energy detectors. These data clearly demonstrate that the observed disk spectrum is broader than a simple 'sum of blackbodies', and relativistic smearing of the emission is strongly required. However, the intrinsic emission should be more complex than a (color-corrected) sum of blackbodies as it should also contain photo-electric absorption edges from the partially ionised disk photosphere. These are broadened by the relativistic smearing, but the models predict ~ 3-5 per cent deviations for 1/3- 1 solar abundance around the edge energies, significantly stronger than observed. This indicate that the models need to include more physical processes such as self-irradiation, bound-bound (line) absorption, and/or emission from recombination continuua and/or lines. Alternatively, if none of these match the data, it may instead require that the accretion disk density and/or emissivity profile with height is different to that assumed. Thus these data demonstrate the feasibility of observational tests of our fundamental understanding of the vertical structure of accretion disks.

preprint2010arXiv

The ASTRO-H Mission

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe by performing high-resolution, high-throughput spectroscopy with moderate angular resolution. ASTRO-H covers very wide energy range from 0.3 keV to 600 keV. ASTRO-H allows a combination of wide band X-ray spectroscopy (5-80 keV) provided by multilayer coating, focusing hard X-ray mirrors and hard X-ray imaging detectors, and high energy-resolution soft X-ray spectroscopy (0.3-12 keV) provided by thin-foil X-ray optics and a micro-calorimeter array. The mission will also carry an X-ray CCD camera as a focal plane detector for a soft X-ray telescope (0.4-12 keV) and a non-focusing soft gamma-ray detector (40-600 keV) . The micro-calorimeter system is developed by an international collaboration led by ISAS/JAXA and NASA. The simultaneous broad bandpass, coupled with high spectral resolution of Delta E ~7 eV provided by the micro-calorimeter will enable a wide variety of important science themes to be pursued.

preprint2009arXiv

An additional soft X-ray component in the dim low/hard state of black hole binaries

We test the truncated disc models using multiwavelength (optical/UV/X-ray) data from the 2005 hard state outburst of the black hole SWIFT J1753.5-0127. This system is both fairly bright and has fairly low interstellar absorption, so gives one of the best datasets to study the weak, cool disc emission in this state. We fit these data using models of an X-ray illuminated disc to constrain the inner disc radius throughout the outburst. Close to the peak, the observed soft X-ray component is consistent with being produced by the inner disc, with its intrinsic emission enhanced in temperature and luminosity by reprocessing of hard X-ray illumination in an overlap region between the disc and corona. This disc emission provides the seed photons for Compton scattering to produce the hard X-ray spectrum, and these hard X-rays also illuminate the outer disc, producing the optical emission by reprocessing. However, the situation is very different as the outburst declines. The optical is probably cyclo-synchrotron radiation, self-generated by the flow, rather than tracing the outer disc. Similarly, limits from reprocessing make it unlikely that the soft X-rays are directly tracing the inner disc radius. This is seen more clearly in a similarly dim low/hard state spectrum from XTE J1118+480. The very small emitting area implied by the relatively high temperature soft X-ray component is completely inconsistent with the much larger, cooler, UV component which is well fit by a truncated disc. We speculate on the origin of this component, but its existence as a clearly separate spectral component from the truncated disc in XTE J1118+480 shows that it does not simply trace the inner disc radius, so cannot constrain the truncated disc models.

preprint2009arXiv

The Ultraluminous State

(Abridged) We revisit the question of the nature of ULXs through a detailed investigation of their spectral shape, using the highest quality X-ray data available in the XMM-Newton public archives. We confirm that simple spectral models commonly used for the analysis and interpretation of ULXs (power-law continuum and multi-colour disc blackbody models) are inadequate in the face of such high quality data. Instead we find two near ubiquitous features in the spectrum: a soft excess and a roll-over in the spectrum at energies above 3keV. We investigate a range of more physical models to describe these data. We find that disc plus Comptonised corona models fit the data well, but the derived corona is cool, and optically thick (tau ~ 5-30). We argue that these observed disc temperatures are not a good indicator of the black hole mass as the powerful, optically thick corona drains energy from the inner disc, and obscures it. We estimate the intrinsic (corona-less) disc temperature, and demonstrate that in most cases it lies in the regime of stellar mass black holes. These objects have spectra which range from those similar to the highest mass accretion rate states in Galactic binaries, to those which clearly have two peaks, one at energies below 1 keV (from the outer, unComptonised disc) and one above 3 keV (from the Comptonised, inner disc). However, a few ULXs have a significantly cooler corrected disc temperature; we suggest that these are the most extreme stellar mass black hole accretors, in which a massive wind completely envelopes the inner disc regions, creating a cool photosphere. We conclude that ULXs provide us with an observational template for the transition between Eddington and super-Eddington accretion flows, with the latter occupying a new ultraluminous accretion state.

preprint2008arXiv

RE J1034+396: The origin of the soft X-ray excess and QPO

The X-ray quasi-periodic oscillation (QPO) seen in RE J1034+396 is so far unique amongst AGN. Here we look at the another unique feature of RE J1034+396, namely its huge soft X-ray excess, to see if this is related in any way to the detection of the QPO. We show that all potential models considered for the soft energy excess can fit the 0.3-10 keV X-ray spectrum, but that the energy dependence of the rapid variability (which is dominated by the QPO) strongly supports a spectral decomposition where the soft excess is from low temperature Comptonization of the disc emission and remains mostly constant, while the rapid variability is produced by the power law tail changing in normalization. The presence of the QPO in the tail rather than in the disc is a common feature in black hole binaries, but low temperature Comptonization of the disc spectrum is not generally seen in these systems. The main exception to this is GRS 1915+105, the only black hole binary which routinely shows super-Eddington luminosities. We speculate that super-Eddington accretion rates lead to a change in disc structure, and that this also triggers the X-ray QPO.

preprint2008arXiv

Suzaku Results on Cygnus X-1 in the Low/Hard State

The black-hole binary Cygnus X-1 was observed for 17 ks with the Suzaku X-ray observatory in 2005 October, while it was in a low/hard state with a 0.7-300 keV luminosity of 4.6 x 10^37 erg/s. The XIS and HXD spectra, spanning 0.7-400 keV, were reproduced successfully incorporating a cool accretion disk and a hot Comptonizing corona. The corona is characterized by an electron temperature of ~100 keV, and two optical depths of ~0.4 and ~1.5 which account for the harder and softer continua, respectively. The disk has the innermost temperature of ~0.2 keV, and is thought to protrude half way into the corona. The disk not only provides seed photons to the Compton cloud, but also produces a soft spectral excess, a mild reflection hump, and a weakly broadened iron line. A comparison with the Suzaku data on GRO J1655-40 reveals several interesting spectral differences, which can mostly be attributed to inclination effects assuming that the disk has a flat geometry while the corona is grossly spherical. An intensity-sorted spectroscopy indicates that the continuum becomes less Comptonized when the source flares up on times scales of 1-200 s, while the underlying disk remains unchanged.

preprint2007arXiv

High Energy spectra of Seyferts and Unification schemes

The Unified Model of AGN predicts the sole difference between Seyfert 1 and Seyfert 2 nuclei is the viewing angle with respect to an obscuring structure around the nucleus. High energy photons above 20 keV are not affected by this absorption if the column is Compton thin, so their 30--100 keV spectra should be the same. However, the observed spectra at high energies appear to show a systematic difference, with Seyfert 1's having $Γ\sim $2.1 whereas Seyfert 2's are harder with $Γ\sim$ 1.9. We estimate the mass and accretion rate of Seyferts detected in these high energy samples and show that they span a wide range in $L/L_{Edd}$. Both black hole binary systems and AGN show a correlation between spectral softness and Eddington fraction, so these samples are probably heterogeneous, spanning a range of intrinsic spectral indices which are hidden in individual objects by poor signal-to-noise. However, the mean Eddington fraction for the Seyfert 1's is higher than for the Seyfert 2's, so the samples are consistent with this being the origin of the softer spectra seen in Seyfert 1's. We stress that high energy spectra alone are not necessarily a clean test of Unification schemes, but that the intrinsic nuclear properties should also change with $L/L_{Edd}$

preprint2006arXiv

Suzaku discovery of iron absorption lines in outburst spectra of the X-ray transient 4U 1630-472

We present the results of six Suzaku observations of the recurrent black hole transient 4U1630-472 during its decline from outburst from February 8 to March 23 in 2006. All observations show the typical high/soft state spectral shape in the 2-50keV band, roughly described by an optically thick disk spectrum in the soft energy band plus a weak power-law tail that becomes dominant only above \~20keV. The disk temperature decreases from 1.4keV to 1.2keV as the flux decreases by a factor 2, consistent with a constant radius as expected for disk-dominated spectra. All the observations reveal significant absorption lines from highly ionized (H-like and He-like) iron Ka at 7.0keV and 6.7keV. The brightest datasets also show significant but weaker absorption structures between 7.8keV and 8.2keV, which we identify as a blend of iron Kb and nickel Ka absorption lines. The energies of these absorption lines suggest a blue shift with an outflow velocity of ~1000km/s.. The H-like iron Ka equivalent width remains approximately constant at ~30 eV over all the observations, while that of the He-like Ka line increases from 7 eV to 20 eV. Thus the ionization state of the material decreases, as expected from the decline in flux. We fit the profiles with Voigt functions (curve of growth) to derive absorbing columns, then use these together with detailed photo-ionization calculations to derive physical parameters of the absorbing material. The data then constrain the velocity dispersion of the absorber to 200-2000 km/s, and the size of the plasma as ~1E10 cm assuming a source distance of 10 kpc.

preprint2004arXiv

Extreme Gravitational Lensing near Rotating Black Holes

We describe a new approach to calculating photon trajectories and gravitational lensing effects in the strong gravitational field of the Kerr black hole. These techniques are applied to explore both the imaging and spectral properties of photons that perform multiple orbits of the central mass before escaping to infinity. Viewed at large inclinations, these higher order photons contribute $\sim 20 %$ of the total luminosity of the system for a Schwarzschild hole, whilst for an extreme Kerr black hole this fraction rises to $\sim 60 %$. In more realistic models these photons will be re-absorbed by the disc at large distances from the hole, but this returning radiation could provide a physical mechanism to resolve the discrepancy between the predicted and observed optical/UV colours in AGN. Conversely, at low inclinations, higher order images re-intercept the disc plane close to the black hole, so need not be absorbed by the disc if this is within the plunging region. These photons form a bright ring carrying approximately 10% of the total disc luminosity for a Schwarzchild black hole. The spatial separation between the inner edge of the disc and the ring is similar to the size of the event horizon. This is resolvable for supermassive black holes with proposed X-ray interferometery missions such as MAXIM, so has the potential to provide an observational test of strong field gravity.

preprint2004arXiv

Observing braneworld black holes

Spacetime in the vicinity of an event horizon can be probed using observations which explore the dynamics of the accretion disc. Many high energy theories of gravity lead to modifications of the near horizon regime, potentially providing a testing ground for these theories. In this paper, we explore the impact of braneworld gravity on this region by formulating a method of deriving the general behaviour of the as yet unknown braneworld black hole solution. We use simple bounds to constrain the solution close to the horizon.

preprint2004arXiv

What is the origin of the soft excess in AGN?

We investigate the nature of the soft excess below 1 keV observed in AGN. We use the XMM-Newton data of the low redshift, optically bright quasar, PG 1211+143, and we compare it with the Narrow Line Seyfert 1 galaxy, 1H 0707-495, which has one of the strongest soft excesses seen. We test various ideas for the origin of the soft X-ray excess, including a separate spectral component (for example a low temperature Comptonized component), a reflection-dominated model, or a complex absorption model. All three can give good fits to the data, and $χ^2$ fitting criteria are not sufficient to discriminate among them. Instead, we favor the complex absorption model on the grounds that it is the most physically plausible.

preprint1997arXiv

Relativistically smeared X-ray reprocessed components in the GINGA spectra of GS 2023+338

We present results of spectral analysis of Ginga data obtained during the decline phase after the 1989 outburst of GS 2023+338 (V404 Cyg). Our analysis includes detailed modelling of the effects of X-ray reflection/reprocessing. We have found that (1) the contribution of the reprocessed component (both continuum and line) corresponds to the solid angle of the reprocessor as seen from the X-ray source of Omega \approx (0.4-0.5) 2π, (2) the reprocessed component (both line and continuum) is broadened (``smeared'') by kinematic and relativistic effects, as expected from the accretion disk reflection. We discuss the constraints these results give on various possible system geometries.

preprint1997arXiv

Simultaneous X-ray and gamma-ray observations of Cyg X-1 in the hard state by Ginga and OSSE

We present four X-ray/gamma-ray spectra of Cyg X-1 observed in the hard ('low') state simultaneously by Ginga and GRO/OSSE on 1991 July 6. The 3-30 keV Ginga spectra are well represented by power laws with an energy spectral index of alpha~0.6 and a Compton reflection component including a fluorescent Fe K-alpha corresponding to the solid angle of the reflector of ~0.3 times 2 Pi. The overall Ginga/OSSE spectra can be modelled by repeated Compton scattering in a mildly-relativistic, tau ~1, plasma. However, the high-energy cutoff is steeper than that due to single-temperature thermal Comptonisation. It can be described by a superposition of dominant optically-thin, thermal emission at kT~140 keV and a Wien-like component from an optically-thick plasma at kT~50 keV. The X-ray spectra do not show the presence of an anisotropy break required if thermal Compton scattering takes place in a corona above a cold disc. Also, the flat spectral index shows that the plasma is soft-photon starved, i.e., the luminosity in incident soft X-ray seed photons is very much less than that in the hard X-rays. Furthermore, the observed solid angle of the reflector is significantly less than 2 pi. These facts taken together strongly rule out a disc-corona geometry. Rather, the observed spectra are consistent with a geometry in which the cold accretion disc (which both supplies the seed soft X-rays and reflects hard X-rays) only exists at large radii, while the Comptonising hot plasma is located in an inner region with no cold disc. This hot plasma consists of either pure electron-positron pairs if the source size is ~5 Schwarzschild radii or it contains also protons if the size is larger.