We present HST/STIS observations and analysis of two prominent nebular structures around the central source of Eta Carinae, the knots C and D. The former is brighter than the latter for emission lines from intermediate or high ionization potential ions. The brightness of lines from intermediate and high ionization potential ions significantly decreases at phases around periastron. We do not see conspicuous changes in the brightness of lines from low ionization potential (<13.6 eV) that the total extinction towards the Weigelt structures is that the total extinction towards the Weigelt structures is AsubV =2/0. that the total extinction towards the Weigelt structures is AV = 2.0. Weigelt C and D are characterized by an electron density of that the total extinction towards the Weigelt structures is AV = 2.0. Weigelt C and D are characterized by an electron density of 10exp6.9 cm-3 that does not significantly change throughout the orbital cycle. The electron temperature varies from 5500 K (around periastron) to 7200 K (around apastron). The relative changes in the brightness of He I lines are well reproduced by the variations in the electron temperature alone. We found that, at phases around periastron, the electron temperature seems to be higher for Weigelt C than that of D. The Weigelt structures are located close to the Homunculus equatorial plane, at a distance of about 1240 AU from the central source. From the analysis of proper motion and age, the Weigelt complex can be associated with the equatorial structure called the Butterfly Nebula surrounding the central binary system.

Aims. In the context of black-hole accretion disks, we aim to compute the plasma-environment effects on the atomic parameters used to model the decay of K-vacancy states in moderately charged iron ions, namely Fe ix - Fe xvi. Methods. We used the fully relativistic multiconfiguration Dirac-Fock (MCDF) method approximating the plasma electron-nucleus and electron-electron screenings with a time-averaged Debye-Huckel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges $10^5$ - $10^7$ K and $10^{18}$ - $10^{22}$ cm$^{-3}$. Conclusions. This study confirms that the high-resolution X-ray spectrometers onboard the future XRISM and ATHENA space missions will be capable of detecting the lowering of the K edges of these ions due to the extreme plasma conditions occurring in accretion disks around compact objects.

Manganese (Mn) is a key Fe-group elements, commonly employed in stellar population and nucleosynthesis studies to explore the role of SN Ia. We have developed a new non-local thermodynamic equilibrium (NLTE) model of Mn, including new photo-ionisation cross-sections and new transition rates caused by collisions with H and H- atoms. We applied the model in combination with 1-dimensional (1D) LTE model atmospheres and 3D hydrodynamical simulations of stellar convection to quantify the impact of NLTE and convection on the line formation. We show that the effects of NLTE are present in Mn I and, to a lesser degree, in Mn II lines, and these increase with metallicity and with effective temperature of a model. Employing 3D NLTE radiative transfer, we derive new abundance of Mn in the Sun, A(Mn)=5.52 +/- 0.03 dex, consistent with the element abundance in C I meteorites. We also apply our methods to the analysis of three metal-poor benchmark stars. We find that 3D NLTE abundances are significantly higher than 1D LTE. For dwarfs, the differences between 1D NLTE and 3D NLTE abundances are typically within 0.15 dex, however, the effects are much larger in the atmospheres of giants owing to their more vigorous convection. We show that 3D NLTE successfully solves the ionisation and excitation balance for the RGB star HD 122563 that cannot be achieved by 1D LTE or 1D NLTE modelling. For HD 84937 and HD 140283, the ionisation balance is satisfied, however, the resonance Mn I triplet lines still show somewhat lower abundances compared to the high-excitation lines. Our results for the benchmark stars confirm that 1D LTE modelling leads to significant systematic biases in Mn abundances across the full wavelength range from the blue to the IR. We also produce a list of Mn lines that are not significantly biased by 3D and can be reliably, within the 0.1 dex uncertainty, modelled in 1D NLTE.

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.

Sunspot catalogs are very useful for studying the solar activity of the recent past. In this context, a catalog covering more than three solar cycles made by the astronomers of the Madrid Astronomical Observatory in Spain (nowadays, the National Astronomical Observatory) from 1952 until 1986 has been recovered. Moreover, a machine-readable version of this catalog has been made available. We have recovered abundant metadata and studied the reliability of this dataset by comparing it with other sunspot catalogs.

A prevailing open problem in planetary nebulae research, and photoionized gaseous nebulae research at large, is the systematic discrepancies in electron temperatures and ionic abundances as derived from recombination and collisionally excited lines. Peimbert (1967) proposed the presence of 'temperature fluctuations' in these nebulae, but the apparent amplitude of such fluctuations, as deduced from spectral diagnostics and/or abundance discrepancy factors, remain unexplained by standard photoionization modeling. While this and other alternative models to explain the temperature and abundance discrepancies remain inconclusive, recent observations seem to point at a connection between nebular abundance discrepancy factors and a binary nature of photoionizing stars. In this paper we show that large amplitude temperature fluctuations are expected to form in planetary nebulae photoionized by short-period binary stars. Resonant temperature fluctuations are first formed along the orbital disk around the binary stars, as the periodically varying ionizing radiation field induces periodic oscillations in the heating-minus-cooling function. Then, the temperatures fluctuations propagate vertically to the disk as thermal waves that later steepen into radiative shocks. The binary period of the ionizing stars is determinant in the formation and propagation of temperature fluctuations, as well as in associated density fluctuations. Fluctuations propagate efficiently only in systems with binary periods significantly shorter than the gas thermalization time, of the order of 10 days. Further, we propose temperature diagnostic line ratios that combine [O III] collisionally excited lines and O II recombination lines to determine the equilibrium temperature and the magnitude of resonant temperature fluctuations in nebulae.

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.

This is the final report of a three-paper series on the K-shell photoabsorption and photoionization of trace elements, namely F, Na, P, Cl, K, Sc, Ti, V, Cr, Mn, Co, Cu and Zn. K lines and edges from such elements are observed in the X-ray spectra of supernova remnants, galaxy clusters and accreting black holes and neutron stars, their diagnostic potential being limited by poor atomic data. We are completing the previously reported radiative datasets with new photoabsorption and photoionization cross sections for isoelectronic sequences with electron number $19\leq N\leq 26$. We are also giving attention to the access, integrity and usability of the whole resulting atomic database. Target representations are obtained with the atomic structure code AUTOSTRUCTURE. Where possible, cross sections for ground-configuration states are computed with the Breit--Pauli $R$-matrix method (BPRM) in either intermediate or $LS$ coupling including damping (radiative and Auger) effects; otherwise and more generally, they are generated in the isolated-resonance approximation with AUTOSTRUCTURE. Cross sections were computed with BPRM only for the K ($N=19$) and Ca ($N=20$) isoelectronic sequences, the latter in $LS$ coupling. For the rest of the sequences ($21\leq N \leq 26$), AUTOSTRUCTURE was run in $LS$-coupling mode taking into account damping effects. Comparisons between these two methods for K-like Zn XII and Ca-like Zn XI show that, to ensure reasonable accuracy, the $LS$ calculations must be performed taking into account the non-fine-structure relativistic corrections.

We have identified two new near-infrared emission lines in the spectra of planetary nebulae (PNe) arising from heavy elements produced by neutron capture reactions: [Te III] 2.1019 $\mu$m and [Br V] 1.6429 $\mu$m. [Te III] was detected in both NGC 7027 and IC 418, while [Br V] was seen in NGC 7027. The observations were obtained with the medium-resolution spectrograph EMIR on the 10.4m Gran Telescopio Canarias at La Palma, and with the high-resolution spectrograph IGRINS on the 2.7m Harlan J. Smith telescope at McDonald Observatory. New calculations of atomic data for these ions, specifically A-values and collision strengths, are presented and used to derive ionic abundances of Te$^{2+}$ and Br$^{4+}$. We also derive ionic abundances of other neutron-capture elements detected in the near-infrared spectra, and estimate total elemental abundances of Se, Br, Kr, Rb, and Te after correcting for unobserved ions. Comparison of our derived enrichments to theoretical predictions from AGB evolutionary models shows reasonable agreement for solar metallicity progenitor stars of $\sim$2 - 4 M$_{\odot}$. The spectrally-isolated [Br V] 1.6429 $\mu$m line has advantages for determining nebular Br abundances over optical [Br III] emission lines that can be blended with other features. Finally, measurements of Te are of special interest because this element lies beyond the first peak of the s-process, and thus provides new leverage on the abundance pattern of trans-iron species produced by AGB stars.

