Accretion Column Research Articles

Discovery of new dipper stars with K2: a window into the inner disc region of T Tauri stars

In recent years a new class of Young Stellar Object has been defined, referred to as dippers, where large transient drops in flux are observed. These dips are too large to be attributed to stellar variability, last from hours to days and can reduce the flux of a star by 10-50\%. This variability has been attributed to occultations by warps or accretion columns near the inner edge of circumstellar disks. Here we present 95 dippers in the Upper Scorpius association and $\rho$ Ophiuchus cloud complex found in K2 Campaign 2 data using supervised machine learning with a Random Forest classifier. We also present 30 YSOs that exhibit brightening events on the order of days, known as bursters. Not all dippers and bursters are known members, but all exhibit infrared excesses and are consistent with belonging to either of the two young star forming regions. We find 21.0 $\pm$ 5.5\% of stars with disks are dippers for both regions combined. Our entire dipper sample consists only of late-type (KM) stars, but we show that biases limit dipper discovery for earlier spectral types. Using the dipper properties as a proxy, we find that the temperature at the inner disk edge is consistent with interferometric results for similar and earlier type stars.

Measuring the masses of intermediate polars with NuSTAR: V709 Cas, NY Lup, and V1223 Sgr

The X-ray spectra of intermediate polars can be modelled to give a direct measurement of white dwarf mass. Here we fit accretion column models to NuSTAR spectra of three intermediate polars; V709 Cas, NY Lup and V1223 Sgr in order to determine their masses. From fits to 3-78 keV spectra, we find masses of $M_{\rm WD}=0.88^{+0.05}_{-0.04}M_{\odot}$, $1.16^{+0.04}_{-0.02}M_{\odot}$ and $0.75\pm0.02M_{\odot}$ for V709 Cas, NY Lup and V1223 Sgr, respectively. Our measurements are generally in agreement with those determined by previous surveys of intermediate polars, but with typically a factor $\sim2$ smaller uncertainties. This work paves the way for an approved NuSTAR Legacy Survey of white dwarf masses in intermediate polars.

A Non-thermal Pulsed X-Ray Emission of AR Scorpii

Abstract We report the analysis result of UV/X-ray emission from AR Scorpii, which is an intermediate polar (IP) composed of a magnetic white dwarf and an M-type star, with the XMM-Newton data. The X-ray/UV emission clearly shows a large variation over the orbit, and their intensity maximum (or minimum) is located at the superior conjunction (or inferior conjunction) of the M star orbit. The hardness ratio of the X-ray emission shows a small variation over the orbital phase and shows no indication of the absorption by an accretion column. These properties are naturally explained by the emission from the M star surface rather than that from the accretion column on the white dwarf’s (WD) star, which is similar to usual IPs. Additionally, the observed X-ray emission also modulates with the WD’s spin with a pulse fraction of ∼14%. The peak position is aligned in the optical/UV/X-ray band. This supports the hypothesis that the electrons in AR Scorpii are accelerated to a relativistic speed and emit non-thermal photons via the synchrotron radiation. In the X-ray bands, evidence of the power-law spectrum is found in the pulsed component, although the observed emission is dominated by the optically thin thermal plasma emissions with several different temperatures. It is considered that the magnetic dissipation/reconnection process on the M star surface heats up the plasma to a temperature of several keV and also accelerates the electrons to the relativistic speed. The relativistic electrons are trapped in the WD’s closed magnetic field lines by the magnetic mirror effect. In this model, the observed pulsed component is explained by the emissions from the first magnetic mirror point.

X-ray reflection from cold white dwarfs in magnetic cataclysmic variables

