Activity and composition of the hyperactive comet 46P/Wirtanen during its close approach in 2018

Context. Hyperactive comets are a small group of comets whose activity is higher than expected. They seem to emit more water than would normally be expected given the size of their nucleus. Comet 46P/Wirtanen (hereafter, 46P) is among these objects of interest. Investigating its activity and composition evolution could provide clues about its origins and its formation region in the Solar nebulae. Aims. Given the exceptional close approach of comet 46P to the Earth in 2018, we aim to study the evolution of its activity and composition as a function of heliocentric distances before and after perihelion. Methods. We used both TRAPPIST telescopes to monitor the comet for almost a year with broad and narrow-band ﬁlters. We derived the production rates of ﬁve gaseous species (OH, NH, CN, C 3 , and C 2 ) using a Haser model as well as the A( θ )f ρ dust proxy parameter. The comet was also observed with the two optical high-resolution spectrographs UVES and ESPRESSO, mounted on the 8-m ESO VLT, to measure the isotopic ratios of C and N, along with the oxygen forbidden-line ratios and the NH 2 ortho-to-para ratios. Results. Over nearly a year, we followed the rise and decline of the production rates of different species as well as the dust activity of 46P on both pre-and post-perihelion. Relative abundances with respect to CN and OH along the orbit of the comet show constant and symmetric abundance ratios and a typical coma composition. We determined the rotation period of the nucleus using high-cadence observations and long series of CN images on several nights. We obtained a value of (9.18 ± 0.05) hours at perihelion. Using the high-resolution spectra of 46P coma, we derived C and N isotopic ratios of 100 ± 20 and 150 ± 30 as well as a green-to-red forbidden oxygen [OI] line ratio of 0.23 ± 0.02. We measured a NH 2 ortho-to-para ratio of 3.31 ± 0.03 and derived an ammonia ratio of 1.19 ± 0.03, corresponding to a spin temperature of 27 ± 1 K. Conclusions. Narrow-band observations show that comet 46P is a hyperactive comet for which 40% of its nucleus surface is active. It has a typical composition, similar to other normal comets; however, an asymmetric behavior with respect to perihelion has been seen in its activity, which is typical of seasonal effects. Photometric measurements show no evidence for a change in the rotation period of the nucleus during this apparition. High-resolution spectra show that 46P has typical NH 2 ortho-to-para, [OI] lines ratios, and C and N isotopic ratios.


Introduction
46P/Wirtanen is a Jupiter-family comet (hereafter, JFC), with an orbital period of 5.5 years, that was discovered on January 17, 1948 by Carl Wirtanen at the Lick Observatory (Jeffers 1948).It has been classified as an hyperactive comet, namely, it is part of a small group of comets whose activity levels are higher than expected based on the sizes of their nuclei.It was an initial target of the Rosetta mission and many studies have been performed to determine its nucleus size, shape, rotation period, albedo, and color, as well as its gas and dust activity (Farnham & Schleicher 1998;Lamy et al. 1998;Schulz et al. 1998;Bertaux et al. 1999;Crovisier et al. 2002).New radar observations have estimated the nucleus radius to be (0.70±0.05) km during its 2018 passage 1 .(McKay et al. 2021;Khan et al. 2021;Roth et al. 2021a), search for icy grains in the coma (Protopapa et al. 2021;Bonev et al. 2021), and coma features from optical imaging (Knight et al. 2021).Millimeter and submillimeter observations of 46P have reported on the detection of several complex organic molecules and showed that the comet is relatively rich in methanol, but relatively depleted in CO, CS, HNC, HNCO, and HCOOH (Biver et al. 2021;Roth et al. 2021b;Coulson et al. 2020;Bergman et al. 2022).Using both TRAPPIST-North (TN) and TRAPPIST-South (TS), we collected photometric observation over almost a year.46P was visible for many hours during the night, allowing for long observing sequences for months during its return, passing near to the Sun.
In this work, we report on the activity, composition, and rotation period of comet 46P during its 2018 perihelion passage.An introduction and historical background of comet 46P is given in Section 1.We describe the observing circumstances and data reduction process in Section 2. In Section 3, we discuss the evolution of the gas and dust activity before and after perihelion and we compare the behavior of 46P to previous apparitions.The relative molecular abundances and their evolution with respect to the heliocentric distances are discussed in the same section.In Section 4, we measure the nucleus' rotation period and its variation around perihelion.In Section 5, we present the C and N isotopic ratios and the NH 2 (and NH 3 ) ortho-to-para ratio (OPR) and the green-to-red [OI] lines (G/R) ratio derived from highresolution spectra.The summary and conclusions of this work are given in Section 6.

