We give an overview of the SImulation of X-ray TElescopes (SIXTE) software package, a generic, mission-independent Monte Carlo simulation toolkit for X-ray astronomical instrumentation. The package is based on a modular approach for the source definition, the description of the optics, and the detector type such that new missions can be easily implemented. The targets to be simulated are stored in a flexible input format called SIMPUT. Based on this source definition, a sample of photons is produced and then propagated through the optics. In order to model the detection process, the software toolkit contains modules for various detector types, ranging from proportional counter and Si-based detectors, to more complex descriptions like transition edge sensor (TES) devices. The implementation of characteristic detector effects and a detailed modeling of the read-out process allow for representative simulations and therefore enable the analysis of characteristic features, such as for example pile-up, and their impact on observations. We present an overview of the implementation of SIXTE from the input source, the imaging, and the detection process, highlighting the modular approach taken by the SIXTE software package. In order to demonstrate the capabilities of the simulation software, we present a selection of representative applications, including the all-sky survey of eROSITA and a study of pile-up effects comparing the currently operating XMM-Newton with the planned Athena-WFI instrument. A simulation of a galaxy cluster with the Athena- X-IFU shows the capability of SIXTE to predict the expected performance of an observation for a complex source with a spatially varying spectrum and our current knowledge of the future instrument.

X-ray observations of galaxy clusters provide insights on the nature of gaseous turbulent motions, their physical scales and on the fundamental processes they are related to. Spatially-resolved, high-resolution spectral measurements of X-ray emission lines provide diagnostics on the nature of turbulent motions in emitting atmospheres. Since they are acting on scales comparable to the size of the objects, the uncertainty on these physical parameters is limited by the number of observational measurements, through sample variance. We propose a different and complementary approach for the computation of sample variance to repeating numerical simulations (i.e. Monte-Carlo sampling) by introducing new analytical developments for lines diagnosis. We consider the model of a "turbulent gas cloud", consisting in isotropic and uniform turbulence described by a universal Kolmogorov power-spectrum with random amplitudes and phases in an optically thin medium. Following a simple prescription for the 4-term correlation of Fourier coefficients, we derive generic expressions for the sample mean and variance of line centroid shift, line broadening and projected velocity structure function. We perform a numerical validation based on Monte-Carlo simulations for two popular models of gas emissivity based on the beta-model. Generic expressions for the sample variance of line centroid shifts and broadening in arbitrary apertures are derived and match the simulations within their range of applicability. Generic expressions for the mean and variance of the structure function are provided and verified against simulations. An application to the Athena/X-IFU and XRISM/Resolve instruments forecasts the potential of sensitive, spatially-resolved spectroscopy to probe the inertial range of turbulent velocity cascades in a Coma-like galaxy cluster.

X-ray observations of the hot gas filling the intra-cluster medium provide a wealth of information on the dynamics of clusters of galaxies. The global equilibrium of the ICM is believed to be partially ensured by non-thermal pressure support, notably the dissipation of energy through turbulent motions. Accurate mapping of turbulence using X-ray emission lines is challenging due to the lack of spatially-resolved spectroscopy. Only future instruments such as the X-ray Integral Field Unit (X-IFU) on Athena will have the spatial and spectral resolution to quantitatively investigate the ICM turbulence at all scales. Powerful diagnostics for these studies are line shift and the line broadening maps, and the second-order structure function. When estimating these quantities, instruments will be limited by uncertainties of their measurements, and by the sample variance (aka cosmic variance) of the observation. We extend here the formalism started in our companion paper I to include the effect of statistical uncertainties in the estimation of these line diagnostics, in particular for structure functions. We demonstrate that statistics contribute to the total variance through different terms, which depend on the geometry of the detector, the spatial binning and the nature of the turbulent field. These terms are important when probing the small scales of the turbulence. An application of these equations is performed for the X-IFU, using synthetic turbulent velocity maps of a Coma-like cluster of galaxies. Results are in excellent agreement with the formulas both for the structure function estimation (<3%) and its variance (<10%). The expressions derived here and in paper I are generic, and ensure an estimation of the total errors in any X-ray measurement of turbulent structure functions. They also open the way for optimisations in the upcoming instrumentation and in observational strategies.

Clusters of galaxies, the largest collapsed structures in the Universe, are located at the intersection of extended filaments of baryons and dark matter. Cosmological accretion onto clusters through large scale filaments adds material at cluster outskirts. Kinetic energy in the form of bulk motions and turbulence due to this accretion provides a form of pressure support against gravity, supplemental to thermal pressure. Significant amount of non-thermal pressure support could bias cluster masses derived assuming hydrostatic equilibrium, the primary proxy for cluster cosmology studies. Sensitive measurements of Doppler broadening and shift of astrophysical lines, and the relative fluctuations in thermodynamical quantities (e.g., density, pressure, and entropy) are primary diagnostic tools. Forthcoming planned and proposed X-ray (with large etendue, throughput, and high spectral resolution) and SZ observatories will provide crucial information on the assembly and virialisation processes of clusters, involving turbulent eddies cascading at various spatial scales and larger gas bulk motions in their external regions to the depth or their potential wells.

