Comet-H: Comparative Comet Observations

• Comet-H is a comparative framework that unifies diverse comet observations—optical, isotopic, and radio—to study composition and dynamics.
• High-resolution studies reveal key metrics such as water D/H ratios, [OI] green-to-red signatures, and enhanced deuterium in organic molecules.
• Multi-scale observations, including CN imaging and radio continuum mapping, provide actionable insights into nucleus rotation and coma evolution.

In these studies, Comet-H is associated with multiple comet-focused observational programs that connect high-resolution optical spectroscopy, isotopic analysis, radio diagnostics, and nucleus-rotation studies. The term is used in connection with the Hanle Echelle Spectrograph program on C/2015 V2 (Johnson), 46P/Wirtanen, 41P/Tuttle–Giacobini–Kresák, and 38P/Stephan–Oterma; with GMRT observations of C/2020 F3 (NEOWISE); and with the CN-morphology study of 103P/Hartley 2. Comparative isotopic context is supplied by Rosetta/ROSINA measurements of 67P/Churyumov–Gerasimenko and Herschel/HIFI constraints on 45P/Honda–Mrkos–Pajdušáková (Aravind et al., 2024, Pal et al., 2021, Samarasinha et al., 2011, Müller et al., 2022, Lis et al., 2013).

1. Observational scope and comparative framework

Across the cited studies, Comet-H spans several measurement regimes. Optical echelle spectroscopy resolves vibrational bands and rotational lines in the coma; in situ mass spectrometry constrains isotopic fractionation in water and simple organics; submillimeter heterodyne spectroscopy measures water isotopologues and production rates; radio interferometry probes centimeter-wave continuum and H I 21-cm absorption; and CN coma imaging tracks rotational dynamics through repeatable jet morphology. This suggests that Comet-H is best understood as a comparative research construct centered on cometary composition, coma physics, and nucleus dynamics rather than as a single narrow observational technique.

Observational context	Facility	Principal output
Optical spectroscopy of several comets	HESP on the 2-m HCT	[OI] green-to-red ratios, NH$_2$ OPR, resolved CN/C$_2$/NH$_2$/C$_3$/CH structure
In situ isotopic analysis of 67P	Rosetta/ROSINA DFMS	Stable water D/H, $^{16}\mathrm{O}/^{17}\mathrm{O}$, and alkane isotope ratios
Submillimeter isotopologues in 45P/HMP	Herschel/HIFI	$\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$ (3$\sigma$)
Radio study of C/2020 F3 (NEOWISE)	GMRT	1050–1450 MHz continuum and H I 21-cm absorption
CN morphology of 103P/Hartley 2	KPNO 2.1-m telescope	Changing effective rotation period and low-excitation NPA rotation

The comparative value of this framework lies in the fact that different diagnostics sample different physical components of a comet. ROSINA sampled the inner coma along Rosetta’s trajectory at 67P, Herschel/HIFI measured line areas and modeled production rates for 45P, HESP isolated optical radicals and forbidden oxygen lines in several comets, GMRT detected a large neutral-hydrogen column and centimeter-wave continuum in NEOWISE, and CN imaging of Hartley 2 tracked jet-driven rotational evolution (Aravind et al., 2024, Müller et al., 2022, Lis et al., 2013, Pal et al., 2021, Samarasinha et al., 2011).

2. High-resolution optical spectroscopy with HESP

Within the HESP program, Comet-H refers specifically to the set of comets observed with the Hanle Echelle Spectrograph on the 2-m Himalayan Chandra Telescope at Hanle. HESP is a bench-mounted, dual fibre-fed echelle spectrograph that uses an R2 echelle grating, two cross-dispersing prisms, and a 4k×4k E2V CCD. Its resolving power is defined as $R = \lambda / \Delta \lambda$, with operation at $R \approx 30{,}000$ in medium resolution and $R \approx 60{,}000$ in high resolution, over an effective wavelength range of 3700–10,000 Å. In star–sky mode used for comets, both fibres are filled by coma emission, so the two fibres simultaneously sample distinct nucleocentric regions separated by approximately 13 arcsec on the sky. Calibration uses ThAr for wavelength and quartz blue for flats; orders are continuum-normalised, frames are co-added to increase S/N, and Doppler corrections for geocentric velocity are applied with IRAF dopcor (Aravind et al., 2024).

