Free Access
Issue
A&A
Volume 591, July 2016
Article Number A102
Number of page(s) 9
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201527497
Published online 22 June 2016

© ESO, 2016

1. Introduction

It is generally accepted that active galactic nuclei (AGN) are powered by the accretion of matter onto supermassive black holes (BHs) situated at the centres of galaxies in the form of a geometrically thin, optically thick disc. Assuming black-body emission of the locally dissipated energy in the disc, it is expected that the optical/UV (opt/UV) spectra of AGN are dominated by a broad quasi-thermal feature with an inner disc radius temperature of the order of ~104−105 K (e.g. Shields 1978). This feature has been detected in the opt/UV part of the broad-band energy distributions of bright quasars (e.g. Elvis et al. 1994), and is commonly referred to as the AGN “big blue bump”. In addition to this emission component, X-ray emission is also a ubiquitous feature of the AGN spectra. Variability arguments and microlensing observations (e.g. Mosquera et al. 2013) indicate that the X-ray source size is less than ~10−20rg (where rg is the gravitational radius), and is probably located very close to the central BH where most of the gravitational energy is liberated. The main process for the production of the X-rays is generally thought to be the Comptonization of soft photons (like the ones emitted by the disc) by electrons of a hot (kTe ~ 100–150 keV) corona.

If the above picture is correct, the opt/UV and X-ray emission should be somehow related in AGN. For example, if the seed photons for Comptonization are those emitted by the disc, their variations should affect the cooling/heating of the X-ray corona, and hence result in X-ray spectral variations. Furthermore, given that the AGN emitted power is provided by matter accretion, one expects variations of the accretion matter rate to affect both the opt/UV and the X-ray region. If the X-ray source is located in the innermost region in these objects, the variations are expected to originate at longer wavelengths (optical) and then move to shorter ones (UV, X-rays). Finally, a fraction of the X-ray emission may be intercepted and reprocessed by the material in the disc, serving as an external heating source for the disc. In this case, the X-ray variations are expected to preceed those observed in the opt/UV band with a delay proportional to the light travel time from the X-ray source to the disc region responsible for the opt/UV emission.

An additional physical component, namely absorbing material which lies along the line of sight towards the central engine of AGN, is revealed by the detection of absorption lines in the UV and soft X-ray spectra of AGN. Approximately 50% of Seyferts and quasars show absorption signatures in the UV band (e.g. Crenshaw et al. 1999) with a similar fraction of AGN showing warm absorbing material signatures in their X-ray spectra (e.g. George et al. 1998; Piconcelli et al. 2005). Both the UV and the soft X-ray band absorbers appear to outflow. The warm absorber velocities are rather low, of the order of a few hundred km s-1, and their kinetic power is less than 1% of the AGN bolometric luminosities (Blustin et al. 2005; Krongold et al. 2007). However, over the last few years highly ionized, high-velocity outflows (with a velocity greater than 104 km s-1) have been detected in a few AGN. They are manifested as blue-shifted K-shell absorption lines from Fe XXV and Fe XXVI at energies E> 7 keV (e.g. Tombesi et al. 2013, and references therein). These ultra-fast outflows may contribute significantly to the AGN feedback to their host galaxy.

We present the results from the analysis of the Swift monitoring data of the quasar PG1211+143. This is a Type I AGN at a redshift of z = 0.0809 (Marziani et al. 1996). Its optical emission lines are narrow; the FWHM of the Hβ line is 1832 ± 81 km s-1 (Kaspi et al. 2000). This has led to the classification of this object as a Narrow Line Seyfert 1 (NLS1) galaxy. PG1211+143 shows a strong big blue bump and a bright X-ray emission. It is one of the first AGN where an ultra-fast ionized outflow with a velocity of ~0.09c has been proposed (Pounds et al. 2003). Kaspi & Behar (2006) have analysed both the CCD and the RGS data of the same 2001 XMM-Newton observation and report the presence of an outflow component of about 3000 km s-1, in contrast with the ultra-fast velocity reported by Pounds et al. (2003). Further evidence of the fast ionized outflow has been reported by Pounds & Reeves (2007, 2009) using additional XMM-Newton data. Zoghbi et al. (2015) did not detect absorption features indicative of a fast outflow in a recent NuSTAR observation. In a recent paper, Pounds (2014) has reported the results from a new analysis of the XMM-Newton RGS data, which suggests the presence of two absorbers with outflow velocities of ~0.07 and 0.14c.

Swift observed PG1211+143 several times over a period of just over two months, in March–May 2007. Bachev et al. (2009) have already studied the Swift data of the source from the monitoring campaign in 2007. They found that the source was highly variable in X-rays but showed little variations in the opt/UV bands. They also found that its X-ray spectrum can be explained by either a partial covering absorber or by X-ray reflection onto the disc, but they were not able to distinguish between the two scenarios.

We model in detail the broad-band (optical/UV/X-ray) spectral energy distribution (SED) of the source. Our main aim is to investigate whether a simple physical model, consisting of an accretion disc, a hot corona, and a warm absorber, can fit the broad-band SED of the source. The PG1211+143 data from the 2007 Swift monitoring are particularly useful for this purpose as the X-ray flux was highly variable during the campaign, while the optical/UV flux remained constant at all flux intervals. This observational fact can put strong constraints on the modelling of the broad-band SED, simultaneously, in all the different X-ray flux intervals of the source.

2. Observations and data analysis

Table 1

Swift observations log.

PG1211+143 was observed by Swift twenty times from 9 March 2007 to 20 May 2007. At first, Swift observed the source daily; after March 19, Swift observed PG1211+143 every 5 to 9 days. A summary of these observations is given in Table 1. Columns 1 and 2 list the number, the date, month, and the start time (UT) of each observation. The next four columns list the exposure time of each UVOT filter and XRT in seconds. We note that the source was not observed in the V-band during the last observation (20 May 2007).

The Swift XRT (Burrows et al. 2005) observations were performed in Photon Counting mode. The event file of the observation was created by using the Swift analysis tool xrtpipeline. Source photons were extracted from a circle with a radius of 20 pixels centred on the source. The background was selected from a nearby, source-free region with a radius of 100 pixels. Only single to quadruple events in the energy range of 0.3–10 keV were selected for further analysis. Source and background spectra were extracted from the event file by using XSELECT version 2.4b. The auxiliary response files were created by the Swift tool xrtmkarf. We used the response matrix version 012 with a grade selection 0–12.

