The ion populations most frequently adopted for diagnostics in collisional plasmas are derived from the density independent, coronal approximation. In higher density, lower temperature conditions, ionisation rates are enhanced once metastable levels become populated, and recombination rates are suppressed if ions recombine into Rydberg levels. As a result, the formation temperatures of ions shift, altering the diagnostics of the plasma. To accurately model the effect of ionisation from metastable levels, new electron impact, ionisation cross sections have been calculated for oxygen, both for direct ionisation and excitation--auto-ionisation of the ground and metastable levels. The results have been incorporated into collisional radiative modelling to show how the ionisation equilibrium of oxygen changes once metastable levels become populated. Suppression of dielectronic recombination has been estimated and also included in the modelling, demonstrating the shifts with density in comparison to the coronal approximation. The final results for the ionisation equilibrium are used in differential emission measure modelling to predict line intensities for many lines emitted by O II-VI in the solar transition region. The predictions show improved agreement by 15-40% for O II, O VI and the inter-combination lines of O III-V, when compared to results from coronal approximation modelling. While there are still discrepancies with observations of these lines, this could, to a large part, be explained by variability in the observations.

Recent improvements to atomic energy-level data allow, for the first time, accurate predictions to be made for the Fe III line emission strengths in the spectra of luminous, $L_\text{bol}=10^{46}-10^{48}$ erg/s, Active Galactic Nuclei. The Fe III emitting gas must be primarily photoionized, consistent with observations of line reverberation. We use CLOUDY models exploring a wide range of parameter space, together with 26,500 rest-frame ultraviolet spectra from the Sloan Digital Sky Survey, to constrain the physical conditions of the line emitting gas. The observed Fe III emission is best accounted for by dense ($n_H=10^{14}$ cm$^{-3}$) gas which is microturbulent, leading to smaller line optical depths and fluorescent excitation. Such high density gas appears to be present in the central regions of the majority of luminous quasars. Using our favoured model, we present theoretical predictions for the relative strengths of the Fe III UV34 $\lambda\lambda$1895,1914,1926 multiplet. This multiplet is blended with the Si III] $\lambda$1892 and C III] $\lambda$1909 emission lines and an accurate subtraction of UV34 is essential when using these lines to infer information about the physics of the broad line region in quasars.

We present new collisional-radiative models (CRMs) for helium in the quiescent solar corona, and predict the emissivities of the He and He$^+$ lines to be observed by DKIST, Solar Orbiter, and Proba-3. We discuss in detail the rates we selected for these models, highlighting several shortcomings we have found in previous work. As no previous complete and self-consistent coronal CRM for helium existed, we have benchmarked our largest model at a density of 10$^{6}$ cm$^{-3}$ and temperature of 20,000 K against recent CRMs developed for photoionised nebulae. We then present results for the outer solar corona, using new dielectronic recombination rates we have calculated, which increase the abundance of neutral helium by about a factor of two. We also find that all the optical triplet He I lines, and in particular the well known He I 10830 and 5876 A lines are strongly affected by both photo-excitation and photo-ionisation from the disk radiation, and that extensive CRM models are required to obtain correct estimates. Close to the Sun, at an electron density of 10$^{8}$ cm$^{-3}$ and temperature of 1 MK, we predict the emissivity of the He I 10830 A to be comparable to that of the strong Fe XIII coronal line at 10798 A. However, we expect the He I emissivity to sharply fall in the outer corona, with respect to Fe XIII. We confirm that the He$^+$ Lyman $\alpha$ at 304 A is also significantly affected by photo-excitation and is expected to be detectable as a strong coronal line up to several solar radii.

Emission and absorption features from C-like ions serve as temperature and density diagnostics of astrophysical plasmas. $R$-matrix electron-impact excitation data sets for C-like ions in the literature merely cover a few ions, and often only for the ground configuration. Our goal is to obtain level-resolved effective collision strength over a wide temperature range for C-like ions from \ion{N}{II} to \ion{Kr}{XXXI} (i.e., N$^{+}$ to Kr$^{30+}$) with a systematic set of $R$-matrix calculations. We also aim to assess their accuracy. For each ion, we included a total of 590 fine-structure levels in both the configuration interaction target and close-coupling collision expansion. These levels arise from 24 configurations $2l^3 nl^{\prime}$ with $n=2-4$, $l=0-1$, and $l^{\prime}=0-3$ plus the three configurations $2s^22p5l$ with $l=0-2$. The AUTOSTRUCTURE code was used to calculate the target structure. Additionally, the $R$-matrix intermediate coupling frame transformation method was used to calculate the collision strengths. We compare the present results of selected ions with archival databases and results in the literature. The comparison covers energy levels, transition rates, and effective collision strengths. We illustrate the impact of using the present results on an \ion{Ar}{xiii} density diagnostic for the solar corona. The electron-impact excitation data is archived according to the Atomic Data and Analysis Structure (ADAS) data class adf04 and will be available in OPEN-ADAS. The data will be incorporated into spectral codes, such as CHIANTI and SPEX, for plasma diagnostics.

At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact parameter calculations or semi-empirical fits and the rate coefficients agree within a factor of two. We show in this paper that these uncertainties cause errors lower than 5% in the emission of radio recombination lines (RRL) of most ionized plasmas of typical nebulae. However, in special circumstances where the transitions between Rydberg levels are amplified by maser effects, the errors can increase up to 20%. We present simulations of the optical depth and H$n\alpha$ line emission of Active Galactic Nuclei (AGN) Broad Line Regions (BLRs) and the Orion Nebula Blister to showcase our findings.

