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.
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.
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.
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.
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.
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.
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.
Abstract
We present significantly improved measurements of turbulent velocities in the hot gaseous haloes of nearby giant elliptical galaxies. Using deep XMM–NewtonReflection Grating Spectrometer ...(RGS) observations and a combination of resonance scattering and direct line broadening methods, we obtain well bounded constraints for 13 galaxies. Assuming that the turbulence is isotropic, we obtain a best-fitting mean 1D turbulent velocity of ∼110 km s−1. This implies a typical 3D Mach number ∼0.45 and a typical non-thermal pressure contribution of ∼6 per cent in the cores of nearby massive galaxies. The intrinsic scatter around these values is modest – consistent with zero, albeit with large statistical uncertainty – hinting at a common and quasi-continuous mechanism sourcing the velocity structure in these objects. Using conservative estimates of the spatial scales associated with the observed turbulent motions, we find that turbulent heating can be sufficient to offset radiative cooling in the inner regions of these galaxies (<10 kpc, typically 2–3 kpc). The full potential of our analysis methods will be enabled by future X-ray micro-calorimeter observations.
Aims. We observed the Seyfert 1 galaxy NGC 985 on several occasions to search for variability in its UV and X-ray absorption features to establish their location and physical properties. Methods. We ...used XMM-Newton to obtain X-ray spectra using the EPIC-pn camera, and the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) to obtain UV spectra. Our observations were simultaneous and span timescales of days to years. Results. We find that the soft X-ray obscuration that absorbed the low energy continuum of NGC 985 in August 2013 diminished greatly by January 2015. The total X-ray column density decreased from 2.1 × 1022 cm-2 to ~6 × 1021 cm-2. We also detect broad, fast UV absorption lines in COS spectra obtained during the 2013 obscuration event. Lines of C iii*, Lyα, Si iv, and C iv with outflow velocities of −5970 km s-1 and a full-width at half-maximum of 1420 km s-1 are prominent in the 2013 spectrum, but have disappeared in all but Lyα in the 2015 spectra. The ionization state and the column density of the UV absorbing gas is compatible with arising in the same gas as that causing the X-ray obscuration. The high velocity of the UV-absorbing gas suggests that the X-ray obscurer and the associated UV outflow are manifestations of an accretion disk wind.
PDF