We model X-ray reflection from white dwarfs (WD) in magnetic cataclysmic variables (mCVs) with a Monte Carlo simulation. A point source with a power-law spectrum or a realistic post-shock accretion column (PSAC) source irradiates a cool and spherical WD. The PSAC source emits thermal spectra of various temperatures stratified along the column according to the PSAC model. In the point source simulation, we confirm (1) a source harder and nearer to the WD enhances the reflection, (2) higher iron abundance enhances the equivalent widths (EWs) of fluorescent iron K${\alpha}_{1,2}$ lines and their Compton shoulder, and increases cut-off energy of a Compton hump, and (3) significant reflection appears from an area that is more than 90 deg apart from the position right under the point X-ray source because of the WD curvature. The PSAC simulation reveals that (1) a more massive WD basically enhances the intensities of the fluorescent iron K${\alpha}_{1,2}$ lines and the Compton hump, except for some specific accretion rate, because the more massive WD makes the hotter PSAC from which higher energy X-rays are preferentially emitted, (2) a larger specific accretion rate monotonically enhances the reflection because it makes the hotter and shorter PSAC, and (3) the intrinsic thermal component hardens by occultation of the cool base of the PSAC by the WD. We quantitatively evaluate influences of the parameters on the EWs and the Compton hump with both types of sources. We also calculate X-ray modulation profiles brought about by the WD spin. They depend on angles of the spin axis from the line of sight and from the PSAC, and whether the two PSACs can be seen. The reflection spectral model and the modulation model involve the fluorescent lines and the Compton hump and can directly be compared to the data which allow us to evaluate these geometrical parameters with unprecedented accuracy.

Looking at A 0535+26 at low luminosities with NuSTAR

We report on two NuSTAR observations of the HMXB A 0535+26 taken toward the end of its normal 2015 outburst at very low $3-50$ keV luminosities of ${\sim}1.4\times10^{36}$ erg/s and ${\sim}5\times10^{35}$ erg/s which are complemented by 9 Swift observations. The data clearly confirm indications seen in earlier data that the source's spectral shape softens as it becomes fainter. The smooth, exponential rollover at high energies present in the first observation evolves to a much more abrupt steepening of the spectrum at $20-30$ keV. The continuum evolution can be well described with emission from a magnetized accretion column, modeled using the compmag model modified by an additional Gaussian emission component for the fainter observation. Between the two observations, the optical depth changes from $0.75\pm0.04$ to $0.56^{+0.01}_{-0.04}$, the electron temperature remains constant, and there is an indication that the column decreases in radius. Since the energy resolved pulse profiles remain virtually unchanged in shape between the two observations, the emission properties of the accretion column, however, reflect the same accretion regime. This conclusion is also confirmed by our result that the energy of the cyclotron resonant scattering feature (CRSF) at ${\sim}45$ keV is independent of the luminosity, implying that the magnetic field in the region in which the observed radiation is produced is the same in both observations. Finally, we also constrain the evolution of the continuum parameters with rotational phase of the neutron star. The width of the CRSF could only be constrained for the brighter observation. Based on Monte-Carlo simulations of CRSF formation in single accretion columns, its pulse phase dependence supports a simplified fan beam emission pattern. The evolution of the CRSF width is very similar to that of the CRSF depth, which is in disagreement with expectations.

MHD simulations of resistive viscous accretion disk around millisecond pulsar

We perform MHD simulations of a thin resistive and viscous accretion disk around a neutron star with the surface dipolar magnetic field of 108 Gauss. The system evolution is followed during the interval of 500 millisecond pulsar rotations. Matter is accreted through a stable accretion column from the disk onto the star. We also show propagation of the stellar wind through the corona. Analysis of the mass accretion flux and torques on the star shows that the disk reaches the quasi-stationary state.

On the radiation beaming of bright X-ray pulsars and constraints on neutron star mass–radius relation

The luminosity of accreting magnetised neutron stars can largely exceed the Eddington value due to appearance of accretion columns. The height of the columns can be comparable to the neutron star radius. The columns produce the X-rays detected by the observer directly and illuminate the stellar surface, which reprocesses the X-rays and causes additional component of the observed flux. The geometry of the column and the illuminated part of the surface determines the radiation beaming. Curved space-time affects the angular flux distribution. We construct a simple model of the beam patterns formed by direct and reflected flux from the column. We take into account the possibility of appearance of accretion columns, whose height is comparable to the neutron star radius. We argue that depending on the compactness of the star the flux from the column can be either strongly amplified due to gravitational lensing, or significantly reduced due to column eclipse by the star. The eclipses of high accretion columns result in specific features in pulse profiles. Their detection can put constraints on the neutron star radius. We speculate that column eclipses are observed in X-ray pulsar V 0332+53, leading us to the conclusion of large neutron star radius in this system ($\sim15\,{\rm km}$ if $M\sim 1.4M_\odot$). We point out that the beam pattern can be strongly affected by scattering in the accretion channel at high luminosity, which has to be taken into account in the models reproducing the pulse profiles.

An estimation of the white dwarf mass in the Dwarf Nova GK Persei with NuSTAR observations of two states