The program analyzed four comets from distinct reservoirs and dynamical classes: the long-period comet C/2015 V2 (Johnson), the Jupiter-family comets 46P/Wirtanen and 41P/Tuttle–Giacobini–Kresák, and the Halley-type comet 38P/Stephan–Oterma. Across 3700–10,000 Å, strong emissions were detected from CN, C$_2$0, NH$_2$1, C$_2$2, and CH, while CH$_2$3 and CO$_2$4 were too weak to be distinguished from noise. CN(B$_2$5–X$_2$6) in the violet, especially the (0–0) band, showed clearly resolved P and R rotational branches, and CN(A$_2$7–X$_2$8) was detected in the red. Rotational lines were cross-identified against laboratory–comet line lists by requiring wavelength coincidence within $_2$9 and $_2$0 significance above the local continuum noise. The Swings effect drives epoch-to-epoch changes in relative line strengths in CN and CH because of solar Fraunhofer line pumping and heliocentric velocity. C$_2$1 Swan bands with $_2$2 were widely detected, and the $_2$3 band head at 5165 Å was strong in 46P and C/2015 V2. High-resolution spectra revealed substructure in features that appear blended at $_2$4, including deblending of NH$_2$5(0–5–0) from parts of the C$_2$6 $_2$7 region (Aravind et al., 2024).

A central result concerns the forbidden oxygen lines [OI] 5577.339 Å, 6300.304 Å, and 6363.776 Å. HESP’s resolution and the comet’s geocentric Doppler shift made it possible to separate cometary [OI] from telluric [OI]; when line separation was marginal, Gaussian deblending was applied. The measured green-to-red ratios were $_2$8 for C/2015 V2, $_2$9 for 46P, and $_3$0 for 41P, indicating H$_3$1O as the dominant parent of [OI] in the sample. The 46P dual-fibre observation on 2018-12-28 is especially diagnostic: the fibre centered on the photocentre, probing within approximately 200 km, yielded $_3$2, while the second fibre about 1000 km away yielded $_3$3. This is the textbook signature of collisional quenching of O($_3$4D) in the dense inner coma. Away from the nucleus, the ratio relaxes to a water-dominated photochemical value. The same dataset yielded an NH$_3$5 ortho-to-para ratio for 46P of approximately $_3$6, an inferred NH$_3$7 OPR of approximately $_3$8, and a spin temperature of $_3$9. The physical meaning of nuclear spin temperatures remains debated, and the paper explicitly notes that a full quantification of the uncertainty in $^{16}\mathrm{O}/^{17}\mathrm{O}$0 would require forward modeling of NH$^{16}\mathrm{O}/^{17}\mathrm{O}$1 populations and radiative transfer (Aravind et al., 2024).

3. In situ isotopic fractionation in 67P/Churyumov–Gerasimenko

Rosetta’s escort of 67P/Churyumov–Gerasimenko provides the most detailed isotopic component of the Comet-H corpus. ROSINA/DFMS, a Nier–Johnson double-focusing mass spectrometer with $^{16}\mathrm{O}/^{17}\mathrm{O}$2 at the 1% level on $^{16}\mathrm{O}/^{17}\mathrm{O}$3 28, ionized neutral gas at 45 eV, transmitted selected $^{16}\mathrm{O}/^{17}\mathrm{O}$4 to a Chevron MCP, and recorded position-resolved counts on the LEDA detector. Isotopic ratios were derived from integrated peak areas of fitted DFMS spectral lines, with each line modeled as a double Gaussian whose relative amplitudes and widths were defined by instrument calibration and held common across peaks within a spectrum. A calibrated mass scale mapped pixel position to $^{16}\mathrm{O}/^{17}\mathrm{O}$5, and overlapping peaks were deconvolved in the fit. For water, back-to-back measurements of $^{16}\mathrm{O}/^{17}\mathrm{O}$6 18 and 19 within one minute enabled robust isotopic ratios under nearly identical conditions. The definitions used were $^{16}\mathrm{O}/^{17}\mathrm{O}$7 and $^{16}\mathrm{O}/^{17}\mathrm{O}$8. For alkanes $^{16}\mathrm{O}/^{17}\mathrm{O}$9, the study used $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$0 and $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$1 (Müller et al., 2022).

