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

We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures 2.8 K < T < 1e10 K and densities spanning the low to high-density limits. The radiation field between photon energies $10^{-8}$ Ry and 100 MeV is considered, along with all atoms and ions of the lightest 30 elements, and ~100 molecules. For ease of maintenance, the data are stored in a format as close as possible to the original data sources. Few data sources include the full range of data we need. We describe how we fill in the gaps in the data or extrapolate rates beyond their tabulated range. We tabulate data sources both for the atomic spectroscopic parameters and for collision data for the next release of Cloudy. This is not intended as a review of the current status of atomic data, but rather a description of the features of the database which we will build upon.

The detection of neutron(n)-capture elements in several planetary nebulae (PNe) has provided a new means of investigating s-process nucleosynthesis in low-mass stars. However, a lack of atomic data has inhibited accurate trans-iron element abundance determinations in astrophysical nebulae. Recently, photoionization and recombination data were determined for Se and Kr, the two most widely detected n-capture elements in nebular spectra. We have incorporated these new data into the photoionization code Cloudy. To test the atomic data, numerical models were computed for 15 PNe that exhibit emission lines from multiple Kr ions. We found systematic discrepancies between the predicted and observed emission lines that are most likely caused by inaccurate photoionization and recombination data. These discrepancies were removed by adjusting the Kr$^+$--Kr$^{3+}$ photoionization cross sections within their cited uncertainties and the dielectronic recombination rate coefficients by slightly larger amounts. From grids of models spanning the physical conditions encountered in PNe, we derive new, broadly applicable ionization correction factor (ICF) formulae for calculating Se and Kr elemental abundances. The ICFs were applied to our previous survey of near-infrared [Kr III] and [Se IV] emission lines in 120 PNe. The revised Se and Kr abundances are 0.1-0.3 dex lower than former estimates, with average values of [Se/(O, Ar)]=0.12$\pm$0.27 and [Kr/(O, Ar)]=0.82$\pm$0.29, but correlations previously found between their abundances and other nebular and stellar properties are unaffected. We also find a tendency for high-velocity PNe that can be associated with the Galactic thick disk to exhibit larger s-process enrichments than low-velocity PNe belonging to the thin disk population.

When modelling an ionised plasma, all spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum. Until now no data set exists that would meet these needs completely. We have therefore produced a table of relativistic Gaunt factors over a much wider range of parameter space than has ever been produced before. We present tables of the thermally averaged Gaunt factor covering the range log10(gamma^2) = -6 to 10 and log10(u) = -16 to 13 for all atomic numbers Z = 1 through 36. The data were calculated using the relativistic Bethe-Heitler-Elwert (BHE) approximation and were subsequently merged with accurate non-relativistic results in those parts of the parameter space where the BHE approximation is not valid. These data will be incorporated in the next major release of the spectral synthesis code Cloudy. We also produced tables of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to 10 for all values of Z between 1 and 36. All the data presented in this paper are available online.

Many massive galaxies at the centres of relaxed galaxy clusters and groups have vast reservoirs of cool (~10,000 K) and cold (~100 K) gas. In many low redshift brightest group and cluster galaxies this gas is lifted into the hot ISM in filamentary structures, which are long lived and are typically not forming stars. Two important questions are how far do these reservoirs cool and if cold gas is abundant what is the cause of the low star formation efficiency? Heating and excitation of the filaments from collisions and mixing of hot particles in the surrounding X-ray gas describes well the optical and near infra-red line ratios observed in the filaments. In this paper we examine the theoretical properties of dense, cold clouds emitting in the far infra-red and submillimeter through the bright lines of [C II]157 \mu m , [O I]63 \mu m and CO, exposed to these energetic ionising particles. While some emission lines may be optically thick we find this is not sufficient to model the emission line ratios. Models where the filaments are supported by thermal pressure support alone also cannot account for the cold gas line ratios but a very modest additional pressure support, either from turbulence or magnetic fields can fit the observed [O I]/[C II] line ratios by decreasing the density of the gas. This may also help stabilise the filaments against collapse leading to the low rates of star formation. Finally we make predictions for the line ratios expected from cold gas under these conditions and present diagnostic diagrams for comparison with further observations. We provide our code as an Appendix.

