Locally Optimally-emitting Clouds and the Narrow Emission Lines in Seyfert Galaxies

The narrow emission line spectra of active galactic nuclei are not accurately described by simple photoionization models of single clouds. Recent Hubble Space Telescope images of Seyfert 2 galaxies show that these objects are rich with ionization cones, knots, filaments, and strands of ionized gas. Here we extend to the narrow line region the ``locally optimally emitting cloud'' (LOC) model, in which the observed spectra are predominantly determined by powerful selection effects. We present a large grid of photoionization models covering a wide range of physical conditions and show the optimal conditions for producing many of the strongest emission lines. We show that the integrated narrow line spectrum can be predicted by an integration of an ensemble of clouds, and we present these results in the form of diagnostic line ratio diagrams making comparisons with observations. We also predict key diagnostic line ratios as a function of distance from the ionizing source, and compare these to observations. The predicted radial dependence of the [O III]/[O II] ratio may be matched to the observed one in NGC4151, if the narrow line clouds see a more intense continuum than we see. The LOC scenario when coupled with a simple Keplerian gravitational velocity field will quite naturally predict the observed line width versus critical density relationship. The influence of dust within the ionized portion of the clouds is discussed and we show that the more neutral gas is likely to be dusty, although a high ionization dust-free region is most likely present too. This argues for a variety of NLR cloud origins.Comment: 29 pages plus 16 figures, accepted for publication in Ap

Classification and analysis of emission-line galaxies using mean field independent component analysis

We present an analysis of the optical spectra of narrow emission-line galaxies, based on mean field independent component analysis (MFICA). Samples of galaxies were drawn from the Sloan Digital Sky Survey (SDSS) and used to generate compact sets of `continuum' and `emission-line' component spectra. These components can be linearly combined to reconstruct the observed spectra of a wider sample of galaxies. Only 10 components - five continuum and five emission line - are required to produce accurate reconstructions of essentially all narrow emission-line galaxies; the median absolute deviations of the reconstructed emission-line fluxes, given the signal-to-noise ratio (S/N) of the observed spectra, are 1.2-1.8 sigma for the strong lines. After applying the MFICA components to a large sample of SDSS galaxies we identify the regions of parameter space that correspond to pure star formation and pure active galactic nucleus (AGN) emission-line spectra, and produce high S/N reconstructions of these spectra. The physical properties of the pure star formation and pure AGN spectra are investigated by means of a series of photoionization models, exploiting the faint emission lines that can be measured in the reconstructions. We are able to recreate the emission line strengths of the most extreme AGN case by assuming the central engine illuminates a large number of individual clouds with radial distance and density distributions, f(r) ~ r^gamma and g(n) ~ n^beta, respectively. The best fit is obtained with gamma = -0.75 and beta = -1.4. From the reconstructed star formation spectra we are able to estimate the starburst ages. These preliminary investigations serve to demonstrate the success of the MFICA-based technique in identifying distinct emission sources, and its potential as a tool for the detailed analysis of the physical properties of galaxies in large-scale surveys.Comment: MNRAS accepted. 29 pages, 24 figures, 3 table

Interpreting the Ionization Sequence in AGN Emission-Line Spectra

We investigate the physical cause of the great range in the ionization level seen in the spectra of narrow lined active galactic nuclei (AGN). Mean field independent component analysis identifies examples of individual SDSS galaxies whose spectra are not dominated by emission due to star formation (SF), which we designate as AGN. We assembled high S/N ratio composite spectra of a sequence of these AGN defined by the ionization level of their narrow-line regions (NLR), extending down to very low-ionization cases. We used a local optimally emitting cloud (LOC) model to fit emission-line ratios in this AGN sequence. These included the weak lines that can be measured only in the co-added spectra, providing consistency checks on strong line diagnostics. After integrating over a wide range of radii and densities our models indicate that the radial extent of the NLR is the major parameter in determining the position of high to moderate ionization AGN along our sequence, providing a physical interpretation for their systematic variation. Higher ionization AGN contain optimally emitting clouds that are more concentrated towards the central continuum source than in lower ionization AGN. Our LOC models indicate that for the objects that lie on our AGN sequence, the ionizing luminosity is anticorrelated with the NLR ionization level, and hence anticorrelated with the radial concentration and physical extent of the NLR. A possible interpretation that deserves further exploration is that the ionization sequence might be an age sequence where low ionization objects are older and have systematically cleared out their central regions by radiation pressure. We consider that our AGN sequence instead represents a mixing curve of SF and AGN spectra, but argue that while many galaxies do have this type of composite spectra, our AGN sequence appears to be a special set of objects with negligible SF excitation.Comment: 57 pages; 18 figures, accepted by MNRA

Two-Component Structure of the Hbeta Broad-Line Region in Quasars. I. Evidence from Spectral Principal Component Analysis

