arxiv: v1 [astro-ph] 22 Sep 2008

Transcription

1 Mon. Not. R. Astron. Soc. 000, 1 8 (2008) Printed 22 June 2009 (MN LATEX style file v2.2) Water in the Near IR spectrum of Comet 8P/Tuttle arxiv: v1 [astro-ph] 22 Sep 2008 R. J. Barber 1, S. Miller 1, N. Dello Russo 2, M. J. Mumma 3, J. Tennyson 1, P. Guio 1 1 Department of Physics & Astronomy, University College London, WC1E 6BT 2 Space Department, The Johns Hopkins University Applied Physics Laboratory, MD , USA 3 Solar System Exploration Division, NASA/GSFC, Greenbelt, MD 20771, USA Accepted XXX. Received XXX; in original form XXX 1 INTRODUCTION Cometary nuclei are the most primitive objects in the solar system. They accreted from the primordial solar nebula at a very early stage and those that are observed today were gravitationally-scattered into one of two distant, dynamical reservoirs, the Oort Cloud and the Kuiper Belt. It is believed that the structure and composition of comet nuclei were altered during and immediately after their early formation, principally as a result of heating of the interior by the radioactive decay of 26 Al, resulting in a stratified distribution of volatiles (Prialink et al., 2008). A second phase of evolution occurs when comet nuclei are perturbed into orbits that bring them close to the Sun. The outer layers are mainly affected by solar radiation which results in the loss of volatiles and structural changes in the region of the surface. Typically there is little activity at heliocentric distances greater than about 4.5 AU (Licandro et al., 2000). However, within about 3 AU of the sun, solar heating causes the surface temperature of the cometary nucleus to rise above 200 K; water ice sublimes and other volatiles and icy dust grains are expelled from the surface. The motion, chemistry and rate of sublimation of molecules from the icy grains are influenced by the solar wind and radiation. Initially the escaping gases form a dynamic, gravitationally-unbound atmosphere, the coma. They are subsequently swept into the comet s plasma and dust tails and eventually dissipate into interplanetary space. The surface of the nucleus is processed during the course of repeated visits to the inner solar system, becoming covered with a rubble blanket of particles too large to rjb@star.ucl.ac.uk ABSTRACT High resolution spectra of Comet 8P/Tuttle were obtained in the frequency range cm 1 on 3 January 2008 UT using CGS4 with echelle grating on UKIRT. In addition to recording strong solar pumped fluorescent (SPF) lines of H 2 O, the long integration time (152 minutes on target) enabled a number of weaker H 2 O features to be assigned, most of which had not previously identified in cometary spectra. These transitions, which are from higher energy upper states, are similar in character to the so-called SH lines recorded in the post Deep Impact spectrum of Comet Tempel 1, suggesting an alternative production mechanism. We derive an H 2 O rotational temperature of 62±5 K and a water production rate of (1.4±0.3) molecules s 1. Key words: 8P/Tuttle, water, SPF lines, SH lines, BT2. be dragged away by escaping gas. During the course of numerous perihelions, the outer surface becomes completely depleted of volatile materials which either escape from the surface or else migrate inwards. This creates an outer mantle of siliceous dust, mainly of an amorphous nature (Huebner, 2008), above a layered structure in which the more volatile species are found at greater depths (Prialnik et al., 2008). This structure results in the great variety that is observed in active comets. Typically cometary nuclei have radii 3 km. However, neither the nucleus, nor the inner coma (a region close to the nuclear where energy transfer is collisionally-dominated) are able to be resolved from Earth, and data obtained with Earth-based instruments are indicative of conditions some distance away from the inner coma. Consequently, our understanding of conditions within the inner coma is indirect, and relies largely on modelling. 2 COMET 8P/TUTTLE 8P/Tuttle is a short period comet, P orb = yr, and is the parent of the Ursid meteor stream (Jenniskens et al., 2002). Current estimates put its radius at 7.5 km, making it the largest of the group of 18 short period comets reported on by Licandro et al. (2000). In the scheme proposed by Levison (1996), comets are classified by reference to the Tisserand invariants. To a first order, this parameter is a constant of the motion of a comet. It is based on the Jacobi integral in the restricted threebody problem (Tisserand, 1894; Moulton, 1947). Because 8P/Tuttle has a Tisserand invariant of less than 2 (T J =

2 2 R.J. Barber et al ), it is classified as a near isotropic comet (NIC) of the Halley type. It is thought that NIC comets were formed in the same vicinity of the primordial nebula as the giant planets, that is within 4 40 AU of the centre (Dones et al., 2004), where the temperature was K. It is likely that this range of temperatures gave rise to variations in the volatile fractions of these comets, depending on the precise region of their formation (Magee-Sauer et al., 2008). Even though NIC comets were formed from remnant planetesimals in the region of the giant planets, gravitational interaction with the giant planets caused them to be ejected to very large distances, leading to the development of the Oort Cloud, a near-spherical region at a distance of around AU (Heisler, 1990) containing possibly comets. If the giant planets had been the only perturbers, the comets would have become unbound from the solar system. It is thought that their lifting into the Oort Cloud was due to tidal torquing by the galaxy and random passing stars (Dones et al., 2004). The formation of the Oort Cloud is thought to have taken place over a protracted period, possibly extending 1 billion years after planetary formation. Those NIC comets that are observed from Earth, have been perturbed in recent times, possibly by interactions with passing stars, into orbits that bring them back into the inner solar system. Because of the great distance of the Oort Cloud, most NIC comets have long orbital periods. Comet 8P/Tuttle is unusual in this respect. Their orbital periods decrease with successive perihelion passages. Eventually the NIC comet impacts with the Sun or else it is ejected from the planetary system. Most short period comets (defined as P orb 200 yr) have a Tisserand invariant 2 T J 3 and in the Levison scheme are referred to as ecliptic comets (EC). It is suggested that ECs are collisional fragments of Kuiper belt objects (Lamy et al., 2004) and that they formed within 35 AU of the centre of the primordial nebula and were subsequently pushed outward by Neptune s 1:2 mean motion resonance during its final phase of migration, into orbits lying between 40 and 50 AU (Levison and Morbidelli, 2003). If it is true that they originated in different parts of the protoplanetary disk, it is possible that EC and NIC comets have different physical properties, (Lamay et al., 2004). However, since many of the identifiable species have had their abundances measured in only a small number of comets, until recently it was not known whether there are any observable differences in the compositions of comets that correlate with their origins or evolution (A Hearn, 2004). Recently evidence has been presented for 3 organic classes based on 13 comets. Surprisingly, members of all three classes are found in each cometary reservoir (Mumma, 2008). Being a short period NIC, the chemistry of 8P/Tuttle is of particular interest. Unlike most short period comets, which are of the EC type, 8P/Tuttle is believed to have been formed in the region of the giant planets. It therefore has the potential to reveal possible differences between EC and NIC objects. Because it has made many approaches to the Sun, 8P/Tuttle s surfaces is highly processed like that of an EC. If then, despite having a highly-processed surface, its spectra reveal a chemistry that is more like that of other, less-well processed NICs, than of ECs, it may be possible to identify differences in the chemistries of various regions of the very early solar system. Even though, with the exception of 1953, 8P/Tuttle had been seen on each orbit since its discovery in 1858, until recently there had been no detailed investigation of its near- IR spectrum. During the 2008 apparition, 8P/Tuttle approached to within 0.25 AU of the Earth and was favourably placed for observing (see for example, Bonev et al., 2008 and Böhnhardt et al., 2008). 3 WATER IN COMETARY SPECTRA In the near IR, the spectra of comets are rich in rovibrational H 2O emission lines which can, in principle, be observed from Earth. However, because of the low temperatures that characterise cometary coma, the strongest water lines are fundamental transitions (that is to say, transitions to ground vibrational states) from low-energy ro-vibrational states. Photons from these transitions are absorbed by water vapour in the Earth s atmosphere, the molecules of which are in ground vibrational states. Bockelée-Morvan et al. (2004), give an overview of the main excitation processes in various regions of the cometary coma. Close to the nucleus the flow of gas in the coma is turbulent and complex, frequent collisions resulting in complex multi-body chemistry (Lisse et al., 2005). In this region, the density is sufficiently great (generally upward of 10 6 molecules cm 3 ) to ensure that the molecules are collisionally thermalised in low J ground vibrational states, with a rotational temperature equal to the local kinetic temperature. There is no precise boundary for the collisionallydominated region, as the breakdown of thermal equilibrium occurs at different distances from the nucleus for each of the allowed rotational transitions between ground vibrational states. Weaver and Mumma (1984) define a disequilibrium distance (D), which is the distance from the cometary nucleus at which the rotational excitation temperature differs from the kinetic temperature by more than 10%. D is dependent on the relative magnitudes of the collisional and the radiative rates, and is lower for transitions with short radiative lifetimes. Assuming an inner coma temperature of 50 K, and a molecular production rate, Q H2 O = s 1, Bockelée-Morvan (1987) calculates D to be km for the (000)[111] (000)[000] transition, but only 33 km for the (000)[330] (000)[303] transition (in this labelling system, vibrational quantum numbers are in round brackets, rotational quantum numbers are in square. Upper states are given first, see Barber et al., 2006). D is also dependent on Q H2 O. At 50 K, D for the (000)[221] (000)[212] transition is km when Q H2 O = s 1, but only 58 km when Q H2 O = s 1. Bockelée-Morvan (1987) demonstrates that at distances of up to several thousand kilometres from the nucleus, many rotational H 2O lines are optically thick and trap photons emitted by nearby molecules. This radiative trapping has the net effect of delaying the radiative decay to lower states, and hence maintaining local thermodynamic equilibrium (LTE) to greater values of D than would be the case in optically thin conditions. This intermediate region has not yet been adequately modelled, and consequently the various excitation and de-excitation processes are not well-quantified. However, Xie and Mumma (1992) showed that at distances of several thousand kilometres, e H 2O collisions play an

