Astrophysical limits on very light axion-like particles from Chandra grating spectroscopy of NGC 1275

Axions and axion-like particles (ALPs) are a well motivated extension of the Standard Model and are generic in String Theory. The X-ray transparency of the magnetized intracluster medium (ICM) in galaxy clusters is a powerful probe of very light ALPs (masses $0<m_a<10^{-11}\,{\rm eV}$); as X-ray photons propagate through the magnetic field of the ICM, they may undergo energy-dependent quantum mechanical conversion into ALPs (and vice versa), imprinting distortions on the observed X-ray spectrum. We present new Chandra data for the active galactic nucleus NGC1275 at the center of the Perseus cluster. Employing the High-Energy Transmission Gratings (HETG) with a 490ks exposure, we obtain a high-quality 1-9keV spectrum free from photon pileup and ICM contamination. Apart from iron-band features, the spectrum is accurately described by a power-law continuum, with any spectral distortions at the $<3$% level. We compute photon survival probabilities as a function of ALP mass $m_a$ and ALP-photon coupling constant $g_{a\gamma}$ for an ensemble of ICM magnetic field models, and then use the NGC1275 spectrum to derive constraints on the $(m_a, g_{a\gamma})$-plane. Marginalizing over the magnetic field realizations, the 99.7% credible region limits the ALP-photon coupling to $g_{a\gamma}<6-8\times 10^{-13}\, {\rm GeV}^{-1}$ (depending upon the magnetic field model) for masses $m_a<1\times 10^{-12}\,{\rm eV}$. These are the most stringent limit to date on $g_{a\gamma}$ for these very light ALPs, and have already reached the sensitivity limits of next-generation helioscopes and light-shining-through-wall experiments. We highlight the potential of these studies with the next-generation X-ray observatories Athena and Lynx, but note the critical importance of advances in relative calibration of these future X-ray spectrometers.


