We present a detailed analysis of the gaseous component of the Si K edge using high-resolution Chandra spectra of low-mass X-ray binaries. We fit the spectra with amodified version of the ISMabs model, including new photoabsorption cross sectionscomputed for all Si ionic species. We estimate column densities for Si i, Si ii, Si iii, Si xii and Si xiii, which trace the warm, intermediate temperature and hot phases of the Galactic interstellar medium. We find that the ionic fractions of the first two phases are similar. This may be due to the physical state of the plasma determined by the temperature or to the presence of absorber material in the close vicinity of the sources. Our findings highlight the need for accurate modeling of the gaseous component before attempting to address the solid component.

Thanks to high-resolution and non-dispersive spectrometers onboard future X-ray missions such as XRISM and Athena, we are finally poised to answer important questions about the formation and evolution of galaxies and large-scale structure. However, we currently lack an adequate understanding of many atomic processes behind the spectral features we will soon observe. Large error bars on parameters as critical as transition energies and atomic cross sections can lead to unacceptable uncertainties in the calculations of e.g., elemental abundance, velocity, and temperature. Unless we address these issues, we risk limiting the full scientific potential of these missions. Laboratory astrophysics, which comprises theoretical and experimental studies of the underlying physics behind observable astrophysical processes, is therefore central to the success of these missions.

About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos.

We present a study of the C K-edge using high-resolution LETGS Chandra spectra of four novae during their super-soft-source (SSS) phase. We identified absorption lines due to C~ii K$\alpha$, C~iii K$\alpha$ and C~ iii K$\beta$ resonances. We used these astronomical observations to perform a benchmarking of the atomic data, which involves wavelength shifts of the resonances and photoionization cross-sections. We used improved atomic data to estimate the C~ii and C~iii column densities. The absence of physical shifts for the absorption lines, the consistence of the column densities between multiple observations and the high temperature required for the SSS nova atmosphere modeling support our conclusion about an ISM origin of the respective absorption lines. Assuming a collisional ionization equilibrium plasma the maximum temperature derived from the ratio of C~ii/C~iii column densities of the absorbers correspond to $T_{max}< 3.05\times10^{4}$ K.

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.

We present an X-ray absorption model for the interstellar medium, to be referred to as ISMabs, that takes into account both neutral and ionized species of cosmically abundant elements, and includes the most accurate atomic data available. Using high-resolution spectra from eight X-ray binaries obtained with the Chandra High Energy Transmission Grating Spectrometer, we proceed to benchmark the atomic data in the model particularly in the neon K-edge region. Compared with previous photoabsorption models, which solely rely on neutral species, the inclusion of ions leads to improved spectral fits. Fit parameters comprise the column densities of abundant contributors that allow direct estimates of the ionization states. ISMabs is provided in the appropriate format to be implemented in widely used X-ray spectral fitting packages such as XSPEC, ISIS and SHERPA.

Chandra high-resolution spectra toward eight low-mass Galactic binaries have been analyzed with a photoionization model that is capable of determining the physical state of the interstellar medium. Particular attention is given to the accuracy of the atomic data. Hydrogen column densities are derived with a broadband fit that takes into account pileup effects, and in general are in good agreement with previous results. The dominant features in the oxygen-edge region are O I and O II K$\alpha$ absorption lines whose simultaneous fits lead to average values of the ionization parameter of $\log\xi=-2.90$ and oxygen abundance of $A_{\rm O}=0.70$. The latter is relative to the standard by Grevesse & Sauval (1998), but a rescaling with the revision by Asplund et al. (2009) would lead to an average abundance value fairly close to solar. The low average oxygen column density ($N_{\rm O}=9.2 \times 10^{17}$ cm$^{-2}$) suggests a correlation with the low ionization parameters, the latter also being in evidence in the column density ratios OII/OI and OIII/OI that are estimated to be less than 0.1. We do not find conclusive evidence for absorption by any other compound but atomic oxygen.

