Output File#

Introduction#

The output file is and HDF5 file and contains the output of a calculation. Overlap with variables in the input file is minimized, i.e., if a variable is present in the input file it may not be written to the output file.

Below we discuss all of the variables in each output group hierarchically (groups are in bold, datasets are in code). Similar to the input file, if the Dim of a variable is [ : ], this means the array has variable dimension, and may change from calculation to calculation.

Note

To further understand the output file, go through this document with the output file open in HDFView, a graphical user interface to HDF5 files: download

Warning

Fortran is a column-major language, which affects reading/writing of multi-dimensional datasets with HDF5. For example, an \(n \times m\) matrix, written to an HDF5 file with Fortran, will be read in to a row-major language (e.g., python) as an \(m \times n\) matrix. Dimensions listed below are the variables dimensions when written with EXCEED-DM

General Groups#

exdm_version - Version number of EXCEED-DM.
dm_model
- FIF_id - String ID of the scattering form factor used.
- mX - Dark matter masses, \(m_\chi\).
  - Units: \(\text{eV}\)
  - Dim: [ : ]
- med_FF - Mediator form factor powers, \(\beta\).
  - Formula: \(\mathcal{F}_\text{med} = \left( \frac{q_0}{q} \right)^\beta\)
  - Dim: [ : ]
- particle_type - Type of DM particle
- rho_X - DM density, \(\rho_\chi\)
  - Units: \(\text{GeV}/\text{cm}^3\)
experiment
- M - Mass of the target.
  - Units: \(\text{kg}\)
- T - Exposure time of the target.
  - Units: \(\text{year}\)
material
- band_gap - Band gap of the target material.
  - Units: \(\text{eV}\)
- name - Name of the target material.
- pc_vol - Primitive cell/unit cell volume.
  - Units: \(\text{Å}^3\)
- rho_T - Target density.
  - Units: \(\text{g}/\text{cm}^3\)
timing
- dt_compute - Time taken to perform sums in the calculation.
  - Units: \(\text{s}\)
- dt_total - Total run time of the calculation.
  - Units: \(\text{s}\)
- start_date - Start year/month/day of the calculation.
  - Dim: [3]
- end_date - End year/month/day of the calculation.
  - Dim: [3]

Calculation Specific Groups#

Binned Scatter Rate#

binned_scatter_rate
- model_<n>/v_e_<v>/mass_<m> - Folder structure for a binned_scatter_rate calculation. Note that if only a single DM model, or Earth velocity vector, is specified these folders will be absent. For example, if only a light mediator with one Earth velocity vector is chosen the folder structure will be binned_scatter_rate/mass_<m> for each mass. However if a heavy and light mediator are calculated the folder structure will be binned_scatter_rate/model_<n>/mass_<m>. Similarly for mulitple Earth velocity vectors.
  - total_binned_scatter_rate - Total number of events in each \(q, \omega\) bin, assuming \(\overline{\sigma}_e = 1\) and an exposure of M_kg x T_year (kg-yr), where M_kg, T_year are specified in the input file.
    
    Note
    
    Note the absence of units on \(\overline{\sigma}_e\) above. To compute the total number of events simply multiply the data by \(\overline{\sigma}_e\) in units of \(\text{cm}^2\). For example to compute the number of events with \(\overline{\sigma}_e = 10^{-40} \, \text{cm}^2\), mulitply the data by \(10^{-40}\).
    
    The event rate, in units of events/kg-year, is then calculated by dividing the entries by the exposure in kg-year. For example, if M_kg = 1, T_year = 1 (the default) in the experiment group in the input file, the entries are both, a) the total number of events in a kg-year and b) the number of events/kg-year. If M_kg = 10, T_year = 1, the entries are the total number of events in 10 kg-year, and the event rate is easily calculated (simply divide by 10 in this example).
    - Dim: [ :, : ], [\(N_q\), \(N_\omega\)]
    - Units: \(\text{cm}^{-2}\)
  - i_<i>/binned_scatter_rate - Total number of events in each \(q, \omega\) bin from each initial state group labelled by i, e.g., the band number for Bloch states. The sum of the entries in each i sum to the total_binned_scatter_rate.
    - Dim: [ :, : ], [\(N_q\), \(N_\omega\)]
    - Units: \(\text{cm}^{-2}\)
numerics_binned_scatter_rate
- E_bin_width - Width of the bins in \(\omega\) space.
  - Units: \(\text{eV}\)
- q_bin_width - Width of the bins in \(q\) space.
  - Units: \(\text{keV}\)

Absorption Rate#

absorption_rate
- width_<i>/mass_<m> - Folder structure for an absorption_rate calculation. Note that if only a single width is specified these folders will be absent. For example, if only a single width is chosen the folder structure will be absorption_rate/mass_<m> for each mass. However if multiple widths are calculated the folder structure will be binned_scatter_rate/width_<i>/mass_<m>.
  - absorption_rate - Total number of events assuming \(g_e = 1\), and an exposure of M_kg x T_year (kg-yr). The event rate, in units of events/kg-year, is then calculated by dividing the entries by the exposure in kg-year. For example, if M_kg = 1, T_year = 1 (the default) in the experiment group in the input file, the entries are both, a) the total number of events in a kg-year and b) the number of events/kg-year. If M_kg = 10, T_year = 1, the entries are the total number of events in 10 kg-year, and the event rate is easily calculated (simply divide by 10 in this example).
numerics_absorption_rate
- widths - List of width parameterizations computed for, [\(a\), \(b\), \(c\)].
  - Formula: \(\delta = \text{min}(c, a + b \omega)\)
  - Units: [ : , { \(\text{eV}\), - , \(\text{eV}\) } ]
  - Dim: [ : , 3]
- smear_type - Defines broadening behavior for the imaginary part of the Greens function.

Dielectric#

dielectric
- width_<i> - Folder containing the dielectric for width parameterization i
  - dielectric_r - Real part of the averaged dielectric, \(\overline{\varepsilon}(\mathbf{q}, \omega)\).
    - Dim: [ : , : , : , : ], [\(N_q\), \(N_\theta\), \(N_\phi\), \(N_\omega\)]
  - dielectric_c - Imaginary part of the averaged dielectric, \(\overline{\varepsilon}(\mathbf{q}, \omega)\).
    - Dim: [ : , : , : , : ], [\(N_q\), \(N_\theta\), \(N_\phi\), \(N_\omega\)]
numerics_dielectric
- E_bin_width - Width of the bins in \(\omega\) space.
  - Units: \(\text{eV}\)
- q_bin_width - Width of the bins in \(q\) space.
  - Units: \(\text{keV}\)
- widths - List of width parameterizations computed for, [\(a\), \(b\), \(c\)].
  - Formula: \(\delta = \text{min}(c, a + b \omega)\)
  - Units: [ : , { \(\text{eV}\), - , \(\text{eV}\) } ]
  - Dim: [ : , 3]
- smear_type - Defines broadening behavior for the imaginary part of the Greens function.