In addition to the X-ray data, we also obtained photometry with the UV/Optical Telescope (UVOT; Roming et al. 2005) in the V,B,U,UVW1, and UVW2 filters. Swift observed the source with the UVM2 filter as well. Since its effective area overlaps the effective area of the UVW1 and UVW2 filters, the respective source’s flux will be heavily correlated with the flux measured with the UVW1 and UVW2 filters. For this reason, we do not consider the data from the UVM2 filter. Photons were extracted from a circular region with r = 5 arcsec, and the background from a circular region with a radius of 50–75 arcsec located next to the source extraction region. The UVOT tool uvotsource was used to determine the magnitudes and fluxes (Poole et al. 2008; Breeveld et al. 2010). The fluxes were corrected for Galactic reddening (E(BV) = 0.035; Schlegel et al. 1998) with the standard reddening correction curves by Cardelli et al. (1989), as described by Eq. (2) in Roming et al. (2009).

3. Optical/UV and X-ray light curves

thumbnail Fig. 1

PG1211+143 UVOT and X-ray light curves (top and bottom panels, respectively) during the 2007 Swift monitoring campaign.

The top panel in Fig. 1 shows a plot of the PG1211+143 light curves in all the UVOT bands (a value of 0.5 × 10-26 erg s-1 cm-2 Hz-1 has been added to the V filter’s fluxes for clarity). The bottom panel shows the X-ray light curves in the 0.3–2 and 2–8 keV band (soft and hard X-rays, respectively). Time on the x-axis is measured in days starting from the date of the first observation.

The UVOT band light curves exhibit variations which are within the observational errors. Application of the χ2-test confirms that the V,B,U, and UVW1 light curves are consistent with the hypothesis of a constant flux (the probability, p, is larger than 0.22 in all cases). The probability of constant flux is somewhat lower for the UVW2 light curve (p = 0.06), but this is still not low enough to claim significant variations. The UVW2 intrinsic fractional root mean square variability amplitude (Frms) is 0.02 ± 0.01 (the errors correspond to the measurement errors only, and are computed as described in Vaughan et al. 2003). The X-ray light curves are highly variable. We observe the typical, erratic variations we observe in the long term light curves of other AGN as well. The observed variations show a soft and hard band min-to-max variability amplitude of the order of ~10 and ~5.5, respectively. We found that Frms = 0.86 ± 0.02 and 0.59 ± 0.03 in the soft and hard X-ray bands, respectively.

Given that we do not detect significant variations in the optical/UV bands, it is meaningless to search for X-ray/opt/UV correlations. On the other hand, the soft and hard band X-ray variations are highly correlated. This is demonstrated in Fig. 2, where we plot the soft vs. the hard band count rate. We fitted a straight line (i.e. y = a + bx) to the data using the subroutine fitexy of Press et al. (1992) in order to take into account the errors on both count rates. A straight line fits the data reasonably well: χ2 = 25.7/18 deg of freedom (pnull = 0.11). However, the best-fit y-intercept turned out to be negative (a = −0.024 ± 0.007), which does not have an obvious physical explanation. A straight line also provides an acceptable fit when we consider the data with a hard count rate smaller than 0.04 counts/s (χ2 = 14.2 for 15 deg of freedom). In this case, the best-fit y-intercept is consistent with zero: −0.015 ± 0.007 counts/s. This is the case (good fit and best-fit y-intercept consistent with zero) irrespective of the upper limit on the 2–8 keV band count rate we consider when we fit the data, as long as it is smaller than 0.04 counts/s.

The red dashed line in Fig. 2 indicates the best-fit line to the data with 2–8 keV count rate smaller than 0.04 counts/s. Clearly, the highest X-ray flux points are not consistent with this line. They lie above the extrapolation of the best-fit line to higher count rates. Our results suggest that the relation between the soft and hard band count rates is linear, but at higher fluxes the soft band varies with a larger amplitude than the hard band. The value of the upper count rate of 0.04 c/s may seem arbitrary; however, as we discuss in Sect. 5.3, this behaviour (i.e. non-linear above a certain X-ray luminosity level) could arise in the case of changes in the opacity of a warm absorber, which would affect mainly the soft band photons.

thumbnail Fig. 2

0.3–2 keV vs. the 2–8 keV band count rate. The dashed line indicates the best-fit model to the data with a hard band count rate less than 0.04 c/s.

4. Flux resolved spectroscopy

In order to study the broad-band SED of the source, we defined three different flux intervals based on the observed 0.3–2 keV light curve: (a) a high flux-interval (HF: 3 Swift observations with a soft band count rate >0.2 ct s-1), (b) a medium flux-interval (MF: 4 Swift observations with count rates in the range 0.07–0.2 ct s-1), and (c) a low flux-interval (LF: 14 Swift observations with count rates <0.07 ct s-1). The horizontal dashed lines in the bottom panel of Fig. 1 indicate the flux limits for these intervals. The three HF observations are those which do not follow the linear trend indicated by the dashed line in Fig. 2. They were taken in the last 25 days of the monitoring campaign. The MF state data correspond to observations that were taken throughout the observing Swift run. The LF state data include most of the observations that were mainly taken during the first intensive monitoring period of the Swift campaign. We note that these intervals refer to distinct X-ray flux intervals of the source and not to its overall flux output which, given the lack of opt/UV variations, remains almost constant throughout the Swift observing run.

To create the total X-ray spectrum for each state, we used mathpha to combine the individual source and background spectra of all the observations, within the respective state. The X-ray spectra were then rebinned with grppha (v. 3.0.1) to have at least 15 photons per bin. We computed the mean flux in each UVOT band using the deredenned mean flux measurements of the individual observations within each state. We then used the flx2xsp command to transform the flux measurements to count rates, and hence to create the respective opt/UV spectrum. For the spectral analysis we used the fitting package Sherpa (Freeman et al. 2001), part of the Chandra Interactive Analysis of Observation software (CIAO; Fruscione et al. 2006).

4.1. Disc plus corona model

