# Hint of an exocomet transit in the CHEOPS light curve of HD 172555

F. Kiefer, V. van Grootel, A. Lecavelier des Etangs, Gy. M. Szabó, A.

Brandeker, C. Broeg, A. Collier Cameron, A. Deline, G. Olofsson, T. G.

Wilson, et al.

► To cite this version:

F. Kiefer, V. van Grootel, A. Lecavelier des Etangs, Gy. M. Szabó, A. Brandeker, et al. Hint of an exocomet transit in the CHEOPS light curve of HD 172555. Astronomy and Astrophysics - A&A, 2023, 671, 10.1051/0004-6361/202245104. insu-04155765

HAL Id: insu-04155765

https://insu.hal.science/insu-04155765

Submitted on 7 Jul 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Hint of an exocomet transit in the CHEOPS light curve of HD 172555*

F. Kiefer¹,², V. Van Grootel³, A. Lecavelier des Etangs⁴, Gy. M. Szabó⁵,⁶, A. Brandeker⁷, C. Broeg⁸,⁹, A. Collier Cameron¹⁰, A. Deline¹¹, G. Olofsson⁷, T. G. Wilson¹⁰, S. G. Sousa¹², D. Gandolfi¹³, G. Hébrard¹,¹⁴, Y. Alibert⁸, R. Alonso¹⁵,¹⁶, G. Anglada¹⁷,¹⁸, T. Bárczy¹⁹, D. Barrado²⁰, S. C. C. Barros¹²,²¹, W. Baumjohann²², M. Beck¹¹, T. Beck⁸, W. Benz⁸,²³, N. Billot¹¹, X. Bonfils²⁴, J. Cabrera²⁵, S. Charnoz²⁶, Sz. Csizmadia²⁵, M. B. Davies²⁷, M. Deleuil²⁸, L. Delrez²⁹,³, O. D. S. Demangeon¹²,²¹, B.-O. Demory²³, D. Ehrenreich³⁰,³¹, A. Erikson²⁵, A. Fortier⁸,²³, L. Fossati²², M. Fridlund³²,³³, M. Gillon²⁹, M. Güdel³⁴, K. Heng²³,³⁵, S. Hoyer²⁸, K. G. Isaak³⁶, L. L. Kiss³⁷,³⁸, J. Laskar³⁹, M. Lendl¹¹, C. Lovis¹¹, D. Magrin⁴⁰, P. F. L. Maxted⁴¹, M. Munari⁴², V. Nascimbeni⁴⁰, R. Ottensamer⁴³, I. Pagano⁴², E. Pallé¹⁵, G. Peter⁴⁴, D. Piazza⁴⁵, G. Piotto⁴⁰,⁴⁶, D. Pollacco³⁵, D. Queloz⁴⁷,⁴⁸, R. Ragazzoni⁴⁰,⁴⁶, N. Rando⁴⁹, F. Ratti⁴⁹, H. Rauer²⁵,⁵⁰,⁵¹, C. Reimers⁴³, I. Ribas¹⁷,¹⁸, N. C. Santos¹²,²¹, G. Scandariato⁴², D. Ségransan¹¹, A. E. Simon⁸, A. M. S. Smith²⁵, M. Steller²², N. Thomas⁸, S. Udry¹¹, I. Walter⁵², and N. A. Walton⁵³

*(Affiliations can be found after the references)*

Received 27 September 2022 / Accepted 10 December 2022

## ABSTRACT

HD 172555 is a young (~20 Myr) A7V star surrounded by a 10 au wide debris disk suspected to be replenished partly by collisions between large planetesimals. Small evaporating transiting bodies, that is exocomets, have also been detected in this system by spectroscopy. After β Pictoris, this is another example of a system possibly witnessing a phase of the heavy bombardment of planetesimals. In such a system, small bodies trace dynamical evolution processes. We aim to constrain their dust content by using transit photometry. We performed a 2-day-long photometric monitoring of HD 172555 with the CHEOPS space telescope in order to detect shallow transits of exocomets with a typical expected duration of a few hours. The large oscillations in the light curve indicate that HD 172555 is a δ Scuti pulsating star. After removing those dominating oscillations, we found a hint of a transient absorption. If fitted with an exocomet transit model, it would correspond to an evaporating body passing near the star at a distance of 6.8±1.4 *R<sub>⋆</sub>* (or 0.05±0.01 au) with a radius of 2.5 km. These properties are comparable to those of the exocomets already found in this system using spectroscopy, as well as those found in the β Pic system. The nuclei of the Solar System's Jupiter family comets, with radii of 2–6 km, are also comparable in size. This is the first piece of evidence of an exocomet photometric transit detection in the young system of HD 172555.

**Key words.** stars: variables: δ Scuti – techniques: photometric – comets: general – circumstellar matter – stars: individual: HD 172555

## 1. Introduction

In the direct aftermath of planet formation, planetary systems can carry a variety of bodies with many sizes and masses: gas and dust particles within disks, small bodies that are more or less icy (asteroids and comets), and planets. How those appear and evolve together – or one after the other – is still an open question. Dynamical instabilities leading to planet orbit reorganisation as well as small body scattering, during or after protoplanetary disk dissipation, are thought to take place early in the life of planetary systems, and they are believed to have occurred in the young Solar System (Morbidelli et al. 2001; Chambers & Lissauer 2002; Tsiganis et al. 2005; Gomes et al. 2005; de Sousa et al. 2020; Liu et al. 2022). In the Solar System, cratering impacts on the Moon are the result of the bombardment of small bodies that occurred several billion years ago (Ryder 1990; Bottke et al. 2012; Morbidelli et al. 2018; Marks et al. 2019; Cartwright et al. 2022). How this phase relates to the configuration and properties

of the individual planets in the Solar System is still under investigation, primarily because the chronology is uncertain due to the possibility of observing the Solar System only now, several billion years after these events. Observing young systems, such as HD 172555 and β Pictoris, that may currently be undergoing a period of heavy bombardment provides a unique window into the dynamical processes at play during the first dozen million years of a planetary system.

The naked-eye star ($V = 4.9$) HD 172555 is a 20 Myr-old A7V dwarf belonging to the β Pic moving group, whose age is well established to be within 15–25 Myr (Barrado y Navascués et al. 1999; Zuckerman et al. 2001; Binks & Jeffries 2014; Mamajek & Bell 2014; Miret-Roig et al. 2018). Its principal properties are summarised in Table 1. It is the member of a wide binary system with a distant M-dwarf companion, CD-64 1208, at ~2000 au, or equivalently 70" (Feigelson et al. 2006; Alonso-Floriano et al. 2015). HD 172555 harbors a cold dust disk with an extent of up to ~1000 au (Nilsson et al. 2009), and a warm 10 au wide dusty debris disk close to edge-on at the inclination $I_{\text{disk}}$ ~ 75° (Smith et al. 2012; Engler et al. 2018) with strong and unexpected

<footnote>* This article uses data from CHEOPS programme CH_PR100010.</footnote>
<table><caption>**Table 1.** Main properties of HD 172555.</caption><thead><tr><th>Parameter</th><th>Unit</th><th>Value</th><th>References</th></tr></thead><tbody><tr><td>RA</td><td>hh:mm:ss</td><td>18:45:26.9019</td><td>1</td></tr><tr><td>Dec</td><td>dd:mm:ss</td><td>-64:52:16.5417</td><td>1</td></tr><tr><td><i>m</i><sub>V</sub></td><td>mag</td><td>4.77</td><td>2</td></tr><tr><td><i>B-V</i></td><td>mag</td><td>0.20</td><td>2</td></tr><tr><td><i>d</i></td><td>pc</td><td>28.8</td><td>1</td></tr><tr><td><i>T</i><sub>eff</sub></td><td>K</td><td>7800</td><td>3,4</td></tr><tr><td>log <i>g</i></td><td>dex</td><td>4.3</td><td>3,4</td></tr><tr><td>[Fe/H]</td><td>dex</td><td>0.07</td><td>3</td></tr><tr><td>Sp. type</td><td></td><td>A7V</td><td>5</td></tr><tr><td>Age</td><td>Myr</td><td>15–25</td><td>6–10</td></tr><tr><td><i>M</i><sub>★</sub></td><td><i>M</i><sub>☉</sub></td><td>1.7</td><td>4</td></tr><tr><td><i>R</i><sub>★</sub></td><td><i>R</i><sub>☉</sub></td><td>1.55</td><td>4, 11<sup>(†)</sup></td></tr><tr><td><i>L</i><sub>★</sub></td><td><i>L</i><sub>☉</sub></td><td>7.8</td><td>4</td></tr></tbody></table>

**References.** (1) Gaia Collaboration (2020); (2) Høg et al. (2000); (3) Erspamer & North (2003); (4) Riviere-Marichalar et al. (2012); (5) Gray et al. (2006); (6) Barrado y Navascués et al. (1999); (7) Zuckerman et al. (2001); (8) Binks & Jeffries 2014; (9) Mamajek & Bell (2014); (10) Miret-Roig et al. (2018); (11) this work. (†) We used a modified version of the infrared flux method (IRFM; Blackwell & Hallis 1977; Schanche et al. 2020), measuring a stellar radius of $R_{\star} = 1.550 \pm 0.012 R_{\odot}$. It agrees with the $R_{\star} = 1.56 R_{\odot}$ derived in Riviere-Marichalar et al. (2012). See Bonfanti et al. (2021) and Wilson et al. (2022) for details and application of the IRFM.

traces of both SiO and CO (Lisse et al. 2009; Schneiderman et al. 2021). Large amounts of oxygen have also been found at a larger distance (Riviere-Marichalar et al. 2012). The main explanation for the presence of gas is a high-speed collision between planetsimals, together with the presence of evaporating icy bodies liberating a mixture of CO, CO$_2$, and H$_2$O out of which oxygen is produced.

Exocomets were discovered in the HD 172555 system through β Pic-like varying spectroscopic signatures of cometary transits in the Ca II doublet of the star in the optical band observed at high resolution (Kiefer et al. 2014a; K14, hereafter). Later, they were also observed in the UV in the CII-CII* lines (Grady et al. 2018). The transiting clouds are made of gas evaporated from inner nuclei on eccentric orbits that cross the star's line of sight with certain transverse speeds and radial velocities (Beust et al. 1990). The main observed events derived from an exomet com transit and already witnessed in a few young systems have the following properties: (i) they are variable, possibly deep and Doppler-shifted, spectroscopic absorptions features due to an extended cloud of sublimed ions covering the solid angle of the stellar surface (Ferlet et al. 1987; Beust et al. 1990; Lagrange-Henri et al. 1992; Vidal-Madjar et al. 1994; Kiefer et al. 2014b); and (ii) they experience photometric transits due to the passage of the head and tail of the dust coma in front of the star (Lecavelier Des Etangs et al. 1999, 2022; Kiefer et al. 2017; Rappaport et al. 2018; Zieba et al. 2019).

Ions that have been observed using the technique of transit spectroscopy place indirect constraints on the dust from which they have sublimated and further ionised (Beust & Tagger 1993; Kiefer et al. 2014b). Photometry has the advantage over spectroscopy in being able to achieve a direct characterisation of the dust evaporating from the cometary nuclei (Lecavelier Des Etangs et al. 1999). The detection of 30 exomet com transits in the Transiting Exomet com Survey Satellite (TESS) photometric light curves of β Pic allowed for the first statistical comparison of the

<table><caption>**Table 2.** Log of CHEOPS observations.</caption><thead><tr><th>Parameter</th><th>Unit</th><th>Value</th></tr></thead><tbody><tr><td>DRP Version</td><td></td><td>13.1.0</td></tr><tr><td>Program ID</td><td></td><td>CHEOPS-10</td></tr><tr><td>Program PI</td><td></td><td>Gyula Szabó</td></tr><tr><td>Obs. ID</td><td></td><td>1514950</td></tr><tr><td>CHEOPS orbit</td><td>h</td><td>1.645</td></tr><tr><td>Aperture radius</td><td>"</td><td>25</td></tr><tr><td>Source magnitude</td><td>CHEOPS band</td><td>4.7</td></tr><tr><td>RA</td><td>hh:mm:ss</td><td>18:45:26.9019</td></tr><tr><td>Dec</td><td>dd:mm:ss</td><td>-64:52:16.5417</td></tr><tr><td>Date start</td><td>MJD</td><td>59 381.445</td></tr><tr><td>Date end</td><td>MJD</td><td>59 383.501</td></tr><tr><td>Total duration</td><td>d</td><td>2.06</td></tr><tr><td>Exposure time</td><td>s</td><td>49.3</td></tr><tr><td>Exposures stacking order</td><td></td><td>24</td></tr><tr><td>Total integration time</td><td>s</td><td>16.8</td></tr><tr><td>Number of frames</td><td></td><td>2272</td></tr><tr><td>Number of flagged frames</td><td></td><td>122</td></tr><tr><td>point-to-point RMS</td><td>ppm</td><td>1581</td></tr><tr><td>CDPP(10 min)</td><td>ppm</td><td>180</td></tr><tr><td>CDPP(3 h)</td><td>ppm</td><td>101</td></tr></tbody></table>

exocomets' population in a young systems with the Solar System comets (Lecavelier Des Etangs et al. 2022). It showed that the inner nuclei of Jupiter family comets (JFC, hereafter) and β Pic's exocomets follow a similar radius distribution with radii 1–10 km and a $dN \propto R^{-3.6 \pm 0.8} dR$. Collisional relaxation in a population of small bodies can produce such a distribution. Our objective is to directly observe the dust produced by the HD 172555 exocomets, and enable the characterisation of the small bodies in this young system.

In this paper, we report on high-precision photometric monitoring of HD 172555 with the CHaracterising ExOPlanets Satellite (CHEOPS; Benz 2021). In Sect. 2 we summarise the instrumental configuration and the CHEOPS data reduction. In Sect. 3 we analyse the δ Scuti variations observed in HD 172555's CHEOPS light curve. In Sect. 4 we report the search for transients in the light curve detrended from the identified δ Scuti oscillations. We compare our result to the K14 exocomets in Sect. 5. We give our conclusions in Sect. 6.

## 2. CHEOPS observations

### 2.1. *Observations' settings and data reduction*

HD 172555 was observed within the frame of the CHEOPS Dusty Debris Disk Guaranteed Time Observation (DDD-GTO) program. It was observed for 30 consecutive CHEOPS orbits (~1.645 h each), for a total duration of 2.1 days. Table 2 gathers details about these CHEOPS observations. The exposure time was set to 0.7 s, leading to a cadence, including readout time, of ~43.2 s. The sub-array frames were automatically processed with the CHEOPS Data Reduction Pipeline (DRP; Hoyer et al. 2020), including smearing, cosmic ray hits, background and stray-light corrections, and finally default aperture photometry extraction.

The CHEOPS point spread function (PSF) is relatively extended because of the defocus design of the instrument. The
<figure><img src="image_2.png" /></figure>

**Fig. 1.** Field of view around HD 172555 showing only background stars with their CHEOPS triangular PSF. East is downward and north is leftward. Gaia identified objects are indicated as red circles and the large black circle represents the optimal aperture of CHEOPS with a size of 25". Apart from HD 172555, the three closest objects within the aperture have *G* magnitudes of 18.5, 15, and 20 in order of increasing distance from the centre of the aperture. The brightest object in the east is the M-dwarf companion CD-64 1208 with *G* = 8.9 mag.

<figure><img src="image_3.png" /></figure>

**Fig. 2.** Flux of the simulated contaminant stars in the aperture (black points) as a function of the roll angle (bottom panel) relative to the flux of the target (in percent).

