A&A 440, 893-899 (2005)
DOI: 10.1051/0004-6361:20052816
C. Henkel1 - N. Jethava1 - A. Kraus1 - K. M. Menten1 - C. L. Carilli2 - M. Grasshoff3 - D. Lubowich4 - M. J. Reid5
1 - Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
2 -
National Radio Astronomy Observatory, Socorro, New Mexico 87801, USA
3 -
Visiting the Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
4 -
Department of Physics and Astronomy, Hofstra University, Hempstead, NY 11549, USA
5 -
Harvard-Smithsonian Center for Astrophysics, 60 Garden St., MS 42, Cambridge, MA 02138, USA
Received 3 February 2005 / Accepted 1 June 2005
Abstract
Using the Effelsberg 100-m telescope, absorption in the (
J,K) = (1,1), (2,2) and (3,3) inversion lines of ammonia
(NH3) was detected at a redshift of
z = 0.6847 toward the gravitational lens system B0218+357. The
cm
absorption peaks at 0.5-1.0% of the continuum level and appears to cover a smaller fraction of the radio continuum background
than lines at millimeter wavelengths. Measured intensities are consistent with a rotation temperature of
35 K, corresponding
to a kinetic temperature of
55 K. The column density toward the core of image A then becomes N(NH
cm-2 and fractional abundance and gas density are of order X(NH
10-8 and n(H
cm-3, respectively. Upper limits are reported for the (2,1) and (4,4) lines of NH3 and for transitions
of the SO, DCN, OCS, SiO, C3N, H2CO, SiC2, HC3N, HC5N, and CH3OH molecules. These limits and the kinetic
temperature indicate that the absorption lines are not arising from a cold dark cloud but from a warm, diffuse, predominantly
molecular medium. The physical parameters of the absorbing molecular complex, seen at a projected distance of
2 kpc to the
center of the lensing galaxy, are quite peculiar when compared with the properties of clouds in the Galaxy or in nearby
extragalactic systems.
Key words: galaxies: abundances - galaxies: ISM - quasars: individual: B0218+357 - quasars: absorption lines - radio lines: galaxies
With an image separation of 334 mas B0218+357 is one of the most compact gravitational lens systems known to date (O'Dea et al. 1992; Patnaik et al. 1993). The lensed source, possibly a BL Lac object (Kemball et al. 2001),
is located at a redshift of
(Cohen et al. 2003) and shows a complex radio structure with two dominating compact sources, A and B, and an Einstein ring (e.g. Biggs et al. 2001a). A time delay of order 10 days has been measured
between the two compact sources (Corbett et al. 1996; Biggs et al. 1999, 2001b; Cohen et al. 2000). Absorption from the lensing galaxy, a face-on spiral likely of type Sa/Sab (York et al. 2005), is observed
at a redshift of z=0.68466 at optical wavelengths (Brown et al. 1993; Stickel & Kühr 1993) and, remarkably, in the
21 cm line of H I (Carilli et al. 1993) and in numerous radio- and millimeter-wavelength lines from a variety of molecules. The absorption is observed against the compact component A which is, at
radio waves,
3 times stronger than component B (e.g. Patnaik et al. 1993). Detected molecular species are CO, HCN, and HCO+ (Wiklind & Combes 1995), H2CO (Menten & Reid 1996), H2O (Combes & Wiklind
1997), CS (Combes et al. 1997), OH (Kanekar et al. 2003) and tentatively also LiH (Combes & Wiklind 1998). The molecular spectra, combined with 21 cm H I absorption profiles, have been used to constrain
the temporal evolution of the fine structure constant (e.g. Carilli et al. 2000; Murphy et al. 2001; Kanekar & Chengalur 2004). The nature of the cloud, however, is poorly understood.
One of the less well known parameters of extragalactic molecular clouds is the kinetic gas temperature. With area filling factors
that are not well known, a thermalized tracer like CO cannot be used to determine
.
Better thermometers are symmetric
top molecules, where relative level populations are determined predominantly by collisions. The two most prominent such molecules
are NH3 and CH3CN. While extragalactic CH3CN was so far only detected in NGC 253 (Mauersberger et al. 1991),
various NH3 "inversion'' lines, from the (
J,K) = (1,1) up to the (6,6) and even the (9,9) line, have now been observed in the
nuclear regions of nearby (
)
galaxies (e.g. Martin & Ho 1986; Henkel et al. 2000; Takano et al.
2000, 2002; Weiß et al. 2001; Mauersberger et al. 2003).
B0218+357 provides a unique view onto a molecular cloud seen at about half a Hubble time in the past (
-cosmology with H0=71 km s-1 Mpc-1,
and
;
Spergel et al. 2003). In an attempt to constrain the physical properties of this cloud at a luminosity distance of
3 Gpc, we have searched for NH3 and
other molecular species.
Table 1: Observational parameters.
Table 2:
rms line to continuum ratios in units of 10
]
(Col. 4) and corresponding channel spacings (Col. 5).
All observations were made with the 100-m telescope at Effelsberg/Germany. In August 2001 and June 2002 we searched, employing a single channel
cm HEMT receiver, for ammonia (NH3), sulfur monoxide (SO), methanol (CH3OH), the SiC2 radical, cyanodiacetylene (HC5N) and carbonylsulfide (OCS). The measurements were carried out in a position
switching mode. In December 2001 we used a dual channel 1.3 cm HEMT receiver to search for cyanoacetylene (HC3N) and silicon monoxide (SiO) in its ground vibrational state (see Tables 1 and 2). These measurements were made in a dual beam switching mode
with a beam throw of 2' and a switching frequency of
1 Hz. In January 2002, we also employed a single-channel 7 mm HEMT receiver to search for formaldehyde (H2CO) and deuterated hydrogen cyanide (DCN). These observations were made in a position switching mode. Most recently, in June 2003, a two channel 5 cm HEMT receiver was used to search for the N=1-0
J=3/2-1/2
transition of the C3N radical.
Frequencies, beamwidths, system temperatures and aperture efficiencies are given in Table 1. For all measurements we employed an "AK 90'' autocorrelator with eight spectrometers, using bandwidths of 40 MHz and 512 channels (20 MHz and 1024 channels and 80 MHz and 256 channels at 5 cm and 7 mm, respectively). Calibration was obtained from measurements of NGC 7027 (Ott et al. 1994). Pointing corrections could be obtained toward B0218+357 itself and were accurate to better than 10''.
![]() |
Figure 1: NH3 lines with a velocity scale relative to z=0.68466, observed towards B0218+357. Channel widths are 3.3 km s-1. The ordinate is in units of the radio continuum flux. In the case of the (1,1) line a two component fit has been added accounting for NH3 hyperfine splitting (Kukolich 1967). The fit demonstrates that the line profile is consistent with two optically thin velocity features. Limited signal-to-noise ratios do not constrain the fit sufficiently to justify an inclusion of its parameters in Table 3. |
| Open with DEXTER | |
Figure 1 shows the measured NH3 profiles. The detected inversion lines, with the ordinate displaying absorption in
units of 1% of the continuum flux, are characterized by a prominent narrow component with a width of a few km s-1,
centered at slightly positive velocities, and a wider and weaker component that is centered at slightly negative velocities
with respect to a redshift of z=0.68466. Tables 2 and 3 display noise levels for all measured spectra and line parameters
for the detected transitions. The weak but wide velocity component has a 2-3 times higher integrated intensity than the more
prominent narrow one. At its peak, the narrow component reaches 0.5-1.0% of the total continuum flux density, while the
broad component peaks at
0.3% in the (1,1) and (2,2) lines and at an even lower level in the (3,3) transition. The
narrow component must be optically thin; otherwise the hyperfine satellites would be visible (for the (1,1) line at
8
and
19 km s-1 with respect to the main feature). Due to a lower signal-to-noise ratio, the optical depth of the wide
component is poorly constrained (see also Sect. 4.1).
At
and 1.3 cm, continuum flux densities were 1.0-1.2 Jy, in good agreement with Patnaik et al. (1993)
and Menten & Reid (1996). At 6 GHz we find 1.45 Jy with an error of
10%, while no flux density was determined
at 7 mm.
Figure 2 shows a Boltzmann plot (rotation diagram) including the four measured "metastable'' (J=K) inversion lines of ammonia. Applying
Table 3: NH3 line parametersa.
![]() |
Figure 2:
Boltzmann plot (rotation diagram) of the normalized NH3 column density (
|
| Open with DEXTER | |
Prior to this study, Menten & Reid (1996) had searched for NH3 absorption towards B0218+357. Their upper limits are higher than the strength of the absorption features shown in Fig. 1. It is tempting to identify the two main features of
the NH3 lines as absorption against the two dominant continuum sources A and B. However, this is extremely unlikely for several reasons. VLBI observations of H I at dm- and of H2CO at cm-wavelengths have shown that the absorption is confined to image A (Carilli et al. 2000; Menten & Reid 1996). While at radio wavelengths
image A is three times as bright as B, image B is brighter at optical wavelengths. This indicates that image A is obscured by dust. Heavy absorption of A is further supported by a large differential rotation measure between A and B (Patnaik et al. 1993) and by the different separations between images A and B at radio and optical wavelengths (
versus
mas, respectively; York et al. 2005). Apparently, radio and optical images of A do not coincide. We further note that images A and B are on opposite sides of the lensing galaxy. Assuming a rotation velocity of order
150 km s-1,
a difference in radial velocity of
5 km s-1 between images A and B would then require an inclination
for the lensing galaxy. While the galaxy is clearly face-on (York et al. 2005), the optical image does not require an inclination that close to zero.
We conclude that the two velocity components are associated with image A that shows a core-jet morphology (Patnaik et al. 1995; Porcas & Patnaik 1995; Biggs et al. 2003). While the stronger narrower velocity component has a linewidth that is well within the range observed in galactic giant molecular clouds (e.g. Combes 1991), the linewidth of the broad component is much larger. Could this be a consequence of hyperfine splitting of the ammonia lines? Or a consequence of differential galactic rotation?
The NH3 hyperfine satellites (for a (
J,K)=(1,1) spectrum, see e.g. Ho & Townes 1983) are symmetrically bracketing
the main feature and are covering
40 km s-1. Thus the partially blended satellite features from both velocity
components can explain the observed linewidths of
50 km s-1. For a clear separation of individual satellite lines
and velocity components, data with higher signal-to-noise ratios would be required.
The jet extends from the core of A over
6 mas toward the center of the lensing foreground galaxy, located at an angular distance
of almost 300 mas (Patnaik et al. 1995; Porcas & Patnaik 1995; Biggs et al. 2003; Wucknitz et al.
2004). Position angles are
68
for the jet (Biggs et al. 2003) and
for the center (Wucknitz et al. 2004). Therefore no significant velocity gradient is expected from the foreground galaxy.
Its low inclination further reduces any such gradient. A wide component with an FWHP linewidth of order 50 km s-1 is not seen
in other molecular lines at cm- and mm-wavelengths (e.g. Wiklind & Combes 1995; Menten & Reid 1996). Thus
hyperfine splitting is likely the main cause for the linewidths observed. Note, that there are indications for the same asymmetry,
i.e. a stronger feature at slightly positive and a weaker feature at slightly negative velocities, in the CO J=2-1 and 3-2,
HCN 2-1 and HCO+ 2-1 lines (Wiklind & Combes 1995; Combes & Wiklind 1997, 1998).
To summarize: While the relatively strong narrow absorption component at slightly positive velocities (see Fig. 1) must be optically thin, the optical depth of the wider and weaker component remains undetermined. Total linewidths are greatly affected by the spacing of the hyperfine features. Gaussian fits to the wide component (Table 3) include the hyperfine satellites of both velocity components.
Molecular lines at mm-wavelengths can show deep absorption (Combes & Wiklind 1995; Wiklind & Combes
1995; Combes & Wiklind 1997). Lines at dm- and cm-wavelengths are generally weaker with respect to the
continuum. This even holds for
cm H I and 18 cm OH, that show wider absorption than ammonia and the mm-wave
lines (Carilli et al. 1993; Kanekar et al. 2003) and that may thus arise from a larger volume. The peak line
to total continuum ratios of the 18 cm OH main lines are, however, similar to those seen in ammonia.
Since the continuum sources are compact and since the measured NH3 linewidths are small by extragalactic standards, the
absorption arises from an area that must be considerably smaller than that studied in emission with single-dish telescopes towards
nearby galaxies. According to Patnaik et al. (1995),
70% of the total 22 GHz flux is associated with image A.
Thus minimum optical depths and source covering factors (
)
of the ammonia lines toward image A are a factor of 1.4 higher
than those suggested for the entire continuum flux by Fig. 1. They reach, in the (1,1) line,
0.014.
Sensitive maps of image A show a source of size
mas2 that is edge-brightened on its south-western side and
tangentially stretched (Biggs et al. 2003). At 15 GHz, the south-eastern region containing the core (components A1
and A2) and the north-western region containing the jet exhibit
% and
% of the observed flux density, respectively
(Patnaik et al. 1995). Although higher frequency VLBI measurements are missing, it would be no surprise if the
contribution of the jet in image A is negligible at 100-150 GHz, thus leaving only the lensed core as the background continuum
source. Assuming for the core component a brightness temperature
and a flux density
not drastically varying
with frequency in the 10-150 GHz frequency range (see e.g. Blandford & Königl 1979), its solid angle would
vary like
(see also Lobanov 1998). According to NED
the mm-wave flux of B0218+357 is about half of that at cm-wavelengths, compatible with a fading jet and a flat
spectrum core.
Adopting this scenario and assuming that the absorption arises exclusively from the core as it is indicated by some of the
deep mm-wave absorption lines (CO should arise from a larger volume than NH3), we obtain for the peak of the ammonia (1,1)
line a minimum optical depth and covering factor of
.
At 150 GHz, this would yield for a cloud
centered near the peaks of image A
,
depending on whether the exponent in the frequency dependence
of the solid angle
is -1 or -2.
High source covering factors at mm-wavelengths can be achieved, if the molecular absorber is elongated along the same position angle
as the continuum, i.e. along a path with roughly constant galactocentric radius, hypothetically defining a spiral arm (see also
Sect. 5). Elongated filaments are common in the Galaxy and are exemplified by the morphology of the Orion giant molecular cloud
or the connection between this GMC and the Monoceros GMC (Maddalena & Thaddeus 1985), reaching linear sizes in
excess of 100 pc. This is more than the few 10 pc that are needed here. Observed
values near unity then suggest
that actual source covering factors (
and
)
are 1-6 times higher than the lower limits estimated above,
not providing an additional constraint for the NH3 source covering factor.
To verify these estimates and to further elucidate the morphology of image A and its foreground molecular cloud at high frequencies, mm-wave VLBI observations are highly desirable.
So far, the nature of the molecular cloud has been an enigma. It was proposed to be either a diffuse or a dark cloud (Combes & Wiklind 1995; Wiklind & Combes 1995; Menten & Reid 1996; Kanekar et al. 2003). The determination
of an NH3 rotation temperature, which is a lower limit to the kinetic temperature, allows us to discriminate between these two scenarios. Provided total gas densities are less than
cm-3 and using the Large Velocity Gradient
(LVG) model of Schilke (1989) with collision rates from Danby et al. (1988), the kinetic temperature is 50-60 K.
For even higher densities the kinetic temperature would be lower but would still be >35 K. Since frequencies of the ammonia lines are quite similar, the resulting temperatures are not significantly affected by potential differences in the morphology of the
background continuum. We therefore conclude that the NH3 absorption does not arise from a cool dark cloud. Since the HCO+ absorption lines are not fully saturated (Wiklind & Combes 1995) and many molecules have not been detected
at levels of <1% of the total continuum flux density (Table 2), the diffuse cloud scenario is preferrable.
The high kinetic temperature of the gas implies that our non-detection of the DCN J=1-0 line does not provide significant
constraints to the cosmic D/H ratio. Adopting a width of the line of order 5 km s-1, we deduce a 5
sensitivity
of 5% of the continuum level (Table 2). Assuming an optical depth of 5 for the opacity of the HCN J=1-0 main line, this
implies that we would have seen a DCN line only if
.
A ratio of 10-2 is, however, not reached at
K (Gerin et al. 1999) and [DCN]/[HCN] ratios at higher temperatures should be even smaller.
Having estimated the kinetic temperature, we can use measurements of molecular excitation to also obtain the density of the gas
toward B0218+357. From the microwave background we obtain a minimum of
(
1+z) = 4.6 K. Combes et al.
(1997) find
K from CS J=1-0 and 3-2. The detection of the 557 GHz and the
non-detection of the 183 GHz line of water vapor (these are rest frequencies) indicate
<30 K (Combes & Wiklind
1997). 13CO J=2-1 and 4-3 spectra suggest
K (Gerin et al. 1997). With
K and n(H
cm-3, LVG calculations for CS (for details, see e.g. Wang et al.
2004) yield an excitation far above that observed. Thus
K. With this kinetic temperature and
K,
cm-3. This is a firm upper limit. If the excitation is close
to 6 K or if the stronger CS line (J=1-0) is optically thick, densities are smaller (
cm-3 in
the former case), while constraints from H2O are less stringent. Again using LVG calculations, we find that a density of
cm-2 is consistent with 13CO excitation temperatures
20 K, if optical depths are smaller
than unity. Since Combes & Wiklind (1995) find that the 13CO J=2-1 line is saturated, however, the
density will be even lower. The upper density limits derived here are also consistent with the upper limit to the NH3 (2,1)
line (Fig. 1 and Table 2), that would only be detectable in a cloud with much higher density (see Mauersberger et al.
1985).
Adopting
K and, less certain,
cm-3 as best guesses, we can also
estimate the excitation temperature within an NH3 inversion doublet:
K (Schilke 1989, his Fig. 6.15).
This can be used to derive total NH3 column densities. Adding up column densities in the three detected lines gives N(NH
cm-2. This is a minimum value. With Eq. (A15) of Ungerechts et al. (1986) and rotation
temperatures of
T12=35-40 K, we obtain
cm-2. The error lies mainly in the assumed density of
the gas. n(H
2)=104 cm-3 implies an excitation temperature of 8 K and a column density of 5 and
cm-2, respectively. So far all given column densities are averaged over the entire continuum source.
Averaging over the core of image A, column densities have to be multiplied by a factor of 2.3.
Typical galactic clouds have fractional abundances of X(NH
3) = 10-8-10-7 (e.g. Hotzel et al. 2004).
Averaged over the entire continuum source, this yields N(H
cm-2 or
cm-2 for the core of image A. The lower limit is too small to allow us to detect
a variety of molecular lines. Therefore the fractional abundance is likely close to 10-8. Our column density range is
slightly lower than what has been proposed by Menten & Reid (1996) on the basis of another molecular cm-wave
transition. Adopting 55 K also for the H I spin temperature, the total background continuum averaged H I column density
becomes
cm-2 (see Carilli et al. 1993). Accounting for an image A core source covering
factor, the column density becomes
cm-2. Thus the line of sight must be mostly molecular.
Column densities averaged over the extent of the continuum source are lower at cm-wavelengths (see Menten & Reid 1996)
then those derived by the higher frequency studies of Gerin et al. (1997),
cm-2, and
Combes & Wiklind (1995),
cm-2. This is expected for a molecular cloud located near the
peak of image A and covering most of the mm-wave but only a fraction of the cm-wave continuum of B0218+357. In spite of this effect
the value proposed by Combes & Wiklind (1995), N(H
cm-2, must be too large. With
such an enormous column density many of the species listed in Table 2 would have been detected. The 18O/17O ratio of
15 suggested by Combes & Wiklind (1995), that is in part responsible for the high column density, is far outside the
range observed so far in interstellar clouds. With
in the LMC (Heikkilä et al. 1998),
4.1 in the local interstellar medium (Wouterloot et al. 2005), 5.5 in the solar system (e.g. Wilson & Rood
1994), and
6.5 in nearby starburst galaxies (Harrison et al. 1999; Wang et al. 2004), new
measurements of rare CO isotopomers would thus be desirable.
To put our data into context with the much more extensively studied interstellar medium of the local universe is not easy. The
kinetic temperature we derive,
55 K, is much higher than the canonical number of 10 K assumed for quiescent Giant
Molecular Cloud gas in the Galaxy. Clouds forming massive stars are warmer and may reach 55 K, but densities tend to be larger
than those estimated in Sect. 4.3 for B0218+357.
The central region of the irregular starburst galaxy M 82 may provide the best agreement with the molecular parameters of
B0218+357. In M 82, NH3 observations indicate
K and
K (Weiß et al.
2001), while observations of CO reveal densities of several 103 cm-3 (e.g. Mao et al. 2000).
The interstellar medium of M 82 is characerized by widespread PDRs (photon dominated regions/photodissociation regions),
where ammonia, being readily destroyed by UV photons, mostly resides in the few remaining well shielded relatively cool
cloud cores. Fractional ammonia abundances are of order
and thus much lower than in B0218+357.
In the inner few 100 pc of nearby spiral galaxies, detected in NH3, temperatures are significantly higher
than in B0218+357 (
K;
K; Mauersberger et al. 2003). This also
holds for clouds within a few 100 pc of the center of our Galaxy (Güsten et al. 1985; Mauersberger et al.
1986; Hüttemeister et al. 1993).
In the Galaxy, few measurements of NH3 absorption from diffuse clouds against radio continuum sources have been reported and
neither this molecule's abundance (if detected) nor its rotation temperature are known in such environments (Nash 1990).
It would thus be easiest to reconcile our NH3 results with a known scenario if we assumed that the B0218+357 absorption line
of sight crossed the central region of the lensing galaxy. This, however, is not the case. Wucknitz et al. (2004) model the
B0218+357 system using the L ENSC LEAN algorithm and find the lens to be centered at a position which is
(
,
mas offset from the A image and (-55,-10) mas from the B image, with
uncertainties of order 5 mas in each coordinate. The best fit size of the lens is
mas with a position angle of
east of north (Wucknitz et al. 2004, Table 4). This implies that image A is at a projected distance
of
2.0 kpc (
)
from the center of the lensing mass distribution, while image B is offset by only 0.4 kpc
(
). Since the starburst in M 82 is, as most starbursts, confined to the nuclear few 100 pc and since its
fractional NH3 abundance is much smaller than that in B0218+357, we do not find, neither in the Galaxy nor in the nearby
extragalactic universe, any analog to the line-of-sight toward image A in B0218+357.
The gravitational lens system B0218+357 allows us to study the details of a molecular cloud observed at about half a Hubble time in the past. Our analysis of the Effelsberg data reveals the following:
Acknowledgements
We wish to thank R. Porcas, C. Böttner, M. Kadler and T. Krichbaum for useful discussions and an anonymous referee for critically reading the text and making a number of important suggestions.