We announce a new facility in the spectral code CLOUDY that enables tracking the evolution of a cooling parcel of gas with time. For gas cooling from temperatures relevant to galaxy clusters, earlier calculations estimated the [Fe XIV] {\lambda}5303 / [Fe X] {\lambda}6375 luminosity ratio, a critical diagnostic of a cooling plasma, to slightly less than unity. By contrast, our calculations predict a ratio ~3. We revisit recent optical coronal line observations along the X-ray cool arc around NGC 4696 by Canning et al. (2011), which detected [Fe X] {\lambda}6375, but not [Fe XIV] {\lambda}5303. We show that these observations are not consistent with predictions of cooling flow models. Differential extinction could in principle account for the observations, but it requires extinction levels (A_V > 3.625) incompatible with previous observations. The non-detection of [Fe XIV] implies a temperature ceiling of 2.1 million K. Assuming cylindrical geometry and transonic turbulent pressure support, we estimate the gas mass at ~1 million solar masses. The coronal gas is cooling isochorically. We propose that the coronal gas has not condensed out of the intracluster medium, but instead is the conductive or mixing interface between the X-ray plume and the optical filaments. We present a number of emission lines that may be pursued to test this hypothesis and constrain the amount of intermediate temperature gas in the system.

Modern spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum emitted by an ionized plasma. Until now no set of data exists that would meet this need in a fully satisfactory way. We have therefore undertaken to produce a table of very accurate non-relativistic Gaunt factors over a much wider range of parameters than has ever been produced before. We first produced a table of non-averaged Gaunt factors, covering the parameter space log10(epsilon_i) = -20 to +10 and log10(w) = -30 to +25. We then continued to produce a table of thermally averaged Gaunt factors covering the parameter space log10(gamma^2) = -6 to +10 and log10(u) = -16 to +13. Finally we produced a table of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to +10. All the data presented in this paper are available online.

This work follows Lykins et al. discussion of classic plasma cooling function at low density and solar metallicity. Here we focus on how the cooling function changes over a wide range of density (n_H<10^12 cm^(-3)) and metallicity (Z<30Z _sun ). We find that high densities enhance the ionization of elements such as hydrogen and helium until they reach local thermodynamic equilibrium. By charge transfer, the metallicity changes the ionization of hydrogen when it is partially ionized. We describe the total cooling function as a sum of four parts: those due to H&He, the heavy elements, electron-electron bremsstrahlung and grains. For the first 3 parts, we provide a low-density limit cooling function, a density dependence function, and a metallicity dependence function. These functions are given with numerical tables and analytical fit functions. For grain cooling, we only discuss in ISM case. We then obtain a total cooling function that depends on density, metallicity and temperature. As expected, collisional de-excitation suppresses the heavy elements cooling. Finally, we provide a function giving the electron fraction, which can be used to convert the cooling function into a cooling rate.

A variety of astronomical environments contain clouds irradiated by a combination of isotropic and beamed radiation fields. For example, molecular clouds may be irradiated by the isotropic cosmic microwave background (CMB), as well as by a nearby active galactic nucleus (AGN). These radiation fields excite atoms and molecules and produce emission in different ways. We revisit the escape probability theorem and derive a novel expression that accounts for the presence of external radiation fields. We show that when the field is isotropic the escape probability is reduced relative to that in the absence of external radiation. This is in agreement with previous results obtained under ad hoc assumptions or with the two-level system, but can be applied to complex many-level models of atoms or molecules. This treatment is in the development version of the spectral synthesis code Cloudy. We examine the spectrum of a Spitzer cloud embedded in the local interstellar radiation field, and show that about 60 percent of its emission lines are sensitive to background subtraction. We argue that this geometric approach could provide an additional tool toward understanding the complex radiation fields of starburst galaxies.

Observations of metal-poor extragalactic H II regions allow the determination of the primordial helium abundance, Y_p. The He I emissivities are the foundation of the model of the H II region's emission. Porter, Ferland, Storey, & Detisch (2012) have recently published updated He I emissivities based on improved photoionization cross-sections. We incorporate these new atomic data and update our recent Markov Chain Monte Carlo analysis of the dataset published by Izotov, Thuan, & Stasinska (2007). As before, cuts are made to promote quality and reliability, and only solutions which fit the data within 95% confidence level are used to determine the primordial He abundance. The previously qualifying dataset is almost entirely retained and with strong concordance between the physical parameters. Overall, an upward bias from the new emissivities leads to a decrease in Y_p. In addition, we find a general trend to larger uncertainties in individual objects (due to changes in the emissivities) and an increased variance (due to additional objects included). From a regression to zero metallicity, we determine Y_p = 0.2465 +/- 0.0097, in good agreement with the Planck result of Y_p = 0.2485 +/- 0.0002. In the future, a better understanding of why a large fraction of spectra are not well fit by the model will be crucial to achieving an increase in the precision of the primordial helium abundance determination.

