LTE Plasma

The LTEPlasma plasma class is the child of BasePlasma but is the first class that actually calculates plasma conditions. After running exactly through the same steps as BasePlasma, LTEPlasma will start calculating the partition functions.

Zi,j=max(k)k=0gk×eEk/(kbT)

where Z is the partition function, g is the degeneracy factor, E the energy of the level and T the temperature of the radiation field.

The next step is to calculate the ionization balance using the Saha ionization equation. and then calculating the Number density of the ions (and an electron number density) in a second step. First ge=(2πmekBTradh2)3/2 is calculated (in LTEPlasma.update_t_rad), followed by calculating the ion fractions (LTEPlasma.calculate_saha).

Ni,j+1×NeNi,j=Φi,(j+1)/jΦi,(j+1)/j=ge×Zi,j+1Zi,jeχjj+1/kBT

In the second step (LTEPlasma.calculate_ionization_balance), we calculate in an iterative process the electron density and the number density for each ion species.

N(X)=N1+N2+N3+N(X)=N1+N2N1N1+N3N2N2N1N1+N4N3N3N2N2N1N1+N(X)=N1(1+N2N1+N3N2N2N1+N4N3N3N2N2N1+)N(X)=N1(1+Φi,2/1Ne+Φi,3/2NeΦi,2/1Ne+Φi,4/3NeΦi,3/2NeΦi,2/1Ne+)αN1=N(X)α

Initially, we set the electron density (Ne) to the sum of all atom number densities. After having calculated the ion species number densities, we recalculate the electron density by weighting the ion species number densities with their ion number (e.g. neutral ion number densities don’t contribute at all to the electron number density, once ionized contribute with a factor of 1, twice ionized contribute with a factor of two, ….).

Finally, we calculate the level populations (LTEPlasma.calculate_level_populations) by using the calculated ion species number densities:

Ni,j,k=gkZi,j×Ni,j×eβradEk

This concludes the calculation of the plasma. In the code, the next step is calculating the τSobolev using the quantities calculated here.

Example Calculations