We consider three recent large-scale calculations for the radiative and electron-impact excitation data of N IV, carried out with different methods and codes. The scattering calculations employed the relativistic Dirac $R$-matrix (DARC) method, the intermediate coupling frame transformation (ICFT) $R$-matrix method, and the B-spline $R$-matrix (BSR) method. These are all large-scale scattering calculations with well-tested and sophisticated codes, which use the same set of target states. One concern raised in previous literature is related to the increasingly large discrepancies in the effective collision strengths between the three sets of calculations for increasingly weak and/or high-lying transitions. We have built three model ions and calculated the intensities of all the main spectral lines in this ion. We have found that, despite such large differences, excellent agreement (to within $\pm$~20\%) exists between all the spectroscopically-relevant line intensities. This provides confidence in the reliability of the calculations for plasma diagnostics. We have used the differences in the radiative and excitation rates amongst the three sets of calculations to obtain a measure of the uncertainty in each rate. Using a Monte Carlo approach, we have shown how these uncertainties affect the main theoretical ratios which are used to measure electron densities and temperatures.

The present debate on the reliability of astrophysical opacities has reached a new climax with the recent measurements of Fe opacities on the Z-machine at the Sandia National Laboratory \citep{Bailey2015}. To understand the differences between theoretical results, on the one hand, and experiments on the other, as well as the differences among the various theoretical results, detailed comparisons are needed. Many ingredients are involved in the calculation of opacities; deconstructing the whole process and comparing the differences at each step are necessary to quantify their importance and impact on the final results. We present here such a comparison using the two main approaches to calculate the required atomic data, the $R$-Matrix and distorted-wave methods, as well as sets of configurations and coupling schemes to quantify the effects on the opacities for the $Fe\ XVII$ and $Ni\ XIV$ ions.

We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully used. The use of such large datasets is not realistic for most calculations due to the time and memory needs, and we describe the limited subset of data we use by default. Despite the fact that we now predict significantly more lines than the previous Cloudy release, this version is faster because of optimization of memory access patterns and other tuning. The size and use of the databases can easily be adjusted in the command-line interface. We give examples of the accuracy limits using small models, and the performance requirements of large complete models. We summarize several advances in the H- and He-like iso-electronic sequences. We use our complete collisional-radiative models of the ionization of these one and two-electron ions to establish the highest density for which the coronal or interstellar medium (ISM) approximation works, and the lowest density where Saha or local thermodynamic equilibrium can be assumed. The coronal approximation fails at surprisingly low densities for collisional ionization equilibrium but is valid to higher densities for photoionized gas clouds. Many other improvements to the physics have been made and are described. These include the treatment of isotropic continuum sources such as the cosmic microwave background (CMB) in the reported output, and the ability to follow the evolution of cooling non-equilibrium clouds.

We compare the results of the semi-classical (SC) and quantum-mechanical (QM) formalisms for angular-momentum changing transitions in Rydberg atom collisions given by Vrinceanu & Flannery, J. Phys. B 34, L1 (2001), and Vrinceanu, Onofrio & Sadeghpour, ApJ 747, 56 (2012), with those of the SC formalism using a modified Monte Carlo realization. We find that this revised SC formalism agrees well with the QM results. This provides further evidence that the rates derived from the QM treatment are appropriate to be used when modelling recombination through Rydberg cascades, an important process in understanding the state of material in the early universe. The rates for $\Delta\ell=\pm1$ derived from the QM formalism diverge when integrated to sufficiently large impact parameter, $b$. Further to the empirical limits to the $b$ integration suggested by Pengelly & Seaton, MNRAS 127, 165 (1964), we suggest that the fundamental issue causing this divergence in the theory is that it does not fully cater for the finite time taken for such distant collisions to complete.

Accurate rates for energy-degenerate l-changing collisions are needed to determine cosmological abundances and recombination. There are now several competing theories for the treatment of this process, and it is not possible to test these experimentally. We show that the H I two-photon continuum produced by astrophysical nebulae is strongly affected by l-changing collisions. We perform an analysis of the different underlying atomic processes and simulate the recombination and two-photon spectrum of a nebula containing H and He. We provide an extended set of effective recombination coefficients and updated l-changing 2s-2p transition rates using several competing theories. In principle, accurate astronomical observations could determine which theory is correct.

For photoionized plasmas, electron energy-loss rates due to radiative recombination (RR) are required for thermal equilib- rium calculations, which assume a local balance between the energy gain and loss. While many calculations of total and/or partial RR rates are available from literature, specific calculations of associated RR electron energy-loss rates are lacking. Here we focus on electron energy-loss rates due to radiative recombination of H-like to Ne-like ions for all the elements up to and including zinc (Z = 30), over a wide temperature range. We use the AUTOSTRUCTURE code to calculate the level-resolved photoionization cross section and modify the ADASRR code so that we can simultaneously obtain level-resolved RR rate coefficients and associated RR electron energy-loss rate coefficients. The total RR rates and electron energy-loss rates of H i and He i are compared with those found in literature. Furthermore, we utilize and parameterize the weighted electron energy-loss factors (dimensionless) to characterize total electron energy-loss rates due to RR. The RR electron energy-loss data are archived according to the Atomic Data and Analysis Structure (ADAS) data class adf48. The RR electron energy-loss data are also incorporated into the SPEX code for detailed modelling of photoionized plamsas.

