# Molecular composition of short-period comets from millimetre-wave spectroscopy: 21P/Giacobini-Zinner, 38P/Stephan-Oterma, 41P/Tuttle-Giacobini-Kresák, and 64P/Swift-Gehrels

N. Biver, D. Bockelée-Morvan, D. C. Lis, D. Despois, R. Moreno, J. Crovisier, P. Colom, J. Boissier, N. Dello Russo, R. J. Vervack, et al.

► To cite this version:

N. Biver, D. Bockelée-Morvan, D. C. Lis, D. Despois, R. Moreno, et al., Molecular composition of short-period comets from millimetre-wave spectroscopy: 21P/Giacobini-Zinner, 38P/Stephan-Oterma, 41P/Tuttle-Giacobini-Kresák, and 64P/Swift-Gehrels. Astronomy and Astrophysics - A&A, 2021, 651, pp.A25. 10.1051/0004-6361/202140765. [hal-03285131](https://hal.science/hal-03285131)

HAL Id: hal-03285131

https://hal.science/hal-03285131

Submitted on 13 Jul 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
# Molecular composition of short-period comets from millimetre-wave spectroscopy: 21P/Giacobini-Zinner, 38P/Stephan-Oterma, 41P/Tuttle-Giacobini-Kresák, and 64P/Swift-Gehrels*,**

N. Biver<fnref n="1" />, D. Bockelée-Morvan<fnref n="1" />, D. C. Lis<fnref n="2" />, D. Despois<fnref n="3" />, R. Moreno<fnref n="1" />, J. Crovisier<fnref n="1" />, P. Colom<fnref n="1" />, J. Boissier<fnref n="4" />
N. Dello Russo<fnref n="5" />, R. J. Vervack Jr.<fnref n="5" />, M. Cordiner<fnref n="6" />,<fnref n="7" />, S. Milam<fnref n="6" />, N. X. Roth<fnref n="6" />, B. P. Bonev<fnref n="8" />, M. A. DiSanti<fnref n="9" />
J. K. Davies<fnref n="10" />, and H. Kawakita<fnref n="11" />

Received 9 March 2021 / Accepted 17 May 2021

## ABSTRACT

We present the results of millimetre-wave spectroscopic observations and spectral surveys of the following short-period comets: 21P/Giacobini-Zinner in September 2018, 41P/Tuttle-Giacobini-Kresák in April 2017, and 64P/Swift-Gehrels and 38P/Stephan-Oterma in December 2018, carried out with the Institut de RadioAstronomie Millimétrique (IRAM) 30-m radio telescope at wavelengths between 1 and 3 mm. Comet 21P was also observed in November 1998 with the IRAM 30-m, *James Clerk Maxwell Telescope*, and the Caltech submillimeter Observatory radio telescopes at wavelengths from 0.8 to 3 mm. The abundances of the following molecules have been determined in those comets: HCN, CH₃OH, CS, H₂CO, CH₃CN, and H₂S in comet 21P; HCN and CH₃OH in 41P; HCN, CH₃OH, and CS in 64P; and CH₃OH in 38P. The last three comets, classified as carbon-chain typical from visible spectro-photometry, are relatively rich in methanol (3.5–5% relative to water). On the other hand, comet 21P, classified as carbon-chain depleted, shows abundances relative to water which are low for methanol (1.7%), very low for H₂S (0.1%), and also relatively low for H₂CO (0.16%) and CO (<2.5%). Observations of comet 21P do not show any change in activity and composition between the 1998 and 2018 per-ihelions. Sensitive upper limits on the abundances of other molecules such as CO, HNCO, HNC, or SO are also reported for these comets.

**Key words.** comets: general – comets: individual: 21P/Giacobini-Zinner – radio lines: planetary systems – submillimeter: planetary systems

## 1. Introduction

Comets are the most pristine remnants of the formation of the Solar System 4.6 billion years ago. They sample some of the oldest and most primitive material in the Solar System, including ices, and they are thus our best window into the volatile composition of the proto-solar disk. Comets may have also played a role in the delivery of water and organic material to the early Earth (see Hartogh et al. 2011, and references therein). The latest

simulations of early Solar System evolution (Brasser & Morbidelli 2013; O'Brien et al. 2014) suggest a more complex scenario. On the one hand, ice-rich bodies that formed beyond Jupiter may have been implanted in the outer asteroid belt and participated in the supply of water to Earth or, on the other hand, current comets coming from either the Oort Cloud or the scattered disk of the Kuiper belt may have formed in the same trans-Neptunian region, sampling the same diversity of formation conditions. Understanding the diversity in the chemical and isotopic composition of cometary material is thus essential in order to assess such scenarios (Altwegg & Bockelée-Morvan 2003; Bockelée-Morvan et al. 2015).

Comet 21P/Giacobini-Zinner is a Jupiter-family comet (JFC) orbiting the Sun every 6.5 yr on a moderately low-inclination (32.0°) orbit. It was first discovered by M. Giacobini in 1900 and

<footnote n="1">¹ LESIA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Université, Université de Paris, 5 place Jules Janssen, 92195 Meudon, France</footnote>

<footnote>e-mail: nicolas.biver@obspm.fr</footnote>

<footnote n="2">² Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA</footnote>

<footnote n="3">³ Laboratoire d'astrophysique de Bordeaux, Université de Bordeaux, CNRS, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France</footnote>

<footnote n="4">⁴ IRAM, 300 rue de la Piscine, 38406 Saint Martin d'Hères, France</footnote>

<footnote n="5">⁵ Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Rd., Laurel, MD 20723, USA</footnote>

<footnote n="6">⁶ Solar System Exploration Division, Astrochemistry Laboratory Code 691, NASA-GSFC, Greenbelt, MD 20771, USA</footnote>

<footnote n="7">⁷ Department of Physics, Catholic University of America, Washington, DC 20064, USA</footnote>

<footnote n="8">⁸ Department of Physics, American University, Washington, DC, USA</footnote>

<footnote n="9">⁹ Solar System Exploration Division, Planetary Systems Laboratory Code 693, NASA-GSFC, Greenbelt, MD 20771, USA</footnote>

<footnote n="10">¹⁰ UK Astronomy Technology Centre, Royal Observatory, Edinburgh, UK</footnote>

<footnote n="11">¹¹ Koyoma Astronomical Observatory, Kyoto Sangyo University Motoyama, Kamingamo, Kita-ku, Kyoto 603-8555, Japan</footnote>

<footnote>* The radio spectra are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/651/A25</footnote>

<footnote>** Based on observations carried out with the IRAM 30-m telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain).</footnote>
rediscovered by E. Zinner in 1913. It has since been observed at 14 out of 16 perihelion passages. Its perihelion distance oscillates in the 0.93–1.10 au range during the 20th and 21st century with best apparitions when perihelion happens in September or October. The 2018 perihelion passage was the most favourable one since 1959 and until 2078, with a minimum geocentric distance (Δ) of 0.392 au (the peak visual magnitude was 7.0). This latest apparition was somewhat similar to the 1985 one, when the ICE spacecraft flew by the comet nucleus at a distance of 7800 km (Cowley 1985). In 1998, the shortest distance from Earth was 0.850 au (the peak visual magnitude was 8.5). The next relatively favourable apparition will be in September 2031 (Δ = 0.55 au). The orbit of comet 21P crosses Earth’s orbit and it is well known as being the parent of the Draconids meteor shower, which reached storm levels in 1933 and 1946.

Comet 41P/Tuttle-Giacobini-Kresák is a JFC orbiting the Sun every 5.4 yr with a low inclination (9°) orbit. It was first discovered by Horace Tuttle in 1858. The comet was observed by Michel Giacobini in 1907, but the link between the two apparitions was not established until 1928. Nevertheless, the comet was not recovered until 1951 by Lubor Kresák when its periodic nature was secured. 41P is a weakly active comet (radius <1 km, Lamy et al. 2004) and it has been mostly observed only during its favourable apparitions. Comet 41P is also well known for its frequent and dramatic outbursts: It displayed surges of brightness over a period of a few days, of 10 magnitudes amplitude twice in 1973, 1 magnitude in 1995, and 4–5 magnitudes twice in 2001. The 2017 return was the most favourable apparition in three centuries. The perihelion (1.05 au) lies just outside Earth’s orbit. The minimum geocentric distance was 0.142 au on 1 April 2017 when the comet reached a peak magnitude of 6.5. Optical observations in 2017 reported a dramatic change and slow-down of its rotation speed with a period increasing from 20 h to 46 h (Bodewits et al. 2018) or from 24 h to 48 h according to Schleicher et al. (2018) in only 1–1.5 months around perihelion. This suggests that the previous outbursts may have been associated with partial rotational break-up.

Comet 64P/Swift-Gehrels is a JFC orbiting the Sun every 9.4 yr on a low-inclination (8.9°) orbit. It was first discovered by L. Swift in 1889, rediscovered by T. Gehrels in 1973 and observed during each perihelion passage since. The most recent apparition perihelion was on 3.9 November 2018 UT at 1.393 au from the Sun. It made its closest approach ever to the Earth on 28 October 2018 at 0.445 au which provided one of the best opportunities for studying this comet remotely (peak magnitude of 9).

Comet 38P/Stephan-Oterma is classified as a Halley-type comet orbiting the Sun every 38 yr, but on a relatively low-inclination orbit (18°). Comet 38P was also first discovered in the 19th century (by E.J.M. Stephan but also J.E. Coggia and E.W.L. Tempel, in January 1867) and rediscovered in 1942 by L. Oterma, after one missed perihelion passage. The perihelion passage of 11 November 2018 at 1.589 au from the Sun was similarly favourable to the previous observed ones, with a perigee of 0.768 au on 17 December 2018 (peak magnitude of 9.5).

Comet 21P was the target of dedicated observing programmes with the Institut de RadioAstronomie Millimétrique (IRAM) 30-m telescope in November 1998 (perihelion on 21.3 November 1998 at 1.033 au from the Sun) and September 2018 (perihelion on 10.3 September 2018 at 1.013 au). It was also observed with the 15-m *James Clerk Maxwell Telescope* (JCMT) and the 10.4-m *Caltech submillimeter Observatory (CSO)* in November 1998. An international observing campaign was also set up for comet 41P and we observed it with the IRAM

30-m after perihelion (observations before perihelion suffered from bad weather). We observed comets 64P and 38P with the IRAM 30-m telescope at a time when the main target (comet 46P/Wirtanen) was unavailable on 12, 13, and 18 December 2018 UT. We searched for several molecules in those comets. We report the detections of HCN, CH$_3$OH, CS, CH$_3$CN, and H$_2$CO and tentatively H$_2$S in 21P; HCN and CH$_3$OH in 41P; HCN, CH$_3$OH, and CS in 64P; and a tentative detection of CH$_3$OH in 38P. Sections 2 and 3 present the observations and frequency coverage. Section 4 shows the spectra of the detected molecules. Section 5 presents the information extracted from the observations for analysis and so as to compute production rates. Section 6 provides the retrieved production rates and abundances or upper limits, which are discussed and compared to other comets in Sect. 7.

## 2. Observations with the IRAM 30-m telescope

Comet 21P was observed under good weather conditions (2 to 5 mm of precipitable water vapour – pwv) daily during the 11–15 November 1998 period. In 2018, during the 14–18 September allocated period, the weather was very poor, with high opacity (9–40 mm of pwv, most of the time between 10 and 16 mm). The observations were thus more concentrated on the low frequency bands of the EMIR receiver (2 and 3 mm bands, Carter et al. 2012). The comet was tracked using orbital elements JPL#K182/7¹ in 2018, where ~4'' ephemeris error was estimated afterwards.

Comet 41P observations were originally planned during two sessions, mid-March and at the end of April 2017, to observe its evolution from before to after perihelion. However, adverse weather conditions limited observations to only two nights (out of the eight scheduled) post perihelion on 26.0 April and 1.0 May 2017 UT, with average weather (3–6 mm pwv). 41P was tracked using the JPL#171/7 orbit solution. The difference in ephemeris computed with the JPL#171/18 solution of 2018 and those used by the IRAM new control system (NCS software) was negligible (<0.2''). A map of HCN(3–2) taken on 26 April yielded a final pointing error of less than 2''.

Comets 64P and 38P were observed under good weather conditions (~2 mm of precipitable water vapour) on 11.8 and 17.9 December 2018, to fair weather (pwv ~5 mm) for 38P on 12.1 and 13.1 December 2018 UT. 64P was tracked using the JPL#183/4 orbit solution (JPL#35 for 38P) to compute the position in real time. Pointing was regularly checked on bright sources (rms <1'') but ephemeris errors around 2'' both in RA and Dec were found afterwards. The two comets were observed for a short time on three different days with different frequency coverage (Table 1).

In 1998 we used 3 mm, 2 mm, and 1 mm 'B100', 'C150', and 'B230' newly installed single-side band (SSB) receivers in parallel, each with a maximum 512 MHz bandwidth, connected to filter banks and correlators. In 2017–2018, we used the EMIR (Carter et al. 2012) 3 mm, 2 mm, and 1 mm band receivers in 2SB mode which offer two instantaneous bandwidths of 8 GHz per polarisation. They were connected to the Fourier Transform Spectrometer (FTS) which covers the full bandwidths with 0.195 MHz sampling (corresponding to 0.22–0.24 km s⁻<fnref n="1" /> at 1 mm wavelength). The central part of the band was also connected to the VESPA autocorrelator. Some lines (e.g. CH$_3$OH($S_2$ – 4$_1$)E at 266 838.123 MHz, HCN(1–0) at 88 631.848 MHz) were also observed with VESPA using its

<footnote n="1">¹ https://ssd.jpl.nasa.gov/horizons.cgi</footnote>
<table><caption>Table 1. Log of millimetre observations.</caption><thead><tr><th>UT date<br>(yyyy/mm/dd.d-dd.d)</th><th>⟨r<sub>h</sub>⟩<br>(au)</th><th>⟨Δ⟩<br>(au)</th><th>Phase<sup>(a)</sup></th><th>Integ.<sup>(b)</sup><br>(min)</th><th>Mode<sup>(c)</sup></th><th>pwv<sup>(d)</sup><br>(mm)</th><th>Freq. range<br>(GHz)</th></tr></thead><tbody><tr><td colspan="8"><strong>21P/Giacobini-Zinner 1998: JCMT 15-m</strong></td></tr><tr><td>1998/10/30.21-30.31</td><td>1.080</td><td>0.925</td><td>58.8°</td><td>92</td><td>F</td><td>0.7</td><td>354.46-354.60,345.75-345.89<sup>(e)</sup></td></tr><tr><td>1998/10/31.22-31.31</td><td>1.076</td><td>0.920</td><td>59.0°</td><td>102</td><td>P</td><td>0.8</td><td>338.31-338.45,338.60-338.74</td></tr><tr><td>1998/11/01.13-01.14</td><td>1.073</td><td>0.916</td><td>59.3°</td><td>9</td><td>F</td><td>0.6</td><td>362.59-362.72,354.46-354.60<sup>(e)</sup></td></tr><tr><td>01.15-01.18</td><td>1.073</td><td>0.916</td><td>59.3°</td><td>40</td><td>P</td><td>0.7</td><td>362.58-362.73,354.46-354.60<sup>(e)</sup></td></tr><tr><td>01.20-01.22</td><td>1.072</td><td>0.916</td><td>59.3°</td><td>28</td><td>P</td><td>0.8</td><td>362.58-362.73</td></tr><tr><td>01.24-01.27</td><td>1.072</td><td>0.916</td><td>59.3°</td><td>17</td><td>F</td><td>0.9</td><td>342.82-342.95</td></tr><tr><td>01.27-01.32</td><td>1.072</td><td>0.916</td><td>59.3°</td><td>48</td><td>P</td><td>0.9</td><td>342.81-342.96</td></tr><tr><td>1998/11/02.11-02.21</td><td>1.069</td><td>0.912</td><td>59.5°</td><td>77</td><td>F</td><td>0.9</td><td>345.73-345.86</td></tr><tr><td>02.22-02.26</td><td>1.068</td><td>0.911</td><td>59.5°</td><td>45</td><td>P</td><td>1.0</td><td>345.72-345.87</td></tr><tr><td>02.27-02.31</td><td>1.068</td><td>0.911</td><td>59.5°</td><td>28</td><td>F</td><td>1.1</td><td>354.46-354.60,345.75-345.89<sup>(e)</sup></td></tr><tr><td>1998/11/03.17-03.29</td><td>1.065</td><td>0.907</td><td>59.8°</td><td>105</td><td>F</td><td>1.3</td><td>351.72-351.85</td></tr><tr><td>1998/11/04.13-04.20</td><td>1.062</td><td>0.903</td><td>60.0°</td><td>80</td><td>P</td><td>2.8</td><td>265.81-265.96,262.81-262.96<sup>(e)</sup></td></tr><tr><td>04.24-04.29</td><td>1.062</td><td>0.902</td><td>60.0°</td><td>44</td><td>F</td><td>1.9</td><td>350.86-350.96</td></tr><tr><td colspan="8"><strong>21P/Giacobini-Zinner 1998: CSO 10.4-m</strong></td></tr><tr><td>1998/11/10.17-10.27</td><td>1.046</td><td>0.879</td><td>61.2°</td><td>90.7</td><td>W</td><td>6.6</td><td>265.7-266.1, 262.7-263.1<sup>(e)</sup></td></tr><tr><td>1998/11/13.17-13.27</td><td>1.040</td><td>0.870</td><td>61.6°</td><td>85.3</td><td>W</td><td>6.3</td><td>265.7-266.1, 262.7-263.1<sup>(e)</sup></td></tr><tr><td>1998/11/14.16-14.23</td><td>1.039</td><td>0.868</td><td>61.8°</td><td>53.3</td><td>W</td><td>0.8</td><td>265.7-266.1, 262.7-263.1<sup>(e)</sup></td></tr><tr><td>11/14.23-14.31</td><td>1.039</td><td>0.867</td><td>61.8°</td><td>74.7</td><td>W</td><td>0.8</td><td>304.0-304.4, 307.0-307.4<sup>(e)</sup></td></tr><tr><td>1998/11/15.17-15.25</td><td>1.037</td><td>0.865</td><td>61.9°</td><td>85.3<sup>(f)</sup></td><td>W</td><td>3.2</td><td>304.0-304.4, 307.0-307.4<sup>(e)</sup></td></tr><tr><td colspan="8"><strong>21P/Giacobini-Zinner 1998: IRAM 30-m</strong></td></tr><tr><td>1998/11/11.60-11.83</td><td>1.043</td><td>0.874</td><td>61.4°</td><td>121.4</td><td>F</td><td>2.0</td><td>88.6, 168.8, 265.9</td></tr><tr><td>1998/11/12.58-12.70</td><td>1.041</td><td>0.872</td><td>61.6°</td><td>90</td><td>W</td><td>2.4</td><td>88.6, 168.8, 265.9</td></tr><tr><td>11/12.71-12.83</td><td>1.041</td><td>0.871</td><td>61.6°</td><td>90</td><td>F</td><td>2.7</td><td>88.6, 168.8, 265.9</td></tr><tr><td>1998/11/13.57-13.67</td><td>1.039</td><td>0.869</td><td>61.7°</td><td>80</td><td>W</td><td>5.0</td><td>88.6, 168.8, 265.9</td></tr><tr><td>11/13.68-13.83</td><td>1.039</td><td>0.869</td><td>61.7°</td><td>85.5</td><td>F</td><td>4.5</td><td>88.6, 168.8, 265.9</td></tr><tr><td>1998/11/14.55-14.81</td><td>1.038</td><td>0.866</td><td>61.8°</td><td>176.5</td><td>F</td><td>4.5</td><td>88.6, 168.8, 265.9</td></tr><tr><td>1998/11/15.55-15.81</td><td>1.037</td><td>0.864</td><td>61.9°</td><td>210</td><td>W</td><td>3.8</td><td>88.6, 145.1, 225.7</td></tr><tr><td colspan="8"><strong>21P/Giacobini-Zinner 2018: IRAM 30-m</strong></td></tr><tr><td>2018/09/14.22-14.28</td><td>1.014</td><td>0.394</td><td>77.6°</td><td>61</td><td>W</td><td>16-37</td><td>83.3-91.1, 141.6-149.4</td></tr><tr><td>14.38-14.40</td><td>1.015</td><td>0.394</td><td>77.5°</td><td>19</td><td>W</td><td>16-40</td><td>83.3-91.1, 141.6-149.4</td></tr><tr><td>2018/09/15.09-15.12</td><td>1.015</td><td>0.395</td><td>77.4°</td><td>28</td><td>W</td><td>15-18</td><td>248.7-256.5, 264.4-272.2</td></tr><tr><td>15.15-15.25</td><td>1.015</td><td>0.395</td><td>77.4°</td><td>98</td><td>W</td><td>19-23</td><td>83.3-91.1, 141.6-149.4</td></tr><tr><td>15.27-15.37</td><td>1.015</td><td>0.395</td><td>77.3°</td><td>79</td><td>W</td><td>19-27</td><td>83.3-91.1, 141.6-149.4</td></tr><tr><td>2018/09/16.10-16.22</td><td>1.016</td><td>0.397</td><td>77.1°</td><td>121</td><td>W</td><td>11-14</td><td>146.9-154.7, 162.6-170.4</td></tr><tr><td>16.39-16.40</td><td>1.017</td><td>0.397</td><td>77.1°</td><td>14</td><td>W</td><td>14</td><td>146.9-154.7, 162.6-170.4</td></tr><tr><td>16.43-16.46</td><td>1.017</td><td>0.397</td><td>77.1°</td><td>32.7</td><td>W</td><td>11</td><td>83.3-91.1, 225.1-232.8</td></tr><tr><td>2018/09/17.09-17.11</td><td>1.018</td><td>0.399</td><td>76.9°</td><td>23.3</td><td>W</td><td>9-10</td><td>83.3-91.1, 141.6-149.4</td></tr><tr><td>17.18-17.22</td><td>1.018</td><td>0.399</td><td>76.9°</td><td>42</td><td>W</td><td>9-10</td><td>83.3-91.1, 141.6-149.4</td></tr><tr><td>17.25-17.29</td><td>1.018</td><td>0.399</td><td>76.9°</td><td>46.7</td><td>W</td><td>9-10</td><td>83.3-91.1, 161.6-169.4</td></tr><tr><td>17.30-17.34</td><td>1.018</td><td>0.399</td><td>76.8°</td><td>46.7</td><td>W</td><td>12-27</td><td>83.3-91.1, 161.6-169.4</td></tr><tr><td>2018/09/18.10-18.22</td><td>1.019</td><td>0.401</td><td>76.6°</td><td>121.3</td><td>W</td><td>9</td><td>83.3-91.1, 161.6-169.4</td></tr><tr><td>18.24-18.35</td><td>1.019</td><td>0.402</td><td>76.5°</td><td>103</td><td>W</td><td>10-12</td><td>83.3-91.1, 218.2-226.0</td></tr><tr><td colspan="8"><strong>41P/Tuttle-Giacobini-Kresák: IRAM 30-m</strong></td></tr><tr><td>2017/04/26.00-26.04</td><td>1.061</td><td>0.179</td><td>67.7°</td><td>60</td><td>W</td><td>3.0</td><td>248.7-256.5, 264.4-272.2</td></tr><tr><td>26.07-26.09</td><td>1.061</td><td>0.179</td><td>67.6°</td><td>28.0</td><td>W</td><td>4.2</td><td>224.9-232.7, 240.5-248.3</td></tr><tr><td>2017/05/01.00-01.03</td><td>1.074</td><td>0.192</td><td>64.8°</td><td>36</td><td>W</td><td>5.5</td><td>248.7-256.5, 264.4-272.2</td></tr><tr><td>01.07-01.09</td><td>1.074</td><td>0.192</td><td>64.7°</td><td>18.7</td><td>W</td><td>6.0</td><td>146.9-154.7, 162.6-170.4</td></tr><tr><td colspan="8"><strong>64P/Swift-Gehrels: IRAM 30-m</strong></td></tr><tr><td>2018/12/11.75-11.76</td><td>1.465</td><td>0.577</td><td>26.6°</td><td>18.7</td><td>W</td><td>1.6</td><td>248.7-256.5, 264.4-272.2</td></tr><tr><td>2018/12/17.94-17.96</td><td>1.490</td><td>0.618</td><td>27.5°</td><td>23.3</td><td>W</td><td>2.2</td><td>240.4-248.1, 256.0-263.8</td></tr><tr><td colspan="8"><strong>38P/Stephan-Oterma: IRAM 30-m</strong></td></tr><tr><td>2018/12/12.11-12.12</td><td>1.636</td><td>0.768</td><td>23.9°</td><td>14.0</td><td>W</td><td>4.5</td><td>209.7-217.5, 225.4-233.1</td></tr><tr><td>2018/12/13.11-13.12</td><td>1.639</td><td>0.767</td><td>23.5°</td><td>10.3</td><td>W</td><td>5.3</td><td>146.9-154.7, 162.6-170.4</td></tr><tr><td>2018/12/17.92-17.93</td><td>1.655</td><td>0.766</td><td>21.5°</td><td>9.3</td><td>W</td><td>2.5</td><td>240.4-248.1, 256.0-263.8</td></tr></tbody></table>

Notes. (a) Phase angle. (b) Total (all setups included) integration time (ON+OFF) on the source. (c) Observing mode: F=frequency switching, W=Wobbler switching (nodding of the secondary mirror), P=position/beam switching. (d) Mean precipitable water vapour in the atmosphere above the telescope. (e) Spectra obtained in DSB mode with the two frequency domains overlapping in the same spectrum. (f) Likely unstable tuning, data not used.
higher spectral resolution (e.g. 20–40 kHz sampling) in a narrow
bandwidth (20–80 MHz).

In order to cancel the sky emission, a wobbler switching
mode at a frequency of 0.5 Hz was used for most observations,
with reference sky positions at ±180" in azimuth. In 1998 some
observations were obtained in frequency switching mode with
frequency throws of 7.2, 7.8, or 15.85 MHz.

3. Observations with JCMT and CSO

Comet 21P was also observed from Maunakea between 30 Octo-
ber and 14 November 1998. The weather was good on 30,
31 October, 1, 2, 3, and 14 November 1998 with less than 1 mm
pwv, average on 4 and 15 November (2 mm pwv) and poor for
Maunakea standards on 10 and 13 November (about 5 mm pwv).
With the exception of some observations of the HCN(3–2) line
at 265.9 GHz, most of these observations used the 0.8 mm band
receivers (304–363 GHz). Some of the JCMT observations were
obtained in the slow frequency switching mode, switching by
±8.1 MHz every two seconds. Other observations were obtained
by beam switching with a reference at 2' in azimuth. All CSO
observations and some of the JCMT observations were obtained
in double-side band (DSB) mode with the two side bands sep-
arated by 3 GHz and overlapping in the resulting spectrum.
Additional information is provided in Table 1. The tracking of
the comet was done by interpolating ephemeris generated with
the JPL-DE406#27 orbit solution.

4. Spectra and line intensities

Some of the spectra of 21P obtained in the frequency switch-
ing mode at JCMT or IRAM in 1998 are contaminated by
stratospheric HCN, H₂CO, or mesospheric CO. The telluric line
however falls at an absolute frequency different from the comet
line thanks to its Doppler shift. However the wings of the broader
stratospheric lines overlap with the comet spectra and we had
to remove a fit to the atmospheric line to obtain a flatter base-
line in the figures presented here and to derive accurate noise
and line intensities. Three lines of HCN, seven lines of CH₃OH,
and one each of CS and H₂CO were detected with the JCMT,
CSO, and IRAM telescopes in 1998. Spectra (average of several
days of observations in some cases) are shown in Figs. 1, 2, 3,
and 4. The daily intensities of the lines observed with JCMT
and CSO are provided in Tables A.1 and A.2, respectively. The
intensities of the lines observed with the IRAM 30-m telescope
in 1998 and 2018 are provided in Table A.3. The mean Doppler
shift of each line is computed as the first moment integral over
the line width (∫Tb.v.dv/∫Tb.dv), excluding the weaker (<4%)
hyperfine components for HCN, and taking the average on the
two strongest hyperfine components for HCN(1–0). In 2018,
lines were stronger because the comet was closer to Earth. The
EMIR receivers at the IRAM 30-m telescope are also more sen-
sitive than the receivers used in 1998, and offer a much broader
bandwidth (Table 1). Spectra are shown in Figs. 5, 6, and 7.

Spectra of HCN and CH3OH detected in comet 41P are
shown in Figs. 8 and 9 and integrated intensities are provided
in Table A.4.

For comet 64P, the HCN(3–2) was clearly detected on
11.8 December (Fig. 10) together with a marginal detection of
methanol lines at 266.838 GHz and around 252 GHz. CH3OH
was best detected via its 242 GHz lines (Fig. 11) on 17.9 Decem-
ber UT. CS(5–4) was marginally seen in this same observation.
Integrated intensities or 3 − σ upper limits are provided in

<figure><img src="image_1.png" /><figcaption>Fig. 1. Lines observed with the JCMT 15-m in comet 21P: HCN(4–3) at 354 505.472 MHz between 30 October and 2 November 1998, and HCN(3–2) at 265 886.434 MHz on 4.2 November 1998 UT. Negative ‘ghost’ lines appear at ±13.7 km s⁻¹ (16.2 MHz) in the HCN(4–3) spectrum owing to the frequency switching technique. The position and relative intensity of the hyper-fine components are drawn below the spectra. The vertical scale is the main beam brightness temperature and the horizontal scale is the Doppler velocity in the comet rest frame.</figcaption></figure>

<caption>Table A.5. No individual line was clearly detected in comet 38P/Stephan-Oterma, owing to very limited integration times.</caption>

5. Data analysis

In order to derive production rates, we assume steady-state out-
gassing at a constant velocity and use the Haser density profile.
Excitation of the molecular rotational lines depends on radiative
and collisional processes. Collisions with neutrals and electrons
are taken into account but depend on the gas temperature and
water production rate. Model parameters are given in Table 3.

5.1. Expansion velocity and outgassing pattern

The expansion velocity is derived from the lines observed with
the best signal-to-noise ratio (S/N). Generally the Doppler veloc-
ity at half maximum intensity (VHM) in the comet frame pro-
vides a good estimate of the expansion velocity: on the observer
side from negative velocities and on the opposite side (gas mov-
ing away seen on the redshifted side of the line). When possible
for a double-peak line we fitted two Gaussians from which we
measure the two *VHMs*. We took the average value of the expan-
sion velocity deduced from the blueshifted and redshifted sides
of the lines.
<figure><img src="image_2.png" /><figcaption>**Fig. 2.** Lines observed with the JCMT 15-m in comet 21P: two CH₃OH lines at 338 GHz on 31.3 October and one line at 350 905.100 MHz on 4.3 November 1998 UT (**left side**); H₂CO at 351 768.639 MHz on 3.2 November and CS at 342 882.850 MHz on 1.3 November 1998 UT (**right side**). The vertical scale is the main beam brightness temperature and the horizontal scale is the Doppler velocity in the comet rest frame. The frequency scale for the CH₃OH lines at 338 GHz is also provided on the upper axis.</figcaption></figure>

For 21P average values are given in Table 3 for the various periods. The lines do not show a strong asymmetry (Figs. 1–7, especially looking at the very symmetric HCN(1–0) line observed with high frequency resolution in 2018 in Fig. 7). The Doppler shift is small or not significant (Table A.1–A.3, Col. 7). Therefore we assumed symmetric outgassing at a mean velocity of 0.7 to 0.73 km s⁻¹ for 21P.

The expansion velocity we assumed for comet 41P was derived from the profile of HCN(3–2) (Fig. 8) and the CH₃OH (2₀ − 1₋₁)E and (5₂ − 4₋₁)E lines, which were observed with the high-resolution backend VESPA. The hyperfine structure of the HCN(3–2) line has a minor effect on the broadening of the main component and was taken into account to derive the expansion velocity. From their line shape, we obtained the average values VHM = −0.8 and +0.6 km s⁻¹, and consequently we assumed a mean expansion velocity of 0.7 km s⁻¹.

For comet 64P, from the profile of the spectra of HCN(3–2) and CH₃OH(5₀ − 4₀)A at 241 791.367 MHz which have the highest *S*/*N* (10 and 6, respectively) we determined a gas velocity of 0.80 ± 0.05 km s⁻¹ on the Earth-facing side, and 0.5 ± 0.1 km s⁻¹ on the other side. Since the phase angle (Table 1) was small (less than 28°), in contrary to the two previous comets, the Earth-facing side was close to the day side and line asymmetry can be associated with day-to-night asymmetry in activity. We adopted

the mean value $v_{\text{exp}} = 0.65$ km s⁻¹. The mean Doppler shift of the lines (Table A.5, Col. 7) was also significantly negative, confirming the preferential outgassing on the day side at a higher velocity.

For comet 38P, which was observed at 1.6 au from the Sun, we assumed $v_{\text{exp}} = 0.6$ km s⁻¹, using the general ~0.8 × $r_h^{-0.5}$-dependence of the expansion velocity observed in comets (Biver et al. 2000, 2006).

## 5.2. Gas temperature

Rotational temperatures were measured for methanol and CH₃CN. In some cases a meaningful rotational temperature of HCN could also be deduced from contemporaneous observation of HCN(1–0) with IRAM 30-m and HCN(3–2) with CSO 10.4-m, which probed the same region of the coma (same beam size of ~26′′). Nevertheless the rotational temperature is not necessarily a measurement of the gas temperature *T*<sub>gas</sub> of the inner coma. As radio observations often sample a region larger than the collision-dominated region of the coma, the population of the rotational levels of the considered molecule evolves radiatively throughout the beam owing to spontaneous decay and infrared pumping via vibrational bands. Some series of methanol lines around 165, 157, 252, or 305 GHz sample rotational levels from
<figure><img src="image_3.png" /></figure>

**Fig. 3.** Lines observed with the CSO 10.4-m in comet 21P: HCN(3–2) at 265 886.434 MHz on 10–14 November 1998 UT, and CH₃OH lines at 307 165.911 MHz (Upper Side Band, left upper axis scale) and 304 208.324 MHz (Image and Lower Side Band, right upper axis scale) on 14.3 November 1998 UT. The receiver tuning was adjusted in order to get those two methanol line in the same spectrum separated by 20 MHz. The vertical scale is the main beam brightness temperature and the horizontal scale is the Doppler velocity in the comet rest frame. The frequency range of 40.7 MHz for the CH₃OH spectrum corresponds to a 40 km s⁻¹ Doppler velocity window as for HCN. A velocity scale centred on each line has been plotted on the lower axis of the methanol spectrum.

the same K-ladder that evolve in a similar manner and keep memory of the thermal equilibrium in the collision-dominated region (Bockelée-Morvan et al. 1994). On the other hand, the rotational levels sampled by methanol lines around 145, 242, or 338 GHz, as well as HCN, are more sensitive to the collision rate and evolve rapidly towards a much colder fluorescence equilibrium. Consequently, their average rotational temperature in the beam is much lower than Tgas. Measured values of Trot and constraints derived on the gas temperature Tgas are provided in Table 2.

Several methanol lines were detected in 1998 and in 2018 in comet 21P. In 1998, however, the constraints remain loose on the gas temperature owing to low S/N (Table 2). We adopted a value of Tgas = 40 K. In 2018 the rotational temperature of the methanol lines at 145 and 165 GHz (Fig. 12) could be more precisely measured and we adopted Tgas = 42 K.

For comet 41P, the individual methanol lines were often only marginally detected (Table A.4), so the uncertainty in the rotational temperature is relatively large: T<sub>rot</sub> = 33<sup>+inf</sup><sub>-19</sub> K, T<sub>rot</sub> = 28 ± 7 K, and T<sub>rot</sub> = 51<sup>+53</sup><sub>-17</sub> K for the 165, 252, and 242 GHz lines respectively. The rotational temperature of the 165 and 252 GHz lines are expected to provide a good estimate of T<sub>gas</sub>,

<figure><img src="image_4.png" /></figure>

**Fig. 4.** HCN(3–2) line in comet 21P observed with the IRAM 30-m on 11.7–14.7 October 1998 and HCN(1–0) observed on 11.7–15.7 October 1998. Scales as in Fig. 1.

<figure><img src="image_5.png" /></figure>

**Fig. 5.** CH₃OH lines at 145 GHz in comet 21P observed with the IRAM 30-m on 14.2–17.2 September 2018 UT. The vertical scale is the main beam brightness temperature and the horizontal scale is the Doppler velocity in the comet rest frame and the rest frequency in MHz is provided on the upper axis.

whereas we expect a lower value for the 242 GHz lines. Hence, we used T<sub>gas</sub> = 45 K, but also tested a value of 35 K, which yields similar production rates (3–8% lower for both HCN and CH₃OH).

We had only weak constraints on the gas temperature from the observations of methanol lines in 64P. The rotational temperature of the lines around 252 GHz is T<sub>rot</sub> = 47<sup>+51</sup><sub>-19</sub> K, which suggests a gas temperature T<sub>gas</sub> in the 33–100 K range. The rotational temperature of the 242 GHz lines is better constrained
<figure><img src="image_6.png" /><figcaption>**Fig. 6.** CH₃OH lines at 165 GHz in comet 21P observed with the IRAM 30-m on 16.1–18.2 September 2018 UT. Scales as in Fig. 5.</figcaption></figure>

(Fig. 13) with $T_{\text{rot}} = 18 \pm 5$ K. Unfortunately it does not bring a precise constraint on the gas temperature: for the assumed collision rate, this rotational temperature is compatible with any value of $T_{\text{gas}}$ in the 10–100 K range. We adopted $T = 40$ K, consistent with the values measured in comets of similar activity and heliocentric distance. In any case, a change of $\pm 10$ K does not result in a change of production rate larger than the uncertainties owing to the limited signal-to-noise of the dataset, and much larger changes in temperature have never been seen in comets of similar activity. We assumed a lower temperature $T_{\text{gas}} = 30$ K for comet 38P, typical for comets of moderate activity observed at 1.6 au from the Sun (e.g. Biver et al. 2012).

### 5.3. Reference water production rate

OH was observed in comet 21P with the Nançay radio telescope in 1985 (Gerard et al. 1988), 1998 and 2018. The 1998 data provided OH production rates in the range $3-4 \times 10^{28}$ molec. s$^{-1}$ for the period of the millimetre observations (Crovisier et al. 2002), whereas Combi et al. (2011) provided slightly higher production rates, with a value around $5 \times 10^{28}$ molec. s$^{-1}$ at perihelion (21 November 1998). In 2018, no reliable OH production rate could be derived from Nançay data during the IRAM observations because of the poor OH maser inversion at that time. Combi et al. (2011) also showed a trend observed at each perihelion passage of comet 21P: the activity peaks about 25 days prior to perihelion, with a steep decrease afterwards. This was also observed for the visual heliocentric magnitude, which peaked 1–2 weeks prior to perihelion in 1998 and 2018. In 2018, Roth et al. (2020) observed OH and H₂O in the infrared yielding a water production rate in the range 3.2–3.0 × $10^{28}$ molec. s$^{-1}$ for the 7–11 September period, just before our IRAM observations. Faggi et al. (2019), also from infrared observations of H₂O, found water production rate in the range 2.6–2.8 × $10^{28}$ molec. s$^{-1}$ for the 8 September–7 October 2018 period. We thus assumed $Q_{\text{H}_2\text{O}} = 3.0 \times 10^{28}$ molec. s$^{-1}$ at the time of IRAM observations (Table 3).

Comet 41P outgassing was monitored by several observatories during its 2017 approach, including the Solar Wind ANistropies (SWAN) instrument on the SOHO spacecraft (31 March to 20 April, Combi et al. 2019), TRAPPIST (11 March to 24 May for OH, Moulane et al. 2018) and with the Nançay radio telescope (4 March–5 June). The comet was not observed at Nançay at the time of the IRAM observations because of poor OH inversion. OH production rates of $\approx 1.1 \times 10^{28}$ and $0.6 \times 10^{28}$ molec. s$^{-1}$ were obtained for the 28 March–3 April

and 18–30 May, respectively. Most of the observations showed a plateau of activity $\pm 20$ days around perihelion with $Q_{\text{H}_2\text{O}} \approx 4, 3$ and $8 \times 10^{27}$ molec. s$^{-1}$ respectively. We adopted a median value of $Q_{\text{H}_2\text{O}} = 4 \times 10^{27}$ molec. s$^{-1}$ for the time of our observations.

For comets 64P and 38P, around the time of our observations, water production rates were estimated from ground-based spectro-photometric observations of OH. D.G. Schleicher (priv. comm.) provided $Q_{\text{H}_2\text{O}} = 7.6 \times 10^{27}$ molec. s$^{-1}$ for 64P on 16 December 2018. For comet 38P his values were $Q_{\text{H}_2\text{O}} = 5 - 4 \times 10^{27}$ molec. s$^{-1}$ on 3 and 16 December 2018. We adopted water production rates of $Q_{\text{H}_2\text{O}} = 0.8 \times 10^{28}$ molec. s$^{-1}$ for 64P and $Q_{\text{H}_2\text{O}} = 0.5 \times 10^{28}$ molec. s$^{-1}$ for 38P, both to compute collision rates and determine molecular abundances. 64P and 38P were only observed in September and October at Nançay.

## 6. Production rates and abundances

Line intensities have been converted to production rates using our time-dependent non-LTE code taking into account collisions with neutrals and electrons as well as solar radiative pumping via vibrational bands for HCN, CH₃OH, CS, H₂CO, CO, and CH₃CN, as used in previous papers (Biver et al. 2000, 2006, 2012, 2014, 2015, 2016). Photolysis constants are also from previous papers and detailed in Biver et al. (2021). Production rates or upper limits are provided in Table 4. For species with multiple lines observed at the same time, the weighted average is provided.

### 6.1. Relative abundances

Abundances relative to water were computed using the average water production rates provided in Table 3. For 41P we assumed a lower water production rate ($0.3 \times 10^{28}$ molec. s$^{-1}$) for the second day (1.0 May 2017) because HCN and CH₃OH showed a significant decrease in production rate and Moulane et al. (2018) also found a decrease in the water production rate between 26.1 April and the 28 April–3 May period. For 38P, abundances relative to water assume $Q_{\text{H}_2\text{O}} = 0.5 \times 10^{28}$ molec. s$^{-1}$. Although no line is detected with a signal-to-noise larger than 2, the combined signal of the 13 strongest methanol lines at 165–170 GHz and 242 GHz yields a formal signal-to-noise of 3.5 and a realistic methanol abundance of $5.0 \pm 1.4\%$. For comet 21P we took the average of 1998 and 2018 measurements as they do not show any difference.

Abundances and a comparison to values measured in other comets are provided in Table 5. Range of abundances measured in the full sample of comets and in other JFCs in the radio is also provided in the last two columns (Biver & Bockelée-Morvan 2019; Biver et al. 2021).

### 6.2. Evolution of production rates of comet 21P

Comet 21P was observed during a two-week period in 1998 and again in 2018 after three full orbits around the Sun. Figure 14 shows both series of production rates measurements relative to the time of perihelion. The 1998 data show a decreasing trend of the production rates of methanol and HCN over the two weeks. The 4-day span of the 2018 observation is too short to draw definitive conclusions, but shows production rates very similar to the 1998 values. The heliocentric distance in 2018 was on average 4% lower than in 1998. If the peak outgassing follows $r_h^{-2}$, it should have been 8% higher in 2018. Taking this change into account, one can infer that the comet activity follows the same trend as in 1998 of continuous decrease of active surface around
<figure><img src="image_7.png" /><figcaption>**Fig. 7.** Observations of comet 21P with the IRAM 30-m in 2018: HCN(1–0) on 14.2–17.8 September, CS(3–2) and CH<sub>3</sub>CN on 14.2–17.2 September and H<sub>2</sub>S(1<sub>10</sub> – 1<sub>01</sub>) on 16.1–18.2 September UT. Scales as in Fig. 1. All the spectra are plotted with the same intensity scale and Doppler velocity range (40 km s<sup>-1</sup>). The rest frequency in MHz is provided on the upper axis of CH<sub>3</sub>CN lines.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 8.** HCN(3–2) line in comet 41P observed with the IRAM 30-m on 26–31 April 2017. Scales as in Fig. 1.</figcaption></figure>

perihelion. This was also observed by Combi et al. (2011) and Roth et al. (2020) for the 1998, 2005, and 2018 perihe- lions. Besides this small trend, production rates and relative abundances were very similar in 1998 and 2018, suggesting no significant secular trend both in activity and composition.

**7. Discussion and conclusion**

Comet 64P/Swift-Gehrels had its most favourable apparition in 2018 and was poorly investigated at previous perihelion

<figure><img src="image_9.png" /><figcaption>**Fig. 9.** CH<sub>3</sub>OH lines at 242 GHz in comet 41P observed with the IRAM 30-m on 26.1 April 2017. The vertical scale is the main beam brightness temperature and the horizontal scale is the Doppler velocity in the comet rest frame and the rest frequency in MHz is provided on the upper axis.</figcaption></figure>

passages. It attracted special attention as it started its pre- perihelion activity in August 2018 with a few outbursts (Kelley et al. 2019). Regarding its composition based on optical spec- troscopy, preliminary TRAPPIST (E. Jehin, priv. comm.) values obtained in January 2019 yield log[(Q(C<sub>2</sub>)/Q(CN)) = -0.02. For comet 38P, the value is +0.10 (TRAPPIST) to +0.01 near its previous perihelion passage (A’Hearn et al. 1995). For 41P,
<figure><img src="image_10.png" /><figcaption>**Fig. 10.** HCN(3–2) line in comet 64P/Swift-Gehrels observed with the IRAM 30-m on 11.8 December 2018. Scales as in Fig. 1.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 11.** CH<sub>3</sub>OH lines at 242 GHz in comet 64P/Swift-Gehrels observed with the IRAM 30-m on 17.9 December 2018. Scales as in Fig. 9.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 12.** Rotational diagram of the methanol lines at 165 GHz in comet 21P. We plot the logarithm of a quantity which should be proportional to the average upper level population in the beam (for optically thin line, Bockelée-Morvan et al. 1994) against the energy of the upper level, for each transition. At thermal equilibrium at *T* all points should align on a line of slope 1/*T*. The black dots with errorbars are the observed values, with upper level of the corresponding transition labelled above. The red line shows the best fit to the observations, with dashed lines for the ±1σ values of the rotational temperature given above the plot. Green circles show the predicted values for a model with a *T*<sub>gas</sub> = 42 K implying a rotational temperature of 41.7 K.</figcaption></figure>

log[Q(C<sub>2</sub>)/Q(CN)] = 0.06 from Moulane et al. (2018). These three comets are thus classified as ‘typical’ according to the classification by A’Hearn et al. (1995) based on optical spectroscopy. Comet 21P, on the other hand, is classified as carbon-chain depleted with log(Q(C<sub>2</sub>)/Q(CN)) = -0.67 (A’Hearn et al. 1995). Our millimetre-wave observations suggest that ‘typical’ comets

<table><caption>**Table 2.** Rotational and inferred gas temperatures.</caption><thead><tr><th>UT date<br>(mm/dd,d)</th><th>Molecule</th><th>Freq.<br>(GHz)</th><th>Lines <sup>(a)</sup></th><th><i>T</i><sub>rot</sub><br>(K)</th><th><i>T</i><sub>gas</sub><br>(K)</th></tr></thead><tbody><tr><td colspan="6"><strong>21P/Giacobini-Zinner (1998)</strong></td></tr><tr><td>10/31.3</td><td>CH<sub>3</sub>OH</td><td>338</td><td>4</td><td>19 ± 9</td><td>&gt;45</td></tr><tr><td>11/14.3</td><td>CH<sub>3</sub>OH</td><td>305</td><td>2</td><td>38<sup>+49</sup><sub>-16</sub></td><td>&gt;45</td></tr><tr><td>11/11–15</td><td>HCN</td><td>89–266</td><td>2</td><td>16.4 ± 2.1</td><td>30–40 <sup>(b)</sup></td></tr><tr><td colspan="6"><strong>21P/Giacobini-Zinner (2018)</strong></td></tr><tr><td>09/14–18</td><td>CH<sub>3</sub>CN</td><td>147</td><td>3</td><td>33 ± 28</td><td>33 ± 28</td></tr><tr><td>09/14–18</td><td>CH<sub>3</sub>OH</td><td>145</td><td>10</td><td>41 ± 9</td><td>&gt;32</td></tr><tr><td>09/14–18</td><td>CH<sub>3</sub>OH</td><td>165</td><td>9</td><td>38 ± 7</td><td>39 ± 7</td></tr><tr><td colspan="6"><strong>41P/Tuttle-Giacobini-Kresák (2017)</strong></td></tr><tr><td>04/26.1</td><td>CH<sub>3</sub>OH</td><td>242</td><td>11</td><td>51<sup>+53</sup><sub>-19</sub></td><td>&gt;40</td></tr><tr><td>04/27.4</td><td>CH<sub>3</sub>OH</td><td>252</td><td>7 × 2</td><td>28 ± 7</td><td>22–38</td></tr><tr><td>05/01.1</td><td>CH<sub>3</sub>OH</td><td>165</td><td>6</td><td>33<sup>+∞</sup><sub>-19</sub></td><td>&gt;15</td></tr><tr><td colspan="6"><strong>64P/Swift-Gehrels (2018)</strong></td></tr><tr><td>12/11.7</td><td>CH<sub>3</sub>OH</td><td>252</td><td>5 × 2</td><td>47 ± 35</td><td>33–100</td></tr><tr><td>12/17.9</td><td>CH<sub>3</sub>OH</td><td>242</td><td>9</td><td>18 ± 5</td><td>Any</td></tr></tbody></table>

**Notes.** <sup>(a)</sup>Number of lines used for the determination of *T*<sub>rot</sub>. <sup>(b)</sup>For an electron density factor *x*<sub>ne</sub> = 0.2 (Zakharov et al. 2007).

<table><caption>**Table 3.** Adopted model parameters.</caption><thead><tr><th>Comet</th><th>Date</th><th><i>Q</i><sub>H<sub>2</sub>O</sub><br>(molec. s<sup>-1</sup>)</th><th><i>v</i><sub>exp</sub><br>(km s<sup>-1</sup>)</th><th><i>T</i><sub>gas</sub><br>(K)</th></tr></thead><tbody><tr><td>21P</td><td>0–4 November 1998</td><td>4 × 10<sup>28</sup></td><td>0.7</td><td>40</td></tr><tr><td>21P</td><td>10–15 November 1998</td><td>3 × 10<sup>28</sup></td><td>0.7</td><td>40</td></tr><tr><td>21P</td><td>14–18 September 2018</td><td>3 × 10<sup>28</sup></td><td>0.73</td><td>42</td></tr><tr><td>41P</td><td>26–31 April 2017</td><td>0.5 × 10<sup>28</sup></td><td>0.7</td><td>35</td></tr><tr><td>64P</td><td>12–18 December 2018</td><td>0.8 × 10<sup>28</sup></td><td>0.65</td><td>40</td></tr><tr><td>38P</td><td>12–18 December 2018</td><td>0.5 × 10<sup>28</sup></td><td>0.6</td><td>30</td></tr></tbody></table>

have a higher methanol-to-water abundance ratio (Table 5), which appears to be often correlated with a high C<sub>2</sub>/OH (and C<sub>2</sub>/CN for comets studied here).

Comet 21P is severely depleted in H<sub>2</sub>S; it has the lowest H<sub>2</sub>S abundance relative to water observed so far, and the same is true for CS. Calmonte et al. (2016) found that the sulphur budget in comet 67P/Churyumov-Gerasimenko led to S/O ~ 1.5% in the volatile phase, with H<sub>2</sub>S as a major contributor. If the ratio is similar in comet 21P, it implies that sulphur species could represent more than 1% relative to water (assuming CO<sub>2</sub>/H<sub>2</sub>O ≤25%, based on e.g. Ootsubo et al. 2012, for the O budget). In most other comets investigated at radio wavelengths (Biver & Bockelée-Morvan 2019) the total of sulphur species amounts to 1–2% relative to water (in number). Unless comet 21P is depleted in sulphur in its volatile phase, one may wonder if OCS, SO(SO<sub>2</sub>), H<sub>2</sub>CS, CH<sub>3</sub>SH, or some other molecule, could be major sulphur species in its coma. This cannot be completely excluded by the upper limits we obtained (see Tables 4 and 5 for the former species and we don’t have stringent upper limits for the others). However Saki et al. (2020) derived OCS/H<sub>2</sub>O = 0.11 ± 0.02% in comet 21P.

Besides CH<sub>3</sub>CN, CS, and HCN, all other molecules (CH<sub>3</sub>OH, H<sub>2</sub>CO, H<sub>2</sub>S, CO) show abundances towards the low
<table><caption>**Table 4.** Production rates.</caption><thead><tr><th>UT date<br>(yyyy/mm/dd)</th><th>Molecule</th><th>r<sub>h</sub><br>(au)</th><th>Production rate<br>(10<sup>25</sup> molec. s<sup>-1</sup>)</th><th>Lines<sup>(a)</sup></th></tr></thead><tbody><tr><td colspan="5"><strong>21P/Giacobini-Zinner</strong></td></tr><tr><td rowspan="6">1998/10/30–35</td><td>HCN</td><td rowspan="6">1.071</td><td>3.5 ± 0.5</td><td>2</td></tr><tr><td>CH<sub>3</sub>OH</td><td>71 ± 10</td><td>5</td></tr><tr><td>CS</td><td>3.3 ± 0.5</td><td>1</td></tr><tr><td>H<sub>2</sub>CO<sup>(b)</sup></td><td>5.5 ± 2.1</td><td>1</td></tr><tr><td>HNC</td><td>&lt;0.45</td><td>(1)</td></tr><tr><td>CO</td><td>&lt;100</td><td>(1)</td></tr><tr><td rowspan="4">1998/11/10–16</td><td>HCN</td><td rowspan="4">1.044</td><td>2.5 ± 0.2</td><td>2</td></tr><tr><td>CH<sub>3</sub>OH</td><td>41 ± 10</td><td>5</td></tr><tr><td>H<sub>2</sub>CO<sup>(b)</sup></td><td>&lt;11.5</td><td>(1)</td></tr><tr><td>H<sub>2</sub>S</td><td>&lt;10.8</td><td>(1)</td></tr><tr><td rowspan="14">2018/09/14–18</td><td>HCN</td><td rowspan="14">1.017</td><td>3.4 ± 0.1</td><td>2</td></tr><tr><td>CH<sub>3</sub>OH</td><td>53.9 ± 3.3</td><td>15</td></tr><tr><td>CH<sub>3</sub>CN</td><td>0.78 ± 0.10</td><td>4</td></tr><tr><td>H<sub>2</sub>CO<sup>(b)</sup></td><td>5.5 ± 2.2</td><td>(4)</td></tr><tr><td>H<sub>2</sub>S</td><td>2.8 ± 1.0</td><td>1</td></tr><tr><td>CS</td><td>2.3 ± 0.2</td><td>1</td></tr><tr><td>CO</td><td>&lt;489</td><td>(1)</td></tr><tr><td>HNC</td><td>&lt;0.34</td><td>(1)</td></tr><tr><td>H<sub>3</sub>N</td><td>&lt;0.49</td><td>(6)</td></tr><tr><td>OCS</td><td>&lt;17.8</td><td>(2)</td></tr><tr><td>SO<sup>(c)</sup></td><td>&lt;13.7</td><td>(1)</td></tr><tr><td>HNCO</td><td>&lt;2.6</td><td>(2)</td></tr><tr><td>HCOOH</td><td>&lt;18.8</td><td>(3)</td></tr><tr><td>NH<sub>2</sub>CHO</td><td>&lt;1.3</td><td>(8)</td></tr><tr><td>CH<sub>3</sub>CHO</td><td>&lt;11.9</td><td>(6)</td></tr><tr><td colspan="5"><strong>41P/Tuttle-Giacobini-Kresák</strong></td></tr><tr><td rowspan="3">2017/04/26.0</td><td>HCN</td><td rowspan="3">1.061</td><td>0.38 ± 0.03</td><td>1</td></tr><tr><td>CH<sub>3</sub>OH</td><td>15.3 ± 1.7</td><td>(22)</td></tr><tr><td>CO</td><td>&lt;85</td><td>(1)</td></tr><tr><td rowspan="5">2017/05/01.0</td><td>HCN</td><td rowspan="5">1.074</td><td>0.29 ± 0.06</td><td>1</td></tr><tr><td>CH<sub>3</sub>OH</td><td>6.6 ± 3.5</td><td>(22)</td></tr><tr><td>CH<sub>3</sub>CN</td><td>&lt;0.35</td><td>(8)</td></tr><tr><td>H<sub>3</sub>N</td><td>&lt;0.39</td><td>(2)</td></tr><tr><td>H<sub>2</sub>S</td><td>&lt;2.2</td><td>(1)</td></tr><tr><td rowspan="10">2017/04/26–31</td><td>HCN</td><td rowspan="10">1.065</td><td>0.36 ± 0.03</td><td>1</td></tr><tr><td>CH<sub>3</sub>OH</td><td>14.0 ± 1.5</td><td>(30)</td></tr><tr><td>HNCO</td><td>&lt;1.0</td><td>(2)</td></tr><tr><td>H<sub>2</sub>CO<sup>(b)</sup></td><td>&lt;5.2</td><td>(2)</td></tr><tr><td>HNC</td><td>&lt;0.074</td><td>(1)</td></tr><tr><td>CS</td><td>&lt;0.48</td><td>(2)</td></tr><tr><td>OCS</td><td>&lt;8.9</td><td>(4)</td></tr><tr><td>SO<sup>(c)</sup></td><td>&lt;8.2</td><td>(1)</td></tr><tr><td>SO<sub>2</sub></td><td>&lt;2.5</td><td>(14)</td></tr><tr><td colspan="5"><strong>64P/Swift-Gehrels</strong></td></tr><tr><td rowspan="2">2018/12/11.75</td><td>HCN</td><td rowspan="2">1.465</td><td>1.20 ± 0.12</td><td>1</td></tr><tr><td>CH<sub>3</sub>OH</td><td>38 ± 6</td><td>11</td></tr><tr><td rowspan="5">2018/12/17.95</td><td>CH<sub>3</sub>OH</td><td rowspan="5">1.490</td><td>28 ± 2</td><td>10</td></tr><tr><td>CS</td><td>0.9 ± 0.3</td><td>1</td></tr><tr><td>SO</td><td>&lt;3.8</td><td>(2)</td></tr><tr><td>SO<sup>(c)</sup></td><td>&lt;11.2</td><td>(2)</td></tr><tr><td>HNCO</td><td>&lt;1.1</td><td>(1)</td></tr><tr><td colspan="5"><strong>38P/Stephan-Oterma</strong></td></tr><tr><td>2018/12/12–18</td><td>CH<sub>3</sub>OH</td><td>1.65</td><td>25 ± 7</td><td>(13)</td></tr><tr><td>2018/12/13.12</td><td>H<sub>2</sub>S</td><td>1.639</td><td>&lt;13.9</td><td>(1)</td></tr><tr><td>2018/12/13–18</td><td>CS</td><td>1.64</td><td>&lt;1.8</td><td>(2)</td></tr><tr><td>2018/12/12–13</td><td>H<sub>2</sub>CO<sup>(d)</sup></td><td>1.637</td><td>&lt;2.9</td><td>(3)</td></tr><tr><td>2018/12/12.12</td><td>CO</td><td>1.636</td><td>&lt;378</td><td>(1)</td></tr></tbody></table>

Notes. <sup>(a)</sup>Number of lines used for the determination of Q, in parenthesis if not all detected individually. <sup>(b)</sup>Assumed parent scale length of 5800 km. <sup>(c)</sup>SO assumed to be the daughter product of SO<sub>2</sub>. <sup>(d)</sup>Assumed parent scale length of 8000 km.

<figure><img src="image_13.png" /><figcaption>**Fig. 13.** Rotational diagram of the methanol lines at 242 GHz in comet 64P/Swift Gehrels. Scales are as is Fig. 12. The green circles show the predicted values for a model with a gas temperature $T_{gas} = 50$ K which implies a lower rotational temperature of 19.5 K.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Fig. 14.** Evolution of the production rates of comet 21P. Empty symbols (left side) are production rates measured in October–November 1998 a few days before perihelion ($r_h = 1.034$ on 21.3 November 1998). Filled symbols correspond to measurements of 2018, shortly after perihelion ($r_h = 1.013$ on 10.3 September 2018). Downward pointing empty triangles are upper limits (for CO in dark blue, H<sub>2</sub>S in red and HNC in blue) for the corresponding period. The vertical scale is the production rate in molec. s<sup>-1</sup>, and the horizontal scale the date of observation relative to the relevant perihelion date.</figcaption></figure>

end observed in comets (Table 5; Bockelée-Morvan & Biver 2017). The relatively low abundance of methanol seen in the radio around perihelion was also derived in the infrared in 1998 (Weaver et al. 1999), 2005 (DiSanti et al. 2013) and 2018 (Faggi et al. 2019), with abundances around 1.3% relative to water.

Comet 21P is also depleted in NH (A’Hearn et al. 1995) and NH<sub>2</sub> (Fink 2009) but showed unidentified emission near 8, 9 and 11 µm, some of which could be attributed to complex organics
<table><caption>**Table 5.** Abundance of organic molecules detected in comets.</caption><thead><tr><th rowspan="2">Molecule</th><th rowspan="2">Name</th><th colspan="6">Abundance relative to water in %</th></tr><tr><th>in 21P</th><th>in 41P</th><th>in 38P</th><th>in 64P</th><th>in comets</th><th>in JFCs</th></tr></thead><tbody><tr><td>HCN</td><td>Hydrogen cyanide</td><td>0.10 ± 0.01</td><td>0.09 ± 0.01</td><td></td><td>0.15 ± 0.02</td><td>0.05–0.25</td><td>0.08–0.20</td></tr><tr><td>CH<sub>3</sub>CN</td><td>Methyl cyanide</td><td>0.026 ± 0.003</td><td>&lt;0.12</td><td>&lt;0.38</td><td>&lt;0.25</td><td>0.008–0.054</td><td>0.017–0.036</td></tr><tr><td>HNCO</td><td>Isocyanic acid</td><td>&lt;0.09</td><td>&lt;0.25</td><td></td><td>&lt;0.14</td><td>0.009–0.080</td><td>&lt;0.009–0.08</td></tr><tr><td>CO</td><td>Carbon monoxide</td><td>&lt;2.5</td><td>&lt;21</td><td>&lt;75</td><td></td><td>0.3–35</td><td>0.3–4</td></tr><tr><td>H<sub>2</sub>CO</td><td>Formaldehyde</td><td>0.16 ± 0.02</td><td>&lt;1.3</td><td>&lt;0.6</td><td></td><td>0.16–1.4</td><td>0.16–0.40</td></tr><tr><td>CH<sub>3</sub>OH</td><td>Methanol</td><td>1.7 ± 0.1</td><td>3.5 ± 0.4</td><td>5.0 ± 1.4</td><td>3.9 ± 0.6</td><td>0.6–6.2</td><td>0.8–5.0</td></tr><tr><td>H<sub>2</sub>S</td><td>Hydrogen sulphide</td><td>0.09 ± 0.03</td><td>&lt;0.73</td><td>&lt;2.8</td><td></td><td>0.09–1.5</td><td>0.09–1.0</td></tr><tr><td>CS</td><td>Carbon monosulphide</td><td>0.08 ± 0.01</td><td>&lt;0.12</td><td>&lt;0.35</td><td>0.11 ± 0.04</td><td>0.05–0.54</td><td>0.02–0.20</td></tr><tr><td>SO(SO<sub>2</sub>)</td><td>Sulphur monoxide</td><td>&lt;0.46</td><td>&lt;2.1(&lt; 0.62)</td><td></td><td>&lt;1.4</td><td>0.04–0.30</td><td>&lt;0.09–0.12</td></tr></tbody></table>

according to Ootsubo et al. (2020). Unfortunately the sensitivity of the 2018 observations was limited and our upper limits on the abundance of complex organics such as NH<sub>2</sub>CHO or CH<sub>3</sub>CHO are higher than values measured in comets (Biver et al. 2021), but nevertheless show that comet 21P is not strongly enriched (by a factor >5) in these molecules. The depletion of sulphur species, carbon-chain organics and NH<sub>2</sub> might be present in the volatile phase of comet 21P, but compensated by an enrichment in heavier refractory species as suggested by Ootsubo et al. (2020). Its observations at two different perihelion passages, separated by 20 yr, do not show any significant change.

**Acknowledgements.** This work is based on observations carried out under projects number 035-98, 086-16, 002-18, 112-18 with the IRAM 30-m telescope and M98BH02 with the JCMT. IRAM is supported by the Institut National des Sciences de l'Univers (INSU) of the French Centre national de la recherche scientifique (CNRS), the Max-Planck-Gesellschaft (MPG, Germany) and the Spanish IGN (Instituto Geográfico Nacional). We gratefully acknowledge the support from the IRAM staff for its support during the observations. The Nançay Radio Observatory is operated by the Paris Observatory, associated with the CNRS and with the University of Orléans. The Caltech Submillimeter Observatory was supported by National Science Foundation grant AST 99-80846. The data were reduced and analysed thanks to the use of the GILDAS, class software (http://www.iram.fr/IRAMFR/GILDAS). This research has been supported by the Programme national de planétologie of INSU. Part of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. Part of this work was supported by a fellowship from the James Clerk Maxwell Telescope at the Institute for Astronomy of the University of Hawaii. The James Clerk Maxwell Telescope has historically been operated by the Joint Astronomy Centre on behalf of the Science and Technology Facilities Council of the United Kingdom, the National Research Council of Canada and the Netherlands Organisation for Scientific Research. We are grateful to Henry E. Matthews for his continuous help and contribution to the observations of comets at JCMT. B. P. Bonev acknowledges support of NSF grant AST-2009 398. N.X. Roth acknowledges support by the NASA Postdoctoral Program at the NASA Goddard Space Flight Center, administered by Universities Space Research Association under contract with NASA. Additional support from the NSF grant AST-1614471 for M.A. Cordiner, and by the Planetary Science Division Internal Scientist Funding Program through the Fundamental Laboratory Research (FLaRe) work package (S. N. Milam & M. A. Cordiner). The authors wish to recognise and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain.

## References

* A'Hearn, M. F., Millis, R. L., Schleicher, D. G., Osip, D. J., & Birch, P. V. 1995, *Icarus*, **118**, 223

* Altwegg, K., & Bockelée-Morvan, D. 2003, *Space Sci. Rev.*, **106**, 139

* Altwegg, K., Balsiger, H., Bar-Nun, A., et al. 2015, *Science*, **347**, 1261952

* Biver, N., & Bockelée-Morvan, D. 2019, *ACS Earth Space Chem.*, **3**, 1550

* Biver, N., Bockelée-Morvan, D., Crovisier, J., et al. 2000, *AJ*, **120**, 1554

* Biver, N., Bockelée-Morvan, D., Crovisier, J., et al. 2002, *Earth Moon Planets*, **90**, 323

* Biver, N., Bockelée-Morvan, D., Crovisier, J., et al. 2006, *A&A*, **449**, 1255

* Biver, N., Crovisier, J., Bockelée-Morvan, D., et al., 2012, *A&A*, **539**, A68

* Biver, N., Bockelée-Morvan, D., Debout, V., et al. 2014, *A&A*, **566**, A5

* Biver, N., Bockelée-Morvan, D., Moreno, R., et al. 2015, *Sci. Adv.*, **1**, e1500863

* Biver, N., Moreno, R., Bockelée-Morvan, D., et al. 2016, *A&A*, **589**, A78

* Biver, N., Bockelée-Morvan, D., Paubert, G., et al. 2018, *A&A*, **619**, A127

* Biver, N., Bockelée-Morvan, D., Boissier, J., et al. 2021, *A&A*, **648**, A49

* Bockelée-Morvan, D., & Biver, N. 2017, *Phil. Trans. R. Soc. A*, **375**, 20160252

* Bockelée-Morvan, D., Crovisier, J., Colom, P., & Despois, D. 1994, *A&A*, **287**, 647

* Bockelée-Morvan, D., Hartogh, P., Crovisier, J., et al. 2010, *A&A*, **518**, L149

* Bockelée-Morvan, D., Calmonte, U., Charnley, S., et al. 2015, *Space Sci. Rev.*, **197**, 47

* Bodewits, D., Farnham, T. L., Kelley, M. S. P., & Knight, M. M. 2018, *Nature*, **553**, 186

* Brasser, R., & Morbidelli, A. 2013, *Icarus*, **225**, 40

* Calmonte, U., Altwegg, K., Balsiger, H., et al. 2016, *MNRAS*, **462**, S253

* Carter, M., Lazareff, B., Maier, D., et al. 2012, *A&A*, **538**, A89

* Combi, M. R., Lee, Y., Patel, T. S., et al. 2011, *AJ*, **141**, 128

* Combi, M. R., Mäkinen, T. T., Bertaux, J.-L., Quémerais, E., & Ferron, S. 2019, *Icarus*, **317**, 610

* Cowley, S. W., 2016, *Nature*, **317**, 381

* Crovisier, J., Colom, P., Gérard, E., Bockelée-Morvan, D., Bourgois, G. 2002, *A&A*, **393**, 1053

* DiSanti, M. A., Bonev, B. P., Villanueva, G. L., & Mumma, M. J. 2013, *ApJ*, **763**, 1

* Faggi, S., Mumma, M. J., Villanueva, G. L., Paganini, L., & Lippi, M. 2019, *AJ*, **158**, 254

* Fink, U. 2009, *Icarus*, **201**, 331

* Gérard, E., Bockelée-Morvan, D., Bourgois, G., Colom, P., & Crovisier, J. 1988, *A&AS*, **74**, 485

* Hartogh, P., Lis, D. C., Bockelée-Morvan, D., et al. 2011, *Nature*, **478**, 218

* Huebner, W. F., Keady, J. J., Lyon, S. P. 1992, *Ap&SS*, **195**, 1

* Kelley, M. S. P., Bodewits, D., Quanzhi Ye, et al. 2019, *Res. Notes Amer. Astron. Soc.*, **3**, A126

* Lamy, P. L., Toth, I., Fernandez, Y. R., & Weaver, H. A. 2004, in *Comets II*, eds. M., Festou, H. U., Keller, & H. A., Weaver (University of Arizona Press), 223

* Lis, D. C., Bockelée-Morvan, D., Güsten, R., et al., 2019, *A&A*, **625**, L5

* Moulane, Y., Jehin, E., Opitom, C., et al. 2018, *A&A*, **619**, A158

* Müller, H. S. P., Schlöder, F., Stutzki, J., & Winnewisser, G. 2005, *J. Mol. Struct.*, **742**, 215

* O'Brien, D. P., Walsh, K. J., Morbidelli, A., Raymond, S. N., & Mandell, A. M. 2014, *Icarus*, **239**, 74

* Ootsubo, T., Kawakita, H., Hamada, S., et al. 2012, *ApJ*, **752**, 15

* Ootsubo, T., Kawakita, H., Shinnaka, Y., et al. 2020, *Icarus*, **338**, 113450

* Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, *J. Quant. Spec. Rad. Transf.*, **60**, 883

* Roth, N. X., Gibb, E. L., Bonev, B. P., et al. 2020, *AJ*, **159**, 42

* Saki, M., Gibb, E. L., Bonev, B. P., et al. 2020, *AJ*, **160**, 184

* Schleicher, D. G., Knight, M. M., Eisner, N. L., & Thirouin, A. 2018, *AJ*, **157**, 108

* Weaver, H. A., Chin, G., Bockelée-Morvan, D., et al. 1999, *Icarus*, **142**, 482

* Zakharov, V., Bockelée-Morvan, D., Biver, N., Crovisier, J., & Lecacheux, A. 2007, *A&A*, **473**, 303
## Appendix A: Tables of observed lines

<table><caption>**Table A.1.** Line intensities from JCMT observations.</caption><thead><tr><th>Date<br>(yyyy/mm/dd.dd)</th><th>Molecule</th><th>Transition</th><th>Frequency<sup>(a)</sup><br>(MHz)</th><th>Beam<br>(")</th><th>Offset<sup>(b)</sup><br>(")</th><th>Intensity<br>(K km s<sup>-1</sup>)</th><th>Doppler shift<br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="8"><strong>21P/Giacobini-Zinner</strong></td></tr><tr><td>1998/10/30.26</td><td>HCN</td><td>4-3</td><td>354 505.472</td><td>13.2</td><td>1.5</td><td>0.367 ± 0.014</td><td>-0.09 ± 0.02</td></tr><tr><td>1998/11/01.15</td><td>HCN</td><td>4-3</td><td>354 505.472</td><td></td><td>2.0</td><td>0.366 ± 0.029</td><td>-0.05 ± 0.04</td></tr><tr><td>1998/11/02.29</td><td>HCN</td><td>4-3</td><td>354 505.472</td><td></td><td>2.0</td><td>0.388 ± 0.035</td><td>-0.11 ± 0.05</td></tr><tr><td>1998/11/04.26</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td>18.5</td><td>3.0</td><td>0.210 ± 0.042</td><td>-0.03 ± 0.09</td></tr><tr><td rowspan="5">1998/10/31.26</td><td rowspan="5">CH<sub>3</sub>OH</td><td>7<sub>-1</sub> - 6<sub>-1</sub>E</td><td>338 344.588</td><td rowspan="5">14.0</td><td rowspan="5">2.0</td><td>0.077 ± 0.018</td><td>-0.09 ± 0.13</td></tr><tr><td>7<sub>0</sub> - 6<sub>0</sub>A<sup>+</sup></td><td>338 408.698</td><td>0.072 ± 0.018</td><td>-0.08 ± 0.12</td></tr><tr><td>7<sub>1</sub> - 6<sub>1</sub>E</td><td>338 614.936</td><td>0.028 ± 0.019</td><td></td></tr><tr><td>7<sub>2</sub> - 6<sub>2</sub>A<sup>+</sup></td><td>338 639.802</td><td></td><td>&lt;0.048</td></tr><tr><td>7<sub>2</sub> - 6<sub>2</sub>E</td><td>338 722.290</td><td>0.035 ± 0.025</td><td>+0.01 ± 0.68</td></tr><tr><td>1998/11/04.17</td><td>CH<sub>3</sub>OH</td><td>1<sub>1</sub> - 0<sub>0</sub>A<sup>+</sup></td><td>350 905.100</td><td>13.4</td><td>1.5</td><td>0.119 ± 0.022</td><td>+0.01 ± 0.10</td></tr><tr><td>1998/11/03.23</td><td>H<sub>2</sub>CO</td><td>5<sub>15</sub> - 4<sub>14</sub></td><td>351 768.639</td><td>13.4</td><td>1.5</td><td>0.027 ± 0.010</td><td>-0.29 ± 0.22</td></tr><tr><td>1998/11/01.28</td><td>CS</td><td>7-6</td><td>342 882.850</td><td>13.7</td><td>1.5</td><td>0.095 ± 0.015</td><td>-0.05 ± 0.04</td></tr><tr><td>1998/11/00.8</td><td>CO</td><td>3-2</td><td>345 795.991</td><td>13.6</td><td>2.7</td><td>0.022 ± 0.007</td><td></td></tr><tr><td>1998/11/01.2</td><td>HNC</td><td>4-3</td><td>362 630.303</td><td>13.0</td><td>2.0</td><td>&lt;0.052</td><td></td></tr></tbody></table>

**Notes.** <sup>(a)</sup>From Müller et al. (2005) and Pickett et al. (1998). <sup>(b)</sup>Half power Gaussian beam width and average pointing offset. On average 13.6" corresponded to 9000 km at the distance of the comet.

<table><caption>**Table A.2.** Line intensities from CSO observations.</caption><thead><tr><th>Date<br>(yyyy/mm/dd.dd)</th><th>Molecule</th><th>Transition</th><th>Frequency<sup>(a)</sup><br>(MHz)</th><th>Beam<br>(")</th><th>Offset<sup>(b)</sup><br>(")</th><th>Intensity<br>(K km s<sup>-1</sup>)</th><th>Doppler shift<br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="8"><strong>21P/Giacobini-Zinner</strong></td></tr><tr><td>1998/11/10.25</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td>27.0</td><td>3.0</td><td>0.198 ± 0.015</td><td>-0.16 ± 0.04</td></tr><tr><td>1998/11/13.2</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td></td><td>3.5</td><td>0.116 ± 0.022</td><td>+0.11 ± 0.11</td></tr><tr><td>1998/11/14.18</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td></td><td>3.0</td><td>0.117 ± 0.012</td><td>-0.08 ± 0.06</td></tr><tr><td rowspan="2">1998/11/14.27</td><td rowspan="2">CH<sub>3</sub>OH</td><td>2<sub>1</sub> - 2<sub>0</sub>A<sup>-+</sup></td><td>304 208.324</td><td>23.0</td><td>3.0</td><td>0.054 ± 0.011</td><td>+0.12 ± 0.11</td></tr><tr><td>4<sub>1</sub> - 4<sub>0</sub>A<sup>-+</sup></td><td>307 165.911</td><td>22.8</td><td>3.0</td><td>0.064 ± 0.012</td><td>-0.01 ± 0.09</td></tr></tbody></table>

**Notes.** <sup>(a)</sup>From Müller et al. (2005) and Pickett et al. (1998). <sup>(b)</sup>Half power Gaussian beam width and average pointing offset. On average 27" corresponded to 17 000 km at the distance of the comet.
<table><caption>Table A.3. Line intensities from IRAM observations.</caption><thead><tr><th>Date<br>(yyyy/mm/dd)</th><th>Molecule</th><th>Transition</th><th>Frequency<sup>(a)</sup><br>(MHz)</th><th>Beam<br>(")</th><th>Offset<sup>(b)</sup><br>(")</th><th>Intensity<br>(K km s<sup>-1</sup>)</th><th>Doppler shift<br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="8" style="text-align: left; font-weight: bold;">21P/Giacobini-Zinner</td></tr><tr><td rowspan="2">1998/11/11.71</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td>9.2</td><td>3.0</td><td>0.466 ± 0.047</td><td>-0.13 ± 0.06</td></tr><tr><td>HCN</td><td>1-0</td><td>88 631.848</td><td>26.4</td><td></td><td>0.035 ± 0.011</td><td>-0.10 ± 0.14</td></tr><tr><td rowspan="2">1998/11/12.71</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td>9.2</td><td>3.5</td><td>0.331 ± 0.048</td><td>-0.01 ± 0.08</td></tr><tr><td>HCN</td><td>1-0</td><td>88 631.848</td><td>26.4</td><td></td><td>0.046 ± 0.012</td><td>-0.15 ± 0.17</td></tr><tr><td rowspan="2">1998/11/13.71</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td>9.2</td><td>2.5</td><td>0.567 ± 0.072</td><td>-0.14 ± 0.07</td></tr><tr><td>HCN</td><td>1-0</td><td>88 631.848</td><td>26.4</td><td></td><td>0.051 ± 0.015</td><td>+0.23 ± 0.18</td></tr><tr><td rowspan="2">1998/11/14.68</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td>9.2</td><td>3.5</td><td>0.414 ± 0.082</td><td>-0.07 ± 0.11</td></tr><tr><td>HCN</td><td>1-0</td><td>88 631.848</td><td>26.4</td><td></td><td>0.048 ± 0.009</td><td>-0.05 ± 0.14</td></tr><tr><td rowspan="5">1998/11/15.68</td><td>HCN</td><td>1-0</td><td>88 631.848</td><td>26.4</td><td>3.0</td><td>0.070 ± 0.010</td><td>-0.04 ± 0.09</td></tr><tr><td>H<sub>2</sub>CO</td><td>3<sub>12</sub> - 2<sub>11</sub></td><td>225 697.772</td><td>10.6</td><td>2.0</td><td></td><td>&lt;0.063</td></tr><tr><td>CH<sub>3</sub>OH</td><td>3<sub>0</sub> - 2<sub>0</sub>E</td><td>145 093.707</td><td>16.2</td><td></td><td>0.014 ± 0.014</td><td></td></tr><tr><td></td><td>3<sub>-1</sub> - 2<sub>-1</sub>E</td><td>145 097.370</td><td></td><td></td><td>0.010 ± 0.014</td><td></td></tr><tr><td></td><td>3<sub>0</sub> - 2<sub>0</sub>A</td><td>145 103.152</td><td></td><td></td><td>0.038 ± 0.014</td><td>+0.37 ± 0.27</td></tr><tr><td rowspan="2">1998/11/13.2</td><td></td><td>3<sub>&pm;2</sub> - 2<sub>&pm;2</sub>E</td><td>145 126.290</td><td></td><td></td><td>0.020 ± 0.016</td><td></td></tr><tr><td>H<sub>2</sub>S</td><td>1<sub>10</sub> - 1<sub>01</sub></td><td>168 762.762</td><td>14.0</td><td>3.0</td><td>0.004 ± 0.015</td><td></td></tr><tr><td rowspan="2">2018/09/14.41</td><td>HCN</td><td>1-0</td><td>88 631.848</td><td>27.5</td><td>3.4</td><td>0.230 ± 0.042</td><td>+0.14 ± 0.15</td></tr><tr><td></td><td></td><td></td><td></td><td>7.1</td><td>0.102 ± 0.021</td><td>+0.31 ± 0.19</td></tr><tr><td rowspan="2">2018/09/15.24</td><td>HCN</td><td>1-0</td><td>88 631.848</td><td></td><td>3.3</td><td>0.133 ± 0.011</td><td>+0.07 ± 0.06</td></tr><tr><td></td><td></td><td></td><td></td><td>6.7</td><td>0.119 ± 0.041</td><td>+0.01 ± 0.25</td></tr><tr><td>2018/09/16.44</td><td>HCN</td><td>1-0</td><td>88 631.848</td><td></td><td>3.9</td><td>0.091 ± 0.026</td><td>+0.00 ± 0.20</td></tr><tr><td>2018/09/17.20</td><td>HCN</td><td>1-0</td><td>88 631.848</td><td></td><td>3.3</td><td>0.124 ± 0.011</td><td>-0.06 ± 0.07</td></tr><tr><td>2018/09/17.29</td><td>HCN</td><td></td><td></td><td></td><td>5.4</td><td>0.111 ± 0.015</td><td>-0.11 ± 0.10</td></tr><tr><td>2018/09/18.20</td><td>HCN</td><td>1-0</td><td>88 631.848</td><td></td><td>3.5</td><td>0.128 ± 0.008</td><td>-0.00 ± 0.05</td></tr><tr><td>2018/09/18.29</td><td></td><td></td><td></td><td></td><td>4.4</td><td>0.088 ± 0.019</td><td>+0.11 ± 0.16</td></tr><tr><td>2018/09/15.11</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td>9.4</td><td>5.0</td><td>0.887 ± 0.475</td><td>+0.10 ± 0.31</td></tr><tr><td>2018/09/17.1</td><td>HNC</td><td>1-0</td><td>90 663.568</td><td>26.9</td><td>3.4</td><td>0.007 ± 0.004</td><td></td></tr><tr><td rowspan="5">2018/09/16.2</td><td>CH<sub>3</sub>CN</td><td>8.0-7.0</td><td>147 174.588</td><td>16.8</td><td>2.9</td><td>0.025 ± 0.009</td><td></td></tr><tr><td></td><td>8.1-7.1</td><td>147 171.752</td><td></td><td></td><td>0.027 ± 0.008</td><td></td></tr><tr><td></td><td>8.2-7.2</td><td>147 163.244</td><td></td><td></td><td>0.014 ± 0.008</td><td></td></tr><tr><td></td><td>8.3-7.3</td><td>147 149.068</td><td></td><td></td><td>0.013 ± 0.009</td><td></td></tr><tr><td>CH<sub>3</sub>CN</td><td>9-8 (4 lines)</td><td>165 569-165 540</td><td>15.0</td><td>3.7</td><td>sum = 0.046 ± 0.026</td><td></td></tr><tr><td>2018/09/17.3</td><td>CH<sub>3</sub>CN</td><td>12-11 (4 lines)</td><td>220 709-220 747</td><td>11.1</td><td>4.9</td><td>sum = 0.105 ± 0.056</td><td></td></tr><tr><td>2018/09/17.1</td><td>CH<sub>3</sub>OH</td><td>5<sub>-1</sub> - 4<sub>0</sub>E</td><td>84 521.169</td><td>28.8</td><td>3.4</td><td>0.019 ± 0.004</td><td>-0.20 ± 0.17</td></tr><tr><td rowspan="10">2018/09/16.1</td><td>CH<sub>3</sub>OH</td><td>3<sub>1</sub> - 2<sub>1</sub>A<sup>+</sup></td><td>143 865.801</td><td>17.2</td><td>2.9</td><td>0.035 ± 0.012</td><td>-0.05 ± 0.20</td></tr><tr><td></td><td>3<sub>0</sub> - 2<sub>0</sub>E</td><td>145 093.707</td><td>17.0</td><td></td><td>0.053 ± 0.013</td><td>+0.03 ± 0.17</td></tr><tr><td></td><td>3<sub>-1</sub> - 2<sub>-1</sub>E</td><td>145 097.370</td><td></td><td></td><td>0.066 ± 0.013</td><td>-0.02 ± 0.13</td></tr><tr><td></td><td>3<sub>0</sub> - 2<sub>0</sub>A</td><td>145 103.152</td><td></td><td></td><td>0.059 ± 0.011</td><td>-0.05 ± 0.12</td></tr><tr><td></td><td>3<sub>2</sub> - 2<sub>2</sub>A<sup>-</sup></td><td>145 124.410</td><td></td><td></td><td>0.016 ± 0.011</td><td>-0.28 ± 0.45</td></tr><tr><td></td><td>3<sub>&pm;2</sub> - 2<sub>&pm;2</sub>E</td><td>145 126.290</td><td></td><td></td><td>0.023 ± 0.013</td><td>+0.48 ± 0.46</td></tr><tr><td></td><td>3<sub>1</sub> - 2<sub>1</sub>E</td><td>145 131.855</td><td></td><td></td><td>0.038 ± 0.012</td><td>+0.18 ± 0.19</td></tr><tr><td></td><td>3<sub>2</sub> - 2<sub>2</sub>A<sup>+</sup></td><td>145 133.460</td><td></td><td></td><td>0.017 ± 0.012</td><td></td></tr><tr><td></td><td>3<sub>1</sub> - 2<sub>1</sub>A<sup>-</sup></td><td>146 368.342</td><td>16.9</td><td></td><td>0.021 ± 0.011</td><td>+0.08 ± 0.30</td></tr><tr><td>CH<sub>3</sub>OH</td><td>9<sub>0</sub> - 8<sub>1</sub>A<sup>+</sup></td><td>146 618.794</td><td>16.8</td><td></td><td>0.026 ± 0.011</td><td>-0.38 ± 0.34</td></tr><tr><td rowspan="12">2018/09/17.3</td><td>CH<sub>3</sub>OH</td><td>1<sub>1</sub> - 1<sub>0</sub>E</td><td>165 050.229</td><td>15.0</td><td>3.7</td><td>0.044 ± 0.013</td><td>+0.17 ± 0.19</td></tr><tr><td></td><td>2<sub>1</sub> - 2<sub>0</sub>E</td><td>165 061.187</td><td></td><td></td><td>0.054 ± 0.015</td><td>+0.22 ± 0.24</td></tr><tr><td></td><td>3<sub>1</sub> - 3<sub>0</sub>E</td><td>165 099.300</td><td></td><td></td><td>0.071 ± 0.015</td><td>-0.21 ± 0.15</td></tr><tr><td></td><td>4<sub>1</sub> - 4<sub>0</sub>E</td><td>165 190.539</td><td></td><td></td><td>0.030 ± 0.014</td><td>+0.25 ± 0.22</td></tr><tr><td></td><td>5<sub>1</sub> - 5<sub>0</sub>E</td><td>165 369.410</td><td></td><td></td><td>0.060 ± 0.013</td><td>+0.21 ± 0.17</td></tr><tr><td></td><td>6<sub>1</sub> - 6<sub>0</sub>E</td><td>165 678.724</td><td></td><td></td><td>0.035 ± 0.011</td><td>+0.10 ± 0.22</td></tr><tr><td></td><td>7<sub>1</sub> - 7<sub>0</sub>E</td><td>166 169.179</td><td>14.9</td><td></td><td>0.053 ± 0.012</td><td>+0.02 ± 0.18</td></tr><tr><td></td><td>8<sub>1</sub> - 8<sub>0</sub>E</td><td>166 898.566</td><td>14.8</td><td></td><td>0.021 ± 0.012</td><td>-0.09 ± 0.43</td></tr><tr><td></td><td>9<sub>1</sub> - 9<sub>0</sub>E</td><td>167 931.056</td><td>14.7</td><td></td><td>0.008 ± 0.011</td><td></td></tr><tr><td>CH<sub>3</sub>OH</td><td>3<sub>2</sub> - 2<sub>1</sub>E</td><td>170 060.581</td><td>14.5</td><td></td><td>0.066 ± 0.031</td><td>+0.25 ± 0.36</td></tr><tr><td>CH<sub>3</sub>OH</td><td>4<sub>2</sub> - 3<sub>1</sub>E</td><td>218 440.050</td><td>11.2</td><td>4.9</td><td>0.082 ± 0.027</td><td>+0.17 ± 0.21</td></tr><tr><td>2018/09/18.3</td><td>CH<sub>3</sub>OH</td><td>3<sub>12</sub> - 2<sub>11</sub></td><td>225 697.772</td><td>10.7</td><td>4.6</td><td>0.069 ± 0.028</td><td>+0.27 ± 0.29</td></tr><tr><td>2018/09/18.1</td><td>H<sub>2</sub>CO</td><td>3<sub>03</sub> - 2<sub>02</sub></td><td>218 222.191</td><td>11.2</td><td>4.9</td><td>0.055 ± 0.029</td><td></td></tr><tr><td>2018/09/18.29</td><td>H<sub>2</sub>CO</td><td>2<sub>02</sub> - 1<sub>01</sub></td><td>145 602.953</td><td>17.0</td><td>2.9</td><td>0.000 ± 0.012</td><td></td></tr><tr><td>2018/09/16.2</td><td>H<sub>2</sub>CO</td><td>2<sub>11</sub> - 1<sub>10</sub></td><td>150 498.334</td><td>16.4</td><td>3.3</td><td>0.004 ± 0.011</td><td></td></tr><tr><td>2018/09/16.44</td><td>CO</td><td>2-1</td><td>230 538.000</td><td>10.6</td><td>2.8</td><td></td><td>&lt;0.267</td></tr><tr><td>2018/09/17.3</td><td>H<sub>2</sub>S</td><td>1<sub>10</sub> - 1<sub>01</sub></td><td>168 762.762</td><td>14.6</td><td>3.7</td><td>0.037 ± 0.013</td><td>+0.12 ± 0.21</td></tr><tr><td>2018/09/16.2</td><td>CS</td><td>3-2</td><td>146 969.029</td><td>16.8</td><td>3.1</td><td>0.076 ± 0.008</td><td>+0.02 ± 0.07</td></tr><tr><td>2018/09/14.9</td><td></td><td></td><td></td><td></td><td>6.1</td><td>0.050 ± 0.021</td><td></td></tr><tr><td>2018/09/16.8</td><td>NH<sub>2</sub>CHO</td><td>7<sub>0.7</sub> - &gt; 8<sub>1.8</sub> (5 li.)</td><td>142 701-167 321</td><td>17-15</td><td>3.5</td><td>(average) &lt; 0.017</td><td></td></tr></tbody></table>

Notes. (a) From Müller et al. (2005) and Pickett et al. (1998). (b) Half power Gaussian beam width and average radial pointing offset. On average 9–27" corresponded to 5700–17 000 km at the distance of the comet in 1998 and 2600–7800 km in 2018. (c) Average Doppler shift of the four lines.
<table><caption>**Table A.4.** Line intensities from IRAM observations.</caption><thead><tr><th>Date<br>(yyyy/mm/dd.dd)</th><th>Molecule</th><th>Transition</th><th>Frequency<sup>(a)</sup><br>(MHz)</th><th>Beam<br>(")</th><th>Offset<sup>(b)</sup><br>(")</th><th>Intensity<br>(K km s<sup>-1</sup>)</th><th>Doppler shift<br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="8"><strong>41P/Tuttle-Giacobini-Kresák</strong></td></tr><tr><td>2017/04/26.02</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td>9.3</td><td>1.2</td><td>0.378 ± 0.038</td><td>-0.11 ± 0.06</td></tr><tr><td>2017/05/01.02</td><td>HCN</td><td>3-2</td><td>265 886.434</td><td></td><td>5.2</td><td>0.438 ± 0.100</td><td>-0.25 ± 0.17</td></tr><tr><td>2017/04/27.4</td><td>HNC</td><td>3-2</td><td>271 981.142</td><td>9.1</td><td>2.2</td><td>&lt;0.080</td><td></td></tr><tr><td>2017/05/01.08</td><td>CH<sub>3</sub>CN</td><td>8-7 (4 lines)</td><td>147 174-147 149</td><td>16.5</td><td>2.0</td><td>(sum) &lt; 0.111</td><td></td></tr><tr><td>2017/05/01.08</td><td>CH<sub>3</sub>CN</td><td>9-8 (4 lines)</td><td>165 569-165 540</td><td>14.5</td><td>2.0</td><td>(sum) &lt; 0.141</td><td></td></tr><tr><td rowspan="12">2017/04/26.08</td><td rowspan="12">CH<sub>3</sub>OH</td><td>5<sub>0</sub> - 4<sub>0</sub>E</td><td>241 700.168</td><td rowspan="12">10.1</td><td rowspan="12">1.2</td><td>0.108 ± 0.033</td><td>-0.10 ± 0.17</td></tr><tr><td>5<sub>-1</sub> - 4<sub>-1</sub>E</td><td>241 767.247</td><td>0.110 ± 0.034</td><td>-0.69 ± 0.31</td></tr><tr><td>5<sub>0</sub> - 4<sub>0</sub> A</td><td>241 791.367</td><td>0.082 ± 0.033</td><td>+0.05 ± 0.23</td></tr><tr><td>5<sub>1</sub> - 4<sub>1</sub>E</td><td>241 879.038</td><td>0.032 ± 0.034</td><td></td></tr><tr><td>5<sub>±2</sub> - 4<sub>±2</sub>E</td><td>241 904.401</td><td>0.125 ± 0.038</td><td>-0.12 ± 0.25</td></tr><tr><td>5<sub>2</sub> - 4<sub>2</sub> A<sup>-</sup></td><td>241 842.287</td><td>0.018 ± 0.034</td><td>-0.40 ± 0.63</td></tr><tr><td>5<sub>2</sub> - 4<sub>2</sub> A<sup>+</sup></td><td>241 887.678</td><td>0.034 ± 0.034</td><td></td></tr><tr><td>5<sub>3</sub> - 4<sub>3</sub> A</td><td>241 832.910</td><td>0.065 ± 0.032</td><td>-0.46 ± 0.36</td></tr><tr><td>5<sub>1</sub> - 4<sub>1</sub> A<sup>-</sup></td><td>243 915.811</td><td>0.070 ± 0.038</td><td>+0.26 ± 0.31</td></tr><tr><td>3<sub>3</sub> - 3<sub>2</sub> A<sup>+-</sup></td><td>251 917.057</td><td rowspan="4">9.7</td><td rowspan="4">2.2</td><td>0.055 ± 0.022</td><td>-0.15 ± 0.30</td></tr><tr><td>3<sub>3</sub> - 3<sub>2</sub> A<sup>-</sup></td><td>251 905.720</td><td>0.035 ± 0.021</td><td>-0.44 ± 0.48</td></tr><tr><td>4<sub>3</sub> - 4<sub>2</sub> A<sup>+-</sup></td><td>251 900.439</td><td>0.035 ± 0.021</td><td>-0.55 ± 0.53</td></tr><tr><td>4<sub>3</sub> - 4<sub>2</sub> A<sup>-</sup></td><td>251 866.510</td><td>0.025 ± 0.022</td><td>-0.35 ± 0.70</td></tr><tr><td>5<sub>3</sub> - 5<sub>2</sub> A<sup>+-</sup></td><td>251 890.868</td><td>0.038 ± 0.022</td><td>-0.08 ± 0.42</td></tr><tr><td>5<sub>3</sub> - 5<sub>2</sub> A<sup>-</sup></td><td>251 811.936</td><td>0.028 ± 0.020</td><td>-0.24 ± 0.52</td></tr><tr><td>6<sub>3</sub> - 6<sub>2</sub> A<sup>+-</sup></td><td>251 895.703</td><td>0.043 ± 0.020</td><td>-0.47 ± 0.35</td></tr><tr><td>6<sub>3</sub> - 6<sub>2</sub> A<sup>-</sup></td><td>251 738.411</td><td>0.045 ± 0.018</td><td>+0.01 ± 0.25</td></tr><tr><td>7<sub>3</sub> - 7<sub>2</sub> A<sup>+-</sup></td><td>251 923.671</td><td>0.004 ± 0.021</td><td></td></tr><tr><td>7<sub>3</sub> - 7<sub>2</sub> A<sup>-</sup></td><td>251 641.754</td><td>-0.001 ± 0.018</td><td></td></tr><tr><td>8<sub>3</sub> - 8<sub>2</sub> A<sup>+-</sup></td><td>251 984.802</td><td>0.009 ± 0.020</td><td></td></tr><tr><td>8<sub>3</sub> - 8<sub>2</sub> A<sup>-</sup></td><td>251 517.269</td><td>0.028 ± 0.020</td><td>-0.31 ± 0.52</td></tr><tr><td>9<sub>3</sub> - 9<sub>2</sub> A<sup>+-</sup></td><td>252 090.369</td><td>0.010 ± 0.021</td><td></td></tr><tr><td>9<sub>3</sub> - 9<sub>2</sub> A<sup>-</sup></td><td>251 359.841</td><td>0.020 ± 0.019</td><td></td></tr><tr><td rowspan="12">2017/05/01.08</td><td rowspan="12">CH<sub>3</sub>OH</td><td>2<sub>0</sub> - 1<sub>-1</sub>E</td><td>254 015.367</td><td>9.6</td><td>2.2</td><td>0.038 ± 0.018</td><td>-0.06 ± 0.21</td></tr><tr><td>5<sub>2</sub> - 4<sub>-1</sub>E</td><td>266 838.123</td><td>9.3</td><td>2.2</td><td>0.090 ± 0.022</td><td>-0.20 ± 0.15</td></tr><tr><td>1<sub>1</sub> - 1<sub>0</sub>E</td><td>165 050.229</td><td rowspan="5">14.6</td><td rowspan="5">2.0</td><td>0.025 ± 0.022</td><td>-0.16 ± 0.25</td></tr><tr><td>2<sub>1</sub> - 2<sub>0</sub>E</td><td>165 061.187</td><td>0.027 ± 0.023</td><td>+0.48 ± 0.38</td></tr><tr><td>3<sub>1</sub> - 3<sub>0</sub>E</td><td>165 099.300</td><td>0.014 ± 0.023</td><td>-0.39 ± 0.49</td></tr><tr><td>4<sub>1</sub> - 4<sub>0</sub>E</td><td>165 190.539</td><td>0.037 ± 0.023</td><td>+0.24 ± 0.18</td></tr><tr><td>5<sub>1</sub> - 5<sub>0</sub>E</td><td>165 369.410</td><td>-0.038 ± 0.024</td><td></td></tr><tr><td>6<sub>1</sub> - 6<sub>0</sub>E</td><td>165 678.724</td><td>0.002 ± 0.024</td><td></td></tr><tr><td>7<sub>1</sub> - 7<sub>0</sub>E</td><td>166 169.179</td><td>0.003 ± 0.025</td><td></td></tr><tr><td>2017/04/26.08</td><td>H<sub>2</sub>CO</td><td>3<sub>1,2</sub> - 2<sub>1,1</sub></td><td>225 697.772</td><td>10.8</td><td>1.2</td><td>&lt;0.090</td><td></td></tr><tr><td>2017/05/01.08</td><td>H<sub>2</sub>CO</td><td>2<sub>1,1</sub> - 1<sub>1,0</sub></td><td>150 498.334</td><td>16.1</td><td>2.0</td><td>&lt;0.059</td><td></td></tr><tr><td>2017/04/26.08</td><td>CO</td><td>2 - 1</td><td>230 538.000</td><td>10.6</td><td>1.2</td><td>&lt;0.085</td><td></td></tr><tr><td>2017/04/26.08</td><td>HNCO</td><td>11<sub>0,11</sub> - 10<sub>0,10</sub></td><td>241 774.032</td><td>10.1</td><td>1.2</td><td>0.062 ± 0.031</td><td></td></tr><tr><td>2017/05/01.08</td><td>H<sub>2</sub>S</td><td>1<sub>1,0</sub> - 1<sub>0,1</sub></td><td>168 762.762</td><td>14.3</td><td>2.0</td><td>&lt;0.085</td><td></td></tr><tr><td>2017/04/26.08</td><td>CS</td><td>5-4</td><td>244 935.557</td><td>10.0</td><td>1.2</td><td>&lt;0.110</td><td></td></tr><tr><td>2017/05/01.08</td><td>CS</td><td>3-2</td><td>146 969.029</td><td>16.5</td><td>2.0</td><td>&lt;0.054</td><td></td></tr><tr><td>2017/05/27.4</td><td>SO</td><td>6.5-5.4</td><td>251 825.770</td><td>9.7</td><td>2.2</td><td>&lt;0.064</td><td></td></tr></tbody></table>

Notes. (a)From Müller et al. (2005) and Pickett et al. (1998). (b)Half power Gaussian beam width and average radial pointing offset. On average 9–16.5" corresponded to 1200–2200 km at the distance of the comet.
<table><caption>**Table A.5.** Line intensities from IRAM observations.</caption><thead><tr><th>Date<br>(yyyy/mm/dd.dd)</th><th>Molecule</th><th>Transition</th><th>Frequency<sup>(a)</sup><br>(MHz)</th><th>Beam<br>(")</th><th>Offset<sup>(b)</sup><br>(")</th><th>Intensity<br>(K km s<sup>-1</sup>)</th><th>Doppler shift<br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="8"><strong>64P/Swift-Gehrels</strong></td></tr><tr><td>2018/12/11.75</td><td>HCN</td><td>3–2</td><td>265 886.434</td><td>9.3</td><td>2.9</td><td>0.349 ± 0.035</td><td>-0.225 ± 0.063</td></tr><tr><td>2018/12/17.95</td><td>CH<sub>3</sub>CN</td><td>14–11 (4 lines)</td><td>257 527–257 448</td><td>9.5</td><td>2.4</td><td>(sum) &lt; 0.147</td><td></td></tr><tr><td rowspan="3">2018/12/11.75</td><td>CH<sub>3</sub>OH</td><td>J<sub>3</sub> – J<sub>2</sub>A<sup>±</sup> (8 1.)</td><td>251 738–251 924</td><td>9.7</td><td rowspan="3">2.9</td><td rowspan="3">sum = 0.158 ± 0.048</td><td rowspan="3">-0.12 ± 0.14</td></tr><tr><td>CH<sub>3</sub>OH</td><td>11<sub>0</sub> – 10<sub>1</sub> A<sup>+</sup></td><td>250 507.016</td><td></td></tr><tr><td>CH<sub>3</sub>OH</td><td>2<sub>0</sub> – 1<sub>-1</sub> E</td><td>254 015.377</td><td>9.6</td></tr><tr><td rowspan="7">2018/12/17.95</td><td>CH<sub>3</sub>OH</td><td>5<sub>2</sub> – 4<sub>1</sub> E</td><td>266 838.123</td><td>9.3</td><td rowspan="7">2.4</td><td>0.008 ± 0.020</td><td></td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>0</sub> – 4<sub>0</sub>E</td><td>241 700.168</td><td rowspan="7">10.0</td><td>0.012 ± 0.019</td><td>-0.14 ± 0.13</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>-1</sub> – 4<sub>-1</sub>E</td><td>241 767.247</td><td>0.074 ± 0.022</td><td>+0.04 ± 0.21</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>0</sub> – 4<sub>0</sub> A</td><td>241 791.367</td><td>0.062 ± 0.020</td><td>-0.27 ± 0.17</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>1</sub> – 4<sub>1</sub>E</td><td>241 879.038</td><td>0.124 ± 0.021</td><td>-0.17 ± 0.09</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>±2</sub> – 4<sub>±2</sub>E</td><td>241 904.401</td><td>0.048 ± 0.021</td><td>+0.03 ± 0.23</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>1</sub> – 4<sub>1</sub> A<sup>-</sup></td><td>243 915.811</td><td>0.045 ± 0.021</td><td>-0.22 ± 0.28</td></tr><tr><td>CH<sub>3</sub>OH</td><td>2<sub>1</sub> – 1<sub>0</sub>E</td><td>261 805.736</td><td></td><td>0.034 ± 0.024</td><td>-0.44 ± 0.50</td></tr><tr><td>2018/12/17.95</td><td>CS</td><td>5–4</td><td>244 935.557</td><td>9.9</td><td>2.4</td><td>0.062 ± 0.022</td><td>-0.12 ± 0.19</td></tr><tr><td>2018/12/17.95</td><td>SO</td><td>6.6–5.5</td><td>258 255.826</td><td>9.5</td><td>2.4</td><td>0.025 ± 0.022</td><td></td></tr><tr><td rowspan="2">2018/12/17.95</td><td rowspan="2">HNCO</td><td>6.7–5.6</td><td>261 843.721</td><td>9.4</td><td rowspan="2">2.4</td><td>-0.058 ± 0.025</td><td></td></tr><tr><td>11<sub>0,11</sub> – 10<sub>0,10</sub></td><td>241 774.032</td><td>10.0</td><td>0.003 ± 0.021</td><td></td></tr><tr><td colspan="8"><strong>38P/Stephan-Oterma</strong></td></tr><tr><td>2018/12/13.12</td><td>CH<sub>3</sub>CN</td><td>8–7 (4 lines)</td><td>147 175–147 149</td><td>16.7</td><td>2.2</td><td>(sum) &lt; 0.162</td><td></td></tr><tr><td>2018/12/13.12</td><td>CH<sub>3</sub>CN</td><td>9–8 (4 lines)</td><td>165 569–165 540</td><td>14.9</td><td>2.2</td><td>(sum) &lt; 0.195</td><td></td></tr><tr><td>2018/12/17.93</td><td>CH<sub>3</sub>CN</td><td>14–11 (4 lines)</td><td>257 527–257 448</td><td>9.5</td><td>2.2</td><td>(sum) &lt; 0.386</td><td></td></tr><tr><td>2018/12/13.12</td><td>CS</td><td>3–2</td><td>146 969.029</td><td>16.7</td><td>2.2</td><td>&lt;0.075</td><td></td></tr><tr><td>2018/12/17.93</td><td>CS</td><td>5–4</td><td>244 935.557</td><td>9.9</td><td>3.1</td><td>&lt;0.136</td><td></td></tr><tr><td>2018/12/12/12</td><td>CH<sub>3</sub>OH</td><td>8<sub>-1</sub> – 7<sub>0</sub>E</td><td>229 758.811</td><td>10.6</td><td>2.2</td><td>&lt;0.147</td><td></td></tr><tr><td>2018/12/13.12</td><td>CH<sub>3</sub>OH</td><td>J<sub>1</sub> – J<sub>0</sub>E (6 lines)</td><td>165 050–166 678</td><td>14.9</td><td>2.2</td><td>(sum) &lt; 0.242</td><td></td></tr><tr><td rowspan="6">2018/12/17.93</td><td>CH<sub>3</sub>OH</td><td>3<sub>2</sub> – 2<sub>1</sub>E</td><td>170 060.581</td><td>14.6</td><td rowspan="6">3.1</td><td>&lt;0.140</td><td></td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>0</sub> – 4<sub>0</sub>E</td><td>241 700.168</td><td rowspan="6">10.0</td><td>0.012 ± 0.044</td><td></td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>-1</sub> – 4<sub>-1</sub>E</td><td>241 767.247</td><td>0.099 ± 0.043</td><td>+0.14 ± 0.22</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>0</sub> – 4<sub>0</sub> A</td><td>241 791.367</td><td>0.089 ± 0.039</td><td>+0.18 ± 0.20</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>1</sub> – 4<sub>1</sub>E</td><td>241 879.038</td><td>0.028 ± 0.045</td><td></td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>±2</sub> – 4<sub>±2</sub>E</td><td>241 904.401</td><td>0.051 ± 0.044</td><td>-0.23 ± 0.47</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5<sub>1</sub> – 4<sub>1</sub> A<sup>-</sup></td><td>243 915.811</td><td></td><td>0.038 ± 0.047</td><td></td></tr><tr><td rowspan="2">2018/12/12.12</td><td>H<sub>2</sub>CO</td><td>3<sub>1,3</sub> – 2<sub>1,2</sub></td><td>211 211.469</td><td>11.4</td><td rowspan="2">2.2</td><td>&lt;0.109</td><td></td></tr><tr><td>H<sub>2</sub>CO</td><td>3<sub>1,2</sub> – 2<sub>1,1</sub></td><td>225 697.772</td><td>10.7</td><td>&lt;0.132</td><td></td></tr><tr><td>2018/12/13.12</td><td>H<sub>2</sub>CO</td><td>2<sub>1,1</sub> – 1<sub>1,0</sub></td><td>150 498.334</td><td>16.3</td><td>2.2</td><td>&lt;0.128</td><td></td></tr><tr><td>2018/12/12.12</td><td>CO</td><td>2–1</td><td>230 538.000</td><td>10.5</td><td>2.2</td><td>&lt;0.117</td><td></td></tr><tr><td>2018/12/13.12</td><td>H<sub>2</sub>S</td><td>1<sub>1,0</sub> – 1<sub>0,1</sub></td><td>168 762.762</td><td>14.7</td><td>2.2</td><td>&lt;0.125</td><td></td></tr></tbody></table>

Notes. (a) From Müller et al. (2005) and Pickett et al. (1998). (b) Half power Gaussian beam width and average radial pointing offset. On average 9.3–10.0" corresponded to 4000–4400 km at the distance of comet 64P and 9.4–16.7" corresponded to 5200–9300 km at the distance of 38P.