johnson-c / ColRadPy

Collisional-radiative python solver
Other
27 stars 10 forks source link

Ionization balance with charge exchange using ADF11 rates #27

Open dmkriete opened 2 months ago

dmkriete commented 2 months ago

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: Figure_1

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: Figure_2

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).