We first modelled the broad-band SED (extending over four orders of magnitude in energy) with a single spectral component, namely xsnthcomp in Sherpa (Zdziarski et al. 1996; Zycki et al. 1999), attenuated (only in the X-rays) by the amount of Galactic neutral absorption along the line of sight to PG 1211+143, NH = 2.7 × 1020 cm-2 (Kalberla et al. 2005). This model provides a description of the continuum shape from thermal Comptonization. The model parameters are (a) the electron temperature, Te, of the hot Comptonizing corona, which determines the high-energy cut-off of the X-ray spectrum; (b) the photon index, Γ, of the Comptonized X-ray powerlaw; (c) the redshift, z, of the emitting source; (d) the model normalization, N; and (e) the soft-photon temperature, TBB. We chose the seed photons to have a multicolour, disc black-body energy distribution in order to simultaneously fit the opt/UV and the X-ray data with the same model.

In practice, we did this by fitting the six opt/UV and X-ray spectra (two spectra for each of the HF, MF, and LF intervals) with six distinct xsnthcomp model components. For each of the three flux-state spectra, all parameters of the three X-ray xsnthcomp components are linked to the respective parameters of the associated UV xsnthcomp components. For all spectra, the high-energy cut-off of the Comptonized X-ray power law (to which our data are insensitive with the current UVOT+XRT data) is frozen to Te = 100 keV1. The only three parameters free to vary in the fit are therefore TBB, Γ, and N. Since we do not observe any significant opt/UV flux variability during the Swift campaign, we force TBB to be the same for all three flux-state spectra, while N and Γ are left free to vary independently in the different flux intervals (hereafter Model A).

Table 2

Best-fitting Model A results (χr2\hbox{$\chi^2_{\rm r}$}/d.o.f. = 2.5/218).

thumbnail Fig. 3

Broad-band spectra of PG1211+143 in its HF (black), MF (red), and LF (blue) flux; super-imposed are the corresponding best-fitting model A lines (top panel) and the residuals (in standard deviations) between the data and their best-fit models (bottom panel).

Table 2 summarizes the best-fitting parameter values and statistics, for Model A. Figure 3 shows the HF, MF, and LF spectra (black, red, and blue points, respectively). The lines above each spectrum indicate the best-fit Models A (top panel). The plot in the bottom panel shows the residuals (in standard deviations) between the data and their best-fit models. The model fits the opt/UV part of the spectra well; however, it fails to accurately model the X-ray data, yielding a global reduced χ2 (χr2\hbox{$\chi^2_{\rm r}$}) of 2.5 for 218 deg of freedom (d.o.f.). For all three flux intervals, the residuals in the X-rays show clear systematic deviations; there are two pronounced excesses between 0.3–0.7 keV and 2–6 keV, and a deficit in between (Fig. 3, bottom panel). This is a clear signature, in low-resolution CCD-like spectra, of the presence of an ionized absorber, attenuating the primary X-ray continuum mainly at energies around 1 keV where the residual opacity due to highly ionized metals is stronger. This is not a surprising result, given the numerous previous reports for the presence of such material in PG1211+143.

4.2. Adding the warm-absorber

Table 3

Best-fitting Models B and C: χr2\hbox{$\chi^2_{\rm r}$}/d.o.f. = 0.9/208 and χr2/d.o.f.=0.9/212\hbox{$\chi^2_{\rm r}/{\rm d.o.f.} = 0.9/212$}, respectively.

We therefore added a dust-free ionized absorber component (i.e. that does not affect the opt/UV data) to our Model A by multiplying the continua of the three flux-state spectra by three distinct phase components (our PHoto-ionized Absorber Spectral Engine model, Krongold et al. 2003, implemented as a user-model in Sherpa), one for each flux interval (hereafter Model B). The parameters of a phase component are (a) the ionization parameter U, defined as the ratio between the volume density of ionizing photons (i.e. with energy 13.6 eV) at the illuminated face of the ionized cloud of gas, nphot, over the electron volume density in the cloud, ne: U = [Qion/ (4πcR2)] /ne = nphot/ne (where Qion is the total number of ionizing photons); (b) the equivalent H column density of the cloud along our line of sight, NH; (c) the inflow/outflow radial velocity of the absorber along the line of sight, υr; and (d) the fraction of ionizing source covered by the absorber along the direction of the observer, fc.

We first left all three of the relevant parameters, U, NH, and fc, free to vary independently between the three phase components, while the radial velocity of the absorber was kept linked to a single value in all three spectra and initially left free to vary in the fit. As before with Model A, TBB was forced to be the same for all three flux-state spectra, while N and Γ were left free to vary independently in the different flux-state spectra.

Table 3 summarizes the best-fit parameter values and statistics for Model B, which provides an excellent description of both the UV and X-ray data in all three flux intervals. The best-fit soft-photon temperature is TBB = 1.76 ± 0.03 eV. In columns 8 and 9 we list the model ionizing flux (i.e. the unabsorbed 0.0136–100 keV continuum flux) and the model unabsorbed X-ray flux in the 2–10 keV band.

The model results indicate that neither the column density nor the covering factor varies (within the errors) between the three flux-state spectra. Moreover, while the sign of the line of sight velocity of the absorber indicates that the gas is outflowing (even at our low resolution we rule out positive velocities at 95% confidence level; see Fig. 4), our low-resolution X-ray spectra cannot constrain this velocity to better than ~70% (see Fig. 4). Our best-fit velocity is consistent with the accurate value of υr = −(3000 ± 500) km s-1 reported by Kaspi & Behar (2006), based on the high-resolution XMM-Newton-RGS spectrum of PG 1211+143, but only marginally consistent (within the 3σ level) with the vr ~ −22 000 km s-1 value reported by Pounds (2014) for the mildly relativistic outflow in this object (see Fig. 4). Our 0.3–6 keV spectra are not sensitive to the even higher velocity vr ≃ 0.14c outflow first reported by Pounds & Page (2006).

thumbnail Fig. 4

Contours at 1, 2, and 3σ confidence levels for the best-fitting redshift of the ionized absorber of Model B, which translates into the radial velocity of the absorber through the relationship υr = (zabszsys) * c, where zabs and zsys are respectively the redshift of the absorber and the systemic redshift of PG1211+143, and c is the speed of light.

