A&A 434, 569-584 (2005)
DOI: 10.1051/0004-6361:20047138

Simultaneous X-ray and UV spectroscopy of the Seyfert galaxy NGC 5548

II. Physical conditions in the X-ray absorber[*]

K. C. Steenbrugge1 - J. S. Kaastra1 - D. M. Crenshaw2 - S. B. Kraemer3,4 - N. Arav5 - I. M. George6,7 - D. A. Liedahl8 - R. L. J. van der Meer1 - F. B. S. Paerels9 - T. J. Turner6,7 - T. Yaqoob7,10


1 - SRON National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands
2 - Department of Physics and Astronomy, Georgia State University, Astronomy Offices, One Park Place South SE, Suite 700, Atlanta, GA 30303, USA
3 - Catholic University of America, USA
4 - Laboratory for Astronomy and Solar Physics, NASA's Goddard Space Flight Center, Code 681, Greenbelt, MD 20771, USA
5 - CASA, University of Colorado, 389 UCB, Boulder, CO 80309-0389, USA
6 - Joint Center for Astrophysics, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA
7 - Laboratory for High Energy Astrophysics, Code 660, NASA's Goddard Space Flight Center, Greenbelt, MD 20771, USA
8 - Physics Department, Lawrence Livermore National Laboratory, PO Box 808, L-41, Livermore, CA 94550, USA
9 - Columbia Astrophysics Laboratory, Columbia University, 538W. 120th Street, New York, NY 10027, USA
10 - Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA

Received 26 January 2004 / Accepted 4 January 2005

Abstract
We present the results from a 500 ks Chandra observation of the Seyfert 1 galaxy NGC 5548. We detect broadened (full width half maximum = 8000 km s-1) emission lines of O VII and C VI in the spectra, similar to those observed in the optical and UV bands. The source was continuously variable, with a 30% increase in luminosity in the second half of the observation. The gradual increase in luminosity occurred over a timescale of $\sim$300 ks. No variability in the warm absorber was detected between the spectra from the first 170 ks and the second part of the observation. The longer wavelength range of the LETGS resulted in the detection of absorption lines from a broad range of ions, in particular of C, N, O, Ne, Mg, Si, S and Fe. The velocity structure of the X-ray absorber is consistent with the velocity structure measured simultaneously in the ultraviolet spectra. We find that the highest velocity outflow component, at -1040 km s-1, becomes increasingly important for higher ionization parameters. This velocity component spans at least three orders of magnitude in ionization parameter, producing both highly ionized X-ray absorption lines (Mg XII, Si XIV) as well as UV absorption lines. A similar conclusion is very probable for the other four velocity components.

Based upon our observations, we argue that the warm absorber probably does not manifest itself in the form of photoionized clumps in pressure equilibrium with a surrounding wind. Instead, a model with a continuous distribution of column density versus ionization parameter gives an excellent fit to our data. From the shape of this distribution and the assumption that the mass loss through the wind should be smaller than the accretion rate onto the black hole, we derive upper limits to the solid angle as small as 10-4 sr. From this we argue that the outflow occurs in density-stratified streamers. The density stratification across the stream then produces the wide range of ionization parameter observed in this source. We determine an upper limit of 0.3 $M_{\odot}$ yr-1 for the mass loss from the galaxy due to the observed outflows.

Key words: galaxies: Seyfert - quasars: individual: NGC 5548 - galaxies: active - X-rays: galaxies

  
1 Introduction

Over half of all Seyfert 1 galaxies exhibit signatures of photoionized outflowing gas in their X-ray and UV spectra. Studying these outflows is important for a better understanding of the enrichment of the Inter Galactic Medium (IGM) as well as the physics of accretion of gas onto a super-massive black hole.

One of the best candidates to study these processes is NGC 5548. A recent XMM-Newton Reflection Grating Spectrometer (RGS) observation of NGC 5548 showed evidence for an X-ray warm absorber that has a nearly continuous distribution of ionization states (Steenbrugge et al. 2003).There is a pronounced increase in column density as a function of ionization parameter. The average outflow velocity measured in the X-ray band is consistent with the range of outflow velocities measured in the UV band, however, the spectral resolution in the X-ray band is not sufficient to resolve the individual velocity components observed in the UV band. The detection of the 1s-2p line of O VI in X-rays allows for a direct comparison with the column density as determined from measurements of the 2s-2p lines by FUSE in the UV band (Brotherton et al. 2002). From this direct comparison, as well as the total column density derived from other lowly ionized ions, Arav et al. (2003) and Steenbrugge et al. (2003) concluded that there is substantially more lowly ionized gas than has been claimed from previous UV observations. It was concluded that the X-ray and UV warm absorbers are different manifestations of the same phenomenon. Mathur et al. (1995) already noted this possibility before the advent of high resolution X-ray spectroscopy. However, their model included only one ionization component and this is unrealistic for the observed NGC 5548 high resolution spectra.

The difference in derived column densities in the UV band versus X-rays can be at least partially explained if the UV absorption lines do not cover the narrow emission line (NEL) region. The column densities derived from the UV lines are then only lower limits due to saturation effects, because these lines then locally absorb almost all the radiation from the continuum and the broad emission line (BEL) region (Arav et al. 2002, 2003).

The discrepancy in the O VI column density mentioned above is based on non-simultaneous X-ray and UV observations and therefore it remains possible that the difference is due to variability in the absorber. The present Chandra observation was proposed to remedy this, by having simultaneous HST STIS, FUSE and Chandra observations. Unfortunately, due to technical problems FUSE was unable to observe at the scheduled time, so that no simultaneous O VI observations were obtained. However, the spectral resolution of Chandra allows us to resolve the -1040 km s-1 component from the four other velocity components detected in the UV, for the strongest lines. This is a rather stringent test of the kinematic relation between both absorbers.

The long wavelength range of the LETGS allows us to map out the column density as a function of ionization parameter. Using only iron to avoid abundance effects we can sample about three orders of magnitude in ionization parameter from Fe VI to Fe XXIV. Due to the high signal-to-noise ratio, we are sensitive to changes in ionization parameter as small as 0.15 in log $\xi $. This is the first Seyfert 1 spectrum where the signal to noise is such that we are able to study in detail the absorption lines between 60-100 Å.

The details of the observations and the data reduction are given in Sect. 2. In Sect. 3 we present the spectral data analysis. The short and long term spectral variability of the warm absorber are discussed in Sect. 4. In Sect. 5 we discuss our results. In Sect. 6 we discuss the ionization structure, the geometry and the mass loss through the outflow discussed in Sect. 7. The continuum time variability is discussed in detail in Kaastra et al. (2004). Limits on the spatial distribution of the X-ray emission are given by Kaastra et al. (2003). The data were taken nearly simultaneously with HST STIS observations detailed by Crenshaw et al. (2003).

  
2 Observation and data reduction

NGC 5548 was observed for a full week with the Chandra High Energy Transmission Grating Spectrometer (HETGS) and the Low Energy Transmission Grating Spectrometer (LETGS) in January 2002. In Table 1 the instrumental setup and the exposure times are listed.

Table 1: The exposure time and the observation details for the Chandra and HST observations of NGC 5548. ACIS stands for the Advanced CCD Imaging Spectrometer and HRC for High Resolution Camera. The E140M grating covers the wavelength range between 1150-1730 Å, the E230M grating between 1607-3119 Å.

The LETGS data were reduced as described by Kaastra et al. (2002a). The LETGS spans a wavelength range from 1-180 Å with a resolution of 0.05 Å (full width half maximum, FWHM). The data were binned to 0.5 FWHM. However, due to Galactic absorption toward NGC 5548, the spectrum is heavily absorbed above 60 Å, and insignificant above 100 Å. In the present analysis we rebinned the data between 60-100 Å by a factor of 2 to the FWHM of the instrument and ignore the longer wavelengths.

