A new algorithm for the solution and analysis of the Collisional Radiative Model equations

A new general algorithm for the solution of the Collisional Radiative Model (CRM) for large atomic models is presented. It is based on the separation of total population, and reduced populations for the individual states of each charge state. The sum of the latter is unity. This leads to a double linearized coupled iterative system, including cycling on all charge states. This algorithm is robust. It gives insight into the quality of the atomic model by allowing the follow up of the charge states populations, the global rates and the average charge Z* as a function of the iteration index. We show results on Xe with 39,683 configuration states.