Errors were discovered in our recent calculations of HeI emissivities (Porter et al. 2012; arXiv:1206.4115). This corrigendum discusses the errors and reports corrected emissivities.

The Fe unresolved transition array (UTAs) produce prominent features in the 15-17?A wavelength range in the spectra of Active Galactic Nuclei (AGN). Here we present new calculations of the energies and oscillator strengths of inner- shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in Active Galactic Nuclei. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration.

This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered.

We discuss recent improvements in the calculation of the radiative cooling in both collisionally and photo ionized plasmas. We are extending the spectral simulation code Cloudy so that as much as possible of the underlying atomic data is taken from external databases, some created by others, some developed by the Cloudy team. This paper focuses on recent changes in the treatment of many stages of ionization of iron, and discusses its extensions to other elements. The H-like and He-like ions are treated in the iso-electronic approach described previously. Fe II is a special case treated with a large model atom. Here we focus on Fe III through Fe XXIV, ions which are important contributors to the radiative cooling of hot, 1e5 to 1e7 K, plasmas and for X-ray spectroscopy. We use the Chianti atomic database to greatly expand the number of transitions in the cooling function. Chianti only includes lines that have atomic data computed by sophisticated methods. This limits the line list to lower excitation, longer wavelength, transitions. We had previously included lines from the Opacity Project database, which tends to include higher energy, shorter wavelength, transitions. These were combined with various forms of the g-bar approximation, a highly approximate method of estimating collision rates. For several iron ions the two databases are almost entirely complementary. We adopt a hybrid approach in which we use Chianti where possible, supplemented by lines from the Opacity Project for shorter wavelength transitions. The total cooling including the lightest thirty elements differs significantly from some previous calculations.

The optical [N I] doublet near 5200 {\AA} is anomalously strong in a variety of emission-line objects. We compute a detailed photoionization model and use it to show that pumping by far-ultraviolet (FUV) stellar radiation previously posited as a general explanation applies to the Orion Nebula (M42) and its companion M43; but, it is unlikely to explain planetary nebulae and supernova remnants. Our models establish that the observed nearly constant equivalent width of [N I] with respect to the dust-scattered stellar continuum depends primarily on three factors: the FUV to visual-band flux ratio of the stellar population; the optical properties of the dust; and the line broadening where the pumping occurs. In contrast, the intensity ratio [N I]/H{\beta} depends primarily on the FUV to extreme-ultraviolet ratio, which varies strongly with the spectral type of the exciting star. This is consistent with the observed difference of a factor of five between M42 and M43, which are excited by an O7 and B0.5 star respectively. We derive a non-thermal broadening of order 5 km/s for the [N I] pumping zone and show that the broadening mechanism must be different from the large-scale turbulent motions that have been suggested to explain the line-widths in this H II region. A mechanism is required that operates at scales of a few astronomical units, which may be driven by thermal instabilities of neutral gas in the range 1000 to 3000 K. In an appendix, we describe how collisional and radiative processes are treated in the detailed model N I atom now included in the Cloudy plasma code.

Far ultraviolet emission has been detected from a knot of Halpha emission in the Horseshoe filament, far out in the NGC 1275 nebula. The flux detected relative to the brightness of the Halpha line in the same spatial region is very close to that expected from Hydrogen two-photon continuum emission in the particle heating model of Ferland et al. (2009) if reddening internal to the filaments is taken into account. We find no need to invoke other sources of far ultraviolet emission such as hot stars or emission lines from CIV in intermediate temperature gas to explain these data.