3 Water in the Near IR spectrum of Comet 8P/Tuttle 3 important role in exciting water molecules. In Section 9 we discuss the possibility that electron collisions may be involved in the formation of the so-called SH water lines that were observed to be present in our 3 January 2008 UT spectrum of 8P/Tuttle. SH lines were first identified in the Deep Impact spectrum of Comet Tempel 1 (Barber et al., 2007), and were subsequently noted to be present in another spectrum of the same comet (Mumma et al., 2005, Figure 2) Their characteristics are discussed in Section 6. At distances greater than a few thousand kilometres from the nucleus the density of the species is sufficiently low, and the mean free path between collisions large enough, for fluorescence equilibrium to exist, that is to say, radiation is due to the balance between solar pumping and subsequent spontaneous radiative decay (Crovisier, 1984). Despite being produced in low temperature regions, cometary spectra contain H 2O lines that, under conditions of local thermodynamic equilibrium (LTE), would only be observed at temperatures of several thousand Kelvin. This situation arises because in cometary comae most photons are emitted in the low density, radiation-dominated region, and not in the collisionally-dominated inner region. Transitions from these excited states to lower energy ro-vibrational states (hot bands) are not absorbed in the Earth s atmosphere. Spectra containing these solar-pumped fluorescent (SPF) transitions enable information to be derived about the rotational temperature, T rot of the inner coma and the water production rate, Q H2 O. A comparison of this latter value with the production rates of trace species provides the abundance ratios that are fundamental to understanding cometary chemistry. Because Boltzmann statistics do not apply in non-lte regions, the derivation of physical parameters from the intensities of the SPF water lines presents difficulties. Each upper excited ro-vibrational state can be populated directly by radiative pumping of rotational levels of the ground vibrational state, as well as by cascade from higher energy states (Crovisier, 1984). The relative probabilities of these various upward and downward routes are temperature-dependent. Models have been developed to compute the g-factors for each SPF line, from which line intensities are able to be calculated as a function of temperature (Crovisier, 1984; Mumma et al., 1995, Dello Russo et al., 2000). These models require a knowledge of the Einstein A coefficients for the downward transitions. Ideally, the upward pumping rates should be computed using the Einstein B coefficients. However, for simplicity, the temperature-dependent vibrational band intensities are frequently used to compute the upward pumping rates. For example, Dello Russo et al. (2004) used Einstein A coefficients from the BT2 synthetic water line list (Barber et al., 2006) to model g-factors for all the H 2O SPF transitions in the: ν 1+ν 3 ν 1, ν 1+ν 3 ν 3, 2ν 1 ν 1, 2ν 1 ν 3, ν 1+ν 2+ν 3 ν 1 ν 2 vibrational bands, from upper states having J 7. The µm ( cm 1 ) spectral region is particularly rich in SPF lines and is also largely devoid of other molecular species. The availability of accurate g- factors has enabled the region to be used in determining cometary rotational temperatures and water production rates (e.g. Dello Russo et al., 2004 and 2006). 4 RECENT DATA Using photometric measurements, several observers estimated rotation times of the nucleus of 8P/Tuttle as being either 5.71±0.04 hr (Schleicher and Woodney, 2007b), 7.7±0.2 hr (Harmon et al., 2008), or multiplicities of 5.7 and hr (Drahus et al., 2008). In addition, radar images obtained on 2-4 January 2008 (Harmon et al., 2008) show a bifurcated nucleus consisting of two lobes each 3 4 km in diameter, suggesting a possible contact binary. Our UKIRT observations took place over a period of more than 2.5 hrs, which based on these estimates is a significant fraction of the rotational period of the bifurcated nucleus. During the course of our observations we noted a change in the relative intensities of some of the strong water features, which may have been due to the two nuclei having different out-gassing characteristics (Bonev et al., 2008). In particular, the SPF line at cm 1 was observed to be stronger than the SPF line at cm 1 in the first half of our observing session, whilst the converse was true in the second half. The intensities of the two lines were similar for the night as a whole, as seen in Figure 1.The observed changes could possibly be explained by a single rotating nucleus with different parts of the cometary surface, and different out-gassing events, coming within our line of sight. However, it is significant that over the temperature range K, the ratio of the g-factors of the lines at cm 1 and cm 1 is constant, also the ratio of the g-factor of the the line at cm 1, which is partly blended with that at cm 1, to those of the other two lines varies by less than 3% over this 30 K range. An alternative explanation is a variation in the intensity over the period of the observations of an SH line at cm 1, (003)[431]-(002)[532]. This line, which has A if = 73.2 s 1 and E upper = cm 1, is a strong candidate for being present in the spectra. UKIRT is unable to resolve this line from the strong SPF line at cm 1. There are several estimates of the H 2O production rates. Two are based on OH production rates obtained using narrow-band photometry. These are: molecules s 1 on 1 Nov 2007 UT, at a heliocentric distance of 1.63 AU (Schleicher, 2007) and molecules s 1 heliocentric distance of 1.30 AU, averaged from observations on 3, 4, 5 Dec 2007 UT, (Schleicher and Woodney, 2008a). Neither of these figures comes with an estimate of error. Bonev et al., 2008 give two values obtained on 2007 December UT, using different Keck NIRSPEC settings, when the comet was at a heliocentric distance of 1.15 AU. These are: (2.38 ± 0.06) and (2.13 ± 0.11) molecules s 1. Böhnhardt et al., 2008, give water production rate of (5.97 ± 0.27) on 27 January 2008 UT, when the heliocentric distance was 1.03 AU. 5 OBSERVATIONS We observed 8P/Tuttle on its approach to the Sun (perihelion being on 27 January UT), on the nights of 3, 4 and 5 January 2008 UT from UKIRT, Hawaii, using the echelle grating on CGS4. We report here on the spectra that we obtained on 3 January UT, in the frequency range cm 1. The heliocentric distance on this date was