In most accreting black-hole systems the copious X-rays commonly observed from the inner-most regions are accompanied by a reflection spectrum. The latter is the signature of energetic photons reprocessed by the optically thick material of an accretion disk. Given their abundance and fluorescence yield, the iron K-shell lines are the most prominent features in the X-ray reflected spectrum. Their line profiles can be grossly broadened and skewed by Doppler effects and gravitational redshift. Consequently, modeling the reflection spectrum provides one of the best methods to measure, among other physical quantities, the black-hole spin. At present the accuracy of the spin estimates is called into question because the data fits require very high iron abundances: typically several times the solar value. Concurrently no plausible physical explanation has been proffered for these black-hole systems to be so iron rich. The most likely explanation for the supersolar iron abundances is model shortfall at very high densities ($>10^{18}$ cm$^{-3}$) due to atomic data shortcomings in this regime. We review the current observational evidence for the iron supersolar abundance in many black-hole systems, and show the effects of high density in state-of-the-art reflection models. We also briefly discuss our current efforts to produce new atomic data for high-density plasmas, which are required to refine the photoionization models.

Availability of reliable atomic and molecular opacity tables is essential in a wide variety of astronomical modeling: the solar and stellar interiors, stellar and planetary atmospheres, stellar evolution, pulsating stars, and protoplanetary disks, to name a few. With the advancement of powerful research techniques such as helio-seismology and asteroseismology, solar neutrino-flux measurements, exoplanet survey satellites, three-dimensional hydrodynamic atmospheric simulations (including non-LTE and granulation effects), high-performance computing of atomic and molecular data, and innovative plasma experiments the accuracy and completeness of opacity tables is being taken to an unprecedented level. The goal of the second Workshop on Astrophysical Opacities was to gather opacity data producers and consumers from both the atomic and molecular sectors to contribute to solving outstanding problems and to develop more effective and integrated interfaces. In this review we attempt to summa- rize the discussion at the workshop and propose future directions for opacity research.

A present prevailing open problem planetary nebulae research, and photoionized gaseous nebulae research at large, is the systematic discrepancies in ionic abundances derived from recombination and collisionally excited lines in many H II regions and planetary nebulae. Peimbert (1967) proposed that these discrepancies were due to 'temperature fluctuations' in the plasma, but the amplitude of such fluctuations remain unexplained by standard phtoionization modeling. In this letter we show that large amplitude temperature oscillations are expected to form in gaseous nebulae photoionized by short-period binary stars. Such stars yield periodically varying ionizing radiation fields, which induce periodic oscilla- tions in the heating-minus-cooling function across the nebula. For flux oscillation periods of a few days any temperature perturbations in the gas with frequencies similar to those of the ionizing source will undergo resonant amplification. In this case, the rate of growth of the perturbations increases with the amplitude of the variations of the ionizing flux and with decreasing nebular equilibrium temperature. We also present a line ratios diagnostic plot that combines [O III] collisional lines and O II recombination lines for diagnosing equilibrium and fluctuation am- plitude temperatures in gaseous nebulae. When applying this diagnostic to the planetary nebula M 1-42 we find an equilibrium temperature of ~6000 K and a resonant temperature fluctuation amplitude (Trtf ) of ~4000 K. This equilibrium temperature is significantly lower than the temperature estimated when temperature perturbations are ignored.

We are concerned with improving the diagnostic potential of the K lines and edges of elements with low cosmic abundances that are observed in the X-ray spectra of supernova remnants, galaxy clusters and accreting black holes and neutron stars. Since accurate photoabsorption and photoionization cross sections are needed in their spectral models, they have been computed for isoelectronic sequences with electron number $12\leq N\leq 18$ using a multi-channel method. Target representations are obtained with the atomic structure code AUTOSTRUCTURE, and ground-state cross sections are computed with the Breit--Pauli $R$-matrix method (BPRM) in intermediate coupling, including damping (radiative and Auger) effects. The contributions from channels associated with the 2s-hole $[2{\rm s}]\mu$ target configurations and those containing 3d orbitals are studied in the Mg and Ar isoelectronic sequences. Cross sections for the latter ions are also calculated in the isolated-resonance approximation as implemented in AUTOSTRUCTURE and compared with BPRM to test their accuracy. It is confirmed that the collisional channels associated with the $[2{\rm s}]\mu$ target configurations must be taken into account owing to significant increases in the monotonic background cross section between the L and K edges. Target configurations with 3d orbitals give rise to fairly conspicuous unresolved transition arrays in the L-edge region, but to a much lesser extent in the K-edge which is our main concern; therefore, they have been neglected throughout owing to their computationally intractable channel inventory, thus allowing the computation of cross sections for all the ions with $12\leq N\leq 18$ in intermediate coupling with BPRM. We find that the isolated-resonance approximations performs satisfactorily.

Our ability to model the shapes and strengths of iron lines in the solar spectrum is a critical test of the accuracy of the solar iron abundance, which sets the absolute zero-point of all stellar metallicities. We use an extensive 463-level Fe atom with new photoionisation cross-sections for FeI as well as quantum mechanical calculations of collisional excitation and charge transfer with neutral hydrogen; the latter effectively remove a free parameter that has hampered all previous line formation studies of Fe in non-local thermodynamic equilibrium (NLTE). For the first time, we use realistic 3D NLTE calculations of Fe for a quantitative comparison to solar observations. We confront our theoretical line profiles with observations taken at different viewing angles across the solar disk with the Swedish 1-m Solar Telescope. We find that 3D modelling well reproduces the observed centre-to-limb behaviour of spectral lines overall, but highlight aspects that may require further work, especially cross-sections for inelastic collisions with electrons. Our inferred solar iron abundance is log(eps(Fe))=7.48+-0.04.

After an extensive assessment of the effective collision strengths available to model the He-like triplet of C V, N VI, O VII, Ne IX, Mg XI and Si XIII in collisionally dominated plasmas, new accurate effective collision strengths are reported for Ne IX. The uncertainty intervals of the density and temperature diagnostics due to the atomic data errors are also determined for both Maxwell-Boltzmann and $\kappa$ electron-energy distributions. It is shown that these uncertainty bands limit the temperature range where the temperature line-ratio diagnostic can be applied and its effectiveness to discern the electron-energy distribution type. These findings are benchmarked with Chandra and XMM-Newton spectra of stellar coronae and with tokamak measurements.

