bradkav / axion-miniclusters Goto Github PK

Add a reference for the value of rho_eq in AMC.py

It would be good to have a piece of code which checks which NSs in the interaction list have a B-field strong enough to allow for conversion (at a given axion mass). Then this could be used to correct the encounter rate to include only converting NSs.

It would be good to add a "verbose" option to all of the functions in the pipeline, so that we can suppress some/all output when we loop over different mass functions.

In the genericAMC branch, the monte carlo fails for NFW profiles. It seems that when the Monte Carlo 'adjusts' the maximum impact parameter, it isn't working correctly, because it ends up sampling more interactions than the maximum N_cut.

We need to implement the cross section (in prepare_distributions.py) according to Sec. II of the Supplementary materials in https://arxiv.org/abs/2109.04582.

The tabulated mass loss for NFW miniclusters needs to be re-calculated, because there was a typo in the calculation of the potential. The correct potential is attached.
Eddington_Distribution_Functions.pdf

The mass loss can then be compared with e.g. Fig. 3 in https://arxiv.org/pdf/2207.11276.pdf.

Two different sets of 'inverse_transform_sampling' functions exist - in both 'perturbations.py' and 'NSencounter.py'

In the genericAMC branch, we need to double-check that the AScut is behaving properly for extended mass functions. Including the AScut reducing the rate by ~5-6 orders of magnitude.