next previous
Up: Ultraviolet variability of

5. Physical conditions of the gas: Photoionization models

The strong correlations obtained between the intensity of the different broad components and the UV flux, confirms that the main excitation mechanism of the gas is photo-ionization by a central source. To obtain better insight in the physical conditions and the spatial relations of these parts of the BLR which appear to constitute separate regions, we have made separate photo-ionization calculations for each of the regions to compare with the results described in the previous section.

We used the photo-ionization code CLOUDY (Ferland 1991: CLOUDY 80.06). The initial parameters needed by CLOUDY are: the shape and intensity of the incident continuum (tex2html_wrap_inline7311), the chemical composition, the distribution of the gas (radial extension and the run of density with radius) and the covering factor. We have assumed solar composition, an open (small radiative covering factor) spherical geometry and constant density. For the ionizing continuum we have used the observed continuum derived in Sect. 3.3 at four levels of brightness (Fig. 5 (click here)), corresponding to the tex2html_wrap_inline7313 of 23.27 (high), 10.04 (high-intermediate), 3.78 (low-intermediate) and 1.31 tex2html_wrap_inline7315 (low), respectively. To connect the frequency gap between the Radio and the FIR range, linear interpolation has been used. To account for the big blue bump, the optical-UV continuum has been parametrized with a power law tex2html_wrap_inline7317, fitting tex2html_wrap_inline7319 to the observed tex2html_wrap_inline7321 for each level. The cut-off energy is a critical factor in the input spectrum. Clavel et al.\ (1990) and Binette et al. (1989) suggested 0.83 and 0.73 ryd, respectively, based on photoionization models for Fairall-9. The simultaneous observation with IUE, ROSAT and GINGA (Fig. 5 (click here)c) shows however that such a cut-off frequency does not correspond with the observed instantaneous spectrum of Fairall-9. Although it might be possible that the spectrum incident on the BLR is different from that which we observe, at this stage we have no convincing arguments to apply a spectrum different from the observed one. Then, a considerably slower drop off is indicated and a cut-off at 3.5 ryd as shown in Fig. 17 (click here) was applied. In the X-rays and higher frequencies we used a power law similar to the UV-optical one, with a spectral index tex2html_wrap_inline7323 and an arbitrary cut-off at 0.1 MeV (consistent with the tex2html_wrap_inline7325 ray observations of other Seyfert 1 nuclei with CGRO). The constant was defined by the observed integrated flux at tex2html_wrap_inline7327. To bring all in an absolute scale we used tex2html_wrap_inline7329 and all fluxes have been corrected for redshift with z=0.0461. We have calculated models with a range of total hydrogen density of the gas, n, between tex2html_wrap_inline7335 and tex2html_wrap_inline7337 and for the total hydrogen column density, tex2html_wrap_inline7339, a range between tex2html_wrap_inline7341. The inner radius of the gas cloud, r, was chosen over a range from 50 to 600 light-days, consistent with the delays obtained for the BLR line components with respect to the UV continuum derived in Sect. 4.3.

With these initial conditions many models have been calculated, obtaining from them the line luminosity for a geometric covering factor of 1 (tex2html_wrap_inline7345). Comparing the results of the models for this value of tex2html_wrap_inline7347 with the observed luminosities, the real covering factor can be obtained. The ionization parameter, U, is determined from the incident continuum luminosity, the density and the distance assumed for each model. So, comparing the relative intensities of the lines obtained from the models with the observed ones, the best fitting model can be found.

  figure1098
Figure 17: Fit to the simultaneous observation IUE-ROSAT-GINGA. : Observation (between the UV continuum and the J band tex2html_wrap_inline7353 has been used). : Fit with a cut-off energy of 3.5 rydberg used as input for the photoionization models. : SED with a 0.63 rydberg cut-off energy. It is clear from these simultaneous observations that the tex2html_wrap_inline7355 cut-off energy proposed by Binette et al. (1989) and Clavel & Santos-Lleó (1990) is incompatible with the data

5.1. Comparison with the observations