TRAPPIST
We used both TN and TS telescopes (Jehin et al. 2011) to observe and follow comet 46P over almost a year.We started monitoring the comet at the beginning of August 2018 (r h =1.88 au, pre-perihelion) until the end of March 2019 (r h =1.70 au, postperihelion).More than 2400 broad and narrow-band images of the comet were collected over 45 nights with TS and 40 nights with TN. 46P reached its perihelion on December 12, 2018 at 1.06 au from the Sun and only at 0.08 au from Earth.For the data reduction, we followed standard procedures using frequently updated master bias, flat, and dark frames.The sky contamination was removed and a flux calibration was performed using standard stars.In order to derive the production rates, we converted the flux of different gaseous species (OH, NH, CN, C 3 , and C 2 ), measured through the HB narrow band cometary filters (Farnham et al. 2000), to the column densities using heliocentric distance and heliocentric velocity-dependent fluorescence efficiencies (Schleicher 2010, and references therein).To convert the column densities into production rates, we used a Haser model (Haser 1957) to fit the coma profile.The model adjustment was performed at a physical distance of 10 000 km from the nucleus.The scale lengths and g-factors of different molecules at 1 au were scaled by r −2 h (A' Hearn et al. 1995;Schleicher 2010).More details about the Haser model and its parameters are given in our previous works (Moulane et al. 2018;Moulane et al. 2020).We derived the water-production rate from our the OH production rates using the empirical formula Q(H 2 O)=1.36 r −0.5 h Q(OH), given in Cochran & Schleicher (1993); Schleicher et al. (1998).We derived the Afρ parameter, a proxy for the dust production (A'Hearn et al. 1984), from the dust profiles using the HB cometary dust continuum BC, GC, and RC filters and the broadband R and I filters.We corrected the A(θ)fρ for the phase angle effect to obtain A(0)fρ using the phase function given in Schleicher (2007).The main uncertainties in the gas-production rates and Afρ come from the absolute flux calibration and the sky background subtraction.For the absolute flux calibration, we used an uncertainty of 5% on the extinction coefficients in all filters as seen in our long-term observations of the standard stars.This is almost negligible compared to the sky level uncertainty at lower airmass, but it becomes significant at high airmasses.We estimated the three-sigma-level uncertainty on the sky background value and use it to computed the error on production rates due to the sky subtraction.As a result, errors given in the following sections are a quadratic combination of sky background and extinction coefficient uncertainties.

