CLOUDY is a large-scale spectral synthesis code designed to simulate fully physical conditions within an astronomical plasma and then predict the emitted spectrum. Here we describe version 90 (C90) of the code, paying particular attention to changes in the atomic database and numerical methods that have affected predictions since the last publicly available version, C84. The computational methods and uncertainties are outlined together with the direction future development will take. The code is freely available and is widely used in the analysis and interpretation of emission-line spectra. Web access to the Fortran source for CLOUDY, its documentation Hazy, and an independent electronic form of the atomic database is also described.

We present a complete set of analytic fits to the non-relativistic photoionization cross sections for the ground states of atoms and ions of elements from H through Si, and S, Ar, Ca, and Fe. Near the ionization thresholds, the fits are based on the Opacity Project theoretical cross sections interpolated and smoothed over resonances. At higher energies, the fits reproduce calculated Hartree-Dirac-Slater photoionization cross sections.

The similarity of quasar line spectra has been taken as an indication that the emission line clouds have preferred parameters, suggesting that the environment is subject to a fine tuning process. We show here that the observed spectrum is a natural consequence of powerful selection effects. We computed a large grid of photoionization models covering the widest possible range of cloud gas density and distance from the central continuum source. For each line only a narrow range of density and distance from the continuum source results in maximum reprocessing efficiency, corresponding to "locally optimally emitting clouds" (LOC). These parameters depend on the ionization and excitation potentials of the line, and its thermalization density. The mean QSO line spectrum can be reproduced by simply adding together the full family of clouds, with an appropriate covering fraction distribution. The observed quasar spectrum is a natural consequence of the ability of various clouds to reprocess the underlying continuum, and can arise in a chaotic environment with no preferred pressure, gas density, or ionization parameter.

We present graphically the results of several thousand photoionization calculations of broad emission line clouds in quasars, spanning seven orders of magnitude in hydrogen ionizing flux and particle density. The equivalent widths of 42 quasar emission lines are presented as contours in the particle density - ionizing flux plane for a typical incident continuum shape, solar chemical abundances, and cloud column density of $N(H) = 10^{23} cm^{-2}$. Results are similarly given for a small subset of emission lines for two other column densities ($10^{22} cm^{-2}$ and $10^{24} cm^{-2}$), five other incident continuum shapes, and a gas metallicity of 5 \Zsun. These graphs should prove useful in the analysis of quasar emission line data and in the detailed modeling of quasar broad emission line regions. The digital results of these emission line grids and many more are available over the Internet.Comment: 16 pages, LaTeX (AASTeX aaspp4.sty); to appear in the 1997 ApJS: full contents of the 9 photoionization grids presented in this paper may be found at http://www.pa.uky.edu/~korista/grids/grids.htm

We present new calculations and analytic fits to the rates of radiative recombination towards H-like, He-like, Li-like and Na-like ions of all elements from H through Zn (Z = 30). The fits are valid over a wide range of temperature, from 3 K to 10 9 K.