An analytical formula is developed to represent accurately the photoabsorption cross section of O I for all energies of interest in X-ray spectral modeling. In the vicinity of the Kedge, a Rydberg series expression is used to fit R-matrix results, including important orbital relaxation effects, that accurately predict the absorption oscillator strengths below threshold and merge consistently and continuously to the above-threshold cross section. Further minor adjustments are made to the threshold energies in order to reliably align the atomic Rydberg resonances after consideration of both experimental and observed line positions. At energies far below or above the K-edge region, the formulation is based on both outer- and inner-shell direct photoionization, including significant shake-up and shake-off processes that result in photoionization-excitation and double photoionization contributions to the total cross section. The ultimate purpose for developing a definitive model for oxygen absorption is to resolve standing discrepancies between the astronomically observed and laboratory measured line positions, and between the inferred atomic and molecular oxygen abundances in the interstellar medium from XSTAR and SPEX spectral models.

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.

K-shell photoabsorption cross sections for the isonuclear C I - C IV ions have been computed using the R-matrix method. Above the K-shell threshold, the present results are in good agreement with the independent-particle results of Reilman & Manson (1979). Below threshold, we also compute the strong 1s -> np absorption resonances with the inclusion of important spectator Auger broadening effects. For the lowest 1s -> 2p, 3p resonances, comparisons to available C II, C III, and C IV experimental results show good agreement in general for the resonance strengths and positions, but unexplained discrepancies exist. Our results also provide detailed information on the C I K-shell photoabsorption cross section including the strong resonance features, since very limited laboratory experimental data exist. The resultant R-matrix cross sections are then used to model the Chandra X-ray absorption spectrum of the blazar Mkn 421.

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 present high-resolution spectroscopy of the neon K-shell and iron L-shell interstellar absorption edges in nine X-ray binaries using the High Energy Transmission Grating Spectrometer (HETGS) onboard the Chandra X-ray Observatory. We found that the iron absorption is well fit by an experimental determination of the cross-section for metallic iron, although with a slight wavelength shift of ~20 mA. The neon edge region is best fit by a model that includes the neutral neon edge and three Gaussian absorption lines. We identify these lines as due to the 1s-2p transitions from Ne II, Ne III, and Ne IX. As we found in our oxygen edge study, the theoretical predictions for neutral and low-ionization lines all require shifts of ~20 mA to match our data. Combined with our earlier oxygen edge study, we find that a best fit O/Ne ratio of 5.4+/-1.6, consistent with standard interstellar abundances. Our best fit Fe/Ne ratio of 0.20+/-0.03 is significantly lower than the interstellar value. We attribute this difference to iron depletion into dust grains in the interstellar medium. We make the first measurement of the neon ionization fraction in the ISM. We find Ne II/Ne I ~ 0.3 and Ne III/Ne I ~ 0.07. These values are larger than is expected given the measured ionization of interstellar helium. For Ne IX, our results confirm the detection of the hot ionized interstellar medium of the Galaxy.

Extensive calculations of the atomic data required for the spectral modelling of the K-shell photoabsorption of oxygen ions have been carried out in a multi-code approach. The present level energies and wavelengths for the highly ionized species (electron occupancies 2 <= N <= 4) are accurate to within 0.5 eV and 0.02 A, respectively. For N > 4, lack of measurements, wide experimental scatter, and discrepancies among theoretical values are handicaps in reliable accuracy assessments. The radiative and Auger rates are expected to be accurate to 10% and 20%, respectively, except for transitions involving strongly mixed levels. Radiative and Auger dampings have been taken into account in the calculation of photoabsorption cross sections in the K-threshold region, leading to overlapping lorentzian shaped resonances of constant widths that cause edge smearing. The behavior of the improved opacities in this region has been studied with the XSTAR modelling code using simple constant density slab models, and is displayed for a range of ionization parameters.

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.