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.

Damped Lyman-alpha (DLA) and sub-DLA absorbers in quasar spectra provide the most sensitive tools for measuring element abundances of distant galaxies. Estimation of abundances from absorption lines depends sensitively on the accuracy of the atomic data used. We have started a project to produce new atomic spectroscopic parameters for optical/UV spectral lines using state-of-the-art computer codes employing very broad configuration interaction basis. Here we report our results for Zn II, an ion used widely in studies of the interstellar medium (ISM) as well as DLA/sub-DLAs. We report new calculations of many energy levels of Zn II, and the line strengths of the resulting radiative transitions. Our calculations use the configuration interaction approach within a numerical Hartree-Fock framework. We use both non-relativistic and quasi-relativistic one-electron radial orbitals. We have incorporated the results of these atomic calculations into the plasma simulation code Cloudy, and applied them to a lab plasma and examples of a DLA and a sub-DLA. Our values of the Zn II {\lambda}{\lambda} 2026, 2062 oscillator strengths are higher than previous values by 0.10 dex. Cloudy calculations for representative absorbers with the revised Zn atomic data imply ionization corrections lower than calculated before by 0.05 dex. The new results imply Zn metallicities should be lower by 0.1 dex for DLAs and by 0.13-0.15 dex for sub-DLAs than in past studies. Our results can be applied to other studies of Zn II in the Galactic and extragalactic ISM.

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.

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.

Absorption-line spectroscopy is a powerful tool used to estimate element abundances in the nearby as well as distant universe. The accuracy of the abundances thus derived is, naturally, limited by the accuracy of the atomic data assumed for the spectral lines. We have recently started a project to perform the new extensive atomic data calculations used for optical/UV spectral lines in the plasma modeling code Cloudy using state-of-the-art quantal calculations. Here we demonstrate our approach by focussing on S II, an ion used to estimate metallicities for Milky Way interstellar clouds as well as distant damped Lyman-alpha (DLA) and sub-DLA absorber galaxies detected in the spectra of quasars and gamma-ray bursts (GRBs). We report new extensive calculations of a large number of energy levels of S II, and the line strengths of the resulting radiative transitions. Our calculations are based on the configuration interaction approach within a numerical Hartree-Fock framework, and utilize both non-ralativistic and quasirelativistic one-electron radial orbitals. The results of these new atomic calculations are then incorporated into Cloudy and applied to a lab plasma, and a typical DLA, for illustrative purposes. The new results imply relatively modest changes (~0.04 dex) to the metallicities estimated from S II in past studies. These results will be readily applicable to other studies of S II in the Milky Way and other galaxies.

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.

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.