We identify [Rb IV] 1.5973 and [Cd IV] 1.7204 micron emission lines in high-resolution (R=40,000) near-infrared spectra of the planetary nebulae (PNe) NGC 7027 and IC 5117, obtained with the IGRINS spectrometer on the 2.7-m telescope at McDonald Observatory. We also identify [Ge VI] 2.1930 $\mu$m in NGC 7027. Alternate identifications for these features are ruled out based on the absence of other multiplet members and/or transitions with the same upper levels. Ge, Rb, and Cd can be enriched in PNe by s-process nucleosynthesis during the asymptotic giant branch (AGB) stage of evolution. To determine ionic abundances, we calculate [Rb IV] collision strengths and use approximations for those of [Cd IV] and [Ge VI]. Our identification of [Rb IV] 1.5973 $\mu$m is supported by the agreement between Rb3+/H+ abundances found from this line and the 5759.55 A feature in NGC 7027. Elemental Rb, Cd, and Ge abundances are derived with ionization corrections based on similarities in ionization potential ranges between the detected ions and O and Ne ionization states. Our analysis indicates abundances 2-4 times solar for Rb and Cd in both nebulae. Ge is subsolar in NGC 7027, but its abundance is uncertain due to the large and uncertain ionization correction. The general consistency of the measured relative s-process enrichments with predictions from models appropriate for these PNe (2.0-2.5 M_sun, [Fe/H]= -0.37) demonstrates the potential of using PN compositions to test s-process nucleosynthesis models.

We conduct a spectroscopic search of quasars observed by the Sloan Digital Sky Survey (SDSS) with broad absorption line (BAL) troughs due to Mg II and troughs due to Fe II that simultaneously exhibit strong Balmer narrow emission lines (NELs). We find that in a redshift range of 0.4 less than or equal to z less than or equal to 0.9 approximately 23 of the 70 Mg II BALs and 4 of a subset of 15 Fe II BALs exhibit strong Balmer emission. We also find significant fractions of Mg II BALs (approximately 23%) and those Mg II BALs with Fe II troughs (approximately 27%) have strong continuum reddening, E(B - V) greater than or equal to 0.1. From measurements of the Balmer decrement in three objects, we find similarly significant reddening of the NEL region in three of the four objects; the NELs in the fourth object are not measurable. We also include one object in this study not taken from the SDSS sample that shows Fe II absorption and strong narrow emission, but due to measurement uncertainty and low continuum reddening the comparison is consistent but inconclusive. We find a trend in both the Mg II and Fe II BAL samples between the NEL reddening and continuum reddening. Because the narrow line reddening is consistent with the continuum reddening in every object in the two SDSS samples, it suggests that the reddening sources in these objects likely exist at larger radial distances than the narrow line regions from the central nucleus.

We present extensive calculations of radiative transition rates and electron impact collision strengths for Fe II. The data sets involve 52 levels from the $3d\,^7$, $3d\,^64s$, and $3d\,^54s^2$ configurations. Computations of $A$-values are carried out with a combination of state-of-the-art multiconfiguration approaches, namely the relativistic Hartree--Fock, Thomas--Fermi--Dirac potential, and Dirac--Fock methods; while the $R$-matrix plus intermediate coupling frame transformation, Breit--Pauli $R$-matrix and Dirac $R$-matrix packages are used to obtain collision strengths. We examine the advantages and shortcomings of each of these methods, and estimate rate uncertainties from the resulting data dispersion. We proceed to construct excitation balance spectral models, and compare the predictions from each data set with observed spectra from various astronomical objects. We are thus able to establish benchmarks in the spectral modeling of [Fe II] emission in the IR and optical regions as well as in the UV Fe II absorption spectra. Finally, we provide diagnostic line ratios and line emissivities for emission spectroscopy as well as column densities for absorption spectroscopy. All atomic data and models are available online and through the AtomPy atomic data curation environment.

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.

We report on spectra of two positions in the XA region of the Cygnus Loop supernova remnant obtained with the InfraRed Spectrograph on the Spitzer Space Telescope. The spectra span the 10-35 micron wavelength range, which contains a number of collisionally excited forbidden lines. These data are supplemented by optical spectra obtained at the Whipple Observatory and an archival UV spectrum from the International Ultraviolet Explorer. Coverage from the UV through the IR provides tests of shock wave models and tight constraints on model parameters. Only lines from high ionization species are detected in the spectrum of a filament on the edge of the remnant. The filament traces a 180 km/s shock that has just begun to cool, and the oxygen to neon abundance ratio lies in the normal range found for Galactic H II regions. Lines from both high and low ionization species are detected in the spectrum of the cusp of a shock-cloud interaction, which lies within the remnant boundary. The spectrum of the cusp region is matched by a shock of about 150 km/s that has cooled and begun to recombine. The post-shock region has a swept-up column density of about 1.3E18 cm^-2, and the gas has reached a temperature of 7000 to 8000 K. The spectrum of the Cusp indicates that roughly half of the refractory silicon and iron atoms have been liberated from the grains. Dust emission is not detected at either position.

The classic optical nebular diagnostics [N II], [O II], [O III], [S II], [S III], and [Ar III] are employed to search for evidence of non-Maxwellian electron distributions, namely $\kappa$ distributions, in a sample of well-observed Galactic H II regions. By computing new effective collision strengths for all these systems and A-values when necessary (e.g. S II), and by comparing with previous collisional and radiative datasets, we have been able to obtain realistic estimates of the electron-temperature dispersion caused by the atomic data, which in most cases are not larger than $\sim 10$%. If the uncertainties due to both observation and atomic data are then taken into account, it is plausible to determine for some nebulae a representative average temperature while in others there are at least two plasma excitation regions. For the latter, it is found that the diagnostic temperature differences in the high-excitation region, e.g. $T_e$(O III), $T_e$(S III), and $T_e$(Ar III), cannot be conciliated by invoking $\kappa$ distributions. For the low excitation region, it is possible in some, but not all, cases to arrive at a common, lower temperature for [N II], [O II], and [S II] with $\kappa\approx 10$, which would then lead to significant abundance enhancements for these ions. An analytic formula is proposed to generate accurate $\kappa$-averaged excitation rate coefficients (better than 10% for $\kappa \geq 5$) from temperature tabulations of the Maxwell-Boltzmann effective collision strengths.

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 study the problem of time-dependent photoionization of low density gaseous nebulae subjected to sudden changes in the intensity of ionizing radiation. To this end, we write a computer code that solves the full time-dependent energy balance, ionization balance, and radiation transfer equations in a self-consistent fashion for a simplified pure hydrogen case. It is shown that changes in the ionizing radiation yield ionization/thermal fronts that propagate through the cloud, but the propagation times and response times to such fronts vary widely and non-linearly from the illuminated face of the cloud to the ionization front (IF). Ionization/thermal fronts are often supersonic, and in slabs initially in pressure equilibrium such fronts yield large pressure imbalances that are likely to produce important dynamical effects in the cloud. Further, we studied the case of periodic variations in the ionizing flux. It is found that the physical conditions of the plasma have complex behaviors that differ from any steady-state solutions. Moreover, even the time average ionization and temperature is different from any steady-state case. This time average is characterized by over-ionization and a broader IF with respect to the steady-state solution for a mean value of the radiation flux. Around the time average of physical conditions there is large dispersion in instantaneous conditions, particularly across the IF, which increases with the period of radiation flux variations. Moreover, the variations in physical conditions are asynchronous along the slab due to the combination of non-linear propagation times for thermal/ionization fronts and equilibration times.

We present detailed analyses of oxygen K absorption in the interstellar medium (ISM) using four high-resolution Chandra spectra towards the X-ray low-mass binary XTE J1817-330. The 11-25 A broadband is described with a simple absorption model that takes into account the pileup effect and results in an estimate of the hydrogen column density. The oxygen K-edge region (21-25 A) is fitted with the physical warmabs model, which is based on a photoionization model grid generated with the xstar code with the most up-to-date atomic database. This approach allows a benchmark of the atomic data which involves wavelength shifts of both the K lines and photoionization cross sections in order to fit the observed spectra accurately. As a result we obtain: a column density of N(H)=1.38+/-0.01\times 10^21 cm^-2; ionization parameter of log(xi)=-2.70+/-0.023; oxygen abundance of A(O)= 0.689^{+0.015}_{-0.010}; and ionization fractions of OI/O = 0.911, OII/O = 0.077, and OIII/O = 0.012 that are in good agreement with previous studies. Since the oxygen abundance in warmabs is given relative to the solar standard of Grevesse et al. (1998), a rescaling with the revision by Asplund et al. (2009) yields A(O)=0.952^{+0.020}_{-0.013}, a value close to solar that reinforces the new standard. We identify several atomic absorption lines Kalpha, Kbeta, and Kgamma in OI and OII; and Kalpha in OIII, OVI, and OVII - last two probably residing in the neighborhood of the source rather than in the ISM. This is the first firm detection of oxygen K resonances with principal quantum numbers n>2 associated to ISM cold absorption.