We therefore linked the column densities and covering factors of the absorber to a single value in all three flux-state spectra, we froze the outflow velocity of the absorber to vr = −3000 km s-1 (Model C), and we repeated the fit. The best-fitting parameters and statistics of Model C are given in Table 3, while Fig. 5 shows the spectra with the corresponding best-fitting models (top panel) and residuals (bottom panel). Model C fits the UV and X-ray data in all three flux intervals as accurately as Model B (χr2/d.o.f.=0.9/212\hbox{$\chi^2_{\rm r}/{\rm d.o.f.} = 0.9/212$}). The best-fit soft-photon temperature is again TBB = 1.76 ± 0.03 eV. The residuals are flat over the entire energy range covered by our data (Fig. 5, bottom panel).

thumbnail Fig. 5

Same as Fig. 3, for the Model C best-fit results.

We note that both the ionization parameter of the photo-ionized outflow and FIon decrease monotonically from the HF to the LF spectrum. Figure 6 shows the measured number of ionizing photons per steradian emitted by the central source, Qion (divided by 4πc), plotted against the best-fit ionization parameter U, for the three flux-interval spectra, in log-log space. Clearly the two quantities are linearly correlated (as shown by the best-fit line in the same figure). The ionized absorber is thus consistent with being always in photo-ionization equilibrium with the ionizing flux during the entire Swift campaign. This allows us to infer an accurate estimate on the quantity (neR2). Using the definition of U, it can be shown that log (Qion/ 4πc) = log U + log (neR2). Therefore, the intercept of the best-fit line shown in Fig. 6 is a measure of the product between the electron volume density in the absorbing gas and the square of the distance of the cloud from the ionizing source. In this way, we find that neR2 = 2 × 1042 cm-1.

thumbnail Fig. 6

Best-fitting number Qion of ionizing photons at the illuminated face of the absorbing cloud of gas (divided by the quantity 4πc), vs. the best-fitting ionization parameter, U, for the three flux-interval spectra HS, MS and LS. Both the x and y-axis are in logarithmic scale. The straight line indicates the best-fitting linear relationship between log(Qion) and log (U), whose intercept provides an estimate for the quantity neR2 (see Sect. 4.2).

4.3. Density, distance, and structure of the outflow

The photo-ionization equilibrium scenario also allows us to remove the degeneracy between ne and R, by setting a lower limit on ne. If the PG1211+143 absorber is in photo-ionization equilibrium during the entire campaign, the absorbing gas must be able to relax to photo-ionization equilibrium over a time scale shorter than the shortest time interval between two consecutive different X-ray flux intervals of the quasar.

The time needed for a cloud of gas to reach equilibrium with its ionizing radiation depends inversely on the gas volume density, both during increasing and decreasing phases of the ionizing radiation. It is typically significantly longer during recombination than during ionization phases (e.g. Nicastro et al. 1999; Krongold et al. 2007). For a 3-dominant-ion condition (i.e. when the ion fractions of a given element are mostly distributed among 3 consecutive species) this equilibration time can be analytically approximated as teqXi,Xi+1(tt+dt)~[1αrec(Xi,Te)eqne][1(αrec(Xi1,Te)αrec(Xi,Te))eq+(nXi+1nXi)eq]t+dt,\begin{eqnarray*} &&t_{\rm eq}^{X^i,X^{i+1}}(t\to t+{\rm d}t) \sim\\ &&{\left[{1 \over {\alpha_{\rm rec}(X^i,T_{\rm e})_{\rm eq} n_{\rm e}}}\right]} {\left[{1 \over {\left({{\alpha_{\rm rec}(X^{i-1},T_{\rm e})} \over {\alpha_{\rm rec}(X^i,T_{\rm e})}}\right)_{\rm eq} + \left({n_{X^{i+1}} \over {n_{X^i}}} \right)_{\rm eq}}}\right]_{t+{\rm d}t}}, \end{eqnarray*}

where teqXi,Xi+1\hbox{$t_{\rm eq}^{X^i,X^{i+1}}$} is the equilibration time scale relative to the transition from the ions i and (i + 1) of the element X, as the source flux changes over the time interval from t to (t + dt); αrec(Xi,Te), ne, and nXi are the radiative-recombination coefficient of the ion Xi at an electron temperature Te, the gas electron volume density, and the relative fraction of the ion Xi, respectively. The subscript eq indicates quantities at equilibrium at the time (t + dt).

In our case, the 3-dominant-ion hypothesis is verified for several light elements. For example, our Model C best-fit results show that the three highest ionization ions of oxygen (He-like, H-like, and fully-stripped) have three-digit approximated fractions of fOVII = 0.277, fOVIII = 0.325, and fOIX = 0.398 (log U = 1.317: LS); fOVII = 0.042, fOVIII = 0.287, and fOIX = 0.671 (log U = 1.404: MS); and fOVII = 0.002, fOVIII = 0.106, and fOIX = 0.892 (log  U = 1.51: HS), and that their sum is equal to unity in all cases. During the Swift campaign, the shortest time interval for a transition from one flux-interval to a another is that between the fourth (flux-interval MS) and fifth (flux-interval LS) observations, separated by Δt = 1 day (86 400 s).

We can then evaluate the quantities in Eq. (1) for the three highest ionization ions of oxygen at the end of the transition, i.e. in the flux-state LS. The temperature of photo-ionized gas with log U(LS)= 1.317 (and log NH = 23.2, in cm-2) is T(LS)= 2.5 × 105 K. This gives αrec(OVIII;2.5 × 105) = 2.5 × 10-12 cm3 s-1, αrec(OVII;2.5 × 105) = 3.8 × 10-12 cm3 s-1 (coefficients from Shull & van Steenberg 1982), and (nOIX/nOVIII) = 1.22. Thus,teqOVIII,OIX(4th5thobservation)~1.5×1011ne-1<86400s\begin{eqnarray} t_{\rm eq}^{\rm OVIII,OIX}(4{\rm th} \to 5{\rm th} observation) \sim 1.5 \times 10^{11} n_{\rm e}^{-1} < 86\,400~s \end{eqnarray}(1)or ne> 1.7 × 106 cm-3. From (neR2) = 2 × 1042 cm-1, we then get R< 1.1 × 1018 cm = 0.35 pc (or 1.1 lt-years). This radius is fully consistent with the broad line region (BLR) radius in this object, which is of the order of ~0.3 lt-years (Kaspi et al. 2000).