We report on X-ray observations of the Dwarf Nova GK Persei performed by {\it NuSTAR} in 2015. GK Persei, behaving also as an Intermediate Polar, exhibited a Dwarf Nova outburst in 2015 March--April. The object was observed with {\sl NuSTAR} during the outburst state, and again in a quiescent state wherein the 15--50 keV flux was 33 times lower. Using a multi-temperature plasma emission and reflection model, the highest plasma temperature in the accretion column was measured as $19.7^{+1.3}_{-1.0}$~keV in outburst and $36.2^{+3.5}_{-3.2}$~keV in quiescence. The significant change of the maximum temperature is considered to reflect an accretion-induced decrease of the inner-disk radius $R_{\rm in}$, where accreting gas is captured by the magnetosphere. Assuming this radius scales as $R_{\rm in} \propto \dot{M}^{-2/7}$ where $\dot{M}$ is the mass accretion rate, we obtain $R_{\rm in} = 1.9 ^{+0.4}_{-0.2}~R_{\rm WD}$ and $R_{\rm in} = 7.4^{+2.1}_{-1.2}~R_{\rm WD}$ in outburst and quiescence respectively, where $R_{\rm WD}$ is the white-dwarf radius of this system. Utilising the measured temperatures and fluxes, as well as the standard mass-radius relation of white dwarfs, we estimate the white-dwarf mass as $M_{\rm WD} = 0.87~\pm~0.08~M_{\rm \odot}$ including typical systematic uncertainties by 7%. The surface magnetic field is also measured as $B \sim 5 \times 10^{5}$~G. These results exemplify a new X-ray method of estimating $M_{\rm WD}$ and $B$ of white dwarfs by using large changes in $\dot{M}$.

Quasi-periodic oscillations from post-shock accretion column of polars

A set of strongly magnetized accreting white dwarfs (polars) shows quasi-periodic oscillations (QPOs) with frequency about a Hz in their optical luminosity. These Hz-frequency QPOs are thought to be generated by intensity variations of the emitted radiation originating at the post-shock accretion column. Thermal instability in the post-shock region, triggered by efficient cooling process at the base, is believed to be the primary reason behind the temporal variability. Here we study the structure and the dynamical properties of the post-shock accretion column including the effects of bremsstrahlung and cyclotron radiation. We find that the presence of significant cyclotron emission in optical band reduces the overall variability of the post-shock region. In the case of a larger post-shock region above the stellar surface, the effects of stratification due to stellar gravity becomes important. An accretion column, influenced by the strong gravity, has a smaller variability as the strength of the thermal instability at the base of the column is reduced. On the other hand, the cool, dense plasma, accumulated just above the stellar surface, may enhance the post-shock variability due to the propagation of magnetic perturbations. These characteristics of the post-shock region are consistent with the observed properties of V834~Cen and in general with Cataclysmic Variable sources that exhibit QPO frequency of about a Hz.

Super-Eddington accretion on to the neutron star NGC 7793 P13: Broad-band X-ray spectroscopy and ultraluminous X-ray sources

We present a detailed, broadband X-ray spectral analysis of the ULX pulsar NGC 7793 P13, a known super-Eddington source, utilizing data from the $XMM$-$Newton$, $NuSTAR$ and $Chandra$ observatories. The broadband $XMM$-$Newton+NuSTAR$ spectrum of P13 is qualitatively similar to the rest of the ULX sample with broadband coverage, suggesting that additional ULXs in the known population may host neutron star accretors. Through time-averaged, phase-resolved and multi-epoch studies, we find that two non-pulsed thermal blackbody components with temperatures $\sim$0.5 and $\sim$1.5 keV are required to fit the data below 10 keV, in addition to a third continuum component which extends to higher energies and is associated with the pulsed emission from the accretion column. The characteristic radii of the thermal components appear to be similar, and are too large to be associated with the neutron star itself, so the need for two components likely indicates the accretion flow outside the magnetosphere is complex. We suggest a scenario in which the thick inner disc expected for super-Eddington accretion begins to form, but is terminated by the neutron star's magnetic field soon after its onset, implying a limit of $B \lesssim 6 \times 10^{12}$ G for the dipolar component of the central neutron star's magnetic field. Evidence of similar termination of the disc in other sources may offer a further means of identifying additional neutron star ULXs. Finally, we examine the spectrum exhibited by P13 during one of its unusual 'off' states. These data require both a hard powerlaw component, suggesting residual accretion onto the neutron star, and emission from a thermal plasma, which we argue is likely associated with the P13 system.

