Context. It has been shown that for the analysis of X-ray spectra the C-statistic, contrary to the χ2-statistic, provides unbiased estimates of the model parameters and their uncertainty ranges. ...Aims. However, it is often stated that the C-statistic cannot be used to carry out statistical tests on the goodness of fit of the model, and therefore several investigations are still based on χ2-statistics. Methods. Here we show that it is straightforward to calculate the expected value and variance of the C-statistic so that it can be used in tests. Results. We provide formulae and simple numerical approximations to evaluate these expected values and variances. We also give examples indicating that tests based on only the expected value and variance of the C-statistic are reliable for spectra even with only ~30 counts. Conclusions. The C-statistic can be used for statistical tests such as assessing the goodness of fit of a spectral model.
Aims. A theoretical framework is developed to estimate the optimal binning of X-ray spectra. Methods. We derived expressions for the optimal bin size for model spectra as well as for observed data ...using different levels of sophistication. Results. It is shown that by taking into account both the number of photons in a given spectral model bin and their average energy over the bin size, the number of model energy bins and the size of the response matrix can be reduced by a factor of 10−100. The response matrix should then contain the response at the bin centre as well as its derivative with respect to the incoming photon energy. We provide practical guidelines for how to construct optimal energy grids as well as how to structure the response matrix. A few examples are presented to illustrate the present methods.
Atomic data and plasma models play a crucial role in the diagnosis and interpretation of astrophysical spectra, thus influencing our understanding of the Universe. In this investigation we present a ...systematic comparison of the leading photoionisation codes to determine how much their intrinsic differences impact X-ray spectroscopic studies of hot plasmas in photoionisation equilibrium. We carry out our computations using the Cloudy, SPEX, and XSTAR photoionisation codes, and compare their derived thermal and ionisation states for various ionising spectral energy distributions. We examine the resulting absorption-line spectra from these codes for the case of ionised outflows in active galactic nuclei. By comparing the ionic abundances as a function of ionisation parameter ξ, we find that on average there is about 30% deviation between the codes in ξ where ionic abundances peak. For H-like to B-like sequence ions alone, this deviation in ξ is smaller at about 10% on average. The comparison of the absorption-line spectra in the X-ray band shows that there is on average about 30% deviation between the codes in the optical depth of the lines produced at log ξ ~ 1 to 2, reducing to about 20% deviation at log ξ ~ 3. We also simulate spectra of the ionised outflows with the current and upcoming high-resolution X-ray spectrometers, on board XMM-Newton, Chandra, Hitomi, and Athena. From these simulations we obtain the deviation on the best-fit model parameters, arising from the use of different photoionisation codes, which is about 10 to 40%. We compare the modelling uncertainties with the observational uncertainties from the simulations. The results highlight the importance of continuous development and enhancement of photoionisation codes for the upcoming era of X-ray astronomy with Athena.
It has been known for decades that the observed number of baryons in the local Universe falls about 30-40 per cent short
of the total number of baryons predicted
by Big Bang nucleosynthesis, as ...inferred
from density fluctuations of the cosmic microwave background and seen during the first 2-3 billion years of the Universe in the so-called 'Lyman α forest'
(a dense series of intervening H I Lyman α absorption lines in the optical spectra of background quasars). A theoretical solution to this paradox locates the missing baryons in the hot and tenuous filamentary gas between galaxies, known as the warm-hot intergalactic medium. However, it is difficult to detect them there because the largest by far constituent of this gas-hydrogen-is mostly ionized and therefore almost invisible in far-ultraviolet spectra with typical signal-to-noise ratios
. Indeed, despite large observational efforts, only a few marginal claims of detection have been made so far
. Here we report observations of two absorbers of highly ionized oxygen (O VII) in the high-signal-to-noise-ratio X-ray spectrum of a quasar at a redshift higher than 0.4. These absorbers show no variability over a two-year timescale and have no associated cold absorption, making the assumption that they originate from the quasar's intrinsic outflow or the host galaxy's interstellar medium implausible. The O VII systems lie in regions characterized by large (four times larger than average
) galaxy overdensities and their number (down to the sensitivity threshold of our data) agrees well with numerical simulation predictions for the long-sought warm-hot intergalactic medium. We conclude that the missing baryons have been found.
We present the results of a search for unidentified emission lines in deep Suzaku X-ray spectra of the central regions of the X-ray brightest galaxy clusters: Perseus, Coma, Virgo and Ophiuchus. We ...analyse an optimized energy range (3.2–5.3 keV) that is relatively free of instrumental features, and a plasma emission model incorporating the abundances of elements with the strongest expected emission lines at these energies (S, Ar, Ca) as free parameters. For the Perseus Cluster core, employing this model, we find evidence for an additional emission feature at an energy
$E=3.51^{+0.02}_{-0.01}$
keV with a flux of
$2.87_{-0.38}^{+0.33}\times 10^{-7}\,{\rm photons}\,{\rm s}^{-1}\,{\rm cm}^{-2}\,{\rm arcmin}^{-2}$
. At slightly larger radii, we detect an emission line at 3.59 ± 0.02 keV with a flux of
$4.8_{-1.4}^{+1.7}\times 10^{-8}\,{\rm photons}\,{\rm s}^{-1}\,{\rm cm}^{-2}\,{\rm arcmin}^{-2}$
. The properties of these features are broadly consistent with previous claims, although the radial variation of the line strength appears in tension with dark matter (DM) decay model predictions. Assuming a decaying DM origin for these features allows us to predict the energies and detected line fluxes for the other clusters. We do not detect an emission feature at the predicted energy and line flux in the Coma, Virgo and Ophiuchus clusters. The formal 99.5 per cent upper limits on the line strengths in each cluster are well below the decaying DM model predictions, disfavouring a decaying DM interpretation. The results of further analysis suggest that systematic effects associated with modelling the spectra for the Perseus Cluster, details of the assumed ionization balance and errors in the predicted spectral line emissivities may be largely responsible for the ∼3.55 keV feature.
Swift intensive accretion disk reverberation mapping of four AGN yielded light curves sampled ∼200-350 times in 0.3-10 keV X-ray and six UV/optical bands. Uniform reduction and cross-correlation ...analysis of these data sets yields three main results: (1) The X-ray/UV correlations are much weaker than those within the UV/optical, posing severe problems for the lamp-post reprocessing model in which variations in a central X-ray corona drive and power those in the surrounding accretion disk. (2) The UV/optical interband lags are generally consistent with as predicted by the centrally illuminated thin accretion disk model. While the average interband lags are somewhat larger than predicted, these results alone are not inconsistent with the thin disk model given the large systematic uncertainties involved. (3) The one exception is the U band lags, which are on average a factor of ∼2.2 larger than predicted from the surrounding band data and fits. This excess appears to be due to diffuse continuum emission from the broad-line region (BLR). The precise mixing of disk and BLR components cannot be determined from these data alone. The lags in different AGN appear to scale with mass or luminosity. We also find that there are systematic differences between the uncertainties derived by JAVELIN versus more standard lag measurement techniques, with JAVELIN reporting smaller uncertainties by a factor of 2.5 on average. In order to be conservative only standard techniques were used in the analyses reported herein.
This work presents a new plasma cooling curve that is calculated using the SPEX package. We compare our cooling rates to those in previous works, and implement the new cooling function in the ...grid-adaptive framework “AMRVAC”. Contributions to the cooling rate by the individual elements are given, to allow for the creation of cooling curves tailored to specific abundance requirements. In some situations, it is important to be able to include radiative losses in the hydrodynamics. The enhanced compression ratio can trigger instabilities (such as the Vishniac thin-shell instability) that would otherwise be absent. For gas with temperatures below 104 K, the cooling time becomes very long and does not affect the gas on the timescales that are generally of interest for hydrodynamical simulations of circumstellar plasmas. However, above this temperature, a significant fraction of the elements is ionised, and the cooling rate increases by a factor 1000 relative to lower temperature plasmas.
The hot intra-cluster medium (ICM) permeating galaxy clusters and groups is not pristine, as it has been continuously enriched by metals synthesised in Type Ia (SNIa) and core-collapse (SNcc) ...supernovae since the major epoch of star formation (z ≃ 2–3). The cluster/group enrichment history and mechanisms responsible for releasing and mixing the metals can be probed via the radial distribution of SNIa and SNcc products within the ICM. In this paper, we use deep XMM-Newton/EPIC observations from a sample of 44 nearby cool-core galaxy clusters, groups, and ellipticals (CHEERS) to constrain the average radial O, Mg, Si, S, Ar, Ca, Fe, and Ni abundance profiles. The radial distributions of all these elements, averaged over a large sample for the first time, represent the best constrained profiles available currently. Specific attention is devoted to a proper modelling of the EPIC spectral components, and to other systematic uncertainties that may affect our results. We find an overall decrease of the Fe abundance with radius out to ~0.9 r500 and ~0.6 r500 for clusters and groups, respectively, in good agreement with predictions from the most recent hydrodynamical simulations. The average radial profiles of all the other elements (X) are also centrally peaked and, when rescaled to their average central X/Fe ratios, follow well the Fe profile out to at least ~0.5 r500. As predicted by recent simulations, we find that the relative contribution of SNIa (SNcc) to the total ICM enrichment is consistent with being uniform at all radii, both for clusters and groups using two sets of SNIa and SNcc yield models that reproduce the X/Fe abundance pattern in the core well. In addition to implying that the central metal peak is balanced between SNIa and SNcc, our results suggest that the enriching SNIa and SNcc products must share the same origin and that the delay between the bulk of the SNIa and SNcc explosions must be shorter than the timescale necessary to diffuse out the metals. Finally, we report an apparent abundance drop in the very core of 14 systems (~32% of the sample). Possible origins of these drops are discussed.
Four decades ago, the firm detection of an Fe-K emission feature in the X-ray spectrum of the Perseus cluster revealed the presence of iron in its hot intracluster medium (ICM). With more advanced ...missions successfully launched over the last 20 years, this discovery has been extended to many other metals and to the hot atmospheres of many other galaxy clusters, groups, and giant elliptical galaxies, as evidence that the elemental bricks of life—synthesized by stars and supernovae—are also found at the largest scales of the Universe. Because the ICM, emitting in X-rays, is in collisional ionisation equilibrium, its elemental abundances can in principle be accurately measured. These abundance measurements, in turn, are valuable to constrain the physics and environmental conditions of the Type Ia and core-collapse supernovae that exploded and enriched the ICM over the entire cluster volume. On the other hand, the spatial distribution of metals across the ICM constitutes a remarkable signature of the chemical history and evolution of clusters, groups, and ellipticals. Here, we summarise the most significant achievements in measuring elemental abundances in the ICM, from the very first attempts up to the era of
XMM-Newton
,
Chandra
, and
Suzaku
and the unprecedented results obtained by
Hitomi
. We also discuss the current systematic limitations of these measurements and how the future missions
XRISM
and
Athena
will further improve our current knowledge of the ICM enrichment.
In 2016 we carried out a Swift monitoring programme to track the X-ray hardness variability of eight type-I AGN over a year. The purpose of this monitoring was to find intense obscuration events in ...AGN, and thereby study them by triggering joint XMM-Newton, NuSTAR, and HST observations. We successfully accomplished this for NGC 3783 in December 2016. We found heavy X-ray absorption produced by an obscuring outflow in this AGN. As a result of this obscuration, interesting absorption features appear in the UV and X-ray spectra, which are not present in the previous epochs. Namely, the obscuration produces broad and blue-shifted UV absorption lines of Lyα, C iv, and N v, together with a new high-ionisation component producing Fe xxv and Fe xxvi absorption lines. In soft X-rays, only narrow emission lines stand out above the diminished continuum as they are not absorbed by the obscurer. Our analysis shows that the obscurer partially covers the central source with a column density of few 1023 cm-2, outflowing with a velocity of few thousand km s-1. The obscuration in NGC 3783 is variable and lasts for about a month. Unlike the commonly seen warm-absorber winds at pc-scale distances from the black hole, the eclipsing wind in NGC 3783 is located at about 10 light days. Our results suggest that the obscuration is produced by an inhomogeneous and clumpy medium, consistent with clouds in the base of a radiatively driven disk wind at the outer broad-line region of the AGN.