L. Binette1,2 - Y. Krongold2
1 - Département de Physique, de Génie Physique et d'Optique, Université Laval, Québec, QC, G1K 7P4, Canada
2 - Instituto de Astronomía, UNAM, Ap. 70-264, 04510 México, DF, México
Received 12 July 2007 / Accepted 20 September 2007
Abstract
Context. Emission lines in quasars are believed to originate from a photoionized plasma. There are, however, some emission features that appear to be collisionally excited, such as the Fe II multiplet bands. Shortward of Ly
,
there are also a few permitted lines of species from low to intermediate ionization.
Aims. Ton 34 (
)
exhibits the steepest far-UV continuum decline known (
)
shortward of 1050 Å. This object also emits unusually strong low to intermediate-excitation permitted lines shortward of the Lyman limit.
Methods. Using archive spectra of Ton 34 from HST, IUE, and Palomar, we measured the fluxes of all the lines present in the spectra and compared their relative intensities with those observed in composite quasar spectra.
Results. Our analysis reveals unusual strengths with respect to Ly
of the following low to intermediate-excitation permitted lines: O II+O III (835 Å), N III+O III (686-703 Å), and N III+N IV (765 Å). We compared the observed line spectrum with both photoionization and shock models.
Conclusions. Photoionization cannot reproduce the strengths of these far-UV lines. Shocks with
turn out to be extremely efficient emitters of these lines and are favored as an excitation mechanism.
Key words: line: identification - line: formation - atomic processes - galaxies: quasars: emission lines - galaxies: quasars: individual: Ton 34
In this work, we analyze the emission lines of an unusual quasar,
Ton 34, which is alternatively named PG 1017+280 or J1019+2745 with
redshift
.
It is severely deficient in ionizing photons,
since its spectral energy distribution ( SED) shows a remarkable
steepening of the continuum in the rest-frame far-UV, shortward of
1100 Å (Binette & Krongold 2007, hereafter BK07; Binette
et al. 2007). If the far-UV is fitted by a power law (
), the index
is as
steep as
.
BK07 suggest that the extreme-UV flux might
undergo a recovery shortward of 450 Å.
While the near-UV emission-line spectrum appears to be "normal'', the far-UV spectrum shows low to intermediate ionization species with unusual strengths. Using the UV SED constructed by BK07 from archive data, we will quantify this statement and present photoionization and shock models for comparison. The aim is to understand how the extreme deficiency of ionizing photons in Ton 34 might be impacting the emission-line spectrum.
The emission-line spectrum of quasar and Seyfert I galaxies is generally believed to originate from gas photoionized by a nuclear UV source. State-of-the-art photoionization models of the broad emission line region (BELR), such as those developed by Baldwin et al. (1995) and dubbed "locally optimally emitting clouds'' (LOC) models, can successfully reproduce most of the emission lines observed in quasars. A grid of these models can be found in Korista et al. (1997, hereafter KO97) and more recently in Casebeer et al. (2006 and references therein). There are, however, a few exceptions to the success of pure photoionization. In particular, photoionization models require micro-turbulences in order to reproduce the shape and intensity of the Fe II UV-band (Baldwin et al. 2004). A possible alternative is that the region producing Fe II is collisionally ionized, as proposed by Grandi (1981, 1982), Joly (1987), Véron-Cetty et al. (2004, 2006), and Joly et al. (2007). In this work, we present evidence that photoionization might not be sustainable in the case of some of the far-UV permitted lines reported in this paper.
Below we summarize the procedure used by BK07 to derive the UV SED of Ton 34.
The current work is based on four archival or bibliographical sources. The 760-1120 Å spectral segment is provided by the dataset Y2IE0A0AT from the HST-FOS archives (grating G270H). To cover the extreme UV region, we borrowed from the IUE archives. The long wavelength segment (LWP) is from Tripp, Bechtold & Green (1994) and corresponds to the dataset LW0P5708. Fluxes longward of 3000 Å (observer-frame) were severely affected by reflected sunlight or moonlight (Lanzetta et al. 1993) and have been discarded. The shorter wavelength IUE segment (SWP) was extracted directly from the archives and corresponds to the dataset SWP28188. To cover the SED behavior longward of the HST segment, we adopted the published optical spectra of Sargent et al. (1988), which were taken at the Palomar 5.08 m Hale Telescope. Both optical spectra lacked absolute flux calibration, although the authors observed standard stars, which allowed them to provide a relative calibration.
We statistically corrected the UV spectral segments for the
cumulated absorption caused by unresolved Ly
forest lines, which
are responsible for the so-called far-UV "Lyman valley'' (Møller &
Jakobsen 1990). For that purpose, we adopted the
mean
transmission function
for
published by Zheng et al. (1997). We also applied a
Galactic reddening correction assuming the Cardelli et al. (1989) extinction curve corresponding to
RV=3.1 and
EB-V = 0.13. The latter value corresponds to
the mean extinction inferred from the 100
maps of Schlegel
et al. (1998) near Ton 34. The blue and red arm segments have been
scaled to overlap smoothly with the HST-FOS segment. Both the LWP
and SWP segments were multiplied by a factor 0.75. This scaling was
necessary so that the LWP segment superimposes the HST-FOS spectrum
as closely as possible. Continuum variability is a possible
explanation for this continuum difference, since the IUE and HST
observations were made in different years. Finally, all the spectral
segments were shifted to rest-frame wavelengths, and
was
multiplied by
.
The IUE spectra have been re-binned by
grouping n pixels together (SWP with n=5 and LWP with n=3) to
improve the limited S/N. The LWP and HST-FOS spectra overlap
significantly in spectral coverage. Both datasets taken nine year
apart confirm the unusual steepness of the UV break in Ton 34.
Shortward of 1100 Å, the continuum of Ton 34 undergoes a sharp
fall off (see Fig. 2 in BK07), which BK07 model as dust absorption
by nanodiamond grains. This resulted in a deep and broad absorption
trough that fits the observed continuum reasonably well. In our
photoionization calculations presented below in
Sect. 3.2.1, we experiment with two ionizing SEDs. The
first is the intrinsic "unabsorbed'' SED, which is assumed to be a
power law of index +0.1 followed by a roll-over centered on
640 Å that extends up to the X-ray domain. Beyond 2 keV,
SED II behaves as a power law of index -1.0, yielding an
of -1.45. This SED is shown in Fig. 1 and, as in
the work of BK07, it is labeled Model II. The second SED used
in photoionization calculations is the dust-absorbed version of the
same SED, which fits the observed UV continuum of Ton 34 between 400 and 1550 Å (labeled Model IV in
Fig. 1). Shortward of 200 Å and longward of
2000 Å, the two distributions are the same. This is because
nanodiamond dust absorbs radiation over a relatively narrow domain
as compared to other grain compositions.
In Fig. 2, we present the continuum subtracted spectrum
of Ton 34, that is, the residual between the observed Ton 34 SED and
our continuum fit represented by Model IV.
![]() |
Figure 1:
Log-log plot of the input spectral energy distributions
used in our photoionization calculations discussed in
Sect. 3.2.1. These ionizing SED s are labeled II and
IV in either |
| Open with DEXTER | |
![]() |
Figure 2: Residuals of the spectral energy distribution of Ton 34 after subtracting our absorbed continuum Model IV from BK07. The different spectral segments have been color-coded as follows, SWP: red, LWP orange, HST-FOS: blue, and Palomar: dark green. Color-coded fiducial marks indicate the position of observed or expected (labeled with symbol "?'') emission lines. Measurements of line intensities and upper limits are given in Table 1. |
| Open with DEXTER | |
The procedure for measuring the flux of the lines was the following:
we first fit a Gaussian to each observed line in the spectra. For
several lines, a narrow component was required, so we added a second
(narrow) Gaussian. In addition, the lines by C IV
,
Si IV
,
and
Ly
show a clear asymmetry in the line profile, with a blue
shoulder (see Fig. 2). For these lines, we included a
third, broader Gaussian. The FWHM of the broad component spans from
3600 to 5300
.
It is interesting to note that the
O II+O III complex at around 835 Å has a significant and strong
red shoulder extending up to
850 Å, which is observed in
both the IUE-LWP and HST-FOS spectra (see Fig. 2). We
could not find any positive identification of this shoulder with any
line from a different ion/transition, so we considered this feature
as part of the O II+O III emission.
The measured line fluxes extracted from Fig. 2, as well
as upper limits of other permitted lines, are listed with respect to
in Col. 5 of Table 1. Note that we give
the total flux under the profile, that is, the integrated flux from
all the Gaussian components required to fit each emission line. A
consistency check was carried out, which showed that the line fluxes
measured over the original spectra or the continuum subtracted
spectra were indistinguishable from each other.
In Col. 5 of Table 1, we show our error estimates,
which we evaluated at a 1
significance level. We assumed an
S/N of 25 for most lines, except for N III+N IV and N III+O III,
where we assume a, S/N of
10. The line upper limits in
Table 1 correspond to a significance of 2
.
As
for the continuum, we estimate the errors to be
10%.
Of all the emission features that we measure in the far-UV, three line systems stand out by their strengths with respect to the composite spectra: the O II+O III lines at 835 Å, the N III+O III lines at 686-703 Å, and the N III+N IV lines at 765 Å.
Many weaker features in the IUE spectrum appear to lie where other
permitted lines of comparable excitation might be expected, such as
O III
508 Å, O IV
554 Å, O V
630 Å, and O IV
609 Å. A few of these have been reported before in other
quasars (Reimers et al. 1998; Laor et al. 1995) or in composite
AGN spectra (Zheng et al. 1997; Telfer et al. 2002; Scott et al. 2004). However, these line systems appear too narrow in the IUE
spectra compared to typical BELR line profiles (see the profile
comparison of Fig. 3). They lack a broad component at
their base. Given the limited S/N of the IUE spectrum at the far-UV
end, we consider it probable that these lines are spurious features
instead. For this reason, we consider these emission-like features
as upper limits rather than real detections. The symbol "?'' denotes
these unconfirmed lines in our various figures.
![]() |
Figure 3: Emission lines extracted from the Ton 34 spectrum plotted in velocity space. The flux scale is arbitrary for each inset. Left panels: near-UV permitted lines, right panel: far-UV permitted lines. Overall, the lines are all consistent with the rest-frame system of Ton 34. Differences in the position of the lines, on the right panel, may be due to absorption by intergalactic gas. The narrow line of O III at 508 Å ( bottom right panel) is severely affected by intergalactic absorption, and better data would be required to confirm its presence. The same applies to the other lines shown as upper limits in Table 1. |
| Open with DEXTER | |
We find little evidence of the high excitation Ne VIII line at
775 Å reported by Telfer et al. (2002) and Scott et al. (2004)
in their respective composite spectrum, and we favor the
identification of O IV
789 Å instead. Because the line
spectrum of Ton 34 has unusually low excitation as shown below in
Sect. 3.1, we do not believe that the high excitation
lines of Mg X and Ne VIII (listed in Table 1) are
present at a detectable level.
As can be gathered from Fig. 2, the strongest emission
features in the far-UV coincide with the position of lines observed
or expected in quasar spectra (Sect. 3.1). However, the
limited quality of the data and the possible coincidence of
absorbers at inconvenient spectral positions prevent us from
deriving incontrovertible conclusions. In the case of the narrower
features (O III
508 Å, O IV
554 Å, O V
630 Å, and
O IV
609 Å), better quality data is required to confirm or
discard their presence, as discussed in Sect. 2.4.
Clearly, new observations are needed in all wave bands down to the
X-rays. In what follows, we take the data at face value and present
photoionization and shock models that attempt to reproduce the
far-UV lines.
Table 1: Comparison of Ton 34 with composite SED s and with models.
We now quantify to what degree the emission lines differ in Ton 34 from the "average'' quasar. To achieve this, we list the line ratios
characterizing the radio-loud (Col. 3) and radio-quiet (Col. 4)
composite spectra of Telfer et al. (2002) in Table 1.
Comparison between Ton 34 and these two sets of ratios requires some
caution, since significant line ratio variations exist among
quasars. For instance, Telfer et al. (2002) report that the RMS
deviation of line fluxes between the different quasars amounts to as
much as 50-70% for the strong lines of C IV
,
O VI
,
and Ly
.
Hence, intrinsic differences of less than a factor two between the
composites and Ton 34 should not be considered significant.
![]() |
Figure 4:
Line flux ratios renormalized to
|
| Open with DEXTER | |
To facilitate the comparison of Ton 34 with the two composites, we
plot their line ratios in Fig. 4. Inspection of the
Table 1 or Fig. 4 reveals that the
commonly strong BELR lines of C IV, N V, and O VI are all present
in Ton 34. As a result, the apparent sharp turndown of the ionizing UV
in the range 650-912 Å is not radically affecting the high
excitation emission lines. In particular, the O VI
line is quite
strong, although not as much as in the two composites. The C IV is
substantially weaker, by more than a factor of six in Ton 34 with
respect to the radio-quiet composite. Also, the line system
C III+N III near 980 Å is noticeably weaker, although the flux
in this line is difficult to measure accurately due to the
uncertainties introduced by the sharp continuum bent and the many
Ly
forest lines.
In the far-UV, we note that the intensity of the O II+O III and
N III+O III systems in Ton 34 are a factor of
14 and 18 brighter, respectively, than in the RLQ composite. There is also
evidence of significant emission of N III and/or N IV at 764 and
765 Å, which are not detected in the composite spectra
either.
The line spectrum of Ton 34 show peculiarities that deserve further
analysis, in particular, O II+O III (835 Å), N III+O III lines
(686-703 Å), and N III+N IV (765 Å), which are measured
with unusual strengths with respect to Ly
.
Are these emission
features necessarily genuine lines? One possibility is that
extinction resonances, unaccounted for in the extinction curve used
to model the deep continuum trough (BK07), may induce features that
looked like broad emission lines. Another possibility is that Ly
absorbers at intervening redshifts might generate spurious emission
features by bracketing narrow continuum regions. Although we cannot
rule out either possibility with the current data, both appear
unlikely to us, on the grounds that the strongest emission features
coincide quite well with the position of plausible atomic
transitions (see Fig. 2). The two strongest line
systems of O II+O III (835 Å) and N III+O III (686-703 Å)
have previously been reported in the RLQ composite, although at a
much reduced flux level. We thus pursue our analysis under the
assumption that the observed features are real and consist of low to
intermediate-excitation permitted lines.
Can photoionization account for the strength of the far-UV permitted lines? We first establish a comparison with published BELR models and then evaluate the impact of a strongly absorbed ionizing continuum.
Baldwin et al. (1995) show that by integrating line fluxes over a
wide range in gas density
and impinging ionizing flux
,
one obtains a much improved fit to quasar line
spectra. Such models were dubbed "locally optimally emitting
clouds'' (LOC). Baldwin et al. (1995) also show that by
preferentially selecting the optimal slab density and impinging flux
for each individual line, one can derive a line spectrum comparable
(within a factor two) to that of a true LOC model. To derive an
approximate LOC model, we proceeded as follows. From the grid of
photoionization models published by Korista et al. (1997; hereafter KO97), we extracted the highest equivalent width found within the
plane
vs.
,
for each line of interest. The
particular grid that we selected was labeled AGN4
. It
assumes solar abundances and an SED that was defined by KO97,
which peaks at 22 eV. It is the closest to our SED II with a
18.5 eV turnover (Fig. 1; see also Haro-Corzo et al.
2007).
The line ratios from this approximated LOC model are shown in
Col. 6 of Table 1. Unfortunately, the N III+O III line system (
686-703 Å) was not part of the AGN4 grid,
nor was the O II
834 Å line. On the other hand, the
N III+N IV system at 765 Å and the O III line at 835 Å were. The N III+N IV system is significantly weaker than observed,
while the
835 Å O III line is predicted an order of magnitude
weaker than the observed O II+O III system. As we consider
unlikely that the
834 Å O II line (absent from the AGN4
grid) is stronger than O III, we conclude that photoionization would
have difficulty in fitting this system. Hence, even locally
optimally emitting clouds would not be able to account for the
intensities of at least some of the far-UV lines observed in Ton 34.
Could the peculiar shape of the Ton 34 SED be responsible for the
unusual strengths of some far-UV lines? Out of curiosity, we
calculated photoionization models with the multipurpose code MAPPINGS Ic (Ferruit et al. 1997; Binette et al. 1989), using SED II to
compare with the absorbed SED IV, characterized by the deep
trough. We assumed solar metallicities (Anders & Grevesse 1989) and
a gas density of
.
The ionization
parameter
was varied until a maximum in the O III]/H
(
1663/
4861) ratio was found, which occurred at
.
The models were truncated at a depth where H is 10%
ionized. These calculations with
using either SED II or IV (both plotted in Fig. 1) are reported in
Cols. 7 and 8 of Table 1, respectively. Because there
are fewer soft ionizing photons in SED IV, we find that the mean
energy of the photoelectrons is twice as high as the one given by
SED II. This must result in a hotter plasma and therefore in
stronger collisionally excited lines. A comparison of the calculated
ratios between the two models and with Ton 34 (Col. 5) reveals that,
although many metal lines in Col. 8 ( SED IV) are often stronger
than in Col. 7 ( SED II), the deep UV trough does not result in a
sufficient increase in the strengths of either the O III+N III lines at 683, 703 Å or of the O II+O III lines at 835 Å. In
conclusion, photoionization predicts far-UV line intensities that
are much too weak in comparison with our measurements. Furthermore,
making drastic changes in the shape of the ionizing continuum does
not alter this conclusion.
In view of the difficulties producing strong permitted lines of O II, O III, and N III in the case of pure photoionization, we are lead to consider whether collisional ionization might not be more appropriate.
To investigate this possibility, we used MAPPINGS Ic to calculate a
sequence of steady-state plane-parallel shock models with a preshock
density of
,
again assuming solar
metallicities. The postshock temperatures of the different models
covered the range
-
K,
corresponding to shock velocities of 75 to 235
.
The
pre-ionization state of the shocked gas was determined
self-consistently by an iterative scheme, using the ionizing
radiation produced within the cooling shock that propagates upstream
(Dopita et al. 1984). The time evolution of the electron
and ion temperatures was followed separately until they equalized,
making use of the equilibration timescale as defined by Spitzer
(1962). Most of the far-UV resonance lines are emitted downstream
in layers of densities in the range 1010.6-1011.3
,
well below the densities of 1016 where collisional
de-excitation would become a concern for many resonance lines. The
elapsed time for the shocked gas to cool to temperatures of 8500 K
is about 10 s. The adiabatic cooling and recombination of the
plasma was followed in time until the ionized fraction reached
2%. Because the integrated columns of the different ions are
modest in shocks, line opacities turn out to be negligible compared
to those of photoionized slabs. For instance, the line-center
opacity of C III
and C IV
are 20 and 1, respectively, for a
100
shock, compared to 105.3 and 104.9 for the
photoionization model of Col. 8.
![]() |
Figure 5:
Line intensities from high-density cooling shocks
renormalized to
|
| Open with DEXTER | |
The intensities of representative far-UV lines are shown in
Fig. 5 as a function of shock velocity. The
calculations show that shocks with gas densities appropriate to the
BELR are very efficient in producing strong lines of O II+O III (
835 Å) and of N III+O III (
686-703 Å), which reach
71% and 29% of Ly
,
respectively. We also computed the
intensities of many other far-UV lines that might be observable in
future observations. Some high-excitation lines such as O IV
554 Å, O IV
789 Å, and O V
630 Å, become intense
for shock velocities exceeding 120
.
By comparing the observed
upper limits for these lines in Table 1 with the
computed intensities of O III
835 Å or O III
703 Å, we
find that velocities on the order of 90-130
produce line
intensities compatible with the estimated line ratios
. To be definite, we adopted the velocity of 100
for the case model
presented in Col. 9 of Table 1.
Shock models by themselves predict far-UV line intensities that are
too strong with respect to Ly
(compare Cols. 9 and 5), creating
a reverse situation to that of photoionization
(Sect. 3.2.1). We are therefore lead to propose a mixed
model, in which we ascribe only a fraction of the luminosity of
Ly
to shock excitation and the complementary fraction to
photoionization. In this mixed model, photoionization would be
responsible for the emission of the strong near-UV (i.e. classical)
lines, while shocks would be contributing about a third of Ly
and (proportionally) all of the far-UV resonance lines shortward of the Lyman limit.
The preshock density
may be significantly higher than assumed
above. We find similar line ratios for preshock densities up to 100
times higher. The luminosity per unit area of the shock model
in this case exceeds that of the photoionization models presented
in Cols. 7 and 8 (see footnote a in Table 1).
Our code includes three-body recombination of H, but not the process
of stimulated emission, which prevents us from going beyond a
preshock density of 1011.6
.
Beyond this limit, we expect
Ly
to be the first line to thermalize, which would further
enhance the strengths of the metal lines with respect to Ly
.
In summary, the far-UV lines observed in Ton 34 shortward of the
Lyman limit are characterized by a much lower excitation energy
than the near-UV lines. For this reason, collisional excitation
(through shocks) at temperatures significantly higher than
typically provided by photoionization is strongly favored.
Calculations with MAPPINGS Ic show that such a temperature regime is
ensured when shock excitation of moderate
takes place. These
shocks would not only account for the far-UV lines, but may also
contribute significantly to the Fe II multiplet lines that have
been proposed as resulting from mechanical heating by Joly et al.
(2007, and references therein).
Acknowledgements
This work was supported by the CONACyT grants J-50296 and J-49594, and the UNAM PAPIIT grant IN118905. Diethild Starkmeth helped us with proofreading.