UVES and ESPRESSO at 8-m ESO VLT
We obtained a high-resolution spectrum of comet 46P with the Ultraviolet-Visual Echelle Spectrograph (UVES) mounted on the Unit 2 telescope (UT2) of the Very Large Telescopes (VLT) at the European Southern Observatory (ESO) on December 9, 2018 (a week after perihelion, r h =1.05 au and ∆=0.09 au).We used the UVES standard setting DIC#1 346+580 covering the range 3030 to 3880 Å in the blue and 4760 to 6840 Å in the red.We used a 0.44 wide slit, providing a resolving power of R∼80000.We took two spectra at UT 0:59:52 (exposure time of 2300 s) and UT 1:19:54 (exposure time of 2300 s) as a start time.
The ESO UVES pipeline (Dekker et al. 2000) was used to reduce the spectra in the extended object mode, in which the spatial information is kept.The spectra were corrected for extinction and flux calibrated using the UVES master response curve provided by ESO.One-dimensional spectra were then extracted by averaging the 2D spectra with simultaneous cosmic ray rejection and then corrected for the Doppler shift due to the velocity of the comet with respect to the Earth.More details about the UVES data reduction are given in the UVES manual 2 .The dustreflected sunlight was finally removed using a reference solar spectrum, BASS2000 3 .A more detailed description of the steps for the UVES data reduction and the solar spectrum subtraction is given in Manfroid et al. (2009) and references therein.
A set of six spectra were also obtained with the new Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO) at the VLT (Pepe et al. 2021) over two nights on December 9 and 10, 2018.Exposure times ranged from 5400 s for the first spectrum to 6300 s for the following five spectra.We used the ultra-high-resolution mode of the instrument, providing a spectral resolution of 190000 between 3800 and 7880 Å.In this mode, the spectrograph is fed by two 0.5 -size fibers: one centered on the object and another 7 away usually used for simultaneous observation of the sky or a wavelength drift reference.However, in the case of the very extended coma of 46P, both fibers contained comet signal and provided an opportunity to probe different parts of the coma.Observations were executed in visitor mode, with good seeing (typically less than 1 arcsec) and the data were reduced using the ESPRESSO pipeline.which is publicly available from the ESO pipeline repository4 .Moulane et al. 2022 -1.55 -1.48 -1.40 -1.32 -1.25 -1.19 -1.14 -1.10 -1.07 -1.06

Activity and composition
3.1.Gas-production rates and ratios We followed the evolution of the production rates of different gaseous species over a year.The derived production rates for each species and the A(0)fρ parameter, as well as their uncertainties are given in Table A.1.Figure 1 shows the logarithmic evolution of OH, NH, CN, C 2 , and C 3 production rates as a function of the distance to the Sun and time to perihelion.The gas activity of the comet was increasing slightly as the comet was approaching the Sun and decreased rapidly after perihelion, and more sharply for NH, OH and the dust.This fall-off could be due to the orientation of the active source regions on the surface of the comet toward the Earth before and after perihelion.The same behavior has been observed in the Lowell Discovery Telescope data (Knight et al. 2021).This explanation is supported by a study of the morphological features seen in 46P coma during this passage, reported by Knight et al. (2021).They show that 46P had large changes in both the sub-solar and sub-Earth latitudes during the apparition as well as changes in the pole orientation of the nucleus, assuming that there is no significant non-principal-axis rotation considered.It is also possible that the ice material present under the surface across the nucleus of the comet has different sublimation efficiency at the same heliocentric distances pre-versus post-perihelion.As mentioned above, we used narrow and broadband dust filters to estimate the dust production using the Afρ parameter introduced by A'Hearn et al. (1984) (see Figure 2).The behavior of the comet dust activity is similar to the gas on both sides of perihelion.This kind of asymmetry has been seen in many LPCs and JFCs.

