We present the design and implementation of a thermal model, developed in COMSOL, aiming to probe the wafer-scale thermal response arising from realistic rates and energies of cosmic rays at L2 impacting the detector wafer of Athena X-IFU. The wafer thermal model is a four-layer 2D model, where 2 layers represent the constituent materials (Si bulk and Si$_{3}$N$_{4}$ membrane), and 2 layers represent the Au metallization layer's phonon and electron temperatures. We base the simulation geometry on the current specifications for the X-IFU detector wafer, and simulate cosmic ray impacts using a simple power injection into the Si bulk. We measure the temperature at the point of the instrument's most central TES detector. By probing the response of the system and pulse characteristics as a function of the thermal input energy and location, we reconstruct cosmic ray pulses in Python. By utilizing this code, along with the results of the GEANT4 simulations produced for X-IFU, we produce realistic time-ordered data (TOD) of the temperature seen by the central TES, which we use to simulate the degradation of the energy resolution of the instrument in space-like conditions on this wafer. We find a degradation to the energy resolution of 7 keV X-rays of $\approx$0.04 eV. By modifying wafer parameters and comparing the simulated TOD, this study is a valuable tool for probing design changes on the thermal background seen by the detectors.

While kilohertz quasi-periodic oscillations (kHz QPOs) have been well studied for decades since their initial discovery, the cause of these signals remains unknown, as no model has been able to accurately predict all of their spectral and timing properties. Separately, X-ray reverberation lags have been detected in AGN and stellar-mass black hole binaries, and reverberation may be expected to occur in neutron star systems as well, producing lags of the same amplitude as the lags measured of the kHz QPOs. Furthermore, the detection of a relativistically reflected Fe K line in the time-averaged spectra of many neutron star systems provides an additional motivation for testing reverberation. While it has been shown that the lag-energy properties of the lower kHz QPOs are unlikely to be produced by X-ray reverberation, the upper kHz QPOs have not yet been explored. We therefore model the upper kHz QPO lag-energy spectra using relativistic ray-tracing functions and apply them to archival RXTE data on 4U 1728-34 where upper kHz QPOs have been detected. By modeling the time-averaged spectra in which upper kHz QPOs had been significantly detected, we determine the reflected flux fraction across all energies and produce a model for the lag-energy spectra from X-ray reverberation. We explore the dependence of the modeled lag properties on several different types of reflection models, but are unable to successfully reproduce the measured lags of 4U 1728-34. We conclude that reverberation alone does not explain the measured time lags detected in upper kHz QPOs.

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.

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) 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.

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.

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.

Kilohertz quasi-periodic oscillations or kHz QPOs occur on the orbital timescale of the inner accretion flow and may carry signatures of the physics of strong gravity (c$^{2}$ ~ GM/R) and possibly clues to constraining the neutron star equation of state (EOS). Both the timing behavior of kHz QPOs and the time-averaged spectra of these systems have been studied extensively, yet no model completely describes all the properties of kHz QPOs. Here, we present a systematic study of spectral-timing products of kHz QPOs from low-mass X-ray binary systems using archival Rossi X-ray Timing Explorer/Proportional Counter Array data. For the lower kHz QPOs in fourteen objects and the upper kHz QPOs in six, we were able to obtain correlated time-lags as a function of QPO frequency and energy, as well as energy-dependent covariance spectra and intrinsic coherence. For the lower kHz QPOs, we find a monotonic decrease in lags with increasing energy, rising covariance to ~12 keV, and near unity coherence at all energies. For the upper kHz QPOs, we find near zero lags, rising covariance to ~12 keV, and less well-constrained coherence at all energies. These results suggest that while kHz QPOs are likely produced by similar mechanisms across the population of LMXBs, the lower kHz QPOs are likely produced by a different mechanism than upper kHz QPOs

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.