To calculate realistic models of objects with Ni in their atmospheres, accurate atomic data for the relevant ionization stages needs to be included in model atmosphere calculations. In the context of white dwarf stars, we investigate the effect of changing the Ni {\sc iv}-{\sc vi} bound-bound and bound-free atomic data has on model atmosphere calculations. Models including PICS calculated with {\sc autostructure} show significant flux attenuation of up to $\sim 80$\% shortward of 180\AA\, in the EUV region compared to a model using hydrogenic PICS. Comparatively, models including a larger set of Ni transitions left the EUV, UV, and optical continua unaffected. We use models calculated with permutations of this atomic data to test for potential changes to measured metal abundances of the hot DA white dwarf G191-B2B. Models including {\sc autostructure} PICS were found to change the abundances of N and O by as much as $\sim 22$\% compared to models using hydrogenic PICS, but heavier species were relatively unaffected. Models including {\sc autostructure} PICS caused the abundances of N/O {\sc iv} and {\sc v} to diverge. This is because the increased opacity in the {\sc autostructure} PICS model causes these charge states to form higher in the atmosphere, moreso for N/O {\sc v}. Models using an extended line list caused significant changes to the Ni {\sc iv}-{\sc v} abundances. While both PICS and an extended line list cause changes in both synthetic spectra and measured abundances, the biggest changes are caused by using {\sc autostructure} PICS for Ni.

Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which focused on H I, we showed that differences in $l$-changing collisional rate coefficients predicted by three different theories can translate into 10% changes in predictions for H I spectra. Here we consider the more complicated case of He atoms, where low-$l$ subshells are not energy degenerate. A criterion for deciding when the energy separation between $l$ subshells is small enough to apply energy-degenerate collisional theories is given. Moreover, for certain conditions, the Bethe approximation originally proposed by Pengelly & Seaton (1964) is not sufficiently accurate. We introduce a simple modification of this theory which leads to rate coefficients which agree well with those obtained from pure quantal calculations using the approach of Vrinceanu et al. (2012). We show that the $l$-changing rate coefficients from the different theoretical approaches lead to differences of $\sim 10$% in He I emissivities in simulations of H II regions using spectral code Cloudy.