We present a method for computing uncertainties in spectral models, i.e. level populations, line emissivities, and emission line ratios, based upon the propagation of uncertainties originating from atomic data. We provide analytic expressions, in the form of linear sets of algebraic equations, for the coupled uncertainties among all levels. These equations can be solved efficiently for any set of physical conditions and uncertainties in the atomic data. We illustrate our method applied to spectral models of O III and Fe II and discuss the impact of the uncertainties on atomic systems under different physical conditions. As to intrinsic uncertainties in theoretical atomic data, we propose that these uncertainties can be estimated from the dispersion in the results from various independent calculations. This technique provides excellent results for the uncertainties in A-values of forbidden transitions in [Fe II].

We investigate statistical equilibrium of neutral and singly-ionized strontium in late-type stellar atmospheres. Particular attention is given to the completeness of the model atom, which includes new energy levels, transition probabilities, photoionization and electron-impact excitation cross-sections computed with the R-matrix method. The NLTE model is applied to the analysis of Sr I and Sr II lines in the spectra of the Sun, Procyon, Arcturus, and HD 122563, showing a significant improvement in the ionization balance compared to LTE line formation calculations, which predict abundance discrepancies of up to 0.5 dex. The solar Sr abundance is log A = 2.93 \pm 0.04 dex, in agreement with the meteorites. A grid of NLTE abundance corrections for Sr I and Sr II lines covering a large range of stellar parameters is presented.

We present a study of Broad Absorption Line (BAL) quasar outflows that show S IV ?1063 and S IV* ?1073 troughs. The fractional abundance of S IV and C IV peak at similar value of the ionization parameter, implying that they arise from the same physical component of the outflow. Detection of the S IV* troughs will allow us to determine the distance to this gas with higher resolution and higher signal-to-noise spectra, therefore providing the distance and energetics of the ubiquitous C IV BAL outflows. In our bright sample of 156 SDSS quasars 14% show C IV and 1.9% S IV troughs, which is consistent with a fainter magnitude sample with twice as many objects. One object in the fainter sample shows evidence of a broad S IV trough without any significant trough present from the excited state line, which implies that this outflow could be at a distance of several kpc. Given the fractions of C IV and S IV, we establish firm limits on the global covering factor on S IV that ranges from 2.8% to 21% (allowing for the k-correction). Comparison of the expected optical depth for these ions with their detected percentage suggests that these species arise from common outflows with a covering factor closer to the latter.

This workshop brought together scientists (including atomic physicists, theoretical astrophysicists and astronomers) concerned with the completeness and accuracy of atomic data for astrophysical applications. The topics covered in the workshop included the evaluation of uncertainties in atomic data, the propagation of such uncertainties in chemical abundances, and the feedback between observations and calculations. On a different level, we also discussed communication issues such as how to ensure that atomic data are correctly understood and used, and which forum is the best one for a fluid interaction between all communities involved in the production and use of atomic data. This paper reports on the discussions held during the workshop and introduces AstroAtom, a blog created as a platform for timely and open discussions on the needs and concerns over atomic data, and their effects on astronomical research. The complete proceedings will be published on http://astroatom.wordpress.com/.

We present R-matrix calculations of photoabsorption and photoionization cross sections across the K edge of the Li-like to Ca-like ions stages of Ni. Level-resolved, Breit-Pauli calculations were performed for the Li-like to Na-like stages. Term-resolved calculations, which include the mass-velocity and Darwin relativistic corrections, were performed for the Mg-like to Ca-like ion stages. This data set is extended up to Fe-like Ni using the distorted wave approximation as implemented by AUTOSTRUCTURE. The R-matrix calculations include the effects of radiative and Auger dampings by means of an optical potential. The damping processes affect the absorption resonances converging to the K thresholds causing them to display symmetric profiles of constant width that smear the otherwise sharp edge at the K-shell photoionization threshold. These data are important for the modeling of features found in photoionized plasmas.

We present absorption line analysis of the outflow in the quasar IRAS F04250-5718. Far-ultraviolet data from the Cosmic Origins Spectrograph onboard the Hubble Space Telescope reveal intrinsic narrow absorption lines from high ionization ions (e.g., C IV, N V, and O VI) as well as low ionization ions (e.g., C II and Si III). We identify three kinematic components with central velocities ranging from ~-50 to ~-230 km/s. Velocity dependent, non-black saturation is evident from the line profiles of the high ionization ions. From the non-detection of absorption from a metastable level of C II, we are able to determine that the electron number density in the main component of the outflow is < 30 per cubic cm. Photoionization analysis yields an ionization parameter log U ~ -1.6 +/- 0.2, which accounts for changes in the metallicity of the outflow and the shape of the incident spectrum. We also consider solutions with two ionization parameters. If the ionization structure of the outflow is due to photoionization by the active galactic nucleus, we determine that the distance to this component from the central source is > 3 kpc. Due to the large distance determined for the main kinematic component, we discuss the possibility that this outflow is part of a galactic wind.

We investigate the absorption structure of the oxygen in the interstellar medium by analyzing {\it XMM}-Newton observations of the low mass X-ray binary Sco X-1. We use simple models based on the O {\sc i} atomic cross section from different sources to fit the data and evaluate the impact of the atomic data in the interpretation of astrophysical observations. We show that relatively small differences in the atomic calculations can yield spurious results. We also show that the most complete and accurate set of atomic cross sections successfully reproduce the observed data in the $21-24.5$ {\AA} wavelength region of the spectrum. Our fits indicate that the absorption is mainly due to neutral gas with an ionization parameter of $\xi=10^{-4}$ erg cm s$^{-1}$, and an oxygen column density of $N_{\mathrm{O}}\approx 8-10\times 10^{17}$ cm$^{-2}$. Our models are able to reproduce both the K edge and the K$\alpha$ absorption line from O {\sc i}, which are the two main features in this region. We find no conclusive evidence for absorption by other than atomic oxygen.

We present $R$-matrix calculations of photoabsorption and photoionization cross sections across the K-edge of Mg, Si, S, Ar, and Ca ions with more than 10 electrons. The calculations include the effects of radiative and Auger damping by means of an optical potential. The wave functions are constructed from single-electron orbital bases obtained using a Thomas--Fermi--Dirac statistical model potential. Configuration interaction is considered among all states up to $n=3$. The damping processes affect the resonances converging to the K-thresholds causing them to display symmetric profiles of constant width that smear the otherwise sharp edge at the photoionization threshold. These data are important for modeling of features found in photoionized plasmas.

We present new atomic data (radiative transitions rates and collision strengths) from large scale calculations and a non-LTE spectral model for Fe III. This model is in very good agreement with observed astronomical emission spectra, in contrast with previous models that yield large discrepancies with observations. The present atomic computations employ a combination of atomic physics methods, e.g. relativistic Hatree-Fock, the Thomas-Fermi-Dirac potential, and Dirac-Fock computation of A-values and R-matrix with intermediate coupling frame transformation and Dirac R-matrix. We study the advantages and shortcomings of each method. It is found that the Dirac R-matrix collision strengths yield excellent agreement with observations, much improved over previously available models. By contrast, the transformation of LS-coupling R-matrix fails to yield accurate effective collision strengths at around 10^4 K, despite using very large configuration expansions, due to the limited treatment of spin-orbit effects in the near threshold resonances of the collision strengths. The present work demonstrates that accurate atomic data for low ionization iron-peak species is now within reach.

