pith. sign in

arxiv: 2604.02613 · v2 · submitted 2026-04-03 · 🌌 astro-ph.EP

Monitoring Volatile Evolution in Disrupting Comet D/2021 A1 (Leonard) with NOEMA and APEX

Pith reviewed 2026-05-13 19:04 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords comet volatile evolutionHCNCScometary disruptionmixing ratiosNOEMAAPEXD/2021 A1
0
0 comments X

The pith

Disruption and fragmentation in comet D/2021 A1 drove changes in volatile mixing ratios beyond what solar distance alone predicts.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents multi-epoch observations of comet D/2021 A1 from November to December 2021 using NOEMA and APEX, tracking HCN and CS emissions as the comet moved inward from 1.3 to 0.8 au. CS mixing ratios relative to water increased by a factor of about five over this period, consistent with production from a distributed source rather than the nucleus. HCN showed smaller changes and more scatter when combined with data from other telescopes. The observations captured outbursts and fragmentation in mid-December, during which both species varied in ways that point to physical breakup releasing additional material. This demonstrates that rapid, species-specific changes in comets require repeated monitoring across instruments to separate insolation effects from disruption processes.

Core claim

Comet D/2021 A1 exhibited a clear rise in CS mixing ratios from 0.02% to 0.10% as heliocentric distance decreased, while HCN remained near 0.07% with added variability once cross-facility data were included; line kinematics showed HCN released near the nucleus and CS from a distributed source, with both species displaying changes tied to mid-December outbursts and fragmentation.

What carries the argument

Mixing ratios of HCN and CS relative to H2O, derived from line intensities and spatial constraints, that separate nucleus versus distributed sources and track changes across epochs.

If this is right

  • CS production increases steadily inward because the molecule forms in the coma from a distributed parent.
  • HCN shows no robust single trend with distance once all data are merged, indicating multiple release mechanisms.
  • Outbursts and fragmentation cause short-term spikes in both species that solar insolation models alone cannot explain.
  • Multi-facility, multi-epoch coverage is required to resolve species-dependent responses on timescales of days.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar monitoring of other fragmenting comets could reveal whether buried ices are preferentially released during breakup events.
  • The observed CS rise may imply that distributed-source chemistry becomes more efficient closer to the Sun in active comets.
  • Accounting for disruption effects could improve predictions of volatile delivery to inner planets during cometary passages.

Load-bearing premise

That HCN measurements from different telescopes can be combined directly without calibration offsets large enough to erase the apparent lack of monotonic trend with distance.

What would settle it

A new campaign using a single calibrated instrument across the same distance range that finds HCN mixing ratios either constant or strictly decreasing with no outburst-linked spikes would undermine the claim that disruption drives the observed variability.

Figures

Figures reproduced from arXiv: 2604.02613 by Boncho P. Bonev, Dominique Bockel\'ee-Morvan, J\'er\'emie Boissier, Martin A. Cordiner, Michael A. DiSanti, Nathan X. Roth, Neil Dello Russo, Nicolas Biver, Stefanie N. Milam, Steven B. Charnley, Timothy N. Proudkii.