INTRODUCTION
Astrophysical observations have great potential to uncover new physics beyond the Standard Model (SM). Indeed, the clearest experimental indications that new Quinn (PQ) mechanism which protects the strong interaction from CP (i.e. time-reversal) violating effects (Peccei & Quinn 1977;Weinberg 1978;Wilczek 1978). This Quantum Chromodynamics (QCD) axion couple to two photons with a coupling strength g aγ that is proportional to their mass (and inversely proportional to the energy scale of PQ symmetry breaking). Furthermore, many extensions of the SM, including in particular string theory, commonly feature very light axion-like particles (ALPs) that do not couple to the strong interactions, but can interact with photons with a strength that is independent of their masses. The coupling constant g aγ and ALP mass m a can, in a natural way, be small compared with other SM mass scales. Axions and ALPs may be produced in the early Universe via non-thermal mechanisms, vacuum realignment and the decay of topological defects, and hence would be produced with very small velocity dispersion. Thus, despite being very light, they can still be a viable candidate for cold dark matter (Preskill et al. 1983).
Given their coupling to photons, the transparency of astrophysical systems can be used to search for the effects of ALPs; as astrophysical photons traverse through cosmic magnetic fields, they are susceptible to conversion to an ALP via the two-photon interaction described by the Lagrangian term L a = −g aγ aE · B. (1) where a is the ALP field, and E and B are the electric and magnetic fields. The transparency of astrophysical systems can be used to set upper limits on g aγ and may be a route to the eventual detection of ALPs. X-ray observations of active galactic nuclei (AGN) in rich clusters of galaxies are particularly suited to ALP searches. Faraday rotation measure (RM) studies demonstrate that the hot intracluster medium (ICM) in such clusters is magnetized, with a ratio of thermalto-magnetic pressure of β ∼ 100 ). Furthermore, we expect significant regions of this field to be coherent on scales 1 − 10 kpc (Vacca et al. 2012). The result is that, if ALPs exist with a sufficiently high coupling to photons, clusters will be efficient converters of X-ray photons into ALPs. The conversion probability will be energy dependent thereby imprinting distortions into the observed spectrum of any embedded (or background) object, with a precise form that depends upon the magnetic field structure as well as the ALP properties. We note that these transparency studies probe the physics of ALPs independently of whether they actually constitute a significant component of the non-baryonic dark matter. Wouters & Brun (2013) used Chandra imaging spectroscopy of the central AGN Hydra-A in the Hydra cluster of galaxies to set an upper limit on spectral distortions from ALPs, showing that any ALP with mass m a < 7 × 10 −12 eV must have g aγ < 8.3 × 10 −12 GeV −1 (95% confidence level). Subsequently, Berg et al. (2017) examined the AGN NGC 1275 at the center of the Perseus cluster, the target of our current study, using Chandra and XMM-Newton imaging spectroscopy. NGC 1275 is almost 100× brighter than Hydra-A in the 2-10 keV band and, furthermore, does not have the heavy intrinsic absorption of Hydra-A. While this dramatically reduces the statistical errors on the spectrum of the NGC 1275, the high source flux creates systematic issues. For the large body of Chandra imaging data (almost 1Ms) for which NGC 1275 is close the optimal aimpoint, the sub-arcsecond focusing leads to very severe photon-pileup; thus there are very strong spectral distortions that are entirely instrumental in origin. For this reason, Berg et al. (2017) focus their attention on the smaller quantities of Chandra data for which NGC 1275 is off-axis and on data from XMM-Newton. Since the poorer focus for these datasets reduces but does not eliminate the effect of pileup, they employ the pileup model of Davis (2001) and conclude that any spectral distortions must be below the 10% (once smoothed to the spectral resolution of the CCD detectors in Chandra and XMM-Newton). The resulting ALP constraint is g aγ < 4 × 10 −12 GeV −1 for massless ALPs. Chen & Conlon (2018) extended this analysis to massive ALPs, and verified that this limit applies to m a 10 −12 eV. 1 To circumvent the photon pileup issue, Marsh et al. (2017) used 370 ks of short (0.4 s) frame time Chandra imaging spectroscopy of the core of M87 in the Virgo cluster. Combined with the fact that this AGN is almost an order of magnitude fainter than NGC1275, the short frame time exposures led to a high-quality spectrum with negligible pileup, and ALP constraints of g aγ < 2.6 × 10 −12 GeV −1 for m a < 10 −12 eV (95% confidence level).
In this paper, we present new Chandra observations of NGC 1275 that set the tightest limits to date on very light ALP conversion in the Perseus cluster. By employing the High-Energy Transmission Gratings (HETG) and investing almost 500 ks of on-source exposure, we obtain a high-quality separation of the AGN emission from the ICM with no discernible photon-pileup. Outside of the astrophysically-rich iron-band (6-7 keV), we find that the resulting 1-9 keV spectrum of the AGN is well described by a power-law continuum (modified only by the effects of absorption by cold gas in our Galaxy) with any remaining spectral distortions below the 3% level.
We proceed to marginalize over a set of representative realisations of the magnetic field to determine the posterior distribution on the (m a , g aγ )-plane. The resulting 99.7% credible region limits the ALP-photon coupling to g aγ < 6 − 8 × 10 −13 GeV −1 (depending upon the magnetic field model) for most masses m a < 1 × 10 −12 eV. This is the most stringent limit to date on the ALPphoton coupling.
For one of our magnetic field models, the posterior peaks at a non-vanishing value of g aγ , and the the 95% credible region pick out a preferred non-zero value of the couple constant for, g aγ ≈ 4 − 8 × 10 −13 GeV −1 , for a range of masses m a < 1 × 10 −12 . We attribute this "detection" to residual, low-level, errors in the Chandra/HETG calibration. This paper is organized as follows. Section 2 presents the new observations and describes some subtitles encountered during the data reduction. After discussing the modeling of the ICM magnetic field and the associated ALP distortions in Section 3, the new constraints on ALPs are given in Section 4. We put these results into context and draw our conclusions in Section 5. When necessary, we assume a standard Planck cosmology with