To compare the models and the observations we have used the observed line ratios for the components for tex2html_wrap_inline7359/tex2html_wrap_inline7361, tex2html_wrap_inline7363/tex2html_wrap_inline7365 and tex2html_wrap_inline7367/tex2html_wrap_inline7369, imposing that one single model fits all three line ratios, at all brightness levels of the continuum. Since the results for all models with tex2html_wrap_inline7371 and tex2html_wrap_inline7373 are the same, we only consider the first. The conclusions for each broad component are:

The central component

The three line ratios are well described by models with the following characteristics distance (light-days); tex2html_wrap_inline7375; tex2html_wrap_inline7377:
tex2html_wrap_inline7379; tex2html_wrap_inline7381; tex2html_wrap_inline7383; tex2html_wrap_inline7385; tex2html_wrap_inline7387.
The results for these best fitting models are shown in Fig. 18 (click here). Since tex2html_wrap_inline7389 is nearly independent of the distance, also models with inner radius up to 600 light-days could be allowed too, but at these distances the other two ratios are too large at some continuum levels, restricting the distance range to tex2html_wrap_inline7391. We see that all models are consistent with the observed values of tex2html_wrap_inline7393 and tex2html_wrap_inline7395, but underestimate tex2html_wrap_inline7397, at bright continuum levels, since they converge to 1. The mean ratio of tex2html_wrap_inline7399, which increases with distance, for distances of 50, 100, 150, 200 and 250 light-days is, respectively, 12.6, 17.0, 21.7, 26.1 and 36.1. The ratio tex2html_wrap_inline7401 for the central component estimated from the data of Wamsteker et al. (1985) is tex2html_wrap_inline7403 (Sect. 4.3.3), corresponding to models with radius of 50 and 100 light-days.

  figure1128