4 4 R.J. Barber et al AU, and the geocentric distance and velocity were 0.25 AU and 3.3 km s 1 respectively; the latter figure corresponds to a red shift of cm 1, which was equivalent to 1 pixel on the array, and less than the minimum resolution of our instrument, which at R = was cm 1 at 3450 cm 1. The total visual magnitude of 8P/Tuttle at the time of our observations was 5.8. However our observations were of the bright inner coma region which presented a point source in our instrument, and required long integration times. In the echelle mode, each pixel on the array is equivalent to 0.86 arcsec on the sky in the spatial direction and 0.41 arcsec perpendicular to the slit. The slit, of length 90 arcsec, was oriented east-west. At the distance of 8P/Tuttle, its 2-pixel width gave a field of view of 149 km perpendicular to the slit. The telescope was dithered to give 2x2 sampling every 40 seconds. Observations were acquired using a standard ABBA sequence, which was achieved by nodding the telescope 22 pixels along the slit. In this mode, the comet signal was present in both the A and B beams (rows 114 and 92 respectively), which compared to nodding to blank sky, increased the signal to noise ratio by a factor of 2, when both signals were added using the (A-B)-(B-A) routine to remove the sky. Each ABBA series represented a total of 160 seconds integration time. Total time on-target was 152 minutes. Spectra of an A0V star, HD 6457, were obtained at the beginning and end of the observing session. Standard star spectra are required in order to adjust for frequency-dependent atmospheric transmission, and also to enable absolute flux calibration. Before dividing by the standard star, we removed high frequency noise from both the 8P/Tuttle and standard star data using Fourier transform smoothing. We rejected all data at frequencies where the atmospheric transmission was 20%. After dividing by the standard star, flux calibration was achieved by multiplying the result by the frequency-dependent flux of HD 6457 over our 22 cm 1 frequency range, calculated relative to Vega as the photometric 0.0 magnitude standard (Hayes, 1985), and assuming a black body function with T eff = 9480 K (Tokunaga, 2000, p151). Vacca et al. (2003), describe a method for correcting near-ir spectra for telluric absorption using an A0V star, which is particularly useful in the I, z, Y and H bands where the numerous higher level lines of the Paschen and Brackett series are located. However, at the short wavelength end of the L band, where we were observing, hydrogen absorption in the atmospheres of A0V stars is not an issue (Tokunaga, 2000, p156), and we therefore used a simple black body function. Because of the lack of arc lines in the cm 1 region, we frequency-calibrated the spectrum using the position of several known SPF lines that were present in the 8P/Tuttle spectrum. This had the effect of transforming the observed spectrum into the rest frame. The fullyreduced, flux-calibrated spectrum is given in Figure 1. 6 DATA ANALYSIS We compared the features in Figure 1 with our BT2 database and our results are detailed in Table 2. Many were identified as SPF lines (see Section 3). However Figure 1 includes a number of other features, whose shape and intensity above the continuum suggest that they are spectral lines, rather than noise, but whose frequencies do not correspond to SPF lines. In order to assign these other features, we employed the methodology detailed in Barber et al. (2007) for the identification of non-spf features in the post Deep Impact spectrum of Comet Tempel 1. We used BT2 to generate synthetic emission spectra for water at a temperature of 3000 K in the cm 1 frequency range, initially restricting the value of J for the upper state to 8, but subsequently to 5 (as none of the observed transitions involved states with higher values of J), and the intensity of the weakest line to be generated was limited to 10 4 of the strongest. The output, which consisted of 1400 individual lines was convolved to match the resolving power of the CGS4 echelle (37 000). Some of the features in this synthetic spectrum were at frequencies corresponding to features in the spectrum of 8P/Tuttle. The relative intensities of these features in the synthetic and observed spectra did not agree, which was to be expected as the observed spectrum was produced in a non-lte region, whereas the synthetic spectra assume Boltzmann statistics. More problematic was the fact that there were many more features in our synthetic spectrum than there were in the spectrum of 8P/Tuttle. Our next step was therefore to examine the individual transitions in the synthetic spectrum. When synthetic spectra are generated using BT2, one of the output files gives the frequency, assignment, intensity, Einstein A coefficient and lower state energy of every water transition in the frequency range, down to a specified intensity cut-off. Using these data, we were able to identify certain characteristics of those lines in our synthetic spectra whose positions matched those of features in the observed spectrum, that differentiated them from those that did not. Before describing our findings, it will be helpful to comment briefly on the ro-vibrational states of water. The H 2O molecule has three vibration modes: symmetric and asymmetric stretch and bend. Different amounts of energy are required to excite each of these modes by one quanta, less incremental energy being required as the total internal energy of the excited molecule rises. In the energy range in which we are interested, each each quanta of symmetric and asymmetric stretch represents cm 1. These values are in the region of twice the energy of one quanta of bend. Hence it is convenient to talk in terms of polyads, where each polyad, designated ν is the equivalent of one quanta of stretch or two of bend. Single bend quanta are denoted by δ. We noted that all the SPF lines in the observed spectrum came from the 2ν polyad, whilst of the non-spf lines that matched features in our synthetic spectra, five were from the 3ν polyad and two from 3ν+δ. Also, none of the lines that we were able to assign had Einstein A coefficients with a value of less than 1.0 s 1. The Einstein A coefficients of the SH lines were generally higher than those of the SPF lines, the average A if being 20.6 and 7.3 s 1 respectively. The SPF lines had upper states whose energies were all in a very narrow range: cm 1, whilst the non-spf lines had upper states with energies in the range cm 1. Rotation excitation, J, also affects which lines are ob-

