# Chemical composition of comets C/2021 A1 (Leonard) and C/2022 E3 (ZTF) from radio spectroscopy and the abundance of HCOOH and HNCO in comets

N. Biver, D. Bockelée-Morvan, B. Handzlik, Aa. Sandqvist, J. Boissier, M. N.

Drozdovskaya, R. Moreno, J. Crovisier, D. C. Lis, M. Cordiner, et al.

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

► To cite this version:

N. Biver, D. Bockelée-Morvan, B. Handzlik, Aa. Sandqvist, J. Boissier, et al. Chemical composition of comets C/2021 A1 (Leonard) and C/2022 E3 (ZTF) from radio spectroscopy and the abundance of HCOOH and HNCO in comets. Astronomy & Astrophysics - A&A, 2024, 690, pp.A271. 10.1051/0004-6361/202450921. [hal-04733412](https://hal.science/hal-04733412v1)

HAL Id: hal-04733412

https://hal.science/hal-04733412v1

Submitted on 12 Oct 2024

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.

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

Distributed under a Creative Commons Attribution 4.0 International License
# Chemical composition of comets C/2021 A1 (Leonard) and C/2022 E3 (ZTF) from radio spectroscopy and the abundance of HCOOH and HNCO in comets*,**

N. Biver<fnref n="1" />,*,* D. Bockelée-Morvan<fnref n="1" />, B. Handzlik<fnref n="2" />, D., Aa. Sandqvist<fnref n="3" />, J. Boissier<fnref n="4" />, M. N. Drozdovskaya<fnref n="5" />, R. Moreno<fnref n="1" />, J. Crovisier<fnref n="1" />, D. C. Lis<fnref n="6" />, M. Cordiner<fnref n="7" />,<fnref n="8" />, S. Milam<fnref n="7" />, D., N. X. Roth<fnref n="7" />,<fnref n="9" />, B. P. Bonev<fnref n="9" />, N. Dello Russo<fnref n="10" />, R. Vervack<fnref n="10" />, D., C. Opitom<fnref n="11" />, and H. Kawakita<fnref n="12" />

Received 30 May 2024 / Accepted 26 July 2024

## ABSTRACT

We present the results of a molecular survey of long period comets C/2021 A1 (Leonard) and C/2022 E3 (ZTF). Comet C/2021 A1 was observed with the Institut de radioastronomie millimétrique (IRAM) 30-m radio telescope in November-December 2021 before perihelion (heliocentric distance 1.22 to 0.76 au) when it was closest to the Earth (≈0.24 au). We observed C/2022 E3 in January-February 2023 with the *Odin* 1-m space telescope and IRAM 30-m, shortly after its perihelion at 1.11 au from the Sun, and when it was closest to the Earth (≈0.30 au). Snapshots were obtained during 12-16 November 2021 period for comet C/2021 A1. Spectral surveys were undertaken over the 8-13 December 2021 period for comet C/2021 A1 (8 GHz bandwidth at 3 mm, 16 GHz at 2 mm, and 61 GHz in the 1 mm window) and over the 3-7 February 2023 period for comet C/2022 E3 (25 GHz at 2 mm and 61 GHz at 1 mm). We report detections of 14 molecular species (HCN, HNC, CH₃CN, HNCO, NH₂CHO, CH₃OH, H₂CO, HCOOH, CH₃CHO, H₂S, CS, OCS, C₂H₅OH and aGg⁺-(CH₂OH)₂) in both comets. In addition, HC₃N, and CH₂OHCHO were marginally detected in C/2021 A1, and CO and H₂O (with *Odin*) were detected in C/2022 E3. The spatial distribution of several species (HCN, HNC, CS, H₂CO, HNCO, HCOOH, NH₂CHO, and CH₃CHO) is investigated. Significant upper limits on the abundances of other molecules and isotopic ratios are also presented. The activity of comet C/2021 A1 did not vary significantly between 13 November and 13 December 2021, when observations stopped, just before it started to exhibit major outbursts seen in the visible and from observations of the OH radical. Short-term variability in the outgassing of comet C/2022 E3 of the order of ±20% is present and possibly linked to its 8h rotation period. Both comets exhibit rather low abundances relative to water for volatile species such as CO (<2%) and H₂S (0.15%). Methanol is also rather depleted in comet C/2021 A1 (0.9%). Following their revised photo-destruction rates, HNCO and HCOOH abundances in comets observed at millimetre wavelengths have been reevaluated. Both molecules are relatively enriched in these two comets (~0.2% relative to water). Since the combined abundance of these two acids (0.1-1%) is close to that of ammonia in comets, we cannot exclude that these species could be produced by the dissociation of ammonium formate and ammonium cyanate if present in comets.

**Key words.** molecular data – comets: general – radio lines: planetary systems – submillimeter: planetary systems – comets: individual: C/2021 A1 (Leonard) – comets: individual: C/2022 E3 (ZTF)

## 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 are thus our best window on the volatile composition of the solar proto-planetary disk. Comets may also have 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 the early Solar System's evolution

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

<footnote n="2">² Astronomical Observatory, Jagiellonian University, ul. Orla 171, 30-244, Kraków, Poland</footnote>

<footnote n="3">³ Stockholm Observatory, Stockholm University, SCFAB-AlbaNova, 10691 Stockholm, Sweden</footnote>

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

<footnote n="5">⁵ Physikalish-Meteorologisches Observatorium Davos und Weltstrahlungszentrum (PMOD/WRC) Dorfstrasse 33, 7260 Davos Dorf, Switzerland</footnote>

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

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

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

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

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

<footnote n="11">¹¹ Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK</footnote>

<footnote n="12">¹² Koyama Astronomical Observatory, Kyoto Sangyo University, Motoyama, Kamitamo, Kita, Kyoto 603-8555, Japan</footnote>

<footnote>* Based on observations carried out with the IRAM-30m telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain).</footnote>

<footnote>** *Odin* is a Swedish-led satellite project funded jointly by the Swedish National Space Board (SNSB), the Canadian Space Agency (CSA), the National Technology Agency of Finland (Tekes), the Centre National d'Études Spatiales (CNES, France), and the European Space Agency (ESA). The Swedish Space Corporation, today OHB Sweden, is the prime contractor, also responsible for *Odin* operations.</footnote>

