Models of stellar spectra are necessary for interpreting light from individual stars, planets, integrated stellar populations, nebulae, and the interstellar medium. We provide a comprehensive and homogeneous collection of synthetic spectra for a wide range of atmospheric parameters and chemical compositions. We compile atomic and molecular data from the literature. We adopt the largest and most recent set of ATLAS9 model atmospheres, and use the radiative code ASS$\epsilon$T.The resulting collection of spectra is made publicly available at medium and high-resolution ($R\equiv\lambda/\delta\lambda$ = 10,000, 100,000 and 300,000 spectral grids, which include variations in effective temperature between 3500 K and 30,000 K, surface gravity ($0\le \log g \le 5$), and metallicity ($-5 \le$[Fe/H]$\le +0. 5$), spanning the wavelength interval 120-6500 nm. A second set of denser grids with additional dimensions, [$\alpha$/Fe] and micro-turbulence, are also provided (covering 200-2500 nm). We compare models with observations for a few representative cases.

Phosphorus abundance is crucial for DNA-based extraterrestrial life in exoplanets. Atomic data for observed spectral lines of P-ions are needed for its accurate determination. We present the first calculations for collision strengths for the forbidden PIII fine structure transition $3s^23p (^2P^o_{1/2-3/2})$ within the ground state at 17.9 $\mu$m, as well as allowed UV transitions in the $3s^23p (^2P^o_{1/2,3/2}) \rightarrow 3s3p^2 (^2D_{3/2,5/2}, ^2S_{1/2}, ^2P_{1/2,3/2})$ multiplets between 915-1345 $\AA$. Collision strengths are computed using the Breit-Pauli R-Matrix method including the first 18 levels, and they exhibit extensive auto-ionizing resonance structures. In particular, the Maxwellian averaged effective collision strength for the FIR 17.9 $\mu$m transition shows a factor 3 temperature variation broadly peaking at typical nebular temperatures. Its theoretical emissivity with solar phosphorus abundance is computed relative to H$\beta$ and found to be similar to observed intensties from planetary nebulae; the abundances derived in earlier works are 3-5 times sub-solar. The results pertain to the reported paucity of phosphorus from preferred production sites in supernovae, and abundances in planetary nebulae and supernova remnants.

The IRON Project, initiated in 1991, aims at two main objectives, i) study the characteristics of and calculate large-scale high accuracy data for atomic radiative and collisional processes, and ii) application in solving astrophysical problems. It focuses on the complex iron and iron-peak elements commonly observed in the spectra of astrophysical plasmas. The present report will illustrate the characteristics of the dominant atomic process of photoionization that have been established under the project and the preceding the Opacity Project and their importance in applications.

Extensive resonance structures are manifest in R-Matrix (RM) calculations. However, there exist a large number of highly excited electronic configurations that may contribute to background non-resonant bound-free opacity in high-temperature plasmas. Since RM calculations are very complex, and not essential for background contributions, the Relativistic Distorted Wave (RDW) method is utilized to complement ("top-up") photoionization cross sections of Fe XVII obtained using Close-Coupling Breit-Pauli R-Matrix (CC-BPRM) method. There is good agreement between RDW and BPRM for background cross sections where resonances are not present, and individual fine structure levels can be correctly matched spectroscopically, though resonances are neglected in the RDW. To ensure completeness, a high energy range up to 500 Ry above the ionization threshold for each level is considered. Interestingly, the hydrogenic Kramer's approximation also shows the same energy behavior as the RDW. Grouping separately, the BPRM configurations consist of 454 bound levels with resonances corresponding to configurations $1s^22s^22p^4nln'l'$ (n $\leq$ 3, n' $\leq$ 10); including RDW configurations there are 51,558 levels in total. The topup contribution results in $\sim$20\% increment, in addition to the 35\% enhancement from BPRM calculations over the Opacity Project value for the Rosseland Mean Opacity at the Z-temperature of 2.11 $\times 10^6$K (Pradhan and Nahar 2017).

A review of a renewed effort to recalculate astrophysical opacities using the R-Matrix method is presented. The computational methods and new extensions are described. Resulting enhancements found in test calculations under stellar interior conditions compared to the Opacity Project could potentially lead to the resolution of the solar abundances problem, as well as discrepancies between recent experimental measurements at the Sandia Z-pinch inertial confinement fusion device and theoretical opacity models. Outstanding issues also discussed are: (i) accuracy, convergence, and completeness of atomic calculations, (ii) improvements in the Equation-of-State of high-temperature-density plasmas, and (iii) redistribution of resonant oscillator strength in the bound-free continuum, and (iv) plasma broadening of auotionizing resonances.

Aimed at solving the outstanding problem of solar opacity, and radiation transport plasma models in general, we report substantial photoabsorption in the high-energy regime due to atomic core photo-excitations not heretofore considered. In extensive R-Matrix calculations of unprecedented complexity for an important iron ion Fe xvii (Fe$^{16+}$), with a wave function expansion of 99 Fe xviii (Fe$^{17+}$) LS core states from $n \leq 4$ complexes (equivalent to 218 fine structure levels), we find: i) up to orders of magnitude enhancement in background photoionization cross sections, in addition to strongly peaked photo-excitation-of-core resonances not considered in current opacity models, and ii) demonstrate convergence with respect to successive core excitations. The resulting increase in the monochromatic continuum, and 35% in the Rosseland Mean Opacity, are compared with the "higher-than-predicted" iron opacity measured at the Sandia Z-pinch fusion device at solar interior conditions.