Molecular abundances and dust-to-gas ratio
Using the production rates above, we derived the ratios with respect to CN and to OH as well as the dust-to-gas ratios.rate ratios with respect to CN and to OH as well as ratios of A(0)fρ-to-Q(OH) and A(0)fρ-to-Q(CN) with time to perihelion.The C 2 /CN ratio decreases when the comet gets far away from the Sun (40 days post-perihelion).The abundances ratios did not show any significant variation on both sides of perihelion.This implies that 46P has an homogeneous composition along its orbit around the Sun.Comet 46P has been observed with many ground-based telescopes over the previous passages using narrow-band photometry, since it was an important target selected initially for the Rosetta mission in 1994.Table 1 summarizes the mean relative molecular abundances in 46P measured in 2018 (this work), in the 1991 apparition A' Hearn et al. (1995), and in the 1997 apparition (Farnham & Schleicher 1998) compared to the database of comets given in (Schleicher 2008).The   (1995).We note that the difference in C 3 abundance in 2018 is due to the use of an updated dust continuum removal factor for the new C 3 filter, which results in an increase of the C 3 flux of about 2.1 times compared to the previous factors used in A'Hearn et al. (1995) and Schleicher et al. (1998).More details are given in Farnham et al. (2000) and Schleicher & Osip (2002).
measurements in the 1991, 1997 and 2018 apparitions show that 46P has a "typical" coma composition and a normal dust-to-gas ratio.Given the uncertainties on the ratios, no significant change is observed in the molecular abundances of 46P over the last four orbits.This shows that even an hyperactive comet, while outgassing a lot, does not change its main surface composition properties after several passages close to the Sun.

Parent and daughter molecules
In order to investigate the origin of the radicals observed in 46P's coma, we compared our daughter species abundances with parent molecular abundances obtained at Infrared wavelengths.Table 2 summarizes the abundances for some radicals observed in the optical compared to their possible parent abundances observed at NIR wavelengths (Roth et al. 2021a;Bonev et al. 2021;McKay et al. 2021;Khan et al. 2021), measured during the same period of time.We compared these abundances at two epochs, perihelion (r h =1.06 au) and post-perihelion (r h =1.12), when the infrared observations are available.Our CN/OH ratios are in agreement with the HCN/H 2 O ratios derived directly from IR observations around the same dates, namely, Dec. 16-23, 2018 (see Table 2).In addition, the HCN production rates (1.06±0.03)×10 25 molec/s obtained from millimeter (Biver et al. 2021) and (1.60±0.10)×10 25 molec/s from sub-millimeter observations (Bergman et al. 2022) show an agreement with our Q(CN) values obtained during the same nights of Dec. 15 and 20.This indicates that the observed CN abundances are consistent with HCN being the dominant parent.The C 2 abundance is much higher than C 2 H 2 which suggests that C 2 could be dissociated from both C 2 H 2 and C 2 H 6 .It has been shown that C 2 is still linked to both molecules even at large heliocentric distances ≥ 3 au in comet Hale-Bopp (Helbert et al. 2005).We looked for a C 3 parent but none has been identified in 46P at infrared wavelengths pointing to a C 3 production coming from chemical reactions in the coma.Based on direct detections of NH 3 lines, the mixing ratio of NH 3 /H 2 O ∼0.50 in 46P (See Table 2) is lower than the NH/OH ∼ 0.72 obtained with TRAPPIST.This result suggests that NH is not dissociated directly from NH 3 , but obviously via NH 2 which has been demonstrated for several comets (Shinnaka et al. 2011(Shinnaka et al. , 2016)).