Following Krongold et al. (2007), and assuming homogeneity in the flow, we can then estimate the line-of-sight thickness of the outflow, ΔR. This is simply equal to the ratio of the equivalent H column density over the H volume density, nH. For a fully ionized gas (which is a good approximation here), nH is equal to ne/ 1.23. Therefore, using the best-fit NH estimate and the upper limit on ne that we estimated above, we find that ΔR = NH/nH ≃ 1.23NH/ne< 1.1 × 1017 cm = 0.035 pc (0.11 lt-years). We can also provide an upper limit on the relative thickness of the out flow as follows: R/R) = (NH/nH) /R = 1.23NH(neR2)− 1/2(ne)− 1/2. Using our neR2 estimate (see Sect. 4.3), our best-fit NH value, and the lower limit on ne, we find that ΔR/R< 0.1.

5. Discussion and conclusions

We present the results of a detailed study of the data from the Swift monitoring campaign of PG1211+143 in 2007. This is a bright quasar, with a reverberation BH mass estimate of ~1.5 × 108M (Peterson et al. 2004). Based on our best-fit Model C results, the 1 μm−100 keV luminosity of PG1211+143 (which effectively measures the bolometric luminosity, Lbol) is of the order of 5 × 1045 erg s-1. This is ~0.4 of the Eddington limit for a 108 solar mass BH.

We studied the broad-band spectra of the source by constructing opt/UV/X-ray spectra in three different X-ray flux intervals. Our main aim was to describe the broad-band SED with the simplest physical model possible, and to use the observationally strong fact of highly variable X-ray emission in the absence of opt/UV variations to constrain the model as much as possible. We found that the three SEDs are modelled accurately by two physically motivated, spectral components. The first accounts for the effects of the warm absorber in this object. The second accounts for the opt/UV/X-ray emission. We find that a single nthcomp model, with a constant temperature of the input soft photons, can fit the opt/UV/X-ray SED well in all three X-ray flux levels. From a physical point of view, this result corresponds to the emission from an accretion disc consisting of multiple black-body components and Comptonization from a hot corona which covers the disc.

The configuration of the inner region in PG1211+143 could be similar but not identical to what has been observed in other AGN as well. For example, Petrucci et al. (2013) presented the results from the study of the broad-band SED of Mrk 509. The physical picture in PG1211+143 could be similar to the picture shown in their Fig. 10. The deep layers of the accretion disc radiate like a black body at TBB ~ 2 eV, and there is a hot layer on top of the disc (we cannot constrain its temperature with the current data; we assumed a temperature of 100 keV), which upscatters the optical-UV photons to form the X-ray spectrm. Basically, the upper disc layer in PG1211+143 appears to be much hotter than the WARM corona in Mrk 509 (as Petrucci et al. called it), and could be the HOT corona itself.

A soft-excess spectral component is not necessary to fit the present data set. In fact, evidence that this component is missing in PG1211+143 is quite strong and is independent of our spectral model fit assumptions. This evidence results from the fact that the soft and the hard band count rates vary proportionally most of the time, and the soft-vs-hard counts plot is well fitted by a straight line with an intercept consistent with zero. Noda et al. (2011, 2013) have detected significant positive offsets in similar soft-vs.-hard “count-coubt” plots in a handful of X-ray quasars. They have interpreted them as evidence of an extra soft-excess component which is less variable than the X-ray continuum. The lack of a similar non-positive offset in our case argues against the presence of an extra, less-variable, soft-component in this source.

Our results strongly suggest that both the soft and hard band X-ray photons are produced by the same component, hence the linear relation between the two. At the highest X-ray flux intervals the lineal relation breaks. The soft X-ray count rates increase by a factor that is larger than in the hard band. As we demonstrate and discuss in Sect. 5.3, this is most probably due to strong opacity changes, caused by higher ionization of the warm absorber, which affect mainly the 0.3–2 keV energy band, and not the X-ray photons at higher energies.

5.1. X-ray and opt/UV variability

We do not detect any significant optical/UV variability during the 80-day observations. The lack of opt/UV variations is consistent with the high BH mass measured in this object, as the various disc time scales are expected to be long in this case. For example, the viscous time scale for an accretion disc around a BH is given by the relation tvisc~105α0.1-1R33/2M8(r/hd)2\hbox{$t_{\rm visc}\sim 10^5\alpha_{0.1}^{-1}R_{3}^{3/2}M_8(r/h_{\rm d})^2$} s, where α0.1 is the viscosity parameter (in units of 0.1), R3 is the disc radius (in units of 3 Schwarzschild radii, RS), M8 is the BH mass (in units of 108M), and r/hd is the ratio of radius over the disc thickness (Czerny 2006). For a thin disc (i.e. a disc where hd/r ~ 0.1), the viscosity time scale is of the order of ~120 days for M8 ~ 1 at a radius equal to 3RS (i.e. the innermost radius for a stable circular orbit in a non-rotating BH) assuming a value of 0.1 for the viscosity parameter. This time scale is longer than the duration of the Swift monitoring campaign of PG1211+143.

On the other hand, the soft and hard X-ray band count rates are highly variable at all time scales. We believe that the simultaneous absence of significant opt/UV variations is not consistent with the hypothesis of inwards propagating disc accretion rate fluctuations. The absence of opt/UV variations argues against the operation of accretion rate fluctuations at larger disc radii. Therefore, the long term X-ray variations cannot be due to propagating fluctuations that happened at larger radii. Of course, the propagation time scale may be longer than the time span of the 2007 Swift observations, and the X-ray variations we observe may be due to fluctuations at larger radii which happened before the start of the observations. It is possible that we were somewhat unlucky and that disc fluctuations at large radii simply did not occur during the 2007 observations, although we consider it rather unlikely.