The HETGS data were reduced using the standard CIAO software version 2.2. The HETGS spectra consist of a High Energy Grating (HEG) spectrum which covers the wavelength range of 1.5-13 Å with a FWHM of 0.012 Å, and a Medium Energy Grating (MEG) spectrum which covers the wavelength range of 1.5-24 Å with a FWHM of 0.023 Å. The MEG and HEG data are binned to 0.5 FWHM. In the reduction of both HETGS datasets the standard ACIS contamination model was applied (see http://cxc.harvard.edu/caldb/about_CALDB/index.html).

The data shortward of 1.5 Å were ignored for all the instruments, as were the data above 11.5 Å for the HEG spectrum and above 24 Å for the MEG spectrum. The higher orders are not subtracted from the LETGS spectrum, but are included in the spectral fitting, i.e. the response matrix used includes the higher orders up to order $\pm$10, with the relative strengths of the different orders taken from ground calibration and in-flight calibrations of Capella. We did not analyze the higher order HETGS spectra, due to their low count rate. All the spectra were analyzed with the SPEX software (Kaastra et al. 2002b), and the quoted errors are the 1$\sigma $ rms uncertainties, thus $\Delta\chi^2$ = 2.

Kaastra et al. (2002a) found a difference in the wavelength scale between the MEG and LETGS instrument of about 170 km s-1. This wavelength difference is within the absolute calibration uncertainty of the MEG instrument. The best line in NGC 5548 to test the wavelength calibration is the O VII forbidden emission line, as it is found to be narrow. However, due to the noise in the MEG spectrum, the wavelength is not well determined, but it is consistent with the wavelength measured from the LETGS spectrum. We also looked for differences in the line centroid of the O VIII Ly$\alpha $ and O VIII Ly$\beta$ absorption line. However, as these lines are broadened, centroids are less accurate. These lines have consistent centroids for both instruments, and no systematic blue- or redshift is detected.

  
3 Spectral data analysis

In Appendix C the spectrum between 1.5 Å-100 Å is shown. The X-ray spectrum is qualitatively similar to the spectra observed before by Chandra (Kaastra et al. 2000, 2002a; Yaqoob et al. 2001) and XMM-Newton (Steenbrugge et al. 2003). A large number of strong and weak absorption lines cover a smooth continuum. In addition, a few narrow and broadened emission lines are visible. Most absorption lines can be identified with lines already observed in these previous spectra or with predicted lines that, due to the increased sensitivity of the present spectra, become visible for the first time.

The spectral model applied here is based on results obtained from earlier Chandra and XMM-Newton spectra of NGC 5548 (Kaastra et al. 2002a; Steenbrugge et al. 2003). The current spectrum is described by the same spectral components as in the earlier observations, albeit with different parameters. In the next subsections we discuss the results for the different components as obtained from our global fit to the spectra. These previous observations as well as the present observation have a continuum that is well fitted by a power-law and modified black body component (Kaastra & Barr 1989) modified by Galactic extinction and cosmological redshift. The Galactic H I column density was frozen to 1.65 $\times$ 1024 m-2 (Nandra et al. 1993), as was the redshift of 0.01676 from the optical emission lines (Crenshaw & Kraemer 1999). The continuum parameters for our best fit model, namely the one with three slab components (Sect. 3.2.1, model A), are listed in Table 2. For comparison, the continuum parameters for the earlier RGS observation are included (Steenbrugge et al. 2003).

The high resolution spectra reveal a warm absorber (WA) modeled by combinations of slab, xabs or warm components (Sect. 3.1), narrow emission lines (Sects. 3.4 and 3.5) and broadened O VIII, O VII and C VI emission lines (Sect. 3.6).

Table 2: The best fit parameters for the continuum for the earlier RGS observation (Steenbrugge et al. 2003) and the present HETGS and LETGS observations of NGC 5548. The parameters for the HETGS and LETGS are for the best fit model including the three  slab components with the outflow velocities frozen to the UV values.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.8cm,clip]{H0138F1.ps}\end{figure} Figure 1: The SED used (solid line) in the present analysis and the analysis of the HST STIS data by Crenshaw et al. (2003). For comparison, the SED used by Kaastra et al. (2002a) and Steenbrugge et al. (2003) in earlier papers on NGC 5548 is plotted as the dot-dash line. The SED adopted by Dumont et al. (1998) is plotted as the dotted line.
Open with DEXTER

  
3.1 Warm absorber models

In our modeling of the warm absorber we use combinations of three different spectral models for the absorber. The most simple is the slab model (Kaastra et al. 2002a), which calculates the transmission of a slab of a given column density and a set of ion concentrations. Both the continuum and the line absorption are taken into account. As in the other absorption models, the lines are modeled using Voigt profiles. Apart from the ionic column densities, the average outflow velocity v and the Gaussian rms velocity broadening $\sigma $ are free parameters of the model. The last two parameters have the same value for each ion. The slab is assumed to be thin in the sense that  ${\rm e}^{-\tau}$ gives a fair description of its transmission.

The second model (xabs) is similar to the slab model, however now the ionic column densities are coupled through a set of photoionization balance calculations (Kaastra et al. 2002a). As a result, instead of the individual ionic concentrations we now have the ionization parameter $\xi $ = L/ $n_{\rm e}r^2$ and the elemental abundances as free parameters. The set of photoionization balance calculations was done using the XSTAR code (Kallman & Krolik 1999), with the SED as shown in Fig. 1. The same SED was used by Crenshaw et al. (2003) in the analysis of the simultaneous HST STIS data. As there is no evidence of a blue bump this was not included in the currently used SED. This is the main difference with the SED used by Kaastra et al. (2002a) and Steenbrugge et al. (2003), which is plotted as the dot-dash line. For the X-ray part of the spectrum we used the continuum model as measured from the present LETGS spectrum, instead of the continuum obtained from the 1999 LETGS spectrum. For the high energy part of the SED we used a power law component as measured in our Chandra spectrum but with an exponential cut-off at 130 keV, consistent with the BeppoSAX data (Nicastro et al. 2000). As reflection components are time variable and the LETGS has insufficient high energy sensitivity to model it properly, we did not include any reflection component in our SED. Between 50 and 1000 Å we used a power-law and at even longer wavelength used the energy index obtained by Dumont et al. (1998) of 2.5.

The third method used is the warm model, which is a model for a continuous distribution of the d$N_{\rm H}$/d$\xi $ as a function of $\xi $. At each value of $\xi $, the transmission is calculated using the xabs model described above. In our implementation of this model  $\xi {\rm d}N_{\rm H}/{\rm d}\xi$ is defined at a few (in our case two) values of $\xi $. At the intermediate points the logarithm of this function is determined by cubic spline interpolation on the log $\xi $ grid. This determines $N_{\rm H}$ as a function of $\xi $. Integration over ${\rm d}N_{\rm H}$/d$\xi $ between the lowest and the highest values of $\xi $ ($\xi_1$ respectively $\xi_2$) yields the total ionic column densities. In our model, as we have taken only two points, the total hydrogen column density is therefore a power-law function of ionization. The warm model correctly incorporates the fact that ions are formed at a range of ionization parameters. Free parameters of this model are $\xi_1$, $\xi_2$, $\xi_1$ d $N_{\rm H}~(\xi_1)$/d$\xi $, $\xi_2$ d $N_{\rm H}~(\xi_2)$/d$\xi $, and similar to the xabs model the outflow velocity v, velocity broadening $\sigma $ and the elemental abundances.

Finally, in all our models (slab, xabs and warm) we changed the wavelength for the O V 1s22s2-1s2s22p 1P1 X-ray line from 22.33 Å (the HULLAC value) to 22.374 Å (Schmidt et al. 2004) as measured at the University of California EBIT-I electron beam ion trap. This recent value was not yet available for our earlier analysis of the RGS data (Steenbrugge et al. 2003).

3.2 Spectral fits with warm absorber models

  
3.2.1 Warm absorber model A: Column densities with UV velocity structure

Our first method (dubbed model A here) allows us to further ascertain that the X-ray warm absorber and the UV absorber are the same phenomenon, a major goal for proposing this observation. We fitted the warm absorber in the X-ray spectra with the outflow velocities frozen to those measured in the UV spectra.

We implemented this by modeling the WA using three slab components, with the outflow velocities frozen. These outflow velocities were frozen to -1040 km s-1, -530 km s-1 and -160 km s-1, corresponding to UV velocity components 1, the average of 2-4, and 5, respectively. Note that the velocities of components 2-4 (-667, -530 and -336 km s-1, respectively) are too close to be separated in our X-ray spectra. In the fit we left the velocity broadening of each of the three components a free parameter. Allowing for three velocity components instead of one, substantially improved the fit from $\chi^2$ = 2916 to 2441 for 2160 degrees of freedom (only LETGS is quoted here). However, the absorption lines are still only partially resolved, leading to strongly correlated errors for the derived column densities of the velocity components. Only for the six ions with the deepest lines (C VI, O V, O VII, O VIII, Ne IX and Si XIII) we measure column densities with meaningful error ranges (Table 3). For these six ions the velocity structure of the X-ray warm absorber closely resembles the UV velocity structure. Figure 2 shows the line profiles as a function of outflow velocity for four lines of the ions listed in Table 3. The HEG spectrum shows most clearly substructure in the Mg XII Ly$\alpha $ absorption lines (Fig. C.1 in Appendix C).

Table 3: Parameters for six ions measured using model A. The first row gives the component number as listed in the UV. The second row lists the outflow velocity v, which was frozen to the UV values, the third row lists the velocity broadening $\sigma _{\rm X-ray}$ as derived from the X-ray spectra. In the fourth row the velocity broadening $\sigma _{\rm UV}$ = FWHM/2.35, where FWHM is the measured and resolved Full Width Half Maximum of the UV absorption lines (Crenshaw et al. 2003). All are in km s-1. For the -530 km s-1 component we added the velocity broadening of components 2 to 4. We list the logarithms of the column densities in m-2. For comparison we also list the C IV and N V column densities as measured by Crenshaw et al. (2003).


  \begin{figure}
\par\includegraphics[angle=-90,width=8.75cm,clip]{H0138F2.ps}\end{figure} Figure 2: The MEG (open circles) and LETGS (filled squares) line profiles for the deepest line of O VIII, O VII, O V and C VI. The five outflow velocities measured in the UV are indicated by the dotted line. No correction was made for possible blending.
Open with DEXTER

  
3.2.2 Warm absorber model B: Column densities with simple velocity structure
Free velocity broadening