Figure 18: Line ratios for the central component and the corresponding photoionization models. The different models are drawn with several types of lines over the observed ratios: distance (light-days) - tex2html_wrap_inline7407 (tex2html_wrap_inline7409 (tex2html_wrap_inline7411 (- - -), tex2html_wrap_inline7413 (tex2html_wrap_inline7415 - tex2html_wrap_inline7417 -), tex2html_wrap_inline7419 (tex2html_wrap_inline7421 tex2html_wrap_inline7423 tex2html_wrap_inline7425 tex2html_wrap_inline7427), tex2html_wrap_inline7429 (tex2html_wrap_inline7431 tex2html_wrap_inline7433 tex2html_wrap_inline7435 - tex2html_wrap_inline7437 tex2html_wrap_inline7439 tex2html_wrap_inline7441) and tex2html_wrap_inline7443 (--)

We show in Fig. 21 (click here) a comparison between the spectrum derived from the models and the individual line components with the observed mean spectrum. The mean -emission lines only- spectrum shown in Fig. 21 (click here) has been corrected by z=0.0461 and corresponds to an observed tex2html_wrap_inline7447 of tex2html_wrap_inline7449 tex2html_wrap_inline7451 and to a tex2html_wrap_inline7453 central component luminosity of tex2html_wrap_inline7455 (tex2html_wrap_inline7457. The comparison in Fig. 21 (click here) is thus for an high-intermediate level of the continuum, with an observed value tex2html_wrap_inline7459 tex2html_wrap_inline7461. For this level the models give U=0.003, 0.009, 0.039, 0.022 and 0.014 and corresponding covering factors of 12%, 5%, 4%, 4% and 3% for distances of 50, 100, 150, 200 and 250 light-days, respectively. These model spectra are similar, except perhaps the 50 light-days model, which produces a weak CIItex2html_wrap_inline7465, which is not observed. The tex2html_wrap_inline7467 and tex2html_wrap_inline7469 intensities are somewhat higher than observed. In Fig. 21 (click here)a we show the mean spectrum together with the model of the central component for 100 light-days, tex2html_wrap_inline7471 and tex2html_wrap_inline7473 model. The photoionization models with CLOUDY confine the region where the central component is produced to a distance from the ionizing source of tex2html_wrap_inline7475, with a column density of tex2html_wrap_inline7477, an hydrogen density of tex2html_wrap_inline7479 and a covering factor of tex2html_wrap_inline7481, with U between 0.003 and 0.039 (Table 14 (click here)). The distance derived from the models is thus consistent with the lower end of the range derived from the delay of this component with respect to the UV continuum, tex2html_wrap_inline7485.

13

  table1166
Table 14: Resulting parameters for the broad components

The blue component

The acceptable models for the observed ratios of this component are similar in their parameters range to the central component ones. They are characterized by (as above in Sect. 5.1.1):
tex2html_wrap_inline7627; tex2html_wrap_inline7629 ; tex2html_wrap_inline7631; tex2html_wrap_inline7633 and tex2html_wrap_inline7635.
Larger distances are excluded because they result in a too high tex2html_wrap_inline7637 ratio. The results for these models are shown with the observed ratios in Fig. 19 (click here). Similarly to the central component, all models reproduce the observed tex2html_wrap_inline7639/tex2html_wrap_inline7641 and tex2html_wrap_inline7643/tex2html_wrap_inline7645 ratios very well, but underestimate tex2html_wrap_inline7647/tex2html_wrap_inline7649. The mean spectrum has in the blue component a Lytex2html_wrap_inline7651 luminosity of tex2html_wrap_inline7653. For the comparable high-intermediate level of the continuum, the models have U=0.035, 0.089, 0.039, 0.022 and 0.014 and predict a covering factor of 5%, 2.9%, 2.5%, 2.3% and 2.1% for 50, 100, 150, 200 and 250 light-days, respectively. Compared with the mean spectrum, the spectra of the five models are very similar, and we show the 100 light-days, tex2html_wrap_inline7657 and tex2html_wrap_inline7659 model for the blue component with the observed mean spectrum (Fig. 21 (click here)b).

  figure1244
Figure 19: Line ratios for the blue component and the corresponding photoionization models. The different models are drawn with several types of lines over the observed ratios: distance (light-days) - tex2html_wrap_inline7663 (tex2html_wrap_inline7665 (tex2html_wrap_inline7667 (- - -), tex2html_wrap_inline7669 (tex2html_wrap_inline7671 - tex2html_wrap_inline7673 -), tex2html_wrap_inline7675 (tex2html_wrap_inline7677 tex2html_wrap_inline7679 tex2html_wrap_inline7681 tex2html_wrap_inline7683), tex2html_wrap_inline7685 (tex2html_wrap_inline7687 tex2html_wrap_inline7689 tex2html_wrap_inline7691 - tex2html_wrap_inline7693 tex2html_wrap_inline7695 tex2html_wrap_inline7697) and tex2html_wrap_inline7699 (--)

In summary, the models indicate that the gas producing the blue component of Lytex2html_wrap_inline7701, CIV, SiIV and MgII is situated at a distance of the ionizing source of tex2html_wrap_inline7703, with a column density of tex2html_wrap_inline7705, an hydrogen density of tex2html_wrap_inline7707 and a covering factor of 5 - 2%, where U is confined between 0.089 and 0.014. This range of parameters is similar to the central component one (Table 14 (click here)). The distance derived from the photoionization model for the blue component is consistent with the mean delay of this component respect to the UV continuum (tex2html_wrap_inline7713).

The red component

The models, which fit best to the observed ratios for this component are characterized by (as above in Sect. 5.1.1):
tex2html_wrap_inline7719; tex2html_wrap_inline7721; tex2html_wrap_inline7723; tex2html_wrap_inline7725 and tex2html_wrap_inline7727. Also here larger distances are excluded because they result in an overestimate the tex2html_wrap_inline7729 ratio. These models are shown in Fig. 20 (click here) with the observed ratios. All models fit perfectly to tex2html_wrap_inline7731/SiIV and tex2html_wrap_inline7733/tex2html_wrap_inline7735, but tex2html_wrap_inline7737/tex2html_wrap_inline7739 is always underestimated.

  figure1280
Figure 20: Lines ratio for the red component and the corresponding photoionization models. The different models are drawn with several types of lines over the observed ratios: distance (light-days) - tex2html_wrap_inline7743 (tex2html_wrap_inline7745 (tex2html_wrap_inline7747 (- - -), tex2html_wrap_inline7749 (tex2html_wrap_inline7751 - tex2html_wrap_inline7753 -), tex2html_wrap_inline7755 (tex2html_wrap_inline7757 tex2html_wrap_inline7759 tex2html_wrap_inline7761 tex2html_wrap_inline7763), tex2html_wrap_inline7765 (tex2html_wrap_inline7767 tex2html_wrap_inline7769 tex2html_wrap_inline7771 - tex2html_wrap_inline7773 tex2html_wrap_inline7775 tex2html_wrap_inline7777) and tex2html_wrap_inline7779 (--)

For this component we have been able to isolate also the tex2html_wrap_inline7781/tex2html_wrap_inline7783 ratio. Its observed value is low (between 2 and 8) which can only be explained with models that, at similar distances, have an order of magnitude lower hydrogen density and a very small column density (tex2html_wrap_inline7785). This might indicate that the BLR is composed of optically thin clouds besides of the optically thick ones considered above, which might also explain better the tex2html_wrap_inline7787 ratio, as have been recently suggested by Shields et al.\ (1995). However, even though the previous models (optically thick) seem to underestimate the tex2html_wrap_inline7789 line strength, we will use them since at the required low column density no consistent solution was feasible for all other line ratios. The mean spectrum has a red component Lytex2html_wrap_inline7791 luminosity of tex2html_wrap_inline7793. For the high-intermediate level of the continuum, the models have U=0.035, 0.0089, 0.039, 0.022 and 0.014 and a covering factor of 6%, 4%, 3.3%, 3.0% and 2.9% for 50, 100, 150, 200 and 250 light-days, respectively. Comparison with the mean spectrum shows that the model spectra are very similar, and in Fig. 21 (click here)c the model shown is for 100 light-days, tex2html_wrap_inline7797 and tex2html_wrap_inline7799. In summary, the photoionization model calculations indicate that the gas producing the red component of tex2html_wrap_inline7801, tex2html_wrap_inline7803, tex2html_wrap_inline7805 and tex2html_wrap_inline7807 is situated at a distance of the ionizing source of tex2html_wrap_inline7809, with tex2html_wrap_inline7811, tex2html_wrap_inline7813 and a covering factor of tex2html_wrap_inline7815 with U between 0.009 and 0.039 (Table 14 (click here)).

The models alone do not allow to distinguish the gas producing the red component from that producing the blue component, because the same parameter range is applicable for both (and very similar to those for the central component as well). On the other hand the delay of the red component of tex2html_wrap_inline7819 with respect to the continuum falls outside the range suggested by the models, in contrast to the results found for the central and blue component. One possible explanation for this is that the gas emitting the red component is situated along the line of sight to the observer or that the regions producing both the red and the blue component are at the same distance falling in toward the central source (see Sect. 7 for a more detailed explanation).

The difficulty associated with the tex2html_wrap_inline7821 ratio for all models can be explained if: (a) The carbon abundance is less than the solar one (see the next section); or (b) The continuum shape in tex2html_wrap_inline7823) between 13.6 and 48 eV (tex2html_wrap_inline7825 and 16, respectively) is not as steep as the exponential cut-off used for the input continuum spectrum (Fig. 17 (click here)), so that the fraction of ionizing photon density to which CIV responds, decreases with respect to the photon density which Lytex2html_wrap_inline7827 responds, keeping in mind that the cut-off frequencies suggested by Clavel & Santos-Lleó (1990) and Binette et al.\ (1989) are inconsistent with the observed soft X-rays excess (Fig. 17 (click here)); or (c) There is not a smooth connection between the UV and the soft X-rays, with an absorption in the continuum between tex2html_wrap_inline7829, resulting in an increase in the tex2html_wrap_inline7831 ratio. Such absorption can be present if optically thin material, transparent at tex2html_wrap_inline7833 and absorbing at tex2html_wrap_inline7835, is present between the ionizing source and the BLR (Ferland et al. 1990); or (d) It is necessary to introduce in the models two populations of clouds: optically thin and optically thick ones, as it has been shown by Shields et al.\ (1995) to be possible for the specific case of F-9.

  figure1323
Figure 21: Observed mean spectrum and the obtained spectra from the photoionization models for each component. a) The observed mean spectrum and the model with 100 light-days, tex2html_wrap_inline7837 and tex2html_wrap_inline7839 for the central component. b) The observed mean spectrum and the model with 100 light-days, tex2html_wrap_inline7841 and tex2html_wrap_inline7843 for the blue component. c) Equal to the central component for the red one. d) The sum of each model for the different components over the observed mean spectrum. Take into account that the narrow component has not been fitted


next previous
Up: Ultraviolet variability of

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr