Coma environment of comet C/2017 K2 around the water ice sublimation boundary observed with VLT/MUSE

We report a new imaging spectroscopic observation of Oort-cloud comet C/2017 K2 (hereafter K2) on its way to perihelion at 2.53 au, around a heliocentric distance where H2O ice begins to play a key role in comet activation. Normalized reflectances over 6 500--8 500 AA for its inner and outer comae are 9.7+/-0.5 and 7.2+/-0.3 % (10^3 AA)^-1, respectively, the latter being consistent with the slope observed when the comet was beyond the orbit of Saturn. The dust coma at the time of observation appears to contain three distinct populations: mm-sized chunks prevailing at <~10^3 km; a 10^5-km steady-state dust envelope; and fresh anti-sunward jet particles. the dust chunks dominate the continuum signal and are distributed over a similar radial distance scale as the coma region with redder dust than nearby. they also appear to be co-spatial with OI1D, suggesting that the chunks may accommodate H2O ice with a fraction (>~1 %) of refractory materials. The jet particles do not colocate with any gas species detected. The outer coma spectrum contains three significant emissions from C2(0,0) Swan band, OI1D, and CN(1,0 red band, with an overall deficiency in NH2. Assuming that all OI1D flux results from H2O dissociation, we compute an upper limit on the water production rate Q_H2O of ~7 x 10^28 molec s^-1 (with an uncertainty of a factor of two). the production ratio log[Q_C2/Q_CN] of K2 suggests that the comet has typical carbon-chain composition, with the value potentially changing with distance from the Sun. Our observations suggest that water ice-containing dust chunks (>0.1 mm) near K2's nucleus emitted beyond 4 au may be responsible for its very low gas rotational temperature and the discrepancy between its optical and infrared lights reported at similar heliocentric distances.Comment: Accepted for publication in Astronomy & Astrophysic

TRAPPIST monitoring of comet C/2012 F6 (Lemmon)

Comet C/2012 F6 is a long-period comet that reached perihelion on March 23, 2012. The unexpected brightness of this comet since December 2012 allowed us to obtain narrowband photometry and to study its chemical composition as well as its rotation

Comets 12^{12}CO+^+ and 13^{13}CO+^+ fluorescence models for measuring the 12^{12}C/13^{13}C isotopic ratio in CO+^+

Context: CO is an abundant species in comets, creating CO$^+$ ion with emission lines that can be observed in the optical spectral range. A good modeling of its fluorescence spectrum is important for a better measurement of the CO$^+$ abundance. Such a species, if abundant enough, can also be used to measure the $^{12}$C/$^{13}$C isotopic ratio. Aims: This study uses the opportunity of a high CO content observed in the comet C/2016 R2 (PanSTARRS), that created bright CO$^{+}$ emission lines in the optical range, to build and test a new fluorescence model of this species and to measure for the first time the $^{12}$C/$^{13}$C isotopic ratio in this chemical species with ground-based observations. Methods: Thanks to laboratory data and theoretical works available in the scientific literature we developed a new fluorescence model both for $^{12}$CO$^+$ and $^{13}$CO$^+$ ions. The $^{13}$CO$^+$ model can be used for coadding faint emission lines and obtain a sufficient signal-to-noise ratio to detect this isotopologue. Results: Our fluorescence model provides a good modeling of the $^{12}$CO$^+$ emission lines, allowing to publish revised fluorescence efficiencies. Based on similar transition probabilities for $^{12}$CO$^+$ and $^{13}$CO$^+$ we derive a $^{12}$C/$^{13}$C isotopic ratio of 73$\pm$20 for CO$^+$ in comet C/2016 R2. This value is in agreement with the solar system ratio of 89$\pm$2 within the error bars, making the possibility that this comet was an interstellar object unlikely.Comment: 11 pages, 8 figure

Isotopic ratios in outbursting comet C/2015 ER61

Isotopic ratios in comets are critical to understanding the origin of cometary material and the physical and chemical conditions in the early solar nebula. Comet C/2015 ER61 (PANSTARRS) underwent an outburst with a total brightness increase of 2 magnitudes on the night of 2017 April 4. The sharp increase in brightness offered a rare opportunity to measure the isotopic ratios of the light elements in the coma of this comet. We obtained two high-resolution spectra of C/2015 ER61 with UVES/VLT on the nights of 2017 April 13 and 17. At the time of our observations, the comet was fading gradually following the outburst. We measured the nitrogen and carbon isotopic ratios from the CN violet (0,0) band and found that $^{12}$C/$^{13}$C=100 $\pm$ 15, $^{14}$N/$^{15}$N=130 $\pm$ 15. In addition, we determined the $^{14}$N/$^{15}$N ratio from four pairs of NH$_2$ isotopolog lines and measured $^{14}$N/$^{15}$N=140 $\pm$ 28. The measured isotopic ratios of C/2015 ER61 do not deviate significantly from those of other comets.Comment: 4 pages, 4 figures, accepted to be published by A&

An inversion method for cometary atmospheres

Remote observation of cometary atmospheres produces a measurement of the cometary emissions integrated along the line of sight. This integration is the so-called Abel transform of the local emission rate. The observation is generally interpreted under the hypothesis of spherical symmetry of the coma. Under that hypothesis, the Abel transform can be inverted. We derive a numerical inversion method adapted to cometary atmospheres using both analytical results and least squares fitting techniques. This method, derived under the usual hypothesis of spherical symmetry, allows us to retrieve the radial distribution of the emission rate of any unabsorbed emission, which is the fundamental, physically meaningful quantity governing the observation. A Tikhonov regularization technique is also applied to reduce the possibly deleterious effects of the noise present in the observation and to warrant that the problem remains well posed. Standard error propagation techniques are included in order to estimate the uncertainties affecting the retrieved emission rate. Several theoretical tests of the inversion techniques are carried out to show its validity and robustness. In particular, we show that the Abel inversion of real data is only weakly sensitive to an offset applied to the input flux, which implies that the method, applied to the study of a cometary atmosphere, is only weakly dependent on uncertainties on the sky background which has to be subtracted from the raw observations of the coma. We apply the method to observations of three different comets observed using the TRAPPIST telescope: 103P/ Hartley 2, F6/ Lemmon and A1/ Siding Spring. We show that the method retrieves realistic emission rates, and that characteristic lengths and production rates can be derived from the emission rate for both CN and C2 molecules. We show that the retrieved characteristic lengths can differ from those obtained from a direct least squares fitting over the observed flux of radiation, and that discrepancies can be reconciled for by correcting this flux by an offset (to which the inverse Abel transform is nearly not sensitive). The A1/Siding Spring observations were obtained very shortly after the comet produced an outburst, and we show that the emission rate derived from the observed flux of CN emission at 387 nm and from the C2 emission at 514.1 nm both present an easily-identifiable shoulder that corresponds to the separation between pre- and post-outburst gas. As a general result, we show that diagnosing properties and features of the coma using the emission rate is easier than directly using the observed flux, because the Abel transform produces a smoothing that blurs the signatures left by features present in the coma. We also determine the parameters of a Haser model fitting the inverted data and fitting the line-of-sight integrated observation, for which we provide the exact analytical expression of the line-of-sight integration of the Haser model

A physico-chemical model to study the ion density distribution in the inner coma of comet C/2016 R2 (Pan-STARRS)

peer reviewedThe recent observations show that comet C/2016 R2 (Pan-Starrs) has a unique and peculiar composition when compared with several other comets observed at 2.8 au heliocentric distance. Assuming solar resonance fluorescence is the only excitation source, the observed ionic emission intensity ratios are used to constrain the corresponding neutral abundances in this comet. We developed a physico-chemical model to study the ion density distribution in the inner coma of this comet by accounting for photon and electron impact ionization of neutrals, charge exchange and proton transfer reactions between ions and neutrals, and electron-ion thermal recombination reactions. Our calculations show that CO[SUB]2[/SUB][SUP]+[/SUP] and CO[SUP]+[/SUP] are the major ions in the inner coma, and close to the surface of nucleus CH[SUB]3[/SUB]OH[SUP]+[/SUP], CH[SUB]3[/SUB]OH[SUB]2[/SUB][SUP]+[/SUP], and O[SUB]2[/SUB][SUP]+[/SUP] are also important ions. By considering various excitation sources, we also studied the emission mechanisms of different excited states of CO[SUP]+[/SUP], CO[SUB]2[/SUB][SUP]+[/SUP], N[SUB]2[/SUB][SUP]+[/SUP], and H[SUB]2[/SUB]O[SUP]+[/SUP]. We found that the photon and electron impact ionization and excitation of corresponding neutrals significantly contribute to the observed ionic emissions for radial distances smaller than 300 km and at larger distances, solar resonance fluorescence is the major excitation source. Our modelled ion emission intensity ratios are consistent with the ground-based observations. Based on the modelled emission processes, we suggest that the observed ion emission intensity ratios can be used to derive the neutral composition in the cometary coma only when the ion densities are significantly controlled by photon and photoelectron impact ionization of neutrals rather than by the ion-neutral chemistry

Cometary science with CUBES

The proposed CUBES spectrograph for ESO's Very Large Telescope will be an exceptionally powerful instrument for the study of comets. The gas coma of a comet contains a large number of emission features in the near-UV range covered by CUBES (305-400 nm), which are diagnostic of the composition of the ices in its nucleus and the chemistry in the coma. Production rates and relative ratios between different species reveal how much ice is present and inform models of the conditions in the early solar system. In particular, CUBES will lead to advances in detection of water from very faint comets, revealing how much ice may be hidden in the main asteroid belt, and in measuring isotopic and molecular composition ratios in a much wider range of comets than currently possible, provide constraints on their formation temperatures. CUBES will also be sensitive to emissions from gaseous metals (e.g., FeI and NiI), which have recently been identified in comets and offer an entirely new area of investigation to understand these enigmatic objects.Comment: Accepted for publication in Experimental Astronom

The N2 Production Rate in C/2016 R2 (PanSTARRS)

Introduction: Radio observations of long-period comet C/2016 R2 (PanSTARRS) revealed that it was remarkably depleted in water (Biver et al. 2018). The spectrum was instead dominated by bands of CO+ and N2+, rarely seen in such abundance in comets before (Cochran & Mckay 2018). Understanding the nature of this comet would allow us to investigate key features in the timeline of planetesimal formation.By measuring the observed emission fluxes of the observed N2+ in C/2016 R2's spectrum, ionic ratios of N2+/CO+ in the coma were estimated to be between 0.06 (Cochran & Mckay. 2018), Opitom et al. 2019) and 0.08 (Biver et al. 2018). This would be the same ratio for N2/CO since ionization efficiencies of N2 and CO are similar at 1 au for quiet Sun (Huebner et al. 1992).C/2016 R2 provides a unique opportunity to set a baseline for identifying N2 in cometary spectra. By using the Ultraviolet-Visual Echelle Spectrograph (UVES) mounted on the 8.2 m UT2 telescope of the European Southern Observatory Very Large Telescope (ESO VLT) observations, we can constrain the properties of N2 in the cometary coma and establish new Haser scalelengths in order to determine the N2 production rate, which we present here. Observations: The observations of C/2016 R2 used in our work were collected on 2018 February 11, 13, and 14 with UVES. All observations were made when the comet was near its perihelion distance of 2.6 au, at 2.76 and 2.75 au. A full description of the observations and data reduction can be found in Opitom et al. (2019). Methods: We aim to fit the observed flux with a Haser profile (Haser (1957)), providing an analytical solution to the column density of parent- and daughter-species in the coma along the line of sight. N2+ being an ion, the Haser model will be restricted to an area near the coma. The UVES slit covers ~6500 km on either side of the nucleus, a narrow region in which ions should be undisturbed by the solar wind. CN scalelengths and Production Rate: We first fit a Haser profile on the CN emissions to ensure scalelengths can properly be determined from our data. We created a synthetic CN model evaluated by interpolation from a spectrum calculated by Zucconi (1985). This model is then convolved by the response of our instrument, with an FWHM of 0.06 Å. For each night of 11, 13, and 14 Feb, the CN lines are identified and summed along the spectroscopic slit. The total flux measured for CN over the entire spectrograph and averaged over the three nights of observation was 2.1x10-15 erg/s/cm2. The flux intensities are then averaged again over their cometocentric distances so as to allow for a proper fit of the Haser model.By using a X2 test, we estimate the best fit of the Haser model to the observed intensity profile and determine the scalelengths of both the parent- (HCN) and daughter- (CN) species in the coma of C/2016 R2. We found lp = 1.3 x 104 km and ld = 2.8 x 105 km (scaled to 1 au using an rh2 law) as shown on Fig. 1. With g = 3.52 x 10-2 photons/s/molecule at 1 au (Schleicher et al. 2010), we estimate a production rate of Q(CN) = (9.8±0.5) x 1024 mol/s. N2+ scalelengths and Production Rate: The production rate was estimated via relative ratios with g =7 x 10-2 photons/ion/s from Lutz et al. (1993) by Wierzchos (2018) as Q(N2) = (2.8 ±0.4) x 1027 mol/s and by McKay (2019) as Q(N2) = (4.8 ± 1.1)x1027 mol/s. It can be inferred from Biver et al. (2018) to be ~8.5 x 1027 mol/s for a Q(CO) = 1.1 x 1029 mol/s. These results are first re-calculated with the most recent g factor from Rousselot et al. 2022. With g =4.90 x 10-3 photons/mol/s at 1 au, prior measurements of the N2+ production rates become Q(N2) = 4.6x1027 mol/s (Wierzchos & M. Womack 2018), =8.0x1027 mol/s (McKay et al. 2019), and 1.4x1028 mol/s (Biver et al. 2018).We limit the identification process to the 3885.5 Å to 3915.0 Å interval to further avoid contamination by the CN emission lines. We explore this interval with the X2 test and find new scalelengths of lp = 2.8 x 106 km and ld = 3.8 x 106 km scaled to 1 au (see Fig. 1). These values are within the expected range estimated from the rate coefficients. However, at this scale, multiple pairs of scalelengths could be selected for N2+ with an equally good fit. We thus have a large uncertainty on the production rate.Using g = 5.41 x 10-3 photons/mol/s (at rh) for the (0,0) band between 3885.5-3915.0 Å and FTOT = 1.0 x 10-14 erg/s/cm2, we find Q(N2)=(8 ±1) x 1027. With Q(CO) ~ 1.1 × 1029 molecules.s-1, N2/CO = 0.07, consistent with observed intensity ratios.Figure 1: The best fit of the Haser model for CN (top, purple, compared to other fits using scalelengths from literature) and N2+ (bottom, blue)

Optical Spectropolarimetry of Binary Asteroid Didymos–Dimorphos before and after the DART Impact

We have monitored the Didymos-Dimorphos binary asteroid in spectropolarimetric mode in the optical range before and after the DART impact. The ultimate goal was to obtain constraints on the characteristics of the ejected dust for modeling purposes. Before impact, Didymos exhibited a linear polarization rapidly increasing with phase angle, reaching a level of similar to 5% in the blue and similar to 4.5% in the red. The shape of the polarization spectrum was anticorrelated with that of its reflectance spectrum, which appeared typical of an S-class asteroid. After impact, the level of polarization dropped by about 1 percentage point (pp) in the blue band and about 0.5 pp in the red band, then continued to linearly increase with phase angle, with a slope similar to that measured prior to impact. The polarization spectra, once normalized by their values at an arbitrary wavelength, show very little or no change over the course of all observations before and after impact. The lack of any remarkable change in the shape of the polarization spectrum after impact suggests that the way in which polarization varies with wavelength depends on the composition of the scattering material, rather than on its structure, be this a surface or a debris cloud.Peer reviewe

Modeling of N2+ and 14N15N+ fluorescence spectrum in comets

1. IntroductionC/2016 R2 (PanSTARRS) was a surprising comet. Detected on September 7, 2016 by Pan-STARRS it showed an unusual composition when it became a bright comet at the end of 2017 and the beginning of 2018. It developed a coma at large (~6 au) heliocentric distance and observations showed that it had a highly unusual composition: no water molecules (or OH radical) could be detected, and the abundances of the usual radicals (CN, C2, C3) were unusually low, with a surprising coma composition dominated by CO, CO2 and N2 molecules with bright CO+ and N2+ emission lines in the visible range. A high CO production rate of about 1029 molecules s-1 was measured (Biver et al. 2018; Wierzchos & Womack 2018) as well as a high CO2 production rate (CO2/CO=1.1 from Opitom et al. 2019), and a high ratio N2/CO varying between 0.06 and 0.09 (Biver et al. 2018; Cochran & McKay 2018a,b; Opitom et al. 2019; Venkataramani et al. 2020).The detection of such bright N2+ emission lines in this comet highlighted the necessity of a good modeling of the N2+ fluorescence spectrum in comets. The high-quality spectra published by Opitom et al. (2019) provided a good opportunity to test such a model. This model also permits to compute the fluorescence spectrum of the 14N15N+ species, leading to the possibility of future measurements of the 14N/15N isotopic ratio in the N2 molecules, one of the main constituant of the solar nebula.2. ObservationsThe spectra used for this work have been obtained with the UVES spectrograph mounted on the ESO 8.2 m UT2 telescope of the VLT. Three different observing nights have been used, corresponding to February 11, 13 and 14, 2018. One single exposure of 4800 s of integration time was obtained during each night and we used a 0.44" wide slit, providing a resolving power R~80,000. The slit length was 8" corresponding to about 14,500 km at the distance of the comet (geocentric distance of 2.4 au). The average heliocentric distance was 2.76 au. Opitom et al. (2019) describe in more details the data processing.From the 2D spectra having a spatial extension of 30 rows, each of them corresponding to a different cometocentric distance, we extracted different 1D spectra for each night. These spectra were then averaged for similar cometocentric distances allowing a detailed comparison of these spectra at different cometocentric distances, the furthest one corresponding to 2x4 rows at the two extremities of the slit (i.e. at a cometocentric distance varying between 4800 and 6600 km).3. Modeling the N2+ fluorescence spectrumWe developed a new fluorescence model for modeling our observational spectra. The transition involved in this spectrum is the first negative group, i.e. the B2Σu+ → X2+Σg+ electronic transition with the (0,0) bandhead appearing near 3914 Å. We considered the first three vibrational levels (v = 0; 1; 2) for both X2+Σg+ and B2Σu+ state, each of them with all the rotational levels from N = 0 to 40.N2+ having no permanent dipole moment, the pure rotational and vibrational transitions are forbidden (or have a very low probability, through quadrupolar transitions, not taken into account in our model). For that reason it takes a long time for this species to reach its fluorescence equilibrium because it needs a few tens of absorption / emission cycles between the X2+Σg+ and B2Σu+ states to reach this equilibrium. A comparison of the spectrum obtained on the nucleus with the one obtained at the edges of the slit revealed clear differences due to different rotational relative populations. For that reason we decided to model the N2+ fluorescence spectrum with a Monte-Carlo simulation. Such a computational method allows to compute a spectrum at different times from an initial relative population distribution. Our model starts with a Boltzmann relative population distribution of 80 K (representing an estimate of the kinetic temperature in the inner coma) and uses 10,000 s of evolution time.We managed to explain satisfactorily the observed N2+ emission spectrum. Fig. 1 presents a close up view around the (0,0) bandhead. This work, presented in more details in Rousselot et al. (2022) also allowed to compute accurate fluorescence efficiencies. Figure 1: Comparison of the observed VLT UVES spectrum of comet C/2016 R2 (blue) obtained at the ends of the slit with our N2+ model (red).4. 14N15N+ fluorescence spectrumOur modeling of the N2+ fluorescence spectrum can be used to compute the 14N15N+ fluorescence spectrum, leading to the possibility of measuring the 14N/15N isotopic ratio in N2 molecules. We will present such a spectrum as well as a search for this isotopologue in the C/2016 R2 spectra. Such comets are rare but future observations will reveal other comets similar in composition to C/2016 R2. With future observing facilities now under construction (such as the ESO ELT) 14N/15N measurements for N2 molecules will probably become possible, leading to new constraints on this isotopic ratio