ColRadPy
ColRadPy copied to clipboard
Ionization balance with charge exchange using ADF11 rates
Currently, it is not possible to run an ionization balance with charge exchange included when ADF11 rates are being used; you have to use ADF04 rates. The reason for this appeared to be that reading CX rates from ADF11 ccd files was not implemented in read_adf11
. I implemented this and now it is possible to run the ionization balance calculation.
I am using carbon as a test case. Here is what the steady-state ionization balance looks like for n_e = 1e19 m^-3 with no CX:
And here is what it looks like with CX included for a neutral hydrogen density and temperature of 5e17 m^-3 and 5 eV, respectively:
Ideally, these results would be compared against an ionization balance where the CX rates are determined from ADF04 files. However, I don't know how to run ionization balances with ADF04 rates. That being said they look qualitatively similar to the results from previous ionization balances with CX included (e.g. Adam McLean's 2021 PSI contribution).