Source author record

Matthew Middleton

Matthew Middleton 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

16works
6topics
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

16 published item(s)

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

Eclipse Mapping of EXO 0748-676: Evidence for a Massive Neutron Star

Determining the maximum possible neutron star (NS) mass places limits on the equation of state (EoS) of ultra-dense matter. The mass of NSs in low mass X-ray binaries can be determined from the binary mass function, providing independent constraints are placed on both the binary inclination and mass ratio. In eclipsing systems, they relate via the totality duration. EXO 0748-676 is an eclipsing NS low mass X-ray binary with a binary mass function estimated using stellar emission lines from the irradiated face of the companion. The NS mass is thus known as a function of mass ratio. Here we model the X-ray eclipses in several energy bands, utilising archival XMM-Newton data. We find a narrow region of absorbing material surrounding the companion star is required to explain the energy-dependent eclipses. Therefore, we suggest the companion may be experiencing ablation of its outer layers and that the system could transition into a redback millisecond pulsar. Our fit returns a mass ratio of $q=0.222^{+0.07}_{-0.08}$ and an inclination $i = 76.5 \pm^{1.4}_{1.1}$. Combining these with the previously measured radial velocity of $410 \pm 5$ km/s, derived from Doppler mapping analysis of H$_α$ emission during quiescence, returns a NS mass of $\sim 2 M_\odot$ even if the line originates as far from the NS as physically possible, favouring hard EoS. The inferred mass increases for a more realistic emission point. However, a $\sim 1.4 M_\odot$ canonical NS mass is possible when considering radial velocity values derived from other emission lines observed both during outburst and quiescence.

preprint2022arXiv

Phase-resolved spectroscopy of a quasi-periodic oscillation in the black hole X-ray binary GRS 1915+105 with NICER and NuSTAR

Quasi-periodic oscillations (QPOs) are often present in the X-ray flux from accreting stellar-mass black holes (BHs). If they are due to relativistic (Lense-Thirring) precession of an inner accretion flow which is misaligned with the disc, the iron emission line caused by irradiation of the disc by the inner flow will rock systematically between red and blue shifted during each QPO cycle. Here we conduct phase-resolved spectroscopy of a $\sim2.2$ Hz type-C QPO from the BH X-ray binary GRS 1915+105, observed simultaneously with NICER and NuSTAR. We apply a tomographic model in order to constrain the QPO phase-dependent illumination profile of the disc. We detect the predicted QPO phase-dependent shifts of the iron line centroid energy, with our best fit featuring an asymmetric illumination profile ($>2σ$ confidence). The observed line energy shifts can alternatively be explained by the spiral density waves of the accretion-ejection instability model. However we additionally measure a significant ($>3σ$) modulation in reflection fraction, strongly favouring a geometric QPO origin. We infer that the disc is misaligned with previously observed jet ejections, which is consistent with the model of a truncated disc with an inner precessing hot flow. However our inferred disc inner radius is small ($r_\text{in}{\sim} 1.4 GM/c^2$). For this disc inner radius, Lense-Thirring precession cannot reproduce the observed QPO frequency. In fact, this disc inner radius is incompatible with the predictions of all well-studied QPO models in the literature.

preprint2022arXiv

X-ray Eclipse Mapping Constrains the Binary Inclination and Mass Ratio of Swift J1858.6-0814

X-ray eclipse mapping is a promising modelling technique, capable of constraining the mass and/or radius of neutron stars (NSs) or black holes (BHs) in eclipsing binaries and probing any structure surrounding the companion star. In eclipsing systems, the binary inclination, $i$, and mass ratio, $q$ relate via the duration of totality, $t_{e}$. The degeneracy between $i$ and $q$ can then be broken through detailed modelling of the eclipse profile. Here we model the eclipses of the NS low-mass X-ray binary Swift J1858.6$-$0814 utilising archival NICER observations taken while the source was in outburst. Analogous to EXO 0748$-$676, we find evidence for irradiation driven ablation of the companion's surface by requiring a layer of stellar material to surround the companion star in our modelling. This material layer extends $\sim 7000 - 14000$ km from the companion's surface and is likely the cause of the extended, energy-dependent and asymmetric ingress and egress that we observe. Our fits return an inclination of $i \sim 81^{\circ}$ and a mass ratio $q \sim 0.14$. Using Kepler's law to relate the mass and radius of the companion star via the orbital period ($\sim$ 21.3 hrs), we subsequently determine the companion to have a low mass in the range $0.183 M_{\odot} \leq M_{cs} \leq 0.372 M_{\odot}$ and a large radius in the range $1.02 R_{\odot} \leq R_{cs} \leq 1.29 R_{\odot}$. Our results, combined with future radial velocity amplitudes measured from stellar absorption/emission lines, can place precise constraints on the component masses in this system.

preprint2016arXiv

A broadband X-ray spectral study of the intermediate-mass black hole candidate M82 X-1 with NuSTAR, Chandra and Swift