Emission lines from highly-excited states (n >= 5) of H- and He-like ions have been detected in astrophysical sources and fusion plasmas. For such excited states, R-matrix or distorted wave calculations for electron-impact excitation are very limited, due to the large size of the atomic basis set needed to describe them. Calculations for n >= 6 are also not generally available. We study the behaviour of the electron-impact excitation collision strengths and effective collision strengths for the most important transitions used to model electron collision dominated astrophysical plasmas, solar, for example. We investigate the dependence on the relevant parameters: the principal quantum number n or the nuclear charge Z. We also estimate the importance of coupling to highly-excited states and the continuum by comparing the results of different sized calculations. We provide analytic formulae to calculate the electron-impact excitation collision strengths and effective collision strengths to highly-excited states (n >= 8) of H- and He-like ions. These extrapolated effective collision strengths can be used to interpret astrophysical and fusion plasma via collisional-radiative modelling.

Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe dipole $l$-changing collisions ($nl\rightarrow nl'=l\pm 1$) between energy-degenerate states within an $n$-shell. The work of Pengelly & Seaton (1964) has, for half-a-century, been considered the definitive study which "solved" the problem. Recent work by Vrinceanu et al.(2012) recommended the use of rate coefficients from a semi-classical approximation which are nearly an order of magnitude smaller than those of Pengelly & Seaton (1964), with the result that significantly higher densities are needed for the $nl$ populations to come into local thermodynamic equilibrium. Here, we compare predicted H~I emissivities from the two works and find widespread differences, of up to $\approx 10$%. This far exceeds the 1% precision required to obtain the primordial He/H abundance ratio from observations so as to constrain Big Bang cosmologies. We recommend using the rate coefficients of Pengelly & Seaton (1964) for $l$-changing collisions, to describe the H recombination spectrum, based-on their quantum mechanical representation of the long-range dipole interaction.

We have carried-out two intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of C-like Fe $^{20+}$, both of which use the same expansions for their configuration interaction (CI) and close-coupling (CC) representations. The first expansion arises from the configurations 2s$^2$ 2p$^2$, 2s 2p$^3$, 2p$^4$, {2s$^2$ 2p, 2s 2p$^2$, 2p$^3$} nl, with n=3,4, for l=0-3, which give rise to 564 CI/CC levels. The second adds configurations 2s$^2$ 2p 5l, for l=0-2, which give rise to 590 CI/CC levels in total. Comparison of oscillator strengths and effective collision strengths from these two calculations demonstrates the lack of convergence in data for n=4 from the smaller one. Comparison of results for the 564 CI/CC level calculation with an earlier ICFT R-matrix calculation which used the exact same CI expansion but truncated the CC expansion to only 200 levels demonstrates the lack of convergence of the earlier data, particularly for n=3 levels. Also, we find that the results of our 590 CC R-matrix calculation are significantly and systematically larger than those of an earlier comparable Distorted Wave-plus-resonances calculation. Thus, it is important still to take note of the (lack of) convergence in both atomic structural and collisional data, even in such a highly-charged ion as Fe $^{20+}$, and to treat resonances non-perturbatively. This is of particular importance for Fe ions given their importance in the spectroscopic diagnostic modelling of astrophysical plasmas.

We have carried-out 98-level configuration-interaction / close-coupling (CI/CC) intermediate coupling frame transformation (ICFT) and Breit-Pauli R-matrix calculations for the electron-impact excitation of Be-like Al$^{9+}$. The close agreement that we find between the two sets of effective collision strengths demonstrates the continued robustness of the ICFT method. On the other hand, a comparison of this data with previous 238-level CI/CC ICFT effective collision strengths shows that the results for excitation up to n=4 levels are systematically and increasingly underestimated over a wide range of temperatures by R-matrix calculations whose close-coupling expansion extends only to n=4 (98-levels). Thus, we find to be false a recent conjecture that the ICFT approach may not be completely robust. The conjecture was based upon a comparison of 98-level CI/CC Dirac R-matrix effective collision strengths for Al$^{9+}$ with those from the 238-level CI/CC ICFT R-matrix calculations. The disagreement found recently is due to a lack of convergence of the close-coupling expansion in the 98-level CI/CC Dirac work. The earlier 238-level CI/CC ICFT work has a superior target to the 98-level CI/CC Dirac one and provides more accurate atomic data. Similar considerations need to be made for other Be-like ions and for other sequences.

Dielectronic recombination (DR) is the dominant recombination process for most heavy elements in photoionized clouds. Accurate DR rates for a species can be predicted when the positions of autoionizing states are known. Unfortunately such data are not available for most third and higher-row elements. This introduces an uncertainty that is especially acute for photoionized clouds, where the low temperatures mean that DR occurs energetically through very low-lying autoionizing states. This paper discusses S$^{2+} \rightarrow$ S$^+$ DR, the process that is largely responsible for establishing the [S~III]/[S~II] ratio in nebulae. We derive an empirical rate coefficient using a novel method for second-row ions, which do have accurate data. Photoionization models are used to reproduce the [O~III] / [O~II] / [O~I] / [Ne~III] intensity ratios in central regions of the Orion Nebula. O and Ne have accurate atomic data and can be used to derive an empirical S$^{2+} \rightarrow$ S$^+$ DR rate coefficient at $\sim 10^{4}$~K. We present new calculations of the DR rate coefficient for S$^{2+} \rightarrow$ S$^+$ and quantify how uncertainties in the autoionizing level positions affect it. The empirical and theoretical results are combined and we derive a simple fit to the resulting rate coefficient at all temperatures for incorporation into spectral synthesis codes. This method can be used to derive empirical DR rates for other ions, provided that good observations of several stages of ionization of O and Ne are available.

Arguably, the best method for determining the effective temperature ($T_{\mathrm{eff}}$) and surface gravity (log $g$) of a DA white dwarf is by fitting the Hydrogen Lyman and Balmer absorption features. However, as has been shown for white dwarfs with $T_{\mathrm{eff}}$>50,000K, the calculated value from the Lyman and Balmer lines are discrepant, which worsens with increasing temperature. Many different solutions have been suggested, ranging from the input physics used to calculate the models, to interstellar reddening. We will focus on the former, and consider three variables. The first is the atomic data used, namely the number of transitions included in line blanketing treatments and the photoionization cross sections. The second is the stark broadening treatment used to synthesise the Lyman and Balmer line profiles, namely the calculations performed by Lemke (1997) and Tremblay & Bergeron (2009). Finally, the third is the atmospheric content. The model grids are calculated with a pure H composition, and a metal polluted composition using the abundances of Preval et al. (2013). We present the preliminary results of our analysis, whereby we have determined the $T_{\mathrm{eff}}$ for a small selection of white dwarfs. We plan to extend our analysis by allowing metallicity to vary in future model grids.

Modeling the spectral emission of low-charge iron group ions enables the diagnostic determination of the local physical conditions of many cool plasma environments such as those found in H II regions, planetary nebulae, active galactic nuclei etc. Electron-impact excitation drives the population of the emitting levels and, hence, their emissivities. By carrying-out Breit-Pauli and intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of Fe$^{2+}$ which both use the exact same atomic structure and the same close-coupling expansion, we demonstrate the validity of the application of the powerful ICFT method to low-charge iron group ions. This is in contradiction to the finding of Bautista et al. [Ap.J.Lett, 718, L189, (2010)] who carried-out ICFT and Dirac R-matrix calculations for the same ion. We discuss possible reasons.

We present electron collision strengths and their thermally averaged values for the forbidden lines of the astronomically abundant doubly-ionized oxygen ion, O^{2+}, in an intermediate coupling scheme using the Breit-Pauli relativistic terms as implemented in an R-matrix atomic scattering code. We use several atomic targets for the R-matrix scattering calculations including one with 72 atomic terms. We also compare with new results obtained using the intermediate coupling frame transformation method. We find spectroscopically significant differences against a recent Breit-Pauli calculation for the excitation of the [O III] lambda 4363 transition but confirm the results of earlier calculations.

Material ejected during (or immediately following) the merger of two neutron stars may assemble into heavy elements by the r-process. The subsequent radioactive decay of the nuclei can power electromagnetic emission similar to, but significantly dimmer than, an ordinary supernova. Identifying such events is an important goal of future transient surveys, offering new perspectives on the origin of r-process nuclei and the astrophysical sources of gravitational waves. Predictions of the transient light curves and spectra, however, have suffered from the uncertain optical properties of heavy ions. Here we consider the opacity of expanding r-process material and argue that it is dominated by line transitions from those ions with the most complex valence electron structure, namely the lanthanides. For a few representative ions, we run atomic structure models to calculate radiative data for tens of millions of lines. We find that the resulting r-process opacities are orders of magnitude larger than that of ordinary (e.g., iron-rich) supernova ejecta. Radiative transport calculations using these new opacities indicate that the transient emission should be dimmer and redder than previously thought. The spectra appear pseudo-blackbody, with broad absorption features, and peak in the infrared (~1 micron). We discuss uncertainties in the opacities and attempt to quantify their impact on the spectral predictions. The results have important implications for observational strategies to find and study the radioactively powered electromagnetic counterparts to compact object mergers.

We have developed a general model for determining density-dependent effective dielectronic recombination (DR) rate coefficients in order to explore finite-density effects on the ionization balance of plasmas. Our model consists of multiplying by a suppression factor those highly-accurate total zero-density DR rate coefficients which have been produced from state-of-the-art theoretical calculations and which have been benchmarked by experiment. The suppression factor is based-upon earlier detailed collision-radiative calculations which were made for a wide range of ions at various densities and temperatures, but used a simplified treatment of DR. A general suppression formula is then developed as a function of isoelectronic sequence, charge, density, and temperature. These density-dependent effective DR rate coefficients are then used in the plasma simulation code Cloudy to compute ionization balance curves for both collisionally ionized and photoionized plasmas at very low (ne = 1 cm^-3) and finite (ne=10^10 cm^-3) densities. We find that the denser case is significantly more ionized due to suppression of DR, warranting further studies of density effects on DR by detailed collisional-radiative calculations which utilize state-of-the-art partial DR rate coefficients. This is expected to impact the predictions of the ionization balance in denser cosmic gases such as those found in nova and supernova shells, accretion disks, and the broad emission line regions in active galactic nuclei.

We have measured electron-ion recombination for Fe XII forming Fe XI using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. The measured merged beams recombination rate coefficient (MBRRC) for collision energies from 0 to 1500 eV is presented. This work uses a new method for determining the absolute MBRRC based on a comparison of the ion beam decay rate with and without the electron beam on. For energies below 75 eV, the spectrum is dominated by dielectronic recombination (DR) resonances associated with 3s-3p and 3p-3d core excitations. At higher energies we observe contributions from 3-N' and 2-N' core excitations DR. We compare our experimental results to state-of-the-art multi-configuration Breit-Pauli (MCBP) calculations and find significant differences, both in resonance energies and strengths. We have extracted the DR contributions from the measured MBRRC data and transformed them into a plasma recombination rate coefficient (PRRC) for temperatures in the range of 10^3 to 10^7 K. We show that the previously recommended DR data for Fe XII significantly underestimate the PRRC at temperatures relevant for both photoionized plasmas (PPs) and collisionaly ionized plasmas (CPs). This is to be contrasted with our MCBP PRRC results which agree with the experiment to within 30% at PP temperatures and even better at CP temperatures. We find this agreement despite the disagreement shown by the detailed comparison between our MCBP and experimental MBRRC results. Lastly, we present a simple parameterized form of the experimentally derived PRRC for easy use in astrophysical modelling codes.

We present results for the electron-impact excitation of highly-charged sulphur ions (S8+ - S11+) obtained using the intermediate-coupling frame transformation R-matrix approach. A detailed comparison of the target structure has been made for the four ions to assess the uncertainty on collision strengths from the target structure. Effective collision strengths (\Upsilon s) are presented at temperatures ranging from 2\times10^2(z+1)^2 K to 2\times10^6(z+1)^2 K (where z is the residual charge of ions. Detailed comparisons for the \Upsilon are made with the results of previous calculations for these ions, which will pose insight on the uncertainty in their usage by astrophysical and fusion modelling codes.

Brown and Mallik (BM) recently showed that, for hot sources, recombination of non-thermal electrons (NTR) onto highly ionised heavy ions is not negligible compared to non-thermal bremsstrahlung (NTB) as a source of flare hard X-rays (HXRs) and so should be included in modelling non-thermal HXR flare emission. In view of major discrepancies between BM results for the THERMAL continua and those of the Chianti code and of RHESSI solar data, we critically re-examine and correct the BM analysis and modify the conclusions concerning the importance of NTR. Although the analytic Kramers expression used by BM is correct for the purely hydrogenic recombination cross section, the heuristic expressions used by BM to extend the Kramers expression beyond the `bare nucleus' case to which it applies had serious errors. BM results have therefore been recalculated using corrected expressions, which have been validated against the results of detailed calculations. At T ~ 10-30 MK the dominant ions are Fe 22+, 23+, 24+ for which BM erroneously overestimated NTR emission by around an order of magnitude. Contrary to the BM claim, NTR in hot flare plasmas does NOT dominate over NTB, although in some cases it can be comparable and so still very important in inversions of photon spectra to derive electron spectra, especially as NTR includes sharp edge features. The BM claim of dominance of NTR over NTB in deka-keV emission is incorrect due to a serious error in their analysis. However, the NTR contribution can still be large enough to demand inclusion in spectral fitting, the spectral edges having potentially serious effects on inversion of HXR spectra to infer fast electron spectra.

Aims. This work reports radiative transition rates and electron impact excitation collision strengths for levels of the 3s23p, 3s3p2, 3s24s, and 3s23d configurations of Siii. Methods. The radiative data were computed using the Thomas-Fermi-Dirac-Amaldi central potential, but with the modifications introduced by Bautista (2008) that account for the effects of electron-electron interactions. We also introduce new schemes for the optimization of the variational parameters of the potential. Additional calculations were carried out with the Relativistic Hartree-Fock and the multiconfiguration Dirac-Fock methods. Collision strengths in LS-coupling were calculated in the close coupling approximation with the R-matrix method. Then, fine structure collision strengths were obtained by means of the intermediate-coupling frame transformation (ICFT) method which accounts for spin-orbit coupling effects. Results. We present extensive comparisons between the results of different approximations and with the most recent calculations and experiment available in the literature. From these comparisons we derive a recommended set of gf- values and radiative transition rates with their corresponding estimated uncertainties. We also study the effects of different approximations in the representation of the target ion on the electron-impact collision strengths. Our most accurate set of collision strengths were integrated over a Maxwellian distribution of electron energies and the resulting effective collision strengths are given for a wide range of temperatures. Our results present significant differences from recent calculations with the B-spline non-orthogonal R-matrix method. We discuss the sources of the differences.

We have measured electron-ion recombination for Fe$^{9+}$ forming Fe$^{8+}$ and for Fe$^{10+}$ forming Fe$^{9+}$ using merged beams at the TSR heavy-ion storage-ring in Heidelberg. The measured merged beams recombination rate coefficients (MBRRC) for relative energies from 0 to 75 eV are presented, covering all dielectronic recombination (DR) resonances associated with 3s->3p and 3p->3d core transitions in the spectroscopic species Fe X and Fe XI, respectively. We compare our experimental results to multi-configuration Breit-Pauli (MCBP) calculations and find significant differences. From the measured MBRRC we have extracted the DR contributions and transform them into plasma recombination rate coefficients (PRRC) for astrophysical plasmas with temperatures from 10^2 to 10^7 K. This spans across the regimes where each ion forms in photoionized or in collisionally ionized plasmas. For both temperature regimes the experimental uncertainties are 25% at a 90% confidence level. The formerly recommended DR data severely underestimated the rate coefficient at temperatures relevant for photoionized gas. At the temperatures relevant for photoionized gas, we find agreement between our experimental results and MCBP theory. At the higher temperatures relevant for collisionally ionized gas, the MCBP calculations yield a Fe XI DR rate coefficent which is significantly larger than the experimentally derived one. We present parameterized fits to our experimentally derived DR PRRC.

Understanding the thermal equilibrium (stability) curve may offer insights into the nature of the warm absorbers often found in active galactic nuclei. Its shape is determined by factors like the spectrum of the ionizing continuum and the chemical composition of the gas. We find that the stability curves obtained under the same set of the above mentioned physical factors, but using recently derived dielectronic recombination rates, give significantly different results, especially in the regions corresponding to warm absorbers, leading to different physical predictions. Using the current rates we find a larger probability of having thermally stable warm absorber at $10^5 \kel$ than previous predictions and also a greater possibility for its multiphase nature. the results obtained with the current dielectronic recombination rate coefficients are more reliable because the warm absorber models along the stability curve have computed coefficient values, whereas previous calculations relied on guessed averages for the same due to lack of available data.

We present results for the inner-shell electron-impact excitation of Fe$^{15+}$ using the intermediate-coupling frame transformation {\it R}-matrix approach in which Auger-plus-radiation damping has been included. The target and close-coupling expansions are both taken to be the 134 levels belonging to the configurations ${\rm 2s^22p^63}l$, ${\rm 2s^22p^53s3}l$, ${\rm 2s^22p^53p^2}$ and ${\rm 2s^22p^53p3d}$. The comparison of Maxwell-averaged effective collision strengths with and without damping shows that the damping reduction is about 30-40% for many transitions at low temperatures, but up to 80% for a few transitions. As a consequence, the results of previous Dirac $R$-matrix calculations (Aggarwal and Keenan, 2008) overestimate the effective collision strengths due to their omission of Auger-plus-radiation damping.

This work reports radiative transition rates and electron impact excitation rate coefficients for levels of the n= 3, 4, 5, 6, 7, 8 configurations of Ca II. The radiative data were computed using the Thomas-Fermi-Dirac central potential method in the frozen core approximation and includes the polarization interaction between the valence electron and the core using a model potential. This method allows for configuration interactions (CI) and relativistic effects in the Breit-Pauli formalism. Collision strengths in LS-coupling were calculated in the close coupling approximation with the R-matrix method. Then, fine structure collision strengths were obtained by means of the intermediate-coupling frame transformation (ICFT) method which accounts for spin-orbit coupling effects. We present extensive comparisons with the most recent calculations and measurements for Ca II as well as a comparison between the core polarization results and the "unpolarized" values. We find that core polarization affects the computed lifetimes by up to 20%. Our results are in very close agreement with recent measurements for the lifetimes of metastable levels. The present collision strengths were integrated over a Maxwellian distribution of electron energies and the resulting effective collision strengths are given for a wide range of temperatures. Our effective collision strengths for the resonance transitions are within ~11% from previous values derived from experimental measurements, but disagree with latter computations using the distorted wave approximation.

We have measured resonance strengths and energies for dielectronic recombination (DR) of Mg-like Fe XV forming Al-like Fe XIV via N=3 -> N' = 3 core excitations in the electron-ion collision energy range 0-45 eV. All measurements were carried out using the heavy-ion Test Storage Ring at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. We have also carried out new multiconfiguration Breit-Pauli (MCBP) calculations using the AUTOSTRUCTURE code. For electron-ion collision energies < 25 eV we find poor agreement between our experimental and theoretical resonance energies and strengths. From 25 to 42 eV we find good agreement between the two for resonance energies. But in this energy range the theoretical resonance strengths are ~ 31% larger than the experimental results. This is larger than our estimated total experimental uncertainty in this energy range of +/- 26% (at a 90% confidence level). Above 42 eV the difference in the shape between the calculated and measured 3s3p(^1P_1)nl DR series limit we attribute partly to the nl dependence of the detection probabilities of high Rydberg states in the experiment. We have used our measurements, supplemented by our AUTOSTRUCTURE calculations, to produce a Maxwellian-averaged 3 -> 3 DR rate coefficient for Fe XV forming Fe XIV. The resulting rate coefficient is estimated to be accurate to better than +/- 29% (at a 90% confidence level) for k_BT_e > 1 eV. At temperatures of k_BT_e ~ 2.5-15 eV, where Fe XV is predicted to form in photoionized plasmas, significant discrepancies are found between our experimentally-derived rate coefficient and previously published theoretical results. Our new MCBP plasma rate coefficient is 19-28% smaller than our experimental results over this temperature range.

We have carried-out multi-configuration Breit-Pauli AUTOSTRUCTURE calculations for the dielectronic recombination (DR) of Fe^{8+} - Fe^{12+} ions. We obtain total DR rate coefficients for the initial ground-level which are an order of magnitude larger than those corresponding to radiative recombination (RR), at temperatures where Fe 3p^q (q=2-6) ions are abundant in photoionized plasmas. The resultant total (DR+RR) rate coefficients are then an order of magnitude larger than those currently in use by photoionized plasma modeling codes such as CLOUDY, ION and XSTAR. These rate coefficients, together with our previous results for q=0 and 1, are critical for determining the ionization balance of the M-shell Fe ions which give rise to the prominent unresolved-transition-array X-ray absorption feature found in the spectrum of many active galactic nuclei. This feature is poorly described by CLOUDY and ION, necessitating an ad hoc modification to the low-temperature DR rate coefficients. Such modifications are no longer necessary and a rigorous approach to such modeling can now take place using these data.

Reliably interpreting spectra from electron-ionized cosmic plasmas requires accurate ionization balance calculations for the plasma in question. However, much of the atomic data needed for these calculations have not been generated using modern theoretical methods and are often highly suspect. This translates directly into the reliability of the collisional ionization equilibrium (CIE) calculations. We make use of state-of-the-art calculations of dielectronic recombination (DR) rate coefficients for the hydrogenic through Na-like ions of all elements from He up to and including Zn. We also make use of state-of-the-art radiative recombination (RR) rate coefficient calculations for the bare through Na-like ions of all elements from H through to Zn. Here we present improved CIE calculations for temperatures from $10^4$ to $10^9$ K using our data and the recommended electron impact ionization data of \citet{Mazz98a} for elements up to and including Ni and Mazzotta (private communication) for Cu and Zn. DR and RR data for ionization stages that have not been updated are also taken from these two additional sources. We compare our calculated fractional ionic abundances using these data with those presented by Mazzotta et al. for all elements from H to Ni. The differences in peak fractional abundance are up to 60%. We also compare with the fractional ionic abundances for Mg, Si, S, Ar, Ca, Fe, and Ni derived from the modern DR calculations of \citet{Gu03a,Gu04a} for the H-like through Na-like ions, and the RR calculations of \citet{Gu03b} for the bare through F-like ions. These results are in better agreement with our work, with differences in peak fractional abundance of less than 10%.

We have carried-out a series of multi-configuration Breit-Pauli AUTOSTRUCTURE calculations for the dielectronic recombination of Fe^{13+}. We present a detailed comparison of the results with the high-energy resolution measurements reported recently from the Heidelberg Test Storage Ring by Schmidt et al. Many Rydberg series contribute significantly from this initial 3s^2 3p M-shell ion, resulting in a complex recombination `spectrum'. While there is much close agreement between theory and experiment, differences of typically 50% in the summed resonance strengths over 0.1-10 eV result in the experimentally based total Maxwellian recombination rate coefficient being a factor of 1.52-1.38 larger than theory over 10^4-10^5 K, which is a typical temperature range of peak abundance for Fe^{13+} in a photoionized plasma. Nevertheless, this theoretical recombination rate coefficient is an order of magnitude larger than that used by modellers to-date. This may help explain the discrepancy between the iron M-shell ionization balance predicted by photoionization modelling codes such as ION and CLOUDY and that deduced from the iron M-shell unresolved-transition-array absorption feature observed in the X-ray spectrum of many active galactic nuclei. Similar data are required for Fe^{8+} through Fe^{12+} to remove the question mark hanging over the atomic data though.

We have calculated partial final-state resolved radiative recombination (RR) rate coefficients from the initial ground and metastable levels of all elements up to and including Zn, plus Kr, Mo, and Xe, for all isoelectronic sequences up to Na-like forming Mg-like. The data are archived according to the Atomic Data and Analysis Structure (ADAS) data class adf48, which spans a temperature range of z^2(10^1 - 10^7)K, where z is the initial ion charge. Fits to total rate coefficients have been determined, for both the ground and metastable levels, and those for the ground are presented here. Comparison is made both with previous RR rate coefficients and with (background) R-matrix photoionization cross sections. This RR database complements a DR database already produced and both are being used to produce updated ionization balances for both (electron) collisionally ionized and photoionized plasmas.

The paper presents an integrated view of the population structure and its role in establishing the ionisation state of light elements in dynamic, finite density, laboratory and astrophysical plasmas. There are four main issues, the generalised collisional-radiative picture for metastables in dynamic plasmas with Maxwellian free electrons and its particularising to light elements, the methods of bundling and projection for manipulating the population equations, the systematic production/use of state selective fundamental collision data in the metastable resolved picture to all levels for collisonal-radiative modelling and the delivery of appropriate derived coefficients for experiment analysis. The ions of carbon, oxygen and neon are used in illustration. The practical implementation of the methods described here is part of the ADAS Project.

Using the code AUTOSTRUCTURE, extensive calculations of inner-shell atomic data have been made for the chemical elements He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Cr, Mn, Fe and Ni. The results are used to obtain up-dated opacities from the Opacity Project, OP. A number of other improvements on earlier work have also been included. Rosseland-mean opacities from OP are compared with those from OPAL. Differences of 5 to 10% occur. OP gives the `Z-bump', at log(T)=5.2, to be shifted to slightly higher temperatures. The opacities from OP, as functions of temperature and density, are smoother than those from OPAL. Extensive tests show that the numerical accuracy of the OP opacities following integration over frequency mesh and interpolation on temperature-density mesh is better than 1%. Prior to a number of recent investigations which have indicated a need for a downward revision in the solar abundances of oxygen and other elements, there was good agreement between properties of the sun deduced from helioseismology and from stellar evolution models calculated using OPAL opacities. The revisions in abundances destroy that agreement. In a recent paper Bahcall et al argue that the agreement would be restored if opacities for the regions of the sun between 0.7 and 0.4 solar radii were larger than those given by OPAL by about 10%. In the region concerned, the present results from OP do not differ from those of OPAL by more than 2.5%.

Monochromatic opacities from the Opacity Project (OP) (Seaton et al.) have been augmented by hitherto missing inner-shell contributions (Badnell & Seaton). OP Rosseland-mean opacities are compared with results from OPAL (Iglesias & Rogers) for the elements H, He, C, O, S and Fe. Overall, there is good agreement between OP and OPAL Rosseland-mean opacities for the 6-elements, but there are some differences. In particular, recent work (Bahcall et al.) has shown that helioseismology measurements give a very accurate value for the depth of the solar convection zone and that solar models give agreement with that value only if opacities there are about 7% larger than OPAL values. For the 6-element mix, here, we obtain Rosseland-mean opacities to be larger than OPAL by 5%.

For high temperatures and densities, stellar opacities obtained from the Opacity Project (OP) were smaller than those obtained from the OPAL project. Iglesias and Rogers [Ap.J. v443, 469 (1995)] suggested that the discrepancy was due to the omission by OP of important atomic inner-shell processes, and considered in detail results for a mixture of 6 elements: H, He, C, O, S and Fe. Extensive new inner-shell data have now been computed using the code AUTOSTRUCTURE. It is shown that the inclusion of these data in the OP work gives opacities for the 6-element mix which are in much closer agreement with those from OPAL. We also discuss a number of problems relating to the calculation of opacities and of equations-of-state in dense plasmas.

A programme is outlined for the assembly of a comprehensive dielectronic recombination database within the generalized collisional--radiative (GCR) framework. It is valid for modelling ions of elements in dynamic finite-density plasmas such as occur in transient astrophysical plasmas such as solar flares and in the divertors and high transport regions of magnetic fusion devices. The resolution and precision of the data are tuned to spectral analysis and so are sufficient for prediction of the dielectronic recombination contributions to individual spectral line emissivities. The fundamental data are structured according to the format prescriptions of the Atomic Data and Analysis Structure (ADAS) and the production of relevant GCR derived data for application is described and implemented following ADAS. The requirements on the dielectronic recombination database are reviewed and the new data are placed in context and evaluated with respect to older and more approximate treatments. Illustrative results validate the new high-resolution zero-density dielectronic recombination data in comparison with measurements made in heavy-ion storage rings utilizing an electron cooler. We also exemplify the role of the dielectronic data on GCR coefficient behaviour for some representative light and medium weight elements.

We have investigated the accuracy of the 1s-vacancy fluorescence data base of Kaastra & Mewe (1993, A&AS, 97, 443) resulting from the initial atomic physics calculations and the subsequent scaling along isoelectronic sequences. In particular, we have focused on the relatively simple Be-like and F-like 1s-vacancy sequences. We find that the earlier atomic physics calculations for the oscillator strengths and autoionization rates of singly-charged B II and Ne II are in sufficient agreement with our present calculations. However, the substantial charge dependence of these quantities along each isoelectronic sequence, the incorrect configuration averaging used for B II, and the neglect of spin-orbit effects (which become important at high-Z) all cast doubt on the reliability of the Kaastra & Mewe data for application to plasma modeling.

We have measured the resonance strengths and energies for dielectronic recombination (DR) of Fe XX forming Fe XIX via N=2 --> N'=2 (Delta_N=0) core excitations. We have also calculated the DR resonance strengths and energies using AUTOSTRUCTURE, HULLAC, MCDF, and R-matrix methods, four different state-of-the-art theoretical techniques. On average the theoretical resonance strengths agree to within <~10% with experiment. However, the 1 sigma standard deviation for the ratios of the theoretical-to-experimental resonance strengths is >~30% which is significantly larger than the estimated relative experimental uncertainty of <~10%. This suggests that similar errors exist in the calculated level populations and line emission spectrum of the recombined ion. We confirm that theoretical methods based on inverse-photoionization calculations (e.g., undamped R-matrix methods) will severely overestimate the strength of the DR process unless they include the effects of radiation damping. We also find that the coupling between the DR and radiative recombination (RR) channels is small. We have used our experimental and theoretical results to produce Maxwellian-averaged rate coefficients for Delta_N=0 DR of Fe XX. For kT>~1 eV, which includes the predicted formation temperatures for Fe XX in an optically thin, low-density photoionized plasma with cosmic abundances, our experimental and theoretical results are in good agreement. We have also used our R-matrix results, topped off using AUTOSTRUCTURE for RR into J>=25 levels, to calculate the rate coefficient for RR of Fe XX. Our RR results are in good agreement with previously published calculations.

In photoionized gases with cosmic abundances, dielectronic recombination (DR) proceeds primarily via nlj --> nl'j' core excitations (Dn=0 DR). We have measured the resonance strengths and energies for Fe XVIII to Fe XVII and Fe XIX to Fe XVIII Dn=0 DR. Using our measurements, we have calculated the Fe XVIII and Fe XIX Dn=0 DR DR rate coefficients. Significant discrepancies exist between our inferred rates and those of published calculations. These calculations overestimate the DR rates by factors of ~2 or underestimate it by factors of ~2 to orders of magnitude, but none are in good agreement with our results. Almost all published DR rates for modeling cosmic plasmas are computed using the same theoretical techniques as the above-mentioned calculations. Hence, our measurements call into question all theoretical Dn=0 DR rates used for ionization balance calculations of cosmic plasmas. At temperatures where the Fe XVIII and Fe XIX fractional abundances are predicted to peak in photoionized gases of cosmic abundances, the theoretical rates underestimate the Fe XVIII DR rate by a factor of ~2 and overestimate the Fe XIX DR rate by a factor of ~1.6. We have carried out new multiconfiguration Dirac-Fock and multiconfiguration Breit-Pauli calculations which agree with our measured resonance strengths and rate coefficients to within typically better than <~30%. We provide a fit to our inferred rate coefficients for use in plasma modeling. Using our DR measurements, we infer a factor of ~2 error in the Fe XX through Fe XXIV Dn=0 DR rates. We investigate the effects of this estimated error for the well-known thermal instability of photoionized gas. We find that errors in these rates cannot remove the instability, but they do dramatically affect the range in parameter space over which it forms.