Improved collisions strengths for the mid-infrared and optical transitions in Ne V are presented. Breit-Pauli R-Matrix calculations for electron impact excitation are carried out with fully resolved near-threshold resonances at very low energies. In particular, the fine structure lines at 14 micron and 24 micron due to transitions among the ground state levels 1s^22s^22p^3 (^3P_{0,1,2}), and the optical/near-UV lines at 2973, 3346 and 3426 Angstrom transitions among the ^3P_{0,1,2}, ^1D_2, ^1S_0 levels are described. Maxwellian averaged collision strengths are tabulated for all forbidden transistion within the ground configuration. Significant differences are found in the low temperature range Te < 10000 K for both the FIR and the opitcal transitions compared to previous results. An analysis of the 14/24 line ratio in low-energy-density (LED) plasma conditions reveals considerable variation; the effective rate coefficient may be dominated by the very low-energy behaviour rather than the maxwellian averaged collision strengths. Computed values suggest a possible solution to the anomalous mid-IR ratios found to be lower than theoretical limits observed from planetary nebulae and Seyfert galaxies. While such LED conditions may be present in infrared sources, they might be inconsistent with photoionization equilibrium models.

Far-infrared and optical [O III] lines are useful temeprature-density diagnostics of nebular as well as dust obscured astrophysical sources. Fine structure transitions among the ground state levels 1s^22s^22p^3 \ ^3P_{0,1,2} give rise to the 52 and 88 micron lines, whereas transitions among the $^3P_{0,1,2}, ,^1D_2, ^1S_0$ levels yield the well-known optical lines 4363, 4959 and 5007 Angstroms. These lines are excited primarily by electron impact excitation. But despite their importance in nebular diagnostics collision strengths for the associated fine structure transitions have not been computed taking full account of relativistic effects. We present Breit-Pauli R-matrix calculations for the collision strengths with highly resolved resonance structures. We find significant differences of up to 20% in the Maxwellian averaged rate coefficients from previous works. We also tabulate these to lower temperatures down to 100 K to enable determination of physical conditions in cold dusty environments such photo-dissociation regions and ultra-luminous infrared galaxies observed with the Herschel space observatory. We also examine the effect of improved collision strengths on temperature and density sensitive line ratios.

A comprehensive study of high-accuracy photoionization cross sections is carried out using the relativistic Breit-Pauli R-matrix (BPRM) method for (hnu + Fe XVII --> Fe XVIII + e). Owing to its importance in high-temperature plasmas the calculations cover a large energy range, particularly the myriad photoexciation-of-core (PEC) resonances including the n = 3 levels not heretofore considered. The calculations employ a close coupling wave function expansion of 60 levels of the core ion Fe XVIII ranging over a wide energy range of nearly 900 eV between the n = 2 and n = 3 levels. Strong coupling effects due to dipole transition arrays 2p^5 --> 2p^4 (3s,3d) manifest themselves as large PEC resonances throughout this range, and enhance the effective photoionization cross sections orders of magnitude above the background. Comparisons with the erstwhile Opacity Project (OP) and other previous calculations shows that the currently available cross sections considerably underestimate the bound-free cross sections. A level-identification scheme is used for spectroscopic designation of the 454 bound fine structure levels of Fe XVII. Level-specific photoionization cross sections are computed for all levels. In addition, partial cross sections for leaving the core ion Fe XVII in the ground state are also obtained. These results should be relevant to modeling of astrophysical and laboratory plasma sources requiring (i) photoionization rates, (ii) extensive non-local-thermodynamic-equilibrium models, (iii) total unified electron-ion recombination rates including radiative and dielectronic recombination, and (iv) plasma opacities. We particularly examine PEC and non-PEC resonance strengths and emphasize their expanded role to incorporate inner-shell excitations for improved opacities, as shown by the computed monochromatic opacity of Fe XVII.

The unified method for total electron-ion recombination is extended to study the dielectronic satellite (DES) lines. These lines, formed from radiative decay of autoionizing states, are highly sensitive temperature diagnostics of astrophysical and laboratory plasma sources. The computation of the unified recombination rates is based on the relativistic Breit-Pauli R-matrix method and close coupling approximation. Extending the theoretical formulation developed earlier we present recombination rate coefficients for the 22 satellite lines of KLL complexes of helium-like Fe XXV and Ni XXVII. The isolated resonance approximation, commonly used throughout plasma modeling, treats these resonances essentially as bound features except for dielectronic capture into, and autoionization out of, these levels. A line profile or cross section shape is often assumed. On the other hand, by including the coupling between the autoionizing and continuum channels, the unified method gives the intrinsic spectrum of DES lines which includes not only the energies and strengths, but also the natural line or cross section shapes. A formulation is presented to derive autoionization rates from unified resonance strengths and enable correspondence with the isolated resonance approximation. While the rates compare very well with existing rates for the strong lines to <20%, the differences for weaker DES lines are larger. We also illustrate the application of the present results to the analysis of K ALPHA complexes observed in high-temperature X-ray emission spectra of Fe XXV and Ni XXVII. There are considerable differences with previous results in the total KLL intensity for Fe XXV at temperatures below the temperature of maximum abundance in coronal equilibrium. (Abbreviated Abstract)

New, more detailed calculations of the emission spectra of the He-like K ALPHA complex of calcium, iron and nickel have been carried out using data from both distorted-wave and R-matrix calculations. The value of the GD ratio (an extended definition of the G ratio that accounts for the effect of resolved and unresolved satellite lines) is significantly enhanced at temperatures below the temperature of He-like maximum abundance. Furthermore it is shown that satellite lines are important contributors to the GD ratio such that GD/G>1 at temperatures well above the temperature of maximum abundance. These new calculations demonstrate, with an improved treatment of the KLn (n>=3) satellite lines, that K ALPHA satellite lines need to be included in models of He like spectra even at relatively high temperatures. The excellent agreement between spectra and line ratios calculated from R-matrix and distorted-wave data also confirms the validity of models based on distorted-wave data for highly charged systems, provided the effect of resonances are taken into account as independent processes.