M82 X-1 is one of the brightest ultraluminous X-ray sources (ULXs) known, which, assuming Eddington-limited accretion and other considerations, makes it one of the best intermediate-mass black hole (IMBH) candidates. However, the ULX may still be explained by super-Eddington accretion onto a stellar-remnant black hole. We present simultaneous NuSTAR, Chandra and Swift/XRT observations during the peak of a flaring episode with the aim of modeling the emission of M82 X-1 and yielding insights into its nature. We find that thin-accretion disk models all require accretion rates at or above the Eddington limit in order to reproduce the spectral shape, given a range of black hole masses and spins. Since at these high Eddington ratios the thin-disk model breaks down due to radial advection in the disk, we discard the results of the thin-disk models as unphysical. We find that the temperature profile as a function of disk radius ($T(r)\propto r^{-p}$) is significantly flatter ($p=0.55^{+ 0.07}_{- 0.04}$) than expected for a standard thin disk ($p=0.75$). A flatter profile is instead characteristic of a slim disk which is highly suggestive of super-Eddington accretion. Furthermore, radiation hydrodynamical simulations of super-Eddington accretion have shown that the predicted spectra of these systems are very similar to what we observe for M82 X-1. We therefore conclude that M82 X-1 is a super-Eddington accretor. Our mass estimates inferred from the inner disk radius imply a stellar-remnant black hole ($M_{\rm BH}=26^{+9}_{-6} M_{\odot}$) when assuming zero spin, or an IMBH ($M_{\rm BH}=125^{+45}_{-30} M_{\odot}$) when assuming maximal spin.

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

Geometrical beaming of stellar mass ULXs

The presence or lack of eclipses in the X-ray lightcurves of ultraluminous X-ray sources (ULXs) can be directly linked to the accreting system geometry. In the case where the compact object is stellar mass and radiates isotropically, we should expect eclipses by a main-sequence to sub-giant secondary star on the recurrence timescale of hours to days. X-ray lightcurves are now available for large numbers of ULXs as a result of the latest XMM-Newton catalogue. We determine the amount of fractional variability that should be injected into an otherwise featureless lightcurve for a given set of system parameters as a result of eclipses and compare this to the available data. We find that the vast majority of sources for which the variability has been measured to be non-zero and for which available observations meet the criteria for eclipse searches, have fractional variabilities which are too low to derive from eclipses and so must be viewed such that theta =< acos(R*/a). This would require that the disc subtends a larger angle than that of the secondary star and is therefore consistent with a conical outflow formed from super-critical accretion rates and implies some level of geometrical beaming in ULXs.

preprint2016arXiv

The optical-UV emissivity of quasars: dependence on black hole mass and radio loudness

We analyzed a large sample of radio-loud and radio-quiet quasar spectra at redshift 1.0 < z < 1.2 to compare the inferred underlying quasar continuum slopes (after removal of the host galaxy contribution) with accretion disk models. The latter predict redder (decreasing) alpha_3000 continuum slopes (L_ν~ν^alpha at 3000Ang) with increasing black hole mass, bluer alpha_3000 with increasing luminosity at 3000Ang, and bluer alpha_3000 with increasing spin of the black hole, when all other parameters are held fixed. We find no clear evidence for any of these predictions in the data. In particular we find that: (i) alpha_3000 shows no significant dependence on black hole mass or luminosity. Dedicated Monte Carlo tests suggest that the substantial observational uncertainties in the black hole virial masses can effectively erase any intrinsic dependence of alpha_3000 on black hole mass, in line with some previous studies. (ii) The mean slope alpha_3000 of radio-loud sources, thought to be produced by rapidly spinning black holes, is comparable to, or even redder than, that of radio-quiet quasars. Indeed, although quasars appear to become more radio loud with decreasing luminosity, we still do not detect any significant dependence of alpha_3000 on radio loudness. The predicted mean alpha_3000 slopes tend to be bluer than in the data. Disk models with high inclinations and dust extinction tend to produce redder slopes closer to empirical estimates. Our mean alpha_3000 values are close to the ones independently inferred at z<0.5 suggesting weak evolution with redshift, at least for moderately luminous quasars.

preprint2014arXiv

A spectral-timing model for ULXs in the super-critical regime

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

preprint2014arXiv

Doppler disc tomography applied to low mass AGN spin

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

preprint2014arXiv

The low or retrograde spin of the first extragalactic microquasar: implications for Blandford-Znajek powering of jets

Transitions to high mass accretion rates in black hole X-ray binaries are associated with the ejection of powerful, relativistically-moving jets. The mechanism powering such events is thought to be linked to tapping of the angular momentum (spin) of the black hole, the rate of accretion through the disc or some combination of the two. We can attempt to discriminate between these possibilities by comparing proxies for jet power with spin estimates. Due to the small number of sources reaching Eddington rates and have therefore been suggested to act as 'standard candles', there has been much recent debate as to whether a significant correlation exists between jet power and spin. We perform continuum fitting to the high-quality, disc-dominated XMM-Newton spectra of the extragalactic microquasar discovered in M31. Assuming prograde spin, we find that, for sensible constraints the spin is always very low (a < 0.15 at 3-sigma). When combined with a proxy for jet power derived from the maximum 5 GHz radio luminosity during a bright flaring event, we find that the source sits well above the previously reported, rising correlation that would indicate that spin tapping is the dominant mechanism for powering the jets. The notable exceptions require the inclination to be improbably small or the jet to be very fast. We investigate whether this could be a by-product of selecting prograde-only spin, finding that the data statistically favour a substantially retrograde spin for the same constraints (a < -0.17 at 3-sigma). Although theoretically improbable, this remarkable finding could be confirmation that retrograde spin can power such jets via spin-tapping, as has been suggested for certain radio quasars. In either case this work demonstrates the value of studying local extragalactic microquasars as a means to better understand the physics of jet launching.

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.

preprint2012arXiv

Future Science Prospects for AMI

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

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.

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.

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}$