<footnote>*** Corresponding author; nicolas.biver@obspm.fr</footnote>
(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 early in the outer asteroid belt and participated in the supply of water to the Earth. On the other hand, current comets are coming from either the Oort Cloud or the scattered disk of the Kuiper belt and may have formed in the same trans-Neptunian region, sampling the same diversity of formation conditions. Understanding the diversity in composition and isotopic ratios of cometary materials is thus essential for assessing such scenarios (Altwegg & Bockelée-Morvan 2003; Bockelée-Morvan et al. 2015).

Comet C/2021 A1 (Leonard) is a long-period dynamically old Oort-Cloud comet (OCC, with an initial semi-major axis of 2000 au and an inclination of 132.7°) that reached perihelion at a heliocentric distance $r_h=0.615$ au on 3.3 January 2022 UT. It came as close as 0.234 au to the Earth on 12.6 December 2021. It was anticipated to become a bright comet visible to the naked eye, but it under-performed until 14 December when it was in solar conjunction. Then monitoring of the outgassing rate via observations of the OH radical (Crovisier et al. 2021) and visual magnitudes showed strong outbursts of activity, repeating on a ~5 day period from 15 December to 7 January 2022¹. Then, it developed spectacular ion and dust tails while being mostly observable from the southern hemisphere. Later on, the comet faded rapidly, and images taken in April-May 2022 suggest that it was disintegrating during this outbursting phase. The derived pre-disintegration nucleus radius was 0.6 ± 0.2 km (Jewitt et al. 2023). We observed comet C/2021 A1 with the Institut de radioastronomie millimétrique (IRAM) 30-m telescope briefly between 12.0 and 16.4 November, and extensively between 8.4 and 13.4 December 2021 UT, before it reached declinations that are too low for the northern observatories.

Comet C/2022 E3 (ZTF) is also a long-period dynamically old OCC (with an initial semi-major axis of 1400 au and an inclination of 109.2°). It reached perihelion at $r_h = 1.112$ au on 12.8 January 2023 UT. It was closest to the Earth on 1.8 February 2023 at 0.284 au, and reached naked eye visibility (total magnitude of 4.8). It attracted public attention and was the target of a worldwide campaign because it was discovered 11 months ahead of its peak brightness, expected to happen at the end of January 2023 when the comet was circumpolar for the northern hemisphere and visible to the naked eye all night<fnref n="2" />. The comet was further advertised in NASA news releases at the end of 2021³. We observed comet C/2022 E3 with the *Odin* space telescope from 19.3 to 20.3 January 2023 (Biver et al. 2023b) and with the IRAM-30m radio telescope between 3.7 and 7.1 February. Due to its high activity and a favourable apparition, this comet was the focus of an international observing campaign, from the radio – OH observed with the Nançay Radio Telescope (NRT) and the Green Bank Telescope – to the infrared (ground-based observations with IRTF and Keck-NIRSPEC and from space with the James Webb Space Telescope on 28 February and 4 March, Milam et al. 2023).

In this paper, we report clear detections of HCN, HNC, CH$_3$CN, HNCO, NH$_2$CHO, CH$_3$OH, H$_2$CO, HCOOH, H$_2$S, and CS in both comets, the more marginal detections of

CO, HC$_3$N, CH$_3$CHO, OCS, (CH$_2$OH)$_2$, CH$_2$OHCHO, and C$_2$H$_5$OH, obtained by averaging several lines, as well as significant upper limits on the abundances of SO, SO$_2$, H$_2$CS, CH$_2$CO, PH$_3$, and other species.

In addition, following the revised photo-destruction rates of several molecules by Hrodmarsson & van Dishoeck (2023), especially of HNCO and HCOOH for which the lifetime has been reduced by nearly one order of magnitude in comparison to previously published or assumed destruction rates (Heays et al. 2017; Huebner & Mukherjee 2015; Biver et al. 2021a), we have reevaluated their abundances in all comets in which they were observed or searched for.

In Sect. 2, we present the observations and spectra of the detected molecules. The information extracted from the observations to analyse the data and compute production rates is provided in Sect. 3. In Sect. 4, we present the retrieved production rates and abundances or upper limits, which are discussed and compared to other comets. The new analysis of HNCO and HCOOH observations are detailed in Sect. 5 followed by the conclusions in Sect. 6.

## 2. Observations

### 2.1. Observations of comet C/2021 A1 with the IRAM-30m

Comet C/2021 A1 (Leonard)<fnref n="4" /> was the focus of a worldwide campaign as it was expected to become very bright at its closest approach to Earth (total visual magnitude $m_1 \sim 4$). It was less active than anticipated until 14 December, when it started to undergo a series of recurrent outbursts bringing it to a maximum brightness of around $m_1 = 3$ between 15 and 25 December 2021.

It was the target of the observing proposal 100-21 scheduled at the IRAM 30-m telescope between 8 and 13 December 2021. Weather conditions were marginal during the first four days (half of the day too icy, windy, or foggy to observe, with precipitable water vapour (*pwv*) in the 3 to 7 mm range otherwise). On the first day (observations were taking place during daytime), when observations resumed after an ice storm, the left-over ice or the impact of de-icing on the antenna reduced the beam efficiency by a factor of 3 (at 1 mm wavelength) to 1.5 (at 2 mm). The last two days offered good weather conditions with *pwv* = 1–3 mm (Table A.1). A few other snapshots (less than 0.5 h of observations) were obtained between 12 and 16 November during observing run 001-21 when comet 67P/Churyumov-Gerasimenko was too high in elevation for tracking (Biver et al. 2023a).

Observations were obtained in wobbler switching mode with the secondary (wobbling) mirror alternating pointing between the ON and OFF positions separated by 180° every 2 seconds. The wobbler was not working on 14–16 November 2021 and we had to use to the position switching mode (PSW). The reference OFF positions in PSW mode were at 300° from the source, alternating ON and OFF every 15 seconds. We used the EMIR (Eight MIR Receiver, Carter et al. 2012) 3 mm, 2 mm, and 1 mm band receivers in 2SB mode connected to the FTS (Fast Fourier Transform Spectrometer) and the VESPA (VERsatile SPectrometer Array) high-resolution spectrometer (Table A.1). The FTS offers an instantaneous bandwidth of 16 GHz in two polarisations with 200 kHz spectral resolution. The VESPA autocorrelator was optimised to provide 4–6 windows of 20–40 kHz spectral resolution on lines of interest in the centre of the Intermediate Frequency windows (6.25 ± 0.25 GHz).

<footnote n="1">¹ http://www.aerith.net/comet/catalog/2021A1/2021A1.html</footnote>

<footnote n="2">² https://doi.org/10.5281/zenodo.12940168</footnote>

<footnote n="3">³ https://science.nasa.gov/resource/whats-up-december-2021/, https://www.nasa.gov/image-article/views-of-comet-leonard-from-two-sun-watching-spacecraft/</footnote>

<footnote n="4">⁴ https://doi.org/10.5281/zenodo.12940168</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Coarse map of the HCN(3–2) line integrated intensity in comet C/2021 A1 on 13.4 December 2021. The beam size is 9.3". Solar phase angle was 158.6°. The direction of the Sun (PA=193.1°) is indicated at the lower right.</figcaption></figure>

Comet C/2021 A1 was tracked with the IRAM-30m control software (NCS) using the latest JPL Horizons<fnref n="5" /> orbital elements available:#15 in November, #17 on 8 December, and #18 on 9–13 December 2021. Offsets (up to 7") were added during the observation to take into account the difference between the position computed by the NCS and the very latest (1–3 days old) astrometric measurements. This was a critical point as for such a new comet relatively close to Earth, positions errors can easily reach the beam size (10" at 240 GHz). The comet position was also checked in real time using coarse (5–9 points) mapping of HCN *J*=3–2 (Fig. 1). Ephemeris offsets were computed afterwards using the JPL#19 orbit solution (yielding values in agreement with the astrometry that was used). Final pointing offsets for each observation were computed after taking into account reconstruction of pointing errors and finally from intensity maps of HCN, that generally yielded residual offsets of less than 1.2".

Representative spectra of comet C/2021 A1 (Leonard) are provided in Figs. 2 and 3 and Figs. A.1–A.3. They show individual lines or averages of several lines for the complex organics from the observations centred on the nucleus. Detailed line intensities are provided in Tables online and for some molecules with useful spatial information in Table B.5 of Sect. 4.3. Spectra showing the full spectral coverage with the wide band FTS are available online.

## 2.2. Observations of comet C/2022 E3 (ZTF) with Odin

Comet C/2022 E3 (ZTF) was favourably placed (solar elongation between 60°and 120°) for *Odin* during one of its yearly astronomy science operations. 15 orbits (~24h) were dedicated to observe the H₂O(1₁₀ – 1₀₁) line at 556.9 GHz in this comet. *Odin* (Frisk et al. 2003) is a small satellite in a polar orbit (period 95 min) equipped with a 1.1 m sub-millimetre radiometer. *Odin* houses four sub-mm receivers covering the 486–504 GHz and 541–581 GHz range with 3 backends: one acousto-optical-spectrometer (AOS) with 1 GHz band width and 1.2 MHz spectral resolution and two autocorrelators (AC1 and AC2), covering 120 MHz with 140 kHz resolution in their highest resolution mode. Two of the receivers, 555B2 and 549A1, were supposed to observe the 556.9 GHz line, but the 549A1

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Millimetre lines observed with IRAM-30m in comet C/2021 A1 during 12–16 November 2021 (average intensity). The CH₃OH line is the sum of the *J*=1 to 4 *J*₁ – *J*₀ E lines at 165050.175, 165061.130, 165099.240, and 165190.475 MHz observed with the high resolution backend VESPA. The vertical axis is main beam brightness temperature in K and the horizontal axis is Doppler velocity in the rest frame of the comet with respect to the main line.</figcaption></figure>

receiver did not lock, and we only obtained data from the 555B2 receiver connected to the AOS and AC2. We aimed at alternating between pointing on the comet position and coarse 3 × 3 points maps spaced by 60 or 120". The beam size was 127". *Odin* achieved successful tracking of the comet on most of 14 of the 15 orbits scheduled from 19.3 to 20.3 January 2023 UTC (Table A.2). Part of the data from the 13th and 14th orbits was lost also due to memory overload. Reconstruction of the attitude of *Odin* did not always converge with good accuracy due to the lack of reference stars in the star-trackers field of view, but most of the pointings were within 15" (≈12% of the beam width) of the comet position. Figure 4 shows the result from the combination of all the maps.

## 2.3. Observations of comet C/2022 E3 (ZTF) with the IRAM-30m

Comet C/2022 E3 (ZTF) was the target of the observing proposal 097-22 scheduled at the IRAM 30-m telescope between 3 and 8 February 2023. Weather conditions were very good to good during the first four days (but too windy followed by a technical issue with the antenna temperature control during 2/3 of the first night). The last night (February 7/8) was lost due to bad weather (snow + wind). Pointing and focus stability was not very good at the beginning of the nights. Very low opacity on 5.0 February (0.2–1 mm *pwv*) enabled a short observation at 177–185 GHz, but the weather conditions degraded and the observing time was too limited to get a useful result on the H₂O line at 183.3 GHz. Table A.2 provides a log of the observations.

Wobbler switching mode with the secondary (wobbling) mirror alternating pointing between the ON and OFF positions separated by 120" every 2 seconds could only be used on the last night as it had again issues as in November 2021. We had to rely on position switching mode (PSW) with reference OFF positions at 240", alternating ON and OFF every 15 seconds. Even using symmetric mode (ON OFF OFF ON sequence), alternating the reference position in +RA and -RA, leaves some residual at the position of the strong ozone atmospheric lines, where the spectra are noisier with poorer baselines, especially around 231.13,

<footnote n="5">⁵ https://ssd.jpl.nasa.gov/horizons.cgi</footnote>
<figure><img src="image_5.png" /></figure>

**Fig. 3.** Individual lines observed with IRAM-30m in comet C/2021 A1 during 8–13 December 2021 (weighted averages). The vertical axis is main beam brightness temperature in K and the horizontal axis is Doppler velocity in the rest frame of the comet with respect to the line. The position and relative intensity of the hyperfine components of the HCN(3–2) line have been drawn below the line.

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

**Fig. 4.** Map of the average spectra of the H₂O(1₁₀–1₀₁) line obtained with *Odin* between 19.45 and 20.16 January 2023 in comet C/2022 E3 (ZTF), as a function of the pointing offset (in arcsec). A baseline in red is plotted on each spectrum. Solar phase angle was 61.6°, and the direction of the Sun (PA=120.2°) is indicated at the upper right with scales for individual AC2 spectra.

237.17, 242.34, 249.80, 249.97, and 267.28 GHz. As for comet C/2021 A1, we used the EMIR receivers and the FTS and VESPA spectrometers (Sect. 2.1).

Comet C/2022 E3 was tracked with the IRAM-30m control software (NCS) using the latest JPL Horizons orbital elements available at the beginning of the observations, JPL#41, adding up offsets (~5" in declination) from the latest astrometric measurements available and confirmed by coarse mapping of the strongest lines in the setup. Final offset in the reduced data take into account a more recent (JPL#45) orbit solution and reconstructed pointing corrections of the antenna to cope with potential approximations during real time observations (Fig. 5).

Representative spectra of comet C/2022 E3 (ZTF) are provided in Fig. 6 and Figs. A.4–A.6. They show individual lines or averages of several lines for the complex organics from the observations centred on the nucleus. Detailed line intensities are provided online, and for some molecules with useful spatial information in Table B.6 of Section 4.3. Spectra showing the full spectral coverage with the FTS are provided online.

Figs. A.3 and A.6 show the average of 3 to 97 lines, weighted according to the noise in each spectrum. We have selected the strongest lines with expected similar S/N (within a factor around four) which are not blended with other lines. This is only for the purpose of highlighting detection and line shapes.

## 2.4. Observations with the Nançay Radio Telescope

The NRT performances and the observing procedure used for comets were described in Crovisier et al. (2002). The integration time is usually about one hour per day. However, from 12
<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Coarse (3×3 points) map of the HCN(3–2) line integrated intensity on 4.67 February 2023 in comet C/2022 E3. Solar phase angle was 45.5°. The direction of the Sun (PA=287°) is indicated at the upper right.</figcaption></figure>

September to 14 December 2021, it was limited to 30 min due to work on the focal track. The beam size at 18-cm wavelength is 3.5 × 18'.

### 2.4.1. Comet C/2021 A1 (Leonard)

C/2021 A1 (Leonard) was observed at the NRT from 1 October 2021 to 14 February 2022 almost every day or every other day<fnref n="6" />.

A succession of outbursts was observed, beginning on 13 December just at the end of the IRAM observations (Crovisier et al. 2021; Jehin et al. 2021). The OH production rate, which was $2.5 \pm 0.1 \times 10^{28}$ molecules s⁻¹ on average for 8–12 December, rose to $6.3 \pm 0.3 \times 10^{28}$ molecules s⁻¹ on 13.5, up to $22.1 \pm 0.2 \times 10^{28}$ molecules s⁻¹ on 15.5 December.

Representative spectra for selected dates are shown in Fig. 7. The evolution of the retrieved OH production rate is plotted in Fig. 8.

### 2.4.2. Comet C/2022 E3 (ZTF)

C/2022 E3 (ZTF) was observed at the NRT from 17 October to 23 December 2022⁷. The retrieved OH production rate rose from $3.3 \times 10^{28}$ molecules s⁻¹. Then the observations were interrupted due to a technical failure on 24 December 2022. They were resumed from 17 to 24 February 2023 in a degraded mode, resulting in an upper limit $Q(\text{OH}) \le 6 \times 10^{28}$ molecules s⁻¹.

## 3. Data analysis

### 3.1. Expansion velocity and outgassing pattern

The lines profiles (Figs. 2–3) of comet C/2021 A1 do not show systematic asymmetry, excepted for small shifts with respect to the rest velocity in the comet frame. In November the mean Doppler shift of the line is slightly negative ($\delta v = -0.03 \pm 0.03$ km s⁻¹ for HCN(3–2)) and positive for most lines in December ($\delta v(\text{HCN}(3–2))= +0.18 \pm 0.01$ km s⁻¹). The two-Gaussian fit (Biver et al. 2021a) used to estimate the expansion

velocity from the velocity at half maximum intensity $VHM$, yields $VHM = -0.75 \pm 0.04$ km s⁻¹ and $VHM = +0.53 \pm 0.03$ km s⁻¹ in November (from HCN) and $VHM = -0.63 \pm 0.01$ km s⁻¹ and $VHM = +0.97 \pm 0.01$ km s⁻¹ in December (all lines). Those asymmetries are expected for a preferential outgassing at a higher rate and velocity on the sunward hemisphere, since in November the phase angle was ~51° (the sunlit hemisphere is mostly facing us) and in December it was between 116° and 159° (Table A.1) so that we were mostly facing the night side of the comet. We could have simulated a two-component outgassing pattern with a higher production rate and expansion velocity on the day side, but this would not change significantly the retrieved production rates for the optically thin lines from using an isotropic model with a constant expansion velocity equal to the mean of the day and night sides. Thus, we assumed isotropic outgassing with $v_{exp} = 0.60$ km s⁻¹ in November and $v_{exp} = 0.73$ km s⁻¹ in December. The actual expansion velocities needed to fit the observed profiles are 5–10% lower than the $VHM$ (that is –0.58 and 0.88 km s⁻¹ in December) due to thermal broadening. An example of simulated profile with asymmetric outgassing is shown in Fig. 9.

In the case of comet C/2022 E3 (ZTF), the full width at half maximum (FWHM) of the water line at 556.9 GHz is of the order of $1.60$ km s⁻¹, but the line is asymmetric and redshifted due to opacity effects (see for example Lecacheux et al. 2003). In addition the uncertainty on the frequency calibration is of the order of $0.05$ km s⁻¹, so it is difficult to derive information on the gas velocity asymmetry from these data, but the average expansion velocity should be of the order of $0.8$ km s⁻¹. IRAM data obtained two weeks later provide more accurate line profiles but also show evidence of variation with time. On the basis of the lines with the best signal-to-noise ratio (HCN, CH₃OH, H₂S, CS, H₂CO, and CH₃CN), the weighted average $VHMs$ are $-0.83 \pm 0.01$ and $+0.64 \pm 0.01$ km s⁻¹, with a mean Doppler shift of the lines of $-0.12 \pm 0.01$ km s⁻¹. Since the phase angle was not too large (45°, Table A.2), we deduce that outgassing was larger on the day-side with a larger expansion velocity. From the $VHMs$ and fits to line profiles (Figs. 10, 11) we estimate that the expansion velocity was of the order of $0.76$ km s⁻¹ on the day side and $0.52$ km s⁻¹ on the night side, with a production rate two times higher on the day side to get the average measured Doppler shift. Nevertheless, assuming isotropic outgassing with the average $v_{exp} = 0.68$ km s⁻¹ yields very similar production rates and we assume isotropic outgassing at this velocity to compute and compare the production rates. Model fits with constant velocities in Figs. 9, 10, and 11 are underestimating the signal around zero velocity. Better fits to the line shapes are obtained when simulating radial acceleration of the expansion velocity, using the parameter $x_{acc} = 4$ from the formula in Biver et al. (2011, 2022) and $v_{exp,0} = 0.9 \times v_{exp}$. This illustrate a possible way to improve the fit, with a slightly lower (≈3–10%) production rate for most molecules. But there are other possible explanations, like variation of the gas temperature and different azimuthal gas distribution. We keep the simpler model with constant velocity and isotropic outgassing for the computation of production rates and abundances of the numerous molecules observed.

### 3.2. Gas temperature

Several species such as methanol are detected through multiple transitions coming from different rotational upper energy levels $E_u$. In the case of local thermal equilibrium (LTE) we

<footnote n="6">⁶ https://lesia.obspm.fr/planeto/cometes/basecom/LD/indexld.html</footnote>

<footnote n="7">⁷ https://lesia.obspm.fr/planeto/cometes/basecom/ZT/indexzt.html</footnote>
<figure><img src="image_8.png" /></figure>

**Fig. 6.** Individual lines observed with IRAM-30m in comet C/2022 E3 during 3–6 February 2023 (average intensity, excluding offset positions). The vertical axis is main beam brightness temperature in K. The horizontal axis is the Doppler velocity in the rest frame of the comet with respect to the line. A galactic contamination of CO has been blanked out. The position and relative intensity of the hyperfine components of the HCN(3–2), HCN(2–1), and H$^{13}$CN(3–2) lines have been drawn below the lines.

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

**Fig. 7.** Selected averages of the OH lines observed at the NRT in comet C/2021 A1 (Leonard) in November–December 2021 (averages of the 1667 and 1665 MHz lines scaled to 1667 MHz).

expect the relative population levels ($p_u$) to follow the Boltzmann law ($p_u \propto \exp(-E_u/kT)$). When we plot the logarithm of the populations $p_u$ versus the upper energy levels $E_u$ (the so-called rotational diagram) the mean slope of the linear fit of the data is $1/T_{rot}$, where $T_{rot}$ is the rotational temperature, equal to the gas kinetic temperature $T$ in LTE. Due to radiative decay and infrared pumping, deviations from LTE can be observed in some series of lines, especially the CH$_3$OH lines at 242 GHz, for which $T_{rot} < T$. In such a case the full non-LTE modelling of the evolution of the population of the rotational level throughout the cometary atmosphere is required to estimate the value of $T$ resulting in the measured $T_{rot}$ within the radio telescope beam. We provide in Figs. C.1–C.13 the rotational diagrams for species for which we observed several transitions with sufficient signal-to-noise ratio (S/N>3) and in Table 1 the measured $T_{rot}$ and inferred $T$.

The average gas kinetic temperature $T$ for comet C/2021 A1 inferred from CH$_3$OH and CH$_3$CN rotational temperatures (Table 1) is $T = 25 \pm 5$ K in November and $T = 61 \pm 3$ K in December. In December the temperature inferred from HNCO lines (Figs. C.4) suggests a lower value of $T$, but we have not taken into account the effect of infrared pumping that could modify the value expected for $T_{rot}$ (HNCO) for a given $T$. The spatial distribution of HNCO is also not well constrained: if it comes
<figure><img src="image_10.png" /></figure>

**Fig. 8.** Time evolution of comet C/2021 A1 (Leonard), single days, and averages of selected periods in November–December 2021. The OH production rates from the NRT are plotted in black. The daily water production rates inferred from SOHO/SWAN observations of the wide H Lyman-α coma, which tends to smooth out short term variations, are plotted in red. (Combi et al. 2023a).

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

**Fig. 9.** Average spectrum of the HCN(3–2) line observed from 8.5 to 13.4 December in comet C/2021 A1, with simulated profiles. The model in blue assumes a production rate of 2.1 × 10<sup>25</sup> molecules s<sup>-1</sup> at 0.88 km s<sup>-1</sup> on the sunward hemisphere and 1.1 × 10<sup>25</sup> molecules s<sup>-1</sup> at 0.58 km s<sup>-1</sup> on the opposite hemisphere mostly facing the observer. The mean phase angle is ~140°, and the tilt of 40°(or 140°) with respect to the comet-observer line of sight is taken into account in this 3D simulation, as depicted in the upper right. The model in purple uses variable velocities in both hemispheres (see text), with production rates of 2.3 × 10<sup>25</sup> molecules s<sup>-1</sup> and 0.9 × 10<sup>25</sup> molecules s<sup>-1</sup>, respectively. The vertical axis is main beam brightness temperature in K. The horizontal axis is the Doppler velocity in the rest frame of the comet.

from a distributed source then the observed rotational temper-
ature will be compatible with a higher value of T. We have
adopted T = 30 K in November and T = 60 K in December.

In the case of comet C/2022 E3 (ZTF), due to higher signal-
to-noise ratios, we have numerous measurements of rotational
temperatures (rotational diagrams in Figs. C.5–C.13). Inferred
gas temperatures *T* are provided in Table 1. From methanol data,
the weighted average and dispersion is *T* = 58 ± 4 K, from
CH<sub>3</sub>CN we get *T* = 68 ± 9 K, and *T* = 53 ± 16 K from the
four other species in Table 1. We adopt the average *T* = 60 K

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

**Fig. 10.** Average FTS spectrum of the HCN(3–2) line observed from 3.7 to 6.7 February in comet C/2022 E3 with simulated profiles. The model in blue assumes a production rate of 3.0 × 10<sup>25</sup> molecules s<sup>-1</sup> at 0.76 km s<sup>-1</sup> on the sunward hemisphere and 1.5 × 10<sup>25</sup> molecules s<sup>-1</sup> at 0.52 km s<sup>-1</sup> on the other hemisphere, as depicted in the upper right. The model in purple uses variable velocities in both hemispheres (see text), with production rates of 3.0 × 10<sup>25</sup> molecules s<sup>-1</sup> and 1.0 × 10<sup>25</sup> molecules s<sup>-1</sup>, respectively. The mean phase angle is ~45°. The vertical axis is main beam brightness temperature in K, and the horizontal axis is Doppler velocity in the rest frame of the comet.

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

**Fig. 11.** Average spectrum of the CH<sub>3</sub>OH line at 266.838 GHz observed from 3.7 to 6.7 February 2023 in comet C/2022 E3 at the same time as HCN(3–2). The model in blue assumes a production rate of 5.6 × 10<sup>26</sup> molecules s<sup>-1</sup> at 0.76 km s<sup>-1</sup> on the sunward hemisphere and 2.8 × 10<sup>26</sup> molecules s<sup>-1</sup> at 0.52 km s<sup>-1</sup> on the other hemisphere, as depicted in the upper right, while the red profile comes from symmetric outgassing at 0.68 km s<sup>-1</sup> and same total production rate (8.4 × 10<sup>26</sup> molecules s<sup>-1</sup>). The model in purple uses the same variable velocities in both hemispheres as in Fig. 10, with production rates of 6.4 × 10<sup>25</sup> molecules s<sup>-1</sup> and 2.0 × 10<sup>25</sup> molecules s<sup>-1</sup>, in sunward and other hemispheres, respectively. The mean phase angle is ~45°. The vertical axis is main beam brightness temperature in K, and the horizontal axis is Doppler velocity in the rest frame of the comet.

to derive production rates of comet C/2022 E3 for all species.
All measurements, given their uncertainties, are within 10 K of
this value (Table 1) but due to likely different collision rates
for the different molecules, small deviations are not surprising.
Also for NH<sub>2</sub>CHO, HCOOH, HNCO, and CH<sub>3</sub>CHO we have
not taken into account infrared pumping that could introduce
more differences between T<sub>rot</sub> and T. We have also investigated
<table><caption>**Table 1.** Rotational temperatures and inferred gas kinetic temperatures.</caption><thead><tr><th>UT<br>(mm/dd.d)</th><th>Molecule</th><th>Freq. range<br>(GHz)</th><th>lines<br>(c)</th><th>off.<sup>(a)</sup><br>(")</th><th><i>T</i><sub>rot</sub><sup>(b)</sup><br>(K)</th><th><i>T</i><sub>gas</sub><br>(K)</th></tr></thead><tbody><tr><td colspan="7"><strong>C/2021 A1 (Leonard)</strong></td></tr><tr><td>11/15.9</td><td>CH<sub>3</sub>OH</td><td>250–254</td><td>12</td><td>1.4</td><td>23 ± 5</td><td>26 ± 6</td></tr><tr><td>11/14.8</td><td>CH<sub>3</sub>OH</td><td>165–169</td><td>6</td><td>1.5</td><td>22 ± 9</td><td>21 ± 9</td></tr><tr><td>12/12.8</td><td>CH<sub>3</sub>OH</td><td>241.8</td><td>13</td><td>1.6</td><td>54.1 ± 7.1</td><td>63 ± 8</td></tr><tr><td>12/10.6</td><td>CH<sub>3</sub>OH</td><td>250–254</td><td>20</td><td>1.7</td><td>65.3 ± 5.2</td><td>70 ± 6</td></tr><tr><td>12/11.1</td><td>CH<sub>3</sub>OH</td><td>165–169</td><td>10</td><td>1.9</td><td>56.0 ± 6.9</td><td>56 ± 7</td></tr><tr><td>12/10.6</td><td>CH<sub>3</sub>OH</td><td>250–267</td><td>5</td><td>1.7</td><td>54.4 ± 7.3</td><td>55 ± 7</td></tr><tr><td>12/12.5</td><td>HNCO</td><td>153–264</td><td>9</td><td>1.5</td><td>30.4 ± 3.6</td><td>37 ± 6</td></tr><tr><td>12/12.5</td><td>HNCO<sub>d</sub><sup>(d)</sup></td><td></td><td></td><td></td><td>37.0 ± 5.4</td><td>48 ± 7</td></tr><tr><td>12/11.1</td><td>CH<sub>3</sub>CN</td><td>147–165</td><td>11</td><td>1.9</td><td>53 ± 11</td><td>53 ± 10</td></tr><tr><td>12/12.8</td><td>CH<sub>3</sub>CN</td><td>220–257</td><td>13</td><td>1.7</td><td>62.9 ± 9.3</td><td>63 ± 9</td></tr><tr><td>12/12.5</td><td>NH<sub>2</sub>CHO</td><td>213–267</td><td>25</td><td>1.5</td><td>116 ± 59</td><td>102 ± 52</td></tr><tr><td colspan="7"><strong>C/2022 E3 (ZTF)</strong></td></tr><tr><td>02/05.2</td><td>CH<sub>3</sub>OH</td><td>165–169</td><td>10</td><td>2.0</td><td>53.0 ± 1.7</td><td>55 ± 2</td></tr><tr><td></td><td></td><td></td><td></td><td>9.1</td><td>51.5 ± 6.3</td><td>57 ± 7</td></tr><tr><td></td><td></td><td></td><td></td><td>13.0</td><td>53.5 ± 7.8</td><td>63 ± 10</td></tr><tr><td>02/05.2</td><td>CH<sub>3</sub>OH</td><td>165–169</td><td>v-<sup>(e)</sup></td><td>2.0</td><td>52.3 ± 1.9</td><td>54 ± 2</td></tr><tr><td></td><td></td><td></td><td>v+</td><td>2.0</td><td>55.0 ± 3.5</td><td>58 ± 4</td></tr><tr><td>02/04.6</td><td>CH<sub>3</sub>OH</td><td>250–254</td><td>28</td><td>1.6</td><td>58.4 ± 1.6</td><td>63 ± 2</td></tr><tr><td></td><td></td><td></td><td>22</td><td>9.8</td><td>48.7 ± 4.6</td><td>59 ± 6</td></tr><tr><td>02/04.6</td><td>CH<sub>3</sub>OH</td><td>250–267</td><td>5</td><td>1.6</td><td>54.3 ± 1.9</td><td>56 ± 2</td></tr><tr><td>02/05.1</td><td>CH<sub>3</sub>OH</td><td>241.8</td><td>14</td><td>1.8</td><td>42.2 ± 1.4</td><td>54 ± 3</td></tr><tr><td>02/05.1</td><td>CH<sub>3</sub>CN</td><td>257</td><td>4</td><td>1.8</td><td>68 ± 16</td><td>68 ± 16</td></tr><tr><td>02/05.2</td><td>CH<sub>3</sub>CN</td><td>147–165</td><td>12</td><td>2.1</td><td>52.3 ± 4.0</td><td>73 ± 8</td></tr><tr><td>02/05.8</td><td>CH<sub>3</sub>CN</td><td>220–239</td><td>9</td><td>2.7</td><td>47.4 ± 9.2</td><td>54 ± 13</td></tr><tr><td>02/05.2</td><td>NH<sub>2</sub>CHO</td><td>149–265</td><td>40</td><td>2.1</td><td>79 ± 14</td><td>84 ± 15</td></tr><tr><td>02/05.2</td><td>HNCO</td><td>153–264</td><td>12</td><td>2.1</td><td>50 ± 9</td><td>76 ± 13</td></tr><tr><td>02/05.2</td><td>HCOOH</td><td>151–270</td><td>29</td><td>2.1</td><td>49 ± 9</td><td>49 ± 9</td></tr><tr><td>02/05.2</td><td>CH<sub>3</sub>CHO</td><td>151–271</td><td>79</td><td>2.1</td><td>43 ± 6</td><td>44 ± 6</td></tr></tbody></table>

**Notes.** Subscript “d” has been added to the molecules for which a daughter Haser density profile is assumed with the parent scale length provided below. (a)Mean pointing offset. (b)Result of non-linear fit with $\chi^2$ minimisation. (c)Number of lines used for the determination of $T_{rot}$. In some cases lines have been averaged by groups of two lines having similar spectral characteristics (same J level, close energy levels and Einstein A coefficients). (d)Assuming a distributed source with a parent scale length of $L_p=5000$ km. (e)Temperatures deduced from the negative (blueshifted) part of the line (respectively positive or redshifted side for v+) sampling mostly the day-side of the coma (respectively night-side).

possible day versus night differences of the coma temperature. The positive parts of the 166 GHz methanol lines only suggest marginally higher temperature on the night side than on the day-side (blueshifted parts of the lines, Table 1).

## 4. Production rates and abundances

Production rates are computed using our excitation and radiative transfer codes, and parameters as in previous papers (Biver et al. 2021a). We assumed isotropic outgassing and a constant velocity and temperature (Sects. 3.1 and 3.2). We provide both daily production rates when the molecules are detected with a sufficiently high signal-to-noise (Tables B.1 and B.2) and averages over the 13–16 November and 8–13 December 2021 periods for comet C/2021 A1 (Leonard) and 3–7 February 2023 for comet C/2022 E3 (ZTF) (Tables B.3 and B.4). For the molecules for

which we observed several lines ‘i’, either detected individually with a S/N >5 or not, the final production rate is the weighted average of all production rates $Q_i$. They are computed on each considered line (even when $Q_i$ is negative) and averaged with weighting according to $1/\sigma(Q_i)^2$ were $\sigma(Q_i)$ is the uncertainty in production rate deduced from the line ‘i’.

### 4.1. Reference water production rate

Water production rates are inferred from the monitoring of OH lines at 18-cm with the NRT (Sect. 2.4) and H$_2$O observations at 557 GHz with *Odin* (Sect. 2.2). Water productions rates for these comets were also obtained from SOHO/SWAN (Combi et al. 2023a,b). For the period of observations with the IRAM-30m radio telescope, and to compute relative abundances, we estimate $Q_{H_2O} = 2, 3$, and $4 \times 10^{28}$ molecules s$^{-1}$ for the 12, 13–16 November, and 8–13 December 2021, respectively, for comet C/2021 A1 (Leonard). For comet C/2022 E3 (ZTF), extrapolation from *Odin* production rates in Table 2, following $1/r_h^2$, yields $Q_{H_2O} = 5 \times 10^{28}$ molecules s$^{-1}$ for the 3–7 February period. These values also follow the longer term trend of brightness evolution. We use these values to determine collisional rates and abundances relative to water.

### 4.2. Temporal variations in the production rates

Evidence for short-term variability in the outgassing is present both in *Odin* and IRAM data for comet C/2022 E3 (ZTF), but not for comet C/2021 A1. Visual inspection of the production rates folded on a single period seems to indicate a periodic pattern of about 0.35 days in $Q_{H_2O}$ and 0.38 days in $Q_{CH_3OH}$, while optical observations have reported a periodic pattern in the CN structures of $0.363 \pm 0.004$ days (Knight et al. 2023) and $0.354 \pm 0.004$ days (Manzini et al. 2023). We fitted (Figs. 12 and 13) a sine-wave variation (for simplicity – a more complex profile would require more parameters) to our water and methanol production rates and found similar periods. Table 2 provides the fitted parameters and their uncertainty based on $\chi^2$ minimisation for a simple sinusoidal pattern. The relative amplitudes are 18% for methanol and 5% for H$_2$O. H$_2$O being observed with a beam about ten times larger (the *Odin* beam radius of 33000 km covers molecules emitted over a full period of 0.38 days), and having optically thick lines, it is expected to display shallower variations. The amplitude of the variations for production rate of methanol are more likely to represent the full extent of variations of outgassing of the nucleus and are similar to the variations ($\pm 20\%$) observed for HCN. This has to be taken into account when deriving precise ratios (e.g. isotopic ratios in Sect. 4.6), but for general abundances averaged over the four days the impact is small, especially as on the 4, 5, and 6 February evenings (Table A.2) the observations covered more than a full period of ~0.36 days.

For other species day-to-day variations of the production rates are provided in Tables B.1 and B.2, and plotted in Figs. 14 and 15, together with some estimates of the $A\rho$ quantity<fnref n="8" /> in order to have a rough idea of the relative evolution of the dust production rate.

### 4.3. Molecules coming from a distributed source

The coarse mapping done during the observations provides some information on the spatial extent of the emission for some

<footnote n="8">⁸ https://www.lesia.obspm.fr/comets</footnote>
<table><caption>**Table 2.** Time periods found to fit comet C/2022 E3 (ZTF) production rates.</caption><thead><tr><th>Molecule</th><th>Q<sup>(a)</sup><br>(molecules s<sup>-1</sup>)</th><th>ΔQ<sup>(b)</sup><br>(molecules s<sup>-1</sup>)</th><th>T<sub>p</sub><br>(days)</th><th>Confidence<br>Level or χ<sup>2</sup></th><th>Method</th></tr></thead><tbody><tr><td>H<sub>2</sub>O</td><td>6.27 ± 0.07 × 10<sup>28</sup></td><td>0</td><td>0</td><td>χ<sub>10</sub><sup>2</sup> = 2.51</td><td>Constant</td></tr><tr><td>H<sub>2</sub>O</td><td>6.23 ± 0.08 × 10<sup>28</sup></td><td>0.31 ± 0.11 × 10<sup>28</sup></td><td>0.349 ± 0.022</td><td>χ<sub>10</sub><sup>2</sup> = 2.44</td><td>Sine first order</td></tr><tr><td>CH<sub>3</sub>OH</td><td>9.12 ± 0.18 × 10<sup>26</sup></td><td>1.69 ± 0.28 × 10<sup>26</sup></td><td>0.3841 ± 0.0029</td><td>χ<sub>8</sub><sup>2</sup> = 2.31</td><td>Sine first order</td></tr></tbody></table>

Notes. <sup>(a)</sup>Mean production rate. <sup>(b)</sup>Amplitude in production rate of the sine fitting.

<figure><img src="image_14.png" /><figcaption>**Fig. 12.** Water production rates (in molecules s⁻¹) of comet C/2022 E3 (ZTF) between 19.3 and 20.3 January 2023 derived from observations with *Odin* (central pointings). A best fit sinusoidal variation is plotted in blue.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 13.** Methanol production rates (in molecules s⁻¹) of comet C/2022 E3 (ZTF) between 3.7 and 7.1 February 2023 derived from observations with IRAM 30m. An average 5% calibration uncertainty has been added to the formal uncertainty due to system temperature noise. A best fit sinusoidal variation is plotted in blue.</figcaption></figure>

molecules, within 14" from the nucleus (2400 km at the distance of comet C/2021 A1 in December 2021) to ~20" (4500 km) for comet C/2022 E3 (ZTF). Tables B.5 and B.6 provide information on the spatial extent of the emission of some molecules (average line intensities as a function of pointing offset) and constraints on the scale length and production from a distributed source. We adjusted a Haser daughter species density profile to the data and provide the result of the χ² minimisation with 1σ uncertainty on the scale length ($L_p$) and production rate ($Q_p$).

<figure><img src="image_16.png" /><figcaption>**Fig. 14.** Production rates of comet C/2021 A1 (Leonard) in November and December 2021. Water production rates are 1.1 × $Q_{OH}$ from Sect. 2.4. $Af\rho$ were measured from images of the comet by N. Biver, not corrected for the phase angle.</figcaption></figure>

Values for other fixed scale lengths ($Q_p$ and corresponding $\chi^2$) are provided. These results have to be taken with caution since: (i) the excitation of the rotational levels is not perfectly known due to variable temperature in the coma, unknown collisional cross-sections, and other poorly estimated processes that can mimic extended emission, (ii) the comet was variable in activity (it underwent a huge outburst on 14 December just after the last observation), (iii) the observations are mostly probing the 1600–5000 km spatial scales (beam size – extent of the mapping), so that information on very different scale lengths is not well constrained.

**HCN:** it has been shown to be mostly released from nucleus ices (Cordiner et al. 2014), but the IR observations (Dello Russo et al. 2016; Lippi et al. 2021) that probe the inner part of the coma also generally find higher production rates of HCN than the radio ones (Biver et al. 2024). The observations of both comets (data were treated day by day to limit time variability effects) suggest
<figure><img src="image_17.png" /><figcaption>**Fig. 16.** Spatial distribution of formaldehyde in comet C/2022 E3 (ZTF). This plot shows for the average of all observations of the twin 211 GHz and 226 GHz lines observed together in comet C/2022 E3 (ZTF) from 3 to 7 February 2023: (i) in black squares (with error-bars) the observed line integrated intensity, (ii) in solid and dotted lines the expected intensity for parent scale length of 0 or 1000–8000 km. For better visibility a normalisation has been applied: at each offset the line intensity is multiplied by the column density at the first offset divided by the one expected at the considered offset point, for a parent distribution. Production rate of the models are adjusted to get the same integrated intensity at the first point. Deviation from a flat line for $L_p = 0$ (black continuous line) is due to excitation effects. Vertical axis is normalised integrated intensity in mK km/s and horizontal axis is pointing offset converted into km at the distance of the comet.</figcaption></figure>

for other lines ($2_{11}$ – $1_{10}$ line at 150.498 GHz) provide looser constraints that fully encompass those values. Fig. 16 shows the combination of all H₂CO observations at 211 and 226 GHz in comet C/2022 E3 versus expected evolution of line intensity as a function of pointing offset for various parent scale lengths – the retrieved $L_p(\text{H}_2\text{CO})$ is 1700 ± 500 km which corresponds to $0.36 \pm 0.11 \times L_d(\text{H}_2\text{CO})$ at 1.18 au from the sun, with $v_{exp} = 0.68$ km s⁻¹.

**NH₂CHO** and **CH₃CHO**: we have also studied the spatial distribution of these two species (only NH₂CHO in C/2021 A1): results are not conclusive and may be affected by excitation effects (especially for the high $K_a$ transitions of formamide). Slightly distributed production (with $L_p \sim 1500$ km) could be favoured for comet C/2022 E3, but a nuclear source is not clearly excluded (within the $2-\sigma$ uncertainty) as for formaldehyde.

In addition, if a significant production of the molecules comes from the sublimation of icy grains in the coma, all species could also show some distributed production. But the least distributed species like HCN, could put an upper limit on the spatial extent of the distributed production. The production of CS and H₂CO remains clearly more distributed and they are thus daughter products.

### 4.4. Relative abundances

Abundances relative to water, assuming $Q_{\text{H}_2\text{O}} = 2 \times 10^{28}$ molecules s⁻¹ in November and $Q_{\text{H}_2\text{O}} = 4 \times 10^{28}$ molecules s⁻¹ in December 2021 for comet C/2021 A1 and $Q_{\text{H}_2\text{O}} = 5 \times 10^{28}$ molecules s⁻¹ in early February 2023 for comet C/2022 E3 are provided in Table 3. Abundances have also been measured for several molecules in comet C/2021 A1 in the infrared by Faggi et al. (2023), but during the outbursting period of the comet in December–January. While these authors

found very similar abundances for HCN, OCS, H₂CO, and compatible with our upper limit on CO, they did not detect CH₃OH. Their inferred CH₃OH/H₂O is very low, down to one order of magnitude lower than our abundance (based on the secure detection of over 42 lines). Either the comet had exhausted its methanol content during the outburst phase at the end of December, or the infrared observations missed part of the methanol emission, following a trend also seen in millimetre investigations (Biver et al. 2011) which show an apparent decrease in methanol abundance as comets get closer to the Sun. The methanol observations of Faggi et al. (2023) were obtained at 0.62 au from the Sun, where the sampled gas coma gets warmer (assumed to be 120 K) and the fraction of methanol in an excited torsional state may also depopulate the ground vibrational levels.

### 4.5. Upper limits

The surveys cover part of the 2 mm wavelength range and most of the 1 mm (210–272 GHz) as shown in spectra online. Most of the known organic molecules have many transitions in these wavelength ranges. We have not noticed clearly unidentified lines (above the $5-\sigma$ level). We checked the sensitivity of the observations for the following molecules: CH₃SH, CH₃OCH₃, CH₃COCH₃, c-C₂H₄O, c-C₃H₂, CH₃NH₂, C₂H₃CN, CH₃Cl, and propanal. Averaging the lines expected to be the strongest did not yield any significant signal. $3-\sigma$ upper limits on abundances and production rates are provided in Tables 3, B.3, and B.4, and for the species not listed the upper limits on abundances are not expected to be better than in previous comets such as 46P/Wirtanen (Biver et al. 2021a).

### 4.6. Isotopic ratios

Isotopic ratios of H, C, N, and S and upper limits on their abundances in some molecules are provided in Table 4. Regarding HCN isotopologues, none are clearly detected but upper limits or marginal signal at the $2-3-\sigma$ level are compatible with values observed in other comets (Biver et al. 2024). Note that for ¹⁵N/¹⁴N, the Earth value of 272 corresponds already to an enrichment compared to the estimated solar value of 450 (Marty et al. 2011), showing that fractionation was important already for Earth nitrogen but even more for the material that was incorporated into cometary ices 4.5 billion years ago.

The ³²Sf/³⁴S could be measured in CS in both comets, but in comet C/2022 E3 it seems lower than the Earth value while compatible with terrestrial value for H₂S: this comet seems enriched in C³⁴S, at the $3-\sigma$ level. However the PSW observing mode used for comet C/2022 E3 is more prone to produce ripples that affect the signal of marginal lines. We also provide upper limits on abundances of deuterated species. None yields a D/H below the Earth VSMOW (Vienna Standard Mean Ocean Water) value, but since some enrichment has been observed in the interstellar medium and in some molecules ejected by comet 67P/Churyumov-Gerasimenko observed by the ROSINA mass spectrometer on board the Rosetta spacecraft (Drozdovskaya et al. 2021; Müller et al. 2022), we provide those upper limits.

## 5. The abundance of HNCO and HCOOH in comets

The molecules HNCO and HCOOH were first detected in comets C/1996 B2 (Hyakutake) and C/1995 O1 (Hale-Bopp) in 1996–1997 (Lis et al. 1997; Bockelée-Morvan et al. 2000). They have been detected in several comets since (Biver et al. 2023a) but
<table><caption>**Table 3.** Molecular abundances.</caption><thead><tr><th rowspan="2">Molecule</th><th rowspan="2">Name</th><th colspan="3">Abundance relative to water in %</th></tr><tr><th>C/2021 A1 (Leonard)</th><th>C/2022 E3 (ZTF)</th><th>All comets</th></tr></thead><tbody><tr><td>HCN</td><td>Hydrogen cyanide</td><td>0.09 ± 0.01</td><td>0.09 ± 0.01</td><td>0.08–0.25</td></tr><tr><td>HNC</td><td>Hydrogen isocyanide</td><td>0.005 ± 0.0005</td><td>0.0015 ± 0.0003</td><td>0.0015–0.035</td></tr><tr><td>HNC<sub>d</sub><sup>(a)</sup></td><td></td><td>0.006 ± 0.0007</td><td>0.0032 ± 0.0006</td><td></td></tr><tr><td>CH<sub>3</sub>CN</td><td>Methyl cyanide</td><td>0.016 ± 0.001</td><td>0.017 ± 0.001</td><td>0.008–0.054</td></tr><tr><td>HC<sub>3</sub>N</td><td>Cyanoacetylene</td><td>0.004 ± 0.001</td><td>&lt;0.0022</td><td>0.002–0.068</td></tr><tr><td>HNCO</td><td>Isocyanic acid</td><td>0.073 ± 0.008</td><td>0.042 ± 0.004</td><td>&lt;0.01–0.62</td></tr><tr><td>NH<sub>2</sub>CHO</td><td>Formamide</td><td>0.023 ± 0.002</td><td>0.019 ± 0.001</td><td>0.015–0.022</td></tr><tr><td>CO</td><td>Carbon monoxide</td><td>&lt;2.0</td><td>0.7 ± 0.3</td><td>0.4–35</td></tr><tr><td>H<sub>2</sub>CO</td><td>Formaldehyde</td><td>0.12 ± 0.01</td><td>0.12 ± 0.01</td><td></td></tr><tr><td>H<sub>2</sub>CO<sub>d</sub><sup>(b)</sup></td><td></td><td>0.26 ± 0.01</td><td>0.22 ± 0.01</td><td>0.13–1.4</td></tr><tr><td>CH<sub>3</sub>OH</td><td>Methanol</td><td>0.88 ± 0.02</td><td>1.76 ± 0.01</td><td>0.7–6.1</td></tr><tr><td>HCOOH</td><td>Formic acid</td><td>0.19 ± 0.02</td><td>0.19 ± 0.01</td><td>&lt;0.04–0.58</td></tr><tr><td>CH<sub>3</sub>CHO</td><td>Acetaldehyde</td><td>0.036 ± 0.009</td><td>0.070 ± 0.004</td><td>0.05–0.08</td></tr><tr><td>(CH<sub>2</sub>OH)<sub>2</sub></td><td>Ethylene-glycol</td><td>0.13 ± 0.03</td><td>0.13 ± 0.02</td><td>0.07–0.35</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>Ethanol</td><td>0.18 ± 0.04</td><td>0.17 ± 0.03</td><td>0.11–0.19</td></tr><tr><td>CH<sub>2</sub>OHCHO</td><td>Glycolaldehyde</td><td>0.051 ± 0.012</td><td>&lt;0.024</td><td>0.016–0.039</td></tr><tr><td>CH<sub>2</sub>CO</td><td>Ketene</td><td>0.033 ± 0.015</td><td>&lt;0.03</td><td>≤0.0078</td></tr><tr><td>HCOOCH<sub>3</sub></td><td>Methyl formate</td><td>&lt;0.15</td><td>&lt;0.10</td><td>0.06–0.08</td></tr><tr><td>c-C<sub>2</sub>H<sub>4</sub>O</td><td>Ethylene-oxide</td><td>&lt;0.024</td><td>&lt;0.017</td><td>&lt;0.029</td></tr><tr><td>H<sub>2</sub>S</td><td>Hydrogen sulphide</td><td>0.15 ± 0.02</td><td>0.18 ± 0.04</td><td>0.09–1.5</td></tr><tr><td>CS</td><td>Carbon monosulphide</td><td>0.06 ± 0.01</td><td>0.049 ± 0.004</td><td>0.05–0.20</td></tr><tr><td>CS<sub>d</sub><sup>(c)</sup></td><td></td><td>0.10 ± 0.01</td><td>0.075 ± 0.003</td><td>0.05–0.20</td></tr><tr><td>OCS</td><td>Carbonyl sulphide</td><td>0.11 ± 0.03</td><td>0.068 ± 0.016</td><td>0.05–0.40</td></tr><tr><td>SO</td><td>Sulphur monoxide</td><td>&lt;0.04</td><td>&lt;0.03</td><td>0.04–0.30</td></tr><tr><td>SO<sub>2</sub></td><td>Sulphur dioxide</td><td>&lt;0.04</td><td>&lt;0.03</td><td>0.03–0.23</td></tr><tr><td>H<sub>2</sub>CS</td><td>Thioformaldehyde</td><td>&lt;0.03</td><td>&lt;0.02</td><td>0.009–0.090</td></tr><tr><td>CH<sub>3</sub>SH</td><td>Methyl mercaptan</td><td>&lt;0.08</td><td>&lt;0.06</td><td>&lt;0.06</td></tr><tr><td>c-C<sub>3</sub>H<sub>2</sub></td><td>Cyclopropenylidene</td><td>&lt;0.024</td><td>&lt;0.011</td><td>&lt;0.009</td></tr><tr><td>I-C<sub>3</sub>H<sub>2</sub></td><td>Propadienylidene</td><td>&lt;0.004</td><td>&lt;0.002</td><td>&lt;0.0026</td></tr><tr><td>PH<sub>3</sub></td><td>Phosphine</td><td>&lt;0.16</td><td>&lt;0.09</td><td>&lt;0.07</td></tr></tbody></table>

Notes. Molecules with a “d” are modelled with a daughter distribution with the following parent scale length: <sup>(a)</sup>Assuming a daughter distribution with $L_p = 1000$ and 2000 km (Cordiner et al. 2014, 2017). <sup>(b)</sup>Assuming a daughter distribution with $L_p = 1500$ and 1700 km (fitted values for C/2021 A1 and C/2022 E3). <sup>(c)</sup>Assuming a daughter distribution with $L_p = 4 \times L(CS_2)$.

<table><caption>**Table 4.** Isotopic ratios.</caption><thead><tr><th>Ratio</th><th>Molecule</th><th>C/2021 A1</th><th>C/2022 E3 <sup>(a)</sup></th><th>Other comets <sup>(b)</sup></th><th>Earth <sup>(c)</sup></th></tr></thead><tbody><tr><td><sup>12</sup>C/<sup>13</sup>C</td><td>HCN</td><td>&gt;47</td><td>135 ± 75</td><td>88–114</td><td>89.4</td></tr><tr><td><sup>14</sup>N/<sup>15</sup>N</td><td>HCN</td><td>131 ± 81</td><td>223 ± 138</td><td>139–205</td><td>272</td></tr><tr><td><sup>32</sup>S/<sup>34</sup>S</td><td>CS</td><td>17.8 ± 5.6</td><td>11.2 ± 3.7</td><td>16–23</td><td>22.7</td></tr><tr><td></td><td>H<sub>2</sub>S</td><td>&gt;3.4</td><td>16.5 ± 7.3</td><td>16–23</td><td>22.7</td></tr><tr><td>D/H</td><td>H<sub>2</sub>O</td><td></td><td>&lt;6.4 × 10<sup>-4</sup></td><td>[1.4–6.4] × 10<sup>-4</sup></td><td>1.6 × 10<sup>-4</sup> <sup>(d)</sup></td></tr><tr><td></td><td>HCN</td><td></td><td>&lt;1.8%</td><td>0.23%</td><td></td></tr><tr><td></td><td>H<sub>2</sub>CO</td><td></td><td>&lt;3.1%</td><td>&lt;0.7%</td><td></td></tr><tr><td></td><td>CH<sub>2</sub>DOH</td><td></td><td>&lt;1.7%</td><td>&lt;0.35%</td><td></td></tr><tr><td></td><td>CH<sub>3</sub>OD</td><td></td><td>&lt;2.5%</td><td>&lt;1.0%</td><td></td></tr></tbody></table>

Notes. <sup>(a)</sup>The lines observed in C/2022 E3 are only present at the 1.5–3σ level. <sup>(b)</sup>Bockelée-Morvan et al. (2015); Biver et al. (2016); Cordiner et al. (2019). <sup>(c)</sup>Šimečková et al. (2006). <sup>(d)</sup>VSMOW value.
<table><caption>**Table 5.** Revised abundances of HNCO and HCOOH in comets.</caption><thead><tr><th>Comet</th><th>Dates<br>(yyyy/mm/dd-dd)</th><th>r<sub>h</sub> range<br>(au)</th><th>Q<sub>HNCO</sub>/Q<sub>H<sub>2</sub>O</sub></th><th>Q<sub>HCOOH</sub>/Q<sub>H<sub>2</sub>O</sub></th></tr></thead><tbody><tr><td>C/1996 B2 (Hyakutake)</td><td>1996/03–04</td><td>1.10–0.65</td><td>0.11 ± 0.02%</td><td>&lt;1.4%</td></tr><tr><td>C/1995 O1 (Hale-Bopp)</td><td>1997/02–05</td><td>0.92–1.28</td><td>0.22 ± 0.01%</td><td>0.44 ± 0.03%</td></tr><tr><td>153P/Ikeya-Zhang</td><td>2002/05/08–12</td><td>1.25</td><td>0.049 ± 0.008%</td><td>&lt;0.12%</td></tr><tr><td>C/2004 Q2 (Machholz)</td><td>2005/01/13–17</td><td>1.21</td><td>0.025 ± 0.003%</td><td>0.09 ± 0.02%</td></tr><tr><td>73P-C/Schwassmann-W.3</td><td>2006/05/10–18</td><td>1.00</td><td>0.066 ± 0.014%</td><td>&lt;0.18%</td></tr><tr><td>73P-B/Schwassmann-W.3</td><td>2006/05/10–18</td><td>1.00</td><td>0.075 ± 0.011%</td><td>&lt;0.15%</td></tr><tr><td>8P/Tuttle</td><td>2007/12/29–31</td><td>1.11</td><td>&lt;0.06%</td><td>0.42 ± 0.11%</td></tr><tr><td>C/2007 N3 (Lulin)</td><td>2009/02/26–27</td><td>1.41</td><td>0.025 ± 0.009%</td><td>&lt;0.13%</td></tr><tr><td>103P/Hartley 2</td><td>2010/10/25–29</td><td>1.06</td><td>0.034 ± 0.012%</td><td>&lt;0.17%</td></tr><tr><td>C/2009 P1 (Garradd)</td><td>2012/02/15–19</td><td>1.73</td><td>0.031 ± 0.010%</td><td>&lt;0.14%</td></tr><tr><td>C/2013 R1 (Lovejoy)</td><td>2013/11/08–46</td><td>0.83–1.13</td><td>0.050 ± 0.007%</td><td>0.17 ± 0.03%</td></tr><tr><td>C/2012 S1 (ISON)</td><td>2013/11/14–15</td><td>0.61–0.64</td><td>0.21 ± 0.06%</td><td>0.58 ± 0.09%</td></tr><tr><td>C/2012 F6 (Lemmon)</td><td>2013/03/14–39</td><td>0.75–0.78</td><td>0.10 ± 0.02%</td><td>0.25 ± 0.06%</td></tr><tr><td>C/2011 L4 (PanSTARRS)</td><td>2013/03/14–39</td><td>0.33–0.83</td><td>0.62 ± 0.12%</td><td>≤ 0.93%</td></tr><tr><td>C/2014 Q2 (Lovejoy)</td><td>2015/01/13–26</td><td>1.30</td><td>0.016 ± 0.001%</td><td>0.055 ± 0.004%</td></tr><tr><td>C/2013 US<sub>10</sub> (Catalina)</td><td>2015/12/30–34</td><td>1.18</td><td>0.051 ± 0.019%</td><td>&lt;0.22%</td></tr><tr><td>46P/Wirtanen</td><td>2018/12/11–18</td><td>1.06</td><td>&lt;0.011%</td><td>&lt;0.043%</td></tr><tr><td>67P/Churyumov-G.</td><td>2021/11/12–16</td><td>1.22</td><td>0.085 ± 0.024%</td><td>&lt;0.33%</td></tr><tr><td>C/2021 A1 (Leonard)</td><td>2021/12/08–13</td><td>0.76–0.81</td><td>0.073 ± 0.008%</td><td>0.19 ± 0.02%</td></tr><tr><td>C/2022 E3 (ZTF)</td><td>2023/02/03–07</td><td>1.18</td><td>0.042 ± 0.004%</td><td>0.19 ± 0.02%</td></tr></tbody></table>

there were large uncertainties on their destruction rates and lifetimes and a simple excitation model was used to derive their abundances. Hrodmarsson & van Dishoeck (2023) have recently published photo-dissociation and photo-ionisation rates for HNCO and HCOOH that significantly differ from those assumed in the past (Biver et al. 2021a, and references therein). The new values ($\beta_0$(HNCO) = $38 \times 10^{-5}$ s$^{-1}$, $\beta_0$(HCOOH) = $54 \times 10^{-5}$ s$^{-1}$, in the solar radiation field at 1 au) are 5 to 10 times higher than previously assumed, with a stated 2-$\sigma$ uncertainty of the order of 20–30% (Hrodmarsson & van Dishoeck 2023). We have decided to revisit all millimetre observations of comets reporting a detection of at least one of these molecules and recomputed production rates and abundances relative to water assuming values of:

$$- \beta_0(\text{HNCO}) = 3.5 \times 10^{-4} \text{ s}^{-1};$$

$$- \beta_0(\text{HCOOH}) = 5.0 \times 10^{-4} \text{ s}^{-1}.$$

The contribution of Solar Lyman-$\alpha$ to their photo-dissociation is minor and their photo-dissociation cross-section extends well into the 130–200 nm range where the solar flux is less affected by solar activity (Huebner et al. 1992). Therefore, we do not expect a large variation with solar activity. An overview of the observations analysed in this work, partly published in previous papers is given in Table A.3.

Table 5 provides the abundances relative to water of HNCO and HCOOH in these comets. In these calculations infrared pumping via the rotational bands is not taken into account, but with those reduced lifetimes, the photo-dissociation of the molecules takes place in a region where collisional excitation dominates (within ~2000 km from the nucleus) and the photo-dissociation process is as efficient as vibrational excitation (g-factors comparable to $\beta$). So the new reduced lifetimes (compared to previous calculations) should also reduce the uncertainties from the neglected infrared pumping.

HNCO is relatively well detected in the two comets via its $K_a$ = 0 lines at 2 mm and 1 mm wavelengths, and marginal signals at an offset of 10" (Tables B.5, B.6) suggest that it may be produced by a distributed source with a scale length

$L_p$(HNCO) larger than 2000 km. Indeed, assuming a parent scale length of 5000 km also provides a better agreement between the observed rotational temperature and the predicted one for $T = 60$ K (Table 1).

HCOOH: In Tables B.5 and B.6 we have selected the lower energy lines that are expected to be the strongest, that is those with $K_a$ = 0 to 2. For comet C/2021 A1, since at the offset position the weighted average intensity is similar to the one at the central position, we do not find a solution for a reasonably small parent scale length, but data are compatible within $1\sigma$ with $L_p$(HCOOH) = 2000 km, and $L_p$(HCOOH) = 0 km cannot be fully excluded. Similarly, for comet C/2022 E3, we derive $L_p$(HCOOH) = 0–5200 km in the $\pm 1\sigma$ interval, so that distributed source as well as nuclear source are not excluded. Hence, if HCOOH is produced by the sublimation and dissociation of ammonium formate, then this must take place within a few thousands of kilometers from the nucleus at 0.78–1.18 au from the Sun.

The revised measured abundances in comets (increased by a factor 1.2 to 5) of these two molecules are plotted in Fig. 17 as a function of heliocentric distance at which they were detected. This plots suggest that their abundance relative to water increases for comets observed closer to the Sun. This needs to be confirmed, but can suggest that the molecules could be at least in part produced by the degradation of parents in the coma such as ammonium salts. Their combined abundance relative to water in comets (~0.1–1%) is comparable to that of NH$_3$(0.2–5%, Biver et al. 2024), the other decomposition product of ammonium salts. Interestingly, a similar trend with heliocentric distance has been observed for the abundance of ammonia (see, e.g. Dello Russo et al. 2016).

## 6. Conclusions

We have undertaken a comprehensive and sensitive spectroscopic survey at 2 mm and 1 mm of two long period comets,
<figure><img src="image_18.png" /><figcaption>**Fig. 17.** Abundances relative to water of HNCO (green) and HCOOH (blue) in comets as a function of the heliocentric distance. Even though there is dispersion between comets, the trend appears to be that these species are more abundant in the coma close to the Sun.</figcaption></figure>

C/2021 A1 (Leonard) and C/2022 E3 (ZTF) when they were relatively close to the Earth (0.2–0.3 au). We collected in addition some useful spatial information from coarse mapping, to probe the spatial distribution of selected molecules. The main results can be summarised as follows:

- Both comets exhibited relatively stable (within ±20%) production rates during the observational periods, although C/2021 A1 underwent a series of disintegration outbursts just after our IRAM-30m observing run;

- We determined the abundance relative to water or a significant upper limit for 26 molecules. Comets C/2021 A1 and C/2022 E3 share similar compositions: they are depleted in hypervolatiles (low CO and H₂S abundances relative to water) and have a relatively low methanol abundance relative to water (0.9 and 1.8% respectively) compared to the mean of other comets;

- Constraints on the presence of distributed sources have been obtained. A slightly distributed source (scale length smaller than 1/3 beam) for HCN fits better the data than direct release from the nucleus. The parent source of formaldehyde was found to have a scale length $L_p(H_2CO) = 0.68$ and $0.36 \pm 0.11$ times the photo-dissociation length of formaldehyde ($L_d(H_2CO)$) in C/2021 A1 and C/2022 E3 respectively. We measured $L_p(CS) \approx 4 \times L_d(CS_2)$, suggesting that the dissociation rate of the parent producing CS is $\beta_0 = 4.5 \times 10^{-4} \text{ s}^{-1}$ at 1 au from the Sun. This is comparable to the photo-dissociation rate of H₂CS, but the H₂CS upper limit is more than 3 times lower than the CS abundance;

- Both comets are relatively depleted in sulphur-bearing species compared to other comets, with H₂S/H₂O in the low range (0.15–0.18%) and the sum of sulphur-bearing molecule abundances below ~0.5%. SO and SO₂ are not detected with abundances below the lowest measured in a comet and OCS is detected in both comets but with a low abundance, too;

- We obtained only shallow constraints on isotopic ratios: ¹²C/¹³C, ¹⁴N/¹⁵N, and ³²S/³⁴S are compatible with values observed in other comets (90, 150, and 23, respectively), although the C³²S/C³⁴S seems marginally twice lower than the terrestrial value in comet C/2022 E3;

- HNCO and HCOOH are well detected in both comets. Because of recently published lifetimes for these molecules

that are significantly shorter than previous estimates, we have revised their inferred abundances in all comets for which detections of these species have been reported. The values found in C/2021 A1 (0.07% and 0.19%) and C/2022 E3 (0.04% and 0.19%) are close to their median abundances relative to water in other comets, 0.05% and 0.25%, respectively. Their total abundance is also about half the median abundance of NH₃ in comets. We cannot exclude that these species are produced in the coma by a distributed source, and their abundance in the comae seems to increase at lower heliocentric distances. This favours the hypothesis that HNCO, HCOOH, and consequently an important fraction of NH₃ could come from the decomposition of ammonium salts.

These observations have provided complementary and new insights regarding the source of several cometary molecules that are unlikely present as such in cometary ices. Further investigations at higher spatial resolution (e.g. interferometric investigations undertaken on those comets under reduction or on future comets), will be very useful to pin down the production processes and the scale lengths of these molecules (CS, H₂CO, HNCO, HCOOH).

## Data availability

The radio spectra are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/690/A271. The tables of line intensities are available at https://doi.org/10.5281/zenodo.13332037. The spectra of comets are available at https://doi.org/10.5281/zenodo.13341404.

*Acknowledgements.* This work is based on observations carried out under projects number 001-21, 100-21, and 097-22 with the IRAM 30-m telescope. 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 data were reduced and analysed thanks to the use of the GILDAS, class software (https://www.iram.fr/IRAMFR/GILDAS). This research has been supported by the Programme national de planétologie de l'Universités des sciences de l'univers (INSU). The Nançay radio telescope (NRT) is part of the Observatoire radioastronomique de Nançay, operated by the Paris Observatory, associated with the French Centre national de la recherche scientifique (CNRS) and with the University of Orléans. 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 (80NM0018D0004). D.C.L. acknowledges financial support from the National Aeronautics and Space Administration (NASA) Astrophysics Data Analysis Program (ADAP). M.N. Drozdovskaya acknowledges the Holcim Foundation Stipend. B.P. Bonev and N. Dello Russo acknowledge support of NSF grant AST-2009398 and NASA grant 80NSSC17K0705, respectively. N.X. Roth was supported by the NASA Postdoctoral Program at the NASA Goddard Space Flight Center, administered by Universities Space Research Association under contract with NASA. M.A. Cordiner was supported in part by the National Science Foundation (under Grant No. AST-1614471). S.N. Milam and M.A. Cordiner acknowledge the Planetary Science Division Internal Scientist Funding Program through the Fundamental Laboratory Research (FLaRe) work package, as well as the NASA Astrobiology Institute through the Goddard Center for Astrobiology (proposal 13-13NAI7-0032).

## 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., Hänni, N., et al. 2020, *Nat. Astron.*, **4**, 533

Biver, N., & Bockelée-Morvan, D. 2019, *ACS Earth Space Chem.*, **3**, 1550
Biver, N., Bockelée-Morvan, D., Crovisier, J., et al. 1999, AJ, 118, 1850

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

Biver N., Bockelée-Morvan D., Colom P., et al. 2011, A&A, 528, A142

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

Biver, N., Bockelée-Morvan, D., Hofstadter, M., et al. 2019, A&A, 630, A19

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

Biver, N., Bockelée-Morvan, D., Lis, D. C., et al. 2021b, A&A, 651, A25

Biver, N., Bockelée-Morvan, D., Boissier, J., et al. 2022, A&A, 668, A171

Biver, N., Bockelée-Morvan, D., Crovisier, J., et al. 2023a, A&A, 670, A170

Biver, N., Bockelée-Morvan, D., Crovisier, J., et al. 2023b, CBET, 5216

Biver N., Dello Russo N., Opitom C., Rubin, M. 2024, arXiv e-prints
[arXiv:2207.04800]

Bockelée-Morvan, D., Lis, D. C., Wink, J. E., et al. 2000, A&A, 353, 1101

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

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

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

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

Combi, M., Mäkinen, T., Bertaux, J.-L., Quemerais, E., & Ferron, S. 2023a, Icarus, 398, 115543

Combi, M., Mäkinen, T., Bertaux, J.-L., Quemerais, E., & Ferron, S. 2023b, DPS-EPSC, 322.04 (on-line poster)

Cordiner, M. A., Remijan, A. J., Boissier, J., et al. 2014, ApJ, 792, L2

Cordiner, M. A., Boissier, J., Charnley, S. B., et al. 2017, ApJ, 838, 147

Cordiner, M. A., Palmer, M. Y., de Val-Borro, M., et al. 2019, ApJ, 870, L26

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

Crovisier, J., Biver, N., & Bockelée-Morvan, D. 2021, CBET, 5087

Dello Russo N., Kawakita, H., Vervack, R. J., & Weaver, H. A. 2016, Icarus, 278, 301

Drozdovskaya, M. N., Schroeder, I I. R. H. G., Rubin, M., et al. 2021, MNRAS, 500, 4901

Faggi, S., Lippi, M., Mumma, M. J., & Villanueva, G. L. 2023, *Planetary Science Journal*, **4**, 8

Frisk, U., Hagström, M., Ala-Laurinaho, J., et al. 2003, A&A, 402, L27

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

Heays, A. N., Bosman, A. D., & van Dishoeck, E. F. 2017, A&A, 602, A105

Hrodmarsson, H. R., & van Dishoeck, E. F. 2023, A&A, 675, A25

Huebner, W. F., & Mukherjee, J. 2015, Planet. Space Sci., 106, 11

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

Jehin, E., Moulane, Y., & Manfroid, J. 2021, CBET, 5087

Jewitt, D., Kim, Y., Mattiazzo, M., et al. 2023, AJ, 165, 122

Knight, M. M., Holt, C. E., Villa, K. M., Skiff, B. A., & Schleicher, D. G. 2023, ATel, 15879

Lecacheux, A., Biver, N., Crovisier, J., et al. 2003, A&A, 402, L55

Lippi, M., Villanueva, G. L., Mumma, M. J., & Faggi, S. 2021, AJ, 162, 74

Lis, D. C., Keene, J., Young, K., et al. 1997, Icarus, 130, 355

Manzini, F., Oldani, V., Ochner, P., Bedin, L. R., & Reguitti, A. 2023, ATel, 15909

Marty, B., Chaussidon, M., Wiens, R. C., Jurewicz, A. J. G., & Burnett, D. S. 2011, Science, 332, 1533

Milam, S. N., Roth, N. X., Villanueva, G. L., et al. 2023, ACM 2023, LPI Contrib. No. 2851

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

Müller, D. R., Altwegg, K., Berthelier, J. J., et al. 2022, A&A, 662, A69

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

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

Roth, N. X., Milam, S. N., Cordiner, M. A., et al. 2021, ApJ, 921, 14

Šimečková, M., Jacquemart, D., Rothman, L. S., Gamache, R. R., & Goldman, A. 2006, J. Quant. Spec. Radiat. Transf., 98, 130
## Appendix A: Additional material

<table><caption>Table A.1. Log of observations of comet C/2021 A1 (Leonard) in 2021.</caption><thead><tr><th>UT date<br>(mm/dd.d–dd.d)</th><th>&lt;r<sub>h</sub> &gt;<br>(au)</th><th>&lt;Δ &gt;<br>(au)</th><th>Phase<sup>(a)</sup></th><th>Tel.</th><th>Integ. time<br>(min)<sup>(c)</sup></th><th>pwv<sup>(b)</sup><br>(mm)</th><th>Mode<sup>(d)</sup></th><th>Freq. range<br>(GHz)</th></tr></thead><tbody><tr><td>11/12.18–12.20</td><td>1.218</td><td>1.192</td><td>48.5°</td><td>IRAM</td><td>28</td><td>4.8</td><td>WSW</td><td>248.7-256.5, 264.4-272.2</td></tr><tr><td>11/13.19–13.22</td><td>1.203</td><td>1.157</td><td>49.5°</td><td>IRAM</td><td>33</td><td>6.6</td><td>WSW</td><td>146.9-154.7, 162.6-170.4</td></tr><tr><td>11/14.16–14.20</td><td>1.188</td><td>1.123</td><td>50.6°</td><td>IRAM</td><td>34</td><td>4.4</td><td>PSW</td><td>248.7-256.5, 264.4-272.2</td></tr><tr><td>14.22–14.24</td><td>1.187</td><td>1.121</td><td>50.7°</td><td>IRAM</td><td>16</td><td>5.0</td><td>PSW</td><td>209.7-217.5, 225.4-233.1</td></tr><tr><td>11/15.18–15.22</td><td>1.172</td><td>1.087</td><td>51.8°</td><td>IRAM</td><td>36</td><td>3.1</td><td>PSW</td><td>146.9-154.7, 162.6-170.4</td></tr><tr><td>11/16.15–16.20</td><td>1.157</td><td>1.053</td><td>52.9°</td><td>IRAM</td><td>44</td><td>1.3</td><td>PSW</td><td>248.7-256.5, 264.4-272.2</td></tr><tr><td>12/08.44–08.47</td><td>0.826</td><td>0.288</td><td>115.5°</td><td>IRAM</td><td>42</td><td>2.0</td><td>WSW</td><td>248.4-256.5, 264.4-272.5</td></tr><tr><td>08.52–08.58</td><td>0.824</td><td>0.286</td><td>116.4°</td><td>IRAM</td><td>64</td><td>2.5</td><td>WSW</td><td>146.6-154.8, 162.6-170.7</td></tr><tr><td>12/09.41–09.58</td><td>0.813</td><td>0.266</td><td>124.0°</td><td>IRAM</td><td>150</td><td>4-6</td><td>WSW</td><td>248.4-256.5, 264.4-272.5</td></tr><tr><td>12/10.29–10.44</td><td>0.800</td><td>0.250</td><td>131.8°</td><td>IRAM</td><td>136</td><td>4-7</td><td>WSW</td><td>248.4-256.5, 264.4-272.5</td></tr><tr><td>10.46–10.54</td><td>0.798</td><td>0.248</td><td>133.0°</td><td>IRAM</td><td>124</td><td>8-19</td><td>WSW</td><td>146.6-154.8, 162.6-170.7</td></tr><tr><td>10.54–10.59</td><td>0.798</td><td>0.248</td><td>133.5°</td><td>IRAM</td><td>124</td><td>4-8</td><td>WSW</td><td>146.6-154.8, 162.6-170.7</td></tr><tr><td>12/11.51–11.53</td><td>0.785</td><td>0.237</td><td>142.7°</td><td>IRAM</td><td>28</td><td>6.0</td><td>WSW</td><td>82.9-90.9, 224.8-232.8</td></tr><tr><td>11.57–11.59</td><td>0.785</td><td>0.237</td><td>143.5°</td><td>IRAM</td><td>28</td><td>5.7</td><td>WSW</td><td>248.4-256.5, 264.4-272.5</td></tr><tr><td>11.59–11.66</td><td>0.784</td><td>0.237</td><td>143.7°</td><td>IRAM</td><td>65</td><td>3-5</td><td>WSW</td><td>248.4-256.5, 264.4-272.5</td></tr><tr><td>12/12.38–12.41</td><td>0.774</td><td>0.233</td><td>150.7°</td><td>IRAM</td><td>9</td><td>2.8</td><td>WSW</td><td>248.4-256.5, 264.4-272.5</td></tr><tr><td>12.42–12.54</td><td>0.773</td><td>0.233</td><td>151.4°</td><td>IRAM</td><td>113</td><td>1.8</td><td>WSW</td><td>240.0-248.1, 256.0-264.1</td></tr><tr><td>12.56–12.65</td><td>0.772</td><td>0.233</td><td>152.5°</td><td>IRAM</td><td>109</td><td>2.0</td><td>WSW</td><td>209.4-217.5, 225.4-233.4</td></tr><tr><td>12/13.38–13.41</td><td>0.762</td><td>0.236</td><td>158.0°</td><td>IRAM</td><td>19</td><td>0.7</td><td>WSW</td><td>248.4-256.5, 264.4-272.5</td></tr><tr><td>13.43–13.50</td><td>0.761</td><td>0.236</td><td>158.4°</td><td>IRAM</td><td>79</td><td>0.4</td><td>WSW</td><td>217.5-225.5, 234.5-241.5</td></tr><tr><td>13.52–13.57</td><td>0.760</td><td>0.237</td><td>158.7°</td><td>IRAM</td><td>56</td><td>1.0</td><td>WSW</td><td>240.0-248.1, 256.0-264.1</td></tr><tr><td>13.60–13.66</td><td>0.759</td><td>0.237</td><td>159.1°</td><td>IRAM</td><td>62</td><td>2.1</td><td>WSW</td><td>146.6-154.8, 162.6-170.7</td></tr></tbody></table>

Notes. (a) Phase angle. (b) Mean precipitable water vapour in the atmosphere above the telescope.

(c) Total (offset positions included) integration time (ON+OFF) on the source. (d) Observing mode: WSW = Wobbler Switching (reference at ±3'); PSW = Position Switching (reference at 5').
<table><caption>**Table A.2.** Log of observations of comet C/2022 E3 (ZTF) in 2023.</caption><thead><tr><th>UT date<br>(mm/dd.d–dd.d)</th><th>&lt;r<sub>h</sub> &gt;<br>(au)</th><th>&lt;Δ &gt;<br>(au)</th><th>Phase<sup>(a)</sup></th><th>Tel.</th><th>Integ. time<br>(min)<sup>(c)</sup></th><th>pwv<sup>(b)</sup><br>(mm)</th><th>Mode<sup>(d)</sup></th><th>Freq. range<br>(GHz)</th></tr></thead><tbody><tr><td>01/19.33–19.42</td><td>1.117</td><td>0.523</td><td>61.7°</td><td>Odin</td><td>81</td><td>0</td><td>DSW</td><td>556.43-557.44</td></tr><tr><td>01/19.45–19.75</td><td>1.118</td><td>0.516</td><td>61.7°</td><td>Odin</td><td>233</td><td>0</td><td>DSW</td><td>556.43-557.44, map 3x3, step 1'</td></tr><tr><td>01/19.78–19.88</td><td>1.118</td><td>0.511</td><td>61.6°</td><td>Odin</td><td>85</td><td>0</td><td>DSW</td><td>556.43-557.44</td></tr><tr><td>01/19.98–20.14</td><td>1.119</td><td>0.504</td><td>61.6°</td><td>Odin</td><td>108</td><td>0</td><td>DSW</td><td>556.43-557.44, map 3x3, step 2'</td></tr><tr><td>01/20.19–20.28</td><td>1.119</td><td>0.501</td><td>61.5°</td><td>Odin</td><td>72</td><td>0</td><td>DSW</td><td>556.43-557.44</td></tr><tr><td>02/03.67–03.78</td><td>1.168</td><td>0.291</td><td>45.4°</td><td>IRAM</td><td>93</td><td>0.8</td><td>PSW</td><td>248.3-256.4, 264.3-272.4</td></tr><tr><td>02/04.65–04.69</td><td>1.173</td><td>0.300</td><td>45.5°</td><td>IRAM</td><td>39</td><td>0.5</td><td>PSW</td><td>248.3-256.4, 264.3-272.4</td></tr><tr><td>04.71–04.78</td><td>1.173</td><td>0.301</td><td>45.5°</td><td>IRAM</td><td>65</td><td>0.2</td><td>PSW</td><td>240.1-248.2, 256.1-264.2</td></tr><tr><td>04.81–04.86</td><td>1.174</td><td>0.302</td><td>45.5°</td><td>IRAM</td><td>56</td><td>0.3</td><td>PSW</td><td>240.1-248.2, 256.1-264.2</td></tr><tr><td>04.88–04.95</td><td>1.174</td><td>0.303</td><td>45.5°</td><td>IRAM</td><td>76</td><td>0.2</td><td>PSW</td><td>146.8-154.9, 162.8-170.9</td></tr><tr><td>04.97–05.02</td><td>1.174</td><td>0.304</td><td>45.6°</td><td>IRAM</td><td>59</td><td>0.2</td><td>PSW</td><td>209.4-217.5, 225.4-233.5</td></tr><tr><td>05.03–05.04</td><td>1.175</td><td>0.304</td><td>45.6°</td><td>IRAM</td><td>16</td><td>1.1</td><td>PSW</td><td>161.0-169.1, 177.0-185.4</td></tr><tr><td>02/05.65–05.69</td><td>1.178</td><td>0.312</td><td>45.8°</td><td>IRAM</td><td>40</td><td>1.8</td><td>PSW</td><td>248.3-256.4, 264.3-272.4</td></tr><tr><td>05.71–05.81</td><td>1.178</td><td>0.314</td><td>45.8°</td><td>IRAM</td><td>94</td><td>0.9</td><td>PSW</td><td>217.5-225.5, 233.5-241.5</td></tr><tr><td>05.84–05.89</td><td>1.179</td><td>0.315</td><td>45.9°</td><td>IRAM</td><td>58</td><td>0.5</td><td>PSW</td><td>146.8-154.9, 162.8-170.9</td></tr><tr><td>05.91–06.00</td><td>1.179</td><td>0.317</td><td>45.9°</td><td>IRAM</td><td>95</td><td>0.4</td><td>PSW</td><td>209.4-217.5, 225.4-233.5</td></tr><tr><td>06.03–06.04</td><td>1.180</td><td>0.318</td><td>46.0°</td><td>IRAM</td><td>17</td><td>0.8</td><td>PSW</td><td>248.3-256.4, 264.3-272.4</td></tr><tr><td>02/06.66–06.69</td><td>1.183</td><td>0.328</td><td>46.0°</td><td>IRAM</td><td>31</td><td>1.7</td><td>WSW</td><td>248.3-256.4, 264.3-272.4</td></tr><tr><td>06.69–06.70</td><td>1.183</td><td>0.328</td><td>46.0°</td><td>IRAM</td><td>8</td><td>1.6</td><td>PSW</td><td>248.3-256.4, 264.3-272.4</td></tr><tr><td>06.71–06.73</td><td>1.183</td><td>0.329</td><td>46.0°</td><td>IRAM</td><td>23</td><td>2.0</td><td>WSW</td><td>248.3-256.4, 264.3-272.4</td></tr><tr><td>06.75–06.80</td><td>1.184</td><td>0.329</td><td>46.0°</td><td>IRAM</td><td>48</td><td>3.5</td><td>WSW</td><td>240.1-248.2, 256.1-264.2</td></tr><tr><td>06.88–06.90</td><td>1.184</td><td>0.331</td><td>46.0°</td><td>IRAM</td><td>23</td><td>1.7</td><td>WSW</td><td>209.4-217.5, 225.4-233.5</td></tr><tr><td>06.92–07.00</td><td>1.185</td><td>0.333</td><td>46.0°</td><td>IRAM</td><td>91</td><td>2.5</td><td>WSW</td><td>209.4-217.5, 225.4-233.5</td></tr><tr><td>07.02–07.04</td><td>1.185</td><td>0.334</td><td>46.0°</td><td>IRAM</td><td>25</td><td>1.9</td><td>WSW</td><td>146.8-154.9, 162.8-170.9</td></tr></tbody></table>

**Notes.** (a) Phase angle. (b) Mean precipitable water vapour in the atmosphere above the telescope.

(c) Total (offset positions included) integration time (ON+OFF) on the source. (d) Observing mode: WSW = Wobbler Switching (reference at ±2'); PSW = Position Switching (reference at 4'), DSW = Dicke Switching (reference at ~42°).
<table><caption>**Table A.3.** Observations of HNCO and HCOOH in comets.</caption><thead><tr><th rowspan="2">Comet</th><th colspan="4">HNCO</th><th colspan="4">HCOOH</th></tr><tr><th>(yyyy/mm/dd-dd)</th><th>r<sub>h</sub> (au)</th><th>telescopes<sup>(a)</sup></th><th>lines<sup>(b)</sup></th><th>(yyyy/mm/dd-dd)</th><th>r<sub>h</sub> (au)</th><th>telescopes<sup>(a)</sup></th><th>lines<sup>(b)</sup></th></tr></thead><tbody><tr><td>C/1996 B2<br>(Hyakutake)</td><td>1996/03–04</td><td>1.10–0.65</td><td>CSO,<br>PdB,IRAM</td><td>1(+2)</td><td>1996/04/11</td><td>0.65</td><td>IRAM</td><td>(2)</td></tr><tr><td>C/1995 O1<br>(Hale-Bopp)</td><td>1997/02–05</td><td>0.92–1.28</td><td>CSO,<br>PdB,IRAM</td><td>4(+1)</td><td>1997/02–05</td><td>0.92–1.28</td><td>CSO,<br>PdB,IRAM</td><td>7(+1)</td></tr><tr><td>153P/Ikeya-Zhang</td><td>2002/05/12</td><td>1.255</td><td>IRAM</td><td>1</td><td>2002/05/10–12</td><td>1.23</td><td>IRAM</td><td>(2)</td></tr><tr><td>C/2004 Q2 (Machholz)</td><td>2005/01/16–18</td><td>1.21</td><td>IRAM</td><td>1</td><td>2005/01/15–18</td><td>1.21</td><td>IRAM</td><td>(3)</td></tr><tr><td>73P-C/S.-W.3</td><td>2006/05/09–18</td><td>1.00</td><td>CSO, IRAM</td><td>1(+1)</td><td>2006/05/13–18</td><td>1.00</td><td>IRAM</td><td>(3)</td></tr><tr><td>73P-B/S.-W.3</td><td>2006/05/09–18</td><td>1.00</td><td>CSO, IRAM</td><td>2</td><td>2006/05/13–18</td><td>1.00</td><td>IRAM</td><td>(3)</td></tr><tr><td>8P/Tuttle</td><td>2007/12/29–31</td><td>1.10</td><td>IRAM</td><td>1</td><td>2007/12/29–31</td><td>1.11</td><td>IRAM</td><td>1(+2)</td></tr><tr><td>C/2007 N3 (Lulin)</td><td>2009/02/26–27</td><td>1.41</td><td>IRAM</td><td>1</td><td>2009/02/25–27</td><td>1.41</td><td>IRAM</td><td>(3)</td></tr><tr><td>103P/Hartley 2</td><td>2010/10/25–29</td><td>1.06</td><td>IRAM, CSO</td><td>1(+3)</td><td>2010/10/25–39</td><td>1.06</td><td>IRAM, CSO</td><td>(9)</td></tr><tr><td>C/2009 P1 (Garradd)</td><td>2012/02/15–19</td><td>1.73</td><td>IRAM</td><td>1(+6)</td><td>2012/02/15–19</td><td>1.73</td><td>IRAM</td><td>(10)</td></tr><tr><td>C/2013 R1 (Lovejoy)</td><td>2013/11/08–16</td><td>1.13</td><td>IRAM</td><td>(5)</td><td>2013/11/08–16</td><td>1.13</td><td>IRAM</td><td>(21)</td></tr><tr><td></td><td>2013/11/27–31</td><td>0.92</td><td>IRAM</td><td>1(+5)</td><td>2013/11/27–31</td><td>0.92</td><td>IRAM</td><td>7(+13)</td></tr><tr><td></td><td>2013/12/09–16</td><td>0.83</td><td>IRAM</td><td>1(+9)</td><td>2013/12/09–16</td><td>0.82</td><td>IRAM</td><td>2(+24)</td></tr><tr><td>C/2012 S1 (ISON)</td><td>2013/11/14–15</td><td>0.62</td><td>IRAM</td><td>1</td><td>2013/11/14–15</td><td>0.62</td><td>IRAM</td><td>2(+6)</td></tr><tr><td>C/2012 F6 (Lemmon)</td><td>2013/03/14–18</td><td>0.75</td><td>IRAM</td><td>(6)</td><td>2013/03/14–39</td><td>0.75</td><td>IRAM</td><td>(16)</td></tr><tr><td></td><td>2013/04/06–08</td><td>0.75</td><td>IRAM</td><td>1(+8)</td><td>2013/04/06–08</td><td>0.78</td><td>IRAM</td><td>(22)</td></tr><tr><td>C/2011 L4<br>(PanSTARRS)</td><td>2013/03/14–18</td><td>0.36</td><td>IRAM</td><td>2(+4)</td><td>2013/03/14–18</td><td>0.36</td><td>IRAM</td><td>1(+21)</td></tr><tr><td>C/2014 Q2 (Lovejoy)</td><td>2013/04/05–08</td><td>0.81</td><td>IRAM</td><td>(9)</td><td>2013/04/05–08</td><td>0.81</td><td>IRAM</td><td>(28)</td></tr><tr><td></td><td>2015/01/13–16</td><td>1.31</td><td>IRAM</td><td>4(+2)</td><td>2015/01/13–16</td><td>1.31</td><td>IRAM</td><td>6(+14)</td></tr><tr><td></td><td>2015/01/23–26</td><td>1.30</td><td>IRAM</td><td>4(+1)</td><td>2015/01/23–26</td><td>1.29</td><td>IRAM</td><td>3(+9)</td></tr><tr><td>C/2013 US<sub>10</sub> (Catalina)</td><td>2015/12/30–34</td><td>1.18</td><td>IRAM</td><td>(6)</td><td>2015/12/30–34</td><td>1.18</td><td>IRAM</td><td>(10)</td></tr><tr><td>46P/Wirtanen</td><td>2018/12/11–18</td><td>1.06</td><td>IRAM</td><td>(10)</td><td>2018/12/11–18</td><td>1.06</td><td>IRAM</td><td>(17)</td></tr><tr><td>67P/Churyumov-G.</td><td>2021/11/12–16</td><td>1.22</td><td>IRAM</td><td>(10)</td><td>2021/11/12–16</td><td>1.22</td><td>IRAM</td><td>(14)</td></tr><tr><td>C/2021 A1 (Leonard)</td><td>2021/12/08–13</td><td>0.78</td><td>IRAM</td><td>6(+6)</td><td>2021/12/08–13</td><td>0.78</td><td>IRAM</td><td>2(+26)</td></tr><tr><td>C/2022 E3 (ZTF)</td><td>2023/02/03–07</td><td>1.18</td><td>IRAM</td><td>3(+9)</td><td>2023/02/03–07</td><td>1.18</td><td>IRAM</td><td>13(+17)</td></tr></tbody></table>

Notes. <sup>(a)</sup>Radio telescopes used: IRAM 30-m, CSO = Caltech Submillimeter Observatory 10.4-m, PdB = Plateau de Bure interferometer in single-dish ON-OFF mode.

<sup>(b)</sup>The number of lines in parenthesis have a S/N < 3, but are used to infer abundances.
<figure><img src="image_19.png" /></figure>

**Fig. A.1.** CH₃CN series of four lines observed with IRAM-30m in comet C/2021 A1 on 8–13 December 2021 (average intensity). The vertical axis is main beam brightness temperature in K, horizontal axis is Doppler velocity in the rest frame of the comet with respect to the (J, 0) – (J − 1, 0) line, with frequencies in the rest frame of the comet given on the upper axis.
<figure><img src="image_20.png" /><figcaption>**Fig. A.2.** CH₃OH series of lines around 166, 242, and 252 GHz observed with IRAM-30m in comet C/2021 A1 on 8–13 December 2021 (average intensity). The vertical axis is main beam brightness temperature in K, horizontal axis is Doppler velocity in the rest frame of the comet with respect to the main line, with frequencies in the rest frame of the comet given on the upper axis.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Fig. A.3.** Spectra of molecules obtained by averaging several lines observed with IRAM-30m in comet C/2021 A1 between 8 and 13 December 2021. The vertical axis is main beam brightness temperature in K. The horizontal axis is the Doppler velocity in the rest frame of the comet. The number of lines averaged is provided for each molecule, either in the 2 mm band (147–153 and 163–171 GHz) or 1 mm band (209–272 GHz).</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>**Fig. A.4.** CH₃CN series of lines at 147, 165, 220, and 257 GHz observed with IRAM-30m in comet C/2022 E3 on 3–6 February 2023 (average intensity). The vertical axis is main beam brightness temperature in K. The horizontal axis is the Doppler velocity in the rest frame of the comet with respect to the main line, with frequencies in the rest frame of the comet given on the upper axis.</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>**Fig. A.5.** CH₃OH series of lines around 165, 242, and 252 GHz observed with IRAM-30m in comet C/2022 E3 on 3–6 February 2023 (average intensity). The vertical axis is main beam brightness temperature in K. The horizontal axis is the Doppler velocity in the rest frame of the comet with respect to the main line, with frequencies in the rest frame of the comet given on the upper axis.</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>**Fig. A.6.** Spectra of molecules obtained by averaging several lines observed with IRAM-30m in comet C/2022 E3 between 3 and 6 February 2023. The vertical axis is main beam brightness temperature in K. The horizontal axis is the Doppler velocity in the rest frame of the comet. The number of lines averaged is provided for each molecule, either in the 2 mm band (147–155, 163–171, and 176–184 GHz) or 1 mm band (209–272 GHz).</figcaption></figure>
## Appendix B: Tables of production rates

<table><caption>**Table B.1.** Daily production rates in C/2021 A1 in November-December 2021.</caption><thead><tr><th>UT date<br>(mm/dd.dd)</th><th>Molecule</th><th>r<sub>h</sub><br>(au)</th><th>Production rate<br>(×10<sup>25</sup> molecules s<sup>-1</sup>)</th><th>Lines<sup>(a)</sup></th></tr></thead><tbody><tr><td>11/12.19</td><td>HCN</td><td>1.218</td><td>1.7 ± 0.5</td><td>1</td></tr><tr><td>11/14.18</td><td>HCN</td><td>1.187</td><td>3.1 ± 0.4</td><td>1</td></tr><tr><td>11/16.18</td><td>HCN</td><td>1.157</td><td>2.9 ± 0.2</td><td>1</td></tr><tr><td>12/08.46</td><td>HCN</td><td>0.826</td><td>3.4 ± 0.2</td><td>1</td></tr><tr><td>12/09.46</td><td>HCN</td><td>0.812</td><td>3.4 ± 0.1</td><td>1</td></tr><tr><td>12/10.32</td><td>HCN</td><td>0.800</td><td>3.2 ± 0.1</td><td>1</td></tr><tr><td>12/11.58</td><td>HCN</td><td>0.785</td><td>4.6 ± 0.2</td><td>2</td></tr><tr><td>12/12.40</td><td>HCN</td><td>0.774</td><td>3.1 ± 0.1</td><td>1</td></tr><tr><td>12/13.40</td><td>HCN</td><td>0.762</td><td>3.5 ± 0.1</td><td>1</td></tr><tr><td>12/09.50</td><td>HNC</td><td>0.812</td><td>0.28 ± 0.05</td><td>1</td></tr><tr><td>12/10.37</td><td>HNC</td><td>0.800</td><td>0.24 ± 0.05</td><td>1</td></tr><tr><td>12/11.62</td><td>HNC</td><td>0.784</td><td>0.51 ± 0.12</td><td>1</td></tr><tr><td>12/12.40</td><td>HNC</td><td>0.774</td><td>0.22 ± 0.08</td><td>1</td></tr><tr><td>12/13.47</td><td>CH<sub>3</sub>CN</td><td>0.761</td><td>0.59 ± 0.06</td><td>7+(2)</td></tr><tr><td>12/08.55</td><td>H<sub>2</sub>S</td><td>0.824</td><td>3.8 ± 1.5</td><td>(1)</td></tr><tr><td>12/10.53</td><td>H<sub>2</sub>S</td><td>0.798</td><td>8.1 ± 1.1</td><td>1</td></tr><tr><td>12/12.60</td><td>H<sub>2</sub>S</td><td>0.772</td><td>6.0 ± 4.6</td><td>(1)</td></tr><tr><td>12/13.63</td><td>H<sub>2</sub>S</td><td>0.759</td><td>4.4 ± 1.2</td><td>1+(1)</td></tr><tr><td>12/08.55</td><td>CS</td><td>0.824</td><td>3.6 ± 0.4</td><td>1</td></tr><tr><td>12/10.53</td><td>CS</td><td>0.798</td><td>3.7 ± 0.2</td><td>1</td></tr><tr><td>12/12.48</td><td>CS</td><td>0.773</td><td>2.6 ± 0.1</td><td>1</td></tr><tr><td>12/13.59</td><td>CS</td><td>0.760</td><td>2.9 ± 0.1</td><td>2</td></tr><tr><td>12/08.55</td><td>H<sub>2</sub>CO<sub>d</sub><sup>(b)</sup></td><td>0.824</td><td>9.5 ± 2.8</td><td>1</td></tr><tr><td>12/10.53</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>0.798</td><td>11.8 ± 1.8</td><td>1</td></tr><tr><td>12/11.52</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>0.785</td><td>14.6 ± 2.5</td><td>1</td></tr><tr><td>12/12.60</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>0.772</td><td>10.2 ± 0.6</td><td>2+</td></tr><tr><td>12/13.55</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>0.760</td><td>8.6 ± 1.1</td><td>2+(2)</td></tr><tr><td>12/08.54</td><td>CH<sub>3</sub>OH</td><td>0.824</td><td>35 ± 5</td><td>4+(8)</td></tr><tr><td>12/09.50</td><td>CH<sub>3</sub>OH</td><td>0.812</td><td>37 ± 7</td><td>1</td></tr><tr><td>12/10.51</td><td>CH<sub>3</sub>OH</td><td>0.798</td><td>44 ± 3</td><td>11</td></tr><tr><td>12/11.60</td><td>CH<sub>3</sub>OH</td><td>0.784</td><td>43 ± 7</td><td>3+(16)</td></tr><tr><td>12/12.49</td><td>CH<sub>3</sub>OH</td><td>0.774</td><td>38 ± 1</td><td>13+(6)</td></tr><tr><td>12/13.53</td><td>CH<sub>3</sub>OH</td><td>0.760</td><td>28 ± 2</td><td>16+(14)</td></tr></tbody></table>

**Notes.** Subscript “d” has been added to the molecules for which a daughter Haser density profile is assumed with the parent scale length provided below.

(a) Number of lines used for the determination of Q, in parentheses the number of lines that are not individually detected.

(b) Where we assume that H₂CO is produced in the coma with a Haser parent scale length of 1500 km in December.

<table><caption>**Table B.2.** Daily production rates in C/2022 E3 in February 2023.</caption><thead><tr><th>UT date<br>(mm/dd.dd)</th><th>Molecule</th><th>r<sub>h</sub><br>(au)</th><th>Production rate<br>(×10<sup>25</sup> molecules s<sup>-1</sup>)</th><th>Lines<sup>(a)</sup></th></tr></thead><tbody><tr><td>02/03.73</td><td>HCN</td><td>1.168</td><td>4.1 ± 0.1</td><td>1+</td></tr><tr><td>02/04.67</td><td>HCN</td><td>1.173</td><td>6.3 ± 0.1</td><td>1+</td></tr><tr><td>02/05.03</td><td>HCN</td><td>1.175</td><td>4.9 ± 0.3</td><td>1+</td></tr><tr><td>02/05.67</td><td>HCN</td><td>1.178</td><td>5.3 ± 0.1</td><td>1+</td></tr><tr><td>02/06.03</td><td>HCN</td><td>1.180</td><td>4.5 ± 0.1</td><td>1+</td></tr><tr><td>02/06.67</td><td>HCN</td><td>1.183</td><td>4.7 ± 0.1</td><td>1+</td></tr><tr><td>02/03.73</td><td>HNC<sub>d</sub><sup>(b)</sup></td><td>1.168</td><td>0.12 ± 0.05</td><td>1</td></tr><tr><td>02/04.67</td><td>HNC<sub>d</sub></td><td>1.173</td><td>0.34 ± 0.11</td><td>1</td></tr><tr><td>02/05.67</td><td>HNC<sub>d</sub></td><td>1.178</td><td>0.25 ± 0.15</td><td>1</td></tr><tr><td>02/06.03</td><td>HNC<sub>d</sub></td><td>1.180</td><td>0.17 ± 0.12</td><td>1</td></tr><tr><td>02/04.91</td><td>CH<sub>3</sub>CN</td><td>1.174</td><td>0.87 ± 0.05</td><td>12</td></tr><tr><td>02/05.76</td><td>CH<sub>3</sub>CN</td><td>1.178</td><td>0.84 ± 0.06</td><td>17</td></tr><tr><td>02/06.90</td><td>CH<sub>3</sub>CN</td><td>1.184</td><td>0.62 ± 0.17</td><td>8+(4)</td></tr><tr><td>02/04.91</td><td>H<sub>2</sub>S</td><td>1.174</td><td>9.4 ± 0.4</td><td>1+</td></tr><tr><td>02/05.02</td><td>H<sub>2</sub>S</td><td>1.175</td><td>11.6 ± 1.1</td><td>2+</td></tr><tr><td>02/05.90</td><td>H<sub>2</sub>S</td><td>1.179</td><td>7.5 ± 0.6</td><td>2+</td></tr><tr><td>02/07.00</td><td>H<sub>2</sub>S</td><td>1.185</td><td>8.9 ± 1.7</td><td>2+</td></tr><tr><td>02/04.84</td><td>CS<sub>d</sub><sup>(c)</sup></td><td>1.174</td><td>3.8 ± 0.1</td><td>2+</td></tr><tr><td>02/05.86</td><td>CS<sub>d</sub></td><td>1.179</td><td>3.7 ± 0.3</td><td>1+</td></tr><tr><td>02/06.90</td><td>CS<sub>d</sub></td><td>1.184</td><td>3.9 ± 0.3</td><td>2+</td></tr><tr><td>02/04.91</td><td>H<sub>2</sub>CO<sub>d</sub><sup>(d)</sup></td><td>1.174</td><td>12.0 ± 0.9</td><td>1+</td></tr><tr><td>02/05.00</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>1.174</td><td>13.3 ± 0.5</td><td>2+</td></tr><tr><td>02/05.76</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>1.178</td><td>13.5 ± 1.1</td><td>2+</td></tr><tr><td>02/05.91</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>1.179</td><td>10.2 ± 0.4</td><td>3+</td></tr><tr><td>02/06.97</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>1.185</td><td>8.3 ± 0.6</td><td>3+</td></tr><tr><td>02/03.73</td><td>CH<sub>3</sub>OH</td><td>1.168</td><td>76 ± 2</td><td>27</td></tr><tr><td>02/04.67</td><td>CH<sub>3</sub>OH</td><td>1.173</td><td>115 ± 4</td><td>27</td></tr><tr><td>02/04.91</td><td>CH<sub>3</sub>OH</td><td>1.174</td><td>87 ± 1</td><td>43</td></tr><tr><td>02/05.76</td><td>CH<sub>3</sub>OH</td><td>1.178</td><td>87 ± 3</td><td>32</td></tr><tr><td>02/06.03</td><td>CH<sub>3</sub>OH</td><td>1.180</td><td>75 ± 6</td><td>15</td></tr><tr><td>02/06.76</td><td>CH<sub>3</sub>OH</td><td>1.184</td><td>80 ± 2</td><td>41</td></tr><tr><td>02/03.73</td><td>NH<sub>2</sub>CHO</td><td>1.168</td><td>1.5 ± 0.4</td><td>1</td></tr><tr><td>02/04.85</td><td>NH<sub>2</sub>CHO</td><td>1.174</td><td>1.0 ± 0.1</td><td>14</td></tr><tr><td>02/05.85</td><td>NH<sub>2</sub>CHO</td><td>1.179</td><td>0.9 ± 0.1</td><td>12</td></tr><tr><td>02/06.91</td><td>NH<sub>2</sub>CHO</td><td>1.184</td><td>1.2 ± 0.3</td><td>12</td></tr><tr><td>02/04.81</td><td>HNCO</td><td>1.174</td><td>2.0 ± 0.3</td><td>3</td></tr><tr><td>02/05.81</td><td>HNCO</td><td>1.178</td><td>1.4 ± 0.6</td><td>2</td></tr><tr><td>02/06.82</td><td>HNCO</td><td>1.184</td><td>1.8 ± 0.9</td><td>3</td></tr></tbody></table>

**Notes.** Subscript “d” has been added to the molecules for which a daughter Haser density profile is assumed with the parent scale length provided below.

(a) Number of lines used for the determination of Q, in parentheses the number of lines that are not individually detected. A “+” indicates that offset data have been taken into account.

(b) Where we assume that HNC is produced in the coma with a Haser parent scale length 2000 km. (c) Where we assume that CS is produced in the coma with a Haser parent scale length 2000 km. (d) Where we assume that H₂CO is produced in the coma with a Haser parent scale length 1700 km.
<table><caption>**Table B.3.** Production rates in C/2021 A1 in November-December 2021 (weekly average).</caption><thead><tr><th>UT date (mm/dd.d)</th><th>Molecule</th><th>r<sub>h</sub> (au)</th><th>Production rate (×10<sup>25</sup> molecules s<sup>-1</sup>)</th><th>Lines<sup>(a)</sup></th></tr></thead><tbody><tr><td colspan="5" style="text-align:center; font-weight:bold;">13-16 November average</td></tr><tr><td>11/15.8</td><td>HCN</td><td>1.17</td><td>3.0 ± 0.1</td><td>1</td></tr><tr><td>11/15.8</td><td>HNC<sub>d</sub><sup>(b)</sup></td><td>1.17</td><td>&lt;0.44</td><td>1</td></tr><tr><td>11/14.8</td><td>CH<sub>3</sub>CN</td><td>1.18</td><td>0.9 ± 0.3</td><td>(8)</td></tr><tr><td>11/14.8</td><td>H<sub>2</sub>S</td><td>1.18</td><td>&lt;11.5</td><td>(1)</td></tr><tr><td>11/14.8</td><td>CS</td><td>1.18</td><td>2.2 ± 0.5</td><td>1</td></tr><tr><td>11/14.8</td><td>H<sub>2</sub>CO<sub>d</sub></td><td>1.18</td><td>7.0 ± 3.7</td><td>(1)</td></tr><tr><td>11/15.3</td><td>CH<sub>3</sub>OH</td><td>1.17</td><td>39 ± 9</td><td>2+(23)</td></tr><tr><td colspan="5" style="text-align:center; font-weight:bold;">08-13 December average</td></tr><tr><td>12/10.8</td><td>HCN</td><td>0.79</td><td>3.4 ± 0.1</td><td>2</td></tr><tr><td>12/12.8</td><td>H<sup>13</sup>CN</td><td>0.77</td><td>&lt;0.072</td><td>1</td></tr><tr><td>12/12.8</td><td>HC<sup>15</sup>N</td><td>0.77</td><td>&lt;0.047</td><td>1</td></tr><tr><td>12/10.8</td><td>HNC<sub>d</sub><sup>(b)</sup></td><td>0.79</td><td>0.24 ± 0.02</td><td>1</td></tr><tr><td>12/12.5</td><td>CH<sub>3</sub>CN</td><td>0.77</td><td>0.62 ± 0.03</td><td>19+(6)</td></tr><tr><td>12/12.4</td><td>HC<sub>3</sub>N</td><td>0.78</td><td>0.14 ± 0.05</td><td>(8)</td></tr><tr><td>12/12.4</td><td>NH<sub>2</sub>CHO</td><td>0.78</td><td>0.91 ± 0.08</td><td>8+(27)</td></tr><tr><td>12/12.5</td><td>HNCO</td><td>0.78</td><td>2.95 ± 0.30</td><td>6+(6)</td></tr><tr><td>12/11.2</td><td>H<sub>2</sub>S</td><td>0.79</td><td>5.9 ± 0.7</td><td>1+(1)</td></tr><tr><td>12/12.0</td><td>CS</td><td>0.78</td><td>2.9 ± 0.1</td><td>2</td></tr><tr><td></td><td>CS<sub>d</sub><sup>(c)</sup></td><td>0.78</td><td>4.1 ± 0.1</td><td>2</td></tr><tr><td>12/13.1</td><td>C<sup>34</sup>S</td><td>0.77</td><td>0.15 ± 0.05</td><td>1</td></tr><tr><td>12/12.4</td><td>SO</td><td>0.78</td><td>&lt;0.6</td><td>(5)</td></tr><tr><td></td><td>SO<sub>d</sub><sup>(d)</sup></td><td>0.78</td><td>&lt;1.5</td><td>(5)</td></tr><tr><td>12/12.4</td><td>SO<sub>2</sub></td><td>0.78</td><td>&lt;1.6</td><td>(20)</td></tr><tr><td>12/12.4</td><td>OCS</td><td>0.78</td><td>4.3 ± 1.1</td><td>(6)</td></tr><tr><td>12/12.4</td><td>H<sub>2</sub>CS</td><td>0.78</td><td>&lt;0.8</td><td>(5)</td></tr><tr><td>12/12.0</td><td>H<sub>2</sub>CO</td><td>0.78</td><td>5.0 ± 0.2</td><td>4+(2)</td></tr><tr><td></td><td>H<sub>2</sub>CO<sub>d</sub><sup>(e)</sup></td><td>0.78</td><td>10.2 ± 0.5</td><td>4+(2)</td></tr><tr><td>12/12.5</td><td>CO</td><td>0.77</td><td>&lt;78</td><td>(1)</td></tr><tr><td>12/12.4</td><td>HCOOH</td><td>0.78</td><td>7.6 ± 0.9</td><td>(28)</td></tr><tr><td>12/11.6</td><td>CH<sub>3</sub>OH</td><td>0.78</td><td>35 ± 1</td><td>42+(14)</td></tr><tr><td>12/12.4</td><td>CH<sub>3</sub>CHO</td><td>0.78</td><td>1.5 ± 0.3</td><td>(53)</td></tr><tr><td>12/12.4</td><td>(CH<sub>2</sub>OH)<sub>2</sub></td><td>0.78</td><td>5.3 ± 1.3</td><td>(39)</td></tr><tr><td>12/12.0</td><td>HCOOCH<sub>3</sub></td><td>0.78</td><td>&lt;5.8</td><td>(6)</td></tr><tr><td>12/12.4</td><td>CH<sub>2</sub>OHCHO</td><td>0.78</td><td>2.0 ± 0.5</td><td>1+(10)</td></tr><tr><td>12/12.4</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>0.78</td><td>7.3 ± 1.5</td><td>(84)</td></tr><tr><td>12/12.4</td><td>CH<sub>2</sub>CO</td><td>0.78</td><td>1.3 ± 0.6</td><td>(10)</td></tr><tr><td>12/11.6</td><td>CH<sub>3</sub>SH</td><td>0.79</td><td>&lt;3.2</td><td>(32)</td></tr><tr><td>12/10.6</td><td>PH<sub>3</sub></td><td>0.80</td><td>&lt;6.4</td><td>(1)</td></tr><tr><td>12/11.6</td><td>c-C<sub>3</sub>H<sub>2</sub></td><td>0.80</td><td>&lt;0.9</td><td>(11)</td></tr><tr><td>12/11.6</td><td>l-C<sub>3</sub>H<sub>2</sub></td><td>0.80</td><td>&lt;0.2</td><td>(13)</td></tr></tbody></table>

**Notes.** Subscript “d” has been added to the molecules for which a daughter Haser density profile is assumed with the parent scale length provided below.

(a) Number of lines used for the determination of Q, in parentheses the number of lines that are not individually detected.

(b) Where we assume that HNC is produced in the coma with a Haser parent scale length of 1000 km and 500 km in December (Cordiner et al. 2017).

(c) Assuming a parent scale length of 1000 km (about 4× longer than the photo-dissociation scale length of CS₂, assumed parent otherwise).

(d) Assuming a parent scale length of 1700 km (SO₂).

(e) Where we assume that H₂CO is produced in the coma with a Haser parent scale length of 2800 km in November and 1500 km in December.

<table><caption>**Table B.4.** Production rates in C/2022 E3 (ZTF) in February 2023.</caption><thead><tr><th>UT date (mm/dd.d)</th><th>Molecule</th><th>r<sub>h</sub> (au)</th><th>Production rate (×10<sup>25</sup> molecules s<sup>-1</sup>)</th><th>Lines<sup>(a)</sup></th></tr></thead><tbody><tr><td colspan="5" style="text-align:center; font-weight:bold;">3-7 February 2023 average</td></tr><tr><td>02/04.6</td><td>HCN</td><td>1.17</td><td>4.6 ± 0.1</td><td>2</td></tr><tr><td>02/05.1</td><td>H<sup>13</sup>CN</td><td>1.18</td><td>&lt;0.065</td><td>(1)</td></tr><tr><td>02/05.1</td><td>HC<sup>15</sup>N</td><td>1.18</td><td>&lt;0.044</td><td>(1)</td></tr><tr><td>02/05.9</td><td>DCN</td><td>1.18</td><td>&lt;0.084</td><td>(1)</td></tr><tr><td>02/04.6</td><td>HNC<sub>d</sub><sup>(b)</sup></td><td>1.17</td><td>0.16 ± 0.03</td><td>1</td></tr><tr><td>02/05.1</td><td>CH<sub>3</sub>CN</td><td>1.18</td><td>0.84 ± 0.03</td><td>20+(4)</td></tr><tr><td>02/05.2</td><td>HC<sub>3</sub>N</td><td>1.18</td><td>&lt;0.11</td><td>(8)</td></tr><tr><td>02/05.3</td><td>NH<sub>2</sub>CHO</td><td>1.18</td><td>0.94 ± 0.05</td><td>14+(25)</td></tr><tr><td>02/05.3</td><td>HNCO</td><td>1.18</td><td>2.08 ± 0.21</td><td>3+(9)</td></tr><tr><td>02/05.5</td><td>H<sub>2</sub>S</td><td>1.18</td><td>9.0 ± 0.2</td><td>2</td></tr><tr><td>02/05.2</td><td>H<sub>2</sub><sup>4</sup>S</td><td>1.18</td><td>&lt;0.7</td><td>(1)</td></tr><tr><td>02/05.2</td><td>CS</td><td>1.18</td><td>2.4 ± 0.1</td><td>2</td></tr><tr><td></td><td>CS<sub>d</sub><sup>(c)</sup></td><td>1.18</td><td>3.7 ± 0.1</td><td>2</td></tr><tr><td>02/05.3</td><td>C<sup>34</sup>S<sub>d</sub></td><td>1.18</td><td>0.03 ± 0.01</td><td>1</td></tr><tr><td>02/05.3</td><td>SO</td><td>1.18</td><td>&lt;0.5</td><td>(5)</td></tr><tr><td></td><td>SO<sub>d</sub><sup>(d)</sup></td><td>1.18</td><td>&lt;1.5</td><td>(5)</td></tr><tr><td>02/05.3</td><td>SO<sub>2</sub></td><td>1.18</td><td>&lt;1.5</td><td>(16)</td></tr><tr><td>02/05.3</td><td>OCS</td><td>1.18</td><td>3.4 ± 0.8</td><td>(6)</td></tr><tr><td>02/05.1</td><td>H<sub>2</sub>CS</td><td>1.18</td><td>&lt;0.9</td><td>(7)</td></tr><tr><td>02/05.2</td><td>H<sub>2</sub>CO</td><td>1.18</td><td>6.1 ± 0.1</td><td>5+(1)</td></tr><tr><td></td><td>H<sub>2</sub>CO<sub>d</sub><sup>(e)</sup></td><td>1.18</td><td>10.9 ± 0.2</td><td>5+(1)</td></tr><tr><td>02/04.9</td><td>HDCO<sub>d</sub></td><td>1.17</td><td>&lt;0.69</td><td>(5)</td></tr><tr><td>02/05.9</td><td>CO</td><td>1.18</td><td>35 ± 16</td><td>1</td></tr><tr><td>02/05.3</td><td>HCOOH</td><td>1.18</td><td>9.0 ± 0.7</td><td>13+(17)</td></tr><tr><td>02/05.1</td><td>CH<sub>3</sub>OH</td><td>1.18</td><td>87.8 ± 0.8</td><td>62+(12)</td></tr><tr><td>02/05.3</td><td>CH<sub>3</sub>CHO</td><td>1.18</td><td>3.4 ± 0.2</td><td>13+(66)</td></tr><tr><td>02/05.3</td><td>(CH<sub>2</sub>OH)<sub>2</sub></td><td>1.18</td><td>5.5 ± 0.7</td><td>6+(108)</td></tr><tr><td>02/05.3</td><td>HCOOCH<sub>3</sub></td><td>1.18</td><td>&lt;4.8</td><td>(6)</td></tr><tr><td>02/05.3</td><td>CH<sub>2</sub>OHCHO</td><td>1.18</td><td>&lt;1.2</td><td>(13)</td></tr><tr><td>02/05.3</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>1.18</td><td>8.7 ± 1.3</td><td>4+(103)</td></tr><tr><td>02/05.3</td><td>CH<sub>2</sub>CO</td><td>1.18</td><td>&lt;1.5</td><td>(7)</td></tr><tr><td>02/05.3</td><td>CH<sub>3</sub>SH</td><td>1.18</td><td>&lt;3.1</td><td>(20)</td></tr><tr><td>02/05.3</td><td>CH<sub>3</sub>NH<sub>2</sub></td><td>1.18</td><td>&lt;4.8</td><td>(21)</td></tr><tr><td>02/05.3</td><td>CH<sub>3</sub>COCH<sub>3</sub></td><td>1.18</td><td>&lt;1.6</td><td>(29)</td></tr><tr><td>02/04.6</td><td>PH<sub>3</sub></td><td>1.17</td><td>&lt;4.4</td><td>(1)</td></tr><tr><td>02/05.3</td><td>c-C<sub>3</sub>H<sub>2</sub></td><td>1.18</td><td>&lt;0.6</td><td>(23)</td></tr><tr><td>02/05.3</td><td>l-C<sub>3</sub>H<sub>2</sub></td><td>1.18</td><td>&lt;0.1</td><td>(13)</td></tr></tbody></table>

**Notes.** Subscript “d” has been added to the molecules for which a daughter Haser density profile is assumed with the parent scale length provided below.

(a) Number of lines used for the determination of Q, in parentheses the number of lines that are not individually detected.

(b) Where we assume that HNC is produced in the coma with a Haser parent scale length of 2000 km (Cordiner et al. 2017).

(c) Assuming a parent scale length of 2000 km (about 4× longer than the photo-dissociation scale length of CS₂, assumed parent otherwise).

(d) Assuming a parent scale length of 3800 km (SO₂).

(e) Where we assume that H₂CO is produced in the coma with a Haser parent scale length of 1700 km.
<table><caption>**Table B.5.** Production from a distributed source based on offset pointings on comet C/2021 A1 (Leonard).</caption><thead><tr><th>UT (mm/dd.d)</th><th>Molecule</th><th>Freq. range (GHz)</th><th>lines <sup>(a)</sup></th><th>off. <sup>(b)</sup></th><th>Line intensity (K km s<sup>-1</sup>)</th><th>L<sub>p</sub> <sup>(c)</sup> (km)</th><th>Q<sub>p</sub> <sup>(d)</sup> (×10<sup>25</sup> molecules s<sup>-1</sup>)</th><th>χ<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="4">12/12.40</td><td rowspan="4">HCN</td><td rowspan="4">265.886</td><td rowspan="4">1</td><td>1.1</td><td>1.714 ± 0.050</td><td>0</td><td>2.9 ± 0.1</td><td>15.2</td></tr><tr><td>4.8</td><td>1.274 ± 0.091</td><td rowspan="3">560<sup>+330</sup><sub>-220</sub></td><td rowspan="3">3.9 ± 0.1</td><td rowspan="3">7.5</td></tr><tr><td>8.4</td><td>0.977 ± 0.092</td></tr><tr><td>11.0</td><td>1.113 ± 0.125</td></tr><tr><td rowspan="3">12/13.40</td><td rowspan="3">HCN</td><td rowspan="3">265.886</td><td rowspan="3">1</td><td>0.4</td><td>2.022 ± 0.045</td><td>0</td><td>3.4 ± 0.1</td><td>4.5</td></tr><tr><td>10.0</td><td>0.939 ± 0.072</td><td rowspan="2">240 ± 100</td><td rowspan="2">3.9 ± 0.1</td><td rowspan="2">0.1</td></tr><tr><td>14.1</td><td>0.686 ± 0.067</td></tr><tr><td>12/10.6</td><td>HNC</td><td>271.981</td><td>1</td><td>1.7</td><td>0.121 ± 0.012</td><td>0 ± 700</td><td>0.19 ± 0.02</td><td>2.7</td></tr><tr><td>12/11.0</td><td></td><td></td><td></td><td>4.7</td><td>0.011 ± 0.049</td><td>1000</td><td>0.30 ± 0.03</td><td>4.1</td></tr><tr><td>12/11.0</td><td></td><td></td><td></td><td>10.3</td><td>0.050 ± 0.027</td><td></td><td></td><td></td></tr><tr><td rowspan="4">12/12.5</td><td rowspan="4">CS</td><td rowspan="4">244.936</td><td rowspan="4">1</td><td>1.8</td><td>0.344 ± 0.013</td><td>0</td><td>2.2 ± 0.1</td><td>3.7</td></tr><tr><td>3.8</td><td>0.274 ± 0.079</td><td>260</td><td>2.6 ± 0.1</td><td>2.2</td></tr><tr><td>8.2</td><td>0.242 ± 0.038</td><td>500<sup>+580</sup><sub>-390</sub></td><td>3.0 ± 0.1</td><td>1.9</td></tr><tr><td>11.8</td><td>0.105 ± 0.043</td><td>1000</td><td>3.8 ± 0.1</td><td>2.7</td></tr><tr><td rowspan="4">12/13.6</td><td rowspan="4">CS</td><td rowspan="4">244.936</td><td rowspan="4">1</td><td>1.1</td><td>0.357 ± 0.020</td><td>0</td><td>2.5 ± 0.1</td><td>27.6</td></tr><tr><td>8.7</td><td>0.170 ± 0.079</td><td>260</td><td>2.9 ± 0.1</td><td>20.0</td></tr><tr><td>10.4</td><td>0.271 ± 0.039</td><td>1000</td><td>4.2 ± 0.2</td><td>8.3</td></tr><tr><td>146.969</td><td>1</td><td>7.000<sup>+36000</sup><sub>-430</sub></td><td>12.6 ± 0.5</td><td>2.4</td></tr><tr><td rowspan="3">12/12.6</td><td rowspan="3">H<sub>2</sub>CO</td><td rowspan="3">211,226</td><td rowspan="3">2</td><td>0.9</td><td>0.185 ± 0.012</td><td>0</td><td>2.2 ± 0.1</td><td>4.8</td></tr><tr><td>10.2</td><td>0.078 ± 0.022</td><td>1500<sup>+3700</sup><sub>-970</sub></td><td>10.2 ± 0.6</td><td>1.4</td></tr><tr><td>14.7</td><td>0.098 ± 0.043</td><td>4000</td><td>18.8 ± 1.2</td><td>2.0</td></tr><tr><td>12/11.1</td><td>HNCO</td><td>153.865</td><td>1</td><td>1.9</td><td>0.040 ± 0.006</td><td>0</td><td>3.2 ± 0.3</td><td>13.5</td></tr><tr><td>12/11.8</td><td></td><td></td><td></td><td>10.3</td><td>0.036 ± 0.021</td><td>2000</td><td>7.1 ± 0.6</td><td>4.1</td></tr><tr><td>12/12.8</td><td></td><td>219–264</td><td>3</td><td>1.4</td><td>0.050 ± 0.005</td><td>10000</td><td>23.0 ± 1.9</td><td>1.8</td></tr><tr><td>12/12.9</td><td></td><td>242,264</td><td>2</td><td>9.9</td><td>0.030 ± 0.014</td><td></td><td></td><td></td></tr><tr><td rowspan="2">12/12.7</td><td rowspan="2">HCOOH</td><td rowspan="2">215–268</td><td>13<sup>(e)</sup></td><td>1.5</td><td>0.0166 ± 0.0024</td><td>0</td><td>7.3 ± 0.9</td><td>4.9</td></tr><tr><td>10</td><td>10.1</td><td>0.0168 ± 0.0074</td><td>2000</td><td>19.4 ± 2.5</td><td>2.7</td></tr><tr><td>12/12.2</td><td>NH<sub>2</sub>CHO</td><td>211–267</td><td>11</td><td>1.4</td><td>0.0321 ± 0.0031</td><td>0 ± 400</td><td>10.8 ± 1.1</td><td></td></tr><tr><td>12/12.7</td><td></td><td>K<sub>a</sub> < 3</td><td>11</td><td>10.2</td><td>0.0081 ± 0.0079</td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** (a) Number of lines used or averaged. (b) Mean pointing offset. (c) Parent scale length in km. The value with uncertainty (or lower limit) is the one obtained from the $\chi^2$ minimisation. (d) Production rates from nuclear or distributed source with given parent scale length.

(e) The $\chi^2$ fitting was computed using all individual lines intensities separately, not their average. (f) This $\chi^2$ for HCN(3-2) includes an additional 10% calibration uncertainty. Without, $\chi^2$ is about two orders of magnitude but leads to the same $L_p$.
<table><caption>**Table B.6.** Production from a distributed source based on offset pointings on comet C/2022 E3 (ZTF).</caption><thead><tr><th>UT<br>(mm/dd.d)</th><th>Molecule</th><th>Freq. range<br>(GHz)</th><th>lines <sup>(a)</sup></th><th>off.<sup>(b)</sup><br>(")</th><th>Line intensity<br>(K km s<sup>-1</sup>)</th><th>L<sub>p</sub> <sup>(c)</sup><br>(km)</th><th>Q<sub>p</sub> <sup>(d)</sup><br>(×10<sup>25</sup> molecules s<sup>-1</sup>)</th><th>χ<sup>2</sup></th></tr></thead><tbody><tr><td>02/04.55</td><td>HCN</td><td>265.886</td><td>1</td><td>1.6<br>5.4<br>9.8<br>14.0</td><td>2.292 ± 0.011<br>1.716 ± 0.029<br>1.106 ± 0.017<br>0.953 ± 0.122</td><td>0<br>600 ± 500</td><td>4.6 ± 0.3<br>5.4 ± 0.3</td><td>2.2 <sup>(e)</sup><br>0.9</td></tr><tr><td>02/05.13</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/04.47</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/04.38</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.03</td><td>HCN</td><td>177.261</td><td>1</td><td>0.9<br>10.1</td><td>1.122 ± 0.061<br>0.731 ± 0.079</td><td>0 ± 800</td><td>4.9 ± 0.2</td><td>0.01</td></tr><tr><td>02/04.55</td><td>HNC</td><td>271.981</td><td>1</td><td>1.6<br>5.4<br>9.8<br>14.0</td><td>0.038 ± 0.009<br>0.077 ± 0.026<br>0.034 ± 0.017<br>0.084 ± 0.033</td><td>0<br>2000<br>>2800</td><td>0.74 ± 0.14<br>1.59 ± 0.29<br>>1.92 ± 0.34</td><td>9.2<br>6.3<br>6.0</td></tr><tr><td>02/05.13</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/04.47</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/04.38</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.23</td><td>CS</td><td>146.969</td><td>1</td><td>2.1<br>8.8<br>13.0</td><td>0.096 ± 0.004<br>0.087 ± 0.010<br>0.059 ± 0.009</td><td>0<br>550<br>1400</td><td>1.95 ± 0.05<br>2.44 ± 0.06<br>3.23 ± 0.08</td><td>30.1<br>11.8<br>3.3</td></tr><tr><td>02/05.41</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.08</td><td>CS</td><td>244.936</td><td>1</td><td>1.8<br>5.1<br>8.4</td><td>0.247 ± 0.008<br>0.196 ± 0.018<br>0.166 ± 0.034</td><td>1900<sup>+900</sup><sub>-470</sub><br>2000</td><td>3.65 ± 0.08<br>3.74 ± 0.09</td><td>2.6<br>2.7</td></tr><tr><td>02/05.18</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.00</td><td>H<sub>2</sub>CO</td><td>211.226</td><td>2</td><td>1.2<br>10.1<br>14.2<br>20.0</td><td>0.282 ± 0.012<br>0.138 ± 0.016<br>0.078 ± 0.019<br>0.080 ± 0.027</td><td>0<br>1000<sup>+500</sup><sub>-330</sub><br>1700</td><td>7.29 ± 0.28<br>10.89 ± 0.41<br>13.25 ± 0.51</td><td>12.4<br>1.7<br>3.4</td></tr><tr><td>02/05.95</td><td>H<sub>2</sub>CO</td><td>211.226</td><td>2</td><td>1.3<br>10.1<br>20.1<br>28.3</td><td>0.209 ± 0.008<br>0.107 ± 0.018<br>0.053 ± 0.020<br>0.032 ± 0.019</td><td>0<br>1700<sup>+1700</sup><sub>-640</sub></td><td>5.59 ± 0.21<br>10.37 ± 0.39</td><td>8.5<br>0.5</td></tr><tr><td>02/06.94</td><td>H<sub>2</sub>CO</td><td>211.226</td><td>2</td><td>1.0<br>4.9<br>10.0</td><td>0.168 ± 0.013<br>0.117 ± 0.026<br>0.100 ± 0.018</td><td>0<br>1700<sup>+3700</sup><sub>-1100</sub></td><td>4.65 ± 0.33<br>8.21 ± 0.57</td><td>5.0<br>1.5</td></tr><tr><td>02/05.2</td><td>H<sub>2</sub>CO</td><td>150.498</td><td>1</td><td>2.1<br>8.8<br>13.0</td><td>0.075 ± 0.004<br>0.045 ± 0.011<br>0.037 ± 0.012</td><td>0 ± 2200<br>1700</td><td>6.97 ± 0.36<br>11.02 ± 0.57</td><td>0.5<br>1.1</td></tr><tr><td>02/05.25</td><td>HNCO</td><td>220-264</td><td>3</td><td>2.1<br>5.1<br>9.1</td><td>0.036 ± 0.005<br>0.048 ± 0.011<br>0.043 ± 0.018</td><td>0<br>>2000<br>6000</td><td>2.06 ± 0.25<br>>3.97 ± 0.46<br>7.95 ± 0.92</td><td>6.2<br>3.4<br>2.7</td></tr><tr><td>02/04.92</td><td></td><td>K<sub>a</sub> = 0</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.42</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.25</td><td>HNCO</td><td>219-265</td><td>6</td><td>2.0<br>5.2<br>9.3</td><td>0.013 ± 0.003<br>0.021 ± 0.007<br>0.023 ± 0.009</td><td>0<br>>1000<br>6000</td><td>2.99 ± 0.57<br>>6.44 ± 1.20<br>24.9 ± 4.6</td><td>8.4<br>7.3<br>6.5</td></tr><tr><td>02/04.92</td><td></td><td>K<sub>a</sub> = 1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.42</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.34</td><td>HCOOH</td><td>215-268</td><td>15</td><td>1.7<br>5.1<br>9.9</td><td>0.0276 ± 0.0017<br>0.0189 ± 0.0057<br>0.0102 ± 0.0046</td><td>0<br>1100<sup>+3100</sup><sub>-1100</sub></td><td>10.39 ± 0.71<br>17.11 ± 1.03</td><td>0.9<br>0.1</td></tr><tr><td>02/05.25</td><td></td><td>K<sub>a</sub> = 0 - 2</td><td>15</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.31</td><td></td><td></td><td>15</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.23</td><td>NH<sub>2</sub>CHO</td><td>147-169</td><td>7</td><td>2.1<br>8.9<br>13.0</td><td>0.0120 ± 0.0014<br>0.0040 ± 0.0039<br>0.0068 ± 0.0036</td><td>0 ± 2500<br>1300</td><td>0.84 ± 0.10<br>1.18 ± 0.14</td><td>1.1<br>1.7</td></tr><tr><td>02/05.23</td><td></td><td>K<sub>a</sub> = 0 - 2</td><td>7</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.41</td><td></td><td>J = 7 - 8</td><td>7</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.44</td><td>NH<sub>2</sub>CHO</td><td>211-267</td><td>15</td><td>1.7<br>5.1<br>9.9</td><td>0.0235 ± 0.0016<br>0.0329 ± 0.0052<br>0.0106 ± 0.0038</td><td>0<br>1300<sup>+2500</sup><sub>-900</sub></td><td>0.93 ± 0.06<br>1.50 ± 0.10</td><td>8.6<br>6.6</td></tr><tr><td>02/05.49</td><td></td><td>J = 10 - 13</td><td>15</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/04.83</td><td></td><td></td><td>6</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.43</td><td>NH<sub>2</sub>CHO</td><td>212-256</td><td>15</td><td>1.6<br>5.2<br>10.0</td><td>0.0131 ± 0.0018<br>0.0158 ± 0.0066<br>0.0155 ± 0.0031</td><td>0<br>1300<br>>3000</td><td>1.07 ± 0.13<br>1.98 ± 0.23<br>>3.24 ± 0.37</td><td>13.5<br>8.5<br>6.6</td></tr><tr><td>02/05.26</td><td></td><td>K<sub>a</sub> = 3 - 5</td><td>15</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.49</td><td></td><td>J = 10 - 12</td><td>15</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.94</td><td>CH<sub>3</sub>CHO</td><td>211-232</td><td>18</td><td>1.2<br>4.9<br>10.1</td><td>0.0115 ± 0.0013<br>0.0061 ± 0.0069<br>0.0043 ± 0.0017</td><td>0<br>1600<sup>+5000</sup><sub>-1000</sub></td><td>3.3 ± 0.4<br>5.9 ± 0.6</td><td>8.0<br>6.5</td></tr><tr><td>02/06.94</td><td></td><td>K<sub>a</sub> = 0 - 3</td><td>18</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.85</td><td></td><td>J=11,12</td><td>18</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.00</td><td></td><td></td><td>18</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/05.62</td><td></td><td></td><td>18</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/04.55</td><td>CH<sub>3</sub>CHO</td><td>251-271</td><td>16</td><td>1.6<br>5.4<br>9.8</td><td>0.0133 ± 0.0017<br>0.0125 ± 0.0044<br>0.0067 ± 0.0029</td><td>0<br>2000<sup>+4000</sup><sub>-1600</sub></td><td>3.5 ± 0.4<br>7.7 ± 0.9</td><td>1.8<br>0.2</td></tr><tr><td>02/05.13</td><td></td><td>K<sub>a</sub> = 1 - 3</td><td>16</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/04.47</td><td></td><td>J=13,14</td><td>16</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>02/04.38</td><td></td><td></td><td>16</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** (b)Mean pointing offset. (a)Number of lines used or averaged. (c)Parent scale length in km. The value with uncertainty (or lower limit) is the one obtained from the $\chi^2$ minimisation. (d)Production rates from nuclear or distributed source with given parent scale length. (e)This $\chi^2$ for HCN(3-2) includes an additional 10% calibration uncertainty. Without, $\chi^2$ is about two orders of magnitude larger but minimisation leads to the same $L_p$.
## Appendix C: Rotational diagrams

The logarithm of a quantity proportional to the line intensity (line area divided by column density) is plotted against the energy of the upper level for each observed rotational transition. The upper levels of the transition are given at the top. The inverse of the slope of a fitted line provides the rotational temperature (given in pink) of the group of lines.

<figure><img src="image_25.png" /><figcaption>Fig. C.1. Rotational diagram of the CH₃OH lines around 166 GHz observed between 8.6 and 13.7 December 2021 in comet C/2021 A1 (Leonard).</figcaption></figure>

<figure><img src="image_26.png" /><figcaption>Fig. C.2. Rotational diagram of all the NH₂CHO lines observed between 8.6 and 13.7 December 2021 in comet C/2021 A1 (Leonard). Lines with same *J* number and *K*<sub>*ij*</sub>=0 to 2 or 3 to 4 are averaged together. Scales as in Fig. C.1). The values in green are those expected from modelling with *T* = 60 K.</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>Fig. C.3. Rotational diagram of the CH₃CN lines at 147, 165, 220, 239, and 257 GHz observed between 8.4 and 13.7 December 2021 in comet C/2021 A1 (Leonard). Scales as in Fig. C.1). The values in green are those expected from modelling with *T* = 60 K.</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>Fig. C.4. Rotational diagram of the HNCO lines at 1 mm and 2 mm observed between 8.4 and 13.7 December 2021 in comet C/2021 A1 (Leonard). Scales as in Fig. C.1.</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>Fig. C.5. Rotational diagram of the CH₃OH lines around 166 GHz observed between 4.9 and 7.0 February 2023 in comet C/2022 E3 (ZTF). The values in green are those expected from modelling with *T* = 60 K.</figcaption></figure>
<figure><img src="image_30.png" /><figcaption>**Fig. C.6.** Rotational diagram of the CH<sub>3</sub>OH lines at 242 GHz observed between 4.7 and 6.8 February 2023 in comet C/2022 E3 (ZTF). Scales as in Fig. C.5.</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>**Fig. C.9.** Rotational diagram of the CH<sub>3</sub>CN lines at 257 GHz observed between 3.7 and 6.7 February 2023 in comet C/2022 E3 (ZTF). Scales as in Fig. C.5.</figcaption></figure>

<figure><img src="image_32.png" /><figcaption>**Fig. C.7.** Rotational diagram of the CH<sub>3</sub>OH lines around 252 GHz observed between 3.7 and 6.7 February 2023 in comet C/2022 E3 (ZTF). Scales as in Fig. C.5.</figcaption></figure>

<figure><img src="image_33.png" /><figcaption>**Fig. C.10.** Rotational diagram of all the NH<sub>2</sub>CHO lines observed between 3.6 and 7.1 February 2023 in comet C/2022 E3 (ZTF) at 1 mm and 2 mm. Some lines with the same *J* and *K*<sub>*a*</sub> levels (and/or similar energy levels and Einstein A coefficients) have been grouped together (*K*<sub>*c*</sub> given as 0 in the third quantum number of the line). Scales as in Fig. C.5.</figcaption></figure>

<figure><img src="image_34.png" /><figcaption>**Fig. C.8.** Rotational diagram of the CH<sub>3</sub>CN lines at 147 and 165 GHz observed between 4.9 and 7.0 February 2023 in comet C/2022 E3 (ZTF). Scales as in Fig. C.5.</figcaption></figure>

<figure><img src="image_35.png" /><figcaption>**Fig. C.11.** Rotational diagram of all the HNCO lines observed between 3.6 and 7.1 February 2023 in comet C/2022 E3 (ZTF). The lines with the same *J* for *K*<sub>*a*</sub> = 1 levels have been grouped together. Scales as in Fig. C.5.</figcaption></figure>
<figure><img src="image_36.png" /><figcaption>**Fig. C.12.** Rotational diagram of all the HCOOH lines observed between 3.6 and 7.1 February 2023 in comet C/2022 E3 (ZTF). Scales and grouping of lines as in Fig. C.10.</figcaption></figure>

<figure><img src="image_37.png" /><figcaption>**Fig. C.13.** Rotational diagram of all the CH<sub>3</sub>CHO lines observed between 3.6 and 7.1 February 2023 in comet C/2022 E3 (ZTF). Scales and grouping of lines as in Fig. C.10.</figcaption></figure>