A New Implementation of the STA Method for the Calculation of Opacities of Local Thermodynamic Equilibrium Plasmas

Abstract: We present opacity calculations with the newly developed STAR code, which implements the Super-Transition-Array (STA), with various improvements. The model is used to calculate and analyze local thermodynamic equilibrium opacities of mid and high Z elements and of the solar interior plasma. We briefly review the underlying computational model and present calculations for iron and neodymium over a wide range of temperature and density.


Introduction
The calculation of atomic opacities from first principles is an important part in the modeling of various astrophysical phenomena, and especially the physics of stellar interior.Opacities quantify the strong coupling between radiation and matter in a hot dense plasma, and therefore, are directly related to the radiative thermal conductivity in stellar interior.
In recent years, a solar composition problem has emerged [1,2], as a result of a revision in the solar photospheric abundances, in which the metallic abundance was significantly decreased [3].Metallic elements have a significant contribution to the opacity in the solar interior, although they only constitute a few percent of the mixture, since these metallic elements are not completely ionized and give rise to strong bound-bound and bound-free absorption.This gives rise to a connection between the solar composition problem and the theoretical uncertainty in the calculation of opacities at stellar interior conditions [4][5][6][7].
We have recently developed [6][7][8][9][10] the atomic code STAR (STA-Revised), for the calculation of opacities of local thermodynamic equilibrium plasmas, by the STA method [11].The code was used to investigate and analyze sensitivities and uncertainties in the calculation of solar interior opacities.
In this paper we describe the underlying computational model and present the computational capabilities of STAR.We calculate the opacity spectra of iron, neodymium and the solar mixture, as well as Rosseland and Planck mean free paths for a wide range of temperature and density.We also present maps for the contributions of the different atomic processes to the Rosseland opacity.