5 Water in the Near IR spectrum of Comet 8P/Tuttle 5 served in cometary spectra, and an examination of the transitions in our synthetic spectrum revealed that none of the features in the spectrum of 8P/Tuttle corresponded to lines whose upper state had J 5. Finally, we noted that none of the emission lines in our K synthetic spectrum that corresponded to features in the observed spectrum had intensities less than 4.9x10 19 Wm 2. In summary, the non-spf lines that we identified in the spectrum of 8P/Tuttle had the same characteristics as those observed in the Deep Impact spectrum of Comet Tempel 1 (Barber et al., 2007), and which these authors refer to as stochastic heating or SH lines. Our assignments of 20 lines in the spectrum of 8P/Tuttle, together with E upper, A if and the type of transition, SPF/SH and O/P (ortho/para) are given in Table 2. The table includes all the strongest SPF features (including one that is a blend of two SPF lines) that we would expect to see in this spectral region based on our derived inner coma temperature of 62 K. There are also 6 unblended SH lines. In addition there is one feature centred at cm 1 which was identified as a blend of an SPF line at cm 1 and an SH line at cm 1. In fact, at 62 K, the g-factor for the SPF component of the blend is , which equates to a line of only 3% of the intensity of the feature at cm 1, and indicates that the SH line is by far the larger of the two components in the blend. Barber et al. (2007) assigned four SH transitions in the spectrum of Comet Tempel 1. Of these, the two at and cm 1 are blended. In the light of our current work, we are inclined to the view that the feature at this frequency is more likely to be the (211)[322] (210)[221] line at cm 1, which has a 3ν+δ polyad upper state of energy cm 1, than to be blended with the (103)[110] (102)[110] line at cm 1, which unlike any of the SH transitions in Table 2, has a 4ν polyad upper state and an upper level of cm 1. In our spectrum of 8P/Tuttle, we were not able to identify a line at cm 1, the position of one of the SH features reported by Barber et al.. This may be because it was blended with the wing of the broad SPF line at cm 1. The fourth SH line assigned by Barber et al. is at cm 1. This corresponds to a weak feature in Figure 1. We examined the spectra of comet C\1999 H1 Lee (Dello Russo et al., 2006), which contained several unidentified features. We noted that the frequencies of three of these features corresponded to SH lines in the spectrum of 8P/Tuttle. These are: (022)[431] (120)[542] at cm 1 (211)[322] (210)[211] at cm 1 (220)[212] (021)[111] at cm 1. The first of these lines is outside the frequency range examined in Barber et al. (2007), but as stated above, the other two lines were also observed in Comet Tempel 1. In addition, we identified several water lines that were unassigned in the high resolution spectrum of 8P/Tuttle obtained at ESO with CRIRES on 27 January 2008 UT, when the heliocentric distance was 1.03 AU (Böhnhardt et al., 2008, Chart C). These include SPF lines at , and cm 1, all of which are also present in Table 2. We identified a number of SH lines, including (013)[330] (012)[413] at cm 1 (also in Table 2) and (201)[551] (101)[624] at cm 1, (032)[514] Table 1. List of comets in which SH water lines have been observed Name P (yr) Type Origin r h T rot K 9P/Tempel JFC Kuiper C/1999 H1(Lee) LPC Oort P/SW JFC Kuiper P/Tuttle NIC Oort JFC: Jupiter Family Comet LPC: Long Period Comet NIC: Near Isotropic Comet 1 artificial impact 2 fragmented 3 bifurcated nucleus Flux (10-15 W m -2 µm -1 ) SH ( ) SPFs ( blend) SPF( ) SH ( ) SPF ( ) SH ( ) SPF ( ) SPF ( ) SH ( ) SPF ( ) SPF ( ) +SH( )? Frequency (cm -1 ) SH ( ) SPF ( ) 3452 SPF ( ) SPF ( ) Figure 1. Spectrum of Comet 8P/Tuttle, UT 3 January Data, which have been Fourier Transform smoothed to remove high frequency noise, have been divided by standard star HD Only frequency regions where atmospheric transmission is greater than 20% are shown (031)[533] at cm 1 and (230)[313] (031)[212] at cm 1, which do not appear in Table 2. Finally, we note that the unassigned feature at cm 1 in the spectrum of Comet 73P/Schwassmann- Wachmann 3C obtained by N.Dello Russo on 15 May 2006 UT (private communication), is the (211)[322] (210)[211] line (also in Table 2). On the basis of our identification of SH lines in comets of different types, we conclude that this class of water lines, which until recently had not been identified in cometary spectra, is probably always present at heliocentric distances of up to 1.5 AU. Moreover, it is clear that the actual lines observed vary from comet to comet, and indeed between spectra of the same comet obtained on different dates. Table 1 includes details of four comets in which we have identified SH transitions in the frequency range cm

