Electron density of elements

[qu...@mail.ru]

Hi, everyone!

I'm confused and I really want to figure it out

Have PyCloudy a total electron density for each chemical element in output  commands? I use ne, but apparently this is the total electron density for all elements.
If you set the abundance and so on. Is it possible to get exactly the electron density of the selected chemical element(for example, all O ions) or the selected ion (for example, OIII) at the output?

And if it is possible, what values are used to calculate such an electron density inside CLOUDY?

Thank you!!
Sincerely Yours,

Masha

[lzc...@gmail.com]

Hi Masha,

As far as I understand it, PyCloudy by default includes several 'save last element' commands  in the input file (as below) to save the ionic structure vs optical depth for specified ions. It will generate several files after finishing calculation, with the corresponding suffixes .ele_H, .ele_He, etc. You can later on access those using the 'get_ionic' command (https://morisset.github.io/pyCloudy_Manual/html/classpy_cloudy_1_1c1d_1_1cloudy__model_1_1_cloudy_model.html#a17019e1a3c0f0580a132e5de9822a020) and calculate the densities of those ions using the ionic fraction, hydrogen density and abundance. A relevant example is in the tutorial here: https://github.com/Morisset/pyCloudy/blob/master/pyCloudy/docs/Using_pyCloudy_1.ipynb

Best,

Zeren

save last element hydrogen ".ele_H"

save last element helium ".ele_He"

save last element carbon ".ele_C"

save last element nitrogen ".ele_N"

save last element oxygen ".ele_O"

save last element argon ".ele_Ar"

save last element neon ".ele_Ne"

save last element sulphur ".ele_S"

save last element chlorin ".ele_Cl"

save last element iron ".ele_Fe"

save last element silicon ".ele_Si"