A wrapper code around ACROPOLIS (https://github.com/hep-mh/acropolis/) to calculate photo- and hadrodisintegration constraints for more specific BSM scenarios that are not part of the core ACROPOLIS package
| Input | Output |
|---|---|
| param_file.dat | bbundance_file.dat |
| cosmo_file.dat | |
| abundance_file.dat |
First install v2(currently in alpha) of ACROPOLIS (preferably in a virtual environment) via
pip3 install git+https://github.com/hep-mh/acropolis.git@v2
Afterwards run
./bin/xena <model_class> <io_directory>
Here, both command-line arguments are mandatory. The argument <model_class> is the name of the model to use (currently only NeutrinoCascadeModel and EmDecayModel are supported), and the argument <io_directory> is the directory where to find cosmo_file.dat, param_file.dat and abundance_file.dat
An example output of running ./bin/xena NeutrinoCascadeModel io/test would be
Xena v1.0 (using ACROPOLIS v2.0.0 [dev])
INFO : Using the model 'NeutrinoCascadeModel'.
INFO : Using 'universal = True' as required by the model.
INFO : Using the directory 'io/test' to read and write data.
INFO : Using parameters from 'io/test/param_file.dat'.
INFO : Using cosmological data from 'io/test/cosmo_file.dat'.
INFO : Using initial abundances from 'io/test/abundance_file.dat'.
INFO : Running photo- and hadrodisintegration...Done!
INFO : The final abundances are:
| mean | high Γ | low Γ | high ξ | low ξ
-------------------------------------------------------------------------------
n | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 [decayed]
p | 7.559015e-01 | 7.559705e-01 | 7.558327e-01 | 7.559047e-01 | 7.558983e-01
H2 | 2.253442e-05 | 2.213301e-05 | 2.294575e-05 | 2.327524e-05 | 2.179358e-05
H3 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 [decayed]
He3 | 1.352834e-05 | 1.358545e-05 | 1.347749e-05 | 1.461604e-05 | 1.244062e-05
He4 | 6.100487e-02 | 6.098777e-02 | 6.102189e-02 | 6.100321e-02 | 6.100654e-02
Li6 | 3.077093e-14 | 5.069457e-14 | 2.249760e-14 | 3.077093e-14 | 3.077093e-14
Li7 | 3.932416e-10 | 4.187589e-10 | 3.698327e-10 | 3.932416e-10 | 3.932416e-10
Be7 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 [decayed]
INFO : The final abundances have been written to 'io/test/bbundance_file.dat'.
Xena expects an abundance-file with name abundance_file.dat in <io_directory> with the same structure as ACROPOLIS
Xena expects a param-file with name param_file.dat in <io_directory> with at least one line that reads eta=<eta>, e.g. eta=6.137e-10
When using either CascadeNeutrinoModel or EmDecayModel, Xena further expects entries for mphi [in MeV] and tau [in s], as well as for Bree and Braa in the latter case.
Xena expects a cosmo-file with name cosmo_file.dat in <io_directory> with the column structure
tin sTin MeVdTdtin MeV²Tnuin MeVHin MeVnb_etafin MeV³
When using the model NeutrinoCascadeModel, Xena aditionally expects six more column with
Semin MeV⁴Shdin MeV⁴dndt_scin MeV⁴Khd_scin MeVdndt_shin MeV⁴Khd_shin MeV
When using the model EmDecayModel, Xenaadditionally expects a seventh column with
nphiin MeV³