Furthermore, X-ray illumination does not seem to significantly affect the disc emission in PG1211+143 (as again, we would expect the opt/UV flux to be variable). This is not surprising, given the small fraction of the X-ray over the bolometric luminosity (which is mainly determined by the disc flux). The ratio between the Fion and F2−10 fluxes (listed in Table 3) is of the order of 18–36 in the HF and LF intervals (given the steep X-ray spectrum, the F2−10 flux is representative of the total X-ray flux in this object). Therefore, energetically, the X-rays are not powerful enough to influence the opt/UV emission of PG1211+143.

5.2. Constraints on the disc and X-ray corona

According to our results the temperature of the inner disc radius, TBB, remains constant and equal to ~2 eV ~2.3 × 104 K. Assuming that the opt/UV emission in PG1211+143 is produced by a multicolour accretion disc, the hottest temperature in the disc should be ~0.5(3GM/8πRin3σ)1/4\hbox{$0.5(3GM\dot{M}/8\pi R_{\rm in}^3\sigma)^{1/4}$} (Pringle 1981). Assuming a BH mass of ~108M, Rin = 3RS, an accretion rate of ~ 0.4Edd, and an efficiency of 0.1, we find that the maximum disc temperature should be of the order of ~9.5 eV, which is larger than the soft photon temperature of ~2 eV that we find. This temperature, for the same physical parameters, corresponds to a disc radius of ~25RS.

The accretion flow at smaller radii may be filled by the hot X-ray corona. The X-ray spectral slope we found (~2.5–2.6) is rather steep, but is in agreement with the measurement of Zoghbi et al. (2015). We cannot constrain the temperature of the X-ray corona in PG1211+143, but we do observe normalization and spectral slope variations. The spectral shape of Comptonization models depends on the temperature and optical depth of the corona. We detect a ΔΓ ~ 0.15 steepening of the spectral slope in the 2–10 keV band, as the flux decreases by a factor of 3, from the high to the low flux state. Interestingly, this “steeper when fainter” is uncommon in Seyferts, where we usually observer a “steeper when brighter” behaviour (i.e. Sobolewska & Papadakis 2009). This is further evidence for the decoupling between the opt/UV and the X-ray emission. A variable opt/UV soft photon input would result in spectral slope – flux relation opposite to what we observe. The nthcomp model does not provide an estimate for the optical depth of the hot plasma. We compared the nthcomp best-fit spectral models with those from the compps model (Poutanen & Svensson 1996), assuming a slab geometry and the same Te and TBB values. We found that the observed spectral slopes correspond to an optical depth of τ ~ 0.15–0.1 (for the HF and LF spectra, respectively).

We note that nthcomp does not include an explicit factor for the relative power between the disc and the corona emission. We can get a rough estimate of the ratio of the energy released in each phase following Appendix B of Petrucci et al. (2013). Using their Eq. (B.4), we find that the soft photon count rate, ns, is roughly equal to the observed photon count rate (since τ is much smaller than unity). The observed count rate can by computed by the best-fit nthcomp model (we assume here an average Γ of 2.5, an average Norm of 3 × 10-3, and kTe = 100 keV). Knowing ns and the temperature of the disc modified black-body spectrum, we can estimate the soft photon luminosity, Ls. Using Eq. (B10) of Petrucci et al. (2013), we can then estimate the heating power, Lh, liberated in the corona to Comptonize the soft photon field. In this way, we found that the amplification ratio, A = Lh/Ls, is roughly equal to 2. This is the case when all the gravitational power is liberated in the hot medium (Haardt & Maraschi 1991). However, if this were the case, we would expect the disc emission to be entirely due to X-ray reprocessing, and hence, to be variable as well, which is not the case. If the X-ray emission region is outflowing, so that the photons Comptonized in the corona are only emitted upward, then LhLobsLs and A ≈ 1, which is the case when the gravitational power is released almost entirely in the disc. Consequently, our results suggest the case of an outflowing corona in PG1211+143.

5.3. Warm-absorber and X-ray emission

The X-ray fractional variability amplitude is of the order of ~0.6–0.85, which is quite large. For comparison, the maximum observed fractional variability amplitude of nearby bright Seyferts is of the order of ~0.4–0.5 over time scales of a year or so (Zhang 2011). This is smaller than the value we observe in PG1211+143 over a period of just ~2 months. This enhanced variability amplitude is probably due to the warm absorber variations. Indeed, our best-fit Model C results predict a 2–10 keV intrinsic flux variation by a factor of ~3. This is not enough to explain the observed max-to-min variability by a factor of ~5 and 10 in the hard and soft band, respectively. This occurs because the X-ray central source is further absorbed, in situ, by a dust-free ionized outflow crossing our line of sight to the quasar and dominating the observed variations in the soft band.

According to our results, the absorber covers almost 90% of the central source. It is outflowing with a velocity of less than 2.3 × 104 km s-1 (upper 3σ limit) and has a column density of log NH ~ 23.2. It is in photo-ionizing equilibrium with the ionizing flux and is located at a distance of less than 0.35 pc from the central source. The relative thickness of the absorber, ΔR/R, is less than 0.1. Therefore, the PG1211+143 outflow, as seen along our line of sight is significantly thinner than its distance from the central ionizing source. This suggests, at least at the location where our line of sight crosses the flow, that the outflow is not propagating outwards radially but rather transversally with respect to our line-of-sight perspective.

thumbnail Fig. 7

Ratios between the best-fitting transmitted models (at high spectral resolution of ΔE = 5 eV) of the MF/HF (black), LF/HF (green), and LF/MF (magenta) spectra in the 0.3–2 keV spectral range (see text for details).

The opacity of an ionized absorber to X-ray photons depends critically on the degree of ionization of the absorber itself. During the Swift campaign, the ionized absorber varies between log U ≃ 1.3–1.5. We expect both the He- and H-like ions of the elements lighter than S, as well as the Fe unresolved transition array (UTA) and the higher ionization Fe L transitions, to play an important role in modifying the transmitted spectral shape of the continuum between E ≃ 0.3–3 keV. This is clearly shown in Fig. 7 where we plot the ratios between the best-fitting transmitted models (at high spectral resolution of ΔE = 5 eV) of the spectra MF/HF (black), LF/HF (green), and LF/MF (magenta) in the 0.3–2 keV spectral range. The straight lines are the identity ratios for the respective cases (the three ratios have been multiplied by factors of 200 (black), 40 (green), and 0.5 (magenta) to separate them and make them clearly visible in a single panel).