OBSERVATIONS, DATA REDUCTION AND INITIAL SPECTRAL FITTING
Chandra observed NGC 1275 in 15 separate segments (ObsIDs) between 24-Oct-2017 and 5-Dec-2017 using the HETG read out on the Advanced CCD Imaging Spectrometer (ACIS) S array. The total on-source exposure time was 490 ks. To ameliorate the risk of modest photon pileup in the event that the source was brighter than expected, we turned off the two outlying ACIS-S chips and used 1/2 sub-arrays on the remaining four chips, resulting in a reduction of the frame readout time to 2.4s with no loss of observing efficiency. While loss of the two outlying chips in principle affects our ability to observe the softest regions of the spectrum, the contaminant that has built up on the ACIS optical blocking filters unavoidably removes those soft photons anyway. Figure 1 shows the image of the ACIS-S array for one of our ObsIDs (20449). This is dominated by zeroth order image of core of the Perseus cluster with its famous cavity system (Fabian et al. 2000). The two two-sided dispersed spectra of the bright central AGN emission, one from the High Energy Grating (HEG) and one from the Medium Energy Grating (MEG), are clearly visible. Although the gratings are slitless and so the cluster light is also dispersed, it is clear that the AGN light is distinct and can be well isolated from the bulk of the ICM emission. Order-sorting, whereby only photons with CCDdetected energies compatible with their spatial position along the dispersion spectrum are accepted, enables further isolation of the AGN spectrum.
Much but not all of our data reduction is standard.
The extraction of the AGN spectra for each ObsIDs then follows the standard CXC science threads 2 with two exceptions. Firstly, we halve the width of the extraction regions (width factor hetg=18) in order to reduce MEG/HEG overlap at the centre of the array and hence access the higher energy band in the HEG. Secondly, the standard algorithm for automatically locating the zeroth order image and hence setting the energy scale of the spectrum failed for most of the ObsIDs, presumably due to the surrounding high-surface brightness and structured ICM. Instead, we force the zeroth order point to be at the known coordinates of NGC1275, and visually confirm for each ObsID that this correctly locates the zeroth order image of the point-like AGN (i.e. that astrometry errors are within one pixel or 0.5 arcsec). The result is four spectra and associated background spectra, response matrices and effective area files for each ObsID, namely the +1 and −1 order spectra for each of the MEG and HEG. In a final step, we combine the spectra to produce a single HEG and a single MEG spectrum (with associated background spectra, response matrices, and effective area files), summing the +/−1 orders from all ObsIDs.
All spectral fitting presented in this paper uses the 1 − 7 keV band for the MEG, the 1.5 − 9 keV band for the HEG, and employs Cash (C) statistic minimization to fit the unbinned spectrum, modified to allow for the subtraction of a background spectrum 3 . Fitting is performed with the XSPECv12.10.1 code (Arnaud 1996). An initial fit of the spectra with a powerlaw continuum modified by Galactic absorption (N H = 1.32×10 21 cm −2 ; Kalberla et al. 2005) finds 10-15% excesses in three energy bands; below 1.3 keV (MEG only), 2.2−2.5 keV (HEG and MEG), and 6−7 keV (principally in the HEG). The 6−7 keV structure corresponds to the well-established iron fluorescent line from cold gas in the vicinity of the AGN (Hitomi Collaboration et al. 2018) and will be the subject of another paper (Reynolds et al., in prep). The other residuals closely mirror structure in the background spectra (which is actually dominated by the core ICM emission), and suggests that the standard grating extraction algorithms which are developed for isolated point sources have underestimated the normalisation of the background for this embedded source. Empirically, we find that these excesses are reduced to below the 3-5% level, and the power-law is a formerly best fit, if the HEG and MEG backgrounds are scaled up by factors of 2.32 and 1.92 respectively. These factors are validated by MARX simulations of a point source embedded in a simplified model of the Perseus cluster.
With these adjusted background normalizations, the AGN photons comprise 80% and 82% of the HEG and MEG spectra, respectively. A joint HEG/MEG fit gives a power-law index of Γ = 1.890 ± 0.009 and normalization at 1keV of A = (8.28 ± 0.06) × 10 −3 ph s −1 keV −1 cm −2 , with C = 4956 for 4874 degrees of freedom (DOF). This corresponds to a 1-9 keV band flux and luminosity of F 1−9 keV = 3.8 × 10 −11 erg cm −2 s −1 and L 1−9 keV = 2.2 × 10 43 erg s −1 , respectively. If we allow the power-law parameters to float freely between the MEG and HEG spectra, the fit improves significantly (∆C = −76) and we find marginally-significant slope and highly-significant nor-malization differences between the two gratings (Table  1). While highly significant, the absolute flux difference is below 10% and hence within the realm of what can reasonably be attributed to an instrumental crosscalibration uncertainty. Hence, for all subsequent fitting in this paper, we will permit normalization and slope offsets between the HEG and MEG spectra. Figure 2 shows this free fit of the absorbed power-law to the HEG and MEG spectra, heavily binned for plotting purposes (but unbinned for fitting purposes). The residuals about the power-law are not entirely random, with a slight broad dip around 1.2 keV, a very subtle broad hump between 3-4 keV, and an obvious feature in the iron band (6.4-7 keV). Still, outside of the iron band, the remaining residuals are less than 5% and mostly less than 3% and so are entirely consistent with the expected level of residual effective area calibration errors. Having probably reached the level dominated by systematic calibration uncertainties, this is the highest quality 1-9 keV band spectrum of this AGN obtained to date and, as we shall see, permits the most sensitive search yet for very light ALPs.