We report the discovery of very high energy gamma-ray emission from the direction of the SNR G54.1+0.3 using the VERITAS ground-based gamma-ray observatory. The TeV signal has an overall significance of 6.8$\sigma$ and appears point-like given the 5$^{arcminute}$ resolution of the instrument. The integral flux above 1 TeV is 2.5% of the Crab Nebula flux and significant emission is measured between 250 GeV and 4 TeV, well described by a power-law energy spectrum dN/dE $\sim$ E$^{-\Gamma}$ with a photon index $\Gamma= 2.39\pm0.23_{stat}\pm0.30_{sys}$. We find no evidence of time variability among observations spanning almost two years. Based on the location, the morphology, the measured spectrum, the lack of variability and a comparison with similar systems previously detected in the TeV band, the most likely counterpart of this new VHE gamma-ray source is the PWN in the SNR G54.1+0.3. The measured X-ray to VHE gamma-ray luminosity ratio is the lowest among all the nebulae supposedly driven by young rotation-powered pulsars, which could indicate a particle-dominated PWN.

We show that time variations in the UV ionizing continuum of quasars, on scales of $\sim$1 year, affect the dynamic structure of the plasmas responsible for low ionization broad absorption lines. Variations of the ionizing continuum produce non-equilibrium photoionization conditions over a significant fraction of the absorbing clouds and supersonically moving ionization fronts. When the flux drops the contraction of the ionized region drives a supersonic cooling front towards the radiation source and a rarefaction wave in the opposite direction. The pressure imbalance is compensated by an increased speed of the cool gas relative to the front. When the flux recovers the cool gas is re-ionized and re-heated by a supersonic ionization front traveling away from the radiation source and a forward shock is created. The reheated clouds equilibrate to a temperature of $\sim 10^4$ K and are observed to have different radial velocities than the main cloud. Such fragmentation seems consistent with the multicomponent structure of troughs seen in some objects. The velocity differences measured among various components in the quasars QSO 2359--1241 and SDSS J0318--0600 can be reproduced by our model if strong magnetic fields ($\sim$10 mG) are present within the clouds.

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

Using high resolution VLT spectra, we study the multi-component outflow systems of two quasars exhibiting intrinsic Fe II absorption (QSO 2359-1241 and SDSS J0318-0600). From the extracted ionic column densities and using photoionization modeling we determine the gas density, total column density, and ionization parameter for several of the components. For each object the largest column density component is also the densest, and all other components have densities of roughly 1/4 of that of the main component. We demonstrate that all the absorbers lie roughly at the same distance from the source. Further, we calculate the total kinetic luminosities and mass outflow rates of all components and show that these quantities are dominated by the main absorption component.

We report on observations of very high-energy gamma rays from the shell-type supernova remnant Cassiopeia A with the VERITAS stereoscopic array of four imaging atmospheric Cherenkov telescopes in Arizona. The total exposure time for these observations is 22 hours, accumulated between September and November of 2007. The gamma-ray source associated with the SNR Cassiopeia A was detected above 200 GeV with a statistical significance of 8.3 s.d. The estimated integral flux for this gamma-ray source is about 3% of the Crab-Nebula flux. The photon spectrum is compatible with a power law dN/dE ~ E^(-Gamma) with an index Gamma = 2.61 +/- 0.24(stat) +/- 0.2(sys). The data are consistent with a point-like source. We provide a detailed description of the analysis results, and discuss physical mechanisms that may be responsible for the observed gamma-ray emission.

We report the detection of very-high-energy (VHE) gamma-ray emission from supernova remnant (SNR) G106.3+2.7. Observations performed in 2008 with the VERITAS atmospheric Cherenkov gamma-ray telescope resolve extended emission overlapping the elongated radio SNR. The 7.3 sigma (pre-trials) detection has a full angular extent of roughly 0.6deg by 0.4deg. Most notably, the centroid of the VHE emission is centered near the peak of the coincident 12CO (J = 1-0) emission, 0.4deg away from the pulsar PSR J2229+6114, situated at the northern end of the SNR. Evidently the current-epoch particles from the pulsar wind nebula are not participating in the gamma-ray production. The VHE energy spectrum measured with VERITAS is well characterized by a power law dN/dE = N_0(E/3 TeV)^{-G} with a differential index of G = 2.29 +/- 0.33stat +/- 0.30sys and a flux of N_0 = (1.15 +/- 0.27stat +/- 0.35sys)x 10^{-13} cm^{-2} s^{-1} TeV^{-1}. The integral flux above 1 TeV corresponds to ~5 percent of the steady Crab Nebula emission above the same energy. We describe the observations and analysis of the object and briefly discuss the implications of the detection in a multiwavelength context.

We present high spectral resolution VLT observations of the BAL quasar SDSS J0318-0600. This high quality data set allows us to extract accurate ionic column densities and determine an electron number density of n_e=10^3.3 +/- 0.2 cm^-3 for the main outflow absorption component. The heavily reddened spectrum of SDSS J0318-0600 requires purely silicate dust with a reddening curve characteristic of predominately large grains, from which we estimate the bolometric luminosity. We carry out photoionization modeling to determine the total column density, ionization parameter and distance of the gas and find that the photionization models suggest abundances greater than solar. Due to the uncertainty in the location of the dust extinction, we arrive at two viable distances for the main ouflow component from the central source, 6 and 18 kpc, where we consider the 6 kpc location as somewhat more physically plausable. Assuming the canonical global covering of 20% for the outflow and a distance of 6 kpc, our analysis yields a mass flux of 120 M_sun yr^-1 and a kinetic luminosity that is ~0.1% of the bolometric luminosity of the object. Should the dust be part of the outflow, then these values are ~4x larger. The large mass flux and kinetic luminosity make this outflow a significant contributor to AGN feedback processes.

We present a detailed analysis of the Astrophysical Research Consortium 3.5 m telescope spectrum of QSO SDSS J0838+2955. The object shows three broad absorption line (BAL) systems at 22,000, 13,000, and 4900 km s^-1 blueshifted from the systemic redshift of z=2.043. Of particular interest is the lowest velocity system that displays absorption from low-ionization species such as Mg II, Al II, Si II, Si II*, Fe II and Fe II*. Accurate column densities were measured for all transitions in this lowest velocity BAL using an inhomogeneous absorber model. The ratio of column densities of Si II* and Fe II* with respect to their ground states gave an electron number density of log n_e (cm^-3) = 3.75 +/- 0.22 for the outflow. Photoionization modeling with careful regards to chemical abundances and the incident spectral energy distribution predicts an ionization parameter of log U_H = -1.93 +/- 0.21 and a hydrogen column density of log N_H (cm^-2) = 20.80 +/- 0.28. This places the outflow at 3.3+1.5-1.0 kpc from the central AGN. Assuming that the fraction of solid angle subtended by the outflow is 0.2, these values yield a kinetic luminosity of (4.5+3.1-1.8) x 10^45 erg s^-1, which is (1.4+1.1-0.6)% the bolometric luminosity of the QSO itself. Such large kinetic luminosity suggests that QSO outflows are a major contributor to AGN feedback mechanisms.

