A&A 463, 445-454 (2007)
DOI: 10.1051/0004-6361:20054194
X. Mazzalay1 -
A. Rodríguez-Ardila2,
1 - IATE, Observatorio Astronómico Córdoba, Laprida 854, X5000BGR,
Córdoba, Argentina
2 -
Laboratório Nacional de Astrofísica, Rua dos Estados Unidos 154, Itajubá, MG, Brazil
Received 13 September 2005 / Accepted 22 August 2006
Abstract
Aims. Mrk 1210 is an outstanding Seyfert 2 galaxy because it displays signatures of recent circumnuclear star formation and a high level of X-ray activity, in addition to the classical spectral characteristics typical of an AGN. Here we investigate the extinction affecting the nuclear and extended emitting gas, the kinematics of the narrow-line region, and the physical properties and conditions of that gas.
Methods. Near-infrared and optical spectra of the nuclear and extended emission region of Mrk 1210 are presented, covering the interval 0.4-2.4 m. Emission and absorption lines were used to infer, respectively, the geometrical extension of the ionized gas and the contribution of the underlying stellar population to the observed integrated continuum. The emission line profiles were employed to study the kinematics in the NLR. The reddening and physical condition of the gas were investigated by means of flux ratios among permitted and forbidden lines.
Results. The NIR nuclear spectrum is dominated by H I and He I recombination lines, as well as [S II], [S III], and [Fe II] forbidden lines. Coronal lines of [S VIII], [S IX], [Si VI], [Si X], and [Ca VIII], in addition to molecular H2 lines, were also detected. The
1.618
m overtone bandhead helped to estimate the contribution of the stellar population to the continuum. It was found that
of the H-band continuum has a stellar origin. It improves previous estimates, which claimed that at least 50% of the observed continuum was attributed to the AGN. Analysis of the emission line profiles, both allowed and forbidden, shows a narrower (
km s-1) line on top of a broad (
km s-1) blue-shifted component. This seems to be associated to a nuclear outflow. This hypothesis is supported by 6 cm VLBI observations, which show a radio ejecta extending up to
30 pc from the nucleus. This result does not require the presence of the hidden BLR claimed to be present in previous NIR observations of this object. Internal extinction, calculated by means of several indicators including [Fe II] flux ratios not previously used before in AGNs, reveals a dusty AGN, while the extended regions are barely affected by dust, if at all. The density and temperature are calculated for the NLR using optical and NIR lines as diagnostic ratios. The results show electronic temperatures from 10 000 K up to 40 000 K and densities between
103-105 cm-3. The higher temperatures show that shocks, most probably related to the radio outflow, must contribute to the line emission.
Key words: galaxies: active - galaxies: individual: Mrk 1210 - galaxies: Seyfert - infrared: galaxies
Additional evidence of the dusty nature of Mrk 1210 is abundant. It is a member
of the sample of galaxies called 60 m
peakers (60PKs, Heisler & De Roberties 1999) because of its "warm''
far-infrared color and the spectral energy distribution peaking
near 60
m. Over the past few years, evidence has been
accumulating indicating that these properties can be attributed to a dust-obscured
active galactic nucleus (e.g., Keel et al. 1994; Hes et al. 1995;
Heisler & Vader 1995). Whether or not this obscured
material is the putative dusty torus, as postulated by the
unified model (e.g., Antonucci 1993), is highly debatable
although evidence from spectropolarimetry points to this scenario.
In fact, Heisler et al. (1997) showed that the Seyfert 2
galaxies for which polarized broad lines have been detected have
warmer far-infrared colors (
f60/f25 < 4) compared to
Seyfert 2 galaxies for which polarized broad lines have not been detected.
The warm IRAS Seyferts are believed to be viewed at an angle that is more
face-on than that of cool IRAS Seyferts. Moreover, in the near-infrared, the
degree of polarization of Mrk 1210 rises toward longer wavelengths,
reaching
5% at 2.15
m (Watanabe et al. 2003). Observational
evidence suggests that this polarization is due to dichroic absorption by
aligned grains. Watanabe et al. (2003) also suggest the existence of
an additional unpolarized component, most
probably associated to electron scattering.
The picture of a highly reddened nucleus derived from polarimetry is also strongly supported by optical multicolor imaging. Indeed, Heisler & Vader (1994) showed that Mrk 1210 displays a smooth global light distribution but morphologically resembles a typical elliptical galaxy. Masked frames display arcs or shells surrounding the main body of the galaxy. More recently, Martini et al. (2003) classified Mrk 1210 as a tightly wound nuclear dust spiral with the individual dust lanes traced for over a full rotation about the nucleus.
In X-rays, Mrk 1210 is peculiar. It is one of the very few cases of an AGN caught in a dramatic transition between a Compton-thick, reflection dominated state, and a Compton-thin state (Guainazzi et al. 2002). This result was possible after comparing XMM Newton observations with ASCA observations made six years earlier. The transition is attributed to either a change in the properties of the absorbers or a "switching-off'' of the nucleus during the ASCA observations. In the former scenario, it may indicate a clumpy or patchy structure of the torus to which the absorbers are associated, while the latter may provide clues to the duty-cycle of the AGNs phenomenon.
All the above means that Mrk 1210 is a key source
where the interplay between an active nucleus and the
presence of strong star formation can be studied in detail. Therefore, it
is a very important object for investigating the nuclear structure
of Seyfert 2 galaxies. In order to achieve this goal, it is necessary to
characterize these two components (i.e., nuclear and stellar)
as completely as possible. The
presence of broad-emission line components, a high degree of polarization,
power law near-infrared colors, warm far-infrared spectrum, high
radio-brightness temperature in the milli-arcsecond central structure,
and finally, hard X-ray emission, all suggest that a significant
contribution comes from the active nucleus. This paper is the first attempt to study
the properties in Mrk 1210 that can be attributed to the AGN,
although not leaving aside the stellar contribution.
Because of the strong evidence of dust obscuration, the NIR region
is most appropriate
because it attenuates only a fraction (1/7)
of the light compared to that in the optical, allowing an
unbiased estimate of the obscured AGN to be obtained.
Here, we study Mrk 1210 by means of JHK spectroscopy
covering the interval 0.8-2.4 m, complemented by
optical spectroscopy aimed at investigating (i) the
extinction affecting the nuclear
and extended emitting gas; (ii) the kinematics of the narrow
line region (NLR); and (iii)
the physical properties and conditions of that gas. To meet the first goal,
a variety of extinction diagnostic ratios based on permitted and
forbidden lines, some of them never used before in an AGN, are used.
This will allow us to make a detailed map of the dust distribution on that object.
For the second goal, a detailed comparison of the emission line profiles,
both in the optical and NIR, are carried out. It will also allow us to
state whether or not a hidden broad-line region is in fact present in
the spectra of this object. For the third goal, density and
temperature-sensitive line ratios are employed to construct a detailed
picture of the ionization structure of emitting gas. In Sect. 2 we
describe the
observations and data reduction process. Section 3 describes the most
important emission features found in the NIR spectra, some of them reported
for the first time in this work. Section 4 discusses the kinematics
of the NLR gas, while Sect. 5 evaluates and
discusses the dust distribution in the nuclear and extended
emission regions of Mrk 1210. Section 6 explores the physical condition
throughout the narrow-line region, and concluding remarks
are found in Sect. 7.
A Hubble constant of 75 km s-1 Mpc-1 is used in this work.
Telluric absorption correction and flux calibration were applied to the individual 1-D spectra by means of the IDL based routine xtellcor (Vacca et al. 2003). Xtellcor makes use of the spectrum of the A0V observed after the source (see above) and a high resolution model of Vega, to construct a telluric correction spectrum that is free of telluric absorption features. The flux calibration is done by knowing the magnitude of the telluric standard and using the reddening curve of Rieke & Lebofsky (1985) in order to estimate the flux as a function of wavelength for the given star. This response function is then applied to the target.
Final reduced spectra were corrected for redshift, as determined from the average
z measured from the positions of [S III] 0.953 m,
He I 1.083
m, [Fe II] 1.256
m,
[C I] 0.985
m,
Br
and Pa
.
A galactic extinction correction
of
E(B-V)=0.030
was applied to the data (Schlegel et al. 1998).
The reduced spectra in the interval 0.8-2.4
m
are plotted, in rest wavelengths, in Figs. 1 and 2
for the nuclear and extended regions, respectively.
In addition to the NIR data, optical spectra covering the wavelength
interval 3900-6900 Å were taken using two different
configurations and detectors. The observations were done with the
Cassegrain spectrograph attached to the 1.6 m
telescope at the Observatório Pico dos Dias, Brazil.
In both cases, a slit width of
,
oriented
east-west was employed. The spectrum containing the H
line (5100-6900 Å,
red spectrum hereafter) was obtained on November 28, 2002 (UT) using
a
CCD and a 600 l/mm grating that provided a spectral
resolution of 4.5 Å. The observation covering H
(3900-5300 Å, blue spectrum hereafter) was taken on May 7, 2003
(UT) with a
CCD
and a 900 l/mm grating. The spectral resolution obtained with the latter
configuration was 2.4 Å. Data were reduced following IRAF standard
techniques, that is, bias subtraction and flat field division. Wavelength
calibration were done using HeNeAr lamps. For flux calibration,
a spectrophotometric standard was observed after the targets. Due to bad seeing (
), the procedure for extracting the blue and red spectra was to
integrate all the signal along the slit.
After reduction and flux calibration, the continuum in the overlap region
(5100-5300 Å) disagreed. In order to bring both spectra to the same
continuum level, a scale factor of 0.74 was applied to
the blue spectrum. This number corresponds to the ratio of the
integrated flux of [N I] 5200 Å for the red and blue spectra.
Note that [N I] was chosen because it is the only line common to
both spectra. A final reduced optical
spectrum is plotted in Fig. 3.
Throughout this article, wavelengths corresponding to the optical data, from the blue up to 8000 Å will be given in Ångstroms. The NIR wavelengths longwards of to 8001 Å will be reported in microns.
Table 1: Line identifications and fluxes for the infrared spectra.
In order to measure the flux of the emission lines in the
different spectra, we assumed that the line profiles
can be represented either by a single or by a sum of Gaussian profiles.
The continuum underneath each line was fit by a low-order polynomium,
usually a straight line.
The LINER routine (Pogge & Owen 1993) was employed in this process.
Typically, most emission lines were well fitted by a single or a
sum of two Gaussian components.
Fluxes for the detected lines or upper limits found
in the nuclear and extended regions in the near-infrared
are listed in Tables 1 and 2.
Table 3 reports the measurements on the optical spectra.
In all cases, errors correspond to above the continuum level and reflect the uncertainty
in the placement of the continuum and in the S/N around the line of interest. When a line identification
was uncertain, the Atomic Line List Database
was consulted to
select a list of candidates. A final identification was
possible according to the presence or absence of additional
lines of the candidate ion. Additional support for line identification
was obtained by comparing the position of the emission lines with those
observed in planetary nebulae, for instance.
Due to the differences in spatial resolution between the optical and NIR spectra, we summed up the fluxes of each line measured in the five apertures of the NIR data for all the cases where the analysis involved the simultaneous use of optical and NIR lines. This procedure was applied for consistency and guarantees that the size of both optical and NIR apertures are representative of the same emitting region.
The nuclear spectrum of Mrk 1210 displays a plethora of emission lines, from those emitted from molecular H2and low and medium ionization species such as [C I], [S II], and [S III], up to very high ionization lines of [S IX] and [Si X]. This implies a very wide range of physical conditions of the emitting gas, even within the central 250 pc, to allow the simultaneous coexistence of molecules, somehow protected from the intense radiation of the central engine, and highly ionized gas, undoubtedly associated to the AGN phenomenon. In addition, absorption lines of Ca II and CO are also observed, indicating the presence of a circumnuclear stellar population.
The off-nuclear spectra of Mrk 1210 are also peculiar. At 250 pc NE
and SW of the nucleus, emission from H I and He I dominates.
Farther out, at 500 pc from the center, both the NE and SW spectra
show emission from He I, [Fe II], and [S III]. Weak
Pa
is detected in the 500 pc NE, suggesting that the extended
emission is dominated by matter-bounded clouds. The
detection of extended [S III] is compatible with the data
presented by Fraquelli et al. (2003) who show extended [O III] 5006 Å up to 1 kpc from the nucleus.
The brightest forbidden lines observed in the nuclear spectrum of Mrk 1210 are
[S III] 0.953 m and [S III] 0.906
m. In
addition, Fig. 2 shows that the former line is not limited to the
nuclear region. It extends up to 500 pc NE and 250 pc SW from the center.
In the AP2 spectrum, a
feature at the expected position of [S III] 0.953
m
is detected but the S/N is low. For this reason, we only report
an upper limit for its flux.
An inspection of Fig. 1 shows that Mrk 1210 displays a very
rich coronal NIR spectrum. Here, we report the first detection of
[S VIII] 0.991 m,
[S IX] 1.252
m,
[Si X] 1.430
m, and
[Ca VIII] 2.320
m on this object. The strongest coronal line
observed is [Si VI]
which is, to our knowledge, the only one
previously reported for this galaxy in the 0.8-2.4
m range
(Rodríguez-Ardila et al. 2004).
At 250 pc from the nucleus, high-ionization lines such
as [S VI] 1.964
m, and
[Si X] 1.430
m are still visible, mostly in the SW direction. This,
however, cannot be taken as a definitive probe of extended emission in the
coronal lines. It is possible that they represent scattered nuclear photons
or even that they are emitted in the inner few tens of parsecs of the
extended region.
In addition to the
features mentioned above, Mrk 1210 is characterized by a conspicuous
[Fe II] spectrum, not previously reported in the literature
for other AGN. Besides
[Fe II] 1.256 m and 1.643
m, which is common in most Seyferts
(Rodríguez-Ardila et al. 2004) and whose flux ratio is an
excellent reddening indicator for the NLR, up to eight more lines are detected
in the J and H-bands. The presence of these lines
is a prime opportunity to use the forbidden iron lines as a
diagnostic tool for the region where they are formed.
In Sect. 5 we will
explore the diagnostic capabilities of the [Fe II] as a probe of
the physical conditions of the NLR.
As already mentioned, we found that [Fe II] 1.256
m
extends up to 500 pc NE and SW from
the center (see Fig. 2).
Table 2: Line identifications and fluxes in the NIR for the AP2 and AP4.
Table 3: Line identifications and fluxes for the optical spectra.
Other interesting features on the NIR forbidden spectrum
are [S II] 1.028 m, 1.032
m, 1.033
m,
and 1.037
m
. These lines, along
with the optical doublet [S II] 4068, 4076 Å arise from
the same upper term. It means that
the flux ratio between the NIR and optical sulfur lines depends
primarily on their
transition probabilities. Therefore, it can be used as an additional
reddening indicator.
Finally, low ionization lines of [N I], [C I], and
[P II] are also detected in the nuclear spectrum. With
exception of the later, these lines are also observed in AP1 and AP3 (see Fig. 2). In
the nucleus, the flux ratio
[P II]/[C I] is 1.3 (see Table 1).
The detection of [P II] lines of considerable strength
is surprising. This is because in the sun, phosphorus is a factor 1000 times less abundant than carbon.
Therefore, if the P/C abundance is even closer to solar, the
[P II] lines should not be present unless carbon and other
strong abundant elements are much more optically thick than
they appear. An alternative to this apparent contradiction is
to assume that phosphorus is over-abundant
by a factor 10-30 compared to its solar abundance, relative
to hydrogen. A similar problem is found in high-redshift
quasars, where broad absorption lines of
P V
1118, 1128 Å are detected
(Hamman 1998; Hamman et al. 2002). The NIR [P II] lines
may probably help to set some constraints to the abundance
of phosphorus in AGN. This subject, however, is outside of the
scope of this paper.
The flux of the H2 (1, 0)S(1) 2.121 m, listed in
Table 1 for the different apertures and corrected for
reddening (see Sect. 5), can be used to derive the mass
of hot molecular gas present in the inner few hundred parsecs
of Mrk 1210. To this purpose, we used the expression:
The values of the H2 mass found are
for the nucleus,
for AP1 and
AP3, totalizing
nearly 300
of hot molecular gas for the inner 500 pc of
Mrk 1210. It is in very good agreement with the nuclear H2
masses reported by Reunanen et al. (2003) for a sample of Seyfert 1
and 2 galaxies.
The small mass of hot H2 contrasts with the cold molecular mass of
found by Raluy et al.
(1998) from CO(
)
observations and using the standard
conversion factor found in Galactic giant molecular clouds,
.
The aperture of the CO data
covers the inner 5.7 kpc. Interestingly, Raluy et al. (1998) found
high star-formation efficiency in Mrk 1210 (derived from the
ratio), comparable to that of
normal galaxies like the Milky Way or M 33. In fact, from a list
of 12 AGNs, Mrk 1210 is the one that shows the second highest
star-formation efficiency. They, however,
called attention to the fact that a significant fraction
of the IR luminosity (up to 50%) could come from the AGN
itself, resulting in overestimation of the star formation
efficiency by up to a factor of 2.
A rough estimate of the contribution of the stellar population
flux to the H-band can be made by measuring the depth of
the
overtone bandhead at 1.618
m (see Fig. 1).
We found that the observed
line depth is
of the continuum. This can be
compared to the
that is typically expected for a
population of GKM supergiants that dominate the H-band
light for stellar populations older that 107 yr
(Schinnerer et al. 1998). It means that
of the total
flux in the H-band continuum is due to GMK supergiants.
Although these stars dominate the stellar continuum emission in the H-band,
Schinnerer et al. (1998) claims that about one third of the total stellar
flux is due to stars of other spectral types.
Therefore,
of the H-band continuum
in a
aperture has a stellar origin.
This result agrees very well with the
contribution
of a stellar population older that 107 yr determined by
Cid Fernandes et al. (2001) for this object, based
on optical spectroscopy. It means that the H-band continuum is essentially dominated
by the stellar population, with only
20% due to the AGN. This new value
refines the claim made by Raluy et al. (1998), who estimated
that at least 50% of the IR could come from the active
nucleus. The above calculations suggest that NIR spectroscopy
is a potential tool for unveiling hidden starburst components, either
in highly obscured objects or in those in which the optical
continuum emission is dominated by the AGN component (as in Seyfert 1
galaxies).
Even after correcting for the increased contribution of the stellar population to the continuum, Mrk 1210 continues to display one of the highest (if not the highest) star formation rates among the group of 12 AGN studied by Raluy et al. (1998).
The strong contribution of stellar light to the Mrk 1210 continuum
in the inner
is also compatible with the analysis of
the H2 lines presented in the beginning of this section. It
further supports the hypothesis
that in this AGN the stellar population contributes significantly
to the excitation of the molecular lines.
![]() |
Figure 4:
Profiles of nuclear lines in Mrk 1210 in velocity space.
The lines are normalized to the same peak flux. The small bars in panel
a) represent the instrumental profile of the blue spectrum (
![]() ![]() |
Open with DEXTER |
The large amount of spectral information at medium resolution, covering the optical and NIR regions, allows us to study the form and width of the emission lines in Mrk 1210. The purpose is to map the velocity field of the NLR, examine the presence of genuine BLR components, and estimate the most probable location of the emitting gas.
Table 4: FWHM* and shifts in the peak centroid of the broad component measured in the lines displaying narrow and broad components.
The results of the Gaussian decomposition carried out to measure the line fluxes (see Sect. 2.3) showed that most lines were well-fitted by a single narrow Gaussian component of
Figure 4 shows a comparison of optical and
NIR emission line profiles in velocity space. The upper left panel (a)
shows [O III] 5006 Å, H
and Br
,
while the upper right panel (b) plots [S III] 0.953
m,
[O III] 5006 Å, and Br
.
It can be seen that
both permitted and forbidden lines in the NIR region are
quite similar in form and width, ruling out the hypothesis of a
broad component associated to the BLR. Moreover,
Fig. 4 shows that the optical lines are systematically
narrower at FWHM than the NIR ones. In contrast, the full width at
zero intensity (FWZI) are rather similar in both regions. In
panel (c) a comparison of H
,
Pa
,
and Br
is shown.
Note, again, that the FHWM of the optical and NIR lines are different
but the FWZI is equivalent.
The difference in FWHM among the narrow components may be
attributed to either selective extinction or to
a difference in the instrumental resolution. In the first scenario,
optical photons emitted by low-velocity clouds
are more absorbed than the NIR ones, allowing
us to more likely see the optical emission originated
in high velocity clouds, which apparently are less subjected
to extinction. In the second scenario, the difference
in FWHM between NIR and optical lines is due to
differences in the instrumental profiles. In fact, the instrumental
FWHM of the NIR data is a factor of about 2.2 larger than that of the blue
spectrum (instrumental
km s-1). Although at first sight
the latter argument would be sufficient to explain the differences in
width, the values of FHWM listed in Table 4 point out that very
likely selective extinction dominates. This is because in
Table 4, the FWHM of the lines have been corrected
for instrumental resolution. Clearly, all the lines are spectroscopically
resolved, and most important, the narrow components of
[O III] and H
are narrower than those of
[S III] and Pa
,
for instance.
Table 5: Flux ratios and E(B-V).
Note that the Pa
profile shows enhanced emission
blueward of the line center. This feature is
visible not only in Mrk 1210, but also in other Seyfert galaxies
(e.g. NGC 2110, Storchi-Bergman et al. 1999; NGC 4151, Thompson 1995). We
attribute this asymmetry primarily to the presence of
[Fe II] 1.278
m, which contaminates the blue wing
of Pa
.
In spite of this feature, the Pa
and Br
line profiles are rather similar, and the difference at FWHM with
the optical lines is clear. In panel (d) we compare
[Fe II] 1.643
m, H2 2.121
m and
[S III] 0.953
m. Note the slight blue
asymmetry in [S III], present neither in
[Fe II] nor in H2. It can also be seen that from
peak intensity to FWHM, [Fe II] and [S III] are
indistinguishable. From FWHM to FWZI, the two profiles diverge,
with [S III] becoming broader. In contrast, H2 is
significantly narrower, almost half the FWHM value of the former
two. Since extinction cannot be invoked to explain these
differences because all three lines are located in the NIR region,
the most plausible explanation is that the bulk of [S III],
[Fe II], and H2 originates in the different volumes of gas,
which are subjected to different dispersion velocity conditions.
The results of the above comparison show that the same broad component
present in the permitted lines is also detected in
some of the forbidden lines. This indicates that the gas responsible
for the broad lines is most likely associated to the NLR itself.
Additional support for this picture is found from the work
of Lutz et al. (2002). They compared the emission line profiles
of [Si IX] 3.934 m
and Br
in Mrk 1210 and decided against the BLR interpretation
because the two profiles were
indistinguishable within the S/N limitations.
We propose that the broad component seen in [S III], [O III], and the permitted lines is associated to outflows, instead of to the hidden BLR suggested by Veilleux. Evidence of this scenario comes from Middelberg et al. (2004), who carried out radio observations of Mrk 1210 at 18 cm and 6 cm with EVN, MERLIN, and VLBA. The EVN and VLBA 18 cm images reveal a bright compact object and a weaker component toward the southeast at a distance of 8.6 pc. They marginally detected a third component at a distance of 30.6 pc from the bright compact source. They suggest that the last component could probably be the continuation of the radio ejecta. Moreover, their 6 cm VLBA data resolves the brightest source into an arc of four components, as well as of the south east feature. Note that most of the galaxy's 6 cm emission comes from scales between 3 pc and 16 pc. This region is unresolved from our data. Middelberg et al. (2004) suggest three scenarios for the radio structure in Mrk 1210. A free-free emitting disk, synchrotron emission from a torus or extended accretion disk, and a system composed of a core and an outflow. Based on physical arguments, the last hypothesis is most favored. The resolved arc of four components could be bright knots due to a shock in an extended, low brightness temperature radio outflow, with similar conditions to those in NGC 5506, NGC 1068 and Mrk 231.
An additional support for this picture is that the peak of the broad component detected in [S III], [O III], H I, and He I is blueshifted relative to the systemic velocity of the galaxy (see Table 4). The shift in the forbidden lines is similar to that of He I and almost half of H I. This broad component is probably associated to the series of four knots reported by Middelberg et al. (2004), which is not spatially resolved in our nuclear spectrum. This also would explain how the line profiles of some forbidden and permitted lines in Fig. 4 are different at FWHM but similar at FWZI. The highest dispersion velocities are associated to the outflow component, with little or no extinction toward the line of sight. At FWHM, the bulk of the line profile is dominated by the emission from the classical NLR clouds. The NIR lines are broader because we are looking deep into the clouds, where the dispersion velocity is higher.
The numerous H I and forbidden lines in the spectra of Mrk 1210, spanning a large interval in wavelength, allowed us to evaluate the intrinsic extinction affecting the NLR within the inner few hundred parsec by means of several indicators. In those cases where the determination of reddening only involved NIR lines, we were able to trace it at different distances from the nucleus, obtaining additional information about the dust distribution.
The extinction was determined from the comparison of the predicted and
observed emission line ratios, assuming the standard extinction law of
Cardelli et al. (1989) for RV=3.1; note, however,
that at these wavelengths,
the extinction law is nearly independent of the assumed RV.
Listed in Table 5 are the line ratios
used for the reddening determinations (first column), the intrinsic
ratio values (second column), and the measured line ratios and estimated
reddening for three infrared spectra (AP1, NUC and AP3). For AP2 and AP4,
no sufficient information is available. The last two
rows in Table 5 list the ratios involving optical lines. These last
values correspond to an emitting region size equivalent to the one
covered by the sum of the five infrared
extractions (500 pc). In all cases, the errors correspond to
and
were determined by assuming that the only quantities that introduce uncertainties
are the measured line fluxes.
For the hydrogen lines, we used the intrinsic ratios
for case B given by Osterbrock (1989).
Our results point out to extinction variations up to a factor
of nearly 5 within the NLR.
Recall, however, that the E(B-V) derived
using Pa
may be misleading. This is because of
the contamination
of the Pa
line with [Fe II] 1.278
m, as
mentioned in the previous section.
This feature increases the Pa
flux, leading to
an overestimation of the reddening if it is determined from the
Pa
/Pa
ratio or to an underestimation, if it is
derived from Br
/Pa
.
A previous determination of the reddening in Mrk 1210 was carried out by Veilleux et al. (1997). They calculated a value of E(B-V)=0.48 using the Balmer decrement reported by Terlevich et al. (1991). Although this value is lower than the one determined here from the optical H I, the discrepancy can probably be explained in terms of seeing, aperture, and position-angle effects. Also, our result is slightly above the one found in a recent statistical work presented by Ho et al. (2003), which showed that Seyfert galaxies tend to have an internal mean reddening of about 0.43 (see also Koski 1978).
The detection of the [S II] infrared lines near 1.03 m and
the optical lines at 4068, 4076 Å, all
arising from the same upper 2P term, allows us to
estimate the interstellar extinction affecting the S+ region.
This method was first pointed out by Miller (1968) and tested by
Wampler (1968, 1971) in several Seyfert galaxies,
and more recently by
Greve et al. (1994) in the Orion nebula.
Due to the large separation in wavelength between the
NIR and optical [S II], the reddening derived from this
ratio should be very accurate.
Using the atomic data from Keenan (1991), we obtained an [S II]
near-infrared to optical theoretical ratio of 0.67.
The observed ratio and the reddening coefficient are given in Table 5.
The E(B-V) of
found for the S+ region is the smallest one
derived for the nuclear gas of Mrk 1210,
indicating that these lines, as expected, are formed in the outer portions
of the NLR.
In addition to the extinction determined from the H I and
[S II] lines, the rich [Fe II] spectrum detected
in Mrk 1210 offers a unique opportunity to study the extinction
affecting the region where it is emitted. This is because
each of the three pairs of lines involved in the [Fe II] ratios listed in Table 5 share
the same upper level. It means that the intensity ratio depends only
on energy differences between the lines and their Einstein A-coefficients,
making these ratios useful probes of reddening.
The scarce detection of such a rich [Fe II] spectrum
in other AGN has prevented its use as diagnostics
of the internal extinction. For this reason, extinction
estimates with the NIR [Fe II] lines are unusual
and, to our knowledge, this is the first time that they have been
applied to a Seyfert nucleus (except for the ratio
[Fe II] 1.256 m/1.643
m). Intrinsic line ratios
were taken from Bautista & Pradhan (1998), calculated from
Nussbaumer & Storey (1988) transition probabilities.
Our results, listed in Table 5, clearly show that high values of extinction affects the [Fe II] emission region. That all three [Fe II] reddening indicators agree within the errors points out that all [Fe II] is emitted in the same dusty environment. Although the uncertainties in some ratios are large, mostly because the lines involved are intrinsically weak, the high extinction derived for the [Fe II] emitting region leads us to propose that the bulk of these lines is formed deep in the NLR, separate from the zone where other low-ionization species are emitted (i.e. [S II] and [O I]). In contrast to the nucleus, the reddening found for AP1 and AP3 shows that the extended emission region are barely affected by dust, if at all (see Table 5).
The picture that emerges from the above analysis points out that the dust distribution in the inner 250 pc of Mrk 1210 is rather inhomogeneous. Much of the optical emission we see should be emitted in the outer portions of the NLR, while the NIR one comes from the deep NLR. This is supported by the analysis of the permitted and forbidden lines made in Sect. 4. At FHWM, the optical narrow components are narrower than the NIR ones. Moreover, that the broad components display a similar width in velocity space, from the visible through NIR, suggests that they are being emitted in a region not affected by dust. Whether the high extinction measured from the [Fe II] and the H I lines is associated to the circumnuclear starburst, to the host galaxy, or even to the torus cannot be easily distinguished from our data. However, we should recall that Martini et al. (2003) used HST images to classify Mrk 1210 as a tightly-wound nuclear dust spiral with the individual dust lanes traced over a full rotation about the nucleus. This view supports the idea that the dust is mainly concentrated in a planar geometry, probably associated to the host galaxy. In this scenario, the broad component arises in the outflow, which is inclined, from the dust plane relative to our line of sight. Further support for this scenario comes from the dramatic transition between a Compton-thick, reflection-dominated state and a Compton-thin state in Mrk 1210 reported by Guainazzi et al. (2002) from X-ray observations. They claimed that this transition can be explained if the Compton-thick and Compton-thin absorbers are different. The former is associated to the torus while the latter, located on larger scales, may be associated with the host galaxy rather than with the nuclear environment.
Electronic densities and temperatures for the nuclear gas of Mrk 1210 can be determined by means of several diagnostic line ratios. But before this can be done, line fluxes need to be corrected in accordance with the value of reddening obtained from lines of the same ion or similar ions. Thus, low-ionization lines such as [O I], [N I], [N II] and [S II] were corrected by E(B-V)=0.28. This value was determined from the [S II] lines. The exception of this low-ionization group is [Fe II], whose lines were corrected by E(B-V)=1.5, as measured from the different indicators employed. The remaining lines were corrected for E(B-V)=0.5. It results from averaging out our determination of internal extinction by means of optical H I lines, the one reported by Veilleux et al. (1997) for Mrk 1210 and the average internal extinction reported by Ho et al. (2003).
In the low-density regime, the [S II] 6716/
6730
flux ratio is a well-known diagnostics of the
electronic density that is relatively insensitive to temperature.
This diagnostic can be applied to regions where the electronic density is
lower than the critical density of the [S II] lines
(
cm-3). Above this density the lines become
collisionally de-excited. The observed [S II] ratio of 0.9
indicates an electron density
cm-3 for
temperatures in the interval 10 000 K to 50 000 K.
On the other hand, the infrared [S II] 1.03
m lines are
observed to
be quite strong relative to the 6716, 6730 Å doublet
in the nuclei of Mrk 1210, indicating that a second,
higher density (
cm-3) S+ zone
can also exist. The [S II] 1.03
m (
6716+
6730) line ratio is
very sensitive to variations in density, in particular for values of
cm-3 where the optical [S II]
diagnostic no longer varies with
.
The observed
[S II] 1.03
m (
6716+
6730) line ratio is
0.59, which indicates an electron density
13 000 cm-3
for an assumed temperature
K (Keenan 1991).
The differences between the densities derived from the two
[S II] line flux
ratios illustrate that the [S II] lines observed in the spectrum
of the central region of Mrk 1210 arise from different zones with
different densities: a denser region (
13 000 cm-3)
that produces the bulk of the
[S II]
4068,
4076 and 1.03
m emission and a lower density
region (
1000 cm-3) that contributes with most of the
6716 and
6730 flux.
Table 6: Density and temperatures found from different diagnostic line ratios.
The density of the Fe+ emitting region can be derived from the observed
[Fe II] 1.533 m/1.643
m and [Fe II] 0.861
m/1.256
m line
ratios (Pradhan & Zhang 1993). The former ratio yields a density
cm-3 for
in the
interval 3000-12 000 K, while from the latter ratio we derive
a higher density,
cm-3, assuming
K. These results, together with the ones from the
previous sections, strongly suggest that the [Fe II] lines are indeed
emitted in a separate
region, characterized with higher densities and higher extinction.
The gas temperature can be estimated from measurements of pairs of
emission line ratios emitted by a single ion from two levels
with considerably different excitation energies. Two good examples
are the [O III] (4958+
5006)/
4363 and
[N II] (
6548+
6583)/
5754
line ratios. Temperature determinations from these ratios were calculated,
using the TEMDEN task included in the STSDAS Version 3.3
package of IRAF. The atomic data employed by the routine are at least as
recent as those given in the compilation by Pradhan & Peng (1995).
The results obtained for a wide range
of densities are plotted in Fig. 5, which shows very high
temperatures in the O+2 and N+ zones, ranging from 18 000 K
up to 50 000 K.
Although the ionization potential required to form the O+2 and
N+ ions (35.1 and 14.5 eV, respectively) are different, the
critical densities for the upper levels of the
lines involved in the
[O III] and [N II] line ratios are
104.9 cm-3.
Thus it is reasonable to suppose that these lines arise from
the same region with
.
Therefore, we set a common temperature of
K,
for an assumed density
cm-3.
Electron temperatures
K are difficult to obtain
in the O+2 and N+ zones in photoionized gas,
requiring a source of energy
input in addition to photoionization. Our working hypothesis
is that this source is associated to shock heating.
![]() |
Figure 5: Electron densities and temperatures found for Mrk 1210. |
Open with DEXTER |
Another well-known temperature diagnostic is the
[S III] (0.906+0.953) m/
6312 line ratio, which is analogous
to the [O III] (
4958+
5006)/
4363 ratio. Nonetheless, it
is not commonly used, mainly due to the problem of deblending
[S III]
6312 from the usually considerably stronger
[O I]
6300. Our medium-resolution optical spectrum allowed us
to separate these two lines without introducing large uncertainties
in the line fluxes.
From the [S III], we
calculated the temperature of the S+2 region using the TEMDEN task
of IRAF STSDAS nebular package.
The results are also plotted in Fig. 5.
Note that the
calculation is relatively insensitive to density, indicating a
gas temperature of about 40 000 K for densities in the interval
10
cm-3.
Osterbrock et al. (1990) determined the S+2 gas temperature
for several Seyfert galaxies, finding temperatures from 10 000 K to
the extreme case of 26 100 K for the high-ionization Seyfert 1
galaxy III Zw 77. The unusually high temperature found for the S+2 emitting region of Mrk 1210, together with the high temperature
derived from the [O III] line ratio, shows further
evidence that shocks contribute
to the ionization structure in the NLR of this galaxy.
Recall that these
lines are the ones with a blueshifted broad component that we
associated to the radio-knots observed in VLBA maps.
Also plotted in Fig. 5 are the electron temperatures
in the N0 and S+ zones, derived from
the temperature sensitive [N I] 5197/1.040
m and
[S II] (
6716+
6730)/
4068 line ratios, respectively.
The ratio [O I] (
6300+
6363)/
5577 also works
as a temperature indicator.
Although the line [O I]
5577
is definitely present in the optical spectrum of Mrk 1210,
the profile is too noisy to measure it accurately. Therefore, we could only
set an upper limit for the line flux, which we used to derive an upper limit
of the [O I] emitting gas temperature of 12 000 K. The low temperature
of the gas,
in addition to the low ionization potential of the ion indicates that
these lines are probably emitted in the same region as
[S II] and [N I].
The results found from the different indicators (see Table 6)
show that the NLR of Mrk 1210
cannot be characterized by a single value of temperature and density, but that
presents a noticeable gradient in the physical conditions.
The temperatures range from about 10 000 K up to the unusual value of
40 000 K determined from the [S III] lines, although the latter
value may be misleading because of the contribution from shocks.
The range of densities
extends at least over three orders of magnitude, from a low density
region (
cm-3) in which the
[S II]
6716,
6730 lines (and probably the [N I] and
[O I] lines) are emitted to a high density region
(
cm-3) where the
[Fe II] lines are emitted. Note that densities that high are
plausible for the [Fe II] given that the critical densities of most
if its transitions are in the interval 104-105 cm-3
(Nisini et al. 2002).
We have carried out a detailed study of the physical properties of the nuclear
and extended emission line regions of the Seyfert 2 galaxy Mrk 1210.
For this purpose, we presented spectroscopic observations
simultaneously covering the interval 0.8-2.4 m.
The spectrum in the range 0.81-1.03
m is the first one published for this
object. For that reason, all lines detected in this region were not previously
identified in this source.
The spectra of Mrk 1210 display a plethora of emission lines from those emitted by molecular H2, prominent mostly in the K-band, through low and medium ionization lines such as [C I], [S II], and [S III] to very high ionization lines of [S IX] and [Si X]. In addition, Mrk 1210 is characterized by a rich [Fe II] spectrum, not previously reported in the literature for other AGN. Extended emission of [S III], [Fe II], and He I is found up to a distance of 500 pc from the nucleus.
The large amount of spectral information at medium resolution, covering the optical and NIR regions allows us to study the physical conditions in the inner few hundred parsecs of this galaxy. Our main results can be summarized as follows:
Acknowledgements
This research was partly supported by the Brazilian agency CNPq (309054/03-6, 473465/2004-3) to ARA and the European Commission's ALFA-II program through its funding of the Latin-American European Network for Astrophysics and Cosmology, LENAC to X.M. The authors thank the referee Enrique Pérez for his useful comments that improve this manuscript. This research made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.
![]() |
Figure 1: Nuclear spectrum of Mrk 1210 in rest wavelength. The labels mark the position of the most important emission lines identified. |
Open with DEXTER |
![]() |
Figure 2:
Spectra of the extended emission region of Mrk 1210 in rest wavelength.
Each panel represents the spectrum in a rectangular
beam, which is
![]() ![]() |
Open with DEXTER |
![]() |
Figure 3: Optical spectrum of Mrk 1210 in rest wavelength. The labels mark the position of the most important emission lines identified. |
Open with DEXTER |