6 6 R.J. Barber et al. 7 RO-VIBRATIONAL EXCITATION OF H 2O MOLECULES The temperatures of the water vapour in the inner coma of comets that are within 2 AU of the sun are generally in the range K (Combi et al., 2004). Close to the nucleus, molecules are collisionally-thermalised in low J pure rotational states. Using the BT2 line list and assuming LTE, it is easy to demonstrate that at the temperature of 8P/Tuttle ( 62 K), all states with J 5 have < 0.1 % probability of being occupied. Absorption of solar photons by gaseous H 2O in these ground vibrational states excites the molecules into more energetic ro-vibrational states. The absorption is governed by selection rules, J = 0, ±1 (J 0-0 forbidden). Because of the low probability that pure rotational states with J 5 will be populated in the inner coma at a temperature of 62 K we would not expect to observe SPF lines whose upper states have J 6. In fact Table 2 does not include any such lines with J upper > 4. It will be seen, however that a high proportion of the lines do in fact have J upper = 4. A rigorous derivation of the relative probability of lines with differing J upper being observed is beyond the scope of this paper. It depends inter alia, upon: the temperature of the inner coma, the water production rate, solar flux, Einstein A coefficients of the various ground vibrational from which the upper levels are pumped, and upper ro-vibrational states of the downward transitions, together with appropriate branching ratios, and also the limiting S/N ratio of the observations. 8 ROTATIONAL TEMPERATURE AND H 2O PRODUCTION RATE The procedure that we used to determine the temperature of the inner coma is based on the fact that the g-factors of some SPF transitions can behave quite differently as a function of temperature: some g-factors increase with temperature, whilst others decrease (Dello Russo et al., 2004). In broad terms, SPF lines whose upper states have been pumped from higher J values become more intense as temperature rises, whilst those pumped from the lower J states show little increase in intensity at higher temperatures, and in the case of those pumped from the lowest rotational states, may actually weaken with increasing temperature. Using the g-factors for the SPF lines in our spectral frequency range (Dello Russo et al., 2004), we generated SPF synthetic spectra, convolved to the resolving power of our instrument. The spectra were produced using the Fortran program, spectra-bt2 (Barber et al., 2006) which is available in electronic form via: and were generated at various temperatures with 5 K increments (the smallest temperature difference that produced measurable differences in the synthetic spectra). The spectral region contains 38 SPF transitions in five vibrational bands, for which g-factors computed by Dello Russo are available. Many of these lines are not detectable at low rotational temperatures, whilst others are at frequencies that fall within regions of low atmospheric transmission. Of the 38, only 12 could be used for temperature diagnostics (two of which appear blended), and these lines corresponded Table 2. Assignments of SPF and previously unobserved SH lines in the post-impact spectrum of 8P/Tuttle. In order, the columns give: the rest frequency in cm 1 ; the assignment, upper level given first (vibrational quantum numbers in round brackets, rotational quantum numbers in square brackets); energy of the upper level in cm 1 ; Einstein A coefficient from BT2; type of line, SPF/SH, O/P (ortho/para) Frequency Identification E upper A if Type (cm 1 ) (see text for notation) (cm 1 ) s (022)[431] - (120)[542] a SH P (101)[111] - (001)[202] SPF O (101)[431] - (100)[532] SPF O (200)[423] - (100)[432] SPF O (022)[432] - (120)[541] SH O (101)[000] - (001)[111] SPF O (003)[432] - (002)[533] SH P (101)[422] - (100)[523] SPF O (200)[303] - (100)[414] SPF O (013)[330] - (012)[431] b SH P (101)[211] - (001)[220] SPF O (211)[322] - (210)[221] d SH O (200)[110] - (100)[221] SPF O (101)[202] - (100)[321] SPF O (220)[212] - (021)[111] ac SH O (101)[413] - (001)[422] SPF O (200)[110] - (001)[111] SPF O (220)[101] - (120)[212] SH O (200)[220] - (001)[221] SPF P (200)[414] - (100)[505] SPF O Blended 1 See Section 6 2 Possibly blended with: (201)[202] - (200)[101] SH O 3 Possibly blended with: (003)[431] - (002)[532] SH O a Also in Lee, Dello Russo et al., 2006 b Also in 8P/Tuttle Böhnhardt et al., 2008 c Also in Tempel 1, Barber et al., 2007 d In abc and SW3 (see Section 3) precisely to the SPF features in our observed spectrum and which are detailed in Table 2. The SPF/SH blend at cm 1 (see Section 6) was not used for temperature diagnostics. On comparing the relative intensities of the synthetic spectra with our observed spectrum, normalised to the cometary continuum level, we observed that, there was excellent agreement between the relative intensities of all the SPF lines at 62 ±5 K, and consequently we conclude that this was the rotational temperature of the inner coma of 8P/Tuttle at UT 3.30 January 2008 (the mid-point of our observations). It should be noted that the g-factors take account of the temperature-dependent ratio of the ortho and para forms of the H 2O molecule. Above 50 K, this is the same as the ratio of the nuclear spin degeneracies of the ortho and para states, viz 3:1. However, because the ground state (no rotation, no vibration) is a para state, and the lowest ortho state for water is at cm 1, equivalent to a temperature of 34 K, lower ratios would have needed to be used if the temperatures had been below 50 K.