We update our prior work on the case B collisional-recombination spectrum of He I to incorporate \textit{ab initio} photoionisation cross-sections. This large set of accurate, self-consistent cross-sections represents a significant improvement in He I emissivity calculations because it largely obviates the piecemeal nature that has marked all modern works. A second, more recent set of \textit{ab initio} cross-sections is also available, but we show that those are less consistent with bound-bound transition probabilities than our adopted set. We compare our new effective recombination coefficients with our prior work and our new emissivities with those by other researchers, and we conclude with brief remarks on the effects of the present work on the He I error budget. Our calculations cover temperatures $5000 \le T_e \le 25000$ K and densities $10^1 \le n_e \le 10^{14}$ cm$^{-3}$. Full results are available online.

We have supplemented existing spectra of Barnard's Loop with high accuracy spectrophotometry of one new position. Cloudy photoionization models were calculated for a variety of ionization parameters and stellar temperatures and compared with the observations. After testing the procedure with recent observations of M43, we establish that Barnard's Loop is photoionized by four candidate ionizing stars, but agreement between the models and observations is only possible if Barnard's Loop is enhanced in heavy elements by about a factor of 1.4. Barnard's Loop is very similar in properties to the brightest components of the Orion-Eridanus Bubble and the Warm Ionized Medium (WIM). We are able to establish models that bound the range populated in low-ionization color-color diagrams (I([SII])/I(H{\alpha}) versus I([NII])/I(H{\alpha})) using only a limited range of ionization parameters and stellar temperatures. Previously established variations in the relative abundance of heavy elements render uncertain the most common method of determining electron temperatures for components of the Orion-Eridanus Bubble and the WIM based on only the I([NII])/I(H{\alpha}) ratio, although we confirm that the lowest surface brightness components of the WIM are on average of higher electron temperature. The electron temperatures for a few high surface brightness WIM components determined by direct methods are comparable to those of classical bright H II regions. In contrast, the low surface brightness HII regions studied by the Wisconsin H{\alpha} Mapper are of lower temperatures than the classical bright HII regions.

We consider an accretion-disc origin for the broad and luminous forbidden-line emission observed in ultraluminous X-ray (ULX) sources CXOJ033831.8-352604 and XMMU 122939.7+075333 in globular clusters hosted by elliptical galaxies NGC 1399 and NGC 4472, respectively. We will refer to the latter by the globular cluster name RZ2109. The first has strong [OIII] and [NII], the second only [OIII]. Both H$\alpha$ and H$\beta$ are very weak or undetected in both objects. We assume that the large line widths are due to Keplerian rotation around a compact object and derive expressions for maximum line luminosities. These idealized models require central masses $\gtrsim100$ and $\gtrsim30000\Msun$ for CXOJ033831.8-352604 and RZ2109, respectively. An independent, bootstrap argument for the total disc mass yields, for both systems, $M_{\mathrm{disc}}\gtrsim10^{-4}\Msun$ for a purely metallic disc (and two orders of magnitude larger for solar metallicities). If Roche-lobe overflow is implicated, viscous time-scales are $\gtrsim300$ yr. Standard disc theory then offers another limit on the central masses. Lobe radii for a $\sim1\Msun$ donor are $\gtrsim10^{13}$ cm. We therefore rule out Roche-lobe overflow of a white dwarf in both systems. Red giants could fill the necessary lobes. Whether they are too metal-poor to produce the strong forbidden lines without strong hydrogen emission is unclear.

We investigate consequences of the discovery that Fe II emission in quasars, one of the spectroscopic signatures of "Eigenvector 1", may originate in infalling clouds. Eigenvector 1 correlates with the Eddington ratio L/L_Edd so that Fe II/Hbeta increases as L/L_Edd increases. We show that the "force multiplier", the ratio of gas opacity to electron scattering opacity, is ~ 10^3 - 10^4 in Fe II-emitting gas. Such gas would be accelerated away from the central object if the radiation force is able to act on the entire cloud. As had previously been deduced, infall requires that the clouds have large column densities so that a substantial amount of shielded gas is present. The critical column density required for infall to occur depends on L/L_Edd, establishing a link between Eigenvector 1 and the Fe II/Hbeta ratio. We see predominantly the shielded face of the infalling clouds rather than the symmetric distribution of emitters that has been assumed. The Fe II spectrum emitted by the shielded face is in good agreement with observations thus solving several long-standing mysteries in quasar emission lines.

