Observing simulated galaxy clusters with PHOX: a novel X-ray photon simulator
We present a novel, virtual X–ray observatory designed to obtain synthetic observations from hydro–numerical simulations, named PHOX. In particular, we provide a description of the code constituting the photon simulator and of the new approach implemented. We apply PHOX to simulated galaxy clusters in order to demonstrate its capabilities. In fact, X–ray observations of clusters of galaxies continue to provide us with an increasingly detailed picture of their structure and of the underlying physical phenomena governing the gaseous component, which dominates their baryonic content. Therefore, it is fundamental to find the most direct and faithful way to compare such observational data with hydrodynamical simulations of cluster–like objects, which can currently include various complex physical processes. Here, we present and analyse synthetic Suzaku observations of two cluster–size haloes obtained by processing with PHOX the hydrodynamical simulation of the large–scale, filament–like region in which they reside. Taking advantage of the simulated data, we test the results inferred from the X–ray analysis of the mock observations against the underlying, known solution. Remarkably, we are able to recover the theoretical temperature distribution of the two haloes by means of the multi–temperature fitting of the synthetic spectra. Moreover, the shapes of the reconstructed distributions allow us to trace the different thermal structure that distinguishes the dynamical state of the two haloes.
keywords:hydrodynamics – methods: numerical – galaxies: clusters: general
As optimal laboratories for both cosmology and astrophysics, galaxy clusters have been thoroughly investigated with large surveys as well as via dedicated single–object observations. In this respect, X–ray measurements still represent one of the best ways to fully study their structure, traced by the hot plasma filling their potential wells. In particular, the most important and difficult quantity to infer is the total gravitating mass, which can be determined from measurements of the intra–cluster medium (ICM) temperature and density profiles, assuming hydrostatic equilibrium to hold. So far, studies of the X–ray–emitting ICM have been mainly limited to the innermost region, where the investigation of the baryonic physics is observationally less challenging than in the outskirts. Although very difficult to map in the X rays, a characterization of the outer region of clusters is in fact crucial to understand the formation and evolution of this structures and to use them as cosmological probes, since it encloses a significant fraction of the cluster volume, where the properties of the accreting gas and of the dark matter halo are still uncertain (e.g., see Ettori & Molendi, 2011). Only very recently, X–ray observations performed with the Japanese satellite Suzaku have reached the virial radius for a few clusters of galaxies (e.g. Fujita et al., 2008; Bautz et al., 2009; George et al., 2009; Reiprich et al., 2009; Hoshino et al., 2010; Kawaharada et al., 2010; Simionescu et al., 2011; Akamatsu et al., 2011), promisingly taking a step forward on the temperature profile debate and indirectly opening new perspectives in controlling the biases on cluster–based cosmological investigations. These results are just one example of the improvements by which X–ray observations continue to provide us with an increasingly detailed picture of galaxy clusters, for which a clear interpretation of the underlying physical processes is very challenging.
On the theoretical side, many interesting issues on cluster formation and evolution are addressed by means of hydro–numerical simulations, which represent a powerful tool to investigate in detail the nature of these complicated astrophysical objects (e.g. Borgani & Kravtsov, 2009, for a recent review). Simulations of large cosmological boxes that include several massive cluster–like haloes can now be performed incorporating not only the dominating dark matter component but also baryonic matter, in the form of gas– and star–like particles (or cells), governed by hydrodynamical processes. With particular concern for galaxy clusters, current simulations are not only able to account for basic gas hydrodynamics, but implement also more complicated models for star formation from multi–phase medium (e.g. Katz, 1992; Katz et al., 1996; Springel & Hernquist, 2003; Marri & White, 2003) and thermal or kinetic feedback from supernovae–driven winds Navarro & White (1993); Scannapieco et al. (2006); Dalla Vecchia & Schaye (2008), chemical enrichment, metal and molecule cooling (e.g. Mosconi et al., 2001; Yoshida et al., 2003; Tornatore et al., 2004, 2007; Scannapieco et al., 2005; Maio et al., 2007, 2010), thermal conduction Cleary & Monaghan (1999); Jubelgas et al. (2004); Ruszkowski et al. (2010), AGN feedback (e.g. Springel et al., 2005; Di Matteo et al., 2005; Sijacki & Springel, 2006; Sijacki et al., 2007, 2008; Puchwein et al., 2008; Fabjan et al., 2010; Dubois et al., 2010; Teyssier et al., 2011), cosmic rays Pfrommer et al. (2007); Jubelgas et al. (2008) and magnetic fields Phillips & Monaghan (1985); Dolag et al. (1999); Brüggen et al. (2005); Price & Monaghan (2005); Dolag & Stasyszyn (2009), to name the most important effects.
The ideal achievement would be the combination of both these sources of information, by directly comparing simulated clusters to X–ray observations of real objects. Devoted to this goal, sophisticated numerical codes, such as X–MAS/X–MAS2 Gardini et al. (2004); Rasia et al. (2008) and, more recently, XIM Heinz & Brüggen (2009); Heinz et al. (2010), have been developed in the last years in order to turn hydro–simulation outputs into mock images with a given X–ray telescope. In particular, the X–MAS virtual telescope is explicitly designed to process outputs obtained from SPH codes like GADGET, whereas XIM is particularly dedicated to grid–based hydrodynamical simulations, e.g. performed with the AMR code FLASH Fryxell et al. (2000), which makes them fundamentally complementary. Similarly, both X–MAS and XIM use a plasma thermal emission code by which they calculate the projected emission associated to the gas component in the simulation. The computational effort, required to calculate this, is usually reduced by interpolating among the model spectra externally stored in a library of templates. Essentially, the emissivity integrated along the chosen line of sight is calculated in terms of flux, depending on the properties of all the gas at a given projected position. Lastly, the convolution with instrumental response is performed and the final photon event file is generated. Such virtual telescopes represent an important step forward with respect to visualization tools used to produce surface brightness maps from simulation outputs, since they properly calculate the emission associated to the gas component accounting for the structure of the simulated source along the line of sight. Moreover, they are able to include accurate calculations required to obtain realistic mock X–ray observations, for instance convolving the simulated spectra with any given instrument response and telescope PSF. In fact, it is vital that the synthetic observations produced by such simulators are as much as possible similar to the standards commonly used by observers, so that the comparison can be the most faithful. As such, this advanced approach easily leads to a challenging increase of the parameters to handle in order to realistically model the plasma emission and the detailed three–dimensional structure of the simulations.
Analogous in the scope, the virtual telescope presented in this paper, PHOX, is also dedicated to convert hydro–numerical outputs into mock X–ray observations. However, the novelty of our simulator lies in the method adopted, by which the spectral emission calculated singularly for each gas element in the simulation is immediately converted into a discrete sample of photons, collected and stored before projecting along any line of sight and convolving with any desired instrument. With this strategy we are able to significantly gain in computational effort, since it requires to process the original simulation only once, independently of the specific synthetic observation to be performed afterwards. Furthermore, the guarantee for high spatial and spectral resolution, preserved without dramatically increasing the computational cost, offers the possibility to anticipate the observational achievements of upcoming X–ray missions, such as IXO/ATHENA, in which the high–resolution spectroscopy will allow us to explore the intrinsic structure of galaxy clusters through the study of their spectral features. In fact, it has been suggested in the last decade that the diagnostics of the broadening of heavy–ion emission lines in highly–resolved X–ray spectra of galaxy clusters could actually provide valuable information about the underlying structure of the ICM and its velocity field (see Pawl et al. (2005), for a preliminary, interesting work in this field), useful for the dynamical classification of these objects. The expectations for such line diagnostics are related in particular to the most prominent emission line in X–ray spectra, namely the line from helium–like iron. In fact, the large atomic mass of the FeXXV ion significantly reduces the thermal line broadening, so that the line width turns out to be definitely more sensitive to turbulent or bulk gas motions Inogamov & Sunyaev (2003); Sunyaev et al. (2003). This kind of studies can eventually help constraining the detectability of non–thermal motions that are likely to establish in the ICM (e.g. Rebusco et al., 2008) and can compromise X–ray mass measurements Rasia et al. (2006); Fang et al. (2009); Lau et al. (2009); Biffi et al. (2011).
The paper is structured as follows: Section 2 is devoted to the description of the method implemented in the photon simulator and the fundamental units of the code. In Section 3 we apply PHOX to the hydrodynamical simulation of a filament–like structure, describing in detail each phase, from the photon generation till the mock Suzaku observation of two massive galaxy clusters residing in the filament and the spectral analysis of the synthetic spectra. Results about the recovering of the ICM emission measure distribution of the two cluster–like haloes from the multi–temperature fitting are presented in Section 4. Finally, we discuss the results obtained, summarise and draw our conclusions in Section 5.
2 Phox: The Technique
The novelty of our photon simulator arises from the method implemented. In this Section we describe the three main parts of the code that constitute the virtual telescope, PHOX
With respect to the structure of the code, the first unit is the most general, essential one, in which the emission associated to the gaseous component in the simulation is generated without including any additional specification, neither about the observation to be simulated nor about the characteristics of the X–ray instrument, but the theoretical emission model. Ideally, the spectral model assumed can be very complicated and account for more realistic descriptions of the emitting gas which are included in up–to–date sophisticated hydro–simulations, such as metallicity, chemical composition and sub–grid turbulent broadening of the emission lines. Conveniently, this first unit can be executed independently from the others just once per simulation output, enormously reducing the computational effort and time. The second unit of the code takes into account the geometry of the problem: the projection along a given line of sight is applied and photon energies are corrected for the Doppler Shift. The order of first and second unit basically reverses the approach of publicly available X–ray virtual telescopes and this has also the main advantage of delaying any limitation on spectral resolution as much as possible, preserving the possibility for any single photon to contribute to the final spectrum.
2.1 Numerical hydro–simulation input
PHOX requires as input the output of a hydrodynamical simulation. For the snapshots of the simulations we use, the redshifts and cosmological parameters are provided as well as the main quantities characterising the gas elements, such as density, temperature, metallicity, position and velocity. Currently, the simulator is best–suited to process simulations performed with SPH, particle–based codes, in particular outputs of the N–body/SPH code GADGET–2 (e.g. Springel et al., 2001; Springel, 2005). Nevertheless, the approach is very general and the code can be easily adapted to process also grid–based simulations, for instance performed with Adaptive Mesh Refinement codes like ART Kravtsov (1999); Kravtsov et al. (2002), ENZO (Norman et al., 2007, and references therein) or FLASH Fryxell et al. (2000). In the description of the technique, we will preserve this generality and refer to the gaseous component in terms of more general emitting volumes, namely “gas elements”, meaning that these could be either particles, as in the application presented in the following Sections 3 and 4, or grid cells.
2.2 Unit 1: generation of the box of photons
The first, essential unit is devoted to convert the simulation snapshot from a box of gas elements into a box of photons. It is crucial to perform this step before the projection is chosen and the emission is integrated along the line of sight in order to achieve and guarantee a very high spatial and spectral precision, accounting for the details of the three–dimensional structure of the source.
Given the density, temperature and metallicity (or abundances), the spectrum associated with each gas element of the simulation input is generated by means of the publicly available XSPEC package (v.12)
The main calculation of this first unit consists of computing the emission associated to the gas component processing element per element, which even makes the code naturally suitable to parallelization. Our simulator directly populates the spectrum of each gas element with a distribution of photons, without storing every spectrum or building a library of template spectra. Specifically, we calculate a cumulative spectrum and perform a Monte–Carlo–like simulation to sample it with a total number of photons, , determined by the total counts expected from the spectral model (in units of ) and the collecting area () and observing time () specified by the user, so that
Ideally, the fiducial values for () and () have to be assumed in a convenient way such that the spectra are largely over–sampled.
The photons are therefore characterised by energy (in the emitting element rest–frame), position and velocity. In order to reduce the amount of memory required to store the photon data, photons are organised in terms of packages, each package being associated to an emitting gas element. In this way, photon energies can be stored in a separate output file with respect to position and velocity, which in this first version of the code are equal for all the photons emitted by the same gas element and therefore stored just once per package. Moreover, the data are therefore naturally compressed, since positions and velocities need to be recorded only for gas elements that indeed emit photons. The total number of photons produced from the whole gaseous component in the snapshot processed by PHOX is meant to be largely over–abundant in order to permit a dramatic reduction in the following units because of geometrical selection, projection and instrumental response.
The data generated from this basic part of PHOX represent a “cube of photons” associated to the input hydro–simulation. From this stage on, the simulator has all the information required to proceed with the synthetic observation and the original input is no longer needed. Therefore, it is worth to remark that this first unit can be processed once per output independently of the specific study to be performed.
2.3 Unit 2: projection and preparation for the observation
PHOX second unit takes into account the geometry of the observation to be simulated. It accepts the output data produced by the first unit and a parameter file with user–specified parameters for the sub–region of the cube that has to be selected, direction corresponding to the line of sight (l.o.s.), real collecting area and time for the mock observation.
After selecting geometrically the photons, the energies have to be corrected for the Doppler Shift in order to account for the line–of–sight velocity of the gas element that originally emitted them. In particular, the correction from the emitted to the observed frame is expressed as
where and is the velocity component along the line of sight. For the sake of simplicity, the line of sight is assumed to be aligned with the –axis. Nevertheless, any different, desired direction only requires a simple rotation of the package positions and velocities before running this second unit of the code.
The sample of photons obtained has also to be further filtered according to the specified collecting area, and a realistic exposure time, chosen to simulate the observation. Accordingly to Eq. 1, this is done by calculating the re–scaling factor between these values and the fiducial quantities assumed in the first unit as
and by assuming it as a probability factor for each photon to be actually observed. Eq. 3 explicitly accounts for the possibility of re–sacling the number of observable photons by the angular–diameter distance to the source, whenever this parameter differs from the fiducial value, adopted in Unit 1.
Lastly, the photon list has to be stored in the most convenient format either to be convolved with the technical response of the specific instrument or to interface with an external software dedicated to simulate specific X–ray satellites.
This second unit already takes into account the characteristics of the study to be performed but no longer requires the original simulation data, processed by Unit 1. Therefore, any change in the specification of the mock observation can be easily included in this post–processing phase. As such, a straightforward consequence is the unique possibility to investigate the same astrophysical source from many different line of sights, handily processing with Unit 2 the cube of virtual photons obtained with the first unit.
2.4 Unit 3: simulating the observation
At the third and last stage the mock X–ray observation is completed by eventually considering a real telescope. The photon list obtained from the second unit is convolved with the technical characteristics of a specific instrument, defined by the redistribution matrix file (RMF) and the ancillary response file (ARF). Such process produces final event files which satisfy the standards of real X–ray observations, so that they can be analysed with the same procedures and tools used for real data.
In the most general perspective this last unit assumes the RMF and ARF files supplied by the user and it is not constructed a priori for a specific instrument. Particular attention, in general, has to be payed to the normalization of the effective area defined by the ARF file of the instrument with respect to the effective area assumed during the projection phase, in order to avoid unphysical overabundance of observed photons.
The third unit of the simulator, being independent of the others, can be conveniently replaced by any desired X–ray instrument simulator, like the sophisticated tool XISSIM
3 SYNTHETIC OBSERVATION OF THE COSMIC WEB: two cases of study
We present here an optimal application of PHOX to the hydro–numerical simulation of a filament–like structure. The large–scale region provides indeed a study case to test the remarkable gain in computational cost and time for the post–processing of large data sets; also, the cluster–like haloes residing in this simulated filament offer various science cases to perform the spectral analysis of the mock X–ray observations.
3.1 The simulated region
The high–resolution simulated region containing the supercluster–like structure extends for about and consists of 27 haloes at present epoch, four of them being massive cluster–size haloes. This filamentary region has been originally extracted from a cosmological, N–body simulation Jenkins et al. (2001); Yoshida et al. (2001) of a CDM universe with , , and , within a box of a side. Using the “zoomed initial condition” (ZIC) technique Tormen et al. (1997), the overdense region has been re–simulated at higher resolution and a hydrodynamical run with these new initial conditions has been performed with the TreePM/SPH code GADGET–2 Springel et al. (2001); Springel (2005), including cooling, star formation and feedback from supernova winds. The simulation output contains approximately gas particles and DM particles, at , and the final resolution is and for the DM and gas particles, respectively. For an extensive and detailed description of the simulation run and analysis we refer the reader to the paper by Dolag et al. (2006).
3.2 Suzaku mock observations with PHOX
We applied PHOX to the whole simulation output corresponding to redshift and obtained from the first unit a list of all the photons emitted by the sufficiently–hot gas particles in the simulated region. In order to compute the model spectra, we adopted an APEC model with fixed metallicity, set for simplicity to for all the particles, and assumed the solar abundances by Anders & Grevesse (1989). The value assumed for the Galactic equivalent column density for the WABS absorption model was The values for fiducial collecting area and exposure time were initially set to and , respectively. This first stage was the most computationally demanding (see Table 1 for time scales and memory requirements referring to the application presented here) and permitted to generate approximately photons, for the whole simulation output.
|Unit 1||Unit 2||Unit 3/XISSIM|
|1h 40min||few sec||few min|
|1.8 Gb||827 Mb||176 Mb (CL2)||100 Mb (CL2)|
|60 Mb (CL1)||35 Mb (CL1)|
The run reported here has been performed by binning in temperature the particles (assuming ) and calculating a spectrum per bin, rescaling afterwards the expected total number of photons by the specific normalization calculated for each particle in the temperature bin. Since no varying metallicities have been considered in this run, this approach has been used to further reduce the computational effort.
The results are obtained for serial runs performed on standard work station (2300 MHz, AMD Opteron).
In Unit 2, we assumed to observe the large–scale region from a physical distance
In the left–hand panel of Fig. 2 we show the filament–like structure at , projected along the –axis (i.e. the observation l.o.s.). The map represents the emission–weighted temperature of the simulated filament enclosed in the high–resolution region. Overplotted in green, we mark the zoom onto the selected haloes, which are shown in the central column of Fig. 2. The two panels contain in fact the maps of the emission–weighted temperature projected along the l.o.s., for CL1 and CL2. For comparison, we overplot also the region (cyan circle) and the XIS field of view (magenta dashed square). Instead, in the right–hand column of Fig. 2 we show the ideal integrated spectra for CL1 and CL2, as obtained from Unit 1 and 2 after selecting the photons coming from the region enclosed within . These spectra are obtained purely through a regular binning of the photon energies with bins of , without including any sensitivity or effective area.
For the purpose of this paper, we present PHOX by taking advantage of the public package XISSIM Ishisaki et al. (2007), designed to simulate observations with the XIS spectrometer on board the Suzaku X–ray satellite. By adapting the second unit of our code to produce the output photon list in a format best–suited for XISSIM, we test our novel technique and demonstrate its capabilities through a simple science application.
Specifically, we use XISSIM to process the simulated photons generated by PHOX with the real responses and calibration characteristics of different XIS CCDs. Given the high statistics offered by our simulation, we decide to consider the back–illuminated (BI: XIS-1) sensor and the two combined, front–illuminated CCDs (FI: XIS-0, -3) separately. The combination of two observations with the XIS-0 and XIS-3 detectors additionally improves the statistics for the final FI spectrum.
The event files obtained for the two cluster–like haloes of the filament, CL1 and CL2, without any addition of physical background emission, were then analysed following the standards of X–ray data analysis.
3.3 Spectral Analysis
In Fig. 3 we show the simulated photon images (left–hand column) of the XIS observations, with the BI and FI sensor, for CL1 and CL2. The images, as well as the spectra, were extracted from the original event files using Xselect (v.2.4) from the FTOOLS
The normalization of an APEC component is defined as where is the angular–diameter distance to the source.
Simulations of cluster–like haloes provide the possibility of precisely knowing the intrinsic dynamical and thermal properties of the ICM. From a visual inspection of the temperature maps shown in Fig. 2 (central column), the different structure of the two haloes studied is already evident. As we will discuss in more detail in Section 4 (see Fig. 4), the further investigation of the emission measure (hereafter, E.M.) distribution as function of temperature, for the gas particles in the simulation that reside in the FOV of Suzaku both in CL1 and CL2, let us unveil the halo intrinsic thermal structures. Especially, CL1 is mainly dominated by one temperature component, although the E.M. distribution is not very narrow, whereas the second halo, CL2, clearly has a complex thermal structure, which cannot be well described by a single temperature.
Therefore, we fit a multi–temperature model to the BI and FI XIS spectra, restricting our analysis to the energy band and requiring a minimum of
For our spectral modelling of the simulated data, we considered a model, describing each component by an APEC model, as assumed in Unit 1 for the emission associated to each gas element. The equivalent hydrogen column density for the Galactic absorption model WABS was frozen to as in the original run of Unit 1. For analogous reasons, we adopted fixed values for redshift, and metallicity always assuming the Anders & Grevesse (1989) solar abundances.
In both haloes, we assigned to the first temperature component, a value fairly close to the emission–weighted temperature estimated from the gas particles in the simulation. Then, following a strategy similar to what is suggested by Peterson et al. (2003) or Kaastra et al. (2004), we assigned to the other components temperatures of and for the cooler ones. The five normalizations of the APEC components where free parameters in the fit.
The spectra and best–fit models are shown in the right–hand column of Fig. 3, for the two haloes and the two XIS sensor kinds. The results of the spectral fit are presented in Table 2. The five temperatures, frozen in the fit, are different for CL1 and CL2 and therefore reported in the table for clarity reasons. Overall, the resulting is always very good, meaning that the five–temperature model provides a fair description of the thermal structure for both clusters in the central observed region.
4 Results: recovering the emission measure distribution
The interesting purpose of our multi–temperature modelling is the reconstruction of the thermal structure of the emitting ICM. This is pursued by directly comparing the distribution of the best–fit normalizations at the corresponding, fixed temperatures with the expected values calculated from the simulation. Dealing with simulated data, we can faithfully test the results of the spectral analysis against the underlying solution in a straightforward way, which is not possible for real data.
By definition, the normalization of each APEC component is
where is the angular–diameter distance to the source and is the redshift. This implies that the normalization is directly proportional to the E.M. of the hot emitting plasma, defined as
where and are the electron and hydrogen number densities of the plasma, respectively. In order to map the distribution of E.M. in the simulated haloes to the corresponding values at fixed temperatures, the theoretical temperature binning has to be determined in the most faithful way possible.
The shape of the E.M. distribution extracted from the simulation is presented for the two haloes in the left–hand panels of Fig. 4. This has been calculated considering all the gas particles of the simulated haloes, for which the projected position on the plane is enclosed in the Suzaku field of view. A lower temperature cut has been artificially set to which is the lower boundary of the energy band used in the spectral fit. Nevertheless, this continuous distribution has to be reduced to the five fixed temperatures adopted for the spectral fit. The most consistent and independent way of choosing the boundaries between the different temperature components, is by requiring that the value of the E.M.–weighted temperature for each bin equals the fixed value adopted for the fitting (we follow here an approach analogous to the one adopted by Kaastra et al. (2004)). Specifically, we iteratively adjust the division boundaries until the match between the E.M.-weigthed temperature in the bin and the fixed value is reached, starting from the lower–temperature component and proceeding towards the higher ones. Finally, the total E.M. calculated in each bin can be converted into an expected value for the APEC normalization via Eq. 4.
The direct comparison between simulation and mock–observation results is presented in the right–hand panels of Fig. 4, for both CL1 and CL2. The original E.M. distribution for the two haloes is plotted in the left–hand panels of Fig. 4, where it is clear that none of them is strictly isothermal, although the CL1 cluster has a narrower thermal distribution with respect to CL2. The theoretical expectations for the fixed temperature components are plotted as black asterisks on the right–hand panels, where the best–fit values are reported as well.
Interestingly, the agreement between simulated and “observed” values is quite good and the overall shape of the thermal distribution in both clusters is recovered reasonably well. This result is found to be true for both the BI and FI XIS sensors. Moreover, the differences between the thermal structures of the two haloes are reflected by a different width and extension of the E.M. distribution over temperature, which is definitely broader in the CL2 case.
We notice that the dominant component at is underestimated by the fit for the halo CL1.
Nonetheless, it is worth to notice that the determination of the most consistent temperature binning of the theoretical E.M. distribution is expected to be sensitive to the definition of temperature assumed to match the temperature of the thermal component in the spectral fit. For instance, it is known that there is a discrepancy between the emission–weighted temperature, as usually calculated in the analysis of numerical simulations, and the projected spectroscopic temperature obtained from X–ray observations (see, for instance, Mazzotta et al., 2004). In this respect, our choice for the theoretical temperature calculated from the simulation, which is weighted by the E.M., is believed to be the most consistent.
However, we also tested the results against a different definition of temperature, namely the emission–weighted
In the right–hand plots of Fig. 4 we also show two alternative five–temperature models adopted to fit the synthetic spectra. Specifically, we perform two spectral fits with a multi–temperature model consisting of five BAPEC components. The BAPEC spectral model is a velocity– and thermally–broadened emission spectrum, where the velocity broadening is parametrized by a gaussian velocity dispersion . In fact, the technique implemented in our simulator allows to preserve the very high energy resolution of the integrated spectrum till the final convolution with the instrument. Since every photon energy carries the imprint of the emitting particle velocity, the spectral features like metal emission lines could include a non-thermal broadening due to peculiar non–thermal motions of the gas particles of the ICM and the velocity–broadened APEC is a valid choice to take this into account. In the two cases explored with the model, every frozen parameter was assumed in the same manner as in Section 3.3 and the five BAPEC normalizations were let be free. In addition, in the first model (BAPEC–1), we fit for another free parameter, which is the velocity dispersion linked across the five components, whereas in the second case (BAPEC–2) both and the redshift are left free in the fit. As shown in the plots, the results for the normalizations are definitely consistent with the standard case, where the five thermal components are described by the APEC model. However, the best–fit values for in particular, turn out to have very large errors, meaning that this additional parameter is not reliably fitted. We interpret this as an evidence that the energy resolution of the Suzaku XIS spectrometer is not sufficiently high to capture the details of the emission lines and therefore a velocity–broadened model like BAPEC is not improving the fit to the synthetic data.
Moreover, the results look generally stable against metallicity, which could in principle compete with normalization in order to reproduce the spectral features associated to metal emission lines. In particular, we tested the fit with five APEC models, both for CL1 and CL2, allowing for the metallicity value to be free, although linked across the five components, in addition to the five normalizations. Generally, this converges to values within few percents from the expected metallicity of used in the simulation of the X–ray observation. Comparing to the results of the standard fit reported in Table 2, the APEC best–fit normalizations, for initial metallicity values set to are found to be consistent with the original values, within the error bars associated to them.
5 Discussion and Conclusions
In this paper, we have presented a novel X–ray photon simulator, named PHOX
The development of PHOX has been strongly motivated by the ultimate aim of comparing the output of hydrodynamical simulations of cluster–like objects to real X–ray observations in the most faithful and reliable way possible. In fact, X–ray observations still provide us with increasingly detailed data that require a deep understanding of the underlying structure of galaxy clusters. Numerical simulations, on the other side, continuously improve in describing not only the dark matter but also the complex physical processes governing the baryonic component. Therefore, a strategy devoted to combine the advances and results in both observations and simulations can be promisingly successful.
In the past, software packages like X–MAS Gardini et al. (2004) or XIM Heinz & Brüggen (2009) have been developed to produce synthetic observations of astrophysical sources simulated with particle– and grid–based codes, respectively. The existence of such tools is definitely complementary to our simulator in order to test the robustness of comparison studies between state–of–the–art hydrodynamical simulations including different physical phenomena and real observations. With respect to other codes, the novelty of our virtual telescope resides in the original method implemented. In particular, the order of the three main units that constitute the code offers a significant gain in computational cost and time, making PHOX extremely flexible and useful to process hydro–numerical simulations of several astrophysical objects, at various spatial scales, as well as simulated light–cones. Moreover, this method uniquely guarantees extremely high spatial and energy resolution. The more demanding phase is indeed the first one, where the output of a hydrodynamical simulation is processed and the X–ray emission of each gas element (either particle or grid cell) is simulated by using the spectral models available in the X–ray software XSPEC. Essentially, we perform the sampling of each spectrum with an expected number of photons, and we immediately store the sample of photons insted of constructing a library of template spectra to be populated afterwards. As such, the spectral model chosen to mock up the emission from the gas elements, regardless of its complexity, can be explicitly specified and suited to represent any desired astrophysical source. Also, the relatively reasonable requirements in terms of computational time and memory, make out of PHOX an ideal tool to process forthcoming large cosmological, hydrodynamical simulations as well as the usual smaller ones.
Concerning galaxy clusters, PHOX currently uses two possible models for the X–ray emission of the ICM (namely, MEKAL and APEC) and accounts for the particular density and temperature of the simulated gas elements. However, the technique implemented offers the chance to account also for different metallicities, or even metal compositions, across the emitting gas elements, without increasing the memory requirements significantly.
Furthermore, PHOX postpones the projection along the line of sight and the convolution with a real instrument response till after the generation of the discrete sample of photons, permitting to each of them, in principle, to contribute to the final integrated spectrum. Requiring a high energy resolution in the model spectrum associated to every gas element, we can ideally combine the potential of PHOX with the characteristics of upcoming X–ray satellites, like IXO/ATHENA or ASTRO–H, in order to predict the detectability of fundamental spectral features. This could be of great help, for instance, in constraining the ICM non–thermal motions (e.g. turbulent, bulk, rotational motions) via the velocity broadening of heavy–ion emission lines. In general, the theoretical modelling of X–ray–emitting sources is a powerful way to predict and optimize the capabilities of future X–ray missions and to help in the interpretation of the data when the instruments are in operation.
As an interesting application of this new virtual telescope, we focused on the hydro–numerical simulation of a large–scale region, originally performed with the TreePM/SPH code GADGET–2 and then re–simulated including cooling, star formation and feedback from supernova winds. This filament–like structure is particularly worthy for presenting the simulator in each part and especially for demonstrating the advantages of our approach. From a purely technical perspective, the greatest gain was achieved in terms of computational cost since the entire, large data cube contained in the simulation output was processed with PHOX Unit 1 only once, providing an output of relatively compacted photon data (see Table 1). In fact, this first phase permitted to obtain the X–ray photons virtually emitted by all the gas particles in the filament. At a later stage, we restricted our analysis to two cluster–like haloes with different thermodynamical structure, CL1 and CL2, which were geometrically selected from two different regions of the simulation box. Conveniently, the projection and observation units of PHOX could be run separately for each halo processing the same output generated with Unit 1 and we finally obtained synthetic Suzaku images and spectra (Fig. 3). The convolution with the response of the X–ray Imaging Spectrometer (XIS) on board Suzaku, via the public tool XISSIM, provided data that could be analysed with the standard procedures of X–ray analysis. An important aim of this application was the reconstruction of the intrinsic emission measure distribution by the multi–temperature fitting of the simulated spectra. The goal of setting spectroscopic constraints on the multi–phase ICM of galaxy clusters was pursued in X–ray observational studies reported by Peterson et al. (2003), Kaastra et al. (2004) and Fabian et al. (2006), although for real clusters it is intrinsically impossible to know the underlying multi–phase structure unambiguously. Instead, numerical simulations uniquely allow to model and predict the theoretical temperature distribution of the ICM (e.g., studied in Kawahara et al., 2007, 2008) and constitute therefore the ideal possibility to test the power of the method used to reconstruct the structure of the mock data. The direct comparison between the X–ray data and the theoretical thermal distribution of the input simulation is found to be overall faithful and robust for both CL1 and CL2 (see Fig. 4). Moreover, the different thermal structure that distinguishes the two selected haloes is remarkably highlighted by the multi–temperature fitting. This result promisingly suggests that the temperature distribution of the ICM can be in principle traced and unveiled by spectral analysis of X–ray observations. From the additional investigation of velocity–broadened fitting models (like the BAPEC model) for the Suzaku spectra, we can also conclude that a higher spectral resolution is required in order to reliably constrain the ICM velocity structure. In general, it is important to notice that the application discussed here does not include other emission components (e.g. a physical X–ray background) in addition to the ICM one and represents therefore an ideal case study.
With regard to technique, the main feature distinguishing PHOX from other virtual X-ray telescopes, such as XMAS and XIM, is the explicit creation and storage of a 3–dimensional box of photons from the raw simulation data cube itself. This step is performed before and without any reference to the observation process and its result is stored as a separate data product, which represents the ideal X–ray emission from all the gas volumes in the simulated region, allowing us to employ Unit 1 of PHOX to create a photon database associated to the simulation. In Unit 2 this photon cube is projected onto a particular observational plane and the photons are sampled according to the chosen aperture and observation time, as described in Eq. 3. With the photon box pre–calculated, this step is computationally trivial and can be performed very efficiently. As such, it could in fact be made available through an on line web service. The result of this step is a photon list and is processed in this paper by means of the Suzaku simulator XISSIM. Similarly, we believe it should be possible for most such virtual telescopes to be adapted to the same format without undue difficulties.
This leads us to an interesting modular design, with a photon simulator producing a list
of photons directly from a hydrodynamical simulation, and a virtual telescope receiving this list,
and convolving it with the characteristics of a real X–ray instrument to produce a photon event file that can
be directly compared in all its aspects to the analogous files resulting from real observations.
This design presents a proper separation of concerns between theorists and observers, where each module is clearly independent from the others and can therefore be implemented and used separately.
In a larger perspective, the well specified model for the data shared and exchanged between these units could represent
a first step towards a standardised, interoperable way to publish
simulation results in a way that is suitable for investigation by observers.
In particular, this is very suitable for the kind of approach pursued by the International Virtual Observatory Alliance
The authors gratefully thank Eugene Churazov, Hajime Kawahara, Stefano Ettori, Elena Rasia, Gayoung Chon and Umberto Maio for useful discussions, that helped improving this work. V.B. acknowledges support from and participation in the International Max Planck Research School on Astrophysics at the Ludwig Maximilian University. K.D. acknowledges the support by the DFG Priority Programme 1177 and additional support by the DFG Cluster of Excellence ”Origin and Structure of the Universe”. The work of G.L. was supported by Advanced Grant 246797 ”GALFORMOD” from the European Research Council.
- pagerange: Observing simulated galaxy clusters with PHOX: a novel X–ray photon simulator–References
- pubyear: …
- The source code of the presented version of PHOX is made available by request to the authors (see http://www.mpa-garching.mpg.de/ kdolag/Phox/). We also plan to make the code available directly for anonymous download in the near future.
- See http://www.gnu.org/s/gsl/.
- See http://heasarc.gsfc.nasa.gov/fitsio/.
- See http://heasarc.gsfc.nasa.gov/xanadu/xspec/.
- See http://heasarc.nasa.gov/docs/suzaku/prop_tools/xissim/
- We set the observer on the positive part of the –axis, assumed to be aligned with the l.o.s..
- See http://heasarc.gsfc.nasa.gov/ftools/.
- The spectra were rebinned using the grppha routine from the FTOOLS package.
- The emission–weighted temperature is defined as where is the gas density and is the cooling function. As usually done, we consider here assuming the gas to emit mainly via Bremsstrahlung.
- The source code of the presented version of PHOX is made available by request to the authors (see http://www.mpa-garching.mpg.de/ kdolag/Phox/). We also plan to make the code available directly for anonymous download in the near future.
- See http://www.ivoa.net.
- Please check http://www.mpa-garching.mpg.de/HydroSims for updates on this web service.
- Akamatsu H., Hoshino A., Ishisaki Y., Ohashi T., Sato K., Takei Y., Ota N., 2011, ArXiv e-prints
- Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta, 53, 197
- Arnaud K. A., 1996, in G. H. Jacoby & J. Barnes ed., Astronomical Data Analysis Software and Systems V Vol. 101 of Astronomical Society of the Pacific Conference Series, XSPEC: The First Ten Years. pp 17–+
- Bautz M. W., Miller E. D., Sanders J. S., Arnaud K. A., Mushotzky R. F., Porter F. S., Hayashida K., Henry J. P., Hughes J. P., Kawaharada M., Makashima K., Sato M., Tamura T., 2009, Pub. Astron. Soc. Japan, 61, 1117
- Biffi V., Dolag K., Böhringer H., 2011, MNRAS, 413, 573
- Blackburn J. K., 1995, in R. A. Shaw, H. E. Payne, & J. J. E. Hayes ed., Astronomical Data Analysis Software and Systems IV Vol. 77 of Astronomical Society of the Pacific Conference Series, FTOOLS: A FITS Data Processing and Analysis Software Package. pp 367–+
- Borgani S., Kravtsov A., 2009, ArXiv e-prints
- Brüggen M., Ruszkowski M., Simionescu A., Hoeft M., Dalla Vecchia C., 2005, ApJ, 631, L21
- Cleary P. W., Monaghan J. J., 1999, Journal of Computational Physics, 148, 227
- Dalla Vecchia C., Schaye J., 2008, MNRAS, 387, 1431
- Di Matteo T., Springel V., Hernquist L., 2005, Nature, 433, 604
- Dolag K., Bartelmann M., Lesch H., 1999, A&A, 348, 351
- Dolag K., Meneghetti M., Moscardini L., Rasia E., Bonaldi A., 2006, MNRAS, 370, 656
- Dolag K., Stasyszyn F., 2009, MNRAS, 398, 1678
- Dubois Y., Devriendt J., Slyz A., Teyssier R., 2010, MNRAS, 409, 985
- Ettori S., Molendi S., 2011, Memorie della Societa Astronomica Italiana Supplementi, 17, 47
- Fabian A. C., Sanders J. S., Taylor G. B., Allen S. W., Crawford C. S., Johnstone R. M., Iwasawa K., 2006, MNRAS, 366, 417
- Fabjan D., Borgani S., Tornatore L., Saro A., Murante G., Dolag K., 2010, MNRAS, 401, 1670
- Fang T., Humphrey P., Buote D., 2009, ApJ, 691, 1648
- Fryxell B., Olson K., Ricker P., Timmes F. X., Zingale M., Lamb D. Q., MacNeice P., Rosner R., Truran J. W., Tufo H., 2000, ApJS, 131, 273
- Fujita Y., Tawa N., Hayashida K., Takizawa M., Matsumoto H., Okabe N., Reiprich T. H., 2008, Pub. Astron. Soc. Japan, 60, 343
- Gardini A., Rasia E., Mazzotta P., Tormen G., De Grandi S., Moscardini L., 2004, MNRAS, 351, 505
- George M. R., Fabian A. C., Sanders J. S., Young A. J., Russell H. R., 2009, MNRAS, 395, 657
- Heinz S., Brüggen M., 2009, ArXiv e-prints
- Heinz S., Brüggen M., Morsony B., 2010, ApJ, 708, 462
- Hoshino A., Henry J. P., Sato K., Akamatsu H., Yokota W., Sasaki S., Ishisaki Y., Ohashi T., Bautz M., Fukazawa Y., Kawano N., Furuzawa A., Hayashida K., Tawa N., Hughes J. P., Kokubun M., Tamura T., 2010, Pub. Astron. Soc. Japan, 62, 371
- Inogamov N. A., Sunyaev R. A., 2003, Astronomy Letters, 29, 791
- Ishisaki Y., Maeda Y., Fujimoto R., Ozaki M., Ebisawa K., Takahashi T., Ueda Y., Ogasaka Y., Ptak A., 25 co–authors 2007, Pub. Astron. Soc. Japan, 59, 113
- Jenkins A., Frenk C. S., White S. D. M., Colberg J. M., Cole S., Evrard A. E., Couchman H. M. P., Yoshida N., 2001, MNRAS, 321, 372
- Jubelgas M., Springel V., Dolag K., 2004, MNRAS, 351, 423
- Jubelgas M., Springel V., Enßlin T., Pfrommer C., 2008, A&A, 481, 33
- Kaastra J. S., Mewe R., 1993, A&A Supp., 97, 443
- Kaastra J. S., Tamura T., Peterson J. R., Bleeker J. A. M., Ferrigno C., Kahn S. M., Paerels F. B. S., Piffaretti R., Branduardi-Raymont G., Böhringer H., 2004, A&A, 413, 415
- Katz N., 1992, ApJ, 391, 502
- Katz N., Weinberg D. H., Hernquist L., 1996, ApJS, 105, 19
- Kawahara H., Reese E. D., Kitayama T., Sasaki S., Suto Y., 2008, ApJ, 687, 936
- Kawahara H., Suto Y., Kitayama T., Sasaki S., Shimizu M., Rasia E., Dolag K., 2007, ApJ, 659, 257
- Kawaharada M., Okabe N., Umetsu K., Takizawa M., Matsushita K., Fukazawa Y., Hamana T., Miyazaki S., Nakazawa K., Ohashi T., 2010, ApJ, 714, 423
- Kravtsov A. V., 1999, PhD thesis, NEW MEXICO STATE UNIVERSITY
- Kravtsov A. V., Klypin A., Hoffman Y., 2002, ApJ, 571, 563
- Lau E. T., Kravtsov A. V., Nagai D., 2009, ApJ, 705, 1129
- Liedahl D. A., Osterheld A. L., Goldstein W. H., 1995, ApJ, 438, L115
- Maio U., Ciardi B., Dolag K., Tornatore L., Khochfar S., 2010, MNRAS, 407, 1003
- Maio U., Dolag K., Ciardi B., Tornatore L., 2007, MNRAS, 379, 963
- Marri S., White S. D. M., 2003, MNRAS, 345, 561
- Mazzotta P., Rasia E., Moscardini L., Tormen G., 2004, MNRAS, 354, 10
- Mewe R., Gronenschild E. H. B. M., van den Oord G. H. J., 1985, A&A Supp., 62, 197
- Morrison R., McCammon D., 1983, ApJ, 270, 119
- Mosconi M. B., Tissera P. B., Lambas D. G., Cora S. A., 2001, MNRAS, 325, 34
- Navarro J. F., White S. D. M., 1993, MNRAS, 265, 271
- Norman M. L., Bryan G. L., Harkness R., Bordner J., Reynolds D., O’Shea B., Wagner R., 2007, ArXiv e-prints
- Pawl A., Evrard A. E., Dupke R. A., 2005, ApJ, 631, 773
- Peterson J. R., Kahn S. M., Paerels F. B. S., Kaastra J. S., Tamura T., Bleeker J. A. M., Ferrigno C., Jernigan J. G., 2003, ApJ, 590, 207
- Pfrommer C., Enßlin T. A., Springel V., Jubelgas M., Dolag K., 2007, MNRAS, 378, 385
- Phillips G. J., Monaghan J. J., 1985, MNRAS, 216, 883
- Price D. J., Monaghan J. J., 2005, MNRAS, 364, 384
- Puchwein E., Sijacki D., Springel V., 2008, ApJ, 687, L53
- Rasia E., Ettori S., Moscardini L., Mazzotta P., Borgani S., Dolag K., Tormen G., Cheng L. M., Diaferio A., 2006, MNRAS, 369, 2013
- Rasia E., Mazzotta P., Bourdin H., Borgani S., Tornatore L., Ettori S., Dolag K., Moscardini L., 2008, ApJ, 674, 728
- Rebusco P., Churazov E., Sunyaev R., Böhringer H., Forman W., 2008, MNRAS, 384, 1511
- Reiprich T. H., Hudson D. S., Zhang Y.-Y., Sato K., Ishisaki Y., Hoshino A., Ohashi T., Ota N., Fujita Y., 2009, A&A, 501, 899
- Ruszkowski M., Lee D., Bruggen M., Parrish I., Oh S. P., 2010, ArXiv e-prints
- Scannapieco C., Tissera P. B., White S. D. M., Springel V., 2005, MNRAS, 364, 552
- Scannapieco C., Tissera P. B., White S. D. M., Springel V., 2006, MNRAS, 371, 1125
- Sijacki D., Pfrommer C., Springel V., Enßlin T. A., 2008, MNRAS, 387, 1403
- Sijacki D., Springel V., 2006, MNRAS, 366, 397
- Sijacki D., Springel V., Di Matteo T., Hernquist L., 2007, MNRAS, 380, 877
- Simionescu A., Allen S. W., Mantz A., Werner N., Takei Y., Morris R. G., Fabian A. C., Sanders J. S., Nulsen P. E. J., George M. R., Taylor G. B., 2011, Science, 331, 1576
- Smith R. K., Brickhouse N. S., Liedahl D. A., Raymond J. C., 2001, ApJ, 556, L91
- Springel V., 2005, MNRAS, 364, 1105
- Springel V., Di Matteo T., Hernquist L., 2005, MNRAS, 361, 776
- Springel V., Hernquist L., 2003, MNRAS, 339, 289
- Springel V., White M., Hernquist L., 2001, ApJ, 549, 681
- Sunyaev R. A., Norman M. L., Bryan G. L., 2003, Astronomy Letters, 29, 783
- Teyssier R., Moore B., Martizzi D., Dubois Y., Mayer L., 2011, MNRAS, 414, 195
- Tormen G., Bouchet F. R., White S. D. M., 1997, MNRAS, 286, 865
- Tornatore L., Borgani S., Dolag K., Matteucci F., 2007, MNRAS, 382, 1050
- Tornatore L., Borgani S., Matteucci F., Recchi S., Tozzi P., 2004, MNRAS, 349, L19
- Yoshida N., Abel T., Hernquist L., Sugiyama N., 2003, ApJ, 592, 645
- Yoshida N., Colberg J., White S. D. M., Evrard A. E., MacFarland T. J., Couchman H. M. P., Jenkins A., Frenk C. S., Pearce F. R., Efstathiou G., Peacock J. A., Thomas P. A., 2001, MNRAS, 325, 803