A Determination of the HDO/H2O Ratio in Comet C/1995 O1 (Hale-Bopp)

See allHide authors and affiliations

Science  06 Feb 1998:
Vol. 279, Issue 5352, pp. 842-844
DOI: 10.1126/science.279.5352.842


Deuterated water (HDO) was detected in comet C/1995 O1 (Hale-Bopp) with the use of the James Clerk Maxwell Telescope on Mauna Kea, Hawaii. The inferred D/H ratio in Hale-Bopp's water is (3.3 ± 0.8) × 10−4. This result is consistent with in situ measurements of comet P/Halley and the value found in C/1996 B2 (Hyakutake). This D/H ratio, higher than that in terrestrial water and more than 10 times the value for protosolar H2, implies that comets cannot be the only source for the oceans on Earth.

It is assumed that comets contain unprocessed matter from the earliest history of the solar system. To pursue this subject further, we measured the abundance of HDO in the coma of comet Hale-Bopp using the James Clerk Maxwell Telescope (JCMT) on Mauna Kea. The resulting D/H ratio helps define the processes leading to comet formation and trace the contribution of unmodified ice grains from the original interstellar cloud to various reservoirs of H2O throughout the solar system.

The HDO ground-state 101-000 transition at 464.92452 GHz was detected on 4.9 April 1997 UT (universal time) (Fig.1), about 5 days after Hale-Bopp passed perihelion (1). At the time of the observations, the comet was 0.917 astronomical units (AU) from the sun and 1.385 AU from Earth. We used the C2 single-channel superconductor-isolator-superconductor (SIS) receiver in a double-sideband mode, which convolves the two sidebands into one spectrum. Our setup enabled simultaneous measurement of two high rotational transitions in methanol (CH3OH) and the HDO line without overlap. In the upper sideband, the CH3OH 92-91 A± transition at 464.835 GHz and the HDO transition yielded line-integrated antenna temperatures ∫T A dv = 2.37 ± 0.11 and 0.64 ± 0.11 K km s−1, respectively. The CH3OH 112-111 A± line at 457.006 GHz was detected in the lower sideband and had an integrated line area of ∫T A dv = 1.71 ± 0.15 K km s−1 (2). At the beginning of the HDO observation, we took a SCUBA image of Hale-Bopp at 850 μm (3). Within the typical root-mean-square mechanical stability of the telescope of ±1.5 arc sec, this map allowed us to check our pointing. The stronger of the two CH3OH lines was then used to monitor our tracking of Hale-Bopp. To achieve continuous sky cancellation, we nodded the secondary mirror at a rate of 1 Hz with a chop throw of 180 arc sec in the azimuthal direction. During its peak-activity, Hale-Bopp was a daytime object, reaching zenith around noon. The solar elongation was 41°. JCMT is the only telescope of its class that is protected against direct solar irradiation with a membrane of woven polytetrafluoroethylene (Gore-Tex), which shades the antenna surface. The membrane transmits over 80% of the incident submillimeter radiation but reflects most of the solar light at shorter wavelengths. We obtained the HDO spectrum during a period of exceptionally dry weather. The zenith opacity at 225 GHz was ≤0.06, which corresponds to a precipitable water column of about 1 mm.

Figure 1

Spectrum of Hale-Bopp showing the HDO 101-000 line together with two methanol transitions (120 min, on and off). The CH3OH 112-111 A± transition originates in the lower sideband (LSB) (top frequency scale); all other lines belong to the upper sideband (USB). Because of different calibrations of the two sidebands, the intensity of the line in the LSB has to be multiplied by a factor of 1.41 (2). Each bin covers 313 kHz and is an average of two adjacent channels. Antenna temperatures have been calibrated against internal standards and pointing sources W3(OH) and N7538 IRS1. The velocity scale for the HDO line is relative to the rest frame of the cometary nucleus. The double-peak structure of the lines is caused by the velocity distribution in the comet.

