fiasco icon indicating copy to clipboard operation
fiasco copied to clipboard

Add calculation of dielectronic level populations

Open jwreep opened this issue 1 year ago • 2 comments

I don't see another issue for this, so making a new one. See the discussion here.

Fiasco does not currently have a method to calculate level populations for a dielectronic ion.

jwreep avatar Feb 22 '24 21:02 jwreep

I think it's worth clarifying what is meant in this issue since I didn't describe it in particular detail. There are no "dielectronic ions" in reality, but CHIANTI creates a separate model ion for dielectronic excitation to calculate satellite lines. The database includes directories like "fe_17" and "fe_17d" to distinguish the two.

This is described in Section 2.2 of the fourth CHIANTI paper:

Screenshot 2024-12-11 at 12 35 24 PM Screenshot 2024-12-11 at 12 35 29 PM

jwreep avatar Dec 11 '24 22:12 jwreep

I agree that this needs to be worked on, but I'm not quite sure how to go about it. Maybe this is because I don't fully understand how it is done in the IDL code.

In v9 and above, dielectronic recombination is included in the level populations calculation through a two-ion model approach. See #254 for more details. However, prior to v9, these so-called satellite lines were treated separately by including a separate ion with a "_d" suffix in the database. This is still true for some ions even in version 9 and above.

What I'm not sure about is how these "_d" data are actually used. In fiasco, if there is a "dielectronic version" of the data, it can be accessed as _dielectronic_<dataset_name>, e.g. fe_15._dielectronic_elvlc would give you the data from the file fe/fe_15d/fe_15d.elvlc. Is it just a matter of making additional Ion objects for these "dielectronic versions" of each ion if they exist?

wtbarnes avatar Mar 27 '25 20:03 wtbarnes