DOI: 10.5281/zenodo.2615235
title: Supplementary Data: Axion global fits with Peccei-Quinn symmetry breaking before inflation using GAMBIT

description: <strong>Description of Supplementary Data</strong>

This record contains the samples used to create the figures (excluding validation and prior dependence plots) and to derive most of the results in Hoof et al., <em>“Axion global fits with Peccei-Quinn symmetry breaking before inflation using GAMBIT”</em> (available on the arXiv). Please contact the authors if you are interested in other samples, YAML files or plotting scripts.<br>
<br>
This record consists of


	21 <code>YAML</code> files (6 for <code>T-Walk</code>, 15 for <code>Diver</code>). Running <code>./gambit -f path/to/YAML/file.yaml</code> in the GAMBIT directory will start the scan. However, most users might want to adjust the output file name and directory as well as the settings for the samplers to their systems.
	21 <code>hdf5</code> files (6 for <code>T-Walk</code>, 15 for <code>Diver</code>). These files contain the actual samples and were compressed using the <code>tar</code> format.
	Two example <code>pip</code> files (<code>2_QCDAxion_10M1.pip</code> for <code>Diver</code> samples, <code>2_QCDAxion_3041.pip</code> for <code>T-Walk</code> samples) for producing plots from the corresponding <code>hdf5</code> files, using <code>pippi</code> and <code>functions.py</code>.


The files follow the naming scheme <code>V_ModelName_[S][C][I][R][E]</code> plus one of the extensions <code>.yaml</code>, <code>.hdf5.tar.gz</code>, or <code>.pip</code>.


	<code>V</code>: This internal version number can be ignored, but should be quoted when asking for help with the plotting scripts
	<code>ModelName</code>: Corresponds to the axion models in the paper (<em>GeneralALP</em>, <em>QCDAxion</em>, <em>DFSZAxion_I</em>, <em>DFSZAxion_II</em>, <em>KSVZAxion</em>)
	<code>S</code>: Scanner (<code>S=1</code>: <code>Diver</code>, <code>S=3</code>: <code>T-Walk</code>)
	<code>C</code>: Switch to include (<code>C=1</code>) or exclude (<code>C=0</code>) the White Dwarf cooling hints
	<code>I</code>: Setting for the initial misalignment angle <em>θ<sub>i</sub></em> (<code>I=4</code>: flat prior on <em>θ<sub>i</sub></em> with values in [-3.1415, 3.1415]). <code>I=M</code> is used to indicate that the file includes merged samples from other scans in addition to the corresponding <code>I=4</code> scan.
	<code>R</code>: Setting for the DM relic density likelihood (<code>R=1</code>: upper limit, <code>R=2</code>: matching the DM density)
	<code>E</code>: Extra digit for the anomaly ratio <em>E/N</em>; only for <em>KSVZAxion</em> models (<code>E=1</code>: 0, <code>E=2</code>, 2/3, <code>E=3</code>: 5/3, <code>E=4</code>: 8/3), <em>DFSZAxion-I</em> models (<code>E=1</code>: 8/3), <em>DFSZAxion-II</em> models (<code>E=2</code>: 2/3), or some <em>GeneralALP</em> files (<code>E=a</code>: “QCD-like setting” with <em>β</em> = 7.94, <em>T<sub>crit</sub></em> = 147 MeV; <code>E=b</code>: “Simple ALP-like setting” with <em>β</em> = 0, <em>T<sub>crit</sub></em> irrelevant)


For convenience, we provide a mapping between the figures in the paper and the <code>hdf5</code> files:


	Fig. 1: none
	Figs 2 - 11: Validation plots
	Figs 12 + 13: 2_GeneralALP_10M2
	Fig. 14: 2_GeneralALP_10M2a, 2_GeneralALP_10M2b
	Fig. 15: 2_QCDAxion_10M1, 2_QCDAxion_10M2
	Fig. 16: 2_QCDAxion_3041, 2_QCDAxion_3042
	Figs 17 + 18: 2_QCDAxion_10M1, 2_QCDAxion_10M2, 2_QCDAxion_30M1, 2_QCDAxion_30M2
	Fig. 19: 3_KSVZAxion_10M11, 3_KSVZAxion_10M12, 3_KSVZAxion_10M13, 3_KSVZAxion_10M14, 3_DFSZAxion_I_10M11, 3_DFSZAxion_II_10M12
	Fig. 20: 2_QCDAxion_10M1, 3_KSVZAxion_10M11, 3_KSVZAxion_10M12, 3_KSVZAxion_10M13, 3_KSVZAxion_10M14, 3_DFSZAxion_I_10M11, 3_DFSZAxion_II_10M12
	Fig. 21: 2_QCDAxion_11M1, 2_QCDAxion_11M2
	Fig. 22: 2_QCDAxion_3141, 2_QCDAxion_3142
	Figs 23 + 24: 2_QCDAxion_3041, 2_QCDAxion_3042, 2_QCDAxion_3141, 2_QCDAxion_3142
	Fig. 25: 2_QCDAxion_3141, 2_QCDAxion_3142
	Fig. 26: 2_QCDAxion_11M1, 3_DFSZAxion_I_11M11, 3_DFSZAxion_II_11M12
	Fig. 27: 2_QCDAxion_3141, 3_DFSZAxion_I_31411, 3_DFSZAxion_II_31412
	Fig. 28: Validation plot
	Fig. 29: Prior dependence plot


A few caveats to keep in mind:


	The YAML files are designed to work with <code>GAMBIT 1.3.0</code>, and the pip files are tested with <code>pippi 2.1</code>, commit 1a08644. They may or may not work with later versions of either software (these working versions/commits can always be obtained via the <code>git</code> history).
	The <code>pip</code> files will produce an approximately complete, but very basic version of plots in the paper. Re-creating all the plots in the paper requires various manual, undocumented interventions such as additions, deletions and combination of the plotting scripts created by <code>pippi</code>. Users wishing to reproduce the more advanced plots in the paper should contact the authors for tips, scripts, or experiment for themselves.
description: The updated version v2 adjusts this Zenodo record to match the peer-reviewed version of the article (and v2 of the arXiv preprint). In particular, the units of the axion-photon coupling are now inverse GeV and some of the parameter names have changed compared to version v1.