A&A 476, 177-192 (2007)
DOI: 10.1051/0004-6361:20078479
J. P. Pérez-Beaupuits1,2 - S. Aalto2 - H. Gerebro2
1 - Kapteyn Astronomical Institute, Rijksuniversiteit Groningen, 9700 AV Groningen, The Netherlands
2 -
Onsala Rymdobservatorium, Chalmers Tekniska Högskola, 439 92 Onsala, Sweden
Received 14 August 2007 / Accepted 27 September 2007
Abstract
Aims. Bright HNC 1-0 emission, rivalling that of HCN 1-0, has been found towards several Seyfert galaxies. This is unexpected since traditionally HNC is a tracer of cold (10 K) gas, and the molecular gas of luminous galaxies like Seyferts is thought to have bulk kinetic temperatures surpassing 50 K. There are four possible explanations for the bright HNC: (a) large masses of hidden cold gas; (b) chemistry dominated by ion-neutral reactions; (c) chemistry dominated by X-ray radiation; and (d) HNC enhanced through mid-IR pumping. In this work, we distinguish the cause of the bright HNC and to model the physical conditions of the HNC and HCN emitting gas.
Methods. We have used SEST, JCMT and IRAM 30 m telescopes to observe HNC 3-2 and HCN 3-2 line emission in a selection of 5 HNC-luminous Seyfert galaxies. We estimate and discuss the excitation conditions of HCN and HNC in NGC 1068, NGC 3079, NGC 2623 and NGC 7469, based on the observed 3-2/1-0 line intensity ratios. We also observed CN 1-0 and 2-1 emission and discuss its role in photon and X-ray dominated regions.
Results. HNC 3-2 was detected in 3 galaxies (NGC 3079, NGC 1068 and NGC 2623). Not detected in NGC 7469. HCN 3-2 was detected in NGC 3079, NGC 1068 and NGC 1365, it was not detected in NGC 2623. The HCN 3-2/1-0 ratio is lower than 0.3 only in NGC 3079, whereas the HNC 3-2/1-0 ratio is larger than 0.3 only in NGC 2623. The HCN/HNC 1-0 and 3-2 line ratios are larger than unity in all the galaxies. The HCN/HNC 3-2 line ratio is lower than unity only in NGC 2623, which makes it comparable to galaxies like Arp 220, Mrk 231 and NGC 4418.
Conclusions. We conclude that in three of the galaxies the HNC emissions emerge from gas of densities
,
where the chemistry is dominated by ion-neutral reactions. The line shapes observed in NGC 1365 and NGC 3079 show that these galaxies have no circumnuclear disk. In NGC 1068 the emission of HNC emerges from lower (<
)
density gas than HCN (>
). Instead, we conclude that the emissions of HNC and HCN emerge from the same gas in NGC 3079. The observed HCN/HNC and CN/HCN line ratios favor a PDR scenario, rather than an XDR one, which is consistent with previous indications of a starburst component in the central regions of these galaxies. However, the
column density ratios obtained for NGC 3079 can be found only in XDR environments.
Key words: galaxies: ISM - galaxies: starburst - galaxies: active - galaxies: Seyfert - radio lines: galaxies - radio lines: ISM
The hydrogene cyanide, HCN molecule, is commonly used as an extragalactic tracer
of molecular gas with densities n(H2) larger than
(e.g. Solomon et al. 1992; Curran et al. 2000;
Kohno 2005). The HCN to CO intensity ratio varies significantly,
from 1/3 to 1/40 in starburst galaxies, and it has not been determined
whether this variation depends on the dense molecular gas content or on
the abundance and/or excitation conditions. In addition, recent results
seems to indicate that HCN may not be an unbiased tracer of the dense
molecular gas content in LIRGs and ULIRGs (Graciá-Carpio et al. 2006).
It is essential, therefore, to use other molecular tracers than HCN, in
order to understand the physical conditions in the dense gas.
A molecule of particular interest, for comparison with HCN, is its isomer
HNC. The detection of interstellar HNC supports the theory of dominant
ion-molecule chemistry in dark molecular clouds. Both species are thought
to be created by the same dissociative recombination of HCNH+. This ion
can produce HCN and HNC, with approximately equal abundances. Models based
only on this scheme would predict then an HNC/HCN ratio 1. However, the
observed HNC/HCN abundance ratios vary significantly between different kinds
of molecular clouds - the ratio ranges from 0.03 to 0.4 in warm cores (
K),
and can be as high as 4.4 in cold cores (
K).
The CN (cyanogen radical) molecule is another tracer of dense gas, with a lower
(by a factor of 5) critical density than HCN. CN is also chemically linked to
HCN and HNC by photodissociations (e.g. Hirota & Yamamoto 1999). Surveys
of the 1-0 transition of CN and HNC have been done in order to trace a cold,
dense phase of the gas in luminous galaxies (Aalto et al. 2002). It was
found that the HNC 1-0 luminosities often rivalled those of HCN 1-0.
These results seem to contradict the idea of warm (
K) gas in the
centers of luminous galaxies (e.g. Wild et al. 1992; Wall et al.
1993) whose IR luminosities were suggested to originate from star formation
rather than AGN activity (Solomon et al. 1992).
Table 1: Sample of galaxiesa.
According to observations in the vicinity of the hot core of Orion KL, experimental
data, chemical steady state and shock models, the HNC/HCN ratio decreases as the
temperature and density increase (e.g. Schilke et al. 1992; Talbi et al.
1996; Tachikawa et al. 2003). If a bright HNC 1-0 transition
line is nevertheless detected under these conditions, it could be due to the following
possible explanations: (a) the presence of large masses of hidden cold gas and dust at
high densities (
); (b) chemistry dominated by ion-molecule reactions
with HCNH+ at low density (
)
in regions where the temperature
dependence of the HNC abundance becomes weaker; (c) enhancement by mid-IR pumping, also
in low density regions where the lines would not be collisionally excited; and (d) the
influence of UV-rays in Photon Dominated Regions (PDRs) and/or X-rays in X-ray Dominated
Regions (XDRs) at densities
and at total column densities
(Meijerink & Spaans 2005).
In the case of CN, observations of its emission towards the Orion A molecular complex (Rodríguez-Franco et al. 1998) suggest that this molecule is also enhanced in PDRs, but particularly in XDRs, where a CN/HCN abundance ratio larger than unity is expected (e.g., Lepp & Dalgarno 1996; Meijerink et al. 2007).
We have observed low and high transition lines of the HCN, HNC and CN molecules in a group of Seyfert galaxies, which are supposed to host both sources of power, AGN and starburst activity, in their central region. Our interest is to assess the excitation conditions of HCN and HNC, distinguish between the above possible causes of the bright HNC, and to explore the relation between the CN emission, XDRs and dense PDRs in these sources.
In Sect. 2 we describe the observations. The results (spectral lines, line intensities and line ratios) are presented in Sect. 3. The interpretation of line shapes and gas distribution in the most relevant cases, as well as the possible explanations for the bright HNC and the modelling of the excitation conditions of HCN and HNC are discussed in Sect. 4. The conclusions and final remarks of this work are presented in Sect. 5.
We have used the James Clerk Maxwell Telescope (JCMT) in 2005 to observe the HNC J=3-2 (271 GHz) and the HCN J=3-2 (267 GHz) lines towards a sample of Seyfert galaxies. Observations of CN and HNC J=1-0 (90 GHz) were made in 2002 using the Swedish ESO Southern Telescope (SEST). HCN J=1-0 (88 GHz) data from literature were also used. The system temperatures ranged between 350 K and 430 K. In the case of CN J=1-0 (226 GHz) the weather conditions were not so good, making the system temperature range between 490 K and 760 K. Pointing was checked regularly on SiO masers and the rms was found to be about 3''.
Table 2: Beamsizes & efficiencies.
![]() |
Figure 1:
Molecular line emissions in NGC 1068.
The velocity resolution was set to 25
![]() ![]() ![]() ![]() |
Open with DEXTER |
For the SEST observations we alternated between a 500 MHz and 1 GHz backend, depending on weather. Simultaneous observations with the 1 and 3 mm receiver were taking place. In addition, the data of the Seyfert galaxy NGC 3079 obtained during the IRAM survey in 2006 is included in this work. We used the software package XS (written by P. Bergman) to reduce the data and fit the Gaussians. Beamsizes and efficiencies are shown in Table 2.
The sample consists of five Seyfert galaxies of which two are Seyfert 1 type and
three are considered mainly Seyfert 2. Table 1 lists the
coordinates of the center positions observed in these galaxies, their
sub-classification as Seyfert galaxies and heliocentric radial velocities,
according to the NASA/IPAC Extragalactic Database (NED)
(http://nedwww.ipac.caltech.edu/). The distances were calculated assuming
a Hubble constant of
Mpc-1
(Ngeow & Kanbur 2006).
The molecular line emissions observed in NGC 1068 are shown in
Fig. 1. A first-order polynomial was used
in most cases to correct the baselines, with the exception of
the HCN 3-2 and CN 2-1 spectra, for which a second-order polynomial was
required. The velocity resolution was set to 15
for CO and CN, whereas a
25
resolution was used for HNC and HCN. These velocity resolutions
represent less than 10% of the line widths. The spectra are centered
with respect to the heliocentric systemic velocity
(from NED).
The J=1-0 and J=2-1 lines of CO show a triple structure where the peaks can be attributed mainly to the bars and the spiral arms, as described by Helfer & Blitz (1995).
We detect two of the main spingroups of CN 1-0: the 1-0
(
J = 3/2-1/2,
F = 5/2-3/2) line at the
center of the spectrum, and the 1-0 (
J = 1/2-1/2,
F = 3/2-3/2) line shifted 856
to the right.
This spectrum shows three components as well, however they are
hard to distinguish due to the blending.
In fact if we freely fit three Gaussian components, the uncertainties
of the center velocity, amplitude and line width are about 100% or
larger. Instead, if we set the line width of the central component
to 154
,
which corresponds to the line width found for the
CN J=2-1 line as described below, we get reasonable values. Only
the first Gaussian component shows high uncertainties in the
amplitude and line width. Since the beam sizes of the CN and CO
molecules are similar at the frequencies of the J=1-0 lines,
we think that both beams pick up emission coming from the spiral
arms and bars. Although, the nuclear region seems to be the predominant
component in the case of CN.
Due to the second spingroup is corrupted by noise in the backend,
we set the amplitude of the Gaussian to 8.3 mK,
which corresponds to the expected factor of about 0.43 times the
amplitude of the main spingroup, according to the National Institute
of Standards and Technology (NIST).
The resulting central velocity of the second component was 1946
(
113.169 GHz), i.e., about 22 MHz shifted from
the expected frequency for this spingroup.
Table 3: NGC 1068 line parameters.
![]() |
Figure 2:
Molecular line emissions in NGC 1365. The spectra
are centered with respect to the heliocentric systemic velocity
![]() ![]() ![]() ![]() |
Open with DEXTER |
In the CN 2-1 line, the two spingroups at 226.8746 GHz (2-1,
J = 5/2-3/2,
F = 7/2-5/2) and 226.6596 GHz
(2-1,
J = 3/2-1/2,
F = 5/2-3/2) are
severely blended since the shift is 300
.
In order to identify the two spingroups we first fit two
Gaussian components to get the line width of the main spingroup,
which is
.
We then set this value to both Gaussian components and
fit again the other parameters.
The second spingroup is, in the optically thin limit, a factor
0.54 weaker than the main spingroup. We get a factor
0.6 between the
intensities obtained from the Gaussian fit.
On the other hand, the resulting center velocity of the second
component is 1289
,
which corresponds to a shift of about 68 MHz with respect to the expected frequency of the second
spingroup. Besides the noise in the data, this shift may also
be produced by the influence of the unresolved spingroups,
(2-1,
J = 3/2-1/2,
F = 1/2-1/2) and
(2-1,
J = 3/2-1/2,
F = 3/2-1/2),
located in between the two main spingroups. Together, these
inter-spingroups would produce an intensity comparable to
that of the second spingroup.
The line shape and intensity of the CN 2-1 spectrum differs from
the one obtained by Usero et al. (2004,
thereafter U04). Besides the lack of baseline coverage observed
in the U04 spectrum, there is a substantial discrepancy between the
estimated source sizes. In U04 the emission of most of the molecules was
assumed to emerge from a
region, which corresponds to
their estimate of the size of the circumnuclear disk (CND),
based on the CO 2-1 high resolution
map presented by Schinnerer et al. (2000).
Table 4: NGC 1365 line parameters.
If we correct the peak antenna temperature (
mK) of our
CN 2-1 spectrum, in order to obtain the main beam brightness temperature
in the same way done in U04 (
),
we get a peak
of about 0.42 K, which is about 0.18 K (or a
factor of
1.8) larger than the peak temperature obtained by U04.
Due to the chemical link between HCN and CN, here we rather estimate
the source size based on the high resolution map of HCN 1-0, published
in Helfer & Blitz (1995). We estimate that the emissions of
HCN, CN (and very likely HNC as well) emerge from a more extended region
of about 10
.
This source size can also be inferred from
the HCN 1-0 map presented by Kohno et al. (2001).
By correcting the peak antenna temperature (as described in Sect. 3.5) we get
a peak
of about 0.12 K, which is a factor 2 smaller than what was found by U04.
The HCN 3-2 line has its center velocity at
which
coincides with the main HCN emission shown in the position-velocity
(p-v) map by Tacconi et al. (1994) (this is discussed in Sect. 4).
The line width and line shape of this spectrum are consistent with HCN 1-0 spectra published in previous work (Nguyen-Q-Rieu et al. 1992;
Curran et al. 2000).
The HNC spectra present a double peak profile in which the main
peak has its center velocity at
in both the J=1-0and J=3-2 transitions. However, the line shape and intensity of the HNC 1-0 spectrum
differs from the one obtained by Hüttemeister et al. (1995).
The HNC 3-2 spectrum seems to be affected by noise, since the amplitude and line
width of the secondary component have uncertainties larger than 50%.
The different line shapes observed between the HNC and HCN spectra would
indicate that they emerge from regions of different kinematics. The line
parameters are summarized in Table 3.
Figure 2 shows the molecular line emissions observed in NGC 1365. Without considering the second spingroups of CN, all the molecular spectra of NGC 1365 present a double peak structure, irrespective of the beam size or line transition.
This double peak structure can also be seen in the CO 2-1 and CO 3-2 spectra showed in Sandqvist et al. (1995) and in Sandqvist (1999), respectively.
The two strongest spingroups of the CN 1-0 line are also
detected in this galaxy, with a double peak structure at the
nuclear region. The strongest peak was used as reference for the main
spingroup. We first fit two Gaussian components to the main group and
then we set the line width of the second spingroup, observed
at around 2400
,
to the value found for the main
spingroup. The proportion found between the amplitude of the spingroups is
about 0.3 instead of the expected value of 0.43, according to NIST.
The data around the second spingroup are strongly affected
by noise, so its central velocity is shifted by about 50
from the
expected value of about 800
.
In the CN 2-1 line, the second spingroup overlaps with the
double structure itself. The line widths of the two spingroups
are set according to the value found for the CN 1-0 line.
Given their combined intrinsic line
strengths, the integrated line intensity ratio of the two
spingroups is expected to be about 1.85 (assuming optically
thin emission). If we also set the line width of the second peak
of the double structure at 212
,
as found for CN 1-0, we get a
proportion of about 0.4 between the two spingroups. Instead, if we
set the amplitude of the second spingroup to 5.7 mK (the expected proportion)
and we fit the line width of the double peak structure, we get a value
of 181
,
which is about 30
less than the width found in the CN 1-0
line. In any case, the center velocity of the second spingroup is shifted by
about 30
more than the expected velocity.
The HNC 1-0 and HCN 3-2 spectra have a structure similar to that of the CN 1-0 spectrum. Note that in CO 1-0 the double-peak is more pronounced than in the CN, HCN or HNC spectra, where the low velocity peak dominates. The line parameters are summarized in Table 4.
The molecular line emissions observed in NGC 3079 are shown in
Fig. 3. The spectra are centered with respect
to the heliocentric systemic velocity
.
The
velocity resolution was set to 20
for HNC and HCN, and to 10
for CO and CN. Four structures are observed in the CO lines, which are
not present in the other spectra. The CO lines are picking up extended
lower density gas, compared to the CN lines (which have similar beam
size). Hence the difference in line shape.
The CN 1-0 line present a double peak structure.
Only the main spingroup (
J = 3/2-1/2,
F = 5/2-3/2) of the CN 1-0 line is
observed since the second spingroup falls beyond the bandwidth of the backend.
In the CN 2-1 line, the two main spingroups, (
J = 5/2-3/2,
F = 7/2-5/2)
and (
J = 3/2-1/2,
F = 5/2-3/2), are detected. But they are severely blended
and merged with the double peak structure. In this case, only two Gaussians
were fitted because a third (middle) component presented large uncertainties in
both the amplitude and line width, due to the noise in the spectrum. The two
components are separated by only 254
,
and the proportion between their
amplitudes is about 0.65. The difference between these and the expected values
is attributed to the noise and the blending of the emission of the second
spingroup.
The HCN spectrum of the J=1-0 and J=3-2 transitions agree with the HCN 1-0 spectrum obtained by Nguyen-Q-Rieu et al. (1992).
The HNC 1-0 spectrum is different from the tentative detection presented by
Hüttemeister et al. (1995), which is broader and more intense than the spectrum
showed here. Besides the similar line shapes, our HCN and HNC spectra
extends from 900
to 1300
.
Instead, the HNC spectrum showed in Hüttemeister
et al. (1995) peaks at around 1300
.
Hence, we believe that our HNC detections are
correct. In this case, because of the similar line shapes and
widths (although the HNC 3-2 line is strongly affected by noise, and its peak
intensities have uncertainties between 40% and 50%.), the HCN and HNC emissions
are likely emerging from the same gas. The line parameters are summarized in
Table 5.
![]() |
Figure 3:
Molecular line emissions in NGC 3079. The
velocity resolution was set to 20
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
Molecular line emissions in NGC 2623 (left)
and NGC 7469 (right). The velocity resolution was set to 25
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The HNC 3-2 line emission observed in NGC 2623 is shown in the
left top panel of Fig. 4. The
velocity resolution was set to 25
.
The spectrum is centered
with respect to the heliocentric systemic velocity
(NED). After a total integration time of 2 h we do not
detect HCN 3-2 emission in this galaxy.
The right top panel of Fig. 4 shows the
CO 2-1 spectrum observed in NGC 7469. The velocity resolution was set
to 15
.
The spectrum is centered with respect to the systemic
velocity
.
The CO spectrum shows a double
peak structure, like the one observed in NGC 1365. We observed the HNC 3-2
line for about 1 h of integration time, but we do not detect any
emission. The bottom panels show the spectra of the not detected lines.
The line parameters of the detected transition lines are
summarized in Table 6.
The velocity integrated intensities are showed in Table 7.
In order to compute the line intensity ratios
were corrected for the different beam sizes obtained with different
frequencies and telescopes, according to the correction factors
defined for compact sources in Rohlfs and Wilson (2003).
In the case of extended sources (
)
the ratios
were corrected for the main beam filling factor, which was approximated
as
.
The beams, as well as
the source structures, were assumed to be Gaussians.
Table 5: NGC 3079 line parameters.
The source sizes reported in Table 1 were estimated from
high resolution maps available in the literature. For NGC 3079 a source size of
was estimated for the CO 1-0 emission, considering only intensities above 15% of the peak integrated intensity of the contour map presented by Koda
et al. (2002). In NGC 1068 most of the CO 1-0 emission emerges from
the two spiral arms, with the largest extension of
40'' (e.g. Helfer & Blitz 1995;
Schinnerer et al. 2000). Considering intensities above 20% of
the peak emission of the high resolution map by Schinnerer et al. (2000) the
source size of the CO 1-0 emission in NGC 1068 was estimated as
.
For NGC 2623 a source size of
was estimated from the CO 1-0 map presented in Bryant et al. (1999),
which agrees well with the estimate made by Casoli et al. (1988).
In the case of NGC 1365 the CO emission is concentrated in the nuclear and bar regions
(Sandqvist et al. 1995). Hence the source size of the CO emission was
estimated as
,
corresponding to intensities above 20% of the peak
emission in the CO 3-2 map by Sandqvist (1999).
From Papadopoulos & Allen (2000) and Davies et al. (2004)
the source size of the CO emission in NGC 7469 was estimated as
,
which
correspond to intensities above 40% of the peak emission in the high resolution map
by Davies et al. (2004). The criteria of selection of the source size of CO varies depending on the gradient of the emission observed in the different sources.
Table 6: Line parameters of NGC 2623 and NGC 7469. The third column indicate the Gaussian components (G.C.) used to fit the spectral lines.
The source size of the high density tracers CN, HCN and HNC were estimated
through HCN maps available in the literature. In all these maps,
we observed that the HCN 1-0 emissions emerge mainly from the nuclear region
of the galaxies. Since there are no published maps of CN nor HNC emissions
for the galaxies studied here, the corresponding source sizes
were considered equal to that of the HCN emission, due to their chemical link.
Since there are no HCN maps for NGC 1365 and NGC 2623, the factor
found in NGC 1068
was used to estimate the source size of the HCN 1-0 emission,
based on their corresponding
.
The estimated
are shown in Table 1.
The line ratios were computed assuming an error of 10% in the
reported beam efficiencies
,
5% of error in the main
beam
(Table 2), and a 10% error in
the estimated source sizes
(Table 1).
The obtained ratios are shown in Table 8. From
these ratios we can conclude that CO, as well as the high density
tracers - CN, HCN and HNC - are subthermally excited.
We compare the spectral shape of the high density tracers with HCN and CO position-velocity maps available in the literature in order to address the location of the dense gas.
Table 7: Velocity-integrated intensitiesa.
Table 8: Line intensity ratios.
In NGC 1068 the two outer peaks of the CO 1-0 spectrum (Fig. 1)
coincide with the maximum double peak emission seen in the CO position-velocity
(p-v) map obtained by Helfer & Blitz (1995), so they can
be attributed to the emission emerging from the spiral arms. The center is
attributed to the emission emerging from the CND, as can be inferred from
the CO 2-1 spectrum (Fig. 1). This
center peak coincides with the maximum HCN emission at 1100
seen in the
HCN p-v map obtained by Tacconi et al. (1994). This p-v
map is shown in the left panel of Fig. 5 for comparison,
along with the corresponding scaled spectra observed in NGC 1068.
The main spingroup of the CN 1-0 line (Fig. 1) shows
a shape similar to that of the CO 1-0 line. So contributions of emission coming
from the CND, as well as from the spiral arms, can also be inferred. The main
spingroup of the CN 2-1 spectrum also coincides with the maximum HCN emission
around 1100
.
The HCN spectra seem to contain two components, although we fit only one Gaussian to the spectrum. The main
component fits the region in the p-v where the strongest HCN emission is coming from.
The secondary component might be attributed to the secondary peak emission
of HCN, observed around 1035
.
Another smaller peak is seen in the
p-v map around 1255
,
but this component is not detected in the
HCN spectra.
In the HNC spectra, instead, two clear components are observed. The main component
has a center velocity of 1073
which lies in between the two main peaks of the
HCN emission observed in the corresponding p-v map. With respect to the CO p-v map (Helfer & Blitz 1995), the main HNC component comes from a region
where the CO emission is faint. Instead, the secondary component of HNC
seems to emerge from a region around 1250
,
which corresponds to the secondary peak
emission in the CO p-v map, and roughly to the third peak
of the HCN emission in the respective map. Although, the latter is uncertain due
to the noise in the spectra.
![]() |
Figure 5: (Top-left panel) Position-velocity ( p-v) map of HCN 1-0 emission in NGC 1068, adapted from Tacconi et al. (1994). The spectra below are the respective emission lines of Fig. 1, re-rescaled to fit the velocity scale of the HCN p-v map. (Top-middle panel) Position-velocity map of CO 3-2 emission in NGC 1365, adapted from Sandqvist (1999). The spectra below are the re-scaled versions of the ones shown in Fig. 2. (Top-right panel) Position-velocity map of CO 1-0 emission in NGC 3079, adapted from Koda et al. (2002). The spectra below are the re-scaled versions of the ones shown in Fig. 3. |
Open with DEXTER |
On top of the middle panel of Fig. 5 we show the CO 3-2
p-v map of the central region of NGC 1365, from Sandqvist (1999).
The spectra below correspond to the re-scaled spectra shown in Fig. 2.
The double peak structure of the spectra coincides fairly well with the double peak
emission observed in the p-v map. The left components of the spectra seem to
emerge from a region around 1530
.
The right components of most of the spectra could emerge from the region around 1710
,
with the exception of the CN 1-0 line. The right component of CN 1-0 seems
to emerge from around 1750
which, in turn, coincides with the maximum emission
level of the right peak in the CO 3-2 p-v map.
The two peaks in the CO 3-2 p-v map are of about the same intensity, which
is reflected in the 12CO 1-0 spectrum, where both components have intensities
K. The 12CO 3-2 and 13CO 1-0 data obtained by Sandqvist
(1999) also have this feature. However, the CN, HCN and HNC spectra,
exhibit a gradient between the intensities of their high and low velocity components.
This can also be seen in the HCN 1-0 and HCO+ 1-0 spectra
obtained by Sandqvist (1999). This intensity gradient could be due to either
a larger abundance of the species, or to a higher excitation in the corresponding
region of the lower velocities.
Since the observed double peak structure does not change with the beam size (one single feature), we think that the nucleus of NGC 1365 lacks a circumnuclear disk, i.e., it is consistent with a Seyfert 1 nucleus.
The four-peak structure observed in NGC 3079 can be identified in the CO 1-0 p-v map of Koda et al. (2002). The top-right panel of Fig. 5 shows an adaptation of the Fig. 10 in Koda et al. (2002) with 0'' offset along the minor axis of NGC 3079. Since this galaxy is almost edge-on, most of the CO emission is probably coming from the spiral arms, besides the nuclear region. Instead, the high density tracers are expected to be mostly nuclear, as in the case of NGC 1068. This can explain the double-peak structure of the line shape of the CN, HCN and HNC molecules, in contrast to the CO lines.
Note that in this case the CN spectra are slightly dominated by the lower-velocity peak, as
in the case of NGC 1068 and NGC 1365, whereas the HCN and HNC spectra of NGC 3079 are dominated
by the higher-velocity peaks. The lower and higher velocity peaks are not perfectly aligned.
The maximum separation between the peaks is 80
for the lower-velocity peaks and
60
for the higher-velocity peaks. The lower-velocity peaks are centered around
1000
and the higher-velocity peaks around 1270
,
which means the high density tracers
tend to avoid the peaks of the CO 2-1 emission.
Since the line shapes of CN, HCN and HNC are similar, and they do not change with the beam size, their emissions likely emerge from the nuclear region in this galaxy, in contrast to the case of NGC 1068, where the line shape of the HCN emission changes with the beam size.
Table 9: Intensity ratios between the high density tracers.
![]() |
Figure 6: Excitation conditions modelled for the 3-2/1-0 line ratios of HCN (top) and HNC (middle) observed in NGC 1068. The conditions required for the HCN and HNC molecules overlap in a narrow region. The relative column densities in the overlap zone of the excitation conditions of these molecules is shown in the bottom plot. The optical depth in the whole region explored for HCN ranges between 0.03 and 10 in the J=1-0 line, and between 0.32 and 30 in the J=3-2 line. In the case of HNC the optical depth ranges between 0.01 and 30 in the J=1-0 line, and between 0.003 and 30 in the the J=3-2 line. The optically thin limit of both molecules and lines is depicted by the right edge of the excitation conditions, whereas the optically thick limit correspond to the left edge of the figures above. |
Open with DEXTER |
The HCN/HNC 1-0 ratio is about 2.0 in NGC 1068
and NGC 3079, and 1.4 in the others. This indicates a brighter
HNC emission in NGC 2623, NGC 1365 and NGC 7469. Luminous HNC in galaxies
may have the following plausible explanations:
The ratio 0.4 found in NGC 2623 is in the limit between case (a) and (b). On the other hand, case (c) and (d) cannot easily be ruled out. More
information is required in order to distinguish between the proposed scenarios,
as suggested by Aalto et al. (2007).
The HCN 3-2/1-0 line ratio is below 0.3 only in NGC 3079, whereas it is
larger than 0.4 in NGC 1068 and NGC 1365. This result is interesting since it
implies that in NGC 3079 both HCN and HNC emission emerge from the same gas,
whereas in NGC 1068 the HNC emission has to emerge from a lower (<
)
density gas than
HCN. It would be interesting to see if this result holds for NGC 1365. We expect
to obtain the HNC 3-2 data for this galaxy in a future project.
The pumping scenario (c) should lead to a HNC 3-2/1-0 line ratio close to unity. This is not observed either in the HNC or in the HCN data we have. However, it is not possible to rule out this scenario since low excitation may be the result of mid-IR pumping of low density gas. Detailed modelling is needed in this case.
We used the radiative transfer code RADEX
(Van der Tak et al. 2007) to explore a wide range of possible excitation conditions
that can lead to the observed line ratios. This code is sensible to the column density of
a molecule per line width, and uses a constant temperature and density of the collision partner,
which in this case is H2. Another limitation of RADEX is that it cannot handle large optical depths
(
). Our analysis is not depth dependent and we assume a homogeneous sphere for the escape
probability approach. Hence, our models aim to reproduce a sort of average cloud that represent
the physical conditions of the emitting gas, which is a well fitted starting model for single dish observations, where all the emissions detected are convolved with the telescope beams.
The grid consists of densities between 104
and
,
temperatures between 4 and 200 K, and column densities per line width
between 1010 and
.
Excitation maps were generated in order to obtain the total
column density per line width (
)
as function of
the kinetic temperature (
K) and the number density of molecular hydrogen
(
). The countour lines of the maps describe the dichotomy
between temperature and density. This means that, for a given column density,
the observed 3-2/1-0 line ratio can be obtained with high temperatures and
low densities, or low temperatures and high densities.
The hyperfine structure of the HCN J=1-0 transition line
is not included in this study. The extrapolated HCN data of the LAMDA
database were used instead (Schöier et al. 2005).
Note that the hyperfine components may
be overlapping and may interact radiatively in AGN-like
environments. This process requires further analysis and modelling that is
not included in RADEX.
We are not able to generate an excitation map of
CN due to the lack of collision data for this molecule.
Although, an study to extrapolate the collision data for CN from
other know molecules (e.g. CS) is ongoing.
Figure 6 shows the excitation condition maps for the HCN and HNC molecules, modeled from the average ratios observed in NGC 1068. These maps show that HCN (top) requires higher density gas and higher column density than HNC (middle) in order to obtain the observed 3-2/1-0 line ratios.
Considering a kinetic temperature
K and a density
,
as estimated by Tacconi et al. (1994), we found that
,
which agrees with the respective value found by Tacconi et al.
According to the recent PDR and XDR models by Meijerink & Spaans (2005)
and Meijerink et al. (2007), under PDR conditions a temperature
K
can be reached at a total column density of
,
if the Habing
flux is
(a Habing flux G0 = 1 corresponds to a far UV flux
erg
s-1),
or at
if
.
At those depths the density
would be slightly higher than
,
though. On the other hand, in an XDR environment
a temperature of 80 K and a density of
can be reached at a slightly larger column density of
,
if the impinging radiation field is
erg
s-1,
or at a much larger depth equivalent to
if the radiation flux is
erg
s-1.
However, at at a temperature of 80 K and density of
there is no solution for HNC. This means that either
the HNC emission arises from a gas with different physical conditions than HCN,
or that both molecules trace a cooler (
K) and lower density gas
(
), if the emission of both molecules arise from the
same gas.
If the HCN and HNC emissions actually trace the same gas,
the range of possible excitation conditions can be constrained to the zone
where the conditions for HCN and HNC overlap. The bottom plot of
Fig. 6 shows the average overlap zone, where the contours
correspond to the ratio between the total column densities of HCN and HNC.
The overlap zone is defined by the optically thin limit (right edge of the excitation maps)
of the HNC transition lines and the optically thick limit (left edge of the excitation maps)
of the HCN lines. The models show that
is larger than
in all the
overlap zone. According to Fig. 10 in Meijerink & Spaans (2005),
column density ratios lower than unity can be found mostly in PDRs,
but also in XDR environments if the total column density
is lower than
.
If we still assume the same temperature (80 K) proposed by Tacconi et al.
(1994), and we consider the contour line where
is about 2.5 orders of magnitude larger than
,
we find that the gas density should be
,
which agrees with the observed HNC 3-2/1-0 line ratio discussed
in Sect. 4.2. However, this density would not be consistent with the result found by Tacconi et al. (1994)
nor with the observed HCN 3-2/1-0 line ratio, which implies that the HNC emissions arise
from a more diffuse gas than HCN.
![]() |
Figure 7:
Excitation conditions modelled for the 3-2/1-0 line ratios
of HCN (top) and HNC (middle) observed in NGC 3079.
The excitation conditions required for these
molecules overlap in most of the range explored. This suggest that
their emissions emerge from gas with the same physical conditions.
The bottom plot shows the overlap zone.
There is a large zone of excitation conditions where the column density of
HNC is between 3 and 10 times larger than that of HCN.
The contour lines of -0.5 depict the zone of the physical conditions for which the
column density of HNC is about 3 times larger than that of HCN.
For both molecules and transition lines the optical depth ranges between 0.001 and 30.
The optically thick limit is basically defined by the limit of convergence of RADEX
(
![]() |
Open with DEXTER |
According with the PDR and XDR models by Meijerink & Spaans (2005),
a column density ratio
larger than unity can be found only in XDR
environments at a total column density
,
with radiation fields
of about 1.6 erg
(or equivalent Habing flux
)
and total density
.
column density ratios larger than unity at
can also be
found with stronger radiation fields (
erg
)
and higher densities
(
). On the other hand,
can be larger than unity
at
if the radiation fields are on the order of 1.6 erg
and the
total density is about
.
In our models we did not explore densities lower than
since the HCN and HNC molecules are
expected to trace higher densities. If we assume a kinetic temperature of 80 K, as in the case of NGC 1068,
and a density
,
the column densities per line width would be about 1012 and
for HCN and HNC, respectively. In the case of HNC both
transition lines are equally optically thick at 80 K in the whole range of densities, whereas HCN
is more optically thin. At a density of
the optical depth of HNC is
1 in both
lines, whereas
is
0.3 and
0.1 in the J=1-0 and J=3-2 lines of HCN, respectively.
A density lower than
is more likely to be the right case in order to be consistent with
the line ratios observed and discussed in Sect. 4.2.
Both NGC 2623 and NGC 1365 have similar excitation conditions for HNC and HCN,
respectively. As in the case of HCN in NGC 1068, both molecules could be embedded in
gas with densities larger than
at 80 K, which is consistent with
line ratios larger than 0.3, as discussed in Sect. 4.2.
The column densities per line width of these molecules would be about
for HNC and about
for HCN, if the gas density is
.
![]() |
Figure 8: Excitation conditions modelled for the 3-2/1-0 line ratio of HNC (top) and HCN (bottom) observed in NGC 2623 and NGC 1365, respectively. For NGC 2623 the optical depth of HNC 3-2 ranges between 0.003 and 10 in the J=1-0 line, and between 0.32 and 30 in the J=3-2 line. In the case of NGC 1365 the optical depth in the HCN 3-2 line ranges between 0.1 and 30, and between 0.01 and 10 in the J=1-0 line. |
Open with DEXTER |
The CN/HNC ratio increases towards the CND in both galaxies, NGC 1068 and NGC 3079. The CN/HNC 1-0 ratio is lower than unity in all the galaxies, with the exception of NGC 1068. This can be interpreted as that HNC is more abundant than CN, assuming the same excitation conditions and that the emission emerge from the same gas and distribution in the galaxy. However, since we do not really know the source size of CN (and basically neither that of HNC) the beam dilution effects could be more (or less) severe than estimated here.
Although the beam dilution uncertainty of CN also applies to the CN/HCN ratio, this case is more interesting because none of the galaxies where we do have data show ratios larger than one, neither for the global ratio (J=1-0) nor for the nuclear part of the galaxies (higher transitions). All of these galaxies are considered active galaxies, so the presence of an AGN increases the chances of finding an X-ray dominated region (XDR) in their nuclear regions, as it seems to be the case of NGC 1068, according to Usero et al. (2004).
In an XDR the CN molecule is expected to be more abundant than HCN (e.g., Lepp & Dalgarno 1996, Meijerink & Spaans 2005; and Meijerink et al. 2007), and hence the CN/HCN intensity ratio could also be larger than unity, especially for the higher transition lines - CN 2-1 and HCN 3-2 - where the beam dilution is less important than for the J=1-0 line. However, this is not seen in any of the sources presented here, nor in the sample of galaxies shown by Aalto et al. (2002), where results for AGN are presented along with starburst galaxies.
According to Meijerink et al. (2007), in a PDR environment the
column density ratio ranges between 2.0 and 0.5 for densities between 104 and
,
respectively, whereas in an XDR this ratio varies from over a 1000 (at
)
to 40 (at
). If the CN intensity
lines are proportional to the column density, and the estimate of the
beam dilution is accurate enough, our results would favor a PDR scenario, rather than an XDR,
with densities <
in the central regions of all these galaxies. However, if
the CN molecule is indeed more abundant than HCN, then the weaker intensity lines
could be also due to a stronger optical depth effect (in the escape probability sense) for CN.
A rigorous modeling and analysis of the CN molecule would be required in order to
understand and predict the intensities of this molecule.
On the other hand, the PDR scenario would also be consistent with the results found for the HNC and HCN molecules
described above. HCN/HNC line intensity ratios larger than unity are usually found
in PDR environments at total column densities
lower than
,
while
this ratio is larger for strong (>10 erg s
)
radiation fields and low (
10
)
densities in an XDR environment (Meijerink et al. 2007).
The lower densities required to observe this ratios
in an XDR environment tend to dismiss this alternative since our simulations in Sect. 4.3
favor densities
10
for the HNC and HCN molecules, as in the case of NGC 1068.
In particular, there is evidence of recent starburst in the dense nuclear
disk of NGC 1068 (Davies et al. 2006). These results are consistent
with a model considering both, AGN and starburst components, required for modelling
the UV to FIR atomic spectrum of NGC 1068 (Spinoglio et al. 2005).
We have used the SEST and JCMT telescopes to carry out a survey of CN 2-1, HCN 3-2 and HNC 3-2 line emission in a sample of 4 Seyfert galaxies, plus NGC 3079 which was observed with the IRAM 30m telescope. The conclusions we draw are as follows:
Modeling of the collision data for the CN molecule would be useful to estimate the
column density ratio, which would complement the
ratio in order to have a more sophisticated tool
to estimate and distinguish the prevalent environment conditions of the high density
gas in the nuclear region of Seyfert galaxies.
The AGN contribution (through XDR effects) is typically of a small angular scale and can be seriously affected by beam dilution at the transition lines studied in this work. On the other hand, the starburst contribution is of a larger angular scale than the AGN, and it effects can be contaminating our observations, and hence leading to the favored PDR scenario found with our models. Hence, our suggested interpretations could change if we zoom in on these sources. Therefore, high resolution maps of HNC and CN molecules are necessary to complement those of HCN, and to do a more accurate estimate of molecular abundances and line intensity ratios, which take source size into account. Observations of the higher transition lines (e.g. J=4-3) can also aid to disentangle the effects of the AGN and the starburst ring, due to the smaller beam size obtained at higher frequencies.
Acknowledgements
We are grateful to the SEST staff for their help during the observing run. We thank S. Curran and A. Polatidis for help with some of the SEST observations. We are grateful to S. Hüttemeister, M. Spaans, and R. Meijerink for discussions. We also thank F. van der Tak and J. Black for their help and discussions about RADEX. Molecular databases that have been helpful include LAMBDA and NIST.