Resonant scattering, metallicity, and redshift (Version 2)

Notes for the Data of Project Resonant Scattering as a function of redshift and metallicity (RSZZ), stored on PURR servers. v.02 The warmer gaseous component of the Large Scale Structure of the Universe emits primarily in the X-ray portion of the electromagnetic spectrum, mostly through thermal emission. Resonant scattering of X-ray background radiation by the Warm-Hot Intergalactic Medium (105 K < T < 107 K), however, can contribute significantly to the amplitude of the resonant emission lines of selected ions (with an increment of more than a factor 30). Tables for the emission (thermal and resonant scattering) at different temperature, density, metallicity, and redshift are generated with Cloudy. The comparison between thermal emission and resonant scattering contributions is performed with Jupyter Notebooks. All the material related to the published paper is stored in four sets: “Cloudy scripts”, “Jupyter scripts”, “Data”, and “Figures”. Additional material, that extends beyond the initial investigation, will be stored at a later time (the description will be updated accordingly). Cloudy Scripts. The Cloudy script “runcloudy” is the basic script to launch cloudy activities. The Cloudy script “runbatchcloudy” is called by “runcloudy” and runs all simulations. The Cloudy scripts “*.in” are called by runbatchcloudy and contain the instructions for each simulation. Input files ending with the tags "01HM", "02HM", "05HM", "HM12", and "KS19" are used to compare the different CXB models. The output is very large, each simulation generating aboutr 30 GB of data (not included in this repository). One simulation runs in about 1 hour. Jupyter scripts. The Jupyter notebook “Scattering_continuum_file_summary” reads the output of the Cloudy simulations and selects only two fields (energy bin and reflected emission) for the band [450-970 eV]. The output is stored in the “Data” set. This effectively reduced the ~30 GB of Cloudy data for one simulation to ~130 MB. The Jupyter notebook “Scattering_emission_maps_short” reads the output of “Scattering_continuum_file_summary” and plots the data for the selected ions (thermal emission, total emission, and ratios between total emission and thermal emission) as a function of temperature and density. The script also saves the raw data used for the plots. The Jupyter notebook “Scattering_emission_graphs” reads the output of “Scattering_emission_maps_short” and plots the data as function of density (or temperature) at selected metallicities and redshifts. The Jupyter notebook “Scattering_continuum_file_summary_bg” is equivalent to "Scattering_continuum_file_summary" but also handles the different CXB model cases. The Jupyter notebook “Scattering_emission_maps_short_bg” is equivalent to "Scattering_emission_maps_short" but also handles the different CXB model cases. The Jupyter notebook “Scattering_emission_graphs_bg” reads the output of “Scattering_emission_maps_short_bg”, it is equivalent to “Scattering_emission_graphs” but also handles the different CXB model cases. The Jupyter notebook “phase_diagrams” generates temperature-density phase diagrams of the output from the Illustris TNG data. The original input data is not provided and must be retrieved from the Illustris TNG simulation repositories. Data. The data folder does not include very large sets, like the Illustris TNG data and the Cloudy output. The Illustris dataset can be retrieved from the Illustris TNG Simulation repository. The Cloudy data can be generated with scripts in the “Cloudy script” set in less than 24 hours. The data folder contains the summary of the Cloudy data (labeled *_s.con) generated by “Scattering_emission_maps_short”. The summary is stored in a set of tar.gz files. If the original Cloudy data is not required, it is likely faster to download the *_s.con files from this repository. The Jupyter notebook copies data from the original continuum file (adopting the same units as the original); the first column (the energy of the spectrum bin) has units of eV, the second column (the “reflected” spectrum) has units of erg cm-2 s-1. The Figures data include the raw data generated by the "Scattering_emission_maps_short" Jupyter notebook. Figures. The "Maps" subset contains the phase diagrams generated by the “phase_diagrams” script (“fig_rho_T*.pdf”) and the emission and ratio maps generated by the script “Scattering_emission_maps_short” (“fig_ioname_*.pdf”). If the original Cloudy data or the summary data is not required, it is faster to download the figures from this repository. The phase diagrams have no units (the pixels show the mass fractions normalized to 1 over the full image). The emission maps have units of erg cm-2 s-1 Hz-1. The ratio maps have no units. The "Graphs" contains one-dimensional cuts of the emission and ratio maps as a function on density, at various temperatures, metallicities, redshifts, and CXB models. Graphs have the same units as the corresponding maps.