7 Water in the Near IR spectrum of Comet 8P/Tuttle 7 In computing the rate of production of gaseous H 2O molecules we used the methodology detailed in Dello Russo et al. (2004). This includes an adjustment for slit losses, which can have the effect of lowering product rates derived from nucleus-centred spectra. The result from our 3 Jan 2008 UT observations of 8P/Tuttle was 1.4 ± molecules s 1. OH fragments rather than excitation transfer (Dunlea and Ravishankara, 2004). Tests of these possible processes will require comparisons of SH in various comets with model predictions, such as dependence on cometary activity and the variations with nucleocentric distance (i.e., local density). However, such comparisons are beyond the scope of this paper. 9 ORIGIN OF SH LINES We comment here on the possible origins of the SH lines that we have identified in cometary spectra and which are characterised by transitions from 3ν or higher polyad states and low J values (typically 5 and below). We consider three mechanisms by which the higher energy states may be populated: direct excitation by electrons, dissociative recombination of H 3O +, and the photodissociation of water and the subsequent quenching of O( 1 D) by an H 2O molecule. Xie and Mumma (1992) demonstrate that in the case of active comets, electron-water collisions can play an important role in populating the rotational states. They base their analysis on Giotto measurements of the ion temperatures in the coma of Comet Halley (Lämmerzahl et al., 1987) and assume the electron temperature is approximately equal to the ion temperature. The temperature of the electrons rises, and hence the cross-section for e H 2O collisions declines (Faure et al., 2004) with distance from the nucleus. In addition, away from the inner coma, electron density is assumed to fall-off as the inverse square of the nuclear distance. Hence the greatest effect of electron excitation will be within the neutral-neutral collisionally-dominated inner coma, where T e 200 K. However, Xie and Mumma (1992) show that the rate of e H 2O collisions is still significant at distances of km from the nucleus, where T e is > K, which represents energies that are able to excite the 3ν polyad in the H 2O molecule. It is important to note that dipole-forbidden electron excitation processes (e.g., pure vibrational excitation) have maximum cross sections near the threshold energy, in this case about cm 1. For pure rotational excitation, monopole-dipole interactions are important. The second process also involves electrons. The recombination of H 3O + (formed as a result of dissociation of water) gives rise to H 2O and H with 25% efficiency, and the fragments carry 6.4 ev of excess energy (Jensen et al., 2000). This exceeds the bond energy of H 2O, but it is very likely that a large fraction of the energy is carried by fast H-atoms, leaving the H 2O in bound states that are highly vibrationally excited. In the last of the three processes, photolysis of H 2O produces O( 1 D) atoms with a quantum yield of 5% under quiet Sun conditions (Huebner et al., 1991). O( 1 D) is long lived ( 110 s), so it can collide with water in the inner coma before radiative relaxation. About 1.97 ev ( cm 1 ) is available upon collision, more than enough to excite the 3ν polyad, and even the 4ν polyad. Subsequent vibrational cascade could feed the observed SH lines, but detailed models are needed to evaluate this. An important issue is the balance between excitation transfer and dissociative reactions. Some laboratory work suggests efficient production of two 10 ACKNOWLEDGMENTS We are grateful to have been granted observing time on UKIRT, and acknowledge the assistance that we have received from the directors and staff of the observatory. 11 REFERENCES A Hearn M.F., 2004, In Comets II, Festu M.C., Keller H.U., Weaver H.A. eds., UA Space Science Series, Arizona, Barber R.J., Tennyson J., Harris G. J., Tolchenov R., 2006, MNRAS, 368, 1087 Barber R.J., Miller S., Stallard T., Tennyson J. Hirst P., Carroll T., Adamson A., 2007, Icarus, 187, 167 Bockelée-Morvan D., 1987, A&A, 181,167 Bockelée-Morvan D., Crovisier J., Mumma M. J., Weaver H.A., 2004, In Comets II, Festu M.C., Keller H.U., Weaver H.A. eds., UA Space Science Series, Arizona, Böhnhardt H., Mumma M. J., Villanueva G. L., DiSanti M.A., et al., 2008, ApJ., 683, L71 Bonev B. P., Mumma M.J., Radeva Y.L., DiSanti M. A., Gibb E. L., Villanueva L., 2008, ApJ., 2008, 680, L61 Combi M.R., Harris W.M., Smyth W.H., 2004, In Comets II, Festu M. C., Keller H.U., Weaver H.A. eds., UA Space Science Series, Arizona, 18 Crovisier J., 1984, A&A, 130, 361 Crovisier J., Encrenaz T., 1983, A&A, 126, 170 Dello Russo N., Mumma M. J., DiSanti M. A., Magee-Sauer K., Navak R., Rettig T.W., 2000, Icarus, 143, 324 Dello Russo N., DiSanti M. A., Magee-Sauer K., Gibb E. L., Mumma M. J., Barber R. J., Tennyson J., 2004, Icarus, 168, 186 Dello Russo N., Mumma M. J., DiSanti M. A., Magee-Sauer K., Gibb E. L., 2006, Icarus, 184, 255 Dones L., Weissman P.R., Levison H.F., Duncan M.J., 2004, ASP Conference Series, Ed., Johnstone D., et al., 323, 371