Figure 1
Figure 1. Figure 1: Spectrally integrated flux maps for NOEMA HCN (J=3-2) detections on November 5 (left) and November 21 (right). In the bottom left of the map, we display the size and orientation of the synthesized beam. The lower right of the map indicates the comet’s illumination, the Sun’s orientation (S), and the dust trail (T), while the upper right panel displays the spectral line extracted from the brightest pixel. F… view at source ↗
Figure 2
Figure 2. Figure 2: The two panels show ∆χ 2 as a function of parent scale length (Lp) for November 5 (left) and 21 (right). The black line connects the discrete ∆χ 2 points to guide the eye. Dotted horizontal lines indicate the ∆χ 2 = 1 threshold (1σ, 68% confidence), while dashed horizontal lines mark the ∆χ 2 = 6.63 threshold (2.6σ, 99% confidence), with χ 2 computed from the difference between the observed interferometric… view at source ↗
Figure 3
Figure 3. Figure 3: (A) HCN (J=3-2) spectrum of A1 on UT 2021 November 5 (black). The top panel shows the ON-OFF spectrum, while the lower two panels display spectra extracted over different baseline ranges (i.e., angular scales). The spectra are shown at a frequency resolution of 125 kHz (velocity resolution 0.15 km s−1 ). The best-fit model is overplotted in red. (B) Real part of the observed visibility amplitude as a funct… view at source ↗
Figure 4
Figure 4. Figure 4: Evolution of the HCN spectral line in comet A1. The black line depicts the extracted spectra, while the red line represents our modeled data. NOEMA spectra shown here are from the ON-OFF observations. in others, and help constrain the processes shaping A1’s activity during its inbound leg. Our stringent 3σ up￾per limit on the CH3OH mixing ratio in A1 (< 1.2%; relative to H2O) is depleted compared to averag… view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of the CS spectral line in comet A1. The black line depicts the extracted spectra, while the red line represents our modeled data. NOEMA spectra shown here are from the ON-OFF observations [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Pre-perihelion time evolution (decreasing rH) of the CS (top) and HCN (bottom) mixing ratios relative to H2O in comet A1. Black open symbols show results from this work. Blue colored, solid symbols show contemporaneous measurements adopted from N. Biver et al. (2024) while the red colored, solid symbol shows measurements from S. Faggi et al. (2023). Shown error bars represent 1σ uncertainties that include … view at source ↗
Figure 7
Figure 7. Figure 7: Pre-perihelion time evolution (decreasing rH) of the CS/HCN mixing ratio in comet A1. Because the ratio does not rely on normalization to H2O, it provides a water-independent comparison of relative abundance variability. Black open symbols show results from this work. Blue colored, solid symbols show contemporaneous measurements adopted from N. Biver et al. (2024). In November, for N. Biver et al. (2024), … view at source ↗
read the original abstract

We report a pre-perihelion survey of volatile emissions from comet D/2021 A1 (Leonard) with the Northern Extended Millimeter Array (NOEMA; UT 2021 Nov. 5, 21, and Dec. 1) and the Atacama Pathfinder Experiment (APEX; UT 2021 Dec. 9-10), spanning heliocentric distances ($r_H$) from 1.3 to 0.80 au. We securely detected HCN and CS and place 3$\sigma$ upper limits on CH$_3$OH, H$_2$CO, and CO abundances. Line kinematics and NOEMA spatial constraints indicate that HCN was released at or near the nucleus (parent scale length $<300$ km), while CS showed higher gas expansion velocities and mixing ratios that increased with decreasing $r_H$ $-$ consistent with production from a distributed source. Across our campaign, CS mixing ratios relative to H$_2$O increased by a factor of $\sim$5, from $0.02 \pm 0.01\%$ at $r_H$ = 1.3 au to $0.10\pm0.02\%$ by $r_H$ = 0.80 au. HCN mixing ratios in our data rose modestly, from $0.04 \pm 0.02\%$ at $r_H$ = 1.3 au to $0.07 \pm 0.02\%$ by $r_H$ = 0.81 au. However, contemporaneous measurements from other facilities placed HCN consistently at a higher absolute level ($\sim\!0.08\%$) with additional variability. Once cross-facility measurements were included, the HCN abundance showed no statistically robust monotonic dependence on $r_H$. Variability in both species during the mid-December outbursts and fragmentation suggests that D/2021 A1's volatile evolution reflected not only solar insolation but also disruption processes, underscoring the value of multi-epoch, multi-instrument monitoring to capture rapid, species-dependent changes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The manuscript reports pre-perihelion millimeter observations of comet D/2021 A1 (Leonard) with NOEMA (Nov 5, 21 and Dec 1 2021) and APEX (Dec 9-10 2021) at heliocentric distances 1.3–0.80 au. Secure detections of HCN and CS are presented together with 3σ upper limits on CH₃OH, H₂CO and CO. Line kinematics and NOEMA spatial constraints indicate HCN release at or near the nucleus (parent scale length <300 km) while CS exhibits higher expansion velocities and increasing mixing ratios with decreasing r_H, consistent with a distributed source. CS/H₂O rises by a factor of ~5 (0.02±0.01% to 0.10±0.02%); the authors’ HCN/H₂O values rise modestly (0.04±0.02% to 0.07±0.02%), but inclusion of contemporaneous external data (~0.08%) removes any monotonic r_H trend and reveals additional variability during mid-December outbursts and fragmentation. The paper concludes that volatile evolution reflects both solar insolation and disruption processes.

Significance. If the cross-facility HCN data can be placed on a common absolute scale, the work provides concrete evidence that fragmentation events drive species-dependent volatile changes on short timescales, beyond the usual r_H dependence. The kinematic and spatial arguments for parent versus distributed sources are a clear strength, and the multi-epoch, multi-instrument approach demonstrates the scientific return from coordinated monitoring of disrupting comets.

major comments (1)
  1. [Abstract] Abstract: the interpretation that mid-December variability in HCN (and CS) indicates disruption-driven evolution beyond insolation rests on combining the authors’ NOEMA/APEX mixing ratios (0.04–0.07%) with external measurements (~0.08%). No section quantifies beam-filling factors, absolute flux calibration offsets, or differences in production-rate modeling assumptions between facilities; if these systematics are comparable to the reported factor-of-2 swings, the physical conclusion is not yet secured.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive review and positive assessment of the manuscript's significance. We address the single major comment below and will incorporate revisions to strengthen the cross-facility comparison.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the interpretation that mid-December variability in HCN (and CS) indicates disruption-driven evolution beyond insolation rests on combining the authors’ NOEMA/APEX mixing ratios (0.04–0.07%) with external measurements (~0.08%). No section quantifies beam-filling factors, absolute flux calibration offsets, or differences in production-rate modeling assumptions between facilities; if these systematics are comparable to the reported factor-of-2 swings, the physical conclusion is not yet secured.

    Authors: We agree that the abstract's interpretation would be strengthened by explicit quantification of cross-facility systematics. In the revised manuscript we will add a new subsection (Section 4.2) that (i) estimates beam-filling factors using the NOEMA synthesized beam and the comet's known angular size at each epoch, (ii) propagates absolute flux calibration uncertainties (typically 10–15% for NOEMA and 20% for APEX) into the mixing ratios, and (iii) compares production-rate modeling assumptions (Haser vs. vectorial models, expansion velocity prescriptions) with those adopted in the external data sets. Preliminary estimates show that the combined uncertainties are ~25–30%, smaller than the observed factor-of-~2 variability in HCN. We will also revise the abstract and conclusions to state that the disruption-driven interpretation is supported once these systematics are accounted for, while noting that contemporaneous multi-facility observations would further reduce residual offsets. These changes directly address the concern and secure the physical conclusion. revision: yes

Circularity Check

0 steps flagged

No significant circularity; claims rest on direct empirical measurements

full rationale

The paper presents observational results from NOEMA and APEX telescope data, reporting detected mixing ratios for HCN and CS, upper limits for other species, and trends with heliocentric distance. These are derived from line detections, kinematics, and spatial constraints without any equations, fitted parameters, or self-citations that reduce a claimed prediction back to the input data by construction. Cross-facility comparisons are included as external data points to contextualize variability, but the paper does not define its own results in terms of those or invoke uniqueness theorems. The central claims about disruption-driven evolution are empirical interpretations of measured abundances, not forced by internal definitions or self-referential loops. This is a standard observational astronomy paper with no load-bearing circular steps.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on standard cometary science assumptions about gas expansion velocities and parent scale lengths used to interpret line profiles and spatial distributions. No new free parameters, ad-hoc axioms, or invented entities are introduced.

axioms (1)
  • domain assumption Standard assumptions about cometary gas expansion velocities and parent scale lengths for interpreting line profiles and spatial constraints.
    Invoked to conclude HCN release at or near the nucleus (scale length <300 km) and CS from a distributed source.

pith-pipeline@v0.9.0 · 5762 in / 1333 out tokens · 72347 ms · 2026-05-13T19:04:32.001805+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages

  1. [1]

    2020, Evidence of ammonium salts in comet 67P as explanation for the nitrogen depletion in cometary comae, NatAs, 4, 533

    Altwegg, K., Balsiger, H., H¨ anni, N., et al. 2020, Evidence of ammonium salts in comet 67P as explanation for the nitrogen depletion in cometary comae, NatAs, 4, 533

  2. [2]

    2022, Emission from hcn and ch3oh in comets-onsala 20-m observations and radiative transfer modelling, Astronomy & Astrophysics, 660, A118

    Bergman, P., Lerner, M., Olofsson, A., et al. 2022, Emission from hcn and ch3oh in comets-onsala 20-m observations and radiative transfer modelling, Astronomy & Astrophysics, 660, A118

  3. [3]

    2011, Molecular investigations of comets C/2002 X5 (Kudo-Fujikawa), C/2002 V1 (NEAT), and C/2006 P1 (McNaught) at small heliocentric distances, Astronomy & Astrophysics, 528, A142

    Biver, N., Bockel´ ee-Morvan, D., Colom, P., et al. 2011, Molecular investigations of comets C/2002 X5 (Kudo-Fujikawa), C/2002 V1 (NEAT), and C/2006 P1 (McNaught) at small heliocentric distances, Astronomy & Astrophysics, 528, A142

  4. [4]

    2022, Observations of comet C/2020 F3 (NEOWISE) with IRAM telescopes, Astronomy & Astrophysics, 668, A171

    Biver, N., Boissier, J., Bockel´ ee-Morvan, D., et al. 2022, Observations of comet C/2020 F3 (NEOWISE) with IRAM telescopes, Astronomy & Astrophysics, 668, A171

  5. [5]

    1999, Spectroscopic Monitoring of Comet C/1996 B2 (Hyakutake) with the JCMT and IRAM Radio

    Biver, N., Bockel´ ee-Morvan, D., Crovisier, J., et al. 1999, Spectroscopic Monitoring of Comet C/1996 B2 (Hyakutake) with the JCMT and IRAM Radio

  6. [6]

    Biver, N., Bockel´ ee-Morvan, D., Handzlik, B., et al. 2024, 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, 690, A271 Bockel´ ee-Morvan, D., & Biver, N. 2017, The composition of cometary ices, PTRSA, 375, 20160252, doi: 10.1098/rsta.201...

  7. [7]

    Weaver, H. A. 2004, in Comets II, ed. H. U. Keller & H. A. Weaver (University of Arizona Press), 391 Bockel´ ee-Morvan, D., Lis, D. C., Wink, J. E., et al. 2000, New molecules found in comet C/1995 O1 (Hale-Bopp). Investigating the link between cometary and interstellar material, A&A, 353, 1101 Bockel´ ee-Morvan, D., Boissier, J., Biver, N., & Crovisier, ...

  8. [8]

    2007, Interferometric imaging of the sulfur-bearing molecules H2S, SO, and CS in comet C/1995 O1 (Hale-Bopp), A&A, 475, 1131

    Boissier, J., Bockel´ ee-Morvan, D., Biver, N., et al. 2007, Interferometric imaging of the sulfur-bearing molecules H2S, SO, and CS in comet C/1995 O1 (Hale-Bopp), A&A, 475, 1131

  9. [9]

    P., Dello Russo, N., DiSanti, M

    Bonev, B. P., Dello Russo, N., DiSanti, M. A., et al. 2021, First Comet Observations with NIRSPEC-2 at Keck: Outgassing Sources of Parent Volatiles and Abundances Based on Alternative Taxonomic Compositional Baselines in 46P/Wirtanen, PSJ, 2, 45

  10. [10]

    D., Witteborn, F

    Bregman, J. D., Witteborn, F. C., Allamandola, L. J., et al. 1987, Airborn and groundbased spectrophotometry of comet P/Halley from 5-13 micrometers, Astronomy and Astrophysics, 187, 616

  11. [11]

    A., Chin, G., et al

    Brooke, T., Weaver, H. A., Chin, G., et al. 2003, Spectroscopy of Comet Hale–Bopp in the infrared, Icarus, 166, 167

  12. [12]

    Ferron, S. 2023, Water production rates from SOHO/SWAN observations of comets C/2020S3 (Erasmus), C/2021 A1 (Leonard) and C/2021 O3 (PanSTARRS), Icarus, 398, 115543, doi: https://doi.org/10.1016/j.icarus.2023.115543

  13. [13]

    A., Remijan, A

    Cordiner, M. A., Remijan, A. J., Boissier, J., et al. 2014, Mapping the Release of Volatiles in the Inner Comae of Comets C/2012 F6 (Lemmon) and C/2012 S1 (ISON) Using the Atacama Large Millimeter/Submillimeter

  14. [14]

    A., Boissier, J., Charnley, S

    Cordiner, M. A., Boissier, J., Charnley, S. B., et al. 2017b, ALMA Mapping of Rapid Gas and Dust Variations in Comet C/2012 S1 (ISON):New Insights into the Origin of Cometary HNC, ApJ, 838, 147

  15. [15]

    A., Coulson, I

    Cordiner, M. A., Coulson, I. M., Garcia-Berrios, E., et al. 2022, A SUBLIME 3D Model for Cometary Coma Emission: The Hypervolatile-rich Comet C/2016 R2 (PanSTARRS), The Astrophysical Journal, 929, 38, doi: 10.3847/1538-4357/ac5893

  16. [16]

    A., Roth, N

    Cordiner, M. A., Roth, N. X., Milam, S. N., et al. 2023, Gas Sources from the Coma and Nucleus of Comet 46P/Wirtanen Observed Using ALMA, https://arxiv.org/abs/2305.04822

  17. [17]

    2008, Distributed Sources in Comets, Space Science Reviews, 138, 179

    Cottin, H., & Fray, N. 2008, Distributed Sources in Comets, Space Science Reviews, 138, 179

  18. [18]

    2021, Central Bureau Electronic Telegrams, 5087, de Val-Borro, M., Milam, S., Cordiner, M., et al

    Crovisier, J., Biver, N., & Bockel´ ee-Morvan, D. 2021, Central Bureau Electronic Telegrams, 5087, de Val-Borro, M., Milam, S., Cordiner, M., et al. 2018, Measuring molecular abundances in comet C/2014 Q2 (Lovejoy) using the APEX telescope, Monthly Notices of the Royal Astronomical Society, 474, 1099 Dello Russo, N., Kawakita, H., Vervack, R. J. J., & Wea...

  19. [19]

    Magee-Sauer, K., & Griep, D. M. 2003, Evidence for a dominant native source of carbon monoxide in Comet C/1996 B2 (Hyakutake), Journal of Geophysical Research: Planets, 108

  20. [20]

    2001, Carbon monoxide production and excitation in Comet C/1995 O1 (Hale-Bopp): Isolation of native and distributed CO sources, Icarus, 153, 361

    Magee-Sauer, K. 2001, Carbon monoxide production and excitation in Comet C/1995 O1 (Hale-Bopp): Isolation of native and distributed CO sources, Icarus, 153, 361

  21. [21]

    A., Mumma, M

    DiSanti, M. A., Mumma, M. J., Russo, N. D., et al. 1999, Identification of two sources of carbon monoxide in Comet Hale–Bopp, Nature, 399, 662

  22. [22]

    A., Bonev, B

    DiSanti, M. A., Bonev, B. P., Gibb, E. L., et al. 2016, En Route to Destruction: The Evolution in Composition of Ices in Comet D/2012 S1 (ISON) between 1.2 and 0.34 AU from the Sun as Revealed at Infrared Wavelengths, ApJ, 820, 20

  23. [23]

    Dubernet, M., & Quintas-S´ anchez, E. 2019, First quantum study of the rotational excitation of HCN by para-H2O: Convergence of quantum results, influence of the potential energy surface, and approximate rate coefficients of interest for cometary atmospheres., Molecular Astrophysics, 16, 100046, doi: https://doi.org/10.1016/j.molap.2019.100046

  24. [24]

    1988, in Exploration of Halley’s Comet (Springer), 481–484

    Eberhardt, P., Krankowsky, D., Schulte, W., et al. 1988, in Exploration of Halley’s Comet (Springer), 481–484

  25. [25]

    J., & Villanueva, G

    Faggi, S., Lippi, M., Mumma, M. J., & Villanueva, G. L. 2023, Strongly Depleted Methanol and Hypervolatiles in Comet C/2021 A1 (Leonard): Signatures of Interstellar Chemistry? The Planetary Science Journal, 4, 8, doi: 10.3847/PSJ/aca64c

  26. [26]

    2010, VizieR Online Data Catalog: Collisional excitation of water in warm media (Faure+, 2008), VizieR Online Data Catalog, J/A+A/492/257

    Faure, A., & Josselin, E. 2010, VizieR Online Data Catalog: Collisional excitation of water in warm media (Faure+, 2008), VizieR Online Data Catalog, J/A+A/492/257

  27. [27]

    D., Cochran, A

    Feldman, P. D., Cochran, A. L., & Combi, M. R. 2004, in Comets II, ed. H. U. Keller & H. A. Weaver (University of Arizona Press), 425

  28. [28]

    Fray, N., B´ enilan, Y., Biver, N., et al. 2006a, Heliocentric evolution of the degradation of polyoxymethylene: Application to the origin of the formaldehyde (H2CO) extended source in Comet C/1995 O1 (Hale–Bopp), Icarus, 184, 239

  29. [29]

    Fray, N., B´ enilan, Y., Biver, N., et al. 2006b, Heliocentric evolution of the degradation of polyoxymethylene: Application to the origin of the formaldehyde (H2CO) extended source in Comet C/1995 O1 (Hale-Bopp), Icarus, 184, 239

  30. [30]

    F., Tsiganis, K., & Morbidelli, A

    Gomes, R., Levison, H. F., Tsiganis, K., & Morbidelli, A. 2005, Origin of the cataclysmic Late Heavy Bombardment period of the terrestrial planets, Nature, 435, 466

  31. [31]

    1957, Distribution d’intensit´ e dans la tˆ ete d’une com` ete, BSRSL, 43, 740

    Haser, L. 1957, Distribution d’intensit´ e dans la tˆ ete d’une com` ete, BSRSL, 43, 740

  32. [32]

    F., & Mukherjee, J

    Huebner, W. F., & Mukherjee, J. 2015, Photoionization and photodissociation rates in solar and blackbody radiation fields, P&SS, 106, 11, doi: 10.1016/j.pss.2014.11.022

  33. [33]

    Irvine, W., Bergman, P., Lowe, T., et al. 2004, Hcn and Hnc in Comets C/2000 Wm1 (Linear) and C/2002 C1 (Ikeya-Zhang), Origins of life and evolution of the biosphere : the journal of the International Society for the Study of the Origin of Life, 33, 609, doi: 10.1023/A:1025791101127

  34. [34]

    M., Bergin, E

    Irvine, W. M., Bergin, E. A., Dickens, J. E., et al. 1998, Chemical processing in the coma as the source of cometary HNC, Nature, 393, 547

  35. [35]

    2022, Outburst of comet C/2021 A1 (Leonard) on January 7, 2021, The Astronomer’s Telegram, 15189, 1

    Jehin, E., Moulane, Y., Manfroid, J., et al. 2022, Outburst of comet C/2021 A1 (Leonard) on January 7, 2021, The Astronomer’s Telegram, 15189, 1

  36. [36]

    2023, Disintegration of Long-period Comet C/2021 A1 (Leonard), The Astronomical Journal, 165, 122, doi: 10.3847/1538-3881/acb53b

    Jewitt, D., Kim, Y., Mattiazzo, M., et al. 2023, Disintegration of Long-period Comet C/2021 A1 (Leonard), The Astronomical Journal, 165, 122, doi: 10.3847/1538-3881/acb53b

  37. [37]

    F., Morbidelli, A., Tsiganis, K., Nesvorny, D., & Gomes, R

    Levison, H. F., Morbidelli, A., Tsiganis, K., Nesvorny, D., & Gomes, R. 2011, Late Orbital Instabilities in the Outer Planets Induced by Interaction with a Self-gravitating Planetesimal Disk, AJ, 142, 152

  38. [38]

    2006, Rotational excitation of carbon monosulfide by collisions with helium, A&A, 451, 1125, doi: 10.1051/0004-6361:20054363

    Lique, F., Spielfiedel, A., & Cernicharo, J. 2006, Rotational excitation of carbon monosulfide by collisions with helium, A&A, 451, 1125, doi: 10.1051/0004-6361:20054363

  39. [39]

    C., Bockel´ ee-Morvan, D., Boissier, J., et al

    Lis, D. C., Bockel´ ee-Morvan, D., Boissier, J., et al. 2008, Hydrogen Isocyanide in Comet 73P/Schwassmann-Wachmann (Fragment B), ApJ, 675, 931

  40. [40]

    A., ¨Oberg, K

    Loomis, R. A., ¨Oberg, K. I., Andrews, S. M., et al. 2018, Detecting Weak Spectral Lines in Interferometric Data through Matched Filtering, The Astronomical Journal, 155, 182, doi: 10.3847/1538-3881/aab604

  41. [41]

    N., & Zolensky, M

    Raymond, S. N., & Zolensky, M. 2024, Comets III (University of Arizona Press)

  42. [42]

    Meier, R., Eberhardt, P., Krankowsky, D., & Hodges, R. R. 1993, The extended formaldehyde source in comet P/Halley, A&A, 277, 677 21

  43. [43]

    N., Remijan, A

    Milam, S. N., Remijan, A. J., Womack, M., et al. 2006, Formaldehyde in Comets C/1995 O1 (Hale-Bopp), C/2002 T7 (LINEAR), and C/2001 Q4 (NEAT): Investigating the Cometary Origin of H2CO, The Astrophysical Journal, 649, 1169, doi: 10.1086/506501

  44. [44]

    F., Tsiganis, K., & Gomes, R

    Morbidelli, A., Levison, H. F., Tsiganis, K., & Gomes, R. 2005, Chaotic capture of Jupiter’s Trojan asteroids in the early Solar System, Nature, 435, 462

  45. [45]

    J., & Charnley, S

    Mumma, M. J., & Charnley, S. B. 2011, The Chemical Composition of Comets - Emerging Taxonomies and Natal Heritage, ARA&A, 49, 471 Nesvorn´ y, D., Vokrouhlick´ y, D., Dones, L., et al. 2017, Origin and Evolution of Short-period Comets, ApJ, 845, 27, doi: 10.3847/1538-4357/aa7cf6

  46. [46]

    B., et al

    Newville, M., Stensitzki, T., Allen, D. B., et al. 2016, LMFIT: Non-linear least-square minimization and curve-fitting for Python, Astrophysics Source Code Library, ascl

  47. [47]

    2005, in SF2A-2005: Semaine de l’Astrophysique Francaise, ed

    Pety, J. 2005, in SF2A-2005: Semaine de l’Astrophysique Francaise, ed. F. Casoli, T. Contini, J. M. Hameury, & L. Pagani, 721

  48. [48]

    Rabli, D., & Flower, D. R. 2010, The rotational excitation of methanol by molecular hydrogen, MNRAS, 406, 95, doi: 10.1111/j.1365-2966.2010.16671.x

  49. [49]

    M., Papadakis I

    Rodgers, S. D., & Charnley, S. B. 2005, Suprathermal chemical reactions driven by fast hydrogen atoms in cometary comae, MNRAS, 356, 1542, doi: 10.1111/j.1365-2966.2004.08606.x

  50. [50]

    X., Milam, S

    Roth, N. X., Milam, S. N., Cordiner, M. A., et al. 2021a, Leveraging the ALMA Atacama Compact Array for Cometary Science: An Interferometric Survey of Comet C/2015 ER61 (PanSTARRS) and Evidence for a Distributed Source of Carbon Monosulfide, The Astrophysical Journal, 921, 14, doi: 10.3847/1538-4357/ac0441

  51. [51]

    Black, J. H. 2005, An atomic and molecular database for analysis of submillimetre line observations, A&A, 432, 369, doi: 10.1051/0004-6361:20041729

  52. [52]

    2013, GILDAS: Grenoble Image and Line Data Analysis Software, Astrophysics Source Code Library, ascl

    Team, G. 2013, GILDAS: Grenoble Image and Line Data Analysis Software, Astrophysics Source Code Library, ascl

  53. [53]

    2013, Rotational quenching of H2CO by molecular hydrogen: cross-sections, rates and pressure broadening, MNRAS, 432, 2573, doi: 10.1093/mnras/stt616

    Wiesenfeld, L., & Faure, A. 2013, Rotational quenching of H2CO by molecular hydrogen: cross-sections, rates and pressure broadening, MNRAS, 432, 2573, doi: 10.1093/mnras/stt616

  54. [54]

    2022, Comets in Context: Comparing Comet Compositions with Protosolar Nebula Models, The Astrophysical Journal, 931, 164, doi: 10.3847/1538-4357/ac67e3

    Willacy, K., Turner, N., Bonev, B., et al. 2022, Comets in Context: Comparing Comet Compositions with Protosolar Nebula Models, The Astrophysical Journal, 931, 164, doi: 10.3847/1538-4357/ac67e3

  55. [55]

    Farnham, T. L. 2021, Preview of Comet C/2021 A1 (Leonard) and Its Encounter with Venus, The Astronomical Journal, 162, 194, doi: 10.3847/1538-3881/ac19ba

  56. [56]

    E., Zega, T

    Zolensky, M. E., Zega, T. J., Yano, H., et al. 2006, Mineralogy and Petrology of Comet 81P/Wild 2 Nucleus