The Model
The monochromatic opacity is given as a sum of four different processes: (i) photon scattering and (ii) free-free (iii) bound-free and (iv) bound-bound absorption.The cross-section for photon-absorption (including stimulated emission) for the ith component of a mixture plasma at temperature T and density ρ, is written as: where k B is Boltzmann's constant.The monochromatic opacity is by definition the absorption coefficient per unit mass: where A i is the atomic mass of the i component and N A is Avogadro's constant.The total monochromatic opacity of the mixture is given by the sum of the individual opacities, weighted by the mass fraction: The Rosseland mean opacity is given by: with the Rosseland weight function: where u = E/k B T.
A simple way to measure fractional contributions, such as the contribution of different elements in a mixture or different atomic processes, to the total Rosseland mean, is the following: given N individual spectral contributions: we define the fractional Rosseland contributions by: where κ i R is the Rosseland mean (4) of the cumulative spectra ∑ i j=1 κ j (E), while for the first contributer δκ 1 = κ 1 R /κ N R .The photon scattering cross section is written as (see Refs. [10,12] and references therein): where Z f is the average number of free electrons, η = µ/k B T is the reduced chemical potential and the Thomson cross-section is with α the fine structure constant and a 0 the Bohr radius.G(u, T ), where T = k B T/m e c 2 , corrects the Klein-Nishina cross-section due to the finite temperature and contains corrections for inelastic scattering.The factor R(η, T ) contains corrections due to the Pauli blocking as a result of partial degeneracy and includes a relativistic electron dispersion relation.The factor f (η, δ) is a correction due to plasma collective effects.
The free-free photoabsorption cross section is calculated via a screened-hydrogenic approximation with a multiplicative degeneracy correction, and is given by: where the Kramers cross-section is: with n e the free electron number density.For the (thermal) averaged free-free Gaunt factor ḡ f f , we use a screened-hydrogenic approximation and a Maxwellian distribution via the tables given in Ref. [13].
In principle, the photoexcitation opacity can result from all possible allowed transitions between pairs of bound-states, from all populated electronic configurations in the plasma.In practice, for mid and high-Z elements, typical configurations may contain a huge number of lines [14,15], so that Detailed-Line-Accounting codes [16][17][18][19][20], and the Unresolved-Transition-Array (UTA) method is used to group this huge number of lines [9,[21][22][23][24][25].However, in many situations, for hot-dense plasmas and for Z 30, there is an intractable number of populated configurations, and a detailed accounting of all UTAs in the spectra is computationally impossible.In Refs.[10,26] the number of populated relativistic configurations was estimated by the formula: where: is the variance of the population of shell s, n s = 1/(e −( s −µ)/k B T + 1) is the Fermi-Dirac distribution, g s = 2j s + 1 is the orbital degeneracy, s is the orbital energy and µ is the chemical potential.It is seen from Equation ( 12) that the occupation numbers of shells whose energies are near the Fermi-Dirac step | s − µ| ≈ k B T fluctuate, while the other shells are either filled or empty.Thus, when the temperature is high enough, such that there are multiple shells that lie within the Fermi-Dirac step, the number of populated configurations can become huge.A calculation of N C for iron (Z = 26) and gold (Z = 79) in a wide range of temperature and density is given in Figure 1.The calculations were performed via the ion-sphere model implemented in STAR.It is evident that for iron the number of populated configurations can become larger than ∼10 15 , while for gold, this number can become larger than ∼10 25 .It is also evident that N C has a single maximum in temperature for each density, and in density for each temperature.This is explained by the fact that shells become ionized for high temperatures or full for low temperatures, and pressure ionized for high densities.For low densities, shell populations become non-degenerate-which reduces the fluctuations in shells occupation numbers.It is also seen, as expected, that the maximum of N C for gold has a larger temperature and density than iron has, due to the higher nuclear attraction and number of bound shells.STAR implements the STA method [11] for the calculation of bound-bound and bound-free photoabsorption spectra.The STA method enables the grouping of the huge number of lines between the enormous number of configurations into large sets, called Super-Transition-Arrays (STAs).The STA moments are calculated analytically via the partition function algebra, and are split into smaller STAs until convergence is achieved.This procedure enables the handling of situations where the number of populated configurations is too large, as shown in Figure 1.STAR also implements several additional advanced capabilities, such as (i) an adaptive integration of resonances in the electronic density of state [27][28][29][30], which has an effect on the bound-free structure and on the self-consistent field average-atom calculation, (ii) ion-sphere and ion-correlation models for the plasma environment [7,27,[31][32][33][34] and (iii) stable recursive calculation of partition functions, as in Refs.[8,[35][36][37].12)) for iron (upper pane) and gold (lower pane).

Opacity Calculations
In Figures 2 and 3, spectral opacities for iron (Z = 26) and neodymium (Z = 60) are presented in detail.It is seen in Figure 3 how for neodymium, lines are "condensed" into transition arrays due to the huge number of lines.Figure 4 shows in detail the opacity spectra of the solar mixture at R = 0.715R (near the convection zone).It is evident that the opacity contribution of the fully ionized hydrogen and helium, which consist more than 99% of the mixture, is only less than 10%.On the other hand, it is seen that most of the opacity is due to bound-bound and bound-free photoabsorption of various mid-Z elements such as oxygen and neon (via K-shell absorption), as well as iron (via L-shell absorption).
In Figures 5 and 6 we present maps of the Rosseland and Planck mean free paths for iron and neodymium, respectively, over a wide range of temperature (100 eV-10 keV) and density (0.001-1000 g/cm 3 ).In Figures 7 and 8 maps for the contributions of the atomic fractions (scattering, free-free, bound-free and bound-bound) to the Rosseland opacity are shown.It is evident that in both cases the bound-bound and bound-free processes are major contributers over most of the temperature-density range, and that the bound-bound contribution is much more dominant for neodymium, due to its higher atomic number.The total (blue), as well as the bound-bound (black), bound-free (red), free-free (magenta) and scattering (lime) opacity spectra are shown (left axis), together with the Rosseland wieght function (green, right axis).The Rosseland mean opacity and the processes contributions to it, are given in the legend.We note that extensive comparisons between STAR and other widely used opacity codes (specifically, OP [16], OPAL [17] and OPAS [38]), were performed in Ref. [10] for the case solar opacities.These comparisons resulted in a very satisfactory agreement in the Rosseland opacity of the solar mixture of about 5% throughout the solar interior.Spectral opacities and ionic populations of several metallic elements were also compared at the thermodynamic conditions of the solar convection zone, and good agreement was reached.

Figure 1 .
Figure 1.Maps on the temperature-density plane of the number of occupied relativistic configurations (see Equation (12)) for iron (upper pane) and gold (lower pane).

3 Figure 2 .
Figure 2. Opacity spectrum, calculated by STAR, for iron (Z = 26) at T = 200 eV and ρ = 0.5 g/cm3 .The total (blue), as well as the bound-bound (black), bound-free (red), free-free (magenta) and scattering (lime) opacity spectra are shown (left axis), together with the Rosseland wieght function (green, right axis).The Rosseland mean opacity and the processes contributions to it, are given in the legend.

Figure 4 .Figure 5 .
Figure 4.The opacity spectrum for the solar mixture at the conditions found at R = 0.715R (near the solar convection-zone) where T = 183.8eV and ρ = 0.1745 g/cm3 .The total opacity due to all element in the solar mixture (solid blue) is compared to the opacity due to hydrogen and helium only (solid red), and the spectra of various individual elements, with a Rosseland opacity contribution that is larger than 1% (dashed lines).The Rosseland mean free path l R and the Rosseland opacity contributions δκ R , are listed in the legend.The Rosseland weight function (solid black) is also shown.

Figure 7 .
Figure 7. Maps on the temperature-density plane of the scattering (upper left pane), free-free (upper right pane), bound-free (lower left pane) and bound-bound (lower right pane) atomic processes contributions to the Rosseland opacity for iron.