Interpretation of the emission line spectra of Seyfert 2 galaxies by multi-component photoionization models

We present multi-component photoionization models allowing for local density inhomogeneities in the NLR to interpret the emission line spectra of Seyfert 2 galaxies. This approach leads to a successful match of a large set of line intensities from the UV to the NIR. In particular, the hitherto elusive NIR features [SIII]9062+9531 as well as high-ionization lines like [FeVII] are consistently fitted. The predictions of CIII] and CIV are considerably improved. From the detailed analysis of single-component photoionization models we derive the minimal radial extent of the NLR and the necessary span in density. Furthermore, we determine constraints on suggestions made about the role of matter-bounded clouds, and on proposed explanations for large [OIII]4363/5007 ratios (the so-called `temperature problem'), and assess the usability of some emission-line ratios as indicators of the ionization parameter. We find that a systematic variation of the cloud column densities in a population of matter-bounded clouds is inconsistent with the trends and correlations exhibited by the emission lines in the diagnostic diagrams. Concerning the temperature problem, the only possibility that leads to an overall consistency with the strengths of all other observed emission lines is subsolar metal abundances (as compared to e.g. the presence of dust, the existence of a high-density component, or matter-bounded clouds). In addition, the consequences of the presence of (Galactic-ISM-like) dust internal to the clouds were investigated. These models alleviate the [OIII]-ratio problem but did not lead to overall consistent fits. In our final model series, the NLR is composed of a mixture of metal-depleted (0.5 x solar) clouds with a radius-independent range in densities (10^2 to 10^5 cm^-3) distributed over a range of distances from the nucleus (galactocentric radii from at least $\sim$ 10^{20} cm to 10^{21.5} cm, for $Q_{tot} = 10^{54}$ s^{-1}). In order to encompass the observed range of each line intensity relative to H$\beta$, it turns out to be necessary to vary the spectral energy distribution incident on the clouds, qualitatively confirming the findings of Ho et al. (1993). We found a successful continuum sequence by adding an increasing contribution of a hot black body ($T \approx 200000$ K) to a steep powerlaw ($\alpha_{uv-x} \approx -2$). These continua imply that low and high-excitation objects differ in the strength but not in the basic shape of the EUV bump.