Across four mission phases—first equinox, perihelion, peak gas production, and second equinox—a total of 150 spectra around $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$2 18–19 yielded a relative overall mean $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$3. The phase means, $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$4, $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$5, $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$6, and $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$7, were all $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$8-compatible with the overall value. Including systematic uncertainties from detector gain and pixel-gain corrections gave an absolute mean $\mathrm{D/H}(\mathrm{H_2O}) < 2.0 \times 10^{-4}$9. The oxygen isotopic ratio from the same dataset was $\sigma$0 relative, or $\sigma$1 absolute when systematics were included, with a relative $\sigma$2 variation of approximately 2.3%. Proper inclusion of the $\sigma$3 contribution at $\sigma$4 19 was essential for a stable oxygen ratio across the mission. Crucially, water D/H showed no dependence on heliocentric distance in the 1.24–2.65 au interval of the core analysis, on water production rate, on observational phase angle, on sub-spacecraft latitude, or on geometry above different hemispheres. The study therefore concluded that the D/H ratio in 67P’s coma is independent of heliocentric distance, level of cometary activity, and spacecraft location with respect to the nucleus (Müller et al., 2022).

The same work extended isotopic analysis to the first four linear alkanes in the coma: methane, ethane, propane, and butane. DFMS detected CH$\sigma$5 in mid-August to early September 2016, C$\sigma$6H$\sigma$7 in October 2014 and March 2016, C$\sigma$8H$\sigma$9 in March 2016, and C$R = \lambda / \Delta \lambda$0H$R = \lambda / \Delta \lambda$1 in March 2016.

Alkane	D/H	$R = \lambda / \Delta \lambda$2
CH$R = \lambda / \Delta \lambda$3	$R = \lambda / \Delta \lambda$4	$R = \lambda / \Delta \lambda$5
C$R = \lambda / \Delta \lambda$6H$R = \lambda / \Delta \lambda$7	$R = \lambda / \Delta \lambda$8	$R = \lambda / \Delta \lambda$9
C$R \approx 30{,}000$0H$R \approx 30{,}000$1	$R \approx 30{,}000$2	$R \approx 30{,}000$3
C$R \approx 30{,}000$4H$R \approx 30{,}000$5	$R \approx 30{,}000$6	$R \approx 30{,}000$7

These alkane D/H values are 4.1–4.8 times higher than D/H in H$R \approx 30{,}000$8O, while the carbon isotopic ratios are compatible, within uncertainties, with previously measured $R \approx 30{,}000$9 in 67P’s CO$R \approx 60{,}000$0 and CO and with typical Solar System values. The paper states that the relatively high D/H ratio in alkanes is in line with other cometary organic molecules and suggests that these organics may be inherited from the presolar molecular cloud from which the Solar System formed. It also notes that 67P’s water, at approximately $R \approx 60{,}000$1, is more than three times terrestrial VSMOW and among the highest measured for a Jupiter-family comet, whereas some hyperactive JFCs such as 103P/Hartley 2 exhibit Earth-like D/H. This comparison underscores isotopic diversity within the comet population (Müller et al., 2022).

4. Remote submillimeter D/H constraints in 45P/Honda–Mrkos–Pajdušáková