Numerical simulations of high-energy flows in accreting magnetic white dwarfs

Some polars show quasi-periodic oscillations (QPO) in their optical light curves which have been interpreted as the result of shock oscillations driven by the cooling instability. Although numerical simulations can recover this physics, they wrongly predict QPOs in the X-ray luminosity and have also failed to reproduce the observed frequencies, at least for the limited range of parameters explored so far. Given the uncertainties on the observed polar parameters, it is still unclear whether simulations can reproduce the observations. The aim of this work is to study QPOs covering all relevant polars showing QPOs. We perform numerical simulations including gravity, cyclotron and bremsstrahlung radiative losses, for a wide range of polar parameters, and compare our results with the astronomical data using synthetic X-ray and optical luminosities.We show that shock oscillations are the result of complex shock dynamics triggered by the interplay of two radiative instabilities. The secondary shock forms at the acoustic horizon in the post-shock region in agreement with our estimates from steady-state solutions. We also demonstrate that the secondary shock is essential to sustain the accretion shock oscillations at the average height predicted by our steady-state accretion model. Finally, in spite of the large explored parameter space, matching the observed QPO parameters requires a combination of parameters inconsistent with the observed ones. This difficulty highlights the limits of one-dimensional simulations, suggesting that multi-dimensional effects are needed to understand the non-linear dynamics of accretion columns in polars and the origins of QPOs.

Studying the accretion geometry of EXO 2030+375 at luminosities close to the propeller regime

The Be X-ray binary EXO 2030+375 was in an extended low luminosity state during most of 2016. We observed this state with NuSTAR and Swift, supported by INTEGRAL observations as well as optical spectroscopy with the NOT. We present a comprehensive spectral and timing analysis of these data here to study the accretion geometry and investigate a possible onset of the propeller effect. The H-alpha data show that the circumstellar disk of the Be-star is still present. We measure equivalent widths similar to values found during more active phases in the past, indicating that the low-luminosity state is not simply triggered by a smaller Be disk. The NuSTAR data, taken at a 3-78 keV luminosity of ~6.8e35 erg/s (for a distance of 7.1 kpc), are well described by standard accreting pulsar models, such as an absorbed power-law with a high-energy cutoff. We find that pulsations are still clearly visible at these luminosities, indicating that accretion is continuing despite the very low mass transfer rate. In phase-resolved spectroscopy we find a peculiar variation of the photon index from ~1.5 to ~2.5 over only about 3% of the rotational period. This variation is similar to that observed with XMM-Newton at much higher luminosities. It may be connected to the accretion column passing through our line of sight. With Swift/XRT we observe luminosities as low as 1e34 erg/s during which the data quality did not allow us to search for pulsations, but the spectrum is much softer and well described by either a blackbody or soft power-law continuum. This softer spectrum might be due to the fact that accretion has been stopped by the propeller effect and we only observe the neutron star surface cooling.

Decade long RXTE monitoring observations of Be/X-ray binary pulsar EXO 2030+375

We present a comprehensive timing and spectral studies of Be/X-ray binary pulsar EXO 2030+375 using extensive Rossi X-ray Timing Explorer observations from 1995 till 2011, covering numerous Type I and 2006 Type II outbursts. Pulse profiles of the pulsar were found to be strongly luminosity dependent. At low luminosity, the pulse profile consisted of a main peak and a minor peak that evolved into a broad structure at high luminosity with a significant phase shift. A narrow and sharp absorption dip, also dependent on energy and luminosity, was detected in the pulse profile. Comparison of pulse profiles showed that the features at a particular luminosity are independent of type of X-ray outbursts. This indicates that the emission geometry is solely a function of mass accretion rate. The broadband energy spectrum was described with a partial covering high energy cutoff model as well as a physical model based on thermal and bulk Comptonization in accretion column. We did not find any signature of cyclotron resonance scattering feature in the spectra obtained from all the observations. A detailed analysis of spectral parameters showed that, depending on source luminosity, the power-law photon index was distributed in three distinct regions. It suggests the phases of spectral transition from sub-critical to super-critical regimes in the pulsar as proposed theoretically. A region with constant photon index was also observed in ~(2-4) x 10^37 erg/s range, indicating critical luminosity regime in EXO 2030+375.