8 8 R.J. Barber et al. Drahus M., Jarchow C., Hartogh P., Waniak W., Bonev T., et al., 2008, CBAT, 1294 Dunlea E. J. and Ravishankara A. R., 2004, Phys. Chem. Chem. Phys., 6, 3333 Faure A., Gorfinkiel J. D., Tennyson J., 2004, MNRAS, 347, 323 Harmon J.K., Nolan M.C., Howell E.S., 2008, CBAT, 8909 Hayes D.S., 1985, IAUS, 111, 225 Heisler J., 1990, Icarus, 88, 104 Huebner W. F., 2008, Space Sci. Rev., 138, 5 Huebner W. F., Keady J. J., Lyon S. P., 1992, Astrophys. Sp. Sci., 195, 1 Jenniskens P., Lyytinen E., de Lignie M. C., Johannink C., et al., 2002, Icarus, 159, 197 Space Sci. Rev., 138, 147 Schleicher D., 2007, CBET, 1113 Schleicher D. and Woodney L., 2007a, IAUC, 8903 Schleicher D. and Woodney L., 2007b, IAUC, 8906 Tokunaga A. T., 2000, In Allen s Astrophysical Quantities, ed. Cox A.N., fourth ed.; Springer, New York Tisserand F., 1894, In Traité de Mécanique Céleste, Vol. III, Gauthier-Villars, Paris Vacca W.D., Cushing M. C., Rayner J.T., 2003, Astron Soc. Pacific, 115, 389 Weaver H.A. and Mumma M. J., 1984, ApJ., 276, 782 Xie X. and Mumma M. J., 1992, ApJ., 386, 720 Jensen M.J., Bilodeau R.C., Safvan C. P., Seiersen K., Andersen L.H., 2000, ApJ., 543, 764 Lämmerzahl P., Krankowsky D., Hodges R. R., Stubbernann U., et al., 1987, A&A, 187, 169 Lamy P.L., Toth I., Fernández Y.R., Weaver H.A., 2004, In Comets II, Festu M.C., Keller H.U., Weaver H.A. eds., UA Space Science Series, Arizona, 224 Levison H. F., 1996, In Completing the inventory of the Solar System, Rettig T.W. and Hahn J. M. eds., 173, ASP Conference Series, Astron. Soc. Pacific Levison H.A. and Morbidelli A., 2003, Nature, 426, 419 Lisse C. M., A Hearn M. F., Farnham T. L., Groussin O., Meech K.J., Fink U., Schleicher D.G., 2005, Space Sci. Rev., 117, 161 Licandro J., Tancredi G., Lindgren M., Rickman H., Hutton G.R., 2000, Icarus, 147, 161 Magee-Sauer K., Mumma M. J., DiSanti M. A., Dello Russo N., Gibb E. L., Bonev B.P., Villanueva G.L., 2008, Icarus, 194, 347 Moulton, F. R., 1947, In An Introduction to Celestial Mechanics 2nd rev. ed., Sect. 159, MacMillan, London Mumma M. J., DiSanti M.A., Tokunaga A., Roettger E. E., 1995, Bull. Am.Astron.Soc., 27, 1144 Mumma M.J., DiSanti M. A., Magee-Sauer K., Bonev B. P., Villanueva G.L., et al., 2005, Science 310, 270 Mumma M. J., 2008, Proceedings of Tenth ACM Conference, Baltimore (in prep.) Prialnik D., Sarid G., Rosenberg E. D., Merk., 2008,