A&A 376, 837-852 (2001)
DOI: 10.1051/0004-6361:20011039
C. Horellou 1 - J. H. Black 1 - J. H. van Gorkom 2 - F. Combes 3 - J. M. van der Hulst 4 - V. Charmandaris 5
1 - Onsala Space Observatory, Centre for Astrophysics & Space Science,
Chalmers University of Technology,
439 92 Onsala, Sweden
2 -
Astronomy Department, Columbia University, 538 W 120 Street, New York, NY 10027,
USA
3 -
Observatoire de Paris, DEMIRM, 61 avenue de l'Observatoire, 75014 Paris, France
4 -
Kapteyn Institute, Department of Astronomy,
Postbus 800, 9700 AV Groningen, The Netherlands
5 -
Cornell University, Astronomy Department,
106 Space Sciences Building,
Ithaca, NY 14853, USA
Received 4 May 2001 / Accepted 12 July 2001
Abstract
We present and interpret
observations of atomic and molecular gas toward the southern elliptical galaxy
NGC 1316 (Fornax A), a strong double-lobe radio source with a disturbed optical morphology
that includes numerous shells and loops.
The 12CO(1-0), 12CO(2-1), and H I observations
were made with SEST and the VLA.
CO emission corresponding to a total molecular hydrogen mass of
was detected towards the central position as well as northwest and
southeast of the nucleus in the regions of the dust patches.
The origin of that gas is likely external
and due to accretion of one or several small gas-rich galaxies.
H I was not detected in the central region of NGC 1316, but
of atomic gas was found
towards the giant H II region
discovered by Schweizer (1980)
located
(or 36.2 kpc) from the nucleus.
H I was also found at three other locations in the outer part of NGC 1316.
The H I distributions and kinematics of the two nearby spiral companions of NGC 1316,
NGC 1317 (a barred galaxy to the north) and
NGC 1310 (to the west) could be studied. Both galaxies have unusually
small H I disks that may
have been affected by ram-pressure stripping.
Key words: galaxies: ellipticals - individual (NGC 1316, NGC 1317, NGC 1310) - galaxies: clusters: Fornax cluster - radio lines: galaxies
![]() |
Figure 1:
Radio continuum contour map at ![]() |
Open with DEXTER |
Over the past fifteen years, evidence
for the existence of a complex, multi-phase
interstellar medium in elliptical galaxies
has accumulated.
Most of the gas in ellipticals is hot
(107 K) and radiates at X-ray wavelengths, but there is also a cooler
component of H I, CO, ionized gas and dust
(e.g. Bregman et al. 1992).
The origin and fate of the gas
are still debated, but it seems that the hot gas is associated with
the spheroid, whereas the atomic and molecular gas as traced by the CO
characterize a prominent disk
(e.g. van Gorkom & Schiminovich 1997).
The detection rates of H Ior CO are significantly
higher among ellipticals with peculiar features
(such as tidal tails, loops, shells, dust lanes) than in the
whole class of ellipticals, which suggests that the gas is
of external origin and due to accretion of smaller companions and/or
to major mergers.
NGC 1316 (Fornax A, Arp 154, PKS 0320-37)
on the outskirts of the Fornax cluster presents several characteristics
of a merger remnant.
The disturbed outer morphology with numerous shells and loops was
noted by Matthews et al. (1964) and Arp (1964). The galaxy
was extensively studied in the optical by Schweizer
(1980,
hereafter S80) and Schweizer (1981) who
classified it as a D-type galaxy with an elliptical-like
spheroid surrounded by a large envelope.
The brightness profile approximately obeys the r1/4 law.
The galaxy has an inner dust lane of length 2'
oriented close to the minor axis.
S80 also reported the discovery of a giant H II region
(hereafter denoted SH2)
from the nucleus, and the existence of an inner disk
of ionized gas probably associated with the dust lane and
extending out to 54''.
The disk is in rapid rotation
(
kms-1).
To explain those features and the lack of a nearby perturber,
S80 suggested that NGC 1316 accreted one or several gas-rich
companion galaxies about
-
years ago.
Mackie & Fabbiano (1998) reported the discovery of an extended
emission-line region (hereafter denoted EELR)
located in a tidal tail
from the nucleus
and suggested on the basis of its position and extreme size that it is related
to the merger process.
NGC 1316 is the
third nearest strong radio galaxy in the sky after NGC 5128 (Centaurus A) and
M 87 (Virgo A) (see Fig. 1).
The radio emission arises in two giant lobes
separated by 33'
and located mainly outside the optical galaxy, at a position angle of
(Wade 1961; Ekers et al. 1983;
Fomalont et al. 1989).
In the central arcmin a bent radio jet
abruptly terminates at the location of a dust lane, which may have deflected
it
(Geldzahler & Fomalont 1984).
The nucleus is a bright UV source, unresolved by HST
(Fabbiano et al. 1994).
A study of the stellar kinematics of NGC 1316 showed that the galaxy is rotating
rapidly around the minor axis at a position angle of
(Bosma et al. 1985).
More recently Arnaboldi et al. (1998)
measured radial velocities of 43 planetary nebulae in the outer parts
of NGC 1316
to determine the distribution of angular momentum.
They derived a constant velocity v= 110 kms-1 along
the major axis of the galaxy in the radial range R=50-200''.
They also probed the kinematics of the inner region (R <80'')
using absorption line spectra.
They confirmed the alignment between the photometric and the kinematic minor axes
but also found evidence for a disturbed kinematical pattern
beyond
.
X-ray emission from NGC 1316 has been observed with several instruments.
The
ratio of NGC 1316 is low compared to that of E and S0 galaxies
(Einstein observations by Fabbiano et al. 1992;
Kim et al. 1992),
but in the same range as that of dynamically young ellipticals which may have
experienced mergers (Fabbiano & Schweizer 1995;
see also Sansom et al. 2000).
In order to place NGC 1316 into the merger sequence as discussed by Hibbard & van Gorkom (1996) it is necessary to know the distribution of the cool gas, both the neutral atomic hydrogen and the molecular hydrogen as traced by the CO. Sage & Galleta (1993) detected 12CO(1-0) line emission at two positions along the inner dust lane and with kinematics similar to that of the ionized gas. In this paper, we present more extensive and more sensitive maps of the 12CO(1-0) and the first 12CO(2-1) observations in the central part of NGC 1316 as well as a search for CO emission toward the outer H II complexes. We also present results from a VLA search for H Igas in NGC 1316 and its surroundings.
The paper is organized as follows: we first describe the H Iand CO observations, then discuss the H Iresults for the companion galaxies NGC 1317 and NGC 1310. We then focus on the H Iand CO results for NGC 1316 and the implications for the merger picture.
The HI observations consist of data collected with the Very Large Array
(VLA)
in its spectral-line mode.
The galaxy was observed in the hybrid DnC configuration:
the 1.3 km D-array configuration with the antennas of the north arm in the
3 km C-array configuration to improve the North-South angular resolution, thus
compensating for the low declination of NGC 1316.
The observations consisted of two runs of 5 hours each.
A 6.25 MHz bandwidth was used, corresponding to a total velocity range
of 1300 kms-1.
On-line Hanning smoothing has been applied to the data after which every other channel
was discarded, leaving a set of 31 independent channels centered at 1774 kms-1
and spaced by 41.7 kms-1.
The observing parameters are summarized in Table 1.
The data were inspected, calibrated, and maps were produced using NRAO's Astronomical
Image Processing System (AIPS).
The continuum was estimated from channels free of line emission and was subtracted
in the UV plane.
Line channel maps were made using uniform weighting with a ROBUST factor of 1 which,
after some experimentation, was found to give the best combination of
sensitivity and spatial resolution.
The channel maps were CLEANed and a correction for primary beam attenuation was applied
to the entire data cube.
Maps of the atomic hydrogen column density and velocity field were produced.
The maps were not smoothed. Only pixels above a flux cut-off of 1.5 mJy/beam
(2.5 )
were included in the summation of the channel maps.
The data were taken in 1985 by van Gorkom, Ekers, Schweizer and Wrobel.
Due to limitations in on-line and off-line
computing power only 25 antennas and one polarization (RR) could be used.
The data were fully reduced by van Gorkom, but never published. Although NGC 1310,
NGC 1317 and the HII region were unambiguously detected, the quality
of the data reduction was poor. This paper shows that reprocessing
of ancient (but venerable) data is very worthwhile, especially when a large
spectral dynamic range is needed, because there have been major improvements
in software such as the UVLIN program within AIPS and far better weighting
schemes in the imaging tasks.
Phase center (B1950) |
![]() ![]() |
Velocity center | 1774 km s-1 (heliocentric) |
Primary beam ( FWHM) | 30
![]() |
Flux calibrator | 3C 48 |
Phase calibrator | 0332-403 |
Bandpass calibrator | 3C 48 |
Array | DnC |
Date | 19, 21 Oct. 1985 |
Synthesized beam | |
- FWHM: major ![]() |
52
![]() |
- Position angle | +2![]() |
Bandwidth | 6.25 MHz |
Number of channels | 31 |
Channel separation | 41.7 km s-1 |
Time on source | 6 h 51 min |
Noise level (1![]() |
|
- Flux density | 0.6 mJy beam-1ch-1 |
- Column density |
![]() |
- H Imass |
![]() ![]() |
Note | NGC 1316 | NGC 1317 | NGC 1310 | SH2 | EELR | |
RA(1950) | (1) | 03![]() ![]() ![]() |
03![]() ![]() ![]() |
03![]() ![]() ![]() |
03![]() ![]() ![]() |
03![]() ![]() ![]() |
DEC(1950) | (1) | -37![]() |
-37![]() |
-37![]() |
-37![]() |
-37![]() |
![]() |
(2) | 1783 | 1948 | 1739 | 1690 | |
D [Mpc] | (3) | 18.6 | 18.6 | 18.6 | 18.6 | 18.6 |
Type | (4) | .LXSOP. | .SXR1.. | .SAS5*. | ||
Size | (5) | 12
![]() |
2
![]() |
2
![]() |
1
![]() |
81
![]() |
PA [![]() |
(5) | 50 | 78 | 95 | ![]() |
|
![]() |
(6) | 9.4 | 11.81 | 12.56 | ||
LB [![]() |
(6) | 9.35 ![]() |
1.02
![]() |
5.1 ![]() |
||
![]() ![]() |
(7) | 2.05 ![]() |
2.45 ![]() |
7.7 ![]() |
||
M(H I) [![]() |
(8) | <108 | 2.7 ![]() |
4.8 ![]() |
1.9 ![]() |
<3 ![]() |
M(H2) [![]() |
(9) | ![]() ![]() |
3.1 ![]() |
<0.2 ![]() |
<7 ![]() |
<7 ![]() |
![]() |
(8) | 1960, 132 | 1827, 144 | 1690, ![]() |
1750, ![]() |
|
![]() |
(8) | <0.1 | 3.3 ![]() |
![]() |
0.2 | 0.1 |
![]() |
(9) | 1948, 84 |
(3) Distance estimate for NGC 1365 based on HST measurements of cepheid variable stars (Madore et al. 1999). At that distance, 1' corresponds to 5.4 kpc. That distance was adopted for all the systems discussed here.
(4) RC3.
(5) For the three galaxies: optical diameters at the 25 mag arcsec-2 isophote taken from RC3; SH2: Schweizer (1980); EELR: Mackie & Fabbiano (1998). (6) RC3.
(7) Derived from the fluxes at 60 and 100 m listed in the IRAS
catalogue using the relation:
where D is the distance in Mpc, f60 and f100 are the fluxes in Jy,
and the far-infrared luminosity
is in solar luminosities.
(8) This work.
(9) NGC 1316, SH2 and EELR: this work. H2 mass of NGC 1317 and NGC 1310: Horellou et al. (1995) after correction for the distance adopted here.
The observations were done with the 15 m Swedish-ESO Submillimeter
Telescope (SEST)
on La Silla, Chile, in two sessions on 1999 December 7-13
and in 2001 January.
We used the IRAM dual receiver to observe simultaneously at the
frequencies of the 12CO(1-0) and 12CO(2-1) lines.
At 115 GHz the half-power beamwidth of the telescope is 43'', which
corresponds to a linear scale of about 3.9 kpc at the adopted distance of
18.6 Mpc.
The CO(2-1) map has a resolution of 22'' (or
2 kpc).
This is the highest resolution map of molecular gas in NGC 1316 available
to date.
A balanced on-off beam-switching mode was used.
The pointing was regularly checked on nearby SiO masers.
The temperature scale used throughout this article is the main-beam
temperature
.
The main-beam efficiency of SEST
is
at 115 GHz and 0.5 at 230 GHz.
The data were analyzed using the software CLASS.
The spectra were smoothed to a final velocity resolution
of
14.5 kms-1.
Only first-order baselines were subtracted.
The integrated HI masses are calculated using:
HI column densities can be derived from the observed integrated intensities by
using the relation:
From the 12CO(1-0) line intensities we can estimate the mass of molecular
gas (H2) within the 43'' SEST beam:
![]() |
Figure 2:
Integrated HI map overlaid on an optical image from the Digitized Sky Survey.
The contour levels are
![]() ![]() |
Open with DEXTER |
Figure 2 shows the HI field of interest around NGC 1316 superimposed on an
optical image from the Digitized Sky Survey.
HI emission was found in the two companion galaxies NGC 1317 (to the north)
and NGC 1310 (to the west).
No HI was detected in the main body of NGC 1316 to a limiting mass of
108
but some was
found at four locations in the outer part of the galaxy.
One of these is associated with the outer HII region SH2.
The other is coincident with
the extended emission line region EELR.
The two other HI concentrations to the west and northwest of NGC 1316
are located at the edge of an optical plateau between the area of X-ray emission
as observed by ROSAT (Mackie & Fabbiano 1998) and the
extended western radio lobe (see Fig. 1).
The mass of each of these HI concentrations is
.
We will defer discussion of these HI features and the lack of HI emission
in the main spheroid of NGC 1316 until later and first discuss the HI properties
of NGC 1317 and NGC 1310.
![]() |
Figure 3:
a) HI integrated intensity map of NGC 1317 overlaid on a grey scale
optical image from the Digitized Sky Survey.
A correction for the primary beam has been applied.
The levels are
(4, 8, 16, 31, 47) ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
The HI line emission of NGC 1317 displayed as channel maps running from
2066 kms-1 (top left panel) to 1857 kms-1 (bottom right panel).
North is up and east is left.
The levels are 1.5, 3, 4, 5, 6, 7, 8, 9 mJy/beam.
The first contour corresponds to the 2![]() |
Open with DEXTER |
![]() |
Figure 5:
a) HI integrated intensity map of NGC 1310 overlaid on a grey scale
optical image from the Digitized Sky Survey.
A correction for the primary beam has been applied.
The levels range from
21 to
![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
The HI line emission of NGC 1310 displayed as channel maps running from
1857 kms-1 (left panel) to 1732 kms-1 (right panel).
North is up and east is left.
The levels are 4, 8 and 12 mJy/beam.
The first contour corresponds to the 2![]() |
Open with DEXTER |
NGC 1317 is a small barred galaxy north of NGC 1316.
The two galaxies have a projected separation on the sky of 6
3
(34 kpc) and a radial velocity difference of +165 kms-1.
An inner bar perpendicular to the primary bar is visible in the optical images
of Schweizer (1980) and Wozniak et al. (1995).
A ring of star formation is seen in the optical images and in
the ultraviolet image obtained with the
Ultraviolet Imaging Telescope.
The primary bar is "fat" (
60''
80'').
The secondary bar has a position angle of
60
and
a total length of
14''.
S80 assigns NGC 1317 a late SB(rs)a or early SB(rs)b type in the
Hubble classification.
Figure 3 presents in form of isocontours
the atomic hydrogen distribution
and the velocity field in NGC 1317 derived from the observations.
Also shown is the global HI spectrum.
The channel maps are shown in Fig. 4.
The HI distribution is more extended
than the beam size and is elongated in the direction of the minor axis of the galaxy.
The gas is clearly rotating, with the south-west side receding.
The velocity and the linewidth of the global HI spectrum agree with those of the
12CO(1-0) line detected
by Horellou et al. (1995) (see Table 2).
NGC 1317 has comparable amounts of molecular and atomic gas
(
).
The HI disk of NGC 1317 is small compared to the optical extent.
The HI isophote corresponding to a column density of 1020 cm-2
crosses the major axis of the galaxy
at a radius of about 55'' whereas the optical isophote
at 25 mag arcsec-2 is found at 82
5. The ratio of the extent of the
HI disk to that of the optical image is therefore 0.7, which places NGC 1317, if
classified as an Sab galaxy,
in the group of galaxies with an HI distribution strongly affected by the
environment (group III of Cayatte et al. 1994 for Virgo galaxies).
Cayatte et al. (1994) showed that those Virgo galaxies with a truncated
HI disk have a normal central HI surface density average as compared to field galaxies
and concluded that the galaxies (all located close to the cluster core) have been
affected by ram-pressure
sweeping. NGC 1317 has a central HI surface density averaged
over half the optical radius of
cm-2,
which is lower than that of Sab galaxies in the comparison sample of Cayatte et al.
(
cm-2), but in the range of values found for galaxies
of earlier types (S0, S0/a, Sa have
cm-2). Given
the uncertainty in the morphological type and the fact that we are discussing
an individual galaxy, it is difficult to draw any firm conclusion about the central HI
density of NGC 1317. Nevertheless, it is clear that the HI disk is small and that
NGC 1317 has a low HI mass for its optical size and morphological type. We calculate
an HI deficiency of .64 in the log, or a factor of 4.4 using the average value for
isolated Sab galaxies observed by Haynes & Giovanelli (1984) and the definition:
HI
where cs stands for
control sample.
Although NGC 1316 is not located in the center of the Fornax cluster,
it does dominate its corner of the cluster and is a source of X-ray emission.
It is possible that NGC 1317
has passed through the X-ray halo of NGC 1316 and been stripped of some of its atomic gas.
NGC 1310 is an almost face-on spiral fortuitously positioned between us and the western
radio lobe of Fornax A. Ionized gas from NGC 1310 produces a Faraday screen that
depolarizes the background radiation (Fomalont et al. 1989;
Schulman & Fomalont 1992).
Although the galaxy is located 20
5 away from the center of Fornax A
on the plane of the sky
and in the outer region of the primary beam of the VLA, HI emission is well detected.
Figures 5 and 6 present the results of the HI observations in the same form as those
for NGC 1317. The approaching gas is located on the eastern side of the galaxy.
If the spiral arms are trailing, then the northern part of NGC 1310 is the near side. The
HI distribution is slightly asymmetric and the emission is stronger on the eastern side.
We estimate a total mass of atomic hydrogen of
for NGC 1310, which is about 20 times higher than the upper limit on the H2 mass
set by Horellou et al. (1995).
Comparing the extent of the HI disk to that of the optical galaxy, one finds that the
isocontour of column density 1020 cm-2 roughly coincides with the blue
isophote at the 25th mag arcsec-2. This is unusual compared to Sc galaxies in
the field, which normally have an HI diameter larger than their D25 diameter
by a factor of
(Cayatte et al. 1994).
However, NGC 1310 is not significantly HI -deficient compared to the isolated Sc
galaxies of Haynes & Giovanelli (1984) (HI
in the log).
The small size of the disk (
2' across)
together with the location of NGC 1310 in the outer part of
the primary beam of the VLA make it difficult to say more about its HI properties.
![]() |
Figure 7:
a)
Contour map of the integrated HI line emission in the HII region SH2.
The levels are 2, 3, 4 and 5 mJy/beam.
The first contour corresponds to the 2![]() |
Open with DEXTER |
![]() |
Figure 8:
Top: CO(1-0) and CO(2-1) spectra obtained in December 1999.
The offsets are in arcseconds.
The x scale ranges from
1000 to 2300 kms-1. The y scale
(main-beam temperature,
![]() |
Open with DEXTER |
The spheroid.
Huchtmeier & Richter (1989)
give an upper limit on the HI integrated line intensity of 4.5 Jy km s-1,
which corresponds
to
at our adopted distance.
We didn't detect any HI in the central region of the galaxy down to a
5
limit of 107
per beam and per channel.
If, however, the HI emission is as
broad in velocity as the CO detected in the central region (see next section)
then the
limit is
.
This could be even higher if
the HI is extended over several beams. Very extended HI, moreover, will
be missed in the present observations because the quality of the data at
the shortest spacings is limited by the presence of 34 Jy of continuum
emission. A safe limit to the HI emission in the main body of NGC 1316
is 108
.
We did not detect HI absorption towards the nuclear continuum source
which has a flux density of 220 mJy at 1412 MHz. The 5
column
density limit is
cm-2, where
is the spin temperature (in K) of the absorbing atomic
hydrogen.
The four HI detections in the outer regions of NGC 1316 are all far
outside the main body of NGC 1316 and lie at or close to the edge of the
faint optical shells and X-ray emission of NGC 1316. Each clump
contains 1 to
of HI. The southern clumps are
associated with optical emission line gas: the HII region SH2
and the EELR.
SH2.
This HII region has a size of
(135 pc
316 pc);
it is located
,
or 36.2 kpc south of the nucleus and has a radial velocity of
kms-1 relative to it (S80).
S80 notes that a "blue fuzz", probably an associated complex of OB stars,
is visible near the HII region.
The HI emission in SH2 appears in one channel centered at
v=1690 kms-1 (our channel width is 41.7 kms-1) (see Fig. 7).
This is in agreement with the velocity of
the ionized gas quoted by Schweizer.
The HI mass is about
.
The distribution is unresolved at our resolution of
(
kpc).
The EELR. This extended emission line region
is located
,
or 28.6 kpc from the nucleus and has a size of
81
,
or 7.3 kpc
2.4 kpc (for our distance of 18.6 Mpc).
On the plane of the sky, it appears on a region of high surface brightness in
the loop called L1 by S80 and near a major dust lane.
The HI associated with the EELR is
and also
unresolved by our beam.
The two northern clumps are located at the outer edge of a plateau in the optical surface brightness distribution (Fig. 2, see also Fig. 7 of S80). The velocities of the HI clumps are (from north to south through west) 1857, 2000, 1750 (EELR) and 1690 (SH2), respectively, following the general rotation pattern on the central CO (next section) and suggesting that the clumps may be the left-overs of a more widespread HI structure, possibly a tidal shred. The location and velocity structure of the HI are reminiscent of other shell galaxies such as Centaurus A (Schiminovich et al. 1994) and NGC 2865 (Schiminovich et al. 1995).
12CO(1-0) | 12CO(2-1) | ||||||||||
X | Y | RA | DEC | v | ![]() |
![]() |
![]() |
v | ![]() |
![]() |
![]() |
'' | '' | '' | '' | km s-1 | km s-1 | K km s-1 | mK | km s-1 | km s-1 | K km s-1 | mK |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | |||
1999 observing | session | ||||||||||
0 | +22 | +13.5 | -17.3 | 1531 ![]() |
129 ![]() |
2.0 ![]() |
5.9 | 1550 ![]() |
88 ![]() |
1.9 ![]() |
4.9 |
0 | +44 | +27.1 | -34.7 | 1557 ![]() |
128 ![]() |
1.1 ![]() |
3.5 | 1525 ![]() |
75 ![]() |
0.9 ![]() |
3.0 |
-22 | -44 | -9.8 | +48.2 | 1809 ![]() |
167 ![]() |
4.3 ![]() |
5.0 | 1818 ![]() |
136 ![]() |
3.3 ![]() |
4.4 |
0 | -44 | -27.1 | +34.7 | 1914 ![]() |
263 ![]() |
4.5 ![]() |
5.6 | 1919 ![]() |
178 ![]() |
3.0 ![]() |
4.7 |
0 | 0 | 0 | 0 | <2.4 | 5.1 | <2.4 | 5.0 | ||||
0 | -22 | -13.5 | +17.3 | 1813 ![]() |
268 ![]() |
3.4 ![]() |
4.4 | 1859 ![]() |
186 ![]() |
2.0 ![]() |
3.0 |
-22 | -22 | +3.8 | +30.9 | 1778 ![]() |
178 ![]() |
3.5 ![]() |
4.8 | 1797 ![]() |
67 ![]() |
1.4 ![]() |
4.3 |
-22 | 0 | +17.3 | +13.5 | <2.4 | 5.1 | <2.1 | 4.2 | ||||
-22 | -66 | -23.3 | +65.6 | 1835 ![]() |
224 ![]() |
3.4 ![]() |
5.8 | <2.1 | 4.7 | ||
-22 | +22 | +30.9 | -3.8 | <3.3 | 7.1 | <2.4 | 5.4 | ||||
-44 | -44 | +7.6 | +61.8 | <2.7 | 5.9 | <2.1 | 4.2 | ||||
0 | -66 | -40.6 | +52 | <2.7 | 5.9 | <2.4 | 5.4 | ||||
-44 | -22 | +21.1 | +44.4 | <2.7 | 6.1 | <3.6 | 7.6 | ||||
0 | +66 | +40.6 | -52 | <3.3 | 7.3 | <2.4 | 5.0 | ||||
-11 | -33 | -11.6 | +32.8 | 1807 ![]() |
206 ![]() |
3.5 ![]() |
7.7 | 1802 ![]() |
202 ![]() |
3.5 ![]() |
4.6 |
-11 | +55 | +42.5 | -36.6 | <2.7 | 5.9 | <4.8 | 10.6 | ||||
-11 | +33 | +29.0 | -19.2 | <2.4 | 5.0 | <3.3 | 7.0 | ||||
2001 observing | session | ||||||||||
0.0 | 0.0 | 0 | 0 | 1634 | 343 | 3.4 | 2.1 | 1653 | 208 | 2.0 | 2.1 |
-27.7 | -35.5 | 0 | +45 | 1829 ![]() |
157 ![]() |
3.2 ![]() |
5.4 | 1795 ![]() |
48 ![]() |
1.7 ![]() |
3.4 |
-10.4 | -49.0 | -22 | +45 | 1846 ![]() |
190 ![]() |
4.3 ![]() |
6.4 | 1910 ![]() |
186 ![]() |
2.9 ![]() |
5.4 |
+7.8 | -63.2 | -45 | +45 | 1900 ![]() |
400 ![]() |
4.7 ![]() |
4.6 | <1.5 | 9 | ||
-45.0 | -21.9 | +22 | +45 | 1737 ![]() |
242 ![]() |
2.6 ![]() |
3.4 | 1697 ![]() |
270 ![]() |
2.6 ![]() |
4.8 |
-59.2 | -40.0 | +22 | +68 | 1671 ![]() |
207 ![]() |
1.4 ![]() |
4.8 | 1578 ![]() |
394 ![]() |
3.3 ![]() |
6.8 |
+3.8 | -30.9 | -22 | +22 | 1898 ![]() |
292 ![]() |
2.7 ![]() |
4 | 1967 ![]() |
41 ![]() |
1.6 ![]() |
9 |
-3.8 | +30.9 | +22 | -22 | 1514 ![]() |
73 ![]() |
1.5 ![]() |
5.3 | 1529 ![]() |
103 ![]() |
2.8 ![]() |
7 |
-13.5 | -17.3 | 0 | +22 | 1805 ![]() |
333 ![]() |
1.1 ![]() |
4.2 | <1.8 | 7.8 | ||
+21.9 | -45.0 | -45 | +22 | <1.8 | 7.1 | <4.5 | 12. | ||||
+13.5 | +17.3 | 0 | -22 | 1506 ![]() |
142 ![]() |
1.1 ![]() |
4.7 | 1500 ![]() |
327 ![]() |
2.7 ![]() |
10. |
-21.9 | +45.0 | +45 | -22 | <0.9 | 4. | <2.1 | 8. | ||||
-7.8 | +63.2 | +45 | -45 | <1.2 | 4.4 | <2.4 | 9. | ||||
-41.9 | -53.6 | 0 | +68 | <1.8 | 5.4 | <4.5 | 12. | ||||
-23.7 | -67.7 | -23 | +68 | 1963 ![]() |
334 ![]() |
3.6 ![]() |
8.7 | 1960 ![]() |
72 ![]() |
1.2 ![]() |
15. |
Combined | data and | outer HII regions | |||||||||
0.0 | 0.0 | 0 | 0 | 1559 | 189 | 2.4 | 2.1 | 1642 | 213 | 1.9 | 2.1 |
EELR | 3.2 | 3.8 | |||||||||
SH2 | 3.1 | 4.1 |
Results from the CO observations are presented in Figs. 8 to 10
and the line parameters are given in Table 3.
The observations were done by two observers in two separate
sessions, and the observing strategy was slightly different.
During the first observing session
(see top panel of Fig. 8),
we mapped the CO emission along the axis with a position angle PA
,
given by S80 as the direction of the maximum extent of the ionized gas
and probably the rotation axis of the spheroid
(perpendicular to the optical major axis).
During the second session (Fig. 8, bottom panel), the grid was in RA/DEC.
Emission is detected north-west and south-east of the nucleus as well as
toward the center.
The intensities in the NW region are about twice as strong as in the SE,
both in CO(1-0) and CO(2-1).
The CO distribution roughly follows that of the dust along the minor
axis and to the NE.
Emission was detected toward the central position.
The averaged central spectra obtained during both observing sessions
are shown in Fig. 9.
The noise level in our CO(1-0) spectrum is about 5 times lower than that in the
spectrum obtained by
Wiklind & Henkel (1989) and the bandwidth is larger, which makes it
easier to see the line structure. The integrated line intensity we infer is about three
times lower than that reported in that work, but it may agree within the not quoted
error bars of Wiklind & Henkel.
Sage & Galletta (1993) had detected CO(1-0) at two positions 30'' to the SE and
50'' to the NW of the nucleus but not toward the center. The fluxes and the limits
quoted agree with ours within the uncertainties.
The CO(1-0) line seems to be very broad (
500 kms-1 at the base)
and to consist of three
components. It is likely that the molecular gas distribution is discrete and that
the low-velocity component arises from emission from the SE that enters
the 43'' beam,
whereas the hint of emission at high velocity (
1900 kms-1) comes from the
NW concentration. All those features are, if present at all, much weaker in the
CO(2-1) spectrum taken with a narrower beam.
We calculate a mass of the H2 gas of
108
within the inner 43'' (3.9 kpc).
This estimate of mass is based upon a standard Galactic conversion factor.
The NW concentration is centered at
,
relative to the central position,
that is, at a projected distance of 3.5 kpc from the nucleus.
The emission is extended on the scale of the 22'' CO(2-1) beam,
and the measured apparent diameter at half maximum intensity
is 35'', which
correspond to a true extent of
,
or a linear size of
2.5 kpc.
This is consistent with the CO(1-0) measurements, which yield an apparent diameter
of the half peak intensity of 50'' (and a true source size of 24'').
Velocities around 1850 kms-1 are found in the NW region,
that is,
70 kms-1 higher than the systemic velocity
of NGC 1316.
A slight east-west velocity gradient is observed, with higher velocities to the
west.
From the apparent extent of the northwest CO concentration measured above and the
CO(1-0) peak intensity, we derive a mass of molecular hydrogen within that
region of
.
In the southeast, the strongest emission is detected 22''
(
2 kpc) away from the nucleus and is half as strong at 44''.
It corresponds to a H2 mass
of
.
The velocities measured there are of the order
of 1550 kms-1, that is,
230 kms-1 lower than the systemic velocity
of NGC 1316.
The CO velocity fields are displayed in Figs. 10c and 10d. They
are fairly regular. The systemic velocity of NGC 1316, however, is likely to be higher
than suggested from those interpolated contour maps. As discussed above
for the central position,
the low velocity wing around 1500 kms-1 could be due to emission
from the NW
concentration that enters the beam.
The centroid velocity of 1559 kms-1 given in Table 2 for
the central position would therefore be too low and affect the
interpolation in the contour map.
![]() |
Figure 9: CO(1-0) detection toward the center of NGC 1316. The line is broad and seems to consist of three components. The y axis shows the main-beam temperature in Kelvin. The width of the window set to subtract a baseline is shown on the x axis: 1444 to 2032 kms-1 for the CO(1-0) spectrum as indicated by the box, and 1460 to 1917 kms-1 for the CO(2-1) spectrum as indicated by the ticks. The CO(2-1) line seems to be narrower than the CO(1-0) one, which is to be expected if the molecular gas is located in the dusty regions offset from the center and located outside the narrower CO(2-1) beam. The spectra have been produced by averaging the data taken during both observing sessions. |
Open with DEXTER |
We have searched for CO emission toward the giant HII region SH2
and the extended emission line region
but have not detected any,
down to the noise levels quoted in Table 3. For SH2 as for EELR, this corresponds
to a 5
limit on the H2 mass of
107
per 43'' beam per
14.5 kms-1 channel.
While CO is clearly detected in the central region of NGC 1316, atomic hydrogen is not. The distribution of the molecular gas follows that of the dust, and its kinematics provides new insight into the dynamical history of the system. HI is detected in the outer part where no CO is found down to a similar mass limit.
![]() |
Figure 10: a) and c) Contour map of the CO(1-0) distribution and corresponding velocity field. The levels of the total intensity map go from 0.5 to 6 K km s-1 in steps of 1 K km s-1. The levels of the velocity map go from 1500 to 2000 kms-1 in steps of 50 kms-1. The circle shows the size of the beam. The black dots indicate the observed positions. The central position is marked by a cross. b) and d) Contour map of the CO(2-1) distribution and corresponding velocity field. The levels are the same as for the CO(1-0). |
Open with DEXTER |
The central part of NGC 1316
The dust patches are very visible in the B-I HST image of NGC 1316 displayed in Fig. 11c. They appear on the NW and the SE side of the nucleus, roughly aligned along the minor axis of the galaxy. The SE lane can be traced to within 1 arcsec from the center, which led Schweizer (1980) to suggest that it lies on the near side of the galaxy. The molecular gas distribution at 22'' resolution follows roughly the distribution of the dust.
Ionized gas
was found by S80 to be associated with the dust and distributed in a rotating disk
with extremal velocities kms-1 at 84'' NW and 35'' SE.
The SE side of the disk is approaching. The extreme range of CO velocities that we
measure on either side of the nucleus, 1500 to 1960 kms-1
at 32'' SE and 72'' NW, is thus smaller
than the range of velocities in the ionized disk.
Shaya et al. (1996) used HST
V and I band images of the central region
(a field of view of 70'') to analyze the three-dimensional
distribution of the dust. They derived a small amount of extinction
toward the very center, AV< 0.4 mag and AI< 0.2 mag, and
in the regions of maximum extinction.
The amount of molecular gas that we infer from the CO observations
correlates with the amount of extinction: we detect only
108
in the center and about twice as much in the dusty NW region.
The radio jet observed by
Geldzahler & Fomalont (1984)
could be deflected by interstellar gas, especially in the NW. Figure 11a
shows a superposition of our CO(2-1) map and the radio jet at 20 cm.
The direction of the central jet changes abruptly in the NW just below the
peak in the CO distribution. There is a hint of a similar bending near the
SE peak in the CO.
The most luminous galaxies in the far-infrared are mergers.
NGC 1316 was detected by IRAS in the four bands but its far-infrared luminosity
is moderate (
).
The
to M(H2) ratio, often used as an indicator of the efficiency
at which molecular gas is converted into new stars, is
5.
This is lower than the average value of
15
that we calculated for
the sample of ellipticals
observed by Wiklind et al. (1995)
(after exclusion
of the two galaxies with the highest values of
ratios,
an S0 and a galaxy larger than the beam size, for which the H2 mass must have been
underestimated).
Nevertheless, this may not mean that the star formation efficiency of
NGC 1316 is
lower
than in other ellipticals since it is known that some of the
far-infrared emission is not related to massive ionizing stars, but to dust
heated by the overall stellar radiation field.
Moreover, the Galactic CO-to-H2 conversion factor may not apply to ellipticals
or it may vary from one elliptical to another. For example, an elliptical galaxy
that has accreted a gas-rich
but metal-poor companion could be CO-poor and its H2 mass would be underestimated.
The outer HII regions
It is interesting that SH2, despite being so much smaller than EELR,
possesses about twice as much atomic gas.
Mackie & Fabbiano (1998) give an H
luminosity
for the EELR, and a mass of ionized gas
(when corrected for our adopted distance).
If the EELR were photoionized by hot stars, this quantity of
ionized gas could be maintained by
40 O-type stars. However,
Mackie and Fabbiano argued that there is no evidence of a population of
young stars associated with the EELR and suggested that the ionized
gas is produced by shocks.
Although Mackie & Fabbiano
do not quote an ionized gas mass for SH2, it appears from their H
image, which has been smoothed to a resolution of
,
that the
H
flux is of the order of
ergss-1cm-2, integrated over the
size reported
by S80. At our adopted distance, such recombination emission can be maintained
by a single O-type star. For comparison, the giant HII region NGC 604 in
the Local Group galaxy M 33 contains
of ionized gas
(Churchwell & Goss 1999) and approximately 200 O-type
stars (González Delgado & Pérez 2000). Although the
H
emission of SH2 may signify recent star-forming activity, the
population of massive stars is not large in comparison with the mass of
neutral gas in the vicinity,
.
In the central region of NGC 1316 the HI 21 cm line is not clearly
detected in emission nor is there any noticeable absorption against the
continuum emission. Our limit on HI emission in the spheroid of NGC 1316
could represent a mass of H as high as 108
.
The limit on
absorption against the central continuum source corresponds to a
column density of H
cm-2,
implying that only gas with spin temperature
K
could reasonably have been detected. The upper limit on HI emission
may at first sight seem surprising in view of the
large amounts of dust and molecular gas in the spheroid. Fornax A
is not unique in this respect: other well known elliptical galaxies
and merger remnants show CO and dust
but no detectable H in their centers
(e.g. NGC 7252, Wang et al. 1992; Hibbard et al. 1994).
In contrast there are smaller elliptical galaxies or merger
remnants that have dust, CO and plenty of H all the way in to their centers.
Examples are NGC 5128, NGC 3656, and Arp 230, all of which have disks
that are viewed nearly edge-on
(Schiminovich et al. 1994;
Eckart et al. 1990;
Balcells et al. 2001;
Sofue et al. 1993;
Schiminovich et al. 2001;
Galletta et al. 1997).
The H-deficient and H-rich systems also
differ in the size of the spheroid and presence of extended X-ray emission.
However, the detectability of HI absorption may also depend on whether
a gaseous disk is present and if so how it is inclined to the observer's
line of sight.
The CO emission from NGC 1316 probably arises in the surfaces of
molecular clouds. Any neutral atomic gas associated with such
molecular clouds may cover too little of the central continuum source to
be detectable in absorption in the 21 cm line
at the angular resolution and sensitivity of the available measurements.
The high pressure of the hot X-ray-emitting gas and the shocks associated
with the merger events may help confine the neutral gas in small clouds,
while the X-radiation keeps the more dilute gas ionized.
Finally, it is interesting to consider briefly how the spin temperature
of H in the central region of NGC 1316 might be affected by the extreme
environment there. There are radiative processes that can boost the
spin temperature of low-density H: absorption and stimulated emission
induced by the continuum radiation in the 21 cm line itself (cf. Bahcall
& Ekers 1969), and absorption and fluorescence in the
HI L
line (Deguchi & Watson 1985). These processes
will not affect the interpretation of optically thin HI emission lines,
but they could make absorption even less readily detectable.
Although Fornax A is a bright continuum source at
21 cm, the
integrated flux density of 0.24 Jy over the central
(Ekers et al. 1983; Geldzahler & Fomalont 1984) corresponds to an
average continuum brightness temperature of 565 K and yields an
absorption rate
of
s-1.
The brightness temperature in the extended central jets is lower,
150 K, which might be a lower limit on the spin temperature of
dilute gas in the center.
Pumping in the HI L
line will have a larger effect on the
spin temperature of HI. The average surface brightness of the
HI L
line can be estimated from the observed surface brightness
in the H
line,
ergs-1cm-2arcsec-2 (Kim et al. 1998), with the assumption that a fraction
of the Balmer-
emissions are followed by emission in L
.
The corresponding pumping rate
is
10
s-1,
averaged over the central arcmin of the galaxy. Since
is likely
to be of the order of unity, this excitation rate is much higher
than the mean rate of continuum pumping in the 21 cm line and is
high enough to dominate over spin-changing collisions if the
density of the HI-emitting gas is
102 cm-3. When
L
pumping controls the spin temperature, this temperature
is probably
104 K (Deguchi & Watson 1985). Thus
the L
pumping could make the HI line undetectable in absorption in
the center of NGC 1316 by ensuring that its spin temperature remains high
even where the kinetic temperature might be much lower than 104 K.
The true magnitude of the L
pumping effect is uncertain, of course,
because the distributions of dilute neutral gas and of L
photons
are not known.
There is a large amount of hot, ionized, X-ray-emitting gas in the central kpc of NGC 1316, which arises at the expense of neutral gas that might otherwise exist there. In the future, it will be valuable to obtain HI observations of greater sensitivity and higher resolution to try to trace the neutral gas and its interactions with its environment in more detail.
![]() |
Figure 11: a) CO(2-1) contour map overlaid on the 20 cm radio jet mapped by Geldzahler & Fomalont (1984). b) CO(2-1) contour map overlaid on Schweizer (1980)'s sketch of the dust lanes (dots). c) HST B-I image of the central region of NGC 1316 (press release STScI-PR99-06, Grillmair et al. 1999). |
Open with DEXTER |
The complex structure of NGC 1316 (five tails or loops of varying morphology)
argues for a rich dynamical history and probably a succession of merging events.
The presence of shells around ellipticals and minor-axis dusty disks
can be interpreted in two different ways: as tracers of the accretion of
a small companion onto an already existing elliptical galaxy (e.g. Quinn 1984),
or as the returning of tidal material during the merger of two disk galaxies
(e.g., Hernquist & Spergel 1992;
Balcells 1997).
There is observational evidence that NGC 1316 formed in a major merger 3 Gyr ago
and has been accreting material more recently, as discussed below.
NGC 1316 has a relaxed morphological appearance within 3',
which indicates that, if the galaxy is the result of a major merger between two disk
galaxies, the event
must have occurred more than
years ago.
The population of globular clusters constrains the age of the system.
Among 24 star clusters associated with NGC 1316, Goudfrooij et al. (2001)
identified four exceptionally luminous ones (
MV < 12.3) with near-solar metallicities,
whose age of
Gyr
could be determined from comparison with single-burst population models.
Those intermediate-age clusters would have formed during a merger of two disk galaxies
that resulted in the formation of the elliptical whereas the older (
14 Gyr old)
metal-poor clusters are likely to have belonged to the progenitor disk galaxies.
As a remnant of the merger of two disk galaxies, NGC 1316 would be at a more advanced
stage than for instance NCC 7252, which is believed to be 0.5-1 Gyr old and still
displays extended tidal tails (Schweizer & Seitzer
1998).
Optical spectra of the inner 3'' of NGC 1316 have revealed a stellar population
of age 2 Gyr (Kuntschner 2000). Those stars would therefore have
formed from the molecular gas that fell toward the center following the merger event
and would have contributed to the depletion of molecular gas observed in the central
region today. Strong nuclear H
absorption has been found
(Schweizer 1981; Bosma et al. 1985).
It might indicate an important component of A-type stars which are part of a
post-starburst population.
The X-ray luminosity increases as a result of gas being shock heated or evaporated by
electrons in the already existing hot plasma
(e.g., Goudfrooij & Trinchieri 1998)
and mass loss from evolving giant stars and supernovae Ia
(e.g., Brighenti & Mathews 1999).
Using ROSAT, Kim et al. (1998)
detected
(at D=18.6 Mpc) of hot
(
107 K) gas, which is small compared to the amount of hot gas found in some
early-type galaxies of comparable optical luminosity. However,
Kim et al. (1998) point out that the amount of hot gas
is comparable with what would be expected from stellar evolution in
1 Gyr. A comparable amount of hot gas
is coincident with the tidal tails in the outskirts
of NGC 1316 (Mackie & Fabbiano 1998) and may be due to more recent events.
The presence of loops and ripples as well as of a disk of ionized gas
is strong evidence in favor of the infall of a small galaxy onto NGC 1316
less than 1 Gyr ago (Schweizer 1980).
In particular, the giant loop called L1 by Schweizer is similar to features seen
in simulations of shell galaxies (e.g., Quinn 1984;
Hernquist & Quinn 1988) and
could be interpreted
as due to the infall of a dynamically colder galaxy 0.5 Gyr ago.
Mackie & Fabbiano (1998) have argued from the temperature of the hot gas detected in the
L1 region that L1 could be the
remnant of a small gas-rich galaxy on a high-velocity encounter
(
380 kms-1) with NGC 1316, and its origin would therefore be different from that
of the other loops which formed in the aftermath of the major, older merger.
The extended emission-line region EELR is found in the area of highest surface
brightness of L1.
The fact that we don't detect any CO from the EELR is not surprising if the infalling galaxy
was metal-poor.
The low level of detectable HI emission from that region may be due to the fact that the
most of the atomic hydrogen has been stripped, or photo-ionized by the newly formed OB associations.
The NW concentration could as well be a recently accreted gas-rich galaxy.
It has a linear dimension of about 2.5 kpc and contains
about
of molecular hydrogen.
This H2 mass is comparable to that of NGC 1317,
the small (
15 kpc diameter) spiral companion of NGC 1316
(Horellou et al. 1995).
This gas concentration is located 40'' (3.5 kpc) from the nucleus,
slightly to the north of the minor axis. The velocity of the molecular gas,
kms-1, is
70 kms-1 higher than the systemic
velocity given by Arnaboldi et al. (1998) and agrees exactly
with the stellar velocity derived from absorption line spectra in that region
(Fig. 3 of Arnaboldi et al. 1998). Interestingly, beyond
30'' from the nucleus along the minor axis, the stellar radial velocities increase
on both sides, producing a U-shaped velocity curve, and the velocity dispersion
profile is irregular. Arnaboldi et al. argue that this kinematical feature is
caused by infalling material within one dynamical time-scale at 40'' (i.e.
less than
years ago), since such features are not expected
to last longer.
It is interesting to compare NGC 1316 with NGC 5128 (Centaurus A).
NGC 5128 lacks tidal tails but it
has a more prominent dust lane and a better defined system
of shells in both of which atomic and molecular gas have been detected
(Eckart et al. 1990;
Schiminovich et al. 1994;
Charmandaris et al.
2000).
NGC 5128 has comparable amounts of HI and H2 in the region of the dust lane
(
)
whereas NGC 1316 has at least four times more
molecular gas than HI.
This difference in HI content may be due to the fact that the HI in NGC 1316
is more affected
by the X-ray emission, since
NGC 1316 is brighter in X-ray than NGC 5128
(factor of 8 in the 0.1-2.4 keV band, ROSAT observations, and
factor of 1.4 in the 0.2-4 keV band, Einstein IPC observations).
Both galaxies have
a double-lobe radio source. The radio distribution, however, has a
different orientation: in Centaurus A it is perpendicular to the dust lane
whereas in Fornax A it is not. Nuclear effects may therefore affect the gas and dust properties
more clearly in Fornax A than in Centaurus A.
In the outer parts of NGC 5128, HI is found in the shells, with CO detected in
the highest column density regions with M(HI)/M(H
.
In NGC 1316, we could not detect any CO from the bright outer HII regions, but
our limit on the H2 mass is comparable to the HI mass.
The possible chain of HI clumps in the outer part of NGC 1316
could be the highest peaks in a continuous shell. More sensitive observations are
required to confirm this.
The HST image of NGC 1316 displays beautifully the distribution of the dust
(Fig. 11c).
Northwest of the nucleus it roughly has a T-shape with a radial extent of
about 40'' and an azimuthal extent of almost 90.
CO is detected in the area of both the radial and the azimuthal dust patches.
Grillmair et al. (1999) point out the similarity of the radial
structure to the photo-evaporating columns in the Eagle Nebula M 16
(Hester et al. 1996), though on a much larger scale, and suggest
that it may be due to Rayleigh-Taylor instabilities or subject to erosion
by particles or radiation from the nucleus of the galaxy.
The azimuthal feature may be a dusty shell.
The azimuthal configuration is expected from recently accreted material
oscillating between two apocenters and forming shells (e.g. Hernquist & Quinn
1988). Gaseous shells could form in the same phase-wrapping process
if the gas is clumpy and not very dissipative
(Charmandaris & Combes 2000).
A detailed dynamical model of the gaseous features in NGC 1316 based on the accretion
of a small companion will be the subject of another paper.
Our new observations of molecular gas (as traced by 12CO(1-0) and 12CO(2-1) line emission) toward the peculiar elliptical galaxy NGC 1316 and a new analysis of H Idata have shown that:
Acknowledgements
We are grateful to the SEST staff for their help and support during the observations and to the VLA analysts for their prompt replies and useful help. This work was in part supported by an NSF grant to Columbia University. We have used the tool DEXTER provided by NASA Astrophysics Data System to extract data from Fig. 11a of Schweizer (1980) and from Fig. 2 of Geldzahler & Fomalont (1984). Those figures and the radio continuum map of Fomalont et al. (1989) shown in Fig. 1 were used with the permission of the authors. We have used images from the Digitized Sky Surveys. The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Government grant NAG W-2166. The images of these surveys are based on photographic data obtained using the Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. We are grateful to the referee for a careful reading of the manuscript and useful comments.