Herschel/HIFI observations of 45P/Honda–Mrkos–Pajdušáková provide a complementary remote-sensing constraint on water deuteration. 45P/HMP is a short-period Jupiter-family comet with orbital period 5.3 yr, eccentricity $R \approx 60{,}000$2, perihelion distance $R \approx 60{,}000$3 AU, and nuclear radius $R \approx 60{,}000$4 km. Observations were obtained on 2011 Aug 13 at heliocentric distance $R \approx 60{,}000$5 AU and geocentric distance $R \approx 60{,}000$6 AU, using frequency-switched single-point spectra interleaved between species and on-the-fly mapping. The targeted transitions were $R \approx 60{,}000$7 $R \approx 60{,}000$8 at 556.939 GHz, $R \approx 60{,}000$9 $_2$00 at 547.676 GHz, HDO $_2$01 at 509.292 GHz, and $_2$02 $_2$03 at 987.927 GHz. Production rates were derived with a non-LTE excitation and radiative transfer framework using $_2$04 km s$_2$05, a gas temperature decreasing from 150 K near the nucleus to 10 K at approximately 80 km, an electron density scaling factor $_2$06, and a contact-surface scaling factor $_2$07 (Lis et al., 2013).

The measured line areas, corrected for main-beam efficiency, were $_2$08 mK km s$_2$09 for $_2$10, $_2$11 mK km s$_2$12 for $_2$13, and $_2$14 mK km s$_2$15 for HDO at $_2$16. These corresponded to $_2$17, $_2$18, and $_2$19. For low D/H in water, the study used $_2$20 and $_2$21. A direct upper limit based on optically thick $_2$22 gave $_2$23 at $_2$24, while the optically thin HDO/$_2$25 ratio with $_2$26 implied $_2$27 at $_2$28. After propagating modeling and isotopic uncertainties, the final conservative constraint was $_2$29 at $_2$30 (Lis et al., 2013).

This upper limit is consistent with terrestrial VSMOW and with the Herschel result for 103P/Hartley 2, $_2$31, while the pre-Herschel canonical Oort-cloud value $_2$32 is excluded for 45P/HMP at the 4.5$_2$33 level. The study therefore concludes that cometary D/H ratios are diverse and that previously observed high values are not universal across the comet population. It further notes that the small sample suggests JFCs may tend to lower D/H than Oort-cloud comets, although the statistical basis remains preliminary. Because the HDO/$_2$34 ratio uses optically thin lines, it is less sensitive to radiative-transfer details and coma anisotropy than estimates based directly on the main isotopologue. The same paper identifies future prospects with CCAT and CHAI, including ground-accessible HDO lines at 509, 464, and 894 GHz and an estimated programmatic yield of approximately 15 HDO detections over a decade (Lis et al., 2013).

5. Radio continuum and H I 21-cm absorption in C/2020 F3 (NEOWISE)

GMRT observations extend the Comet-H framework into the centimeter-radio regime. C/2020 F3 (NEOWISE) was observed on 2020-07-30 from 08:43 to 11:34 UTC with the Giant Metrewave Radio Telescope, using thirty 45 m dishes with baselines of approximately 0.1–25 km and a dual-polarization L-band receiver. The comet was at 0.79 AU from the Sun and 0.76 AU from Earth. The frequency range 1050–1450 MHz included the H I 21-cm hyperfine line at rest frequency 1420.405 MHz. The spectral setup used 8192 channels at native resolution 0.0488 MHz; for continuum imaging, channels were averaged to 2.0 MHz bins. Calibration used 3C 286 for flux and bandpass on the Perley–Butler 2017 flux-density scale and 1227+365 for phase. Data reduction in CASA included extensive flagging, continuum imaging with tclean using w-projection and multifrequency synthesis, self-calibration with gaincal/applycal, primary-beam correction with wpbgmrt, and uv-plane continuum subtraction with uvcontsub (Pal et al., 2021).

The continuum emission was detected at flux-density levels of $_2$35 to $_2$36 mJy across 1050–1450 MHz in a 91″ aperture. A single power-law fit $_2$37 gave $_2$38, indicating a rising continuum with frequency over the observed band. At 1378 and 1430 MHz, the emission is elongated sunward, opposite the comet’s motion, with the peak continuum aligned with the JPL Horizons ephemeris position. The paper attributes the L-band continuum to a large Icy Grains Halo comprising volatile-rich particles surrounding the nucleus. Because dust continuum is efficient when particle size $_2$39, detection at $_2$40 cm implies unusually large grains or icy pebbles in the coma or halo. Thermal dust emission from cm-to-dm sized icy grains is described as consistent with the positive $_2$41, but $_2$42 is lower than the canonical Rayleigh–Jeans slope $_2$43, so the mechanism cannot be pinned down definitively; the paper states that a contribution from partially optically thick free–free emission is plausible, whereas optically thin free–free emission with $_2$44 is disfavored (Pal et al., 2021).