SMC X-3: the closest ultraluminous X-ray source powered by a neutron star with non-dipole magnetic field

Magnetic field of accreting neutron stars determines their overall behaviour including the maximum possible luminosity. Some models require an above-average magnetic field strength (> 10^13 G) in order to explain super-Eddington mass accretion rate in the recently discovered class of pulsating ultraluminous X-ray sources (ULX). The peak luminosity of SMC X-3 during its major outburst in 2016-2017 reached ~2.5x10^39 erg/s comparable to that in ULXs thus making this source the nearest ULX-pulsar. SMC X-3 belongs to the class of transient X-ray pulsars with Be optical companions, and exhibited a giant outburst in July 2016 - February 2017. The source has been observed during the entire outburst with the Swift/XRT and Fermi/GBM telescopes, as well as the NuSTAR observatory. Collected data allowed us to estimate the magnetic field strength of the neutron star in SMC X-3 using several independent methods. Spin evolution of the source during and between the outbursts and the luminosity of the transition to so-called propeller regime in the range of (0.3 - 7)x10^35 erg/s imply relatively weak dipole field of (1 - 5)x10^12 G. On the other hand, there is also evidence for much stronger field in the immediate vicinity of the neutron star surface. In particular, transition from super- to sub-critical accretion regime associated with cease of the accretion column, absence of cyclotron absorption features in the broadband X-ray spectrum of the source obtained with NuSTAR and very high peak luminosity favor an order of magnitude stronger field. This discrepancy makes SMC X-3 a good candidate to posses significant non-dipolar components of the field, and an intermediate source between classical X-ray pulsars and accreting magnetars which may constitute an appreciable fraction of ULX population.

General Relativistic Radiation MHD Simulations of Supercritical Accretion onto a Magnetized Neutron Star: Modeling of Ultraluminous X-Ray Pulsars

Abstract By performing 2.5-dimensional general relativistic radiation magnetohydrodynamic simulations, we demonstrate supercritical accretion onto a non-rotating, magnetized neutron star, where the magnetic field strength of dipole fields is 1010 G on the star surface. We found the supercritical accretion flow consists of two parts: the accretion columns and the truncated accretion disk. The supercritical accretion disk, which appears far from the neutron star, is truncated at around ≃3 R * (R * = 106 cm is the neutron star radius), where the magnetic pressure via the dipole magnetic fields balances with the radiation pressure of the disks. The angular momentum of the disk around the truncation radius is effectively transported inward through magnetic torque by dipole fields, inducing the spin up of a neutron star. The evaluated spin-up rate, ∼−10−11 s s−1, is consistent with the recent observations of the ultraluminous X-ray pulsars. Within the truncation radius, the gas falls onto a neutron star along the dipole fields, which results in a formation of accretion columns onto the northern and southern hemispheres. The net accretion rate and the luminosity of the column are ≃66 L Edd/c 2 and ≲10 L Edd, where L Edd is the Eddington luminosity and c is the light speed. Our simulations support a hypothesis whereby the ultraluminous X-ray pulsars are powered by the supercritical accretion onto the magnetized neutron stars.

Features of the accretion in the EX Hydrae system: Results of numerical simulation

A two dimensional numerical model in the axisymmetric approximation that describes the flow structure in the magnetosphere of the white dwarf in the EX Hya system has been developed. Results of simulations show that the accretion in EX Hya proceeds via accretion columns, that are not closed and have curtain-like shapes. The thickness of the accretion curtains depends only weakly on the thickness of the accretion disk. This thickness developed in the simulations does not agree with observations. It is concluded that the main reason for the formation of thick accretion curtains in the used model is the assumption that the magnetic field penetrates fully into the plasma of the disk. An analysis based on simple estimates shows that a diamagnetic disk that fully or partially shields the magnetic field of the star may be a more attractive explanation for the observed features of the accretion in EX Hya.

Super-Eddington accretion on to a magnetized neutron star