To obtain an accurate D/H ratio in water, one would ideally want to compare HDO with an optically thin line of the rare isotopic molecules H2 18O or H2 17O using the same telescope; however, all of the stronger transitions of H2 18O and H2 17O are inaccessible to JCMT. Hence, we compared our observations with production rates measured elsewhere at about the same time. For the first week of April, a water production rateQ(H2O) = (1.0 ± 0.2) × 1031s−1 was derived from OH observations at the Nançay radio telescope (4). An inspection of the spatial distribution of H2O (5) and ground-based OH data in the near ultraviolet (6), however, revealed that derived production rates may depend on aperture size. Follow-up measurements of HDO on 5 April indicated that the absolute calibration of the C2 receiver was not fully reliable at the time of the HDO experiment. As a check, we compared the CH3OH production rate derived from the two CH3OH lines in the HDO spectrum of 4 April (Fig. 1) with those deduced from Caltech Submillimeter Observatory (CSO, Mauna Kea) and Institut de Radio Astronomie Millimétrique (IRAM, Spain) observations during the same period of time (4). Measurements at IRAM and CSO yielded Q(CH3OH) = 2.3 × 1029 s−1, a factor 1.53 higher than the value of 1.5 × 1029 s−1 derived from the JCMT data (7). Accordingly, we adjusted the C2 temperature scale by a factor of 1.53 (8). In the millimeter and submillimeter spectral range, numerous rotational lines of CH3OH can be used to probe the gas temperature in the coma. At the time we measured HDO, the beam-averaged kinetic temperature of Hale-Bopp reached ∼120 K (4, 9). This temperature was measured at IRAM and CSO and is consistent with the temperature derived from the two CH3OH lines in the HDO spectrum (Fig. 1).

The D/H ratio for different excitation models and a gas temperature of 120 K were determined (Table 1). All models assumed a Gaussian beam shape with a full width at half power of 10.3 ± 0.5 arc sec at 464.9 GHz. In model 1, we assumed a thermal distribution of the rotational levels at a constant equilibrium temperature. A Haser (10) model was used to describe the gas distribution in the coma. For HDO, we adopted a total photodissociation lifetime of 7.7 × 104 s at 1 AU (11) and an expansion velocity of 1.1 km s−1, in agreement with the line shapes. Model 1 is the simplest model and is discussed here because it is widely used and yields an upper limit and rough estimate for Q(HDO). Model 2 considers excitation of rotational levels by vibrational levels in the infrared and by collisions with neutral water molecules (12). Unlike in model 1, no equilibrium is assumed. Instead, the model accounts for the changing physical environment as the coma expands radially outward. In the first few thousand kilometers from the nucleus, collisions determine the population of the rotational levels, and at larger distances, radiative processes prevail. Model 2 depends on the number of collision partners given by Q(H2O) = 1 × 1031s−1 and on the adopted collision cross section σcoll = 5 × 10−14 cm2. Model 3 is an improved version of model 2 that includes collisions with electrons (13). The effects of optical depth are negligible for HDO. The more elaborate models 2 and 3 are less sensitive to the gas temperature than model 1 because a sizable fraction of the total line flux originates from areas in the beam where fluorescence excitation dominates over collisional excitation. Model 3 with the preferred temperature gives an HDO production rate of Q(HDO) = (6.5 ± 1.1) × 1027 s−1 or a D/H ratio in water of (D/H)H2 O = (3.3 ± 0.8) × 10−4. The 1σ error for (D/H)H2 Ois the sum of the formal errors of Q(HDO) andQ(H2O).

Table 1

D/H of water in Hale-Bopp.T gas is the kinetic gas temperature (results for 80 K are shown to illustrate the temperature dependence); ∫T B dv is the calculated main beam efficiency-corrected line area for Q(HDO) = 1.0 × 1027 s−1; 〈N〉 is the beam-averaged column density for the measured value ∫T B dv = 1.85 K km s−1 (2); and Q(HDO) is the nuclear production rate of HDO. Note that (D/H)H2O= 0.5[Q(HDO)/Q(H2O)].