The same dataset yielded an H I 21-cm absorption detection with SNR approximately 5.7. A single-Gaussian MCMC fit gave a peak optical depth $_2$45 and FWHM $_2$46 km s$_2$47. The spectrum was extracted by disk-averaging over a 91″ circular aperture centered on the comet. Using the standard absorption relation and assuming spin temperature $_2$48 K and filling factor $_2$49, the reported H I column density was $_2$50. The paper does not provide the fitted line-center velocity, so velocity centroids relative to the comet and the LSR cannot be quoted from this dataset. The linewidth is much broader than contemporaneous OH widths of approximately 10 km s$_2$51, which the paper attributes to fast H atoms produced by water and OH photodissociation at characteristic speeds of approximately 18–20 km s$_2$52, combined with additional broadening in the distributed hydrogen coma. The result is presented as the first interferometric high-resolution detection of a comet with GMRT and as a rare radio constraint on a large neutral hydrogen coma and its associated halo (Pal et al., 2021).

6. CN coma morphology and the rotational dynamics of 103P/Hartley 2

The rotational component of Comet-H is anchored by the ground-based CN monitoring of 103P/Hartley 2 during the 2010 apparition. The coma was imaged on twenty nights spanning four observing runs with the Kitt Peak National Observatory 2.1-m telescope: 2010 September 1–3, September 30–October 4, November 2–8, and December 11–15. These dates bracketed perigee on Oct 20, perihelion on Oct 28, and the EPOXI flyby on Nov 4. Narrowband CN images used a filter centered near 3869 Å with FWHM approximately 56 Å. After standard bias subtraction and flat-fielding, the images were enhanced by azimuthal profile division using the Schleicher and Farnham technique. Because Hartley 2 had a low dust-to-CN ratio, dust contamination in the CN band was regarded as negligible. The CN violet (0–0) fluorescence band isolated gas emitted from active regions, and evolving jets mapped the rotational state of the nucleus through spirals, arcs, and opposing features (Samarasinha et al., 2011).

The morphology changed systematically across the apparition. In September, the CN images showed a clear clockwise progression of an Archimedean spiral-like jet, with a periodicity near 17.1 h. In October the spiral persisted, a weaker partially overlapping feature appeared to the east or northeast, and the best repeatability occurred at $_2$53 h rather than 17.1 h. In November the morphology switched to two nearly opposite CN features with no spiral-like pattern, and the repeatability period grew to approximately 18.8 h. In December, the southern feature moved east $_2$54 south $_2$55 west, indicating counter-clockwise motion. The authors interpret the cycle-to-cycle differences at identical rotational phase as evidence for a low-excitation non-principal-axis rotational state. Synodic corrections were $_2$56 h and cannot explain the observed period changes (Samarasinha et al., 2011).

Using $_2$57, the paper reports an angular-speed decrease from approximately $_2$58 in September to $_2$59 in November, corresponding over approximately 60 days to an average spin-down rate $_2$60. In period space, this is approximately $_2$61, with stepwise increases of about $_2$62 from September to October and about $_2$63 from October to November. The angular-momentum relations used are $_2$64 and $_2$65. With $_2$66, $_2$67, $_2$68, $_2$69, and $_2$70, the characteristic period-change timescale $_2$71 is approximately 6 months, consistent with the observed month-scale increase in $_2$72. The corresponding outgassing torque estimate is $_2$73, while the observation-based estimate $_2$74 gives approximately $_2$75, an order-of-magnitude consistency that the paper regards as reasonable given uncertainties in lever arm and mass flux (Samarasinha et al., 2011).