MODELING THE ALP SIGNATURES
Our modeling of ALP spectral distortions follows that of Marsh et al. (2017). We solve the linearized Schrödingerlike equation that describes the quantum mechanical oscillations between photons and ALPs as they traverse through the cluster magnetic field towards the observer.  The survival probability of photons emitted from the nucleus that exit the cluster as photons, rather than ALPs, depends on the plasma density and the the cluster magnetic field. For the electron density, we use the simple analytic approximation derived from XMM-Newton observations of the Perseus cluster observations by Churazov et al. (2003): n e (r) = 3.9 × 10 −2 (1 + (r/80 kpc)) 1.8 + 4.05 × 10 −3 (1 + (r/280 kpc)) 0.87 cm −3 . (2) For the magnetic field, we use two stochastic models.
Model A: We first consider a slight modification to the magnetic fields used in Berg et al. (2017). This model is motivated by VLBA observations of the nucleus of NGC 1275 by Taylor et al. (2006), which found Our model A modifies this description in two ways. First, we note that the large value of B 0 combined with the moderate central plasma density of eq. (2) leads to non-negligible ALP-photon oscillations from the innermost region of the cluster. However, the simple analytic model of Churazov et al. (2003) does not apply to small radii, r < 10 kpc, where it underestimates the electron density, and where the spherically symmetric approximation is not justified. Applying eq. (2) to this region leads to an overestimate of the ALP-photon conversion probability. In our work, we conservatively exclude the central region, and simulate the ALP-photon oscillations from 10 kpc out to the virial radius, R vir = 1.8 Mpc. Second, we note that the cluster itself acts as a Faraday screen and sources RMs in addition to those arising from the central region. Since this model attributes the observed RMs to the Faraday screen close to the center of the cluster, we consistently select only those magnetic field configurations in which the cluster contribution is subleading: RM cluster ≤ 2000 rad m −2 . However, we have found that this restriction has no statistically significant impact on the typical conversion probabilities.

Model B:
We furthermore consider a model in which the ratio of the thermal-to-magnetic pressure is fixed to β = 100 throughout the cluster. We use the Perseus pressure profile of Fabian et al. (2006) to derive a magnetic field strength of B 25 ≈ 7.5 µG at r = 25 kpc.
Approximating the cluster as isothermal, the magnetic field decreases with radius as ∼ n e (r), where we again use (2) for the electron density. 4 The coherence lengths of the magnetic field can be expected to grow with distance from the centre. We model this effect by drawing the coherence lengths randomly from (1 + r/50 kpc) × 3.5 kpc to (1 + r/50 kpc) × 10 kpc, with a power-law fall-off as ∼ L −1.2 . This model produces Faraday RMs of the same order as those observed by Taylor et al. (2006), with the cluster as the Faraday screen.
For each of our two field models, we generate 500 RM-acceptable magnetic field configurations and solve the Schrödinger-like equation in order to calculate photon survival probabilities across a grid of m a and g aγ . Our models sample the (m a , g aγ )-plane, spanning the range log 10 (m a / eV) ∈ [−13.6, −11.1] and log 10 (g aγ /GeV −1 ) ∈ [−13, −10.7].
The result is a library of approximately 260,000 energy-dependent photon survival probability curves for each of our two magnetic field models (Model-A and Model-B), p γ (E; m a , g aγ , i A/B ), where i A/B indexes the 500 RMacceptable magnetic field realizations for that given magnetic field model. Some representative photon survival probability curves, and their functional dependence on m a and g aγ , are shown in Figure 3.