View this table:

We now have three values of D/H in cometary H2O (Fig.2). Of these, only the Halley value was an in situ measurement; hence, its precision is correspondingly greater (14). Nevertheless, it appears that the (D/H)H2 O values in these three comets are similar to each other and higher than the value for standard mean ocean water (SMOW), which is (D/H)H2 O = 1.56 × 10−4 (15). Although the three comets vary in their gas-to-dust ratio, size, and orbital period, the retrograde orbit of Halley and the long orbital periods and high inclinations of Hyakutake and Hale-Bopp (16) suggest that the Oort cloud is the source region for all three. It seems for the moment reasonable to accept the Halley value as characteristic of ices in Oort cloud comets, on the basis of this sample. Note that we do not yet have a value of D/H in a Jupiter family comet, that is, a short-period comet that is in a low-inclination direct orbit, indicating an origin in the Kuiper Belt (17). As all three of the comets appear to have come from the Oort cloud, statistics favor their origin in the Uranus-Neptune region of the solar nebula (18).

Figure 2

D/H ratio of Hale-Bopp (HB), Halley (14), and Hyakutake (Hyak) (12) as a function of the (current) inverse semimajor axis a. The wavy line shows the SMOW value (15), and the dashed-dotted line, the in situ value for Halley (14). Depicted are 1σ standard deviations. The error bar for Halley is about the size of the dot. The reciprocal semimajor axis 1/a is directly related to the orbital period Ta 2/3and is a crude measure of the dynamical age of the comet. For none of the comets on our list was it possible to determine an inverse semimajor axis that is back-extrapolated in time, (1/a)0, a quantity which would have provided a more realistic measure for the dynamical age (16).

The relatively high cometary (D/H)H2 Oratio helps us understand processes in the primordial nebula. Evidently, radial mixing was not very thorough because the mixture of ices that formed Earth's oceans is not the same as the mixture we find in the comets (14). We find high (D/H)H2 O values (1 × 10−4 to 5 × 10−4) in the cores of hot (∼200 K) interstellar clouds (19). In hot cores, the high values of HDO/H2O may reflect the composition of icy grain mantles that formed in an earlier cold phase (19). Millaret al. (20) have calculated that a steady-state value of 6.4 × 10−4 for HDO/H2O will be achieved in an interstellar cloud at a temperature of about 35 ± 10 K as a result of gas-phase ion-molecule reactions. Brown and Millar (21) have shown that still higher values of D/H could occur as a result of interactions with grain surfaces. The absence of ionization and the low ambient temperature (50 ± 20 K) in the Uranus-Neptune region of the solar nebula, the source region for Oort-cloud comets, ensure that little change in the interstellar medium values of D/H will occur during comet formation even if the icy grain mantles sublimate as the particles fall into the solar nebula mid-plane, to recondense on their refractory cores at the local temperature of the nebula (22). Our (D/H)H2 O measurement therefore strengthens the case for the primitive nature of cometary material. It invites further study of D/H in other H-bearing species in comets, as experience with interstellar clouds dictates that these values will not be the same as (D/H)H2 O in SMOW (23). Although we must still await determinations of HDO in classical short period comets to be certain, it appears that comets by themselves cannot be the only source of water for Earth's oceans (14). This result was anticipated from the high value of O/C in Earth's crust (24), but D/H provides a much stronger argument (14). Mixing of cometary water with water absorbed from the solar nebula by the rocky grains that accreted to form Earth may provide the observed value of (D/H)SMOW(25), as water vapor in the inner solar nebula could have D/H < 1 × 10−4 (26).

  • * Visiting astronomer at the James Clerk Maxwell Telescope, 660 North A'Ohoku Place, Hilo, HI 96720, USA.


View Abstract

Navigate This Article