A further result concerns damping. Hartley 2 remained in an excited, evolving NPA state over four months without evident damping toward principal-axis rotation. Using $_2$76, $_2$77, and $_2$78, the absence of damping implies $_2$79, where $_2$80 is the shear modulus and $_2$81 is the quality factor. The study interprets this as evidence that Hartley 2 is mechanically weak, flexible, and elastic relative to asteroids. Combined with the small size, low density, elongated shape, and activity concentrated near the ends of the long axis, the CN morphology and torque analysis yield a coherent picture of a porous nucleus in a low-excitation tumbling state whose rotation rate and spin axis evolve measurably under jet-induced torques (Samarasinha et al., 2011).

7. Synthesis, interpretive limits, and open problems

Several themes recur across the Comet-H literature. First, water is repeatedly identified as a dominant parent species in observable coma diagnostics: low [OI] green-to-red ratios away from the nucleus indicate H$_2$82O-dominated oxygen emissions in the HESP sample, while in 67P the water D/H ratio remains stable through perihelion and peak outgassing, and in 45P the low upper limit on D/H is consistent with terrestrial-like values in at least some JFCs. Second, cometary isotopic diversity is explicit rather than exceptional. 67P’s water has $_2$83, whereas 45P has $_2$84 and 103P/Hartley 2 is quoted in comparison at $_2$85. Third, organics need not track water isotopically: in 67P, simple linear alkanes are enriched in deuterium by factors of 4.1–4.8 relative to H$_2$86O, while their carbon isotopes remain compatible with Solar System norms. Fourth, coma structure is multiscale: HESP resolves radial and spatial changes in [OI] quenching and NH$_2$87 OPR, GMRT detects a large neutral hydrogen column and centimeter-wave continuum from an Icy Grains Halo, and CN imaging of Hartley 2 shows that coma morphology can directly encode non-principal-axis rotation and evolving torques (Aravind et al., 2024, Müller et al., 2022, Lis et al., 2013, Pal et al., 2021, Samarasinha et al., 2011).

The literature also places clear limits on interpretation. ROSINA/DFMS sampled the inner coma along Rosetta’s trajectory, so localized isotopic heterogeneity could remain undersampled or mixed away in the coma. Herschel/HIFI relied on spherical-symmetry assumptions, an adopted $_2$88, and a simplified temperature profile, although the optically thin HDO/$_2$89 ratio is relatively robust. HESP ratios are robust against the absence of absolute flux calibration, but residual blends, variable S/N across orders, and nearby C$_2$90 contamination near 5577 Å remain possible systematics. The NH$_2$91 spin temperature inferred from NH$_2$92 OPR is explicitly described as physically debated. In NEOWISE, the 1050–1450 MHz slope alone cannot definitively separate thermal dust from free–free emission, and the paper does not provide the H I line-center velocity. In Hartley 2, December’s period was not robustly measured because of limited nightly coverage (Aravind et al., 2024, Lis et al., 2013, Pal et al., 2021, Samarasinha et al., 2011).

The forward program implied by these studies is correspondingly broad. The 67P paper calls for long-term, multi-comet in situ isotopic monitoring, high-resolution remote spectroscopy across perihelion, and ultimately cryogenic sample return from multiple nucleus locales. The 45P paper identifies CCAT and ALMA as essential for expanding the HDO database. The HESP paper points to future work on isotopic lines, detailed C$_2$93/NH$_2$94 deblending, and comprehensive modeling of [OI] photochemistry and spin temperatures as non-sidereal tracking and pipeline tools mature. The GMRT study calls for deeper integrations, broader frequency coverage, and interferometric mapping of both continuum and H I absorption at multiple epochs. Together these directions suggest that Comet-H, understood as a comparative observational framework, is primarily valuable because it links isotopic inheritance, coma excitation, photochemical processing, grain-halo structure, and rotational dynamics within a single technically coherent view of cometary physics (Müller et al., 2022, Lis et al., 2013, Aravind et al., 2024, Pal et al., 2021, Samarasinha et al., 2011).