It is shown that above the temperature of maximum abundance, recombination rates into the excited states of He-like ions that are calculated using earlier, more approximate methods differ markedly from rates obtained from recent distorted-wave and R-Matrix calculations (unified recombination rate coefficients) for Ca, Fe, and Ni. The present rates lead to G ratios that are greatly lower than those resulting from the more approximate rates in previous works, by up to a factor of six at high electron temperatures. Excellent agreement between the distorted-wave and the R-Matrix rates, as well as excellent agreement in the G ratios calculated from them, provides support for the accuracy of these new values which have a broad applicability to the modelling and interpreting of X-ray spectra from a variety of astrophysical and laboratory sources.

Based on new calculations we reconfirm the low and high density limits on the forbidden fine structure line ratio [O II] I(3729)/I(3726): lim_{N_ e} --> 0} = 1.5 and lim_{N_ e} --> \infty} = 0.35. Employing [O II] collision strengths calculated using the Breit-Pauli R-matrix method we rule out any significant deviation due to relativistic effects from these canonical values. The present results are in substantial agreement with older calculations by Pradhan (1976) and validate the extensive observational analysis of gaseous nebulae by Copetti and Writzel (2002) and Wang et al (2004) that reach the same conclusions. The present theoretical results and the recent observational analyses differ significantly from the calculations by MacLaughlin and Bell (1998) and Keenan et al (1999). The new maxwellian averaged effective collision strengths are presented for the 10 transitions among the first 5 levels to enable computations of [O II] line ratios.