A number of recent works in astronomy and cosmology have relied upon theoretical He I emissivities, but we know of no effort to quantify the uncertainties in the atomic data. We analyze and assign uncertainties to all relevant atomic data, perform Monte Carlo analyses, and report standard deviations in the line emissivities. We consider two sets of errors, which we call "optimistic" and "pessimistic." We also consider three different conditions, corresponding to prototypical Galactic and extragalactic H II regions and the epoch of cosmological recombination. In the extragalactic H II case, the errors we obtain are comparable to or larger than the errors in some recent $Y_p$ calculations, including those derived from CMB observations. We demonstrate a systematic effect on primordial abundance calculations; this effect cannot be reduced by observing a large number of objects. In the cosmological recombination case, the errors are comparable to many of the effects considered in recent calculations.

(abridged) Non-ionizing stellar continua are a source of photons for continuum pumping in the hydrogen Lyman transitions. In the environments where these transitions are optically thick, deexcitation occurs through higher series lines, so that the flux in these lines has a fluorescent contribution in addition to recombination; in particular, Balmer emissivities are systematically enhanced above case B. The effectiveness of such mechanism in HII regions and the adequacy of photoionization models as a tool to study it are the two main focuses of this work. We find that photoionization models of H II regions illuminated by low-resolution population synthesis models significantly overpredict the fluorescent contribution to the Balmer lines. Conversely, photoionization models in which the non-ionizing part of the continuum is omitted or is not transferred underpredict the fluorescent contribution to the Balmer lines, producing a bias of similar amplitude in the opposite direction. In this paper, we carry out realistic estimations of the fluorescent Balmer intensity and discuss the variations to be expected as the simulated observational setup and the stellar population's parameters are varied. In all the cases explored, we find that fluorescent excitation provides a significant contribution. We also show that differential fluorescent enhancement may produce line-of-sight differences in the Balmer decrement, mimicking interstellar extinction. Fluorescent excitation emerges from our study as a small but important mechanism for the enhancement of Balmer lines, which should be taken into account in the abundance analysis of photoionized regions, particularly in the case of high-precision applications such as the determination of primordial helium.

[Abridged] Photoionization, whether by starlight or other sources, has difficulty in accounting for the observed spectra of the optical filaments that often surround central galaxies in large clusters. Our first paper examined whether heating by energetic particles or dissipative MHD wave can account for the observations. Here we include atomic and low-ionization regions. The model of the hydrogen atom, along with all elements of the H-like iso-electronic sequence, is now fully nl-resolved. We show how the predicted HI spectrum differs from the pure recombination case. The second update is to the rates for H^0 - H2 inelastic collisions. We now use the values computed by Wrathmall et al. The rates are often much larger. We calculate the chemistry, ionization, temperature, gas pressure, and emission-line spectrum for a wide range of gas densities and collisional heating rates. We assume that the filaments are magnetically confined and free to move along field lines so that the gas pressure is equal to that of the surrounding hot gas. A mix of clouds, some being dense and cold and others hot and tenuous, can exist. The observed spectrum will be the integrated emission from clouds with different densities and temperatures but the same pressure. We assume that the gas filling factor is given by a power law in density. The power-law index is set by matching the observed intensities of IR H2 lines relative to optical HI lines. We conclude that the filaments are heated by ionizing particles, either conducted in from surrounding regions or produced in situ by processes related to MHD waves.

The optical filaments found in many cooling flows in galaxy clusters consist of low density ($\sim 10^3 \pcc$) cool ($\sim 10^3$ K) gas surrounded by significant amounts of cosmic-ray and magnetic-field energy. Their spectra show anomalously strong low-ionization and molecular emission lines when compared with galactic molecular clouds exposed to ionizing radiation such as the Orion complex. Previous studies have shown that the spectra cannot be produced by O-star photoionization. Here we calculate the physical conditions in dusty gas that is well shielded from external sources of ionizing photons and is energized either by cosmic rays or dissipative MHD waves. Strong molecular hydrogen lines, with relative intensities similar to those observed, are produced. Selection effects introduced by the microphysics produce a correlation between the \htwo line upper level energy and the population temperature. These selection effects allow a purely collisional gas to produce \htwo emission that masquerades as starlight-pumped \htwo but with intensities that are far stronger. This physics may find application to any environment where a broad range of gas densities or heating rates occur.