We report on a spectral principal component analysis (SPCA) of a sample of 816 quasars, selected to have small Fe II velocity shifts with spectral coverage in the rest wavelength range 3500--5500 \AA. The sample is explicitly designed to mitigate spurious effects on SPCA induced by Fe II velocity shifts. We improve the algorithm of SPCA in the literature and introduce a new quantity, \emph{the fractional-contribution spectrum}, that effectively identifies the emission features encoded in each eigenspectrum. The first eigenspectrum clearly records the power-law continuum and very broad Balmer emission lines. Narrow emission lines dominate the second eigenspectrum. The third eigenspectrum represents the Fe II emission and a component of the Balmer lines with kinematically similar intermediate velocity widths. Correlations between the weights of the eigenspectra and parametric measurements of line strength and continuum slope confirm the above interpretation for the eigenspectra. Monte Carlo simulations demonstrate the validity of our method to recognize cross talk in SPCA and firmly rule out a single-component model for broad Hbeta. We also present the results of SPCA for four other samples that contain quasars in bins of larger Fe II velocity shift; similar eigenspectra are obtained. We propose that the Hbeta-emitting region has two kinematically distinct components: one with very large velocities whose strength correlates with the continuum shape, and another with more modest, intermediate velocities that is closely coupled to the gas that gives rise to Fe II emission.Comment: 22 pages, 17 figures, accepted for publication in The Astrophysical Journa

Interpreting the Ionization Sequence in Star-Forming Galaxy Emission-Line Spectra

High ionization star forming (SF) galaxies are easily identified with strong emission line techniques such as the BPT diagram, and form an obvious ionization sequence on such diagrams. We use a locally optimally emitting cloud model to fit emission line ratios that constrain the excitation mechanism, spectral energy distribution, abundances and physical conditions along the star-formation ionization sequence. Our analysis takes advantage of the identification of a sample of pure star-forming galaxies, to define the ionization sequence, via mean field independent component analysis. Previous work has suggested that the major parameter controlling the ionization level in SF galaxies is the metallicity. Here we show that the observed SF- sequence could alternatively be interpreted primarily as a sequence in the distribution of the ionizing flux incident on gas spread throughout a galaxy. Metallicity variations remain necessary to model the SF-sequence, however, our best models indicate that galaxies with the highest and lowest observed ionization levels (outside the range -0.37 < log [O III]/H\b{eta} < -0.09) require the variation of an additional physical parameter other than metallicity, which we determine to be the distribution of ionizing flux in the galaxy.Comment: 41 pages, 17 figures, 9 tables, accepted to MNRA

Outflows from active galactic nuclei: The BLR-NLR metallicity correlation

The metallicity of active galactic nuclei (AGNs), which can be measured by emission line ratios in their broad and narrow line regions (BLRs and NLRs), provides invaluable information about the physical connection between the different components of AGNs. From the archival databases of the International Ultraviolet Explorer, the Hubble Space Telescope and the Sloan Digital Sky Survey, we have assembled the largest sample available of AGNs which have adequate spectra in both the optical and ultraviolet bands to measure the narrow line ratio [N II]/H{\alpha} and also, in the same objects, the broad-line N V/C IV ratio. These permit the measurement of the metallicities in the NLRs and BLRs in the same objects. We find that neither the BLR nor the NLR metallicity correlate with black hole masses or Eddington ratios, but there is a strong correlation between NLR and BLR metallicities. This metallicity correlation implies that outflows from BLRs carry metal-rich gas to NLRs at characteristic radial distances of ~ 1.0 kiloparsec. This chemical connection provides evidence for a kinetic feedback of the outflows to their hosts. Metals transported into the NLR enhance the cooling of the ISM in this region, leading to local star formation after the AGNs turn to narrow line LINERs. This post-AGN star formation is predicted to be observable as an excess continuum emission from the host galaxies in the near infrared and ultraviolet, which needs to be further explored.Comment: 19 pages, 13 figures, 3 tables. Accepted for publication at MNRA

On the size of the Fe II emitting region in the AGN Akn 120

We present a reverberation analysis of the strong, variable optical Fe II emission bands in the spectrum of Akn 120, a low-redshift AGN which is one of the best candidates for such a study. On time scales of several years the Fe II line strengths follow the variations in the continuum strength. However, we are unable to measure a clear reverberation lag time for these Fe II lines on any time scale. This is due to the very broad and flat-topped nature of the Fe II cross correlation functions, as compared to the H-beta response which is much more sharply localized in time. Although there is some suggestion in the light curve of a 300-day response time, our statistical analysis does not pick up such a feature. We conclude that the optical Fe II emission does not come from a photoionization-powered region similar in size to the H-beta emitting region, but we cannot say for sure where it does come from. Our results are generally consistent either with emission from a photoionized region several times larger than the H-beta zone, or with emission from gas heated by some other means, perhaps responding only indirectly to the continuum variations.Comment: Accepted for publication in the Ap