Model catalogues and histograms of KSVZ axion models with multiple heavy quarks

This record contains the files pertaining to the results mentioned in the linked work Plakkot & Hoof, “Anomaly Ratio Distributions of KSVZ Axion Models with Multiple Heavy Quarks.” The contents are:

4 histogram files

2 compressed tar.gz files containing detailed catalogues,

2 Python scripts to extract information from the files (requires Python 3 and the numpy, matplotlib and h5py packages)

Histogram files

The columns of the histogram files (listed below) correspond to the numerator of E, denominator of E, numerator of N, denominator of N, and the number of models (frequency) with that E/N ratio. The file histogram_complete_NQ_1_to_9.txt contains additional columns showing the number of models per entry for each NQ.

histogram_all_LP_allowed_models.txt for all LP-allowed models

histogram_additive_LP_allowed_models.txt for LP-allowed additive models

histogram_same_reps_LP_allowed_models.txt for LP-allowed additive models where all new quarks live on the same representation

histogram_complete_NQ_1_to_9.txt for all possible models with NQ ≤ 9, regardless of the LP criterion

Catalogues

The catalogues are contained in the two compressed HDF5 files listed below. The catalogues contain groups for different NQ, each with subgroups for numerators and denominators of E and N, model representation (as a list). The file catalogue_additive_models_NQ_1_to_9.tar.gz contains additionally the energy scale at which the first LP appears. The integers in the model lists represent the quark representations (integer m for the representation rm in the text), and the negative signs indicate opposite PQ charge.

catalogue_additive_models_NQ_1_to_9.tar.gz contains the catalogue for additive models with NQ ≤ 9 (9 groups, 6 subgroups)

catalogue_all_LP_allowed_models.tar.gz Contains catalogues for all LP-allowed models (28 groups, 5 subgroups)

Scripts

The two scripts are:

print_histogram_info.py is a sample script to extract information from the histogram files. The default choice for NQ can be adjusted by e.g. invoking python print_histogram_info.py 5 for NQ = 5

read_catalogue.py is a sample script to extract data from the catalogue files (which need to be unpacked first). The default settings can be overwritten by e.g. invoking python read_catalogue.py 5 file.h5 for NQ = 5 and the catalogues contained in file.h5