Table 4: The best fit column densities using model B for the ions for which the outflow velocity and velocity broadening are well constrained. The outflow velocity v and the rms velocity broadening $\sigma $ are listed as well as the ionization parameters $\xi $ and U for which the ion has its maximum column density. The column densities are from the fit including the broad emission lines (Sect. 3.6).

Using model A we can only derive accurate column densities for six ions. In order to obtain column densities for other ions we need to reduce the number of free parameters. To do so we use a model with one outflow velocity component, leaving the velocity broadening as well as the outflow velocity free. This allows us to obtain column densities, outflow velocities and velocity broadening for the ions with stronger lines. We first made a fit with a single slab component. This slab component has a single outflow velocity and velocity broadening for all the ions. After this fit is obtained, we freeze all the parameters in this slab component.

A second slab component is then added containing only a single ion, leaving the column density, outflow velocity and velocity broadening as free parameters. The column density for this particular ion is frozen to zero in the first slab component, and the spectra are fit again in order to fine tune the column density, outflow velocity and velocity broadening for each ion. Each ion was cycled through the second slab component in this manner, and the results for those ions which have well determined outflow velocity and velocity broadening are listed in Table 4. In this table we also list the ionization parameter for which the column density of the particular ion peaks. We list $\xi $ = L/ $n_{\rm e}r^2$ expressed in 10-9 W m, where L is the energy luminosity, ne is the electron density, r the distance. We also list U = $Q_{\rm tot}$/(4$\pi$r2c$n_{\rm e}$), with $Q_{\rm tot}$/(4$\pi$r2) the total photon flux for photons with an energy greater than 13.6 eV. This table indicates a weak correlation between outflow velocity and the ionization of the ion, which is shown in Fig. 3. As a side effect, leaving the outflow velocity free potentially allows to correct for any remaining wavelength inaccuracies.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.7cm,clip]{H0138F3.ps}\end{figure} Figure 3: The measured outflow velocity versus the logarithm of the ionization parameter. The dotted line indicates the -530 km s-1 outflow velocity, the dominant component in the UV absorber. A trend toward higher outflow velocities for higher ionized ionization parameters is observed.
Open with DEXTER

Fixed velocity broadening

As it is not possible to determine the outflow velocity and broadening accurately for all ions, we also fitted the spectra using a fixed velocity broadening to further reduce the number of free parameters. The velocity broadening was fixed to 140 km s-1 (average $\sigma $ for O VII and O VIII) and 70 km s-1 (average $\sigma $ for Si XIII and Si XIV). Those two values were chosen to study the influence of saturation and velocity broadening on the determined column density. In this fit the outflow velocity was fixed to either the derived value (see Table 4) or a value between -400 and -800 km s-1, depending on the ionization parameter of the ion. This method allows us to obtain an accurate column density for the less abundant ions. The results of this fit are listed in Appendix A as Table A.1. As in Table 4 we list the ionization parameter for which the column density of the particular ion peaks. For several ions we find a higher column density and larger velocity broadening in Table 4, which can be explained by the fact that a velocity broadening of 140 km s-1 does not include the full blend due to the five different velocity components.

  
3.2.3 Warm absorber model C: Multiple separate ionization components

In our previous models A and B we determined ionic column densities independently from the photoionization balance. A disadvantage of this method is that the column densities of the many ions with small abundances are poorly constrained, while the combined effect of these ions may still be noticeable in the spectrum. Therefore we have tried spectral fits with three xabs components.

The fit with three xabs components is not statistically acceptable. As in previous papers the iron abundance is measured to be overabundant for lower ionization ions (Kaastra et al. 2002a; Blustin et al. 2002; Steenbrugge et al. 2003; Netzer et al. 2003). A possible explanation was given by Netzer et al. (2003). Namely, the ionization balance for M-shell iron in low temperature, photoionized plasmas depends strongly on the dielectronic recombination rates. These rates are sometimes poorly known due to resonance effects at low energies, and are thus not correctly incorporated in the photoionization balance codes that are used to predict the ionic column densities. In recent articles Kraemer et al. (2003) and Netzer (2003) estimate that the ionization parameter for the iron ions that produce the UTA can change by as much as a factor of two or more. An increase of the effective ionization parameter $\xi $ by a factor $\sim$2 for Fe VII-Fe XIII indeed would bring these ions more in line with the trend observed for the other elements.

Fitting iron simultaneously with the other elements will bias the ionization parameter measured. A solution is to decouple iron from the other elements. Since however we need at least three xabs components, this would imply doubling of the number of free parameters to six xabs models. Another option is to fit instead iron separately with a slab model. The results of this fit are listed in Table 5. However, as iron has the largest ionization range and is the only element that samples both the highest ionized and lowest ionized absorber, both options severely limit the ionization range we are able to detect unbiased. The error bars derived for the fit without iron are much larger; and the ionization parameters measured should not be compared with those listed in our earlier RGS analysis, as we included iron in the fit. Forcing a large velocity broadening to fit all the outflow velocity components, and fitting the iron ions separately with a slab model, we find a decent fit for three ionization components ($\chi^2$ = 3435 for 2758 degrees of freedom). Adding more xabs components leads to highly correlated errors and fitting results. We list the results of this fit in Table 5.

Table 5: The best fit results for a model with three xabs components fitting all elements but iron, which was fit separately with a  slab model (model C). To fit the full blend we froze the velocity broadening to 200 km s-1 and the outflow velocity to -530 km s-1.

Considering the error bars, due to the exclusion of iron in the fit, the model with different ionization parameters is only constrained for two ionization components (B and C) using the full observation. In order to further reduce the number of free parameters, and thus to obtain a better constrained fit as well as to test whether the ionization parameter distribution is continuous, we next test a model with a continuous ionization structure.

  
3.2.4 Warm absorber model D: Continuous ionization structure

Steenbrugge et al. (2003) found a power law like distribution of deduced hydrogen column density versus ionization parameter. In order to test whether this apparent power law distribution can be produced indeed by an underlying power law, we have fitted the spectrum directly to such a distribution using the warm model (Sect. 3.1). Similar to model C, we get a poor fit if iron is included. Modeling all iron ions using a slab model and all remaining elements with the warm model, we obtain a good fit to the data, namely $\chi^2$ = 3305 for 2754 degrees of freedom, assuming solar abundances. Table 6 lists the best fit parameters using the warm model. In this fit the best fit outflow velocity is -530 km s-1 and the measured velocity broadening is 140 km s-1. The derived column density distribution is in excellent agreement with the results obtained from model B (see Fig. 4).


  \begin{figure}
\par\includegraphics[angle=-90,width=8.75cm,clip]{H0138F4.ps}\end{figure} Figure 4: Comparison between the ionic column densities calculated with the warm model (x-axis is model D) and the ionic column densities obtained with the slab model (y-axis is model B). The solid line represents the solution where the predicted warm and measured slab column densities are equal.
Open with DEXTER

Table 6: The best fit values for model D, using one warm component and fitting iron separately as a  slab component.

In Fig. 4 we compare the ionic column densities predicted by model D with the column densities measured by model B. The correlation between the measured and predicted column density is rather tight, with a measured slope of 1.10 $\pm$ 0.07. However, one also notes that the measured column densities using model B are higher than for model D. This is due to the fact that with method B we optimized the outflow velocity. This is not possible with model D, where the outflow velocity and velocity broadening are tied for all ions.

  
3.3 Reliability of column density estimates

In Table A.1, in Appendix A we list all the ions that we fitted with the slab model (model B). As the slab model fits both for the absorption lines and the corresponding edges, some column density determinations could be dominated by an edge measurement, rather than by absorption lines. Measurements dominated by an edge, although resulting in quite accurate column densities are highly dependent on the correct calibration of the instrument and the continuum fit, and thus potentially have a large systematic error. A further uncertainty is introduced by unknown blending at the longer wavelengths (i.e. above 60 Å), as the atomic data for lines produced at this wavelength region are sometimes poorly known. In Appendix B we discuss those ions for which the column density maybe affected by the above effects.

  
3.4 Narrow emission lines

Table 7: Narrow emission lines. The Equivalent Width ( EW) as measured with the LETGS, except for the Mg XI, Al XIII and Si XIII forbidden lines (MEG), and flux are listed. In the last column the ionization parameter, $\xi $, where the ion has the highest column density is given. Forbidden lines are indicated by f, intercombination lines by i. Below the line the C V and O VII RRC are listed. The rest wavelengths are taken from Drake (1988).

Our spectrum shows the presence of a few narrow emission lines. Table 7 lists the measured strength or upper limits of forbidden and intercombination lines of several ions as well as the Radiative Recombination Continua (RRC) observed in the spectra. For all narrow emission lines, with the exception of the O VII and Ne IX forbidden lines, the wavelength was frozen to the value in the restframe of NGC 5548. For the O VII and Ne IX forbidden lines we determined blueshifts of -0.009 $\pm$ 0.004 Å and -0.005 $\pm$ 0.007 Å, corresponding to an outflow velocity of -150 km s-1 and -175 km s-1, respectively. These blueshifts, however, become negligible if instead, the redshift of 0.01717 based upon the 21 cm line (Crenshaw & Kraemer 1999) is assumed.