default aperture is defined with a radius of 25 pixels (or equivalently 25"). The spacecraft rolls with a period of 1.645 h. While it rolls, background objects enter and leave the extended photometric aperture, leading to periodic contamination effects. The closest neighbours which PSF may enter the default aperture around HD172555 have a magnitude difference of 4.2 mag with HD 172555 (see Fig. 1). The DRP provides an estimate of the contamination effects due to background objects. The average contamination level that is found for HD 172555's light curve is 155 ppm, which is on the same order of magnitude as the combined differential photometric precision (CDPP; Hoyer et al. 2020) of the CHEOPS flux at 10 min cadence, or ~180 ppm. However, the full amplitude of the correlation between the roll angle and contamination is 50 ppm (Fig. 2), thus it is 3.6 times smaller than the CDPP.

## 2.2. The HD 172555 light curve

The CHEOPS light curve of HD 172555 is shown in Fig. 3. The observations are not continuous and have a large number of interruptions or gaps, about 45% of the total monitoring duration. Those are mainly due to the passage of the star behind the Earth, and crossings of the spacecraft over the South Atlantic Anomaly.

As summarised in Table 2, the root mean square (RMS) of the relative flux in the light curve is 1581 ppm, but the CDPP that ignores large-scale variations is much lower, 180 ppm at a 10 min cadence and 100 ppm at 3 h cadence. This provides a better estimate of the short-cadence flux dispersion. The light curve is indeed strongly dominated by large variations that are of stellar origin.

The time scales of the variability can be clearly seen in the Lomb-Scargle periodogram (LSP), where the variations dominate in the range ~200–800 µHz (Fig. 4). We attributed these strong periodic variations to stellar oscillations, which, given the position of the star ($T_{\text{eff}}$ = 7800 K, $\log g$ = 4.3; Table 1) in a $T_{\text{eff}}$ – $\log g$ diagram such as Fig. 10 of Uytterhoeven et al. (2011) and the oscillation frequencies that have been identified, are of a $\delta$ Scuti nature.

We aim to search for weaker signals, such as transits, that are possibly present in the light curve. At this stage of the reduction, those would still be hidden behind the strong oscillations. Our purpose now is thus to remove this variability.

# 3. HD 172555's $\delta$ Scuti oscillations

First, we stress that from the frequency extraction carried out below, there is no intention to use the stellar oscillations to attempt to model the star with asteroseismology. The observation duration of about *T* = 2.1 days is too short for this purpose, as it results in a resolution of only 1/*T* = 5.6 µHz. Instead, our aim is to clean the light curve as much as possible, from any periodic variability of stellar or instrumental origin, in order to search for transient signatures.

To remove the periodic oscillations, we applied the standard pre-whitening technique (Deeming 1975) using the tool Felix (Frequency Extraction for Light curve Exploration; Charpinet et al. 2010; Zong et al. 2016). In short, in the LSP, we identified the frequency and amplitude of the highest-amplitude peak, which are used as initial guesses in a subsequent non-linear least square (NLLS) fit of a cosine wave in a time domain using the Levenberg-Marquardt algorithm. The fitted wave of derived frequency, amplitude, and phase was then subtracted from the light curve. The operation was repeated as long as there was a peak above a pre-defined threshold, defined as a given level of the signal-to-noise ratio (S/N). The *S/N* = 1 level – the noise – is defined locally as the median of the points within a gliding window (centred on each point of the LSP) of ~300 times the resolution of the data. This median was re-evaluated at each step of the pre-whitening, that is each time a peak was removed.

We detrended the light curve from the oscillations at different levels of significance (as determined by the S/N), from relatively aggressive to conservative. The minimum significance required for identifying a peak to be removed is 4$\sigma$, that is a false alarm probability of $3.2 \times 10^{-5}$. This 4$\sigma$ significance level can be converted into an $S/N = x$ level. To determine the value of $x$, we used the approach developed in Zong et al. (2016): using the same time sampling (i.e. a cadence of 43.2s) and the same window (i.e. accounting for gaps) as in the real light curve, we simulated 10 000 pure Gaussian white-noise light curves. For a given S/N threshold, we then searched for the number of times that at least one peak in the LSP of these artificial light curves (that are, by construction, just noise) happen to be above this threshold. We obtained the false alarm probability by dividing by the number of tests (10 000 here). We found that the threshold corresponding to a 4$\sigma$ significance (false alarm probability of $3.2 \times 10^{-5}$) is $S/N = 4.8$. Table A.1 presents the properties
<figure><img src="image_4.png" /><figcaption>Fig. 3. CHEOPS light curve of HD 172555, with $t_0 = 0$ at MJD-59 381.445.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Fig. 4. LSP of the CHEOPS light curve of HD 172555. The red line corresponds to a 4σ significance. The Nyquist frequency is at 11 579 µHz. The insert shows a zoomed window on the low frequency domain <1,100 µHz.</figcaption></figure>

of the periodic variations extracted in the CHEOPS light curve down to $S/N = 4.8$. Frequencies are detected from 103.4 µHz almost up to the Nyquist frequency of 11 579 µHz. From $f_{34}$ (1011.23 µHz) and beyond, it should be noticed that all peaks correspond to multiples of the CHEOPS orbital frequency (orbital period of 98.9 min). Thus, contaminations correlated to the roll angle of the spacecraft with amplitudes <200 ppm were still present in the light curve, which FELIX succeeded to identify and remove. More conservative thresholds that we applied for removing periodic variabilities are $S/N > 6, 8, \& 10$, and finally, most conservatively, removing only the ten dominant frequencies. The Fig. 5 shows the effect of the different detrending schemes.

# 4. Transient signature search

## 4.1. Hint of a transit-like event in the binned residual light curves

Our purpose is now to analyse the detrended light curves of HD 172555 in the search for any transient events of circumstellar origin. We expect to witness them in such a young inclined system. Some short-scale (<1 h) variabilities are still present in the detrended light curves, which cannot be of circumstellar origin because a central transit of the star at a distance >1 $R_\star$ is at

least 1.2 h long. We thus averaged out this variability by applying a mean filter to the data, that is equivalent to binning. A window size of 1.645 h, equal to the CHEOPS orbital period, offers a good compromise between decreasing noise, removing excess variability, synchronising the bins and the gaps, and enhancing the prominence of >2 h variations.

Before calculating the binned light curve, we linearly interpolated the data on a grid of epochs with a constant step of 43 s. Epochs within gaps were treated as missing data. For each bin, the binned flux is the mean among all epochs within the binning window. The equation that we used to calculate error bars is:

$$ \sigma_i^{\text{bin}} = \sqrt{\sum_{j=1}^{N} \frac{\sigma_{i,j}^2}{N^2} + \frac{s_i^2}{N}}, \qquad (1) $$

where i is the index of a binned point, j is a sub-index of the binning window going from 1 to N, $\sigma_{i,j}$ is the intrinsic (CHEOPS) flux error bar of the considered epoch in the light curve, and $s_i$ is the standard deviation of flux within bin i. The leftmost term is the quadratic sum of the flux errors of all epochs within bin i divided by $N^2$; this is the variance on the mean of a set of random variables. Missing data (~45% of the time span) are taken into account by attributing an uncertainty to them that is equal to one standard deviation of the full light curve $\sigma_{lc}$. The rightmost
<figure><img src="image_6.png" /></figure>

**Fig. 5.** HD 172555's light curve after the removal of δ Scuti variations with different degrees of precision (see text). From top to bottom: raw light curve (black), the ten dominant frequencies subtracted (blue), and the oscillations subtracted with $S/N > 10$ (orange), $S/N > 8$ (green), $S/N > 6$ (red), and $S/N > 4.8$ (purple). The vertical scale is common to all curves, allowing for the effect of the different δ Scuti removal schemes on the residuals' amplitude to be compared.

term is the square of the standard deviation of flux within bin *i*, divided by *N*; this is the estimator of the variance of the mean within a sample of values (here bin *i*).

Figure 6 shows the effect of binning on the different detrended light curves obtained after subtracting δ Scuti oscillations and other periodic variations at various thresholds (see also Sect. 3). In the binned detrended light curves, there is one prominent feature, an absorption at about *t* = 0.3 days with an apparent depth 300–500 ppm, common to every detrending scheme. Another series of dips attract attention near the edge within 1.6–1.8 days. However, they are less prominent than the *t* = 0.3-days feature, and with a short duration, they are similar to the bin period (1.645 h), and thus more likely to originate from noise. For this reason, in this work, we only focus on the *t* = 0.3-days feature, our transient candidate.

In all the detrended light curves, residual δ Scuti oscillations dominate over other stochastic noise. Even though it is apparently not periodic in the 2-days long data, we cannot fully exclude that such a transient signal is due to the constructive interference of several δ Scuti oscillations. We have checked, however, that the frequencies involved in a transient with such a duration (~0.4 days) rather stand below 100 μHz (or *T* = 0.1 days). As shown in the light curve LSP (Fig. 4), a signal at these frequencies with an amplitude of >300 ppm is well separated from the majority of frequencies (500–800 μHz) attributed to the δ Scuti pulsator. This suggests that the identified transient is not an artefact of the stellar oscillations.

<figure><img src="image_7.png" /></figure>

**Fig. 6.** Binning of the raw and residual light curves with a 1.645 h bin. The colour code is the same as in Fig. 5. The dots show, in comparison, the original, non-binned, light curve.

Furthermore, in Sect. 4.2, we are able to exclude that this transient is due to instrumental or reduction systematics. In Sect. 4.3, assuming this transient might be produced by the passage of an exocomet in front of HD 172555, we fitted it with an exocomet transit model with a good match. We then show in Sect. 4.4 that an injected exocomet transit of a few hundred parts-per-million is expected to survive the δ Scuti oscillations' removal process applied above, down to $S/N \sim 6$.

## 4.2. Systematics

We first compared the detrended light curve with the systematics of CHEOPS: thermal ramp, PSF principal component coefficients, and centroid variations (Fig. 7). The variations were binned, as the light curves above, with a 1 h timestamp. We used the reduction pipeline called PIPE to extract the PSF photometry (Brandeker et al., in prep.; see also descriptions in Szabó et al. 2021; Brandeker et al. 2022), including the PSF principal component coefficients and the centroid variations.

We identified a long-term variation with an exponential-decay shape. This variation is indeed present in the data as well and commonly identified as being due to thermal relaxation of the telescope tube, as recorded by the thermal front sensor of CHEOPS (Deline et al. 2022). However, we did not observe short timescale transient-like signals in these diagnostics. This suggests that the transient signal is of astrophysical origin, and not of instrumental origin.
<figure><img src="image_8.png" /><figcaption>Fig. 7. Systematics variations, from top to bottom: the recording from one of the thermal front sensors (in arbitrary units); the PSF principal components' coefficients $u_1, u_2, u_3$, and $u_4$ (in arbitrary units); and the PSF centroid $(x_c, y_c)$ (in pixels).</figcaption></figure>

Moreover, the PIPE PSF photometry gives an independent photometric extraction with potentially different systematics from the aperture photometry of the DRP. In the present case, PIPE gives very similar results to the DRP with only marginal improvements in the mean absolute deviation (163 ppm for PIPE compared to 168 ppm for the DRP). The transient signal is still present. In the rest of the analysis, we use only the DRP reduction.

### 4.3. Transit signature fit

HD 172555 is known to be surrounded by a debris disk and transiting exocomets (K14; Grady et al. 2018). Given this information, we make the hypothesis that the identified transient absorption feature is due to the transit of an exomet in front of HD 172555. Such exomet's photometric transit has already been witnessed in $\beta$ Pictoris (Zieba et al. 2019; Pavlenko et al. 2022; Lecavelier Des Etangs et al. 2022).

We could fit this feature with a photometric exomet transit model of the relative flux variation $\Delta F/F$. We followed Lecavelier Des Etangs et al. (2022) and used a 1D empirical model of an exomet transit based on thorough numerical simulations of transit shapes of an evaporating comet by Lecavelier Des Etangs et al. (1999):

$$ \frac{\Delta F(t)}{F(t)} = K \left( e^{-\beta(t-t_0)} H(t-t_0) - e^{-\beta(t-t_0-\Delta t)} H(t-t_0-\Delta t) \right), \quad (2) $$

with $H(x)$ being the Heaviside function. In this model with four parameters, $t_0$ is the time of the beginning of transit of the head of the exomet; $\beta$ is the speed of the transit of one scale length of the cometary tail; $\Delta t$ is the duration of the transit of the comet nucleus; and $K$ is a scalar amplitude related to absorption depth AD – that is, the deepest relative flux variation during the transit – through AD = $K(1 - \exp(-\beta\Delta t))$. We used $\beta$ from 0 to 100 days$^{-1}$ and $\Delta t$ from 0.04 to 0.4 days, corresponding to periastron distances ranging from about 0.01 to 0.13 au.

More importantly, from these quantities, we could estimate physical properties of the exomet nucleus using a scaling relation based on the well-constrained properties of the bright solar system's Hale-Bopp comet (Jewitt & Matthews 1999). It has a radius of ~30 km with a dust production rate of $2 \times 10^6$ kg s$^{-1}$ at 1 au that can be translated into the following scaling relation (Lecavelier Des Etangs et al. 2022):

$$ \dot{M}_{\text{dust, 1 au}} = 2 \times 10^6 \text{ kg s}^{-1} \left( \frac{r_{\text{nucl}}}{30 \text{ km}} \right)^2 \left( \frac{L_{\star}}{L_{\odot}} \right), \quad (3) $$

with $\dot{M}_{\text{dust, 1 au}}$ being the dust-mass production rate of the exomet if it were at 1 au from the central star, $r_{\text{nucl}}$ being the radius of the solid comet nucleus, and $q$ being the periastron distance to the central star. The Hale-Bopp dust production rates, in agreement with Lecavelier Des Etangs et al. (1999)'s model assumptions, follow a $\dot{M} \propto 1/q^2$ law, with a decrease by a factor ~10 when distance increases from 1 to 3 au (Jewitt & Matthews 1999).

We also used equations relating the model parameters to physical properties of the exocomets. They are based on the following identities (Lecavelier Des Etangs et al. 2022):

$$ AD = 5 \times 10^{-5} \left( \frac{\dot{M}_{1 \text{ au}}}{10^5 \text{ kg s}^{-1}} \right) \left( \frac{q}{1 \text{ au}} \right)^{-1/2} \left( \frac{M_{\star}}{M_{\odot}} \right) \quad (4) $$

$$ \Delta t_{\text{transit}} = \frac{\bar{L}_{\text{path}}}{v_{\text{transit}}} \quad (5) $$

$$ \text{with } v_{\text{transit}} \approx \sqrt{\frac{2GM_{\star}}{q}} \quad \text{and} \quad \bar{L}_{\text{path}} = \frac{\pi R_{\star}}{2}. $$

The duration of transit $\Delta t_{\text{transit}}$ is given by the average path length through the stellar disk divided by the transit velocity $v_{\text{transit}}$. Using the stellar parameters of HD 172555 given in Table 1, these identities translate into the following:

$$ q \approx 0.013 \left( \frac{\Delta t}{1 \text{ h}} \right)^2 \text{ au} \quad (6) $$

$$ \dot{M}_{1 \text{ au}} \approx 3.7 \times 10^5 \left( \frac{AD}{1000 \text{ ppm}} \right) \sqrt{\frac{q}{0.1 \text{ au}}} \text{ kg s}^{-1} \quad (7) $$

$$ r_{\text{nucl}} \approx 2.4 \left( \frac{\dot{M}_{1 \text{ au}}}{10^5 \text{ kg s}^{-1}} \right)^{1/2} \text{ km.} \quad (8) $$

Previous exomet detections in Ca II spectra of the HD 172555 system have a radial velocity within ±20 km s$^{-1}$ from the star's systemic velocity, and a distance to the star ≤0.3 au such that Calcium can sublimate at a rate enabling the ion cloud to cover a significant portion of the line of sight (Beust et al. 1990). Their transit should thus happen shortly before or after the periastron passage along their near-parabolic orbit. In the above formula, we thus assumed that the exomet is transiting the stellar disk at a locus close to the periastron, with a distance to the central star ~q.
<table><caption>**Table 3.** Fit results with the exocomet transit model of Eqs. (2)–(8).</caption><thead><tr><th>Parameters</th><th>Unit</th><th>10 dom. freq.</th><th>S/N &gt; 10</th><th>S/N &gt; 8</th><th>S/N &gt; 6</th><th>S/N &gt; 4.8</th></tr></thead><tbody><tr><td>t<sub>0</sub></td><td>BJD-t<sub>min</sub></td><td>0.273 ± 0.015</td><td>0.271 ± 0.012</td><td>0.271 ± 0.011</td><td>0.270 ± 0.010</td><td>0.310 ± 0.020</td></tr><tr><td>β</td><td>day<sup>-1</sup></td><td>43 ± 38</td><td>8.6 ± 2.3</td><td>14.3 ± 4.5</td><td>14.3 ± 4.3</td><td>44.0 ± 29.4</td></tr><tr><td>Δt</td><td>h</td><td>3.75 ± 0.40</td><td>1.69 ± 0.59</td><td>1.90 ± 0.54</td><td>1.93 ± 0.50</td><td>1.00 ± 0.57</td></tr><tr><td>K</td><td>ppm</td><td>337 ± 73</td><td>920 ± 380</td><td>497 ± 187</td><td>435 ± 148</td><td>224 ± 11</td></tr><tr><td>χ<sup>2</sup></td><td>DoF = 2212</td><td>284 776</td><td>11 9804</td><td>78 834</td><td>56 445</td><td>31 623</td></tr><tr><td>χ<sup>2</sup><sub>null</sub></td><td>DoF = 2216</td><td>289 145</td><td>123 559</td><td>81 087</td><td>58 240</td><td>32 009</td></tr><tr><td>F-test</td><td></td><td>8.5</td><td>17.3</td><td>15.8</td><td>17.6</td><td>6.8</td></tr><tr><td>p-value</td><td></td><td>8.6 × 10<sup>-7</sup></td><td>5.2 × 10<sup>-14</sup></td><td>9.2 × 10<sup>-13</sup></td><td>3.2 × 10<sup>-14</sup></td><td>2.1 × 10<sup>-5</sup></td></tr><tr><td>AD</td><td>ppm</td><td>336 ± 73</td><td>417 ± 219</td><td>336 ± 148</td><td>297 ± 119</td><td>188 ± 110</td></tr><tr><td>q</td><td>au</td><td>0.186 ± 0.014</td><td>0.038 ± 0.010</td><td>0.047 ± 0.010</td><td>0.049 ± 0.009</td><td>0.013 ± 0.005</td></tr><tr><td>M&#770;<sub>dust, 1 au</sub></td><td>10<sup>5</sup> kg s<sup>-1</sup></td><td>1.81 ± 0.40</td><td>1.01 ± 0.55</td><td>0.92 ± 0.41</td><td>0.82 ± 0.34</td><td>0.27 ± 0.17</td></tr><tr><td>r<sub>nucl</sub></td><td>km</td><td>3.38 ± 0.37</td><td>2.53 ± 0.68</td><td>2.40 ± 0.54</td><td>2.28 ± 0.47</td><td>1.30 ± 0.41</td></tr></tbody></table>

We fitted all differently detrended light curves with this exocomet transit model and a 2-degree polynomial to account for the long-term trend identified as a thermal relaxation. The results of the fit are summarised in Table 3 and the model compared to the binned light curves is shown in Fig. 8. We derived the $\chi^2$ and an F test, comparing the cases of fitting an exocomet transit profile plus a 2-degree polynomial, and as the null hypothesis, a 2-degree polynomial only.

The depth of the signature indeed decreases for detrending deeper than $S/N > 8$. The shape of the transit profile also seems strongly affected by residual $\delta$ Scuti oscillations if only the ten dominant frequencies are subtracted. It can be seen in the null hypothesis $\chi^2_{null}$ that oscillations still dominate the error budget of the detrended light curve, even if they are removed down to $S/N = 4.8$. Nevertheless, the adjunction of the transit signature to the fitted model always significantly improves the $\chi^2$, with an F test >6.8 ($p$-value < $10^{-5}$).

It is unclear, however, at what S/N threshold of the $\delta$ Scuti detrending schemes would the derived parameters be the closest to the true transiting exocomet, if any. To understand how $\delta$ Scuti removal may also remove transit signal, we have to test transit injection recovery.

### 4.4. Transit injection recovery

Here, we check the following: that an exocomet transit signal could indeed survive the $\delta$ Scuti removal process; at what threshold the transit parameters are best retrieved; and, especially, below what threshold such a transit signal is identified by FELIX and removed as if it were of $\delta$ Scuti origin. To do so, we injected a transient into the raw light curve of HD 172555 prior to the detrending scheme presented in Sect. 3. We injected a signal at $t - t_0 = 1.2$ days similar to the one potentially detected here (Sect. 4.3), with a transit depth of 544 ppm and a full width at half maximum of $\Delta t = 0.25$ days. All parameters of the injected model are given in Table 4.

Figure 9 shows the detrended light curves when removing the periodic variabilities down to $S/N > 4.8$, 6, 8, and 10. The injected signal is clearly detected in all the detrending schemes. For the '$S/N > 4.8$' scheme, the injected signal is strongly affected by the removal of oscillations and the absorption depth is strongly damped to below 544 ppm (as injected), while the $t = 0.3$-day candidate transient has disappeared.

In Table 4 we added the parameters fitted at the position of the injected transient, showing that the transit time, $\Delta t$, and K

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Transient signature detected at $t - t_0 = 0.3$ day and modelled by a transiting exocomet with the model given in Eqs. (2)–(8), compared to the 1-h binned detrended light curves. From top to bottom, the colour codes identify the same datasets as in Fig. 6.</figcaption></figure>

parameters derived for the detrended light curves with $S/N > 8–10$ are almost not affected by removal of periodic variabilities. The value derived for $\beta$ is strongly discrepant with the injected value, showing the difficulty to have a good estimate as to the real signal. The most accurate estimation of periastron distance $q$, dust production rate $M_{dust,1 au}$, and nucleus radius $r_{nucl}$ were obtained by the fit of the injected transit signature in the '$S/N > 8$' detrending scheme.

### 4.5. Final adopted parameters of the candidate exomet

In this scheme, the candidate exocomet transit model detected in this work, fitted in Sect. 4 with the results shown in Table 3,
<table><caption>**Table 4.** Injection model and fits (see text for explanations).</caption><thead><tr><th>Parameters</th><th>Unit</th><th>Injected</th><th>S/N &gt; 10</th><th>S/N &gt; 8</th><th>S/N &gt; 6</th><th>S/N &gt; 4.8</th></tr></thead><tbody><tr><td>t<sub>0</sub></td><td>BJD-t<sub>min</sub></td><td>1.2</td><td>1.226 ± 0.005</td><td>1.232 ± 0.003</td><td>1.228 ± 0.004</td><td>1.250 ± 0.003</td></tr><tr><td>β</td><td>day<sup>-1</sup></td><td>15</td><td>43.0 ± 13.0</td><td>98.1 ± 31.8</td><td>58.3 ± 18.0</td><td>45.9 ± 12.4</td></tr><tr><td>Δt</td><td>h</td><td>2.4</td><td>2.11 ± 0.15</td><td>2.18 ± 0.08</td><td>2.02 ± 0.12</td><td>1.00 ± 0.18</td></tr><tr><td>K</td><td>ppm</td><td>700</td><td>697 ± 85</td><td>582 ± 44</td><td>460 ± 51</td><td>601 ± 146</td></tr><tr><td>χ<sup>2</sup></td><td>DoF = 2211</td><td></td><td>121 972</td><td>90 049</td><td>56 594</td><td>32 371</td></tr><tr><td>χ<sup>2</sup><sub>null</sub></td><td>DoF = 2215</td><td></td><td>133 317</td><td>100 737</td><td>61 623</td><td>34 136</td></tr><tr><td>F-test</td><td></td><td></td><td>51.4</td><td>65.6</td><td>49.1</td><td>30.1</td></tr><tr><td>AD</td><td>ppm</td><td>544</td><td>681 ± 85</td><td>582 ± 44</td><td>456 ± 51</td><td>513 ± 136</td></tr><tr><td>q</td><td>au</td><td>0.037</td><td>0.029 ± 0.002</td><td>0.031 ± 0.001</td><td>0.026 ± 0.002</td><td>0.007 ± 0.001</td></tr><tr><td>M&#775;<sub>dust, 1 au</sub></td><td>10<sup>5</sup> kg s<sup>-1</sup></td><td>1.24</td><td>1.36 ± 0.18</td><td>1.20 ± 0.09</td><td>0.87 ± 0.10</td><td>0.49 ± 0.14</td></tr><tr><td>r<sub>nucl</sub></td><td>km</td><td>2.82</td><td>2.96 ± 0.18</td><td>2.78 ± 0.11</td><td>2.37 ± 0.13</td><td>1.77 ± 0.23</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Recovery of an injected transient signal (green filled area) at about $t - t_0 = 1.2$ days from detrended light curves where the removed $\delta$ Scuti oscillations have the same S/N levels as analysed here with the same colour code as in Fig. 5. The raw light curve with the injected signal is shown on top (black line). The adopted bin step is also 1.645 h.</figcaption></figure>

has $\beta = 14 \pm 5$ day$^{-1}$, $\Delta t = 1.9 \pm 0.5$ h, and $K = 500 \pm 200$ ppm. This leads to an absorption depth AD ~ 340±150 ppm, a periastron distance $q \sim 0.05 \pm 0.01$ au, a dust production rate $\dot{M} \sim (0.9 \pm 0.4) \times 10^5$ kg s$^{-1}$, and a nucleus size of 2.5±0.5 km. This is our best estimate of the true profile, assuming the transient is indeed due to an exocomet transit.

**5. Comparing this candidate exocomet to K14 spectroscopic detections**

It is interesting to make a comparison between the present candidate detection and those reported in K14. It can allow us to understand if they belong to similar or different populations.

Some parameters of their detections are reported in Table 5. Using the formulas from Kiefer et al. (2014b), we could calculate the evaporation efficiency of all those detections, with the evaporation efficiency defined as the log ratio between the energy spent to evaporate gas and dust from the exocomet nucleus and the incoming energy flux. For HD 1725555, it is related to the dust production rate, with $d$ being the distance to the star in $R_\star$ as

$$ \eta = \log(\dot{M}_{\text{dust}} d^2) - 1.86. \qquad (9) $$

For variable absorption lines in Ca II K & H spectra due to the occultation of the stellar surface by a cloud of atomic gas evaporated from an exocomet, the depth of the variable lines in Ca II K & H are $\alpha(1 - e^{-A})$ and $\alpha(1 - e^{-A/2})$, respectively, with the oscillation strength ratio $f_H/f_K = 0.5$ (K14; Kiefer et al. 2014b). We introduced the cloud-to-star surface ratio $\alpha$ and the optical absorption depth $A$. The evaporation efficiency can be determined from $\alpha$ and $A$ (see Kiefer et al. 2014b for details) as

$$ \eta = 9.2 + \log[\alpha(1 - e^{-A} - e^{-A/2})]. \qquad (10) $$

Another useful formula from Kiefer et al. (2014b) allows for the distance to the star to be expressed with respect to the surface ratio $\alpha$ and the dust production rate:

$$ \alpha = 5 \times 10^{-14} \dot{M}^{4/3} d^{4/3}. \qquad (11) $$

Combining this equation with Eq. (9) leads to an expression of the distance, in stellar radius, with respect to the evaporation efficiency and the surface ratio

$$ d = 7.7 \times 10^{-9} 10^{\eta} \alpha^{-3/4}. \qquad (12) $$

It is then straightforward to derive the dust production rate at 1 au, $\dot{M}_{\text{dust, 1 au}}$, again assuming that the exocomet locus is close to its orbit periastron, such that $d \sim q$ is the periastron distance. Further, applying Eqs. (6)–(8) allows one to determine, for each K14 detection, an estimation of transit duration $\Delta t$ and the photometric transit absorption depth. All derived parameters are reported in Table 5. The error bars were obtained by: (i) sampling split-normal distributions of $\alpha$ and $A$ by drawing 10 000 samples; (ii) applying all formulas given above and in Sect. 4.3 to determine $\eta$, $q$, $\dot{M}_{\text{dust, 1 au}}$, AD, $\Delta t$, and $r_{\text{nucl}}$; and (iii) calculating the 16th, 50th, and 84th percentiles of the samples to determine the main values and the $1\sigma$ uncertainties given in Table 5.
<table><caption>**Table 5.** Summary of some parameters from K14's detections.</caption><thead><tr><th>Date</th><th>Δt<sub>min</sub></th><th>α</th><th>A</th><th>η</th><th>q<br>(au)</th><th>M&#x0307;<sub>dust, 1 au</sub><br>(10<sup>5</sup> kg s<sup>-1</sup>)</th><th>AD<br>(ppm)</th><th>Δt<br>(hour)</th><th>r<sub>nucl</sub><br>(km)</th></tr></thead><tbody><tr><td>22/09/2004</td><td>3.5 h</td><td>≥0.84</td><td>0.07 ± 0.005</td><td>8.17 ± 0.03</td><td>0.0088 ± 0.0006</td><td>5.6 ± 0.5</td><td>5050 ± 250</td><td>0.81 ± 0.03</td><td>5.66 ± 0.23</td></tr><tr><td>21/08/2005</td><td>3 min</td><td>0.04<sup>+0.04</sup><sub>-0.01</sub></td><td>1.7 ± 1.1</td><td>8.05 ± 0.30</td><td>0.059<sup>+0.014</sup><sub>-0.023</sub></td><td>4.3 ± 2.5</td><td>1530 ± 680</td><td>2.12<sup>+0.24</sup><sub>-0.46</sub></td><td>5.0 ± 1.5</td></tr><tr><td>08/07/2010</td><td>3 min</td><td>≥0.024</td><td>≤10</td><td>9.16<sup>+0.22</sup><sub>-0.57</sub></td><td>0.132<sup>+0.026</sup><sub>-0.048</sub></td><td>54 ± 38</td><td>12 800 ± 7500</td><td>3.16<sup>+0.30</sup><sub>-0.64</sub></td><td>17.6<sup>+5.1</sup><sub>-8.5</sub></td></tr><tr><td>11/06/2011</td><td>1.9 h</td><td>0.04 ± 0.01</td><td>1.7 ± 0.5</td><td>7.93 ± 0.14</td><td>0.054<sup>+0.009</sup><sub>-0.010</sub></td><td>3.23 ± 0.95</td><td>1210 ± 280</td><td>2.02 ± 0.18</td><td>4.30 ± 0.65</td></tr></tbody></table>

**Notes.** For comparison, the evaporation efficiency obtained for the transient reported here, if due to an exocomet, is $η=7.39±0.14$ for $q = 0.05 ± 0.01$ au and $M&#x0307;<sub>dust, 1 au</sub> = 0.92 ± 0.41 × 10<sup>5</sup> kg s<sup>-1</sup>.

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Periastron distance and evaporation efficiency for the exo- comets detected in K14 (blue circles) and the candidate exocomet reported in this work (orange star).</figcaption></figure>

The evaporation efficiency of the CHEOPS exocomet transit reported here can be calculated similarly to Eq. (9) by using the measured $M&#x0307;<sub>dust, 1 au</sub>$ and fixing $d = 1$ au. We find $η = 7.39±0.14$ and an estimated periastron distance of $0.05±0.01$ au, or equivalently $6.8±1.4 R_★$. It is graphically compared to the K14 detections in Fig. 10. The CHEOPS transit event and the 2004, 2005, and 2011 K14 events have comparable, though slightly larger, evaporation efficiencies with $η_{K14} ≈ 8.05±0.10$, but similar peri- astron distance $q_{K14} ≈ 0.04±0.02$ au. The small difference in evaporation efficiency translates to about four times larger evap- oration rates $M&#x0307;<sub>dust, 1 au</sub> ≈ (4.4±1.0)×10<sup>5</sup> kg s<sup>-1</sup> and twice as large nucleus radii ≈5.0±0.6 km. In comparison, the 2010 K14 event has radically different properties with $η > 8.5$, $q > 0.08$ au, $M&#x0307;<sub>dust, 1 au</sub> > 16 × 10<sup>5</sup> kg s<sup>-1</sup>$, and $r_{nucl} > 9$ km. The 2004, 2005, and 2011 K14 spectroscopic detections and the candidate pre- sented in this work hence most likely belong to a common class of exocomets in this system to which the 2010 K14 event would conversely not belong.

The calculated AD and $Δt$ for all K14 detections show that if they were observed in photometry, they would have an even larger amplitude in the light curve. This implies that the HD 172555 exocomets that are observed transiting in spec- troscopy would produce a detectable photometric signature as well. Given the higher evaporation rates of K14 detections if they are due to a selection effect, the reverse, however, might not always be true. It, nevertheless, opens the opportunity to observe the HD 172555 exocomets simultaneously with a photometer and

a spectrograph in order to detect and measure both the gas and dust evaporated from a single exocomet at once. This could put strong constraints on the relative abundances of chemical ele- ments in exocomets nucleus, that is to say volatiles in ices and refractories in dust grains.

## 6. Conclusion

With CHEOPS photometric monitoring of HD 172555, we have shown that this young star is a $δ$ Scuti pulsator. Once the $δ$ Scuti oscillations were removed, we discovered a candidate transient event that nicely matched an exocomet transit signature. Fitting the transient light curve with an exocomet transit model, we mea- sured a dust production rate at 1 au of $M&#x0307;<sub>dust, 1 au</sub> ~ 10<sup>5</sup> kg s<sup>-1</sup>. Using the scaling relationship between the dust production rate and nucleus radius (Lecavelier Des Etangs et al. 2022), we deduced that the candidate HD 172555 exocomet hinted at in this work would have a nucleus size of 2.5±0.5 km. This radius is similar to those measured for exocomets in the $β$ Pic sys- tem with radii from 1.5 to 6.7 km (Lecavelier Des Etangs et al. 2022), and for JFC in the Solar System with typical radii of 2–5.5 km (Tancredi et al. 2006).

More observations on longer uninterrupted time periods of HD 172555, typically longer than 3 days, are needed to better characterise the $δ$ Scuti variations and more exquisitely remove them from the measured light curve. This is a crucial step to uncover transient absorptions and exclude constructive interfer- ences from residual, low-amplitude $δ$ Scuti oscillations. More- over, if the reported transient is indeed a transiting exocomet, we would anticipate with more than 89% probability a new detection if HD 172555 is observed on a time interval longer than 4 days, and a probability of detecting at least two transiting exocomets of at least 73%. After $β$ Pictoris, the young system of HD 172555 is therefore one of the best targets for investigating the statis- tics and composition of exocomets, using both spectroscopy and photometry.

**Acknowledgements.** CHEOPS is an European Space Agency (ESA) mission in partnership with Switzerland with important contributions to the payload and the ground segment from Austria, Belgium, France, Germany, Hungary, Italy, Portugal, Spain, Sweden, and the United Kingdom. The CHEOPS Consortium would like to gratefully acknowledge the support received by all the agencies, offices, universities, and industries involved. Their flexibility and willingness to explore new approaches were essential to the success of this mission. F.K., A.L. and G.H. acknowledge funding from the Centre National d'Etudes Spatiales and from the French National Research Agency (ANR) under contract number ANR-18-CE31-0019 (SPIaSH). F.K. also acknowledges support from the Université Paris Sciences et Lettres under the DIM-ACAV program Origines et conditions d'apparition de la vie. V.V.G. is an F.R.S-FNRS Research Associate. G.M.S. acknowledges the support of the Hungarian National Research, Development and Innovation Office (NKFIH) grant K-125015, a PRODEX Experiment Agreement No. 4000137122, the Lendület LP2018-7/2021
grant of the Hungarian Academy of Science and the support of the city of Szombathely. A.B. was supported by the SNSA. A.C.C. acknowledges support from STFC consolidated grant numbers ST/R000824/1 and ST/V000861/1, and UKSA grant number ST/R003203/1. A.C.C. and T.W. acknowledge support from STFC consolidated grant numbers ST/R000824/1 and ST/V000861/1, and UKSA grant number ST/R003203/1. S.G.S. acknowledges support from FCT through FCT contract nr. CEECIND/00826/2018 and POPH/FSE (EC). D.G. gratefully acknowledges financial support from the CRT foundation under Grant No. 2018.2323 "Gaseous rocky? Unveiling the nature of small worlds". Y.A. and M.J.H. acknowledge the support of the Swiss National Fund under grant 200020_172746. We acknowledge support from the Spanish Ministry of Science and Innovation and the European Regional Development Fund through grants ESP2016-80435-C2-1-R, ESP2016-80435-C2-2-R, PGC2018-098153-B-C33, PGC2018-098153-B-C31, ESP2017-87676-C5-1-R, MDM-2017-0737 Unidad de Excelencia Maria de Maeztu-Centro de Astrobiología (INTA-CSIC), as well as the support of the Generalitat de Catalunya/CERCA programme. The MOC activities have been supported by the ESA contract No. 4000124370. S.C.C.B. acknowledges support from FCT through FCT contracts nr. IF/01312/2014/CP1215/CT0004. X.B., S.C., D.G., M.F. and J.L. acknowledge their role as ESA-appointed CHEOPS science team members. This project was supported by the CNES. The Belgian participation to CHEOPS has been supported by the Belgian Federal Science Policy Office (BELSPO) in the framework of the PRODEX Program, and by the University of Liège through an ARC grant for Concerted Research Actions financed by the Wallonia-Brussels Federation. L.D. is an F.R.S.-FNRS Postdoctoral Researcher. This work was supported by FCT - Fundação para a Ciência e a Tecnologia through national funds and by FEDER through COMPETE2020 - Programa Operacional Competitividade e Internacionalização by these grants: UID/FIS/04434/2019, UIDB/04434/2020, UIDP/04434/2020, PTDC/FIS-AST/32113/2017 & POCI-01-0145-FEDER-032113, PTDC/FIS-AST/28953/2017 & POCI-01-0145-FEDER-028953, PTDC/FIS-AST/28987/2017 & POCI-01-0145-FEDER-028987, O.D.S.D. is supported in the form of work contract (DL 57/2016/CP1364/CT0004) funded by national funds through FCT. B.-O.D. acknowledges support from the Swiss National Science Foundation (PP00P2-190080). This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (project FOUR ACES. grant agreement No 724427). It has also been carried out in the frame of the National Centre for Competence in Research PlanetS supported by the Swiss National Science Foundation (SNSF). D.E. acknowledges financial support from the Swiss National Science Foundation for project 200021_200726. M.F. and C.M.P. gratefully acknowledge the support of the Swedish National Space Agency (DNR 65/19, 174/18). M.G. is an F.R.S.-FNRS Senior Research Associate. S.H. gratefully acknowledges CNES funding through the grant 837319. K.G.I. is the ESA CHEOPS Project Scientist and is responsible for the ESA CHEOPS Guest Observers Programme. She does not participate in, or contribute to, the definition of the Guaranteed Time Programme of the CHEOPS mission through which observations described in this paper have been taken, nor to any aspect of target selection for the programme. This work was granted access to the HPC resources of MesoPSL financed by the Region Ile de France and the project Equip@Meso (reference ANR-10-EQPX-29-01) of the programme Investissements d'Avenir supervised by the Agence Nationale pour la Recherche. M.L. acknowledges support of the Swiss National Science Foundation under grant number PCEFP2_194576. P.M. acknowledges support from STFC research grant number ST/M001040/1. L.B., G.B., V.N., I.P., G.P., R.R., G.S., V.S., and T.Z. acknowledge support from CHEOPS ASI-NAF agreement n. 2019-29-HH.0. This work was also partially supported by a grant from the Simons Foundation (PI Queloz, grant number 327127). I.R.I. acknowledges support from the Spanish Ministry of Science and Innovation and the European Regional Development Fund through grant PGC2018-098153-B-C33, as well as the support of the Generalitat de Catalunya/CERCA programme. N.A.W. acknowledges UKSA grant ST/R004838/1.

## References

Alonso-Floriano, F. J., Caballero, J. A., Cortés-Contreras, M., Solano, E., & Montes, D. 2015, A&A, 583, A85

Barrado y Navascués, D., Stauffer, J. R., Song, I., & Caillault, J. P. 1999, ApJ, 520, L123

Benz, W. 2021, COSPAR Sci. Assembly, 43, 502

Beust, H., & Tagger, M. 1993, Icarus, 106, 42

Beust, H., Lagrange-Henri, A. M., Vidal-Madjar, A., & Ferlet, R. 1990, A&A, 236, 202

Binks, A. S., & Jeffries, R. D. 2014, MNRAS, 438, L11

Blackwell, D. E., & Shallis, M. J. 1977, MNRAS, 180, 177

Bonfanti, A., Delrez, L., Hooton, M. J., et al. 2021, A&A, 646, A157

Bottke, W. F., Vokrouhlický, D., Minton, D., et al. 2012, Nature, 485, 78

Brandeker, A., Heng, K., Lendl, M., et al. 2022, A&A, 659, L4

Cartwright, J. A., Hodges, K. V., & Wadhwa, M. 2022, Earth Planet. Sci. Lett., 590, 117576

Chambers, J. E., & Lissauer, J. J. 2002, Lunar Planet. Sci. Conf., 1093

Charpinet, S., Green, E. M., Baglin, A., et al. 2010, A&A, 516, L6

Deeming, T. J. 1975, Ap&SS, 36, 137

Deline, A., Hooton, M. J., Lendl, M., et al. 2022, A&A, 659, A74

de Sousa, R. R., Morbidelli, A., Raymond, S. N., et al. 2020, Icarus, 339, 113605

Engler, N., Schmid, H. M., Quanz, S. P., Avenhaus, H., & Bazzon, A. 2018, A&A, 618, A151

Erspamer, D., & North, P. 2003, A&A, 398, 1121

Feigelson, E. D., Lawson, W. A., Stark, M., Townsley, L., & Garmire, G. P. 2006, AJ, 131, 1730

Ferlet, R., Hobbs, L. M., & Vidal-Madjar, A. 1987, A&A, 185, 267

Gaia Collaboration 2020, VizieR Online Data Catalog: I/350

Gomes, R., Levison, H. F., Tsiganis, K., & Morbidelli, A. 2005, Nature, 435, 466

Grady, C. A., Brown, A., Welsh, B., et al. 2018, AJ, 155, 242

Gray, R. O., Corbally, C. J., Garrison, R. F., et al. 2006, AJ, 132, 161

Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27

Hoyer, S., Guterman, P., Demangeon, O., et al. 2020, A&A, 635, A24

Jewitt, D., & Matthews, H. 1999, AJ, 117, 1056

Kiefer, F., Lecavelier des Etangs, A., Augereau, J. C., et al. 2014a, A&A, 561, L10

Kiefer, F., Lecavelier des Etangs, A., Boissier, J., et al. 2014b, Nature, 514, 462

Kiefer, F., Lecavelier des Étangs, A., Vidal-Madjar, A., et al. 2017, A&A, 608, A132

Lagrange-Henri, A. M., Gosset, E., Beust, H., Ferlet, R., & Vidal-Madjar, A. 1992, A&A, 264, 637

Lecavelier Des Etangs, A., Vidal-Madjar, A., & Ferlet, R. 1999, A&A, 343, 916

Lecavelier des Etangs, A., Cros, L., Hébrard, G., et al. 2022, Sci. Rep., 12, 5855

Lisse, C. M., Chen, C. H., Wyatt, M. C., et al. 2009, ApJ, 701, 2019

Liu, B., Raymond, S. N., & Jacobson, S. A. 2022, Nature, 604, 643

Mamajek, E. E., & Bell, C. P. M. 2014, MNRAS, 445, 2169

Marks, N. E., Borg, L. E., Shearer, C. K., & Cassata, W. S. 2019, J. Geophys. Res. Planets, 124, 2465

Miret-Roig, N., Antoja, T., Romero-Gómez, M., & Figueras, F. 2018, A&A, 615, A51

Morbidelli, A., Petit, J. M., Gladman, B., & Chambers, J. 2001, Meteor. Planet. Sci., 36, 371

Morbidelli, A., Nesvorny, D., Laurenz, V., et al. 2018, Icarus, 305, 262

Nilsson, R., Liseau, R., Brandeker, A., et al. 2009, A&A, 508, 1057

Pavlenko, Y., Kulyk, I., Shubina, O., et al. 2022, A&A, 660, A49

Rappaport, S., Vanderburg, A., Jacobs, T., et al. 2018, MNRAS, 474, 1453

Riviere-Marichalar, P., Barrado, D., Augereau, J. C., et al. 2012, A&A, 546, L8

Ryder, G. 1990, EOS Transactions, 71, 313

Schanche, N., Hébrard, G., Collier Cameron, A., et al. 2020, MNRAS, 499, 428

Schneiderman, T., Matrà, L., Jackson, A. P., et al. 2021, Nature, 598, 425

Smith, R., Wyatt, M. C., & Haniff, C. A. 2012, MNRAS, 422, 2560

Szabó, G. M., Gandolfi, D., Brandeker, A., et al. 2021, A&A, 654, A159

Tancredi, G., Fernández, J. A., Rickman, H., & Licandro, J. 2006, Icarus, 182, 527

Tsiganis, K., Gomes, R., Morbidelli, A., & Levison, H. F. 2005, Nature, 435, 459

Uytterhoeven, K., Moya, A., Grigahcène, A., et al. 2011, A&A, 534, A125

Vidal-Madjar, A., Lagrange-Henri, A. M., Feldman, P. D., et al. 1994, A&A, 290, 245

Wilson, T. G., Goffo, E., Alibert, Y., et al. 2022, MNRAS, 511, 1043

Zieba, S., Zwintz, K., Kenworthy, M. A., & Kennedy, G. M. 2019, A&A, 625, L13

Zong, W., Charpinet, S., & Vauclair, G. 2016, A&A, 594, A46

Zuckerman, B., Song, I., Bessell, M. S., & Webb, R. A. 2001, ApJ, 562, L87

<footnote n="1">¹ Sorbonne Université, CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France</footnote>

<footnote n="2">² LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, 5 place Jules Janssen, 92195 Meudon, France
e-mail: flavien.kiefer@obspm.fr</footnote>

<footnote n="3">³ Space sciences, Technologies and Astrophysics Research (STAR) Institute, Université de Liège, Allée du 6 Août 19C, 4000 Liège, Belgium</footnote>

<footnote n="4">⁴ Institut d'astrophysique de Paris, UMR7095 CNRS, Université Pierre & Marie Curie, 98bis blvd. Arago, 75014 Paris, France</footnote>

<footnote n="5">⁵ ELTE Eötvös Lorandand University, Gothard Astrophysical Observatory, 9700 Szombathely, Szent Imre h. u. 112, Hungary</footnote>

<footnote n="6">⁶ MTA-ELTE Exoplanet Research Group, 9700 Szombathely, Szent Imre h. u. 112, Hungary</footnote>
⁷ Department of Astronomy, Stockholm University, AlbaNova University Center, 10691 Stockholm, Sweden

⁸ Physikalisches Institut, University of Bern, Sidlerstrasse 5, 3012 Bern, Switzerland

⁹ Center for Space and Habitability, Gesellsschaftstrasse 6, 3012 Bern, Switzerland

¹⁰ Centre for Exoplanet Science, SUPA School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews KY16 9SS, UK

¹¹ Observatoire Astronomique de l'Université de Genève, Chemin Pegasi 51, 1290 Versoix, Switzerland

¹² Instituto de Astrofísica e Ciencias do Espaco, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal

¹³ Dipartimento di Fisica, Universita degli Studi di Torino, via Pietro Giuria 1, I-10125, Torino, Italy

¹⁴ Observatoire de Haute-Provence, CNRS, Université d'Aix-Marseille, 04870 Saint-Michel-l'Observatoire, France

¹⁵ Instituto de Astrofísica de Canarias, 38200 La Laguna, Tenerife, Spain

¹⁶ Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain

¹⁷ Institut de Ciencies de l'Espai (ICE, CSIC), Campus UAB, Can Magrans s/n, 08193 Bellaterra, Spain

¹⁸ Institut d'Estudis Espacials de Catalunya (IEEC), 08034 Barcelona, Spain

¹⁹ Admatis, 5. Kandó Kálmán Street, 3534 Miskolc, Hungary

²⁰ Depto. de Astrofísica, Centro de Astrobiología (CSIC-INTA), ESAC campus, 28692 Villanueva de la Cañada (Madrid), Spain

²¹ Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, 4169-007 Porto, Portugal

²² Space Research Institute, Austrian Academy of Sciences, Schmiedlstrasse 6, 8042 Graz, Austria

²³ Center for Space and Habitability, University of Bern, Gesellschaftsstrasse 6, 3012 Bern, Switzerland

²⁴ Université Grenoble Alpes, CNRS, IPAG, 414 rue de la Piscine, 38400 St-Martin d'Hères, France

²⁵ Institute of Planetary Research, German Aerospace Center (DLR), Rutherfordstrasse 2, 12489 Berlin, Germany

²⁶ Université de Paris, Institut de physique du globe de Paris, CNRS, 75005 Paris, France

²⁷ Centre for Mathematical Sciences, Lund University, Box 118, 221 00 Lund, Sweden

²⁸ Aix-Marseille Univ, CNRS, CNES, LAM, 38 rue Frédéric Joliot-Curie, 13388 Marseille, France

²⁹ Astrobiology Research Unit, Université de Liège, Allée du 6 Août 19C, 4000 Liège, Belgium

³⁰ Observatoire Astronomique de l'Université de Genève, Chemin Pegasi 51, 1290 Versoix, Switzerland

³¹ Centre Vie dans l'Université des sciences, Université de Genève, Quai Ernest-Ansermet 30, 1211 Genève 4, Switzerland

³² Leiden Observatory, University of Leiden, PO Box 9513, 2300 RA Leiden, The Netherlands

³³ Department of Space, Earth and Environment, Chalmers University of Technology, Onsala Space Observatory, 439 92 Onsala, Sweden

³⁴ University of Vienna, Department of Astrophysics, Türkenschanzstrasse 17, 1180 Vienna, Austria

³⁵ Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK

³⁶ Science and Operations Department - Science Division (SCI-SC), Directorate of Science, European Space Agency (ESA), European Space Research and Technology Centre (ESTEC), Keplerlaan 1, 2201-AZ Noordwijk, The Netherlands

³⁷ Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, 1121 Budapest, Konkoly Thege Miklós út 15-17, Hungary

³⁸ ELTE Eötvös Loránd University, Institute of Physics, Pázmány Péter sétány 1/A, 1117, Budapest, Hungary

³⁹ IMCCE, UMR8028 CNRS, Observatoire de Paris, PSL Univ., Sorbonne Univ., 77 av. Denfert-Rochereau, 75014 Paris, France

⁴⁰ INAF, Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 35122 Padova, Italy

⁴¹ Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK

⁴² INAF, Osservatorio Astrofisico di Catania, Via S. Sofia 78, 95123 Catania, Italy

⁴³ Department of Astrophysics, University of Vienna, Tuerkenschanzstrasse 17, 1180 Vienna, Austria

⁴⁴ Institute of Optical Sensor Systems, German Aerospace Center (DLR), Rutherfordstrasse 2, 12489 Berlin, Germany

⁴⁵ Physikalisches Institut, University of Bern, Sidlerstrasse 5, 3012 Bern, Switzerland

⁴⁶ Dipartimento di Fisica e Astronomia "Galileo Galilei", Universita degli Studi di Padova, Vicolo dell'Osservatorio 3, 35122 Padova, Italy

⁴⁷ ETH Zurich, Department of Physics, Wolfgang-Pauli-Strasse 2, 8093 Zurich, Switzerland

⁴⁸ Cavendish Laboratory, JJ Thomson Avenue, Cambridge CB3 0HE, UK

⁴⁹ ESTEC, European Space Agency, Keplerlaan 1, 2201AZ, Noordwijk, The Netherlands

⁵⁰ Zentrum für Astronomie und Astrophysik, Technische Universität Berlin, Hardenbergstr. 36, 10623 Berlin, Germany

⁵¹ Institut für Geologische Wissenschaften, Freie Universität Berlin, 12249 Berlin, Germany

⁵² German Aerospace Center (DLR), Institute of Optical Sensor Systems, Rutherfordstraße 2, 12489 Berlin, Germany

⁵³ Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, UK
Appendix A: List of periodic variabilities in the CHEOPS light curve of HD 172555

<table><caption>**Table A.1.** List of frequencies $f_n$ extracted in the CHEOPS light curve of HD 172555 down to S/N=4.8 and their fitted parameters. The CHEOPS orbital frequency and its harmonics are indicated in the 'Comments' column.</caption><thead><tr><th>Id.</th><th>Frequency (μHz)</th><th>σ<sub>f</sub> (μHz)</th><th>Period (s)</th><th>σ<sub>P</sub> (s)</th><th>Amplitude (%)</th><th>σ<sub>A</sub> (%)</th><th>Phase</th><th>σ<sub>Ph</sub></th><th>S/N</th><th>Comments</th></tr></thead><tbody><tr><td>f<sub>39</sub></td><td>103.41</td><td>0.58</td><td>9670.24</td><td>55</td><td>0.0060</td><td>0.0011</td><td>0.740</td><td>0.060</td><td>5.3</td><td></td></tr><tr><td>f<sub>14</sub></td><td>169.48</td><td>0.17</td><td>5900.40</td><td>6.0</td><td>0.0198</td><td>0.0011</td><td>0.215</td><td>0.018</td><td>18.0</td><td>~ f<sub>orb</sub></td></tr><tr><td>f<sub>29</sub></td><td>193.29</td><td>0.36</td><td>5173.57</td><td>9.7</td><td>0.0093</td><td>0.0011</td><td>0.469</td><td>0.037</td><td>8.5</td><td></td></tr><tr><td>f<sub>7</sub></td><td>310.588</td><td>0.070</td><td>3219.70</td><td>0.73</td><td>0.0460</td><td>0.0010</td><td>0.8054</td><td>0.0072</td><td>44.1</td><td></td></tr><tr><td>f<sub>21</sub></td><td>358.91</td><td>0.27</td><td>2786.21</td><td>2.1</td><td>0.0118</td><td>0.0010</td><td>0.960</td><td>0.028</td><td>11.4</td><td></td></tr><tr><td>f<sub>4</sub></td><td>387.405</td><td>0.050</td><td>2581.28</td><td>0.33</td><td>0.0641</td><td>0.0010</td><td>0.7624</td><td>0.0051</td><td>62.5</td><td></td></tr><tr><td>f<sub>10</sub></td><td>391.669</td><td>0.094</td><td>2553.17</td><td>0.61</td><td>0.0339</td><td>0.0010</td><td>0.6527</td><td>0.0096</td><td>33.2</td><td></td></tr><tr><td>f<sub>24</sub></td><td>396.21</td><td>0.31</td><td>2523.90</td><td>2.0</td><td>0.0103</td><td>0.0010</td><td>0.862</td><td>0.032</td><td>10.1</td><td></td></tr><tr><td>f<sub>8</sub></td><td>409.295</td><td>0.092</td><td>2443.22</td><td>0.55</td><td>0.0341</td><td>0.0010</td><td>0.5535</td><td>0.0094</td><td>33.7</td><td></td></tr><tr><td>f<sub>1</sub></td><td>423.084</td><td>0.023</td><td>2363.60</td><td>0.13</td><td>0.1340</td><td>0.0010</td><td>0.7407</td><td>0.0024</td><td>133.0</td><td></td></tr><tr><td>f<sub>13</sub></td><td>439.82</td><td>0.15</td><td>2273.64</td><td>0.76</td><td>0.0212</td><td>0.0010</td><td>0.504</td><td>0.015</td><td>21.1</td><td></td></tr><tr><td>f<sub>20</sub></td><td>452.69</td><td>0.23</td><td>2209.00</td><td>1.1</td><td>0.01344</td><td>0.00100</td><td>0.869</td><td>0.024</td><td>13.5</td><td></td></tr><tr><td>f<sub>3</sub></td><td>458.138</td><td>0.044</td><td>2182.75</td><td>0.21</td><td>0.07081</td><td>0.00100</td><td>0.1308</td><td>0.0045</td><td>71.1</td><td></td></tr><tr><td>f<sub>11</sub></td><td>465.88</td><td>0.10</td><td>2146.48</td><td>0.47</td><td>0.03035</td><td>0.00099</td><td>0.139</td><td>0.010</td><td>30.6</td><td></td></tr><tr><td>f<sub>25</sub></td><td>474.33</td><td>0.30</td><td>2108.2</td><td>1.3</td><td>0.01021</td><td>0.00099</td><td>0.649</td><td>0.031</td><td>10.3</td><td></td></tr><tr><td>f<sub>6</sub></td><td>480.561</td><td>0.060</td><td>2080.90</td><td>0.26</td><td>0.05146</td><td>0.00099</td><td>0.4813</td><td>0.0061</td><td>52.1</td><td></td></tr><tr><td>f<sub>12</sub></td><td>486.26</td><td>0.10</td><td>2056.52</td><td>0.43</td><td>0.03021</td><td>0.00099</td><td>0.763</td><td>0.010</td><td>30.6</td><td></td></tr><tr><td>f<sub>22</sub></td><td>490.93</td><td>0.26</td><td>2036.90</td><td>1.1</td><td>0.01159</td><td>0.00098</td><td>0.089</td><td>0.027</td><td>11.8</td><td></td></tr><tr><td>f<sub>5</sub></td><td>501.348</td><td>0.058</td><td>1994.62</td><td>0.23</td><td>0.05218</td><td>0.00098</td><td>0.2455</td><td>0.0060</td><td>53.3</td><td></td></tr><tr><td>f<sub>15</sub></td><td>520.98</td><td>0.16</td><td>1919.44</td><td>0.59</td><td>0.01888</td><td>0.00097</td><td>0.450</td><td>0.016</td><td>19.4</td><td></td></tr><tr><td>f<sub>2</sub></td><td>536.714</td><td>0.035</td><td>1863.19</td><td>0.12</td><td>0.08470</td><td>0.00096</td><td>0.2276</td><td>0.0036</td><td>87.9</td><td></td></tr><tr><td>f<sub>9</sub></td><td>548.173</td><td>0.088</td><td>1824.24</td><td>0.29</td><td>0.03397</td><td>0.00096</td><td>0.4118</td><td>0.0090</td><td>35.3</td><td></td></tr><tr><td>f<sub>27</sub></td><td>574.38</td><td>0.31</td><td>1741.01</td><td>0.95</td><td>0.00947</td><td>0.00096</td><td>0.811</td><td>0.032</td><td>9.9</td><td></td></tr><tr><td>f<sub>19</sub></td><td>583.23</td><td>0.20</td><td>1714.58</td><td>0.58</td><td>0.01501</td><td>0.00095</td><td>0.354</td><td>0.020</td><td>15.8</td><td></td></tr><tr><td>f<sub>30</sub></td><td>597.04</td><td>0.32</td><td>1674.92</td><td>0.90</td><td>0.00916</td><td>0.00095</td><td>0.130</td><td>0.033</td><td>9.7</td><td></td></tr><tr><td>f<sub>26</sub></td><td>609.30</td><td>0.30</td><td>1641.22</td><td>0.81</td><td>0.00969</td><td>0.00094</td><td>0.039</td><td>0.031</td><td>10.3</td><td></td></tr><tr><td>f<sub>31</sub></td><td>663.47</td><td>0.33</td><td>1507.24</td><td>0.75</td><td>0.00870</td><td>0.00093</td><td>0.964</td><td>0.034</td><td>9.4</td><td></td></tr><tr><td>f<sub>18</sub></td><td>673.06</td><td>0.19</td><td>1485.76</td><td>0.42</td><td>0.01508</td><td>0.00092</td><td>0.969</td><td>0.019</td><td>16.4</td><td>~ 4 * f<sub>orb</sub></td></tr><tr><td>f<sub>40</sub></td><td>676.16</td><td>0.49</td><td>1478.90</td><td>1.1</td><td>0.00580</td><td>0.00092</td><td>0.056</td><td>0.050</td><td>6.3</td><td></td></tr><tr><td>f<sub>33</sub></td><td>691.33</td><td>0.33</td><td>1446.49</td><td>0.70</td><td>0.00852</td><td>0.00092</td><td>0.904</td><td>0.034</td><td>9.3</td><td></td></tr><tr><td>f<sub>32</sub></td><td>700.44</td><td>0.33</td><td>1427.67</td><td>0.68</td><td>0.00852</td><td>0.00091</td><td>0.054</td><td>0.034</td><td>9.4</td><td></td></tr><tr><td>f<sub>35</sub></td><td>783.29</td><td>0.36</td><td>1276.67</td><td>0.58</td><td>0.00770</td><td>0.00089</td><td>0.718</td><td>0.037</td><td>8.7</td><td></td></tr><tr><td>f<sub>43</sub></td><td>815.85</td><td>0.52</td><td>1225.72</td><td>0.79</td><td>0.00519</td><td>0.00088</td><td>0.739</td><td>0.054</td><td>5.9</td><td></td></tr><tr><td>f<sub>42</sub></td><td>894.11</td><td>0.50</td><td>1118.44</td><td>0.63</td><td>0.00531</td><td>0.00086</td><td>0.326</td><td>0.052</td><td>6.2</td><td></td></tr><tr><td>f<sub>34</sub></td><td>1011.23</td><td>0.31</td><td>988.89</td><td>0.30</td><td>0.00832</td><td>0.00083</td><td>0.782</td><td>0.032</td><td>10.0</td><td>6 * f<sub>orb</sub></td></tr><tr><td>f<sub>44</sub></td><td>1180.09</td><td>0.61</td><td>847.39</td><td>0.44</td><td>0.00399</td><td>0.00079</td><td>0.174</td><td>0.063</td><td>5.1</td><td>7 * f<sub>orb</sub></td></tr><tr><td>f<sub>41</sub></td><td>1685.15</td><td>0.41</td><td>593.42</td><td>0.14</td><td>0.00540</td><td>0.00071</td><td>0.746</td><td>0.042</td><td>7.6</td><td>10 * f<sub>orb</sub></td></tr><tr><td>f<sub>38</sub></td><td>2190.70</td><td>0.30</td><td>456.48</td><td>0.062</td><td>0.00628</td><td>0.00060</td><td>0.122</td><td>0.030</td><td>10.5</td><td>13 * f<sub>orb</sub></td></tr><tr><td>f<sub>36</sub></td><td>2527.68</td><td>0.25</td><td>395.62</td><td>0.039</td><td>0.00670</td><td>0.00054</td><td>0.293</td><td>0.026</td><td>12.3</td><td>15 * f<sub>orb</sub></td></tr><tr><td>f<sub>28</sub></td><td>2696.19</td><td>0.17</td><td>370.89</td><td>0.024</td><td>0.00939</td><td>0.00052</td><td>0.522</td><td>0.018</td><td>18.1</td><td>16 * f<sub>orb</sub></td></tr><tr><td>f<sub>23</sub></td><td>3033.26</td><td>0.13</td><td>329.68</td><td>0.014</td><td>0.01158</td><td>0.00049</td><td>0.670</td><td>0.013</td><td>23.7</td><td>18 * f<sub>orb</sub></td></tr><tr><td>f<sub>16</sub></td><td>3201.824</td><td>0.081</td><td>312.32</td><td>0.0079</td><td>0.01819</td><td>0.00048</td><td>0.9213</td><td>0.0084</td><td>38.1</td><td>19 * f<sub>orb</sub></td></tr><tr><td>f<sub>17</sub></td><td>3370.411</td><td>0.092</td><td>296.70</td><td>0.0081</td><td>0.01562</td><td>0.00046</td><td>0.2217</td><td>0.0094</td><td>33.7</td><td>20 * f<sub>orb</sub></td></tr><tr><td>f<sub>37</sub></td><td>3539.16</td><td>0.22</td><td>282.55</td><td>0.017</td><td>0.00660</td><td>0.00046</td><td>0.482</td><td>0.022</td><td>14.4</td><td>21 * f<sub>orb</sub></td></tr><tr><td>f<sub>54</sub></td><td>3706.50</td><td>1.1</td><td>269.80</td><td>0.077</td><td>0.00124</td><td>0.00044</td><td>0.54</td><td>0.11</td><td>2.8</td><td>22 * f<sub>orb</sub></td></tr><tr><td>f<sub>48</sub></td><td>4886.75</td><td>0.53</td><td>204.63</td><td>0.022</td><td>0.00244</td><td>0.00042</td><td>0.851</td><td>0.055</td><td>5.8</td><td>29 * f<sub>orb</sub></td></tr><tr><td>f<sub>45</sub></td><td>5392.15</td><td>0.44</td><td>185.45</td><td>0.015</td><td>0.00293</td><td>0.00042</td><td>0.188</td><td>0.045</td><td>7.1</td><td>32 * f<sub>orb</sub></td></tr><tr><td>f<sub>50</sub></td><td>5897.95</td><td>0.58</td><td>169.55</td><td>0.017</td><td>0.00216</td><td>0.00040</td><td>0.495</td><td>0.060</td><td>5.3</td><td>35 * f<sub>orb</sub></td></tr><tr><td>f<sub>47</sub></td><td>6740.91</td><td>0.47</td><td>148.35</td><td>0.010</td><td>0.00258</td><td>0.00039</td><td>0.869</td><td>0.048</td><td>6.6</td><td>40 * f<sub>orb</sub></td></tr><tr><td>f<sub>46</sub></td><td>7245.77</td><td>0.42</td><td>138.01</td><td>0.0080</td><td>0.00285</td><td>0.00039</td><td>0.279</td><td>0.043</td><td>7.3</td><td>43 * f<sub>orb</sub></td></tr><tr><td>f<sub>49</sub></td><td>7751.10</td><td>0.52</td><td>129.01</td><td>0.0086</td><td>0.00228</td><td>0.00038</td><td>0.618</td><td>0.053</td><td>6.0</td><td>46 * f<sub>orb</sub></td></tr><tr><td>f<sub>53</sub></td><td>9099.72</td><td>0.64</td><td>109.89</td><td>0.0077</td><td>0.00185</td><td>0.00038</td><td>0.304</td><td>0.066</td><td>4.8</td><td>54 * f<sub>orb</sub></td></tr></tbody></table>