Recent calculations of collisionless, Case-B, He I emissivities were performed by Bauman et al. (2005). The source code used in the calculation has been freely available online since that paper was published. A number of changes have been made to simplify the use of the code by third parties. Here I provide details on how to obtain, compile, and execute the program and interpret the results.

A complete model of helium-like line and continuum emission has been incorporated into the plasma simulation code Cloudy. All elements between He and Zn are treated, any number of levels can be considered, and radiative and collisional processes are included. This includes photoionization from all levels, line transfer including continuum pumping and destruction by background opacities, scattering, and collisional processes. The model is calculated self-consistently along with the ionization and thermal structure of the surrounding nebula. The result is a complete line and continuum spectrum of the plasma. Here we focus on the ions of the He I sequence and reconsider the standard helium-like X-ray diagnostics. We first consider semi-analytical predictions and compare these with previous work in the low-density, optically-thin limit. We then perform numerical calculations of helium-like X-ray emission (such as is observed in some regions of Seyferts) and predict line ratios as a function of ionizing flux, hydrogen density, and column density. In particular, we demonstrate that, in photoionized plasmas, the $R$-ratio, a density indicator in a collisional plasma, depends on the ionization fraction and is strongly affected by optical depth for large column densities. We also introduce the notion that the $R$-ratio is a measure of the incident continuum at UV wavelengths. The $G$-ratio, which is temperature-sensitive in a collisional plasma, is also discussed, and shown to be strongly affected by continuum pumping and optical depth as well. These distinguish a photoionized plasma from the more commonly studied collisional case.

We apply a recently developed theoretical model of helium emission to observations of both the Orion Nebula and a sample of extragalactic H II regions. In the Orion analysis, we eliminate some weak and blended lines and compare theory and observation for our reduced line list. With our best theoretical model we find an average difference between theoretical and observed intensities $<I_{predicted}/I_{observed}-1> = 6.5%$. We argue that both the red and blue ends of the spectrum may have been inadequately corrected for reddening. For the 22 highest quality lines, with $3499 {\AA}\le\lambda\le6678 {\AA}$, our best model predicts observations to an average of 3.8%. We also perform an analysis of the reported observational errors and conclude they have been underestimated. In the extragalactic analysis, we demonstrate the likelihood of a large systematic error in the reported data and discuss possible causes. This systematic error is at least as large as the errors associated with nearly all attempts to calculate the primordial helium abundance from such observations. Our Orion analysis suggests that the problem does not lie in the theoretical models. We demonstrate a correlation between equivalent width and apparent helium abundance of lines from extragalactic sources that is most likely due to underlying stellar absorption. Finally, we present fits to collisionless case-B He I emissivities as well as the relative contributions due to collisional excitations out of the metastable $2s {}^{3}S$ term.

We discuss new functionality of the spectral simulation code CLOUDY which allows the user to calculate grids with one or more initial parameters varied and formats the predicted spectra in the standard FITS format. These files can then be imported into the x-ray spectral analysis software XSPEC and used as theoretical models for observations. We present and verify a test case. Finally, we consider a few observations and discuss our results.

Determinations of the primordial helium abundance are used in precision cosmological tests. These require highly accurate He I recombination rate coefficients. Here we reconsider the formation of He I recombination lines in the low-density limit. This is the simplest case and it forms the basis for the more complex situation where collisions are important. The formation of a recombination line is a two-step process, beginning with the capture of a continuum electron into a bound state and followed by radiative cascade to ground. The rate coefficient for capture from the continuum is obtained from photoionization cross sections and detailed balancing, while radiative transition probabilities determine the cascades. We have made every effort to use today's best atomic data. Radiative decay rates are from Drake's variational calculations, which include QED, fine structure, and singlet-triplet mixing. Certain high-$L$ fine-structure levels do not have a singlet-triplet distinction and the singlets and triplets are free to mix in dipole-allowed radiative decays. We use quantum defect or hydrogenic approximations to include levels higher than those treated in the variational calculations. [Abstract clipped]

We calculate the He I case B recombination cascade spectrum using improved radiative and collisional data. We present new emissivities over a range of electron temperatures and densities. The differences between our results and the current standard are large enough to have a significant effect not only on the interpretation of observed spectra of a wide variety of objects but also on determinations of the primordial helium abundance.