The intercombination lines as well as the RRC's are undetectable or very weak. In particular the O VII intercombination line is not detected, likely due to blending by the absorption line of O VI at 21.79 Å rest wavelength. The forbidden O VII line seems double-peaked in the MEG spectrum (see Appendix C, Fig. C.1), and is poorly fit by a single Gaussian line. However, the LETGS spectrum and the earlier MEG spectrum of 2000 do not show this line profile. None of the other forbidden lines show any broadening or a double-peak profile. We thus conclude that the MEG O VII forbidden line shape is due to noise.

  
3.5 Fe K$\alpha $

The narrow Fe K$\alpha $ emission line is clearly seen in the HEG spectrum (see Appendix C, Fig. C.1). The flux of the line is (0.24 $\pm$ 0.08) ph m-2 s-1. Our results are in good agreement with the EPIC results presented by Pounds et al. (2003) and the earlier HEG data presented by Yaqoob et al. (2001) (Table 8). The Fe K$\alpha $ emission line is also discussed by Yaqoob & Padmanabhan (2004). There is no evidence for a broadened Fe K$\alpha $ emission line, consistent with the results obtained by Yaqoob et al. (2001) and Pounds et al. (2003).

Table 8: Parameters of the Fe K$\alpha $ line from the present HEG spectrum (2002), previous HEG spectrum (2000; Yaqoob et al. 2001), and EPIC spectrum (2001; Pounds et al. 2003). The FWHM is given in km s-1, the flux in ph m-2 s-1.

  
3.6 Broad emission lines

Kaastra et al. (2002a) detected a 3$\sigma $ significant relativistically broadened emission line for N VII Ly$\alpha $ and O VIII Ly$\alpha $, and a non-relativistically broadened emission line for C VI Ly$\alpha $. Ogle et al. (2004) detect a broadened C VI Ly$\alpha $ line with a FWHM of 1100 km s-1 in NGC 4051. No relativistically or non-relativistically broadened emission lines were detected in the high flux state RGS spectrum (Steenbrugge et al. 2003). The LETGS and MEG spectra show clear, broad excess emission centered on some of the deeper absorption lines, although none show the asymmetric line profile of a relativistically broadened line (Fig. 5). We fit these broad excesses with Gaussians. Due to the noise in the data, the broadness of these features, and the rather low contrast with the continuum, we fixed the width of these Gaussians to 8000 km s-1 (Arav et al. 2002) as observed for the broadest component of the C IV and Ly$\alpha $ broad emission lines in the UV. The wavelength was fixed to the rest wavelength, i.e. assuming no outflow velocity. Table 9 lists the broad emission line fluxes for the LETGS and MEG spectra. A detail of the LETGS spectrum of the O VII triplet broadened line is shown in Fig. 5.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.75cm,clip]{H0138F5.ps}\end{figure} Figure 5: Detail of the LETGS spectrum showing the fit with (thick line) and without (thin line) a broad emission line for the O VII resonance line. The profile of the broadened emission line is also plotted.
Open with DEXTER

To ascertain the significance of these broad emission lines, we checked the lines at shorter wavelengths separately in the LETGS and MEG spectra. The line fluxes are consistent within the error bars, between the MEG and LETGS observations. Inclusion of these broad emission lines in the model alters the derived column densities for the different ions in the warm absorber. Therefore the differences between the model with and without a broad emission line is in the flux of this line (see Fig. 5). As a result the column densities quoted in this paper, which are based on the fit with the broad emission lines included, should not be compared directly with column densities derived from earlier observations as small continuum mismatches are likely to remain in those spectra (e.g. Kaastra et al. 2002a; Steenbrugge et al. 2003) in a search for variability of the warm absorber. As an example, the logarithm of the derived column density for O VII (using one slab component) in a model without a broad emission line is 21.72 m-2, compared to 22.18 m-2 if the broad emission lines are included. The broad emission lines are more easily detected than in previous observations, this partly results from the low continuum flux level, and thus the larger contrast during the present observation. Further discussion on these broad emission lines will be given in a forthcoming paper (Steenbrugge et al. 2005).

Table 9: Flux of the broad emission lines, from the simultaneous LETGS, MEG and HEG fit. The wavelength was frozen to the rest wavelength of the line, while the FWHM was frozen to 8000 km s-1. The lines for which upper limits are detected were not used in the further analysis of the data.

4 Timing analysis

Kaastra et al. (2004) give a detailed analysis of the lightcurve in different energy bands. Here we focus on the possible differences in the warm absorber due to flux variations.

  
4.1 Short term variability

First we studied possible differences in the HETG and LETG spectrum, as the average flux level differed by 30% between both observations on a timescale of about 300 ks. In our first approach, we have taken the fluxed MEG spectrum and folded this through the LETGS response matrix. This is possible because of the two times higher spectral resolution of the MEG as compared to the LETGS. We adjusted the MEG continuum to match the higher LETGS continuum. Restricting our analysis to the 1.5-24 Å band, where both instruments have sufficient sensitivity, we find good agreement between both spectra.

We find no significant difference in the equivalent width of the strongest absorption lines such as the Ly$\alpha $ and Ly$\beta$ lines of Ne X and O VIII and the resonance lines of Ne IX and O VII. There may be a small enhancement of the most highly ionized lines in our spectrum (Si XIV and Mg XII Ly$\alpha $ and the Si XIII resonance line), corresponding to a $\sim$50% increase in ionic column density in the LETGS observation. For the Si XIII line this is uncertain because the line coincides with an instrumental edge in the MEG spectrum. Moreover, at the high ionization parameter where these ions are formed ( $\log\xi\sim 2.5$) the opacity of the O VIII Ly$\alpha $ line is even higher than for the Mg and Si lines. Since we do not observe an enhanced column density for O VIII, we conclude that the evidence for enhanced opacity at high ionization parameter during the LETGS observation is uncertain. A formal fit to the difference spectrum shows that any additional high ionization component emerging during the LETGS observation has a hydrogen column density less than 1024 m-2, or an order of magnitude less than the persistent outflow.

In a different approach, we have fitted the warm absorber in the LETGS spectrum between 1.5 and 24 Å using a combination of three xabs components, the results are listed in Table 10. Fitting the HETG spectra with the same ionization parameters and column densities, but fixing the continuum parameters to the best fit parameters (see Table 2) we find an excellent fit. Allowing the ionization and column densities to be free parameters, we obtain the results listed in Table 10. The improvement for allowing the ionization and column densities to be free is $\Delta\chi^2$ = 10 for 4242 degrees of freedom. The maximum difference between the two data sets occurs for the lowest ionization component 3. For the HETG spectra this component has log $\xi $ = -0.1 $\pm$ 2.5, this clearly indicates the lack of sensitivity at the longer wavelengths in the HETG spectra. Figure 6 shows the transmission for the best fit model for LETGS (solid line) and the HETGS (dotted line). Note that the largest discrepancy is around 20 Å, where the effective area of the MEG is rather low, and the lines are for the lowest ionized component 3. The flux in the soft band increased by $\sim$50% between the HETGS and LETGS observation. Therefore, an increase of $\xi $ by 50% is excluded for component 2. However, our data cannot rule out that the ionization parameter responds linearly to continuum flux enhancements for components 1 and 3. We conclude that the difference in ionization parameter of the warm absorber between the HETG and the LETGS observation is small and below the sensitivity in log $\xi $ of $\sim$0.15.

Table 10: Comparison of the best fit parameters for the LETGS and HETGS spectra between 1.5 and 24 Å. The labeling of the different components follows Kaastra et al. (2002a).


  \begin{figure}
\par\includegraphics[angle=-90,width=8.8cm,clip]{H0138F6.ps}\end{figure} Figure 6: The transmission for the best fit model with three xabs components for the LETGS (thin solid line) and the HETG (dotted line) spectra. Note that the difference is largest, 10% in continuum transmission and 35% for the largest difference in absorption lines, for lowly ionized ions which have absorption lines around 20 Å. Here the MEG effective area is already rather low, and the errors on this lowest ionization parameter are large.
Open with DEXTER

A similar situation holds for the time variability during the LETGS observations. Despite the significant flux increase during the observation (see Kaastra et al. 2004), the signal to noise ratio of our data is insufficient to rule out or confirm a significant response of the warm absorber to the change in ionizing flux.

4.2 Long term spectral variations


  \begin{figure}
\par\includegraphics[angle=-90,width=8.75cm,clip]{H0138F7.ps}\end{figure} Figure 7: Fit residuals of the 1999 LETGS spectrum minus the fit residuals of the 2002 spectrum, normalized as described in the text. For clarity, the residuals have been rebinned by a factor of 8. The dashed lines indicate the 2.5$\sigma $ significance level used in the analysis.
Open with DEXTER