Although Galactic cosmic rays (protons and nuclei) are widely believed to be dominantly accelerated by the winds and supernovae of massive stars, definitive evidence of this origin remains elusive nearly a century after their discovery [1]. The active regions of starburst galaxies have exceptionally high rates of star formation, and their large size, more than 50 times the diameter of similar Galactic regions, uniquely enables reliable calorimetric measurements of their potentially high cosmic-ray density [2]. The cosmic rays produced in the formation, life, and death of their massive stars are expected to eventually produce diffuse gamma-ray emission via their interactions with interstellar gas and radiation. M 82, the prototype small starburst galaxy, is predicted to be the brightest starburst galaxy in gamma rays [3, 4]. Here we report the detection of >700 GeV gamma rays from M 82. From these data we determine a cosmic-ray density of 250 eV cm-3 in the starburst core of M 82, or about 500 times the average Galactic density. This result strongly supports that cosmic-ray acceleration is tied to star formation activity, and that supernovae and massive-star winds are the dominant accelerators.

The recent detection by the Fermi gamma-ray space telescope of high-energy gamma-rays from the radio galaxy NGC 1275 makes the observation of the very high energy (VHE: E > 100 GeV) part of its broadband spectrum particularly interesting, especially for the understanding of active galactic nuclei (AGN) with misaligned multi-structured jets. The radio galaxy NGC 1275 was recently observed by VERITAS at energies above 100 GeV for about 8 hours. No VHE gamma-ray emission was detected by VERITAS from NGC 1275. A 99% confidence level upper limit of 2.1% of the Crab Nebula flux level is obtained at the decorrelation energy of approximately 340 GeV, corresponding to 19% of the power-law extrapolation of the Fermi Large Area Telescope (LAT) result.

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.

Reliable atomic data have been computed for the spectral modeling of the nitrogen K lines, which may lead to useful astrophysical diagnostics. Data sets comprise valence and K-vacancy level energies, wavelengths, Einstein $A$-coefficients, radiative and Auger widths and K-edge photoionization cross sections. An important issue is the lack of measurements which are usually employed to fine-tune calculations so as to attain spectroscopic accuracy. In order to estimate data quality, several atomic structure codes are used and extensive comparisons with previous theoretical data have been carried out. In the calculation of K photoabsorption with the Breit--Pauli $R$-matrix method, both radiation and Auger damping, which cause the smearing of the K edge, are taken into account. This work is part of a wider project to compute atomic data in the X-ray regime to be included in the database of the popular {\sc xstar} modeling code.

The spectrum from the black hole X-ray transient GRO J1655-40. obtained using the $Chandra$ High Energy Transmission Grating (HETG) in 2005 is notable as a laboratory for the study of warm absorbers, and for the presence of many lines from odd-$Z$ elements between Na and Co (and Ti and Cr) not previously observed in X-rays. We present synthetic spectral models which can be used to constrain these element abundances and other parameters describing the outflow from the warm absorber in this object. We present results of fitting to the spectrum using various tools and techniques, including automated line fitting, phenomenological models, and photoionization modeling. We show that the behavior of the curves of growth of lines from H-like and Li-like ions indicate that the lines are either saturated or affected by filling-in from scattered or a partially covered continuum source. We confirm the conclusion of previous work by \cite{Mill06} and \cite{Mill08} which shows that the ionization conditions are not consistent with wind driving due to thermal expansion. The spectrum provides the opportunity to measure abundances for several elements not typically observable in the X-ray band. These show a pattern of enhancement for iron peak elements, and solar or sub-solar values for elements lighter than calcium. Models show that this is consistent with enrichment by a core-collapse supernova. We discuss the implications these values for the evolutionary history of this system.

We present evidence that the very-high-energy (VHE, E > 100 GeV) gamma-ray emission coincident with the supernova remnant IC 443 is extended. IC 443 contains one of the best-studied sites of supernova remnant/molecular cloud interaction and the pulsar wind nebula CXOU J061705.3+222127, both of which are important targets for VHE observations. VERITAS observed IC 443 for 37.9 hours during 2007 and detected emission above 300 GeV with an excess of 247 events, resulting in a significance of 8.3 standard deviations (sigma) before trials and 7.5 sigma after trials in a point-source search. The emission is centered at 06 16 51 +22 30 11 (J2000) +- 0.03_stat +- 0.08_sys degrees, with an intrinsic extension of 0.16 +- 0.03_stat +- 0.04_sys degrees. The VHE spectrum is well fit by a power law (dN/dE = N_0 * (E/TeV)^-Gamma) with a photon index of 2.99 +- 0.38_stat +- 0.3_sys and an integral flux above 300 GeV of (4.63 +- 0.90_stat +- 0.93_sys) * 10^-12 cm^-2 s^-1. These results are discussed in the context of existing models for gamma-ray production in IC 443.

We present results from a long-term monitoring campaign on the TeV binary LSI +61 303 with VERITAS at energies above 500 GeV, and in the 2-10 keV hard X-ray bands with RXTE and Swift, sampling nine 26.5 day orbital cycles between September 2006 and February 2008. The binary was observed by VERITAS to be variable, with all integrated observations resulting in a detection at the 8.8 sigma (2006/2007) and 7.3 sigma (2007/2008) significance level for emission above 500 GeV. The source was detected during active periods with flux values ranging from 5 to 20% of the Crab Nebula, varying over the course of a single orbital cycle. Additionally, the observations conducted in the 2007-2008 observing season show marginal evidence (at the 3.6 sigma significance level) for TeV emission outside of the apastron passage of the compact object around the Be star. Contemporaneous hard X-ray observations with RXTE and Swift show large variability with flux values typically varying between 0.5 and 3.0*10^-11 ergs cm^-2 s^-1 over a single orbital cycle. The contemporaneous X-ray and TeV data are examined and it is shown that the TeV sampling is not dense enough to detect a correlation between the two bands.

We present results of deep echelle spectrophotometry of the brightest knot of the HH202 in the Orion Nebula --HH202-S-- using the ultraviolet Visual Echelle Spectrograph (UVES). The high spectral resolution has permitted to separate the component associated with the ambient gas from that associated with the gas flow. We derive electron densities and temperatures for both components, as well as the chemical abundances of several ions and elements from collisionally excited lines, including the first determinations of Ca^{+} and Cr^{+} abundances in the Orion Nebula. We also calculate the He^{+}, C^{2+}, O^{+} and O^{2+} abundances from recombination lines. The difference between the O^{2+} abundances determined from collisionally excited and recombination lines --the so-called abundance discrepancy factor-- is 0.35 dex and 0.11 dex for the shock and nebular components, respectively. Assuming that the abundance discrepancy is produced by spatial variations in the electron temperature, we derive values of the temperature fluctuation parameter, t^2, of 0.050 and 0.016, for the shock and nebular components, respectively. Interestingly, we obtain almost coincident t^2 values for both components from the analysis of the intensity ratios of He I lines. We find significant departures from case B predictions in the Balmer and Paschen flux ratios of lines of high principal quantum number n. We analyze the ionization structure of HH202-S, finding enough evidence to conclude that the flow of HH202-S has compressed the ambient gas inside the nebula trapping the ionization front. We measure a strong increase of the total abundances of nickel and iron in the shock component, the abundance pattern and the results of photoionization models for both components are consistent with the partial destruction of dust after the passage of the shock wave in HH202-S.

Spectroscopy of planetary nebulae (PNe) provides the means to investigate s-process enrichments of neutron(n)-capture elements that cannot be detected in asymptotic giant branch (AGB) stars. However, accurate abundance determinations of these elements present a challenge. Corrections for unobserved ions can be large and uncertain, since in many PNe only one ion of a given n-capture element has been detected. Furthermore, the atomic data governing the ionization balance of these species are not well-determined, inhibiting the derivation of accurate ionization corrections. We present initial results of a program that addresses these challenges. Deep high resolution optical spectroscopy of ~20 PNe has been performed to detect emission lines from trans-iron species including Se, Br, Kr, Rb, and Xe. The optical spectral region provides access to multiple ions of these elements, which reduces the magnitude and importance of uncertainties in the ionization corrections. In addition, experimental and theoretical efforts are providing determinations of the photoionization cross-sections and recombination rate coefficients of Se, Kr, and Xe ions. These new atomic data will make it possible to derive robust ionization corrections for these elements. Together, our observational and atomic data results will enable n-capture element abundances to be determined with unprecedented accuracy in ionized nebulae.