The changes in opacity between the three flux-state affect the entire 0.3–2 keV spectral range, and clearly extend even at energies higher than 2 keV. The strongest variations are between the HF and the LF spectra (green curves), where the high-ionization Fe L opacity changes mimic a flattening of the spectral slope at E ≳ 1.2 keV (at the low resolution of our CCD spectra), while the dramatic changes in the relative fraction of iron UTA ions and in the C-Ne K region result in the deepening of the trough between 0.4–0.8 keV during the LF. The latter is clearly also visible in the transition between the MF and the LF (magenta curves), while a smoother change in the overall spectral shape between 0.5 and 2 keV is induced by the change of opacity of the absorber between the HF and the MF intervals (black curves). These opacity changes are most probably responsible for the deviation from the linear soft vs. hard X-ray count rate relation at high X-ray fluxes (Fig. 2).

5.4. Implications for feed-back

Based on our best-fitting Model C estimates of the physical (log (NH) = 23.20), kinematic (vr = −3000 km s-1), and geometrical (thin flow with R< 0.35 pc and ΔR/R< 0.1) parameters of the ionized outflow of PG 1211+143, we can estimate the mass outflow rate of this moderate-velocity outflow. We assume the same conical disc-wind geometry that Krongold et al. (2007) adopted for the low-velocity outflow of the low-luminosity Seyfert 1 galaxy NGC 4051, where all the evidence also pointed to a transverse geometrically thin absorber. Following Krongold et al. (2007), the mass outflow rate for this geometry can be estimated as out ≃ 1.2πmpNHvrR< 3.2 × 1026 g s-1 ≃ 5 M yr-1 (assuming a vertical disc wind and an inclination angle of 30 deg). This upper limit is about 2.3 times larger than the Eddington accretion rate for PG1211+143 (which is about 2.2 M yr-1, assuming Lbol = 5 × 1045 ergs/s, and a 10% radiative efficiency). If the mass outflow rate is indeed that high, then outflows like this are probably short-lived episodes in the quasar lifetime.

Finally, for the kinetic power of the outflow we obtain Pout=1/2vr2out<1.4×1043\hbox{$P_{\rm out} = 1/2 \dot{M}_{\rm out} v_{\rm r}^2 < 1.4 \times 10^{43}$} erg s-1, i.e. only about 7% of the radiative power of the quasar. At this power, it would require about 2.3 Gyr for the outflow to deploy enough mechanical energy (1060 ergs) into the galaxy and the surrounding inter-galactic medium (IGM) to control their evolution. This is probably far too long for a plausible wind lifetime. However, it would only require about 3 Myr for such an outflow to deposit enough energy (about 1057 ergs, assuming an average ISM density of nH = 0.01 cm-3; e.g. Krongold et al. 2007 and references therein) in the 100 kpc radius, 1 kpc thick disc of its host galaxy, and heat the ISM to its evaporation temperature of 107 K. This is probably sufficient to produce a fast decline in the star formation rate of the galaxy.


1

Zoghbi et al. (2015) have recently reported a lower limit of 124 keV on Te. Our results remain unaffected even if we choose a value as large as 200 keV for Te.

Acknowledgments

This work was supported by the “AGNQUEST” project, which is implemented under the “Aristeia II” Action of the “Education and lifelong Learning” operational programme of the GSRT, Greece.