We have tested for the presence of long term variations in the warm absorber of NGC 5548 by comparing our best fit model for the present data with the LETGS spectrum of December 1999 (Kaastra et al. 2000, 2002a). To do this, we have fit the 1999 spectrum keeping all parameters of the warm absorber and all fluxes of narrow or broadened emission lines fixed to the values of our 2002 spectrum. Only the parameters of the power law and modified blackbody component were allowed to vary. Then we subtracted the relative fit residuals for our 2002 spectrum from the fit residuals of the 1999 spectrum (Fig. 7). The relative fit residuals r are defined by r = (o - m)/m where o is the observed spectrum, m the best fit model spectrum. Subtracting the fit residuals instead of dividing the spectra allows us to take out continuum variations. Further, as the higher orders are not subtracted any changes in the continuum will produce broad band features in the ratio of two spectra. In Fig. 7 we normalized the difference of the fit residuals by the standard deviation of this quantity. This erases all remaining errors due to either shortcomings in the effective area calibration or in the spectral modeling.

Table 11: Features in the difference spectrum of the 1999 and 2002 LETGS data. We list Gaussian centroids $\lambda $ (in the restframe of NGC 5548), $\sigma $ and peak value p; peaks are expressed as a fraction of the 2002 spectrum at the same energy.

The figure shows clear evidence for five features with more than $2.5\sigma$ significance. There is evidence for an increasing flux toward smaller wavelengths starting around 5 Å and resulting in about a 20% excess at 2 Å in the 1999 spectrum. This corresponds to a difference in the continuum and not in the warm absorber. The remaining features are not due to continuum variations (Table 11). As the optical and UV broad emission lines are known to be variable on timescales as short as a day, we identify the broader features with changes in the broad emission lines. We identify the 34.13 Å feature as the red wing of a broad C VI Ly$\alpha $ line (restframe wavelength 33.74 Å, hence velocity of +3500 $\pm$ 700 km s-1). Similarly, the 27.18 Å feature can be identified as the red wing of C VI Ly$\gamma$ (restframe wavelength 26.99 Å, hence velocity of +2100 $\pm$ 800 km s-1). Apparently, the red wing of the C VI Ly$\alpha $ line has decreased while the red wing of Ly$\gamma$ increased and no change in the Ly$\beta$ line. However, the process that would cause the Ly$\alpha $ line to decrease while the Ly$\gamma$ line increases is unknown. The intrinsic width, $\sigma $, of the variable parts of both lines corresponds to about 900-1000 km s-1, this is much smaller than the assumed FWHM of 8000 km s-1. The 16.10 Å feature can be identified by a red wing of O VIII Ly$\beta$ (restframe wavelength 16.01 Å, hence velocity of +1700 $\pm$ 1100 km s-1). A change in the strength of the broad emission line due to the O VII triplet can also explain the dip at 21-22 Å.

Finally, the 22.36 Å feature most likely has a different origin. First, it is intrinsically narrow and poorly resolved (the upper limit to its Gaussian width corresponds to 700 km s-1). Its wavelength coincides within the error bars with the wavelength of the strongest O V absorption line (22.374 $\pm$ 0.003 Å). The 2002 spectrum shows a deep O V absorption line, which is much weaker in the 1999 spectrum. This demonstrates that the O V column is variable on the time scale of a few years. Note that the O V line is by far the strongest absorption line in the LETGS spectrum for the ionization parameter log $\xi $ = -0.2, the ionization parameter where O V reaches it maximum concentration. The Fe M-shell lines are formed at the same ionization parameter, but are much weaker.

  
5 Discussion

  
5.1 Comparison between the warm absorber models

In our model B we determined the column density, the average outflow velocity, and the velocity broadening for all observed ions separately (Table 4). The total column densities determined using this method, for the six ions C VI, O V, O VII, O VIII, Ne IX, and Si XIII are consistent with the results obtained using method A, with fixed velocity structure. This gives confidence that the column densities as listed in Tables 4 and A.1 are reliable, even if there is substructure to the lines. The average outflow velocity determined for these six ions is consistent with the outflow velocity of the deepest UV components. Comparing our measured column densities with those obtained in the UV for the five outflow velocity components, we conclude that the column density measured for the -160 km s-1 component is in fact the blend of UV components 4 and 5. Our -530 km s-1 component is a blend of UV components 2 and 3. Components 2 through 5 form one unresolvable blend in X-rays, and only component 1 is clearly separated.

  
5.2 Outflow velocity

Our method B does not resolve the full velocity structure of the outflow. Therefore one should take care in interpreting the relation between the ionization parameter and the average outflow velocity (Fig. 3). Namely, all line profiles probably are a blend of five outflow velocities, but due to limited signal to noise ratio and spectral resolution, components 3 and 4 of the UV dominate the blend.

In general, there is a trend that more highly ionized ions have a larger average outflow velocity. This could indicate that the -1040 km s-1 outflow velocity component becomes more prominent for the more highly ionized ions. The -530 km s-1 outflow velocity component dominates for lower ionized ions observed in X-rays and the UV band. To further study this effect we plot the total column density as a function of outflow velocity (Fig. 8) for the six ions of Table 3. We also added C IV and N V taken from Crenshaw et al. (2003) and O VI taken from the non-simultaneous FUSE data by Brotherton et al. (2002). For the UV data we added the two components with the lowest outflow velocity to compare it with the -160 km s-1 X-ray component. The two middle outflow velocity components were added to represent the -530 km s-1 component. Finally, for O VI we indicate the total column density as measured in the X-rays. From Fig. 8 we note that the ionic column density of the -1040 km s-1 velocity component is the smallest of the five velocity components for low $\xi $, while it becomes the largest for highly ionized gas. There is thus a clear difference in ionization structure between the -1040 km s-1 outflow component and the lower outflow velocity components. A possible explanation for this difference in ionization structure is that the outflowing wind at -1040 km s-1 is less dense, leading to a higher overall ionization of the gas, while overall the column density is smallest.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.75cm,clip]{H0138F8.ps}\end{figure} Figure 8: The ionic column density versus the logarithm of the ionization parameter for the three outflow velocity components as measured for six ions in X-rays, 2 ions (C IV and N V) measured simultaneously with HST STIS (Crenshaw et al. 2003) and the lower limits for O VI from non-simultaneous FUSE data (Brotherton et al. 2002, taking their preferred uncovered model; Arav et al. 2003). In calculating the column densities for the UV data we kept the highest outflow value separate, added the two middle ones and the two lower ones (see Sect. 5.2). For O VII the ionization parameter was decreased by 0.05 for easier identification. The open circle indicates the total column density for O VI obtained from the LETGS and MEG spectra.
Open with DEXTER

  
5.3 Velocity broadening

In Table 3 we list the measured Gaussian velocity broadening $\sigma $ of the three velocity components as well as the velocity broadening of the UV components 1, 2-4, and 5. For components 2-4 we find a slightly smaller X-ray width, while for component 5 we find a larger width as compared to the UV lines. We attribute this to the blending of (a part of) component 4 into component 5. We note that UV component 4 has $\sigma $ = 62 km s-1 compared to 68 km s-1 in component 3 and only 18 km s-1 in component 2, hence we need approximately half of component 4 to blend into component 5 in order to explain the difference.

Interestingly, we find a much smaller velocity broadening for UV component 1, the only one we can resolve. Modeling the O VIII Ly$\alpha $ line with a column density of 1022.2 m-2 and a velocity dispersion $\sigma $ of 40 km s-1, we find that the line is heavily saturated, and effectively produces a FWHM of 240 km s-1. This is similar to the 222 $\pm$ 18 km s-1 value measured from the UV spectrum. We are, however, able to obtain the velocity broadening $\sigma $ from the line ratios of the non-saturated Lyman series of O VIII and the other five ions.

In Table 4 the velocity broadening $\sigma $ is listed based upon model B. It is small compared to the velocity range of the components measured from the UV spectra. This indicates that we do not fully resolve the blend of the five outflow velocity components. The velocity broadening in our case is mainly determined from equivalent width ratios of absorption lines from the same ion. The velocity broadening for the majority of absorption lines is less than 220 km s-1. This is consistent with detecting mainly the -667 km s-1, -530 km s-1 and -336 km s-1 outflow components, which are the dominant velocity components in the UV. These three velocity components form an unresolved blend in the X-rays. Plotting the velocity broadening versus the ionization parameter produces a scatter plot.

  
6 Ionization structure

An important question about the physical nature of these outflows is whether they occur as clumps in an outflow or in a more homogeneous wind. In the first case these clumps should be in pressure equilibrium with the outflowing less dense wind, with a finite number of solutions for temperature and ionization $\xi $ for the same constant pressure ionization $\Xi $ value. In the second case one expects a continuous ionization distribution, which depends on the density of the particular part of the wind. Thus less dense material will be higher ionized, while more dense regions will be lower ionized. In this section we elaborate on both scenario's, confronting them with our observations.

  
6.1 Clumps in a wind

In several AGN outflow scenario's, instabilities in the wind may cause the formation of clumps. Such clumps can survive for longer times if they are in pressure equilibrium with the surrounding wind and are on the stable branch of the $T(\Xi)$ curve. Where the ionization parameter for constant pressure, $\Xi $, is given by $\Xi = L/(4\pi c r^2 P)$ = 0.961 $\times$ 104 $\xi/T$, with L the luminosity, P the pressure, r the distance from the ionizing source, c the speed of light and T the temperature. Ionization codes, like XSTAR (Kallman & Krolik 1999) predict a specific ionization range for which pressure equilibrium can occur (Krolik & Kriss 2001). Krongold et al. (2003) claim that there are only two ionization components in such an equilibrium in the Seyfert 1 galaxy NGC 3783. Netzer et al. (2003) needed three ionization components in pressure equilibrium to fit the same NGC 3783 spectrum. Ogle et al. (2004) find that in the case of NGC 4051 the different ionization components are not in pressure equilibrium. Another possibility is that the clumps are magnetically confined.