The high-frequency-peaked BL-Lacertae object \objectname{1ES 0806+524}, at redshift z=0.138, was observed in the very-high-energy (VHE) gamma-ray regime by VERITAS between November 2006 and April 2008. These data encompass the two-, and three-telescope commissioning phases, as well as observations with the full four-telescope array. \objectname{1ES 0806+524} is detected with a statistical significance of 6.3 standard deviations from 245 excess events. Little or no measurable variability on monthly time scales is found. The photon spectrum for the period November 2007 to April 2008 can be characterized by a power law with photon index $3.6 \pm 1.0_{\mathrm{stat}} \pm 0.3_{\mathrm{sys}}$ between $\sim$300 GeV and $\sim$700 GeV. The integral flux above 300 GeV is $(2.2\pm0.5_{\mathrm{stat}}\pm0.4_{\mathrm{sys}})\times10^{-12}\:\mathrm{cm}^{2}\:\mathrm{s}^{-1}$ which corresponds to 1.8% of the Crab Nebula flux. Non contemporaneous multiwavelength observations are combined with the VHE data to produce a broadband spectral energy distribution that can be reasonably described using a synchrotron-self Compton model.

We analyze the physical conditions of the outflow seen in QSO 2359-1241 (NVSS J235953-124148), based on high resolution spectroscopic VLT observations. This object was previously studied using Keck/HIRES data. The main improvement over the HIRES results is our ability to accurately determine the number density of the outflow. For the major absorption component, level population from five different Fe II excited level yields n_H=10^4.4 cm^-3 with less than 20% scatter. We find that the Fe ii absorption arises from a region with roughly constant conditions and temperature greater than 9000 K, before the ionization front where temperature and electron density drop. Further, we model the observed spectra and investigate the effects of varying gas metalicities and the spectral energy distribution of the incident ionizing radiation field. The accurately measured column densities allow us to determine the ionization parameter log(U) = -2.4 and total column density of the outflow (log(N_H) = 20.6 cm^-2). Combined with the number density finding, these are stepping stones towards determining the mass flux and kinetic luminosity of the outflow, and therefore its importance to AGN feedback processes.

Radiative and Auger decay data have been calculated for modelling the K lines in ions of the nickel isonuclear sequence, from Ni$^+$ up to Ni$^{27+}$. Level energies, transition wavelengths, radiative transition probabilities, and radiative and Auger widths have been determined using Cowan's Hartree--Fock with Relativistic corrections (HFR) method. Auger widths for the third-row ions (Ni$^+$--Ni$^{10+}$) have been computed using single-configuration average (SCA) compact formulae. Results are compared with data sets computed with the AUTOSTRUCTURE and MCDF atomic structure codes and with available experimental and theoretical values, mainly in highly ionized ions and in the solid state.

The HFR and AUTOSTRUCTURE atomic structure codes are used to compute complete data sets of level energies, wavelengths, A-values, and radiative and Auger widths for K-vacancy states of the Ne, Mg, Si, S, Ar, and Ca isonuclear sequences. Ions with electron number N>9 are treated for the first time. Detailed comparisons with previous measurements and theoretical data for ions with N<=9 are carried out in order to estimate reliable accuracy ratings.

Spectroscopic studies of Planetary Nebulae (PNe) and H {\sc ii} regions have driven much development in atomic physics. In the last few years the combination of a generation of powerful observatories, the development of ever more sophisticated spectral modeling codes, and large efforts on mass production of high quality atomic data have led to important progress in our understanding of the atomic spectra of such astronomical objects. In this paper I review such progress, including evaluations of atomic data by comparisons with nebular spectra, detection of spectral lines from most iron-peak elements and n-capture elements, observations of hyperfine emission lines and analysis of isotopic abundances, fluorescent processes, and new techniques for diagnosing physical conditions based on recombination spectra. The review is directed toward atomic physicists and spectroscopists trying to establish the current status of the atomic data and models and to know the main standing issues.

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 study the nature of the [TiII] and [NiII] emission from the so-called strontium filament found in the ejecta of eta Carinae. To this purpose we employ multilevel models of the TiII and NiII systems which are used to investigate the physical condition of the filament and the excitation mechanisms of the observed lines. For the TiII ion, for which no atomic data was previously available, we carry out ab initio calculations of radiative transition rates and electron impact excitation rate coefficients. It is found that the observed spectrum is consistent with the lines being excited in a mostly neutral region with an electron density of the order of $10^7$ cm$^{-3}$ and a temperature around 6000 K. In analyzing three observations with different slit orientations recorded between March~2000 and November~2001 we find line ratios that change among various observations, in a way consistent with changes of up to an order of magnitude in the strength of the continuum radiation field. These changes result from different samplings of the extended filament, due to the different slit orientations used for each observation, and yield clues on the spatial extent and optical depth of the filament. The observed emission indicates a large Ti/Ni abundance ratio relative to solar abundances. It is suggested that the observed high Ti/Ni ratio in gas is caused by dust-gas fractionation processes and does not reflect the absolute Ti/Ni ratio in the ejecta of \etacar. We study the condensation chemistry of Ti, Ni and Fe within the filament and suggest that the observed gas phase overabundance of Ti

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.

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%.

We calculate the efficiency of iron K line emission and iron K absorption in photoionized models using a new set of atomic data. These data are more comprehensive than those previously applied to the modeling of iron K lines from photoionized gases, and allow us to systematically examine the behavior of the properties of line emission and absorption as a function of the ionization parameter, density and column density of model constant density clouds. We show that, for example, the net fluorescence yield for the highly charged ions is sensitive to the level population distribution produced by photoionization, and these yields are generally smaller than those predicted assuming the population is according to statistical weight. We demonstrate that the effects of the many strongly damped resonances below the K ionization thresholds conspire to smear the edge, thereby potentially affecting the astrophysical interpretation of absorption features in the 7-9 keV energy band. We show that the centroid of the ensemble of K$\alpha$ lines, the K$\beta$ energy, and the ratio of the K$\alpha_1$ to K$\alpha_2$ components are all diagnostics of the ionization parameter of our model slabs

We report extensive calculations of the decay properties of fine-structure K-vacancy levels in Fe X-Fe XVII. A large set of level energies, wavelengths, radiative and Auger rates, and fluorescence yields has been computed using three different standard atomic codes, namely Cowan's HFR, AUTOSTRUCTURE and the Breit-Pauli R-matrix package. This multi-code approach is used to the study the effects of core relaxation, configuration interaction and the Breit interaction, and enables the estimate of statistical accuracy ratings. The K-alpha and KLL Auger widths have been found to be nearly independent of both the outer-electron configuration and electron occupancy keeping a constant ratio of 1.53+/-0.06. By comparing with previous theoretical and measured wavelengths, the accuracy of the present set is determined to be within 2 mA. Also, the good agreement found between the different radiative and Auger data sets that have been computed allow us to propose with confidence an accuracy rating of 20% for the line fluorescence yields greater than 0.01. Emission and absorption spectral features are predicted finding good correlation with measurements in both laboratory and astrophysical plasmas.