Soft lags from the emission of the lower kilohertz quasi-periodic oscillations (kHz QPOs) of neutron star low mass X-ray binaries have been reported from 4U1608-522 and 4U1636-536. Those lags hold prospects for constraining the origin of the QPO emission. In this paper, we investigate the spectral-timing properties of both the lower and upper kHz QPOs from the neutron star binary 4U1728-34, using the entire Rossi X-ray Timing Explorer archive on this source. We show that the lag-energy spectra of the two QPOs are systematically different: while the lower kHz QPO shows soft lags, the upper kHz QPO shows either a flat lag-energy spectrum or hard variations lagging softer variations. This suggests two different QPO-generation mechanisms. We also performed the first spectral deconvolution of the covariance spectra of both kHz QPOs. The QPO spectra are consistent with Comptonized blackbody emission, similar to the one found in the time-averaged spectrum, but with a higher seed-photon temperature, suggesting that a more compact inner region of the Comptonization layer (boundary/spreading layer, corona) is responsible for the QPO emission. Considering our results together with other recent findings, this leads us to the hypothesis that the lower kHz QPO signal is generated by coherent oscillations of the compact boundary layer region itself. The upper kHz QPO signal may then be linked to less-coherent accretion-rate variations produced in the inner accretion disk, being detected when they reach the boundary layer.

The intense radiation flux of Type I X-ray bursts is expected to interact with the accretion flow around neutron stars. High frequency quasiperiodic oscillations (kHz QPOs), observed at frequencies matching orbital frequencies at tens of gravitational radii, offer a unique probe of the innermost disk regions. In this paper, we follow the lower kHz QPOs, in response to Type I X-ray bursts, in two prototypical QPO sources, namely 4U 1636-536 and 4U 1608-522, as observed by the Proportional Counter Array of the Rossi X-ray Timing Explorer. We have selected a sample of 15 bursts for which the kHz QPO frequency can be tracked on timescales commensurable with the burst durations (tens of seconds). We find evidence that the QPOs are affected for over ~200 s during one exceptionally long burst and ~100 s during two others (although at a less significant level), while the burst emission has already decayed to a level that would enable the pre-burst QPO to be detected. On the other hand, for most of our burst-kHz QPO sample, we show that the QPO is detected as soon as the statistics allow and in the best cases, we are able to set an upper limit of ~20 s on the recovery time of the QPO. This diversity of behavior cannot be related to differences in burst peak luminosity. We discuss these results in the framework of recent findings that accretion onto the neutron star may be enhanced during Type I X-ray bursts. The subsequent disk depletion could explain the disappearance of the QPO for ~100 s, as possibly observed in two events. However, alternative scenarios would have to be invoked for explaining the short recovery timescales inferred from most bursts. Clearly the combination of fast timing and spectral information of Type I X-ray bursts holds great potential in the study of the dynamics of the inner accretion flow around neutron stars.

The Athena+ mission concept is designed to implement the Hot and Energetic Universe science theme submitted to the European Space Agency in response to the call for White Papers for the definition of the L2 and L3 missions of its science program. The Athena+ science payload consists of a large aperture high angular resolution X-ray optics and twelve meters away, two interchangeable focal plane instruments: the X-ray Integral Field Unit (X-IFU) and the Wide Field Imager (WFI). The X-IFU is a cryogenic X-ray spectrometer, based on a large array of Transition Edge Sensors (TES), offering 2.5 eV spectral resolution, with ~5" pixels, over a field of view of 5 arc minutes in diameter. In this paper, we briefly describe the Athena+ mission concept and the X-IFU performance requirements. We then present the X-IFU detector and readout electronics principles, the current design of the focal plane assembly, the cooling chain and review the global architecture design. Finally, we describe the current performance estimates, in terms of effective area, particle background rejection, count rate capability and velocity measurements. Finally, we emphasize on the latest technology developments concerning TES array fabrication, spectral resolution and readout performance achieved to show that significant progresses are being accomplished towards the demanding X-IFU requirements.