CONSTRAINTS ON ALP PARAMETERS
Equipped with the library of photon survival probability curves, we can now use our HETG spectra of NGC 1275 to determine the allowed regions of the (m a , g aγ )-plane.
In order to determine the probability of the parameters given this data, appropriately marginalised over the unknown cluster magnetic field configuration, we follow the Bayesian procedure of Marsh et al. (2017). We assume flat priors on ln m a and ln g aγ in the range log 10 (m a / eV) ∈ [−30, −11.1] and log 10 (g aγ /GeV −1 ) ∈ [−19, −10.7]. We will find that our results are insensitive to the minimum allowed mass, and the particular choice considered here corresponds to cosmologically large Compton wavelengths of the ALPs. The minimum allowed coupling constant corresponds to the inverse Planck mass, below which quantum gravitational corrections are expected to become important. We also assume flat priors on the randomly generated magnetic field configurations, labelled by i A/B . Motivated by the initial fitting presented in Section 2, our baseline spectral model for NGC 1275 consists of a power-law continuum modified by the effects of Galactic absorption (N H = 1.32 × 10 21 cm −2 ; Kalberla et al. 2005).
For a given magnetic field model (Model-A and Model-B), we take each of our photon survival probability curves (indexed by m a , g aγ and the magnetic field realization i A/B ), multiply by the absorbed power-law spectrum, and then fit to the unbinned HEG/MEG spectra, minimizing the C-statistic over the HEG/MEG photon indices and HEG/MEG normalizations. The lowest masses in our model library (log 10 (m a / eV) = −13.6) yield fits that are indistinguishable from the massless case, and hence these model fits are used as proxies for the very low-mass region of parameter space. Similarly, the smallest coupling constant in our model library (log 10 (g aγ /GeV −1 ) = −13) are indistinguishable from the zero coupling case and hence these model fits are used as proxies for the very small coupling region of parameter space.
From the resulting set of C(m a , g aγ , i A/B ), we form posterior probabilities P A/B (m a , g aγ , i) normalized such that ma,gaγ ,i A/B P A/B (m a , g aγ , i A/B ) = 1 . ( To obtain the posterior over the ALP-parameters alone (and account for the 'look-elsewhere effect' associated with the unknown magnetic field), we then marginalize over the magnetic field configurations, Contours of equal P A/B (m a , g aγ ) provide the boundaries of the credible regions of the ALP parameters. Figure 4 shows the main result of this paper, our new constraints on the (m a , g aγ )-plane. The posterior over the ALP parameters has a complex structure.
At the 95% level and for model-B, the marginalization process picks out a preferred value of g aγ ≈ 4 − 8 × 10 −13 GeV −1 . 5 An examination of a sample of the photon survival probability curves for the most probable models shows that they all share a modest dip at 1.2 keV and a broad hump at 3-4 keV, structure that is visually apparent in the spectra (Fig. 2). At these levels, however, such spectral structures could easily result from remaining calibration errors in the HETG energydependence effective area and thus we cannot claim even a tentative detection of ALPs on the basis of the enclosed 95% confidence contour.

DISCUSSION AND CONCLUSIONS
A deep exposure of NGC 1275 at the center of the Perseus cluster of galaxies with the Chandra/HETG has allowed us to obtain the highest quality spectrum of this AGN free from strong contamination by the surrounding cluster emission and free from the effects of photon pileup. Apart from subtle structure in the iron-band (6-7 keV), we find that the 1-9 keV spectrum of NGC 1275 is accurately described by a power-law continuum form modified by the effects of modest Galactic absorption; deviations from the power-law are at the ±3% level. Taking this to be the intrinsic spectrum of the AGN, we proceed to use these data to constrain models for ALPphoton oscillations in the magnetic field of the Perseus cluster. We have obtained the most stringent limit yet on the ALP-photon coupling constant of very light ALPs 5 We note that the strength of the preference for a non-vanishing value of the coupling gaγ is prior-dependent: should we have restricted the prior range to only the sampled region of log 10 (gaγ /GeV −1 ) ∈ [−13, −10.7], the 99.7% confidence regions obtained from either of the magnetic field models would exclude a vanishing ALP-photon coupling. In contrast, the upper limit on gaγ is highly prior-independent, and the large prior range considered in this paper is therefore conservative. Figure 4. Constraints on the (ma, gaγ)-plane from this study for magnetic field Model-A (top) and Model-B (bottom). We show 99.7% confidence limit (heavy black line bounding the excluded, light blue shaded region), the 95% confidence region (allowed region with yellow shading), and best-fit parameter values (red dots). Also shown are the previous most stringent 95% limit from the Marsh et al. (2017) study of M87 (blue line), the limit from the absence of a γ-ray burst associated with SN1987 (Payez et al. 2015) (black line), the limit from FERMI (gray line bounding blue excluded region at the 95% confidence level) as well as the projected sensitivity from the next generation helioscope IAXO. For other, weaker limits obtained using X-ray astronomy in the same region of parameter space, see also Wouters & Brun (2013); Berg et al. (2017); .
with masses m a < 1 × 10 −12 eV: g aγ < 7.9 × 10 −13 GeV −1 , (99.7% confidence), (5) for magnetic field model-B, with an even more stringent limit of g aγ < 6.3 × 10 −13 GeV −1 for model-A. Even at this much higher level of confidence, our limits are 3−4× stronger than those obtained by Marsh et al. (2017) for M87/Virgo, and over 5× stronger than the those found by Berg et al. (2017) with non-dispersive spectroscopy of NGC 1275/Perseus. For one of our two magnetic field models, the ALP conversion models pick out some remaining structure in the HEG/MEG spectra and hence, at the 95% level (but not the 99.7% level) lock onto a preferred non-zero coupling constant of g aγ ≈ 4 − 8 × 10 −13 GeV −1 . However, we acknowledge that ±3% is very plausibly at the level of residual calibration uncertainties in the HETG energydependent effective area and hence this result must be viewed with extreme caution.
For these very light ALPs (m a < 1 × 10 −11 eV) our astrophysical limits are now tighter than those obtained by the fact that SN1987A did not generate a gammaray burst via the ALP-mediated escape of gamma-rays from the collapsing stellar core (Payez et al. 2015). Our limits have also exceeded the projected sensitivity of the next-generation helioscope, the International Axion Observatory (IAXO) (Armengaud et al. 2014), as well as the next generation light-shining-through-walls experiments such as ALPS-III 6 .
While the focus of this paper are the constraints on ALPs, we note that our data allow a strong test of the  fluorescent dark matter model. These authors note that there are hints of an absorption feature at 3.5 keV in the XMM-Newton/EPIC-pn of NGC 1275. In an attempt to reconcile the claim of a 3.5 keV dark matter emission line in the XMM-Newton spectrum of the Perseus ICM (Bulbul et al. 2014;Boyarsky et al. 2014) with the non-detection of any such feature in the Hitomi spectrum of the the system (Hitomi Collaboration et al. 2017),  proceed to formulate a two-level dark matter model in which a dark matter absorption line in the AGN spectrum offsets the dark matter emission line from the cluster in system-integrated spectra such as that produced by Hitomi. This hypothesis requires that the dark matter imprints an absorption line close to 3.5 keV with an equivalent width of 15 eV. Our HETG current spectrum sets an upper limit of 4 eV with 99% confidence, allowing us 6 https://indico.cern.ch/event/523655/contributions/ 2246871/attachments/1332644/2003496/LSW.pdf to rule out this version of the fluorescent dark matter model. Following on from Wouters & Brun (2013), Berg et al. (2017) and Marsh et al. (2017), our work is just the latest to highlight the power of X-ray transparency studies of galaxy clusters to probe the physics of ALPs. Equipped with high-resolution, non-dispersion, high count rate micro-calorimater arrays at the focus of X-ray telescopes with A > 1 m 2 , the next generation of X-ray flagship missions Athena and Lynx will dramatically advance the quality of the spectrum possible for an embedded AGN such as NGC 1275.
But the current study also points to the significant challenges that need to be overcome if we are extend this technique further and dig deeper into the parameter space. Even if photon pileup/deadtime and contamination by ICM emission is rendered negligible, the relative calibration of the energy-dependent effective area of the X-ray spectrometer will set a floor on the sensitivity of any studies based on spectral distortions irrespective of the photon statistics. For example, for a representative Perseus magnetic field model, a very light ALP with g aγ = 3 × 10 −13 GeV −1 produces distortions at the 1.5% (3%) level in the E ∼ 2 − 4 keV(5 − 8 keV) band. For g aγ = 2 × 10 −13 GeV −1 , the corresponding distortion is 0.75% (1.5%). Thus, to significantly further this technique for searching for ALPs, future X-ray spectrometers must achieve relative effective area calibrations of 1% or better. It is clear that in order to push the ALP constraints significantly further than we have done here (or actually obtain a robust detection of very light ALPs) will require advances in broad-band, on-orbit, relative calibration of future X-ray spectrometers.
Astronomy, where this work was concluded. S.V. also acknowledges support by the Science and Technology Facilities Council (STFC) and by the Kavli Institute for Cosmology, Cambridge. This work used, in part, the Cambridge COSMOS SMP system, part of the STFC DiRAC HPC Facility supported by BIS NeI capital grant ST/J005673/1 and STFC grants ST/H008586/1, ST/K00333X/1.