References

  1. Bachev, R., Grupe, D., Boeva, S., et al. 2009, MNRAS, 399, 750 [NASA ADS] [CrossRef] [Google Scholar]
  2. Blustin, A. J., Page, M. J, Fuerst, S. V., Branduardi-Raymont, G., & Ashton, C. E. 2005, A&A, 431, 111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Breeveld, A. A., Curran, P. A., Hoversten, E. A., et al. 2010, MNRAS, 406, 168 [Google Scholar]
  4. Burrows, D., Hill, J. E., Nousek, J. A., et al. 2005, Space Sci. Rev., 120, 165 [NASA ADS] [CrossRef] [Google Scholar]
  5. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
  6. Czerny, B. 2006, ASP Conf. Ser., 360, 265 [Google Scholar]
  7. Crenshaw, D. M., Kraemer, S. B., Boggess, A., et al. 1999, ApJ, 516, 750 [NASA ADS] [CrossRef] [Google Scholar]
  8. George, I. M., Turner, T. J., Netzer, H., et al. 1998, ApJS, 114, 73 [NASA ADS] [CrossRef] [Google Scholar]
  9. Haard, F., & Maraschi, L. 1991, ApJ, 380, L51 [NASA ADS] [CrossRef] [Google Scholar]
  10. Elvis, M., Wilkes, B., McDowell, J. C., et al. 1994, ApJS, 95, 1 [NASA ADS] [CrossRef] [Google Scholar]
  11. Freeman, P. E., Doe, S., Siemiginowska, A. 2001, SPIE Proc., 4477, 76 [NASA ADS] [CrossRef] [Google Scholar]
  12. Fruscione, A., McDowell, J. C., Allen, G. E., et al. 2006, SPIE Proc. 6270, eds. D. R. Silvia, & R. E. Doxsey, 62701V [Google Scholar]
  13. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Kaspi, S., & Behar, E. 2006, ApJ, 636, 674 [NASA ADS] [CrossRef] [Google Scholar]
  15. Kaspi, S., Smith, P. S., Netzer, H., et al. 2000, ApJ, 533, 631 [NASA ADS] [CrossRef] [Google Scholar]
  16. Krongold, Y., Nicastro, F., Brickhouse, N. S., et al. 2003, ApJ, 597, 823 [Google Scholar]
  17. Krongold, Y., Nicastro, F., Elvis, M., et al. 2007, ApJ, 659, 1022 [NASA ADS] [CrossRef] [Google Scholar]
  18. Marziani, P., Sulentic, J. W., Dultzin-Hacyan, D., Calvani, M., & Moles, M. 1996, ApJS, 104, 37 [Google Scholar]
  19. Mosquera, A. M., Kochanek, C. S., Chen, B., et al. 2013, ApJ, 769, 53 [NASA ADS] [CrossRef] [Google Scholar]
  20. Nicastro, F., Fiore, F., Perola, G. C., & Elvis, M. 1999, ApJ, 512, 184 [NASA ADS] [CrossRef] [Google Scholar]
  21. Noda, H., Makishima, K., Yamada, S., et al. 2011, 2011, PASJ, 63, 925 [Google Scholar]
  22. Noda, H., Makishima, K., Nakazawa, K., et al. 2013, PASJ, 65, 4 [NASA ADS] [Google Scholar]
  23. Peterson, B. M., Ferrarese, L., Gilbert, K. M., et al. 2004, ApJ, 613, 682 [NASA ADS] [CrossRef] [Google Scholar]
  24. Petrucci, P.-O., Paltani, S., Malzac, J., et al. 2013, A&A, 549, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Piconcelli, E., Guainazzi, M., Cappi, M., et al. 2005, A&A, 432, 835 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Poole, T. S., Breeveld, A. A., Page, M. J., et al. 2008, MNRAS, 383, 627 [NASA ADS] [CrossRef] [Google Scholar]
  27. Pounds, K. A. 2014, MNRAS, 437, 322 [NASA ADS] [CrossRef] [Google Scholar]
  28. Pounds, K. A., & Page, K. L. 2006, MNRAS, 360, 1123 [Google Scholar]
  29. Pounds, K. A., & Reeves, J. N. 2007, MNRAS, 374, 823 [NASA ADS] [CrossRef] [Google Scholar]
  30. Pounds, K. A., & Reeves, J. N. 2009, MNRAS, 397, 249 [NASA ADS] [CrossRef] [Google Scholar]
  31. Pounds, K. A., Reeves, J. N., King, A. R., et al. 2003, MNRAS, 345, 705 [NASA ADS] [CrossRef] [Google Scholar]
  32. Poutanen, J., & Svensson, R. 1996, ApJ, 470, 249 [NASA ADS] [CrossRef] [Google Scholar]
  33. Pringle, J. E. 1981, ARA&A, 19, 137 [NASA ADS] [CrossRef] [Google Scholar]
  34. Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, The Art of Scientific Computing (Cambrige, UK: Cambridge University Press) [Google Scholar]
  35. Shields, G. A. 1978, Nature, 272, 706 [NASA ADS] [CrossRef] [Google Scholar]
  36. Roming, P. W. A., Kennedy, T. E., Mason, K. O., et al. 2005, Space Sci. Rev., 120, 95 [NASA ADS] [CrossRef] [Google Scholar]
  37. Roming, P. W. A., Koch, T. S., Oates, S. R., et al. 2009, ApJ, 690, 163 [NASA ADS] [CrossRef] [Google Scholar]
  38. Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [NASA ADS] [CrossRef] [Google Scholar]
  39. Shull, J. M., & van Steenberg, M. 1982, ApJS, 48, 95 [NASA ADS] [CrossRef] [Google Scholar]
  40. Sobolewska, M. A., & Papadakis, I. E. 2009, MNRAS, 399, 1597 [NASA ADS] [CrossRef] [Google Scholar]
  41. Tombesi, F., Cappi, M., Reeves, J. N., et al. 2013, MNRAS, 430, 1102 [Google Scholar]
  42. Vaughan, S., Edelson, R., Warwick, R. S., & Uttley, P. 2003, MNRAS, 345, 1271 [NASA ADS] [CrossRef] [Google Scholar]
  43. Zdziarski, A. A, Johnson, W. N., & Magdziarz, P. 1996, MNRAS, 283, 193 [NASA ADS] [CrossRef] [Google Scholar]
  44. Zhang, Y.-H. 2011, ApJ, 726, 21 [NASA ADS] [CrossRef] [Google Scholar]
  45. Zoghbi, A., Miller, J. M., Walton, D. J., et al. 2015, ApJ, 799, 24 [Google Scholar]
  46. Zycki, P. T., Done, C., & Smith, D. A. 1999, MNRAS, 309, 561 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Swift observations log.

Table 2

Best-fitting Model A results (χr2\hbox{$\chi^2_{\rm r}$}/d.o.f. = 2.5/218).

Table 3

Best-fitting Models B and C: χr2\hbox{$\chi^2_{\rm r}$}/d.o.f. = 0.9/208 and χr2/d.o.f.=0.9/212\hbox{$\chi^2_{\rm r}/{\rm d.o.f.} = 0.9/212$}, respectively.

All Figures

thumbnail Fig. 1

PG1211+143 UVOT and X-ray light curves (top and bottom panels, respectively) during the 2007 Swift monitoring campaign.

In the text
thumbnail Fig. 2

0.3–2 keV vs. the 2–8 keV band count rate. The dashed line indicates the best-fit model to the data with a hard band count rate less than 0.04 c/s.

In the text
thumbnail Fig. 3

Broad-band spectra of PG1211+143 in its HF (black), MF (red), and LF (blue) flux; super-imposed are the corresponding best-fitting model A lines (top panel) and the residuals (in standard deviations) between the data and their best-fit models (bottom panel).

In the text
thumbnail Fig. 4

Contours at 1, 2, and 3σ confidence levels for the best-fitting redshift of the ionized absorber of Model B, which translates into the radial velocity of the absorber through the relationship υr = (zabszsys) * c, where zabs and zsys are respectively the redshift of the absorber and the systemic redshift of PG1211+143, and c is the speed of light.

In the text
thumbnail Fig. 5

Same as Fig. 3, for the Model C best-fit results.

In the text
thumbnail Fig. 6

Best-fitting number Qion of ionizing photons at the illuminated face of the absorbing cloud of gas (divided by the quantity 4πc), vs. the best-fitting ionization parameter, U, for the three flux-interval spectra HS, MS and LS. Both the x and y-axis are in logarithmic scale. The straight line indicates the best-fitting linear relationship between log(Qion) and log (U), whose intercept provides an estimate for the quantity neR2 (see Sect. 4.2).

In the text
thumbnail Fig. 7

Ratios between the best-fitting transmitted models (at high spectral resolution of ΔE = 5 eV) of the MF/HF (black), LF/HF (green), and LF/MF (magenta) spectra in the 0.3–2 keV spectral range (see text for details).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.