Answers to the metal production of the Universe can be found in galaxy clusters, notably within their Intra-Cluster Medium (ICM). The X-ray Integral Field Unit (X-IFU) on board the next-generation European X-ray observatory Athena (2030s) will provide the necessary leap forward in spatially-resolved spectroscopy required to disentangle the intricate mechanisms responsible for this chemical enrichment. In this paper, we investigate the future capabilities of the X-IFU in probing the hot gas within galaxy clusters. From a test sample of four clusters extracted from cosmological hydrodynamical simulations, we present comprehensive synthetic observations of these clusters at different redshifts (up to z = 2) and within the scaled radius R500 performed using the instrument simulator SIXTE. Through 100 ks exposures, we demonstrate that the X-IFU will provide spatially-resolved mapping of the ICM physical properties with little to no biases (<5%) and well within statistical uncertainties. The detailed study of abundance profiles and abundance ratios within R500 also highlights the power of the X-IFU in providing constraints on the various enrichment models. From synthetic observations out to z = 2, we also quantify its ability to track the chemical elements across cosmic time with excellent accuracy, and thereby to investigate the evolution of metal production mechanisms as well as the link to the stellar initial mass-function. Our study demonstrates the unprecedented capabilities of the X-IFU in unveiling the properties of the ICM but also stresses the data analysis challenges faced by future high-resolution X-ray missions such as Athena.

The X-ray Integral Field Unit (X-IFU) that will be on board the Athena telescope will provide an unprecedented view of the intracluster medium (ICM) kinematics through the observation of gas velocity, $v$, and velocity dispersion, $w$, via centroid-shift and broadening of emission lines, respectively. The improvement of data quality and quantity requires an assessment of the systematics associated with this new data analysis, namely biases, statistical and systematic errors, and possible correlations between the different measured quantities. We have developed an end-to-end X-IFU simulator that mimics a full X-ray spectral fitting analysis on a set of mock event lists, obtained using SIXTE. We have applied it to three hydrodynamical simulations of a Coma-like cluster that include the injection of turbulence. This allowed us to assess the ability of X-IFU to map five physical quantities in the cluster core: emission measure, temperature, metal abundance, velocity and velocity dispersion. Finally, starting from our measurements maps, we computed the 2D structure function (SF) of emission measure fluctuations, $v$ and $w$ and compared them with those derived directly from the simulations. All quantities match with the input projected values without bias; the systematic errors were below 5%, except for velocity dispersion whose error reaches about 15%. Moreover, all measurements prove to be statistically independent, indicating the robustness of the fitting method. Most importantly, we recover the slope of the SFs in the inertial regime with excellent accuracy, but we observe a systematic excess in the normalization of both SF$_v$ and SF$_w$ ascribed to the simplistic assumption of uniform and (bi-)Gaussian measurement errors. Our work highlights the excellent capabilities of Athena X-IFU in probing the thermodynamic and kinematic properties of the ICM. (abridged)

The X-ray Integral Field Unit (X-IFU) is the high resolution X-ray spectrometer of the ESA Athena X-ray observatory. Over a field of view of 5' equivalent diameter, it will deliver X-ray spectra from 0.2 to 12 keV with a spectral resolution of 2.5 eV up to 7 keV on ~5 arcsecond pixels. The X-IFU is based on a large format array of super-conducting molybdenum-gold Transition Edge Sensors cooled at about 90 mK, each coupled with an absorber made of gold and bismuth with a pitch of 249 microns. A cryogenic anti-coincidence detector located underneath the prime TES array enables the non X-ray background to be reduced. A bath temperature of about 50 mK is obtained by a series of mechanical coolers combining 15K Pulse Tubes, 4K and 2K Joule-Thomson coolers which pre-cool a sub Kelvin cooler made of a 3He sorption cooler coupled with an Adiabatic Demagnetization Refrigerator. Frequency domain multiplexing enables to read out 40 pixels in one single channel. A photon interacting with an absorber leads to a current pulse, amplified by the readout electronics and whose shape is reconstructed on board to recover its energy with high accuracy. The defocusing capability offered by the Athena movable mirror assembly enables the X-IFU to observe the brightest X-ray sources of the sky (up to Crab-like intensities) by spreading the telescope point spread function over hundreds of pixels. Thus the X-IFU delivers low pile-up, high throughput (>50%), and typically 10 eV spectral resolution at 1 Crab intensities, i.e. a factor of 10 or more better than Silicon based X-ray detectors. In this paper, the current X-IFU baseline is presented, together with an assessment of its anticipated performance in terms of spectral resolution, background, and count rate capability. The X-IFU baseline configuration will be subject to a preliminary requirement review that is scheduled at the end of 2018.