Water-production rates and active area
We derived the water-production rate from OH using an empirical relationship (see Section 2) proposed by Cochran & Schleicher (1993) and Schleicher et al. (1998).Figure 4 shows a comparison of our water-production rates for the 2018 apparition with other measurements derived at different wavelengths using various techniques.First, we have a very good agreement between the TRAPPIST data (red circles) and the Lowell Observatory data (black squares) on both sides of perihelion, with both using the same technique based on narrow-band filters and the same model to derive the water-production rates.We also  2020) who derived the water-production rates from the hydrogen Lyman-α emission observed by the SWAN instrument on board SOHO show an offset on both sides of perihelion.This disagreement is probably the result of using different techniques and models.
It is also tricky to compare these results as SWAN/SOHO instrument threshold only allows for Lyman-α emission to be obtained within a small range of heliocentric distances.SWAN's values have rather large error-bars and dispersion among their measurements during this passage and the same for the previous ones (see Figure 5).We have seen this systematic difference also for comet 21P/Giacobini-Zinner (Moulane et al. 2020).
As 46P was the original target selected for Rosetta mission, many observations were performed in the previous passages in order to characterize the comet.Here, we compare the  water-production rates obtained in the previous passages 1991 and 1997 with those obtained in 2008 and 2018 (see Figure 5).Besides the offset between narrow-band and spectroscopic data, the comet's activity was similar in 1991 and 1997 but 30% lower in 2008-2018 apparitions.We did not include the 2002 observations in this comparison as they were dominated by large outbursts5 after perihelion (Combi et al. 2020).The maximum water-production rate in 2018 was measured ten days Comet 46P is considered a hyperactive comet as it is emitting more water than it would be expected from its nucleus radius.Based on the standard water vaporization model (Cowan & A'Hearn 1979)  6 , we estimate the evolution of the active area of 46P during our observational campaign in 2018 by modeling the water-production rate as a function of the heliocentric distances.We assumed a bond albedo of 5%, 100% infrared emissivity, and the rotational pole pointed at the Sun (see, e.g., A'Hearn et al. 1989;McKay et al. 2018McKay et al. , 2019)).We found that the active area of 46P varied from ∼0.9 km 2 at 1.52 au pre-perihelion, reached a maximum of ∼2.5 km 2 at 1.06 au during the perihelion passage, and decreased to ∼0.7 km 2 at 1.53 au post-perihelion.Based on a radius of ∼0.7 km, the active fraction of the nucleus was ∼14%, ∼40% and ∼11%, respectively.These values are much larger than the active fraction of less than 3% derived for the majority of JFCs that have radius measurements (A'Hearn et al. 1995).Since the comet had been observed in its previous passages, we compared the active fraction of its nucleus at perihelion in the past.Based on the water-production rates derived from the empirical conversion of Haser OH production rates (Schleicher et al. 1998;Farnham et al. 2007), we estimated the active fractions of its nucleus to be ∼72% at r h =1.11 au in 1991, ∼65% at r h =1.07 au in 1997, ∼37% at r h =1.05 au in 2008 and ∼40% at r h =1.06 au in 2018.It is clear that there is a decrease in the active fraction of the nucleus with time.Comet 103P/Hartley 2 was the first known "hyperactive" comet and it was visited by the Deep Impact eXtended Investigation mission, demonstrating that this hyperactivity was produced by sublimation of icy grains in the coma (A'Hearn et al. 2011;Protopapa et al. 2014).Several authors reported on the presence of an extended source of water ice sublimation in 46P coma in the form of small grains or large chunks.Bonev et al. (2021) provided indirect arguments for the presence of an extended source of water vapor in the coma of 46P at NIR spectra, but they did not provide information on the properties of the ice that produces the water vapor, such as particle size and ice-to-dust ratio, while Protopapa et al. (2021) and 6 https://pdssbn.astro.umd.edu/tools/ma-evap/index.shtml Kareta et al. (private communication) did not detect any water ice absorption features in their NIR spectra.

Rotation period and CN coma morphology
We took the opportunity provided by this very close approach of 46P to investigate the rotation period of its nucleus, using long CN series of images collected during the same night and at different epochs on both sides of the perihelion.Thanks to the visibility of the comet in both hemispheres, we collected long series of CN images for many hours on several nights with both telescopes.These series were taken over 12 nights with TS and 8 nights with TN.We determined the rotation period from CN light curves.We tested different size apertures and we finally used 20 which gave the largest amplitude in the light-curve.We phased multiple nights data to construct more extensive light curves and searched for the best alignment of the overlapping segments.We derived a rotation period of 9.12±0.05hr in late November and early December (pre-perihelion), 9.18±0.05hr between December 7 and 10 (perihelion time) and 9.00±0.04hr in mid January 2019 (post-perihelion).These results are in agreement (within the error bars) with those derived by Farnham et al. (2021) from the CN coma morphology and photometric series (see Table 3).Given the uncertainties on the measurements, we did not detect any significant change in the rotation period of the nucleus on both sides of the perihelion.Figure 6 shows the CN light curves at different epochs, showing two asymmetric maxima and minima, which could be due to two active areas on the surface of the nucleus.We searched for CN coma features using the rotational filter technique as described in our previous work (Moulane et al. 2018).Figure 7 shows the detection and evolution of a CN jet at different epochs and on both sides of perihelion.We did not detect any clear rotation of this jet with time in our images for many hours during the same night.This could be due to an orientation of the jet towards the Earth.Around November 16, another small jet appeared in the CN images but it disappeared in the following nights.This could suggest that the nucleus might be in a non-principal-axis (NPA) rotational state, as was previously mentioned by Samarasinha et al. (1996).The orientation of the jet remains the same before perihelion, while it changes by ∼ 95 degrees counterclockwise after perihelion because of the viewing geometry of the comet with respect to Earth.Farnham et al. (2021) used the 4.3-m Lowell Discovery Telescope to investigate the CN coma morphology of 46P.They were able to detect two spiral jets in the coma.The first one appears to have been active (at varying levels) throughout most of a rotation, corresponding to the bright jet detected in our data.The second seems to turn on and off with time and it could correspond to the small jet that appears around mid November in our images.Compared to the previous apparitions, Lamy et al. (1998) derived a rotation period of the nucleus of 6.0±0.3 using data from HST, while Meech et al. (1997) found a possible rotation of 7.6 hr and a very small amplitude variation of 0.09 mag using ground-based telescope during the 1997 apparition.Based on these results, we conclude that the rotation of the 46P nucleus might have slowed down by 16% to 45% over the last four orbits.These changes might be due to some unusual activity in the previous apparitions.There is some evidence for potentially significant outburst activity during the 2002 apparition (Combi et al. 2020) that could explain these changes.

Optical high-resolution spectra
5.1.The 12 C/ 13 C and 14 N/ 15 N isotopic ratios We obtained high-resolution observations with UVES/VLT (Dekker et al. 2000) of comet 46P on December 9, 2018 (r h =1.05 au and ∆=0.09 au).Using the 12 C 14 N B-X(0,0) band, we estimated the 12 C/ 13 C and 14 N/ 15 N isotopic ratios of 46P.We used a CN fluorescence model to create synthetic spectra of 13 C 14 N, 12 C 15 N, and 12 C 14 N.More details of the model are given in Manfroid et al. (2009).Figure 8 shows the observed CN spectrum compared to the synthetic one made under the same observing conditions.The ratios found for 12 C/ 13 C and 14 N/ 15 N are 100±20 and 150±30, respectively.These values are consistent with those of about 20 comets with different dynamical origins, 91.0±3.6 and 147.8±5.7 for 12 C/ 13 C and 14 N/ 15 N, respectively (Manfroid et al. 2009;Bockelée-Morvan et al. 2015).
5.2.The NH 2 and NH 3 ortho-to-para ratios From the UVES spectra, we measured the ortho-to-para (OPR, hereafter) abundance ratio of NH 2 using the three robivronic emission bands (0,7,0), (0,8,0), and (0,9,0) following the method described in Shinnaka et al. (2011).The derived OPRs of NH 2 and of its parent molecule NH 3 are listed for each band in  We would like to point out that recent laboratory experiments show that the OPR of water does not retain the memory of its formation temperature in the natal molecular cloud or in the solar nebula 4.6 Gyrs ago (Hama et al. 2011(Hama et al. , 2016;;Hama & Watanabe 2013).The OPR of water just after its sublimation from the solid phase is set to a statistical weight ratio (i.e., 3 for water).It is likely that this is also the case for ammonia (but the statistical weights ratio is unity for ammonia) although no experimental studies have been reported so far.The OPRs of cometary volatiles might have been modified by some kind of ortho-topara conversion processes in the inner coma such as water clusters, or by other catalyst activities of dust crust surfaces of the nucleus (Shinnaka et al. 2016).Thus, the observed OPRs might be diagnostic of the physico-chemical conditions in the innermost coma or beneath the surface.10) and used them to measure the ratio between the green line (557.73nm) and the sum of the two red lines (630.03 and 636.38 nm), commonly referred to as the G/R ratio.Due to the low geocentric velocity of the comet at the time of our observations (between -5.4 km/s and -4.5 km/s), the cometary and telluric lines were slightly blended, even at the very high resolution of ESPRESSO.We used two Gaussian to fit both lines simultaneously.We performed the measurements for the center and the sky fiber separately for the six spectra and subsequently averaged the values for each fiber.The uncertainty of the measurement is the standard deviation of the measurements made of six spectra.We measured a G/R ratio of 0.23±0.02 in the central fiber and of 0.05±0.01only in the sky fiber (7 away).Both measurements are in agreement with what has been measured in several comets at the same nucleocentric projected distances (Decock et al. 2015).The difference of G/R ratio at different distances from the nucleus is mainly due to the collisional quenching of O( 1 S) and O( 1 D) by water molecules in the inner coma.This fits with the general behavior observed in several comets, where the G/R ratio decreases monotonically with the projected distance to the nucleus (Decock et al. 2015).The intensity of these lines can be used to derive the CO/CO 2 /H 2 O abundances ratios and, hence, the C/O ratio, which can provide a diagnostic of solar nebula chemistry in the comet-forming region (Öberg et al. 2011).The detail analysis and modeling of the [OI] lines are out of the scope of this paper.

Summary and conclusions
We conducted an extensive monitoring of the Jupiter Family hyperactive comet 46P on both sides of perihelion with TRAP-PIST telescopes.These observations led us to an overview of the comet activity along its journey around the Sun in 2018-2019.
Using narrow and broadband filters with the TRAPPIST telescopes, we derived different quantities related to the activity of the comet such as the gas-production rates, the A(0)fρ parameter, the rotation period of the nucleus, and the morphology of the coma.Comet 46P shows an asymmetric activity with respect to perihelion, steeper decrease after perihelion for OH, NH, and dust (Afρ parameter).This might be due to the spin axis orientation or the distribution of activity on the comet's surface.The comet has exhibited a decrease in its activity from the 1997 to 2018 apparition by about 30%.This implies a decrease of the nucleus's active fraction to ∼40% which confirms its hyperactivity compared to most of JFCs.The molecular abundances relative to CN and OH remained the same as a function of the heliocentric distance, showing that 46P has a typical composition similar to what has been found in previous apparitions.These derived abundances in the optical showed a good correlation with their possible mother molecules measured in the IR (HCN, C 2 H 2 , and NH 3 ) and sub-millimeter (HCN), confirming the link between different molecules observed in several comets.Given the close approach of the comet, we obtained high spatial and high resolution spectra for the comet around the perihelion.Using UVES spectra, we measured typical 12 C/ 13 C and 14 N/ 15 N isotopic ratios of 100±20 and 150±30, and average NH 2 OPR of 3.31±0.03and 1.16±0.03for NH 3 , with a nuclear spin temperature of 28±1 K.The ammonia OPR was found equal to 1.19±0.03,corresponding to a spin temperature of 27±1 K. Using very-high-resolution ESPRESSO spectra, we measured a forbidden oxygen lines ratio G/R of 0.23±0.02.These measurements are in agreement with those found for several comets of different dynamical types and origins and do not reveal any peculiarity with regard to the hyperactive comet 46P.

Fig. 2 .
Fig. 2. A(0)fρ parameter from the broad and narrow-band filters as a function of time to perihelion in days.Triangle symbols indicate data acquired with TN.Filled circles indicate TS.

Fig. 3 .
Fig. 3. Logarithm of the production rate and dust-to-gas ratios of comet 46P with respect to CN (Log10[Q(xx)/Q(CN)]) and to OH (Log10[Q(xx)/Q(OH)]) as a function of time to perihelion.Pre-perihelion values are represented with open circles and post-perihelion values with filled circles.The dust-to-gas ratio is expressed in cm second per molecule.The horizontal dashed line indicates the mean ratios of typical comets given in Schleicher (2008).

Fig. 4 .
Fig. 4. Water-production rate of comet 46P as a function of time from perihelion during the 2018 apparition.Our measurements and those from Knight et al. (2021) were derived from the OH production rates.Combi et al. (2020) measurements were derived from H Lyman-α spectra using the SWAN instrument on SOHO.iSHELL/IRTF and (McKay et al. 2021) data were derived from NIR spectra.Lis et al. (2019) data were derived from far-infrared spectra using the GREAT spectrometer aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA).

