Complex pseudo-partition functions in the Configurationally-Resolved Super-Transition-Array approach for radiative opacity

A few years ago, Kurzweil and Hazak developed the Configurationally Resolved Super-Transition-Arrays (CRSTA) method for the computation of hot-plasma radiative opacity. Their approach, based on a temporal integration, is an important refinement of the standard Super-Transition-Arrays (STA) approach, which enables one to recover the underlying structure of the STAs, made of unresolved transition arrays. The CRSTA formalism relies on the use of complex pseudo partition functions, depending on the considered one-electron jump. In this article, we find that, despite the imaginary part, the doubly-recursive relation which was introduced in the original STA method to avoid problems due to alternating-sign terms in partition functions, is still applicable, robust, efficient, and exempt of numerical instabilities. This was rather unexpected, in particular because of the occurrence of trigonometric functions, or Chebyshev polynomials, which can be either positive or negative. We also show that, in the complex case, the recursion relation can be presented in a form where the vector of real and imaginary parts at a given iteration is therefore obtained by a sum of the rotated previous ones.