In order to test the presence of a finite number of ionization components, each with its own value for $\xi $ and $N_{\rm H}$, we fitted the ionic column densities obtained with model B (Table A.1) to a model with a finite number of ionization components. Since there are potential problems with the iron ionization balance (see Sect. 3.2.3), we did our analysis separately for iron and the other elements. The results of these fits are listed in Tables 12 and 13.

Using the measured column densities for all ions, except iron, we notice that the fit does not improve if more than three ionization components are fitted. The program even prefers to split up one component rather than add another in the case we fitted for five ionization components. The fit is never statistically acceptable, which is unlikely due to abundance effects as most elements have abundances consistent with solar. As an extra test we decided to fit oxygen, silicon and sulfur separately. However, due to the smaller span in ionization range and the fewer points (maximum of five), all were well fit with three ionization components.

Table 12: Fits with a finite number of ionization components to the column densities of all elements but iron derived with model B, assuming solar abundances. The number of ionization components is indicated by N. In the last column we give the significance according to an F-test of the added component.

Table 13: Fits with a finite number of ionization components to the iron column densities for the RGS data (Steenbrugge et al. 2003). The number of ionization components is indicated by N. In the last column we give the significance according to an F-test of the added component.

For iron we decided to use the column densities measured by the RGS (Steenbrugge et al. 2003), as here several more ions have a measured column density and the errors are somewhat smaller. For Fe XII and Fe XXIV, which have the same measured column density in the RGS and Chandra spectra, we used the smaller Chandra error bars. To fit these column densities adequately, one needs at least five ionization components (log $\xi $ = 3.10, 2.48, 1.56, 0.65, -0.03). In Fig. 9 we show the fit to the measured column densities with four and five ionization components. The main problem with the fit for four ionization components is Fe XXIV, which is severely under-predicted. As the RGS and Chandra spectra measure the same column density for this ion, this is a rather certain measurement.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.3cm,clip]{H0138F9.ps}\end{figure} Figure 9: The column densities for iron as measured by the RGS instrument, with the fit for a model with four (dash-dot line) and five (thick solid line) ionization components. The fit with four ionization components clearly underpredict the Fe XXIV column density.
Open with DEXTER

We conclude that we need at least five different ionization components to explain the measured ionic column densities in NGC 5548. These components span a wide range in ionization parameter (log $\xi $ between -0.03 and 3.10). Are these ionization components in pressure equilibrium? In Fig. 10 the ionization parameter if pressure is held constant ($\Xi $) versus temperature is shown for the SED used in the current analysis and the earlier RGS and LETGS analysis. The marginally stable part of the curve indicates a constant nT value, therefore components with a different $\xi $ are in pressure equilibrium on the marginally stable branch from log $\xi $ = 1.3 to 2.7. The curve has nowhere a negative slope, the smallest positive slope is between log  $\xi = 1.4{-}1.5$. It should be noted that the difference between the results produced by different codes can be large, but substituting CLOUDY 95.6 (Ferland 2002) for XSTAR we find that still not all components are in pressure equilibrium.

For the three different ionization components as observed in the earlier LETGS and RGS spectra of NGC 5548 (Kaastra et al. 2002a; Steenbrugge et al. 2003), the lowest ionization component has log $\xi $ = 0.4, which cannot be in equilibrium with the other two components at log $\xi $ = 1.98 and 2.67. If we take the five ionization components from the iron analysis, then the highest and the lowest ionization component cannot be in pressure equilibrium. This clearly suggests that the current photo-ionization balance models are either too simple, or that at least the lowest ionization component is confined by another process, possibly magnetic. As the lowest ionization component does have the same kinematics as the higher ionized absorbers, all absorbers must co-exist in a single confined outflow. The result that not all ionization components are in pressure equilibrium is rather robust, even accounting for inaccuracies in the photoionization balance codes, in particular for lowly ionized iron, or the inclusion of the blue bump.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.3cm,clip]{H0138F10.ps}\end{figure} Figure 10: The temperature versus ionization parameter for constant pressure $\Xi $ diagram for the SED used in the present analysis (thick line) and for the SED used in earlier LETGS (Kaastra et al. 2002a) and RGS (Steenbrugge et al. 2003) analysis (thin line). The $\xi $ values corresponding to the currently used SED are indicated on the x-axis on top. The five components derived from fitting iron in the RGS spectrum are indicated. The dotted lines indicate the boundaries of the marginally stable branch for the SED used in analyzing the RGS data; the dash-dotted lines are the same boundaries, but for the current used SED.
Open with DEXTER

Another possible method to distinguish between a continuous and discrete ionization is to look for spectral variability. No spectral variability of the warm absorber, with the exception of O V, was detected although the source was observed at a 2-10 keV luminosity of 4.9 $\times$ 1036 W (1999 LETGS, Kaastra et al. 2002a), 5.7 $\times$ 1036 W (2001 RGS, Steenbrugge et al. 2003) and 4.0 $\times$ 1036 W (present LETGS spectrum). However, if the ionization parameter varies linearly with luminosity, then no variation in ionization is expected to be detectable. Further complicating the possible detection of changes in the ionization parameters is the fact that in the UV the different velocity components show opposite column density variations for the same ion. The net effect will be hard to discern in the X-rays for velocity components 2-5. Component 1, which can be resolved with Chandra is the most variable component, and should be studied in the future for signs of variability. A further lack of variation in the ionization structure of component 1 possibly indicates that the ionization distribution is continuous. In the case of two or three ionization components, a luminosity change should result in a shift in the ionization parameters. In a continuous model one would expect a decrease in column density for the lowest ionization states and an increase in highly ionized material for an increase in luminosity. However, as the warm absorber is detected up to Fe XXIV, and our sensitivity for Fe XXV and Fe XXVI is insufficient, this is hard to measure with this data set. However, for the long NGC 3783 XMM-Newton observation this behaviour was indeed detected. Behar et al. (2003) studying the RGS spectra detected no spectral variability in the lowly ionized absorber, and in particular no change in the iron UTA. Reeves et al. (2004) however did detect spectral variability of the highly ionized absorber with the EPIC pn instrument.

  
6.2 Continuous ionization distribution

We decided to derive the total hydrogen column density as a function of $\xi $ for each ion, similar to Steenbrugge et al. (2003). This in contrast to the warm model where we determined the total column density at the two extrema. Contrary to the previous analysis we did not assume that the ion column densities are predominantly determined by the ionization parameter at which the relative concentration is at its maximum. In the present analysis we took into account that each ion is formed over a range of ionization parameter similar to the warm model. We implemented this by assuming a power-law distribution, as this function well describes the RGS data (Steenbrugge et al. 2003), and appears to be sufficient in the warm model:

 \begin{displaymath}%
\xi {\rm d}N_{\rm H} (\xi)/ {\rm d}\xi = A \xi^{\alpha}
\end{displaymath} (1)

with a given normalization A and index $\alpha $. We calculate the predicted ionic column densities  $N_{\rm i,p}$ for this distribution similar to the warm model. As the ionization parameter range in reality will be finite, this power law needs to have a cut-off value, which was chosen to be log $\xi $ = 4. We will use the results from model B to determine A and $\alpha $. If  $f_{\rm i} (\xi)$ is the ionic concentration relative to hydrogen at a given ionization parameter $\xi $, then  $N_{\rm i,p}$ is given by  $N_{\rm i,p} = \int_0^\infty f_{\rm i} (\xi) A \xi^{\alpha-1}$ d$\xi $. For each ion the effective ionization parameter  $\xi_{\rm i}$ is defined as

 \begin{displaymath}%
\xi_{\rm i} \equiv \frac {\int_0^\infty f_{\rm i}(\xi) \xi ...
...}{\int_0^\infty f_{\rm i}(\xi) \xi^{\alpha-1} {\rm d}\xi}\cdot
\end{displaymath} (2)

This value of $\xi_{\rm i}$ is typically 25% larger than the value of $\xi $ for which  $f_{\rm i} (\xi)$ reaches its maximum. This increase is solely due to our assumption of a power-law distribution. This can be used to derive the equivalent hydrogen column density  $N_{\rm H,i}$  for each ion (not to be confused with the integrated column density over all $\xi $ values) from the observed ionic column densities  $N_{\rm i,obs}$, since obviously we should have $N_{\rm i,obs}$ = $N_{\rm i,p}$. This yields Eq. (3):

 \begin{displaymath}%
