FeI and NiI in cometary atmospheres. Connections between the NiI/FeI abundance ratio and chemical characteristics of Jupiter-family and Oort-cloud comets

FeI and NiI emission lines have recently been found in the spectra of 17 Solar System comets observed at heliocentric distances between 0.68 and 3.25 au and in the interstellar comet 2I/Borisov. The blackbody equilibrium temperature at the nucleus surface is too low to vaporize the refractory dust grains that contain metals, making the presence of iron and nickel atoms in cometary atmospheres a puzzling observation. Moreover, the measured NiI/FeI abundance ratio is on average one order of magnitude larger than the solar photosphere value. We report new measurements of FeI and NiI production rates and abundance ratios for the Jupiter-family comet (JFC) 46P/Wirtanen in its 2018 apparition and from archival data of the Oort-cloud comet (OCC) C/1996 B2 (Hyakutake). The comets were at geocentric distances of 0.09 au and 0.11 au, respectively. The emission line surface brightness was found to be inversely proportional to the projected distance to the nucleus, confirming that FeI and NiI atoms are ejected from the surface of the nucleus or originate from a short-lived parent. Considering the full sample of 20 comets, we find that the range of NiI/FeI abundance ratios is significantly larger in JFCs than in OCCs. We also unveil significant correlations between NiI/FeI and C$_2$/CN, C$_2$H$_6$/H$_2$O, and NH/CN. Carbon-chain- and NH-depleted comets show the highest NiI/FeI ratios. The existence of such relations suggests that the diversity of NiI/FeI abundance ratios in comets could be related to the cometary formation rather than to subsequent processes~in~the~coma.Comment: Accepted for publication in A&A Letter

Dust modelling and a dynamical study of comet 41P/Tuttle-Giacobini-Kresak during its 2017 perihelion passage

Thanks to the Rosetta mission, our understanding of comets has greatly improved. A very good opportunity to apply this knowledge appeared in early 2017 with the appearance of the Jupiter family comet 41P/TGK. We performed an observational campaign with the TRAPPIST telescopes that covered almost the entire period of time when the comet was active. In this work we present a comprehensive study of the evolution of the dust environment of 41P based on observational data from January to July, 2017. Also, we performed numerical simulations to constrain its origin and dynamical nature. To model the observational data set we used a Monte Carlo dust tail model, which allowed us to derive the dust parameters that best describe its dust environment as a function of heliocentric distance. In order to study its dynamical evolution, we completed several experiments to evaluate the degree of stability of its orbit, its life time in its current region close to Earth, and its future behaviour. From the dust analysis, we found that comet 41P has a complex emission pattern that shifted from full isotropic to anisotropic ejection sometime during February 24-March 14 in 2017, and then from anisotropic to full isotropic again between June 7-28. During the anisotropic period, the emission was controlled by two strongly active areas, where one was located in the southern and one in the northern hemisphere of the nucleus. The total dust mass loss is estimated to be $\sim7.5\times10^{8}$ kg. From the dynamical simulations we estimate that $\sim$3600 yr is the period of time during which 41P will remain in a similar orbit. Taking into account the estimated mass loss per orbit, after 3600 yr, the nucleus may lose about 30$\%$ of its mass. However, based on its observed dust-to-water mass ratio and its propensity to outbursts, the lifetime of this comet could be much shorter.Comment: 14 pages, 13 figures. Accepted for its publication in Astronomy & Astrophysic

TRAPPIST comet production rates: 19P/Borrelly, 67P/C-G, 104P/Kowal 2, C/2019 L3 (ATLAS), and C/2021 A1 (Leonard)

E. Jehin, Y. Moulane, J. Manfroid, M. Vander Donckt, F. Pozuelos, M. Ferrais, and D. Hutsemekers (STAR Institute, University of Liege) report that they obtained from TRAPPIST-North (code=Z53, Morocco) and TRAPPIST-South (code=I40, Chile) robotic telescopes (Jehin el al. 2011) recent observations using cometary HB narrowband filters (Farnham et al. 2000) for the following comets and computed preliminary production rates at 10.000 km using a Haser Model (Vp=Vd=1km/s) (Haser 1957)

Outburst of comet C/2021 A1 (Leonard) on January 7, 2021

We report about an outburst of comet C/2021 A1 (Leonard) that happened between January 6 and 8, 00:00 UT 2021. Observations were carried on with TRAPPIST-South (code=I40, Chile) robotic telescopes about every night since Dec 20 until Jan 25, using broad band and cometary HB narrowband filters (Jehin el al. 2011)

2I/Borisov: A C2 depleted interstellar comet

Aims. The discovery of the first active interstellar object 2I/Borisov provides an unprecedented opportunity to study planetary formation processes in another planetary system. In particular, spectroscopic observations of 2I allow us to constrain the composition of its nuclear ices. Methods. We obtained optical spectra of 2I with the 4.2 m William Herschel and 2.5 m Isaac Newton telescopes between 2019 September 30 and October 13, when the comet was between 2.5 au and 2.4 au from the Sun. We also imaged the comet with broadband filters on 15 nights from September 11 to October 17, as well as with a CN narrow-band filter on October 18 and 20, with the TRAPPIST-North telescope. Results. Broadband imaging confirms that the dust coma colours (B − V = 0.82 ± 0.02, V − R = 0.46 ± 0.03, R − I = 0.44 ± 0.03, B − R = 1.28 ± 0.03) are the same as for Solar System comets. We detect CN emission in all spectra and in the TRAPPIST narrow-band images with production rates between 1.6 × 1024 and 2.1 × 1024 molec/s. No other species are detected. We determine three-sigma upper limits for C2, C3, and OH production rates of 6 × 1023 molec/s, 2 × 1023 molec/s and 2 × 1027 molec/s, respectively, on October 02. There is no significant increase of the CN production rate or A(0)fρ during our observing period. Finally, we place a three-sigma upper limit on the Q(C2)/Q(CN) ratio of 0.3 (on October 13). From this, we conclude that 2I is highly depleted in C2, and may have a composition similar to Solar System carbon-chain depleted comets

Pluto's lower atmosphere and pressure evolution from ground-based stellar occultations, 1988-2016

Context. The tenuous nitrogen (N2) atmosphere on Pluto undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has recently (July 2015) been observed by the New Horizons spacecraft. Aims. The main goals of this study are (i) to construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) to constrain the structure of the lower atmosphere using a central flash observed in 2015. Methods. Eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between altitude levels of ~5 and ~380 km (i.e. pressures from ~ 10 μbar to 10 nbar). Results. (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived. (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia; and/or (b) hazes with tangential optical depth of ~0.3 are present at 4–7 km altitude levels; and/or (c) the nominal REX density values are overestimated by an implausibly large factor of ~20%; and/or (d) higher terrains block part of the flash in the Charon facing hemisphere

Theoretical study of ArH+ dissociative recombination and electron-impact vibrational excitation

Cross-sections are presented for dissociative recombination and electron-impact vibrational excitation of the ArH+ molecular ion at electron energies appropriate for the interstellar environment. The R-matrix method is employed to determine the molecular structure data, i.e. the position and width of the resonance states. The cross-sections and the corresponding Maxwellian rate coefficients are computed using a method based on the Multichannel Quantum Defect Theory. The main result of the paper is the very low dissociative recombination rate found at temperatures below 1000K. This is in agreement with the previous upper limit measurement in merged beams and offers a realistic explanation for the presence of ArH+ in exotic interstellar conditions