Code and results for the disruption of axion miniclusters (AMCs) in the Milky Way, as well as radio signals from encounters between AMCs and neutrons stars.
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.
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.
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
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.