Fig. 5 .
Fig. 5. Water-production rates of comet 46P as a function of time from perihelion during 1991 (orange), 1997 (black), 2008 (blue), and 2018 (in red) apparitions.Fink et al. (1998) derived the water-production rate from the [OI]( 1 D) forbidden line doublet emission in the optical, Bertaux et al. (1999) values were derived from the H Ly-α emission observed by the SWAN instrument on board SOHO and Crovisier et al. (2002) values were derived from the OH emission observed at the radio wavelengths.

Fig. 6 .
Fig. 6.Phased CN light curves for the nucleus rotation period of comet 46P at different epochs.Time origin for the phase is JD 2458440 on Nov. 17, 2018.Data combined for (a) Nov. 23, 28 and Dec. 1, 2018 preperihelion, (b) Dec. 7-10, 2018 at perihelion and for (c) Jan. 12 and 15, 2019 post-perihelion.The red curve shows the best fit of the data.

Fig. 7 .
Fig. 7. Evolution of the CN jets of comet 46P along its orbit around the Sun.The images are oriented north (N) up and east (E) right.

Fig. 9 .
Fig. 9. NH 3 T spin in comets of various dynamical origin.The red filled circle is the value of 46P obtained with UVES at VLT (this work).

Fig. 10 .
Fig. 10.Three forbidden oxygen [OI] lines in comet 46P, from ESPRESSO observations, at 557.73 nm (left) for the green line, and 630.03 nm (middle) and 636.38 nm (right) for the red doublet.The cometary oxygen lines are identified by tick marks and the telluric lines marked by dashed tick.The y-axis is in arbitrary units for all three panels.The red line is a fit for one of the six spectra taken in the centered fiber.

Table 1 .
Mean relative molecular abundances of comet 46P for the perihelion 1991 to 2018 compared to typical and depleted comets.
Notes.The Afρ/Q values for typical and depleted comets are from A'Hearn et al.

Table 2 .
Abundances ratios derived from the optical and infrared data of comet 46P in its 2018 apparition.

Table 3 .
Rotation period of 46P nucleus at different epochs.

Table 4 ,
and they have average values of 3.31 ± 0.03 and 1.19 ± 0.02, respectively.A nuclear spin temperature (T spin ) for ammonia of 28 ± 1 K was derived.The 46P value is consistent with typical values measured in comets (see Fig.9), which is a possible cosmogonic indicator linked to the formation temperature of the molecule.

Table 4 .
(Pepe et al. 2021)H 3 OPRs in comet 46P.Upper and lower values for each band are the results for the two spectra taken at UT 0:59:52 and UT 1:19:54 on Dec. 9, respectively.Using the very high-resolution ESPRESSO spectra(Pepe et al. 2021), we were able to resolve the telluric and cometary forbidden oxygen lines even with a Doppler shift as small as 5 km/s.We detected the cometary [OI] lines at 557.73, 630.03, and 636.38 nm in all six spectra (See Figure