Most of ultraluminous X-ray sources are thought to be objects accreting above their Eddington limits. In the recently identified class of ultraluminous X-ray pulsars, accretor is a neutron star and thus has a fairly small mass with a small Eddington limit. The accretion disc structure around such an object affects important observables such as equilibrium period, period derivative and the size of the magnetosphere. We propose a model of a nearly-standard accretion disc interacting with the magnetosphere only in a thin layer near the inner disc rim. Our calculations show that the size of the magnetosphere may be represented as the classical Alfv\'en radius times a dimensionless factor $\xi$ which depends on the disc thickness only. In the case of radiation-pressure-dominated disc, the size of the magnetosphere does not depend on the mass accretion rate. In general, increasing the disc thickness leads to a larger magnetosphere size in units of the Alfv\'en radius. For large enough mass accretion rates and magnetic moments, it is important to take into account not only the pressure of the magnetic field and the radiation pressure inside the disc, but also the pressure of the radiation produced close to the surface of the neutron star in accretion column. The magnetospheric size may increase by up to factor of two as a result of the effects related to the disc thickness and the irradiation from the central source. Accounting for these effects reduces the estimate of the neutron star magnetic moment by a factor of several.

Post-flare formation of the accretion stream and a dip in pulse profiles of LMC X–4

We report here a pulse profile evolution study of an accreting X-ray pulsar LMC X–4 during and after the large X-ray flares using data from the two observatories XMM-Newton and RXTE. During the flares, the pulse profiles were found to have a significant phase offset in the range of 0.2–0.5 compared to the pulse profiles immediately before or after the flare. Investigating the pulse profiles for about 105 s after the flares, it was found that it takes about 2000–4000 s for the modified accretion column to return to its normal structure and formation of an accretion stream that causes a dip in the pulse profile of LMC X–4. We have also carried out pulse phase resolved spectroscopy of LMC X–4 in narrow phase bins using data from EPIC-pn and spectroscopically confirmed the pulsating nature of the soft spectral component, having a pulse fraction and phase different from that of the power-law component.

A broad-band self-consistent modelling of the X-ray spectrum of 4U 1626−67

The accretion-powered X-ray pulsar 4U 1626--67 is one of the few highly magnetized pulsars that accretes through Roche-lobe overflow from a low-mass companion. The characteristics of its broadband spectrum are similar to those of X-ray pulsars hosted in a high-mass X-ray binary systems, with a broad resonant cyclotron scattering feature (CRSF) at $\sim$37 keV. In this work, we examine the pulse-resolved and the pulse-averaged broadband spectrum using data from NuSTAR and Swift. We use the Becker & Wolff model of bulk+thermal Comptonization to infer key physical parameters of the accretion column flow and a broadband model for the disk reflected spectrum. In the softer X-ray band, we need to add a soft black-body component with $kT_{\rm bb}$ $\sim$0.5 keV, whose characteristics indicate a possible origin from the neutron star surface. Residuals suggest that the shape of the cyclotron line could be more satisfactorily fitted using a narrow core and broader wings and, at higher energies, a second harmonic could be present at $\sim$61 keV.

X-ray constraints on the fraction of obscured active galactic nuclei at high accretion luminosities

The wide-area XMM-XXL X-ray survey is used to explore the fraction of obscured AGN at high accretion luminosities, $L_X (\rm 2-10 \, keV) > 10^{44} \, erg \,s ^{-1}$, and out to redshift $z\approx1.5$. The sample covers an area of about $\rm14\,deg^2$ and provides constraints on the space density of powerful AGN over a wide range of neutral hydrogen column densities extending beyond the Compton-thick limit, $\rm N_H\approx10^{24}\,cm^{-2}$. The fraction of obscured Compton-thin ($\rm N_H=10^{22}-10^{24}\,cm^{-2}$) AGN is estimated to be $\approx0.35$ for luminosities $L_X(\rm 2-10\,keV)>10^{44}\,erg\,s^{-1}$ independent of redshift. For less luminous sources the fraction of obscured Compton-thin AGN increases from $0.45\pm0.10$ at $z=0.25$ to $0.75\pm0.05$ at $z=1.25$. Studies that select AGN in the infrared via template fits to the observed Spectral Energy Distribution of extragalactic sources estimate space densities at high accretion luminosities consistent with the XMM-XXL constraints. There is no evidence for a large population of AGN (e.g. heavily obscured) identified in the infrared and missed at X-ray wavelengths. We further explore the mid-infrared colours of XMM-XXL AGN as a function of accretion luminosity, column density and redshift. The fraction of XMM-XXL sources that lie within the mid-infrared colour wedges defined in the literature to select AGN is primarily a function of redshift. This fraction increases from about 20-30% at z=0.25 to about 50-70% at $z=1.5$.