N_{\rm H,i} \equiv \xi_{\rm i} \frac {{\rm d}N_{\rm H}(\xi_...
...\frac {N_{\rm i,obs}}{N_{\rm i,p}} A \xi_{\rm i}^{\alpha}\cdot
\end{displaymath} (3)

The calculated total hydrogen column density using this method, instead of assuming that the ion is formed at the ionization parameter where  $f_{\rm i} (\xi)$ has its maximum, is larger by at most a factor of 3. Ion concentrations were determined using a set of runs with XSTAR (Kallman & Krolik 1999) for a thin photoionized layer with solar abundances (Anders & Grevesse 1989).

The derived hydrogen column densities $N_{\rm H,i}$ versus the ionization parameter  $\xi_{\rm i}$ are presented in Fig. 11. For lower ionized ions (log $\xi $ < - 0.7) we can only determine upper limits to the column density, in contrast to the earlier RGS results (Steenbrugge et al. 2003). This results from the lower effective area of the LETGS compared to the RGS in the iron UTA region. However, due to the longer wavelength band, we do detect many previously unobservable medium ionized ions. This strengthens the relationship between the hydrogen column density and ionization parameter in the range log $\xi $ = 0 to 1. The column density increases by nearly 2 decades for a 3 decades increase in $\xi $, consistent with Steenbrugge et al. (2003).

A formal fit of the column densities of all ions except iron yields $A = 10^{24.4 \pm 0.1}$ (at log $\xi $ = 0) and $\alpha $ = 0.40 $\pm$ 0.05 with $\chi^2$ = 107 for 30 d.o.f. A similar fit to the iron ions only yields A = 10 $^{25.0 \pm 0.2}$ (at log $\xi $ = 0) and $\alpha $ = 0.2 $\pm$ 0.1, with a $\chi^2$ = 16 for 13 d.o.f.

The relationship obtained is rather tight, and the assumption of solar abundances seems justified, if we do not include those ions which may suffer from systematic bias in their column density determinations (Sect. 3.3). If the same analysis is repeated for sodium, aluminum, argon and nickel, we find that these elements are overabundant compared to solar, except for sodium. This strengthens our conclusion that these ions may suffer from strong systematic errors in their column density estimates. These ions are therefore not plotted in Fig. 11.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.3cm,clip]{H0138F11.ps}\end{figure} Figure 11: The total hydrogen column density $N_{\rm H}$, assuming solar abundances (Anders & Grevesse 1989) derived using Eq. (3), plotted versus ionization parameter. The ionic column densities were taken from Tables 4 and A.1 assuming a velocity broadening of 140 km s-1. For clarity, no upper limits have been plotted. The best fit results for model D are plotted as the two crosses connected by a dotted line.
Open with DEXTER

7 The outflow

  
7.1 Outflow geometry

In trying to get a physical picture of the outflows in AGN it is important that we know the geometry, as well as the ionization structure. By estimating the opening angle of the outflow, we can distinguish between a spherical outflow model or a model with localized streams. The model by Elvis (2000) prefers a localized stream. The model by Urry & Padovani (1995), on the other hand predicts a less collimated outflow.

Similar to Steenbrugge et al. (2003) we can obtain an upper limit to the opening angle for the outflow (Eq. (4)). We use the mass conservation formula ( $\dot{M}_{\rm loss}$ = $\Omega m_{\rm p} n_{\rm e} r_0^2 v$), with r0 the distance from the ionizing source, $n_{\rm e}$ the electron density, $m_{\rm p}$ the proton mass, v the outflow velocity and $\Omega$ the opening angle. Further we assume that the system is stationary and that the mass loss rate is equal or less than the mass accretion rate $\dot{M}_{\rm acc} = L/\eta c^2$, with L the bolometric luminosity and $\eta$ the accretion efficiency. The outflow velocities of the wind are taken from the UV measurement. This leads to the following constraint on the solid angle $\Omega$:

 \begin{displaymath}%
\Omega \leq \frac{\xi}{\eta c^2 m_{\rm p} v}
\end{displaymath} (4)

where we used:

 \begin{displaymath}%
\xi = \frac{L}{n_{\rm e} r^2}\cdot
\end{displaymath} (5)

We assume a Schwarzschild black hole, with an accretion efficiency $\eta$ = 0.057 (Thorne 1974). For a Kerr black hole, with an efficiency of 0.30 (Thorne 1974), the opening angle is reduced by a factor of 5.4. Even higher efficiencies will further reduce the opening angle. Assuming that only 50% of the mass that is accreted is lost through an outflow halves the upper limit to the opening angle.

As the ionization parameter, as well as the outflow velocities measured in this analysis closely resemble those found with the RGS, also the resulting upper limits for $\Omega$ have a similar range (but now corrected for the extra $4 \pi$). Table 14 lists the upper limits to the opening angle for the five different outflow velocity components observed in the UV and the ionization range observed in the X-rays.

We thus find that the opening angle of the outflow is very small, down to 10-4 sr for the lowest ionization component. This implies that the wind occurs in the form of narrow streamers. The range in ionization parameters at a given outflow velocity component then implies density stratification perpendicular to the outflow.

A further quantitative diagnostic of the thickness of the outflow can be derived from Fig. 11. There is a clear power-law relation between the total hydrogen column density and the ionization parameter. This power-law relation is also established by the fact that one can fit the spectra with the warm model, using only two points, the lowest and highest ionization state observed. From Fig. 11 we find a power-law slope $\alpha $ = 0.40 $\pm$ 0.05, and a normalization log A = 24.4 $\pm$ 0.1 at log $\xi $ = 0 (see Sect. 6.2). Similar to Steenbrugge et al. (2003) we combine Eqs. (1) with (5) and use d$N_{\rm H}$ = $n_{\rm e}$ dr with r the distance from the axis of the streamer, resulting in a simple differential equation for d$\xi $/dr with formal solution n/n0 = (1 + r/r0)$^{-\beta}$ where n0 is the density at the axis and r0 = $A \xi_0^{\alpha + 1}$/ $n_0 (\alpha + 1)$ and $\beta$ = 1/($\alpha $ + 1). As we do not know the distance r0 of the absorber to the nucleus of the AGN, we do not know n0; however we know r0 n0 = 1.0 $\times$ 1024 m-2 at $\xi $ = 0 and $\beta$ = 0.71 $\pm$ 0.01. Accordingly, at large distance from the axis of the streamer (low density, high ionization parameter), the density scales as $n \sim r^{-0.71}$. Steenbrugge et al. (2003) found a similar result, but due to the detection and overabundance of the iron UTA with the RGS, the power-law slope was less well determined.

7.2 Mass loss through outflow

Outflows from AGN are potentially important contributors to the enrichment of the IGM. A critical factor is the amount of mass that escapes from the host galaxy through the observed outflows. Here we estimate an upper limit to the mass loss and discuss briefly the difficulties associated with this estimate.

Table 14: The upper limit to the opening angle as calculated from Eq. (4) for the five outflow velocities measured in the UV, over the ionization range observed in the Chandra spectra. The angles are given in sr.

First the outflowing gas has to escape the potential of the massive black hole. Matter can escape if its radial velocity vr is larger than the escape velocity  $v_{\rm e}=\sqrt{2GM/r_0}$. Inserting a black hole mass of 6.8 $\times$ $10^{7}~M_{\odot}$ (Wandel 2002) we find that matter escapes only if r0 > (0.58/vr2) pc with vr expressed in units of 1000 km s-1. The measured outflow velocities are only lower limits, namely the velocity component along our line of sight, and range between -1041 km s-1 and -166 km s-1. Taking the measured outflow velocities, then r0 needed for the outflow to escape corresponds to 0.5-20 pc. The location of the wind is thus important in determining whether the gas will escape or not. Crenshaw et al. (2003) argue that the wind should cover at least a part of the inner, high ionization narrow line region (NLR), while Arav et al. (2002) argue that the narrow line region is not covered. The NLR spans at least the distance range of 1 pc (highly ionized) to 70 pc (lowly ionized) (Kraemer et al. 1998). If the argument of Crenshaw et al. (2003) holds, then in particular the highest velocity gas may be able to escape. If the vr is significantly higher than the measured v, the lower velocity components will also escape.

After escaping the immediate environment of the central black hole, the gas must escape from the bulge and halo of the galaxy. The stellar velocity dispersion of the bulge is 180 km s-1 (Ferrarese et al. 2001), and hence it is likely that most of the gas can escape the bulge provided it has not been decelerated. The effects of the disk and halo of NGC 5548 are expected to be of equal or less importance.

We neglect the possibility that ram pressure in the host galaxy terminates the outflow. Wilson & Ulvestad (1982) note the possibility that the linear radio structure observed in the host galaxy of NGC 5548 is due to material ejected from the nucleus and stopped due to ram pressure. This seems to indicate that the outflows as observed in the UV and the X-rays can escape the nucleus, whether or not the host galaxy.

If all gas escapes, then our earlier argument that the total mass loss through the wind must be smaller than the accretion rate through the disk leads to an upper limit of the mass loss of about 0.3 $M_{\odot}$ yr-1. NGC 5548 had a major interaction with another galaxy some 0.6-1.0 Gyr ago (Tyson et al. 1998); assuming that the AGN phase lasts for at least 0.6 Gyr with a steady mass loss rate, the upper limit to the mass enrichment of the IGM is 2 $\times$ $10^8~M_{\odot}$.

8 Summary and conclusions

Our long Chandra observation of NGC 5548 allowed us to obtain a unique, high signal to noise, high spectral resolution, broadband spectrum of this Seyfert 1 galaxy.

The fluxes measured for the Fe-K line and the narrow emission lines are within 1$\sigma $ consistent with earlier Chandra and XMM-Newton spectra. There is thus no evidence for long-term variability. The only narrow feature that shows column density variations is O V.

We found clear evidence for broad emission lines, similar to the optical and UV Broad Line Region lines. In particular C VI, O VII and O VIII showed significant broad lines. From a comparison of our spectrum taken in 2002 with the earlier LETGS observation of 1999, we find significant changes in the red wing (at 1700-3500 km s-1) of C VI lines and the O VIII Ly$\beta$ line. A more extensive discussion of these broad lines will be given in a forthcoming paper (Steenbrugge et al. 2005). The presence of these broad lines affects our estimates of column densities in the warm wind, by factors up to 3, a similar situation as has been found in the UV band (Arav et al. 2002).

We find evidence for the presence of an inner-shell X-ray absorption line of N V at 29.42 Å. The long exposure time and the combination of LETGS and HETGS spectra allows us to fit the warm absorber using three different spectral models. Based upon these models, we conclude that there is a good agreement between the properties of the outflow as measured through UV absorption lines and X-ray absorption lines, although in the X-ray band we do not fully resolve the spectral lines. But line centroids and derived line width are in good agreement. We find that the highest velocity component 1 at -1040 km s-1 has a different ionization structure than the other components 2-5.

We have compared our results with models for the ionization structure of a wind. We conclude that our data are not in agreement with a model with discrete clumps that are in gas pressure equilibrium with the surrounding medium. This conclusion is based upon the following arguments:

1.
the need for at least five discrete ionization components to span the broad range (at least 3 orders of magnitude) of ionization parameter $\xi $;
2.
we find, assuming our SED, a limited range in $\xi $ for which multiple solutions of constant $\Xi $ can co-exist. The lowest ionization parameter determined from the data is incompatible with this range;
3.
the fact that there are five discrete velocity components, instead of a continuous velocity range. At least the -1040 km s-1 component spans the observed range in ionization.
Our data are in agreement with a model with a continuous distribution of $N_{\rm H}$ as a function of $\xi $. Comparing the mass outflow through the wind with the accretion rate onto the black hole, we derived upper limits to the solid angle sustained by the outflowing wind. These upper limits are of the order of 10-4 sr for the lowest ionized gas in velocity component 1. We conclude that the wind occurs in the form of narrow streams with density variations perpendicular to the flow velocity. Most likely these streamers are caused by accretion disk instabilities and cross our line of sight when being radiatively accelerated. Assuming density stratification, the dependency of $N_{\rm H}$ upon $\xi $ gives the density profile across these streamers. At large distances the density is proportional to r-0.71, with r the distance from the densest part of the streamer.

Unfortunately the LETGS is not sensitive enough to detect significant variations of the warm absorber as a response to the continuum flare occurring during the LETGS observation. Reverberation studies with more sensitive X-ray instruments will allow us to determine the densities and hence the location of the wind. However regular monitoring of NGC 5548 in both X-rays and UV will allow us to derive important constraints on the long term variability of the wind. A clear example is the variability of O V between 1999 and 2002 that we found with our observations.

Acknowledgements
SRON National Institute for Space Research is supported financially by NWO, the Netherlands Organization for Scientific Research. We thank the referees for their constructive comments.

References

 

  
Online Material

Appendix A

Table A.1: The best fit column densities (in m-2) as measured using model B. The outflow velocity (in km s-1) was taken from Table 4 or a function of ionization and frozen. The velocity broadening (in km s-1) was also frozen during the fit. In the last two columns we list the ionization parameter (in 10-9 W m) for which the ion has its maximum column density. The column densities quoted are for those measured in the fit including the broad emission lines (Sect. 3.6). All ions with uncertain column densities are indicated by * (see Sect. 3.3).

Appendix B

In this Appendix we discuss those ions for which the column density may be affected by the effects described in Sect. 3.3. Carbon, Oxygen and Neon: For C IV the strongest line at 41.42 Å coincides with the instrumental C I edge, which is deep and susceptible to calibration uncertainties. For O IV the wavelength of the dominant line in the X-ray range is not well known and the absorption line is weak. For Ne VIII the deepest line at 88.10 Å is at a rather noisy part of the spectrum, the other line at 67.38 Å is blended with an Al VIII line.

Sodium and Aluminum: both Na VII and Na VIII are only detected through lines above 60 Å. The column density for Al VII is mainly determined from the 2p1/2 edge at 43.56 Å and several lines above 60 Å. A similar situation occurs for Al VIII, where the 2p1/2 and 2p3/2 edges at 51.36 Å dominate in the column density determination.

Silicon and Sulfur: the deepest lines for Si VIII and Si IX at 61.04 Å and 55.30 Å, respectively are saturated. For Si VIII this line is also blended. The saturation of these well detected lines explains the difference in column density obtained for the different velocity broadening assumed. For S X the strongest lines form a blend in the instrumental C I edge.

Argon: the column density of Ar IX is determined from an absorption line which coincides with the instrumental C I edge and the C V forbidden line. Ar X has only one detectable line at 37.45 Å, which is at a rather noisy part of the spectrum. For Ar XI the deepest line is blended with the -1040 km s-1 component of S X at 34.29 Å. Ar XII has only one weak line in the spectrum at 31.37 Å, and for Ar XIII we detect only two weak lines at 29.27 Å and 27.44 Å. The deepest Ar XV absorption line is blended with the N VII Ly$\alpha $ line. The blending is further complicated by the possible substructure in the N VII Ly$\alpha $ line due to the different outflow velocities. Ar XVI is detected from only two edges in the spectrum at 15 and 25 Å.

Calcium: Ca XIV has only one unblended line at 21.13 Å. However, the fit is driven by only one data point of the LETGS spectrum, making the column density uncertain. For all the other calcium ions the detection is very sensitive to the velocity broadening. If we leave $\sigma $ free, we mostly find upper limits to the velocity broadening of 30 km s-1, implying highly saturated lines and therefore unrealistically high column densities. The lower ionized states of calcium (Ca I to Ca VIII) have edges in the spectra, which dominate the fit.

Nickel: the detection of all nickel ions is uncertain. Nickel only produces weak absorption lines, most of them are blended with stronger iron absorption lines. Small errors in wavelength for these iron lines could be compensated in our model by absorption from nickel. For Ni I the column density is determined from the observed short wavelength tail of the edge at 105.97 Å and a smaller edge at 14.42 Å. For Ni II the 3s edge at 87.93 Å dominates the column density determination. The strongest Ni VI and Ni VII lines are blended by Fe XVI at 14.30 Å and Fe XVIII at 14.26 Å, respectively. The strongest Ni IX line is blended by an Fe XX line at 13.73 Å. The column density for Ni XI was determined mainly from the 3p-edge at 38.62 Å, close to the instrumental C I edge, and in a noisy part of the spectrum. The two strongest lines of Ni XII are blended with the Ne IX resonance line at 13.45 Å. Ni XIII causes a slight depression in the spectrum due to many weak absorption lines. However, the detection is uncertain as the continuum spectrum of the MEG around the relevant wavelength of 13.27 Å differs by about 5% from the LETGS spectrum, a similar level as the expected depression. The deepest line for Ni XVII is blended with an Fe XX at 12.90 Å; Ni XVIII has several weaker lines that are blended with Fe XX at 12.60 Å.

Due to the uncertainties in the determined column densities for the above listed ions, we decided not to use them in the further analysis.

Appendix C

Here we present the LETGS, MEG and HEG spectra of NGC 5548. The dotted lines indicate the zero flux level for the MEG and HEG spectra. The HEG spectrum was rebinned by a factor of 2 for plotting purposes, the MEG data was rebinned by a factor of 2 above 19 Å. The Fe K$\alpha $ line appears slightly shifted in wavelength between the HEG and the MEG spectra, due to calibration uncertainties or noise in the MEG spectrum. Between 4 Å and 6.2 Å there is a calibration mismatch between the HEG and the MEG spectra. The thick line is for the model with broadened emission lines, the thin line is the same model but without broadened emission lines.


  \begin{figure}
\par\includegraphics[width=6.5cm,height=8cm,angle=-90]{H0138F12.p...
...4cm}
\includegraphics[width=6.5cm,height=8cm,angle=-90]{H0138F17.ps}\end{figure} Figure C.1: Detail of the HEG, MEG and LETGS spectra. The dotted lines indicate the zero flux level for the HEG and MEG spectra.


  \begin{figure}\par\includegraphics[width=6.5cm,height=8cm,angle=-90]{H0138F18.ps...
...4cm}
\includegraphics[width=6.5cm,height=8cm,angle=-90]{H0138F23.ps}\end{figure} Figure C.1: continued.


  \begin{figure}\par\includegraphics[width=6.5cm,height=8cm,angle=-90]{H0138F24.ps...
...4cm}
\includegraphics[width=6.5cm,height=8cm,angle=-90]{H0138F29.ps}\end{figure} Figure C.1: continued.



Copyright ESO 2005