The inverse processes of photoionization and electron-ion recombination of (hnu + Ni XXVI --> Ni XXVII + e and (hnu + Ni XXVII --> Ni XXVIII + e) are studied using the unified method for the total recombination. The method subsumes both the radiative and di-electronic recombination processes and enables self-consistent sets of results for photoionization and electron-ion recombination by using the same wavefunction expansion. Photoionization cross sections (sigma_PI), recombination cross sections (sigma_{RC}), recombination collision strengths (Omega_{RC}), and recombination rate coefficients (alpha_{RC}) are obtained for ionization balance and spectral analysis of UV and X-ray lines. Level-specific photoionization cross sections and recombination rates are presented to enable accurate computation of recombination-cascade matrices for all fine structure levels n(SLJ) up to n <= 10: 98 bound fine structure levels of Ni XXVI with 0 <= l <= 9, 0 <= L <= 11, 1/2 <= J <= 17/2, and 198 levels of Ni XXVII with 0 <= l <= 9, 0 <= L <= 14, 0 <= J <= 10. Total alpha_{RC} for Ni XXVI and Ni XXVII are compared with the existing values with very good agreement. Total recombination rate coefficients for the hydrogen-like recombined ion, Ni XXVIII, are also presented. The calculations are carried out in relativistic Breit-Pauli R-matrix (BPRM) approximation with inclusion of radiation damping of resonances. With consideration of all details of the processes, the results, which include the level specific sigma_{PI} and alpha_{R} calculated for the first time, should be the most accurate for these ions.

Ar XIII is found to be unique with new features in electron-ion recombination not seen in any other ion. The ion has been studied with the unified method which provides a theoretically self-consistent set of atomic parameters for the inverse processes of photoionization and total electron-ion recombination. The unified method subsumes both the radiative recombination (RR) and dielectronic recombination (DR) within the framework of close-coupling formulations using the R-matrix method. A set of four DR "bumps", two in the low and two in the high temperature regions, is found to exist in the recombination rates of Ar XIV + e $\to$ Ar XIII. This is in contrast to two typical DR "bumps", one at high temperature common for most ions and one at low temperature depending on the presence of near threshold autoionizing resonances in the bound-free process. The large scale ab initio calculations for Ar XIII employ eigenfunction expansion of 37 core Ar XIV states from n = 2 and 3 complexes. This enables core excitations of type $\Delta n$ = 0 and 1. Ar XIII exhibits more extensive resonant structures and wider PEC (photoexitation-of-core) resonances for n = 3 core states than those of n = 2 states. Hence the high energy region of photoionization and recombination are dominated by these structures. Total and partial photoionization cross sections of all bound states, state-specific recombinaion rates of 561 bound states and total recombination rate coefficients are presented.

We present new calculations for transition probabilities of Fe IV, with much more extensive datasets than heretofore available. The large-scale close coupling R-matrix calculations yield 1,798 LS bound states with n $\leq$ 11 and $l \leq$ 9, and corresponding transition probabilities for 138,121 dipole allowed transitions in the form of oscillator strengths $f$, line strengths $S$, and A-coeffficients for a variety of applications. This represents the largest R-matrix dataset in LS coupling for any ion under either the Opacity Project or the Iron Project. Through algebraic transformation of the LS multiplets, a total of 712,120 dipole allowed fine structure transitions for Fe IV are obtained. Observed transition energies are used together with the energy independent line strengths to derive the $f$- and the $A$-values; the adopted algorithm used calculated energies for the remainder. Present results show significantly better accuracy for the important low-lying states than previous calculations. Monochromatic and mean opacities for Fe IV are computed and compared with those obtained using the Opacity Project data. We find differences which could have important consequences for several astrophysical applications involving low ionization stages of iron.

We investigate the effect of resonances in photoionization and recombination cross sections computed using the R-matrix method. Photoionization and recombination rates derived from high-resolution cross sections for oxygen ions are compared with earlier works with less resolution and accuracy, such as in the widely used Opacity Project data. We find significant differences in photoionization rates for O II metastable states, averaged over Planck functions corresponding to ionizing radiation fields, with respect to the intrinsic accuracy of the calculations and improved resolution. Furthermore, for highly charged ions other physical effects are also important. Recombination rate coefficients, averaged over a Maxwellian distribution, are extremely sensitive to the position and resolution of near-threshold resonances, and radiation damping, in (e + O VII) --> O VI + hnu. Surprisingly however, the effect on the monochromatic and the mean Rosseland and Planck bound-free opacities is relatively small, but may be potentially significant.

We present theoretical iron emission line strengths for physical conditions typical of Active Galactic Nuclei with Broad-Line Regions. The non-local thermodynamic equilibrium (NLTE) models include a new and extensive treatment of radiative transfer in the FeIII ion, complementing the FeII emission line strengths predicted in our earlier works. We also briefly present preliminary results for the FeI emission from AGN using a reduced atom model. We can satisfactorily reproduce the empirical UV FeIII emission line template of Vestergaard & Wilkes (2001) for the prototypical narrow-line Seyfert 1 galaxy I Zw 1, both in terms of the general FeIII flux distribution and the relative strength of the FeIII and FeII emission. However, a number of detailed features are still not matched; the most prominent example is the strongest single FeIII feature observed in the I Zw 1 spectrum, UV47: it is predicted to be strong only in models suppressing Fe-H charge exchange reactions. We examine the role of variations in cloud turbulent velocity and iron abundance and carry out Monte Carlo simulations to demonstrate the effect of uncertainties in atomic data on the computed spectra.

For crucial tests between theory and experiment, ab initio close coupling calculations are carried out for photoionization of O II, O III, O IV, O V. The relativistic fine structure and resonance effects are studied using the R-matrix and its relativistic variant the Breit Pauli R-matrix (BPRM) approximation. Detailed comparison is made with high resolution experimental measurements carried out in three different set-ups: Advanced Light Source at Berkeley, and synchrotron radiation experiments at University of Aarhus and University of Paris-Sud. The comparisons illustrate physical effects in photoionization such as (i) fine structure, (ii) resolution, and (iii) metastable components. Photoionization cross sections sigma{PI} of the ground and a few low lying excited states of these ions obtained in the experimental spectrum include combined features of these states. Theoretically calculated resonances need to be resolved with extremely fine energy mesh for precise comparison. In addition, prominent resonant features are observed in the measured spectra from transitions allowed with relativistic fine structure, but not in LS coupling. The sigma_{PI} are obtained for ground and metastable (i) 2s^22p^3(^4S^o, ^2D^o, ^2P^o) states of O II, (ii) 2s^22p^2(^3P,^1D,^1S) and 2s2p^3(^5S^o) states of O III, (iii) 2s^22p(^2P^o_J) and 2s2p^2(^4P_J) levels of O IV, and (iv) 2s^2(^1S) and 2s2p(^3P^o,^1P^o) states of O V. It is found that resonances in ground and metastable cross sections can be a diagnostic of experimental beam composition, with potential ap plications to astrophysical and laboratory plasma environments.

A unified scheme using the R-matrix method has been developed for electron-ion recombination subsuming heretofore separate treatments of radiative and dielectronic recombination (RR and DR). The ab initio coupled channel approach unifies resonant and non-resonant phenomena, and enables a general and self-consistent treatment of photoionization and electron-ion recombination employing idential wavefunction expansion. Detailed balance takes account of interference effects due to resonances in cross sections, calculated explicitly for a large number of recombined (e+ion) bound levels over extended energy regions. The theory of DR by Bell and Seaton is adapted for high-n resonances in the region below series limits. The R-matrix method is employed for (A) partial and total photoionization and photorecombination cross sections of (e+ion) bound levels, and (B) DR and (e+ion) scattering cross sections. Relativistic effects and fine structure are considered in the Breit-Pauli approximation. Effects such as radiation damping may be taken into account where necessary. Unfiied recombination cross sections are in excellent agreement with measurements on ion storage rings to about 10-20%. In addition to high accuracy, the strengths of the method are: (I) both total and level-specific cross sections and rate coefficients are obtained, and (II) a single (e+ion) recombination rate coefficient for any given atom or ion is obtained over the entire temperature range of practical importance in laboratory and astrophysical plasmas, (III) self-consistent results are obtained for photoionization and recombination; comprehensive datasets have been computed for over 50 atoms and ions. Selected data are presented for iron ions.

Aimed at ionization balance and spectral analysis of UV and X-ray sources, we present self-consistent sets of photoionization cross sections, recombination cross sections, and rate coefficients for Li-like O VI and He-like O VII. Relativistic fine structure is considered through the Breit-Pauli R-matrix (BPRM) method in the close coupling approximation, implementing the unified treatment for total electron-ion recombination subsuming both radiative and di-electronic recombination processes. Self-consistency is ensured by using an identical wavefunction expansion for the inverse processes of photoionization and photo-recombination. Radiation damping of resonances, important for H-like and He-like core ions, is included. Compared to previous LS coupling results without radiative decay of low-n (<= 10) resonances, the presents results show significant reduction in O VI recombination rates at high temperatures. In addition to the total rates, level-specific photoionization cross sections and recombination rates are presented for all fine structure levels n (lSLJ) up to n <= 10, to enable accurate computation of recombination-cascade matrices and spectral formation of prominent UV and X-ray lines such as the 1032,1038 A doublet of O VI, and the `triplet' forbidden, intercombination, and resonance X-ray lines of O VII at 22.1, 21.8, and 21.6 \ang respectively. Altogether, atomic parameters for 98 levels of O VI and 116 fine structure levels of O VII are theoretically computed. These data should provide a reasonably complete set of photoionization and recombination rates in collisional or radiative equilibrium.

An extensive set of fine structure levels and corresponding transition probabilities for allowed and forbidden transitions in Fe XVII is presented. A total of 490 bound energy levels of Fe XVII of total angular momenta 0 <= J <= 7 of even and odd parities with 2 <= n <= 10, 0 <= l <= 8, 0 <= L <= 8, and singlet and triplet multiplicities, are obtained. They translate to over 2.6 x 10^4 allowed (E1) transitions that are of dipole and intercombination type, and about 3000 forbidden transitions that include electric quadrupole (E2), magnetic dipole (M1), electric octopole (E3), and magnetic quadrupole (M2) type representing the most detailed calculations to date for the ion. Oscillator strengths f, line strengths S, and coefficients A of spontaneous emission for the E1 type transitions are obtained in the relativistic Breit-Pauli R-matrix approximation. A valus for the forbidden transitions are obtained from atomic structure calculations using codes SUPERSTRUCTURE and GRASP. The energy levels are identified in spectroscopic notation with the help of a newly developed level identification algorithm. Nearly all 52 spectroscopically observed levels have been identified, their binding energies agreeing within 1% with our calculation. Computed transition probabilities are compared with other calculations and measurement. The effect of 2-body magnetic terms and other interactions is discussed. Present data set enhances by more than an order of magnitude the heretofore available data for the transition probabilities of Fe XVII.

The K-alpha resonance complexes in oxygen ions O I - O VI are theoretically computed and resonance oscillator strengths and wavelengths are presented. The highly resolved photoionization cross sections, with relativistic fine structure, are computed in the coupled channel approximation using the Breit-Pauli R-matrix method. A number of strong K-alpha resonances are found to be appreciable, with resonance oscillator strengths f_r > 0.1. The K-alpha resonance wavelengths of O I-O VI lie in a relatively narrow wavelength range 22 - 23.5 A, and the X-ray opacity in this region should therefore be significantly affected by K --> L transitions in oxygen. The results should be useful in the interpretation of soft X-ray spectra observed from Chandra and XMM-Newton.

Ab initio calculations including relativistic effects employing the Breit-Pauli R-matrix (BPRM) method are reported for fine structure energy levels and oscillator strengths upto n = 10 and 0.leq. l .leq.9 for 15 Li-like ions: C IV, N V, O VI, F VII, Ne VIII, Na IX, Mg X, Al XI, Si XII, S XIV, Ar XVI, Ca XIII, Ti XX, Cr XXII, and Ni XXVI. About one hundred bound fine structure energy levels of total angular momenta, 1/2 .leq. J .leq. 17/2 of even and odd parities, total orbital angular momentum, 0 .leq L .leq. 9 and spin multiplicity (2S+1) = 2, 4 are considered for each ion. The levels provide almost 900 dipole allowed and intercombination bound-bound transitions. The BPRM method enables consideration of large set of transitions with uniform accuracy compared to the best available theoretical methods. The CC eigenfunction expansion for each ion includes the lowest 17 fine structure energy levels of the core configurations 1s^2, 1s2s, 1s2p, 1s3s, 1s3p, and 1s3d. The calculated energies of the ions agree with the measured values to within 1% for most levels. The transition probabilities show good agreement with the best available calculated values. The results provide the largest sets of energy levels and transition rates for the ions and are expected to be useful in the analysis of X-ray and EUV spectra from astrophysical sources.

Latest developments in theoretical computations since the international Opacity Project (OP), under the new the Iron Project (IP) and extensions, are described for applications to a variety of objects such as stellar atmospheres, nebulae, and active galactic nuclei. The primary atomic processes are: electron impact excitation (EIE), photoionization, electron-ion recombination, and bound-bound transitions, all considered using the accurate and powerful R-matrix method including relativistic effects. As an extension of the OP and the IP, a self-consistent and unified theoretical treatment of photoionization and recombination has been developed. Both the radiative and the dielectronic recombination (RR and DR) processes are considered in a unified manner. Photoionization and recombination cross sections are computed with identical wavefunction expansions, thus ensuring self-consistency in an ab initio manner. The new unified results differ from the sum of previous results for RR and DR by up to a factor of 4 for the important but complex atomic systems such as Fe I - V. The fundamental differences are due to quantum mechanical intereference and coupling effects neglected in simpler approximations that unphysically treat RR and DR separately, which can not be independently measured or observed. I will also describe the electronic, web-interactive, database, TIPTOPBASE, to archive the OP/IP data in a readily accessible manner. TIPTOPBASE would also include electron-ion recombination data and new fine structure transition probabilities. Efficient codes developed by M.J. Seaton to calculate `customized' mixture opacities and radiative accelerations ('levitation') in stars will also be available.

Most astrophysical plasmas entail a balance between ionization and recombination. We present new results from the powerful R-matrix method which yields in an ab initio manner: (I) self-consistent photoionization and recombination cross sections using identical wavefunction expansions, (II) unified e-ion recombination rates at all temperatures of interest, incorporating non-resonant and resonant processes, radiative and dielectronic recombination (RR and DR), and (III) level-specific recombination rates for many excited atomic levels. RR and DR processes can not be measured or observed independently. In contrast to the shortcomings of simple (but computationally easy) approximations that unphysically treat RR and DR separately with different methods, emphasizing marginal effects for selected ions over small energy ranges, the R-matrix method naturally and compeletely accounts for e-ion recombination for all atomic systems. Photoionization and recombination cross sections are compared with state-of-the-art experiments on synchrotron radiation sources and ion storage rings. Overall agreement between theory and experiments is within 10-20 %. For photoionization the comparison includes not only the ground state but also the metastable states, with highly resolved resonance structures. The recent experiments therefore support the estimated accuracy of the vast amount of photoionization data computed under the Opacity Project (OP), the Iron Project (IP), and related works using the R-matrix method.

Theoretical differential oscillator strengths related to monochromatic opacity due to K- and L-shell absorption from oxygen and iron ions are directly compared with the Chandra and Xmm-Newton spectra of Seyfert~1 galaxy MCG--6-30-15. We compute the highly resolved continuum with resonances due to O I - VI and Fe XVI. It is found that the KLn (n --> infty) resonance series limits in O VI, and the lowest LMM resonance in Fe XVI, lie at the prominent 17.5 A break in the observed spectra. We also calculate and identify, for the first time, the observed gap in spectral flux at 22-23 A in the Chandra spectra due to K --> L resonant absorption features from all O-ions O I - VI, potentially leading to abundance determination of several or all ionization stages; similar signature gaps may be predicted due to other elements. The precise atomic parameters are computed in the relativistic close coupling approximation using the Breit-Pauli R-matrix method. The new X-ray opacities might possibly distinguish between models of a dusty warm absorber, and/or gravitational redshift and broadening due to a massive black hole in MCG--6-30-15.

A hot central star illuminating the surrounding ionized H II region usually produces very rich atomic spectra resulting from basic atomic processes: photoionization, electron-ion recombination, bound-bound radiative transitions, and collisional excitation of ions. Precise diagnostics of nebular spectra depend on accurate atomic parameters for these processes. Latest developments in theoretical computations are described, especially under two international collaborations known as the Opacity Project (OP) and the Iron Project (IP), that have yielded accurate and large-scale data for photoionization cross sections, transition probabilities, and collision strengths for electron impact excitation of most astrophysically abundant ions. As an extension of the two projects, a self-consistent and unified theoretical treatment of photoionization and electron-ion recombination has been developed where both the radiative and the dielectronic recombination processes are considered in an unified manner. Results from the Ohio State atomic-astrophysics group, and from the OP and IP collaborations, are presented. A description of the electronic web-interactive database, TIPTOPBASE, with the OP and the IP data, and a compilation of recommended data for effective collision strengths, is given.

Relativistic close coupling calculations are reported for unified electronic recombination of (e + C IV) including non-resonant and resonant recombination processes, radiative and dielectronic recombination (RR and DR). Detailed comparison of the theoretical unified results with two recent experiments on ion storage rings (Mannervik etal. [1] and Schippers etal. [2]) shows very good agreement in the entire measured energy region 2s -- 2p with 2pnl resonances. The results benchmark theory and experiments to uncertainties of ~15%, and show that the resonant and the background cross sections are not an incoherent sum of separate RR and DR contributions. The limiting values of the DR cross sections, as n ---> infinity, are shown to correspond to those due to electron impact excitation (EIE) at the 2Po_(1/2,3/2) fine structure thresholds, delineated for the first time. The near-threshold 2s 2S_(1/2) - 2p 2Po_(1/2,3/2) EIE cross sections are also compared with recent experimental measurements. The demonstrated threshold fine structure and resonance effects should be of general importance in excitation and recombination of positive ions.

High resolution measurements of photoionization cross sections for atomic ions are now being made on synchrotron radiation sources. The recent measurements by Kjeldsen etal. (1999) showed good agreement between the observed resonance features and the the theoretical calculations in the close coupling approximation (Nahar 1995). However, there were several observed resonances that were missing in the theoretical predictions. The earlier theoretical calculation was carried out in LS coupling where the relativistic effects were not included. Present work reports photoionization cross sections including the relativistic effects in Breit-Pauli R-matrix (BPRM) approximation. The configuration interaction eigenfunction expansion for the core ion C III consists of 20 fine structure levels dominated by the configurations from 1s^22s^2 to 1s^22s3d. Detailed features in the calculated cross sections exhibit the missing resonances due to fine structure. The results benchmark the accuracy of BPRM photoionization cross sections as needed for recent and ongoing experiments.

Recently X-ray photoabsorption in KLL resonances of O VI was predicted [Pradhan, Astrophys.J. Lett. 545, L165 (2000)], and detected by the Chandra X-ray Observatory [Lee et al, Astrophys. J. {\it Lett.}, submitted]. The required resonance oscillator strengths f_r, are evaluated in terms of the differential oscillator strength df/de that relates bound and continuum absorption. We present the f_r values from radiatively damped and undamped photoionization cross sections for Li-like C,O, and Fe calculated using relativistic close coupling Breit-Pauli R-matrix method. The KLL resonances of interest here are: 1s2p (^3P^o) 2s [^4P^o_{1/2,3/2}, ^2P^o_{1/2,3/2}] and 1s2p (^1P^o) 2s [^2P^o_{1/2,3/2}]. The KLL photoabsorption resonances in Fe XXIV are fully resolved up to natural autoionization profiles for the first time. It is demonstrated that the undamped f_r independently yield the resonance radiative decay rates, and thereby provide a precise check on the resolution of photoionization calculations in general. The predicted photoabsorption features should be detectable by the X-ray space observatories and enable column densities in highly ionized astrophysical plasmas to be determined from the calculated f_r. The dielectronic satellites may appear as redward broadening of resonances lines in emission and absorption.

Relativistic and channel coupling effects in photoionization and unified electronic recombination of Fe XVII are demonstrated with an extensive 60-level close coupling calculation using the Breit-Pauli R-matrix method. Photoionization and (e + ion) recombination calculations are carried out for the total and the level-specific cross sections, including the ground and several hundred excited bound levels of Fe XVII (up to fine structure levels with n = 10). The unified (e + ion) recombination calculations for (e + Fe XVIII --> Fe XVII) include both the non-resonant and resonant recombination (`radiative' and `dielectronic recombination' -- RR and DR). The low-energy and the high energy cross sections are compared from: (i) a 3-level calculation with 2s^2p^5 (^2P^o_{1/2,3/2}) and 2s2p^6 (^2S_{1/2}), and (ii) the first 60-level calculation with \Delta n > 0 coupled channels with spectroscopic 2s^2p^5, 2s2p^6, 2s^22p^4 3s, 3p, 3d, configurations, and a number of correlation configurations. Strong channel coupling effects are demonstrated throughout the energy ranges considered, in particular via giant photoexcitation-of-core (PEC) resonances due to L-M shell dipole transition arrays 2p^5 --> 2p^4 3s, 3d in Fe XIII that enhance effective cross sections by orders of magnitude. Comparison is made with previous theoretical and experimental works on photoionization and recombination that considered the relatively small low-energy region (i), and the weaker \Delta n = 0 couplings. While the 3-level results are inadequate, the present 60-level results should provide reasonably complete and accurate datasets for both photoionization and (e + ion) recombination of Fe~XVII in laboratory and astrophysical plasmas.

We present the first detailed ab initio quantum mechanical calculations for total and state-specific recombination rate coefficients for e + Ni III --> Ni II. These rates are obtained using a unified treatment for total electron-ion recombination that treats the nonresonant radiative recombination and the resonant dielectronic recombination in a self-consistent unified manner in the close coupling approximation. Large-scale calculations are carried out using a 49-state wavefunction expansion from core configurations 3d^8, 3d^74s, and 3d^64p that permits the inclusion of prominent dipole allowed core transitions. These extensive calculations for the recombination rates of Ni II required hundreds of CPU hours on the Cray T90. The total recombination rate coefficients are provided for a wide range of temperature. The state-specific recombination rates for 532 bound states of doublet and quartet symmetries, and the corresponding photoionization cross sections for leaving the core in the ground state, are presented. Present total recombination rate coefficients differ considerably from the currently used data in astrophysical models.

Photoionization and recombination cross sections and rate coefficients are calculated for Li-like Fe XXIV and He-like Fe XXV using the Breit-Pauli R-matrix (BPRM) method. A complete set of total and level-specific parameters is obtained to enable X-ray photoionization and spectral modeling. The ab initio calculations for the unified (e + ion) recombination rate coefficients include both the non-resonant and the resonant recombination (radiative and di-electronic recombination, RR and DR, respectively) for (e + Fe XXV) -> Fe XXIV and (e + Fe XXVI) -> Fe XXV. The level specific rates are computed for all fine structure levels up to n = 10, enabling accurate computation of recombination-cascade matrices and effective rates for the X-ray lines. The total recombination rate coefficients for both Fe XXIV and Fe XXV differ considerably, by several factors, from the sum of RR and DR rates currently used to compute ionization fractions in astrophysical models. As the photoionization/recombination calculations are carried out using an identical eigenfunction expansion, the cross sections for both processes are theoretically self-consistent; the overall uncertainty is estimated to be about 10-20%. All data for Fe XXIV and Fe XXV (and also for H-like Fe XXVI, included for completeness) are available electronically.

Unified recombination cross sections and rates are computed for (e + Fe XVIII) --> Fe XVII including non-resonant and resonant (radiative and di-electronic recombination, RR and DR) processes in an ab initio manner with relativistic fine structure. The highly resolved theoretical cross sections exhibit considerably more resonance structures than observed in the heavy-ion storage ring measurements at Heidelberg, Germany. Nonetheless, the detailed resonance complexes agree well with experiment, and the unfiied rates agree with the sum of experimentally derived DR and theoretical RR rates to ~ 20%, within experimental or theoretical uncertainties. Theoretical results may provide estimates of field ionization of rydberg levels close to the DR peak, and non-resonant background contirbutions, particularly close to the RR peak as E --> 0. More generally, the unified results avoid the physical and practical problems in astrophysical models inherent in the separation of electronic recombination into RR and DR on the one hand, and further subdivision into low-energy Delta n = 0 DR and high-energy Delta n > 0 DR in photoionized and collisionally ionized X-ray plasmas.

The Breit-Pauli R-matrix method developed under the Iron Project has been used to obtain extensive sets of oscillator strengths and transition probabilities for dipole allowed and intercombination fine structure transitions in carbon like ions, Ar XIII and Fe XXI. The complete set consists of 1274 fine structure bound energy levels and 198,259 oscillator strengths for Ar XIII, and 1611 bound levels and 300,079 oscillator strengths for Fe XXI. These correspond to levels of total angular momenta of 0 <= J <= 7 of even and odd parities formed from total spin of 2S+1=5,3,1, and orbital angular momenta 0 <= L <= 9 with n <= 10, 0 <= l <= 9 for each ion. The relativistic effects are included in the Breit-Pauli approximation. The close coupling wavefunction expansion for each ion is represented by the lowest 15 fine structure levels of target configuations, 2s2.2p, 2s.2p2 and 2p3. The energy levels are identified spectroscopically using a newly developed identification procedure. The procedure also makes a correspondence between the fine strucure energy levels and LS terms. This provides the check for completeness of the calculated levels. Comparison is made of the present energies and the f-values with the available observed and theoretical values. Present transition probabilitis agree very well with the relativistic atomic structure calculations of Mendoza et al. for the intercombination transitions, 2s.2p3(5So)2 - 2s2.2p2(3P){1,2},(1D)2. This further indicates that the importance of the neglected Breit interaction decreases with ion charge and constrains the uncertainty in the present calculations to within 15% even for the weak transitions.

The first complete set of unified cross sections and rate coefficients are calculated for photoionization and recombination of He- and Li-like ions using the relativistic Breit-Pauli R-matrix method. We present total, unified (e + ion) recombination rate coefficients for (e + C VI ---> C V) and (e + C V \longrightarrow C IV) including fine structure. Level-specific recombination rate coefficients up to the n = 10 levels are also obtained for the first time; these differ considerably from the approximate rates currently available. Applications to recombination-cascade coefficients in X-ray spectral models of K-alpha emission from the important He-like ions is pointed out. The overall uncertainty in the total recombination rates should not exceed 10-20%. Ionization fractions for Carbon are recomputed in the coronal approximation using the new rates. The present (e + ion) recombination rate coefficients are compared with several sets of available data, including previous LS coupling results, and `experimentally derived' rate coefficients. The role of relativistic fine structure, resolution of resonances, radiation damping, and interference effects is discussed. Two general features of recombination rates are noted: (i) the non-resonant (radiative recombination) peak as E,T ---> 0, and the (ii) the high-T resonant (di-electronic recombination) peak.

Theoretically self-consistent calculations for photoionization and (e + ion) recombination are described. An identical eigenfunction expansion for the ion is employed in coupled channel calculations for both processes, thus ensuring consistently accurate cross sections and rates in an ab initio manner. The theoretical treatment of (e + ion) recombination subsumes both the non-resonant recombination ('radiative recombination'), and the resonant recombination ('di-electronic recombination') processes in a unified scheme. In addition to the total, unified recombination rates, level-specific recombination rates and photoionization cross sections are obtained for a large number of atomic levels. Both relativistic Breit-Pauli, and non-relativistic LS coupling, calculations are carried out in the close coupling approximation using the R-matrix method. Although the calculations are computationally intensive, they yield nearly all photoionization and recombination parameters needed for astrophysical photoionization models with higher precision than hitherto possible, estimated at about 10-20 % from comparison with experimentally available data (including 'experimentally derived DR rates'). Results are electronically available for over 40 atoms and ions. Photoionization and recombination of He-, and Li-like C and Fe are described for X-ray modeling. The unified method yields total and complete (e+ion) recombination rate coeffcients, that can not otherwise be obtained theoretically or experimentally. "temp" 23 lines, 1524 characters

An extensive set of dipole-allowed, intercombination, and forbidden transition probabilities for Fe V is presented. The Breit-Pauli R-matrix (BPRM) method is used to calculate 1.46 x 10^6 oscillator strengths for the allowed and intercombination E1 transitions among 3,865 fine-structure levels dominated by configuration complexes with n <= 10 and l <= 9. These data are complemented by an atomic structure configuration interaction (CI) calculation using the SUPERSTRUCTURE program for 362 relativistic quadrupole (E2) and magnetic dipole (M1) transitions among 65 low-lying levels dominated by the 3d^4 and 3d^ 4s configurations. Procedures have been developed for the identification of the large number of fine-structure levels and transitions obtained through the BPRM calculations. The target ion Fe VI is represented by an eigenfunction expansion of 19 fine-structure levels of 3d^3 and a set of correlation configurations. Fe V bound levels are obtained with angular and spin symmetries SL\pi and J\pi of the (e + Fe VI) system such that 2S+1 = 5,3,1, L <= 10, J <= 8 of even and odd parities. The completeness of the calculated dataset is verified in terms of all possible bound levels belonging to relevant LS terms and transitions in correspondence with the LS terms. The fine-structure averaged relativistic values are compared with previous Opacity Project LS coupling data and other works. The 362 forbidden transition probabilities considerably extend the available data for the E2 and M1 transtions, and are in good agreement with those computed by Garstang for the 3d^4 transitions.

Ab initio theoretical calculations are reported for the electric (E1) dipole allowed and intercombination fine structure transitions in Fe V using the Breit-Pauli R-matrix (BPRM) method. We obtain 3865 bound fine structure levels of Fe V and $1.46 x 10^6$ oscillator strengths, Einstein A-coefficients and line strengths. In addition to the relativistic effects, the intermediate coupling calculations include extensive electron correlation effects that represent the complex configuration interaction (CI). Fe V bound levels are obtained with angular and spin symmetries $SL\pi$ and $J\pi$ of the (e + Fe VI) system such that $2S+1$ = 5,3,1, $L \leq$ 10, $J \leq 8$. The bound levels are obtained as solutions of the Breit-Pauli (e + ion) Hamiltonian for each $J\pi$, and are designated according to the `collision' channel quantum numbers. A major task has been the identification of these large number of bound fine structure levels in terms of standard spectroscopic designations. A new scheme, based on the analysis of quantum defects and channel wavefunctions, has been developed. The identification scheme aims particularly to determine the completeness of the results in terms of all possible bound levels for applications to analysis of experimental measurements and plasma modeling. An uncertainty of 10-20% for most transitions is estimated.