◀ Back

Escape fraction catalogs

For all snapshots we produce ionizing escape fraction catalogs directly mirroring the SUBFIND subhalo catalogs. These post-processing files provide supplemental data for each subhalo with at least one star particle and four gas cells within the virial radius, corresponding to identifications of dark matter halos. There is a single HDF5 file for each snapshot named ion_[file].hdf5 containing simulation attributes and subhalo datasets. Full details are given in Section 2 of Yeh et al. (2022).

We calculate the radius \(R_{\rm 200}\) at which the average density reaches 200 times the density of the universe, and adopt this same definition for the halo mass \(M_{\rm halo}\), stellar mass \(M_{\rm stars}\), instantaneous star-formation rate (SFR), and other quantities. This ensures a consistent accounting of properties contributing to the escape fraction, which we define as the ratio between the number of ionizing photons that escape and are intrinsically emitted from all sources within the virial radius, i.e. $$f_{\rm esc} \equiv \dot{N}_{\rm esc} / \dot{N}_{\rm int} \, .$$ The escape fractions for star particles outside the virial radius of any halo are assumed to be maximal (i.e. the birth cloud value \(f_{\rm esc}^{\rm cloud}\)) as these directly contribute to reionizing the IGM. The intrinsic photon rate for each subhalo (including stars outside \(R_{200}\) is provided in the Lyα emission catalog (Subhalo_NoSFR HDF5 group) described in Smith et al. (2022a), which is calculated from the Binary Population and Spectral Synthesis models (BPASS version 2.2.1; Eldridge et al. 2017) with a Chabrier initial mass function (IMF; Chabrier 2003). Passive Monte Carlo radiative transfer (MCRT) calculations were performed with the Cosmic Lyα Transfer (colt) code, with the details of the photo-ionization physics and MCRT implementation most recently described in Smith et al. (2022b).

We now provide the full list of ionizing escape fraction quantities. The main simulation attributes are provided in the HDF5 file as:
Attributes ▾
Attribute Dimensions Units Description
BoxSize 1 \({\rm ckpc} / h\) Spatial extent of the periodic box (in comoving units).
EscapeFraction 1 - Birth cloud escape fraction.
HubbleParam 1 \(100\,{\rm km/s/Mpc}\) Hubble constant (little \(h\) in standard units).
Ngroups_Total 1 - Total number of FoF groups in the entire catalog.
Nsubgroups_Total 1 - Number of SUBFIND subgroups in the entire catalog.
Omega0 1 - The cosmological density parameter for matter.
OmegaBaryon 1 - The cosmological density parameter for baryonic matter.
OmegaLambda 1 - The cosmological density parameter for the cosmological constant.
Redshift 1 - The redshift corresponding to the current output.
Time 1 - The scale factor \(a=1/(1+z)\) corresponding to the current output.
UnitLength_in_cm 1 - Code length units in cm.
UnitMass_in_g 1 - Code mass units in g.
UnitVelocity_in_cm_per_s 1 - Code velocity units in cm/s.
Subhalo properties are saved as HDF5 datasets:
Datasets ▾
Attribute Dimensions Units Description
L_int N \({\rm erg/s}\) Total luminosity of stars within \(R_{200}\).
M_200 N \({\rm M}_{\odot}\) Halo mass within \(R_{200}\).
M_gas N \({\rm M}_{\odot}\) Gas mass within \(R_{200}\).
M_stars N \({\rm M}_{\odot}\) Stellar mass within \(R_{200}\).
Ndot_int N \({\rm photons/s}\) Total intrinsic photon rate of stars within \(R_{200}\).
R_200 N \({\rm kpc}\) Radius of the halo at which its density reaches 200 times the critical density.
SFR_200 N \({\rm M}_{\odot} / {\rm yr}\) Sum of the instantaneous star-formation rates in gas cells within \(R_{200}\).
bin_L_int N \({\rm erg/s}\) Total luminosity of stars within \(R_{200}\), divided into three frequency bins \(({\rm HI, HeI, HeII})\) with energy edges at \((13.6, 24.59, 54.42, \infty)\,{\rm eV}\).
bin_Ndot_int N \({\rm photons/s}\) Total intrinsic photon rate of stars within \(R_{200}\), divided into three frequency bins as above.
bin_f_esc N - Halo escape fraction within \(R_{200}\), divided into three frequency bins as above.
f_HeI N - HeI absorption fraction within \(R_{200}\).
f_HeII N - HeII absorption fraction within \(R_{200}\).
f_esc N - Halo escape fraction within \(R_{200}\).
l_abs N \({\rm cm}\) Absorption distance: \(\langle \ell \rangle \equiv \ell^{-1} \int_0^\ell \ell' e^{-k_\text{a} \ell'} \text{d}\ell'\).
subhalo_id N - Subhalo ID in the SUBFIND subhalo catalog.

Note: \(N\) denotes the number of selected halos. The helium absorption fraction is simply the sum \(f_{\rm He} = f_{\rm HeI} + f_{\rm HeII}\) and the hydrogen absorption fraction can be calculated as \(f_{\rm H} = f_{\rm esc}^{\rm cloud} - f_{\rm esc} - f_{\rm He} - f_{\rm abs}\).