A detailed analysis of the radiative and Auger de-excitation channels of K-shell vacancy states in Fe II-Fe IX has been carried out. Level energies, wavelengths, A-values, Auger rates and fluorescence yields have been calculated for the lowest fine-structure levels populated by photoionization of the ground state of the parent ion. Different branching ratios, namely K-alpha_2/K-alpha_1, K-beta/K-alpha, KLM/KLL, KMM/KLL, and the total K-shell fluorescence yields, omega_K, obtained in the present work have been compared with other theoretical data and solid-state measurements, finding good general agreement with the latter. The K-alpha_2/K-alpha_1 ratio is found to be sensitive to the excitation mechanism. From these comparisons it has been possible to estimate an accuracy of ~10% for the present transition probabilities.

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 discuss implementation of TOPbase and Iron Project opacities for stellar spectral codes. We use a technique employed by Peach, where a Boltzmann-averaged cross section is calculated for selected temperatures, and the opacity obtained from double interpolation in temperature and wavelength. It is straightforward to include {\it all} levels for which cross sections have been calculated. Boltzmann-averaged cross sections for Fe II show a local maximum between 1700 and 2000[A]. We suggest this feature arises from 3d^5 4snl to 3d^5 4pnl transitions within Fe II. IUE spectra of iron-rich CP stars show local minima in this region. Theoretical calculations of a representative stellar continuum demonstrate that Fe II photoionization contributes significantly to the observed minima.

As part of a project to compute improved atomic data for the spectral modeling of iron K lines, we report extensive calculations and comparisons of atomic data for K-vacancy states in Fe XXIV. The data sets include: (i) energy levels, line wavelengths, radiative and Auger rates; (ii) inner-shell electron impact excitation rates and (iii) fine structure inner-shell photoionization cross sections. The calculations of energy levels and radiative and Auger rates have involved a detailed study of orbital representations, core relaxation, configuration interaction, relativistic corrections, cancellation effects and semi-empirical corrections. It is shown that a formal treatment of the Breit interaction is essential to render the important magnetic correlations that take part in the decay pathways of this ion. As a result, the accuracy of the present A-values is firmly ranked at better than 10% while that of the Auger rates at only 15%. The calculations of collisional excitation and photoionization cross sections take into account the effects of radiation and spectator Auger dampings. In the former, these effects cause significant attenuation of resonances leading to a good agreement with a simpler method where resonances are excluded. In the latter, resonances converging to the K threshold display symmetric profiles of constant width that causes edge smearing.

It is shown that the commonly held view of a sharp Fe K edge must be modified if the decay pathways of the series of resonances converging to the K thresholds are adequately taken into account. These resonances display damped Lorentzian profiles of nearly constant widths that are smeared to impose continuity across the threshold. By modeling the effects of K damping on opacities, it is found that the broadening of the K edge grows with the ionization level of the plasma and that the appearance at high ionization of a localized absorption feature at 7.2 keV is identified as the K-beta unresolved transition array.

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.

We calculate the recombination spectra of the He-like ions He~I, C~V, N~VI, O~VII, Ne~IX, Mg~XI, Si~XIII, S~XV, Ar~XVII, Ca~XIX, and Fe~XXV. We include the following physical processes: radiative recombination, dielectronic recombination, three-body recombination, electron impact ionization, and collisional excitation by electrons, protons, and $\alpha$-particles. The calculations account for the effects of lowering of the continuum at high densities and high density corrections to dielectronic recombination. Then we construct models for He-like ions for fast computation of their spectra. Every model includes 29 bound levels up to n=5 and 6 doubly excited levels that account for the most important satellite lines. The models are constructed in a way that allows for proper approach to LTE under appropriate conditions. These models can simultaneously solve for the H/He-like ionization balance in photoionized or collisionally ionized plasmas and compute emission spectra including the combined effects of radiative and dielectronic recombination, collisional excitation, photoionization from excited levels, fluorescence, and line trapping. The models can be used for any temperature between 100 and $10^9$K and electron densities of up to $10^{18}$ \cm3. The models can be easily used within spectral modeling codes or as stand-alone tools for spectral analysis. We present comparisons between the results of the present models and previous work. Significant differences are found between the present effective recombination rate coefficients to the $n=2$ and those of previous estimates. Later, we study various emission line ratio diagnostics under collisional ionization and photoionized conditions.

Several different approximations and techniques have been developed for the calculation of atomic structure, ionization, and excitation of atoms and ions. These techniques have been used to compute large amounts of spectroscopic data of various levels of accuracy. This paper presents a review of these theoretical methods to help non-experts in atomic physics to better understand the qualities and limitations of various data sources and assess how reliable are spectral models based on those data.

The UV emission lines of Hercules X-1, resolved with the HST GHRS and STIS, can be divided into broad (FWHM 750 km/s) and narrow (FWHM 150 km/s) components. The broad lines can be unambiguously identified with emission from an accretion disk which rotates prograde with the orbit. The narrow lines, previously identified with the X-ray illuminated atmosphere of the companion star, are blueshifted at both phi=0.2 and phi=0.8 and the line flux at phi=0.2 is 0.2 of the flux at phi=0.8. Line ratio diagnostics show that the density of the narrow line region is log n=13.4+/-0.2 and the temperature is T=1.0+/-0.2x10^5 K. The symmetry of the eclipse ingress suggests that the line emission on the surface of the disk is left-right symmetric relative to the orbit. Model fits to the O V, Si IV, and He II line profiles agree with this result, but fits to the N V lines suggest that the receding side of the disk is brighter. We note that there are narrow absorption components in the N V lines with blueshifts of 500 km/s.

We present ground state photoionization cross sections of atoms and ions averaged over resonance structures for photoionization modeling of astrophysical sources. The detailed cross sections calculated in the close-coupling approximation using the R-matrix method, with resonances delineated at thousands of energies, are taken from the Opacity Project database TOPbase and the Iron Project, including new data for the low ionization stages of iron Fe I--V. The resonance-averaged cross sections are obtained by convolving the detailed cross sections with a Gaussian distribution over the autoionizing resonances. This procedure is expected to minimize errors in the derived ionization rates that could result from small uncertainties in computed positions of resonances, while preserving the overall resonant contribution to the cross sections in the important near threshold regions. The detailed photoionization cross sections at low photon energies are complemented by new relativistic distorted-wave calculations for Z<= 12, and from central-field calculations for Z>12 at high energies, including inner-shell ionization. The effective cross sections are then represented by a small number of points that can be readily interpolated linearly for practical applications; a Fortran subroutine and data are available. The present numerically averaged cross sections are compared with analytic fits that do not accurately represent the effective cross sections in regions dominated by resonances.

The emission spectra and the ionization structure of the low ionization stages of iron, Fe I--IV, in gaseous nebulae are studied. This work includes: (i) new atomic data: photoionization cross sections, total e-ion recombination rates, excitation collision strengths, and transition probabilities; (ii) detailed study of excitation mechanisms for the [Fe II], [Fe III], and [Fe IV] emission, and spectroscopic analysis of the observed IR, optical, and UV spectra; (iii) study of the physical structure and kinematics of the nebulae and their ionization fronts. Spectral analysis of the well observed Orion nebula is carried out as a test case, using extensive collisional-radiative and photoionization models. It is shown that the [Fe II] emission from the Orion nebula is predominantly excited via electron collisions in high density partially ionized zones; radiative fluorescence is relatively less effective. Further evidence for high density zones is derived from the [O I] and [Ni II] spectral lines, as well as from the kinematic measurements of ionic species in the nebula. The ionization structure of iron in Orion is modeled using the newly calculated atomic data, showing some significant differences from previous models. The new model suggests a fully ionized H II region at densities on the order of $10^3$ cm$^{-3}$, and a dynamic partially ionized H II/H I region at densities of $10^5-10^7$ \cm3. Photoionization models also indicate that the optical [O I] and [Fe II] emission originates in high density partially ionized regions within ionization fronts. The gas phase iron abundance in Orion is estimated from observed spectra.