The Athena X-Ray Integral Field Unit (X-IFU) will provide spatially resolved high-resolution spectroscopy (2.5 eV FWHM up to 7 keV) over the 0.2 to 12 keV energy band. It will comprise an array of 3840 superconducting Transition Edge Sensors (TESs) operated at 90 mK, with an absolute energy scale accuracy of 0.4 eV. Slight changes in the TES operating environment can cause significant variations in its energy response function, which may result in degradation of the detector's energy resolution, and eventually in systematic errors in the absolute energy scale if not properly corrected. These changes will be monitored via an onboard Modulated X-ray Source (MXS) and the energy scale will be corrected accordingly using a multi-parameter interpolation of gain curves obtained during ground calibration. Assuming realistic MXS configurations and using the instrument end-to-end simulator SIXTE, we investigate here both statistical and systematic effects on the X-IFU energy scale, occurring either during ground measurements or in-flight. The corresponding impacts on the energy resolution and means of accounting for these errors are also addressed. We notably demonstrate that a multi-parameter gain correction, using both the pulse-height estimate and the baseline of a pulse, can accurately recover systematic effects on the gain due to realistic changes in TES operating conditions within 0.4 eV. Optimisations of this technique with respect to the MXS line configuration and correction time, as well as to the energy scale parametrization are also show promising results to improve the accuracy of the correction.

With its array of 3840 Transition Edge Sensors (TESs) operated at 90 mK, the X-Ray Integral Field Unit (X-IFU) on board the ESA L2 mission Athena will provide spatially resolved high-resolution spectroscopy (2.5 eV FWHM up to 7 keV) over the 0.2 to 12 keV bandpass. The in-flight performance of the X-IFU will be strongly affected by the calibration of the instrument. Uncertainties in the knowledge of the overall system, from the filter transmission to the energy scale, may introduce systematic errors in the data, which could potentially compromise science objectives - notably those involving line characterisation e.g. turbulence velocity measurements - if not properly accounted for. Defining and validating calibration requirements is therefore of paramount importance. In this paper, we put forward a simulation tool based on the most up-to-date configurations of the various subsystems (e.g. filters, detector absorbers) which allows us to estimate systematic errors related to uncertainties in the instrumental response. Notably, the effect of uncertainties in the energy resolution and of the instrumental quantum efficiency on X-IFU observations is assessed, by taking as a test case the measurements of the iron K complex in the hot gas surrounding clusters of galaxies. In-flight and ground calibration of the energy resolution and the quantum efficiency is also addressed. We demonstrate that provided an accurate calibration of the instrument, such effects should be low in both cases with respect to statistics during observations.

The X-ray Integral Field Unit (X-IFU) is the cryogenic imaging spectrometer on board the future X-ray observatory \textsl{Athena}. With a hexagonal array of 3840 AC-biased Transition Edge Sensors (TES), it will provide narrow-field observations (5$^{\prime}$ equivalent diameter) with unprecedented high spectral resolution (2.5 eV up to 7 keV) over the 0.2 - 12 keV bandpass. Throughout its observations, the X-IFU will face various sources of X-ray background. Specifically, the so-called Non-X-ray Background (NXB) caused by the interaction of high-energy cosmic rays with the instrument, may lead to a degradation of its sensitivity in the observation of faint extended sources (e.g. galaxy clusters outskirts). To limit this effect, a cryogenic anti-coincidence detector (CryoAC) will be placed below the detector plane to lower the NXB level down to the required level of $5 \times 10^{-3}$ cts/s/cm$^{2}$/keV over 2 - 10 keV. In this contribution, we investigate ways to accurately monitor the NXB and ensure the highest reproducibility in-flight. Using the limiting science case of the background-dominated observation of galaxy clusters outskirts, we demonstrate that a reproducibility of 2\% on the absolute knowledge of the background is required to perform driving science objectives, such as measuring abundances and turbulence in the outskirts. Monitoring of the NXB in-flight through closed observations, the detector's CryoAC or the companion instrument (Wide Field Imager) will be used to meet this requirement.

The X-ray Integral Field Unit (X-IFU) on board the Advanced Telescope for High-ENergy Astrophysics (Athena) will provide spatially resolved high-resolution X-ray spectroscopy from 0.2 to 12 keV, with 5 arc second pixels over a field of view of 5 arc minute equivalent diameter and a spectral resolution of 2.5 eV up to 7 keV. In this paper, we first review the core scientific objectives of Athena, driving the main performance parameters of the X-IFU, namely the spectral resolution, the field of view, the effective area, the count rate capabilities, the instrumental background. We also illustrate the breakthrough potential of the X-IFU for some observatory science goals. Then we briefly describe the X-IFU design as defined at the time of the mission consolidation review concluded in May 2016, and report on its predicted performance. Finally, we discuss some options to improve the instrument performance while not increasing its complexity and resource demands (e.g. count rate capability, spectral resolution). The X-IFU will be provided by an international consortium led by France, The Netherlands and Italy, with further ESA member state contributions from Belgium, Finland, Germany, Poland, Spain, Switzerland and two international partners from the United States and Japan.