# Glancing through the debris disk: Photometric analysis of DE Boo with CHEOPS

Á. Boldog, Gy. M. Szabó, L. Kriskovics, A. Brandeker, F. Kiefer, A.

Bekkelien, P. Guterman, G. Olofsson, A. E. Simon, D. Gandolfi, et al.

► To cite this version:

Á. Boldog, Gy. M. Szabó, L. Kriskovics, A. Brandeker, F. Kiefer, et al. Glancing through the debris disk: Photometric analysis of DE Boo with CHEOPS. Astronomy and Astrophysics - A&A, 2023, 671, 10.1051/0004-6361/202245101. insu-04155764

HAL Id: insu-04155764

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

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
# Glancing through the debris disk: Photometric analysis of DE Boo with CHEOPS*

Á. Boldog¹,²,³ Ⓓ, Gy. M. Szabó⁴,³, L. Kriskovics¹,² Ⓓ, A. Brandeker⁵, F. Kiefer⁶,⁷ Ⓓ, A. Bekkelien⁸ Ⓓ, P. Guterman⁹,¹⁰, G. Olofsson⁵, A. E. Simon¹¹ Ⓓ, D. Gandolfi¹² Ⓓ, L. M. Serrano¹² Ⓓ, T. G. Wilson¹³ Ⓓ, S. G. Sousa¹⁴ Ⓓ, A. Lecavelier des Etangs¹⁵ Ⓓ, Y. Alibert¹¹ Ⓓ, R. Alonso¹⁶,¹⁷ Ⓓ, G. Anglada¹⁸,¹⁹, T. Bandy²⁰, T. Bárczy²¹ Ⓓ, D. Barrado²² Ⓓ, S. C. C. Barros¹⁴,²³ Ⓓ, W. Baumjohann²⁴ Ⓓ, M. Beck⁸, T. Beck¹¹, W. Benz¹¹,²⁵, N. Billot⁸, X. Bonfils²⁶, C. Broeg¹¹,²⁷, M. Buder²⁸, J. Cabrera²⁹, S. Charnoz³⁰, A. Collier Cameron¹³, C. Corral van Damme³¹, Sz. Csizmadia²⁹ Ⓓ, M. B. Davies³², A. Deline⁸, M. Deleuil⁹, L. Delrez³³,³⁴, O. D. S. Demangeon¹⁴,²³, B.-O. Demory²⁵ Ⓓ, D. Ehrenreich⁸ Ⓓ, A. Erikson²⁹, J. Farinato³⁵, 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. Mecina⁴⁵ Ⓓ, V. Nascimbeni⁴³, R. Ottensamer⁴⁵ Ⓓ, I. Pagano³⁵ Ⓓ, E. Pallé¹⁶ Ⓓ, G. Peter²⁸, G. Piotto⁴³,⁴⁶, D. Pollacco³⁹, D. Queloz⁴⁷,⁴⁸, R. Ragazzoni⁴³,⁴⁶, N. Rando³¹, H. Rauer²⁹,⁴⁹,⁵⁰ Ⓓ, I. Ribas¹⁸,¹⁹ Ⓓ, N. C. Santos¹⁴,²³ Ⓓ, G. Scandariato³⁵ Ⓓ, D. Ségransan⁸, A. M. S. Smith²⁹ Ⓓ, M. Steller²⁴, N. Thomas¹¹, S. Udry⁸ Ⓓ, V. Van Grootel³⁴ Ⓓ, and N. A. Walton⁵¹ Ⓓ

*(Affiliations can be found after the references)*

Received 30 September 2022 / Accepted 20 January 2023

## ABSTRACT

**Aims.** DE Boo is a unique system, with an edge-on view through the debris disk around the star. The disk, which is analogous to the Kuiper belt in the Solar System, was reported to extend from 74 to 84 AU from the central star. The high photometric precision of the Characterising Exoplanet Satellite (CHEOPS) provided an exceptional opportunity to observe small variations in the light curve due to transiting material in the disk. This is a unique chance to investigate processes in the debris disk.

**Methods.** Photometric observations of DE Boo of a total of four days were carried out with CHEOPS. Photometric variations due to spots on the stellar surface were subtracted from the light curves by applying a two-spot model and a fourth-order polynomial. The photometric observations were accompanied by spectroscopic measurements with the 1m RCC telescope at Piszkéstető and with the SOPHIE spectrograph in order to refine the astrophysical parameters of DE Boo.

**Results.** We present a detailed analysis of the photometric observation of DE Boo. We report the presence of nonperiodic transient features in the residual light curves with a transit duration of 0.3–0.8 days. We calculated the maximum distance of the material responsible for these variations to be 2.47 AU from the central star, much closer than most of the mass of the debris disk. Furthermore, we report the first observation of flaring events in this system.

**Conclusions.** We interpreted the transient features as the result of scattering in an inner debris disk around DE Boo. The processes responsible for these variations were investigated in the context of interactions between planetesimals in the system.

**Key words.** techniques: photometric – methods: data analysis – stars: flare – starspots – circumstellar matter

## 1. Introduction

Debris disks have been observed outside of our Solar System since the 1980s (Aumann et al. 1984; Smith & Terrile 1984). In contrast to primordial disks of gas and dust, which dissipate after a few million years of the formation of the star, dusty debris disks can be observed for a longer period even around K type stars (Moór et al. 2006). These disks are made of dust that experiences frequent collisions with larger bodies such as comets or planetesimals along their orbits. The semi major axes of debris disks differ from system to system. Some debris disks are found as close to the central star as around 1 AU, such as in the HD 69830 system (Lisse et al. 2007), while others can stretch out several hundreds of AU (Smith & Terrile 1984; Kalas et al. 2007). Systems with both an inner and an outer debris disk have also been reported (Lisse et al. 2008; Stark et al. 2009).

High-sensitivity photometry can provide unique information on the structure of an edge-on debris disk. When disk inhomogeneities pass in front of the star, they produce variations in the star's extinction, dimming the stellar light. These variations are interpreted as transits of structures in the circumstellar disk (Lecavelier Des Etangs et al. 1999).

During its nominal mission, TESS also observed quasi-periodic and quasi-stochastic transients in “dipper” stars, which are surrounded by debris disks (Gaidos et al. 2019). The observed features show a quasi-continuous spectrum of separated, non-recurrent random transits to quasi-periodic events with a period of 1–4 days. In the case of HD 240779 for example, the dimming vanished over 2 yr between two TESS visits, suggesting that the origin of the transient was disrupted or underwent efficient evaporation (Gaidos 2022). Most of these dippers are early-type stars, while K-M stars can also be found between them. Another famous example for an edge-on disk is the β Pic system

<footnote>* This article uses data from CHEOPS programme CH_PR100010.</footnote>
<table><caption>**Table 1.** Stellar parameters of the K0V primary (A) and the M dwarf companion (B) from the literature and this work which were taken as inputs to the present analysis.</caption><thead><tr><th>Parameter</th><th>Value</th><th>Reference</th></tr></thead><tbody><tr><td><i>M</i><sub>A</sub></td><td>0.84 <i>M</i><sub>☉</sub></td><td>Marshall et al. (2014)</td></tr><tr><td><i>R</i><sub>A</sub></td><td>0.86 <i>R</i><sub>☉</sub></td><td>Marshall et al. (2014)</td></tr><tr><td><i>L</i><sub>A</sub></td><td>0.498 <i>L</i><sub>☉</sub></td><td>Marshall et al. (2014)</td></tr><tr><td><i>P</i><sub>rot,A</sub></td><td>10.39 ± 0.03 days</td><td>Henry et al. (1995)</td></tr><tr><td><i>M</i><sub>B</sub></td><td>0.45 <i>M</i><sub>☉</sub></td><td>Kennedy (2014)</td></tr><tr><td><i>R</i><sub>A</sub></td><td>0.870 ± 0.021 <i>R</i><sub>☉</sub></td><td>This work</td></tr><tr><td><i>R</i><sub>B</sub></td><td>0.37 ± 0.17 <i>R</i><sub>☉</sub></td><td>This work</td></tr><tr><td><i>T</i><sub>eff,B</sub></td><td>3596 ± 124 K</td><td>This work</td></tr><tr><td><i>L</i><sub>B</sub></td><td>0.021<sup>+0.029</sup><sub>-0.019</sub> <i>L</i><sub>☉</sub></td><td>This work</td></tr></tbody></table>

**Notes.** For a detailed description of the estimation of stellar parameters calculated in this present analysis, see Sect. 3.1.

(Smith & Terrile 1984), where transient dips have been observed (Zieba et al. 2019; Pavlenko et al. 2022; Lecavelier des Etangs et al. 2022), with a flux drop at the level of 0.03% at maximum and a duration of less than 1 day, and with no periodicity seen in the sectors investigated. Unlike the dipper stars, the transients in the β Pic system were attributed to exocomets.

Ground-based facilities are limited in their ability to observe these variations because of the Earth's atmosphere and the diurnal observation gaps. We started a space-based photometry program with the Characterising Exoplanet Satellite (CHEOPS) to detect these tiny variations in a few selected systems. The CHEOPS space telescope launched in 2019 and designed to detect and characterise the transits of small exoplanets (Benz et al. 2021), is well placed to perform this search because of its high photometric precision (Maxted et al. 2021).

DE Bootis (HD 131511) is a chromospherically active spectroscopic binary (SB 1) dwarf RS CVn system as identified by Eker et al. (2008), with a K0V primary component (Gray et al. 2003). The stellar fundamental parameters for the primary and the M dwarf companion are shown in Table 1. The binary has a separation of 0.19±0.03 AU, an orbital period of 125.396 days with the eccentricity of the orbit being *e* = 0.51 ± 0.01 (Kennedy 2014). The age of the system was calculated by Gray et al. (2015) to be less than 797 Myr, while Marshall et al. (2014) found it to be between 570 and 700 Myr. DE Boo exhibits signs of magnetic activity, such as emission in the Ca II H and K lines and large stellar spots. These spots can be observed in the rotational modulation of the light curve with a period of 10.39 days (Henry et al. 1995). A radially narrow debris disk at around 70 AU from the primary component has been reported by Marshall et al. (2014). DE Boo is among the closest systems with a debris disk, lying at 37.928±0.243 light years from Earth (Gaia Collaboration 2020), while its Gaia DR2 magnitude is 5.766, leading to a 19.92 ppm internal precision per hour with CHEOPS. The disk, analogous to the Solar System's Kuiper-belt, is aligned with the orbit of the host binary and has a near edge-on geometry, as derived by Jancart et al. (2005) using HIPPARCOS data. Furthermore, Marshall et al. (2014) showed that the inclination of the disk is at least 84°. Due to this, the debris disk is always in the line of sight, continuously causing dips in the light curve.

The level of the expected variations can be estimated via the value of $f = L_{disk}/L_{star}$, which implies that $f = 4 \times 10^{-6}$ of the starlight is intercepted by the dust. The radial optical depth is $\tau = f/\sin(\theta/2)$, where $\theta$ is the opening angle

(Kennedy et al. 2014). The resulting $\tau \approx 1 \times 10^{-5}$ is roughly the level of variation expected if there are radial holes in the dust distribution. Because of the frequent collisions between the smallest bodies inside the disk and the Keplerian shear, such holes are unexpected. However, a similar process was suggested in the case of dippers (Tajiri et al. 2020). The transit of planetsimals that are collecting dense dust clumps inside their Hill radius can still result in detectable light variations. The viewing geometry of DE Boo along with the high expected S/N make this system an ideal target for the Dusty Debris Disk (DDD) survey with CHEOPS, and we have chosen this system as one of the targets of our program.

The co-planarity of the debris disk with the orbit of the binary was studied by Kennedy (2014). The favorable orientation allows us to examine the behavior of the material in the debris disk – such as extrasolar collisions – early in the evolution of the system. The outstanding precision of CHEOPS makes DE Boo a prime target for seeking light-curve anomalies from an edge-on debris disk.

Here we present photometric observations carried out with the CHEOPS space telescope along with spectroscopic measurements with the Hungarian Im telescope at Piszkéstető Mountain Station and with the SOPHIE spectrograph at the Observatoire de Haute-Provence (OHP). In Sect. 2, we describe our observations and the data-processing methods we used. In Sect. 3, we analyze the light curve and present our results. In Sect. 4, we derive our conclusions and provide a possible explanation of the approximately half-day transient features seen in the residual light curve of DE Boo.

## 2. Observations and data processing

During the 2021 visibility window, DE Boo was observed twice with CHEOPS, in May and June. Simultaneous spectroscopic measurements were also gathered from Piszkéstető Mountain Station, Hungary.

### 2.1. Photometric observations with CHEOPS

Observation windows with a length of 30 CHEOPS orbits (1 orbit = 98.77 min) were set up for both visits in the CHEOPS Proposal Handling Tool (PHT). The total time of the observation was four days. Observation logs are shown in Table 2. As DE Boo is an active star with a brightness of $G = 5^{m}766$ in the Gaia G-band (Gaia Collaboration 2020), we used a short exposure time of 3 s in order to be able to resolve possible flares and to reduce the effect of saturation. The nominal efficiency of the observations was 71% and 68% during the first and second visits to DE Boo. The following ephemeris was used to phase the observations:

$$JD = 2459335.084 + 10.39 \times E, \qquad (1)$$

where 10.39 days was the rotational period of DE Boo.

Photometry was performed using the *imagettes* for each exposure. The imagettes are small images centered on the target with a radius of 30 pixels. The advantage of these images is that because they are smaller in size, they do not have to be co-added on board as in subarrays, and can be downloaded individually. Using the imagettes enabled us to achieve better time resolution, which is crucial in order to adequately deal with flares. Aperture photometry, on the other hand, which is provided by the CHEOPS Data Reduction Pipeline (DRP, Hoyer et al. 2020) for subarrays, becomes more difficult in the case of imagettes
<table><caption>**Table 2.** Logs of 2021 CHEOPS observations of DE Boo.</caption><thead><tr><th>Visit #</th><th>Start date</th><th>End date</th><th>File key</th><th>CHEOPS product</th><th>Integ. time (s)</th><th>Num. of frames</th></tr></thead><tbody><tr><td>1</td><td>2021-05-31 17:56:17</td><td>2021-06-02 18:53:04</td><td>PR100010_TG001101</td><td>Imagettes</td><td>3</td><td>2803</td></tr><tr><td>2</td><td>2021-06-04 20:52:07</td><td>2021-06-06 21:57:03</td><td>PR100010_TG001102</td><td>Imagettes</td><td>3</td><td>2692</td></tr></tbody></table>

**Notes.** The time notation follows the ISO-8601 convention. The File Key supports the fast identification of the observations in the CHEOPS archive.

because of the small sizes. To overcome this issue, we used PIPE<fnref n="1" /> (PSF imagette photometric extraction; Brandeker et al. in prep.; see also descriptions in Szabó et al. 2021; Morris et al. 2021b), a tool specifically developed for photometric extraction of imagettes using point-spread function (PSF) photometry. This allowed us to derive photometry consistent with the DRP and with comparable S/N, account for instrumental effects such as smearing, and take advantage of the shorter cadence of the imagettes, thus providing help with analyzing flares.

To deal with contamination, frames with recognized issues (e.g., strong background, cosmic rays) were dismissed. Only images unaffected by these events were further used for data analysis.

## 2.2. Spectroscopic observations and analysis

A total of 95 spectra were gathered with the echelle spectrograph (R = 21 000) mounted on the 1 m RCC telescope at Piszkéstető Mountain Station, Hungary, between 29 May and 2 June 2021, in parallel with the CHEOPS observations, with exposure times of 300 s and 600 s. The former yielded an average S/N of 83, while with the latter an average S/N of 160 was reached at 6400 Å. Data reduction was carried out with the regular IRAF<fnref n="2" /> echelle tasks. Wavelength calibration was carried out using ThAr calibration spectra taken between observations. We also had five observations from the SOPHIE spectrograph mounted on the 1.93 m telescope at the Observatoire de Haute-Provence, taken between 31 May and 5 June with an exposure time of 600 s and a spectral resolving power of R = 75 000, yielding an average S/N of 255.

In order to determine precise astrophysical parameters, spectral synthesis was carried out using SME (Piskunov & Valenti 2017) on both datasets. During the synthesis, MARCS models were used (Gustafsson et al. 2008). Atomic line parameters were taken from the VALD database (Kupka et al. 1999). Macroturbulence was estimated using the following equation (Valenti & Fischer 2005):

$$v_{\text{mac}} = \left( 3.98 - \frac{T_{\text{eff}} - 5770 \text{ K}}{650 \text{ K}} \right) \text{km s}^{-1}. \quad (2)$$

The whole process is described in Kriskovics et al. (2019). The resulting astrophysical parameters for both the Piszkéstető and OHP datasets are summarized in Table 3, and are in relatively good agreement. The slight difference in the $v \sin i$ values could be attributed to the fact that lower resolution results in wider instrumental profiles, which are not always easy to disentangle. We also note that based on the flux ratio of the two components (see Table 1), the contribution of the secondary component to the spectra is negligible.

<table><caption>**Table 3.** Astrophysical parameters from spectral synthesis for the Piszkéstető and OHP datasets.</caption><thead><tr><th></th><th>Piszkéstető</th><th>OHP</th></tr></thead><tbody><tr><td>$T_{\text{eff}}$</td><td>5290 ± 93 K</td><td>5245 ± 40 K</td></tr><tr><td>$\log g$</td><td>4.47 ± 0.15</td><td>4.5 ± 0.05</td></tr><tr><td>[Fe/H]</td><td>0.09 ± 0.1</td><td>0.14 ± 0.07</td></tr><tr><td>$v_{\text{mic}}$</td><td>1.3 ± 0.9 km s<sup>-1</sup></td><td>1.0 ± 0.1 km s<sup>-1</sup></td></tr><tr><td>$v_{\text{mac}}$</td><td>4.7 km s<sup>-1</sup></td><td>4.7 km s<sup>-1</sup></td></tr><tr><td>$v \sin i$</td><td>7 ± 1 km s<sup>-1</sup></td><td>4.5 ± 0.5 km s<sup>-1</sup></td></tr></tbody></table>

## 2.3. Light-curve pre-processing

Further corrections were carried out in the photometric data to account for the effects of the systematic errors on roll angle. These errors were corrected following a nonparametric estimation from the data, whereby photometric points were phased according to the roll angle. The data points were smeared in the time domain first, in a boxcar of ten data points in length. After phasing the data in the roll-angle domain, a 2.5 sigma-clipping was applied to omit the outliers, such as possible flares and light-curve transients. The phased light-curve points were then averaged in roll-angle bins of 3° in width, and a prediction of the roll-angle effect was made for all measured points with a linear interpolation from the binned pattern.

During the first visit, the star showed flaring activity on two occasions. We plotted these flares in magnification, and find that the noise in the individual data points is compatible with the full amplitude of the flares. This can hide the internal structure of flares if there are multiple flares in succession. To be able to resolve the flare complexes, the light curve was resampled with a cadence of 15 s, which meant binning five data points together. This gave the most informative resolution of the flare structures. We also kept this binning in the spot modeling and the analysis of the residuals.

## 2.4. Stellar spot model

In order to compensate for the rotational modulation, an analytic model of two spots was fitted with SpotModel (Ribárik et al. 2003), following Budding (1977). The code assumes homogeneous, circular spots with radii, longitudes, and latitudes treated as free parameters, and iteratively changes these parameters until the computed synthetic light curve corresponding to a given spot configuration adequately fits the observed light curve. Assuming $\Delta T \approx 1500$ K (which can be a typical value for cool dwarfs; e.g. Solanki 2003) for the spots and using the effective temperature from the spectral synthesis, a spot intensity of 0.25 was used. The rotational period of 10.39 days was kept constant during the spot modeling. The fit is shown in the upper panel of Fig. 1,

<footnote n="1">¹ https://github.com/alphapsa/PIPE</footnote>

<footnote n="2">² iraf.net</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Analytical spot model fitted to our CHEOPS data. Upper panel: CHEOPS light curve (black dots) with the fitted analytic spot model (blue solid line). Dashed lines denote the longitudes of the fitted spots while dash-dotted lines indicate the radii of these spots. Middle panel: residual after subtraction of the spot model. Lower panel: residual after subtraction of the low-order polynomial. Blue ticks denote the positions of the two observed flares. The light curve is phased following Eq. (1) with $P_{rot}$ = 10.39 days.</figcaption></figure>

and the resulting spot parameters are summarized in Table 4. As spot latitudes, sizes, and intensities can mutually compensate for each other, latitudes and sizes should be considered nominal. Due to this degeneracy, different spot latitudes, sizes, and temperatures (and therefore spot intensities) can yield virtually the same fit, but as our primary goal is to remove rotational modulation caused by the spots, and because this does not affect the residuals, the fitted model can definitely be used to compensate for the rotational modulation. In the case of high-precision photometry, it is often hard to perfectly fit the light curve with a relatively low number of analytic spots, and therefore it is possible that a slight trend remains (middle panel in Fig. 1). As it is clearly seen in Fig. 1 that the trends are the product of a slight misfit, a further subtraction of a fourth-order polynomial from the residual data was carried out. To ensure no higher harmonics of the stellar rotation were present after the subtraction of the spot model, we studied the periodogram of the data before fitting the polynomial. The resulting periodogram can be seen in Fig. A.1.

In order to choose the most appropriate model, we performed a Bayesian information criterion (BIC) check over several

<table><caption>**Table 4.** Spot parameters from the analytical spot model. *l*, *b*, and *r* denote spot longitudes, latitudes, and radii, respectively.</caption><thead><tr><th></th><th>Spot 1</th><th>Spot 2</th></tr></thead><tbody><tr><td><i>l</i> (deg)</td><td>12.17 ± 0.05</td><td>160.34 ± 0.11</td></tr><tr><td><i>b</i> (deg)</td><td>-10.16 ± 5.30</td><td>83.60 ± 0.05</td></tr><tr><td><i>r</i> (deg)</td><td>11.42 ± 0.21</td><td>40.12 ± 0.03</td></tr></tbody></table>

polynomial models. The fourth-order polynomial model had the lowest BIC value and was used for fitting and further subtraction. During this step, we removed the variations on longer timescales (≥1 day characteristic period), and the order of the applied polynomial has a direct and significant influence on this timescale. In the middle panel of Fig. 1, we indeed see variations in this timescale, and these features are very likely not instrumental. CHEOPS is known to be satisfactorily stable in this sense (see e.g., Delrez et al. 2021 or Morris et al. 2021a for the various instrumental effects of CHEOPS). Also, there are
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Residual light curves of the first (upper panel) and second (lower panel) visit to DE Boo after the subtraction of a fourth-order polynomial from the spot-model-subtracted data. Red dots show the orbital average flux for each orbit of CHEOPS. Red curves emphasise the change in the brightness of the target.</figcaption></figure>

no telemeric covectors similar to the observed pattern. However, this timescale does include the light variations caused by the fine structure of the spots on the rotating surface. As we do not have a sufficiently detailed spot model, we cannot make conclusions as to the origin of the long-timescale residuals, or decipher whether they are of a transient nature or are mostly residuals from the rotational signal. Therefore, we decided not to interpret this variation, and removed this signal to reveal the variations on shorter timescales. Here, our aim is not the exact reconstruction of all signals from the disk, but to determine the signal in the data that clearly cannot be of stellar or instrumental origin. The applied process supported this aim by removing the longer variations seen in the light curve with characteristic periods on the order of days, without affecting the variations on shorter timescales (lower panel of Fig. 1), which we discuss in the following sections. Residual light curves after the subtraction of the polynomial are shown in Fig. 2.

## 2.5. Estimating flare energies

During the first visit to DE Boo, two flares were detected by CHEOPS. This is the first time that flares have been photometrically observed and reported for DE Boo. This unprecedented detection justified a detailed analysis of these phenomena. In order to do so, we needed imagettes of sufficiently short cadence to resolve the finer structures of the flares. However, the resampling of these data for 15 s was also necessary to increase the S/N. The flares in the re-sampled imagette data are shown in Fig. 3.

In order to obtain the fraction of DE Boo's quiescent luminosity in the CHEOPS band, a BT-NextGen model spectrum for a K0V-type star with $T_{\text{eff}} = 5300$ K and $\log g = 4.5$ was convolved with the transmission curve of CHEOPS (Deline et al. 2020). Flare energies were then calculated by integrating the

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Two flares that occurred during the first visit to DE Boo, after re-sampling the imagette data to a cadence of 15 s.</figcaption></figure>

total flux for the duration of the flaring events, and multiplying them with the bolometric luminosity of DE Boo and the fraction of the luminosity in the CHEOPS band. Though we assumed that the flares originate from the primary, in reality they could arise from the M dwarf companion as well. However, as both the quiescent luminosity and the stellar spectrum are dominated by the K0V primary component, our approach is expected to provide a reasonable estimation of the flare energies. We note that the order of magnitude of the energies would not be different even if they originated from the companion.

# 3. Results

## 3.1. Improved radius of DE Boo

We improved existing measurements of the radius of DE Boo using a Markov-chain Monte Carlo infrared flux method (MCMC IRFM; Blackwell & Shallis 1977; Schanche et al. 2020). This approach uses known relations between the stellar angular diameter, effective temperature, and apparent bolometric flux. For DE Boo, we computed the bolometric flux by constructing a spectral energy distribution (SED) of two stellar components: one of the primary using the spectral parameters derived above with the ATLAS catalogs (Castelli & Kurucz 2003) and a second of the known M-dwarf companion with parameters taken from the PHOENIX catalogs (Allard et al. 2011). Empirical relations from Baraffe et al. (2015) were used to estimate the effective temperature and the radius of the companion. This was necessary in order to derive an accurate stellar radius for DE Boo using the derived parameters and the offset-corrected Gaia EDR3 parallax (Lindegren et al. 2021). We fitted the combined models to observed data taken from the most recent data releases for the following bandpasses: Gaia G, $G_{\text{BP}}$, and $G_{\text{RP}}$, 2MASS J, H, and K, and WISE W1 and W2 (Skrutskie et al. 2006; Wright et al. 2010; Gaia Collaboration 2021). We used this fit to derive $R_{\star} = 0.870 \pm 0.021 R_{\odot}$.
<table><caption>**Table 5.** Estimated energies and durations of the two flares during the first visit to DE Boo.</caption><thead><tr><th>Start date<br>JD-2400000</th><th>End date<br>JD-2400000</th><th>Duration<br>(min)</th><th>Estimated energy<br>(erg)</th></tr></thead><tbody><tr><td>59366.273541</td><td>59366.287429</td><td>20</td><td>3.35 × 10<sup>32</sup></td></tr><tr><td>59366.830820</td><td>59366.834061</td><td>4.7</td><td>7.43 × 10<sup>31</sup></td></tr></tbody></table>

## 3.2. Properties of flare activity

We estimated the energies of the two flares observed during the first visit. The first, more energetic flare, which lasted for 20 min, occurred at the beginning of the observation period and had an estimated energy of $3.35×10^{32}$ erg. In contrast, the second flare had a much shorter duration of 5 min and an energy of about $7.43×10^{31}$ erg and was only identifiable on the imagettes. The detection of two flares with such energies under the total observation time of four days is not out of the ordinary on a fairly active late G or early K dwarf (e.g. Leitinger et al. 2020, and references therein). The energy levels are relatively typical as well (Vida et al. 2021). Assuming that these flares originated from the companion, their energies would still be well within the boundaries of what is expected in the case of an M dwarf star. However, we note that two flares during a four-day period is not enough to make any quantitative conclusions about the flare activity. Derived parameters of the two flares along with their dates of occurrence are shown in Table 5.

## 3.3. Possible photometric transients in the light-curve residuals

After masking out the flares, subtracting a spot model, and removing slowly varying systematic errors due to possible structures on the stellar surface that were not accounted for in our simple two-spot model, we noticed the presence of moderately slowly varying light-curve residuals on a timescale of 0.3–0.8 days. These residuals have a shorter timescale than what could be originated by the rotation of spots (and the timescale of spot emergence or decay is much longer than this; see e.g. Namekata et al. 2019), while the timescale is too long to be due to flares or other activity transients. Other typical light curve changes, such as ellipticity in closer binaries, or other tracers of magnetic activity typical of RS CVn systems are all related to rotation, and in our case would cause changes on longer timescales (see e.g. Berdyugina 2005, and references therein). Therefore, we interpret these residuals as transients from the edge-on disk in the following.

Let us assume a cylindrically symmetric disk consisting of a light scattering material that follows some $\sigma(r, \phi)$ distribution. Here, $\sigma$ is the scattering cross section of the disk locally, that is, the fraction of starlight that is occulted or scattered by the dust particles in the elementary volume; $r$ is the radius vector, and $\phi$ is the angular coordinate. If the dust follows a circular orbit, the $\omega$ angular velocity can be expressed as $\omega = \sqrt{2GMr^{-3/2}}$, according to Kepler's Third Law. Hence, during a dt time, we can write

$$d\phi = \omega dt = \sqrt{2GM} r^{-\frac{3}{2}} dt. \qquad (3)$$

The variation of the disk scattering along the line of sight is then

$$\frac{d\sigma}{dt} = \frac{d\sigma}{dr} \frac{dr}{dt} + \frac{d\sigma}{d\phi} \frac{d\phi}{dt} = \frac{d\sigma}{d\phi} \sqrt{2GM} r^{-\frac{3}{2}}, \qquad (4)$$

because $\frac{dr}{dt} = 0$. The variation of the total scattering cross section along the line of sight $\Sigma = \int_{r1}^{r2} \sigma dr$ can then be expressed as

$$\frac{d\Sigma}{dt} = \int_{r1}^{r2} \frac{d\sigma}{d\phi} \sqrt{2GM} r^{-\frac{3}{2}}, dr. \qquad (5)$$

This formula represents the amount of material along a distinct line of sight. The light-curve transients are related to the transits of the individual grains along the stellar disk. If a material at $r$ radius with $\rho$ density scatters out a given fraction of the light $S$, the contribution to the light curve can be written as

$$S(r, t) = [\rho(r, \tau) * Lc(r, \tau)](t), \qquad (6)$$

where $Lc$ is the transit light curve of a compact transiter at $r$ radius, * denotes convolution, and $\tau$ is a running time-like variable, which was introduced because the transit light curve is a time-dependent variable itself. This expression shows that the shape of the transit light curve can be derived as the shapes of all individual overdensities at their $r$ radius and $\tau$ transit times combined appropriately and evaluated at $t$. If the configuration is very thin, and the scattered light is on the order of one part per thousand at most, the combination can be simplified to a summation. In this case, the combined effect of all contributions to the light curve shape will be written as

$$S(t) = \int_{R1}^{R2} [\rho(r, \tau) * Lc(r, \tau)](t) dr. \qquad (7)$$

Here, the convolution shows that the frequency spectrum of the transit light curve will be the $F(\sigma)$ frequency spectrum belonging to the material distribution multiplied by the $F(T)$ frequency spectra of the transit light curves at different radii. If the cloud of grains is compact (the transiting structure is much smaller than the stellar disk) and $F(\sigma)$ includes components at very large frequencies, the boundary frequency in $F(\text{signal})$ will be the highest frequency in $F(T)$. If the transiting feature is compact, for example a small planet, the duration of the transient will be close to the transit duration. If the transiting feature is extended, for example a loose cloud of large asteroids or a giant comet, where the ingress and egress duration are compatible with (or even exceed) the transit duration, the experienced timescale of the transients will be significantly longer than the transit duration calculated at the orbital semi-major axis of the transiting feature.

If the transiting cloud is extended far beyond the stellar radius, ($\sigma$) has only low-frequency components, and the observed light-curve anomaly will be very slow, on the order of $\Delta\phi/\sqrt{(2GM)r^{-3/2}}$. In this case, the smoothness of the distribution of the material and the orbital radius of the disk features become degenerate. Therefore, giving a lower limit for the orbital radius of the excess material, we also have to give an upper limit for its angular extension. Similar descriptions of these processes can be found in Kennedy et al. (2017).

In summary, the transit of the grains that generate a light-curve anomaly can be at most as long as the anomaly itself. Applying Kepler's Third Law, this gives us an estimate, or more precisely, an upper limit of the radius where the claimed transiting features can be found, which we can write as

$$a_{\text{max}} = \frac{GM_*D^2}{4(1-b^2)R_*^2}, \qquad (8)$$
where *G* is the gravitational constant, *M*<sub>★</sub> and *R*<sub>★</sub> are the mass and radius of DE Boo, respectively, *D* is the duration of the light-curve anomaly itself, measured in units of time, and *b* is the impact parameter, the projected distance of the transit chord from the center point of the star, measured in units of *R*<sub>★</sub>. If the shape of the transient is compatible with the transit of a compact, planet-like object, *D* is related to the relative duration *W*, as *D* = *W* × *P*<sub>orb</sub>, where *W* = *R*<sub>★</sub>/a√(1 + *k*)<sup>2</sup> − *b*<sup>2</sup>/π. If the light variation is more fuzzy and we cannot define a major transit bin, *D* is to be understood as simply the time from the beginning to the end of the negative transit anomalies.

The residuals of the DE Boo observations are plotted in Fig. 2. The residuals of the two visits differ from each other. The first visit is compatible with a stable residual, with the dips suggesting one or two transients passing in front of the star. The most stable part of the residual is observed between 2459367.3 and 2459368.3. The second CHEOPS visit is merely composed of varying residuals and suggests density variations around the star rather than one transiting object. The most significant deviations from a constant residual are observed between 2459367.0 and 2459367.3, 2459370.8 and 2459371.6, and 2459371.75 and 2459372.35. Therefore, the experienced timescale of the suspected transients in the DE Boo light curve is 0.3–0.8 days.

Using the stellar parameters given in Table 1 and inserting 0.3 and 0.8 days for the transit duration into Eq. (8) and assuming *b* = 0, we get 0.35 and 2.47 AU for the maximum orbital distance of the transiting feature. This places the transiting object much closer to the star than the debris disk, to indicatively 1 AU from DE Boo. The suspected structure is consistent with the criteria of dynamical stability according to Holman & Wiegert (1999). This criterion tells us that circumbinary structures farther from the system center than the *a*<sub>c</sub> critical radius are stable in the long term; whereas values of *a*<sub>c</sub> were found by means of numerical simulations and as a function of the *µ* mass ratio and the eccentricity of the binary. Table 7 in Holman & Wiegert (1999) shows us that near *µ* = 0.4–0.5, *e* = 0.5 and *a*<sub>c</sub> = 3.6–3.7 if measured in the unit of the primary-to-secondary separation. As this separation is 0.19 AU, debris structures behind the critical radius of *a*<sub>c</sub> = 0.703 AU can be dynamically stable. Therefore, a transient at 1 AU from the barycenter is dynamically plausible. We also note that Trilling et al. (2007) review several other dynamical criteria, but these authors find the Holman & Wiegert (1999) criterion to be the strictest of these, and therefore we applied the Holman & Wiegert (1999) criterion here.

**4. Discussion and conclusions**

We carried out photometric observations of DE Boo using the CHEOPS space telescope with a total duration of four days. In order to properly remove the effect of stellar variation, a stellar spot model and a fourth-order polynomial were fitted and subtracted from the light curves, respectively. We observed two flares during the observation period. This was the first time that flaring activity of DE Boo was reported. We found the estimated energies of the flares to be consistent with the predicted levels of activity of a magnetically active K dwarf star. We conclude that further measurements are required to obtain quantitative results about the stellar activity of DE Boo.

Local thermodynamic equilibrium synthetic spectra were
fitted to the observations taken with the echelle spectro-
graph mounted on the 1m RCC telescope at Piszkéstető,
Hungary, and with SOPHIE installed on the 1.93 m tele-
scope at the Observatoire de Haute-Provence, France, in

order to refine the astrophysical parameters of the primary.
The two fits yielded the same results within their respective
error bars.

Nonperiodic transiting features with a duration of 0.3–0.8 days were identified in the residual light curves of DE Boo. Based on their transit duration, the semi-major axes of these features were estimated. Our calculations show that these features had a maximum orbital distance of 2.47 AU. This is a much smaller value than the indicated inner radius of the debris disk that is seen in the infrared wavelengths. As the duration of our observations was less than the expected transit duration of the material passing in the debris disk at 70 AU according to Eq. (8) (around 4.5 days), it is not possible to detect their effects on the investigated light curves. Space telescopes with longer observational baselines (such as TESS) are better suited to this task. Moreover, we note that any effect with the duration of a few days was removed during the polynomial fitting. However, the goal of this study was not simply to identify features of an already known debris disk. As it is not out of the ordinary that an inner disk is found in a system with outer debris disks, it cannot be excluded that, for DE Boo, transiting material is observed much closer to the star than the most massive parts of the debris disk. In the case of β Pic, we see a similar geometry: although the visible debris disk extends to 100 (secondary disk) and 200 (primary disk) AU, the observed transients last for less than a day, indicating that the presumed exocomets orbit much closer than the debris disks. Several systems with warm debris disks have also been known to exhibit quasi-periodic dippers with a timescale of 1–4 days; these were found to be due to transits of overdensities in the disk, of surprisingly large amplitude on some occasions (Tajiri et al. 2020). These systems are more similar to DE Boo in terms of the origin of the light variations. In these systems, the overdensities also appear in a structure that is close to the star (as can be concluded from their quasi-period), while the major warm debris disk lies much farther out.

A possible explanation for the observed light variations is the formation of an inner disk region, indicating ongoing dust replenishment Wahhaj et al. (2003), which can have an origin in the planetesimal belts Okamoto et al. (2004). Warm debris disks around evolved stars cannot be primary remnants of planet formation, but formed quite recently (Wyatt et al. 2007), most likely in a late heavy bombardment-like scenario, or from a very massive reservoir of planetesimals that survived planet formation (such as the Oort Cloud in the Solar System). In the presence of warm dust, excess radiation is expected in mid-infrared. However, no excess was detected at 24 µm by the Multiband Imaging Photometer for Spitzer (MIPS) for DE Boo (Gáspár et al. 2013). This may be because the detection limit of MIPS for debris disks at a distance of a few AU around late-type stars is above $f \sim 10^{-5}$. However, variations in the light curve of DE Boo imply a lower fractional luminosity, which is therefore below what may be detectable with MIPS. This suggests that the presence of a faint warm debris disk cannot be excluded.

The current observations covered 2×2 days of DE Boo, and
show that DE Boo is an important example of a system with an
edge-on debris disk that we are looking at through the debris disk
itself. Further observations with TESS in 2022, covering 28 days,
will be an ideal opportunity to confirm our interpretations.

Acknowledgements. CHEOPS is an 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. We acknowl-
edge the support of the PRODEX Experiment Agreement No 4000137122
between the ELTE Eötvös Loránd University and the European Space Agency
(ESA-D/SCI-LE-2021-0025). Gy.MSz. 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. L.K. is supported by the Hungarian
National Research, Development and Innovation Office grants PD-134784 and
K-131508. L.K. is a Bolyai János Research Fellow. ABr was supported by the
SNSA. F.K. acknowledges funding from the Centre National d'Etudes Spatiales,
from the French National Research Agency (ANR) under contract number
ANR-18-CE31-0019 (SPlaSH), and from the Université Paris Sciences et Lettres
under the DIM-ACAV program Origines et conditions d'apparition de la vie.
D.G. gratefully acknowledges financial support from the CRT foundation under
Grant No. 2018.2323 "Gaseous or rocky? Unveiling the nature of small worlds".
L.M.S. gratefully acknowledges financial support from the CRT foundation
under Grant No. 2018.2323 'Gaseous or rocky? Unveiling the nature of small
worlds'. 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. acknowledge support from FCT through FCT contract nr.
CEECIND/00826/2018 and POPH/FSE (EC). Y.A. and M.J.H. acknowledge the
support of the Swiss National Fund under grant 200020_172746. We acknowl-
edge 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 no. 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. A.C.C. acknowledges support from STFC
consolidated grant numbers ST/R000824/1 and ST/V000861/1, and UKSA
grant number ST/R003203/1. 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. and A.D. acknowledge
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. SH gratefully acknowledges CNES funding through
the grant 837319. KGI 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
gratefully acknowledged to the HPC resources of MesoPSL financed by the Region Ile
de France and the project Equip@Meso (reference ANR-10-EQXP-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.Bo., G.Br., V.Na., I.Pa.,
G.Pi., R.Ra., G.Sc., V.Si., and T.Zi. acknowledge support from CHEOPS
ASI-INAF 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.
V.V.G. is an F.R.S-FNRS Research Associate. NAW acknowledges UKSA
grant ST/R004838/1.

**References**

Allard, F., Homeier, D., & Freytag, B. 2011, in 16th Cambridge Workshop on
Cool Stars, Stellar Systems, and the Sun, eds. C. Johns-Krull, M. K. Brown-
ing, & A. A. West, Astronomical Society of the Pacific Conference Series,
448, 91

Aumann, H. H., Gillett, F. C., Beichman, C. A., et al. 1984, ApJ, 278, L23

Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42

Benz, W., Broeg, C., Fortier, A., et al. 2021, Exp. Astron., 51, 109

Berdyugina, S. V. 2005, Living Rev. Solar Phys., 2, 8

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

Budding, E. 1977, Ap&SS, 48, 207

Castelli, F., & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, eds.
N. Piskunov, W. W. Weiss, & D. F. Gray, IAU Symposium, 210, A20

Deline, A., Queloz, D., Chazelas, B., et al. 2020, A&A, 635, A22

Delrez, L., Ehrenreich, D., Alibert, Y., et al. 2021, Nat. Astron., 5, 775

Eker, Z., Ak, N. F., Bilir, S., et al. 2008, MNRAS, 389, 1722

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

Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1

Gaidos, E. 2022, Res. Notes AAS, 6, 49

Gaidos, E., Jacobs, T., LaCourse, D., et al. 2019, MNRAS, 488

Gáspár, A., Rieke, G. H., & Balog, Z. 2013, ApJ, 768, 25

Gray, R. O., Corbally, C. J., Garrison, R. F., McFadden, M. T., & Robinson, P. E.
2003, AJ, 126, 2048

Gray, R. O., Saken, J. M., Corbally, C. J., et al. 2015, AJ, 150, 203

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

Henry, G. W., Fekel, F. C., & Hall, D. S. 1995, AJ, 110, 2926

Holman, M. J., & Wiegert, P. A. 1999, AJ, 117, 621

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

Jancart, S., Jorissen, A., Babusiaux, C., & Pourbaix, D. 2005, A&A, 442, 365

Kalas, P., Fitzgerald, M. P., & Graham, J. R. 2007, ApJ, 661, L85

Kennedy, G. M. 2014, MNRAS, 447, L75

Kennedy, G. M., Murphy, S. J., Lisse, C. M., et al. 2014, MNRAS, 438, 3299

Kennedy, G. M., Kenworthy, M. A., Pepper, J., et al. 2017, Roy. Soc. Open Sci.,
4, 160652

Kriskovics, L., Kővári, Z., Vida, K., et al. 2019, A&A, 627, A52

Kupca, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C., & Weiss, W. W.
1999, A&AS, 138, 119

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

Leitzinger, M., Odert, P., Greimel, R., et al. 2020, MNRAS, 493, 4570

Lindgren, L., Bastian, U., Biermann, M., et al. 2021, A&A, 649, A4

Lisse, C. M., Beichman, C. A., Bryden, G., & Wyatt, M. C. 2007, AJ, 658, 584

Lisse, C. M., Chen, C. H., Wyatt, M. C., & Morlok, A. 2008, ApJ, 673, 1106

Marshall, J. P., Kirchschlager, F., Ertel, S., et al. 2014, A&A, 570, A114

Maxted, P. F. L., Ehrenreich, D., Wilson, T. G., et al. 2021, MNRAS, 514, 77

Moór, A., Ábrahám, P., Derekas, A., et al. 2006, ApJ, 644, 525

Morris, B. M., Delrez, L., Brandeker, A., et al. 2021a, A&A, 653, A173

Morris, B. M., Heng, K., Brandeker, A., Swan, A., & Lendl, M. 2021b, A&A,
651, A12

Namekata, K., Maehara, H., Notsu, Y., et al. 2019, ApJ, 871, 187

Okamoto, Y. K., Kataza, H., Honda, M., et al. 2004, Nature, 431, 660

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

Piskunov, N., & Valenti, J. A. 2017, A&A, 597, A16

Ribárik, G., Oláh, K., & Strassmeier, K. G. 2003, Astron. Nachr., 324, 202

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

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Smith, B. A., & Terrile, R. J. 1984, Science, 226, 1421

Solanki, S. K. 2003, A&ARv, 11, 153

Stark, C. C., Kuchner, M. J., Traub, W. A., et al. 2009, ApJ, 703, 1188

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

Trilling, D. E., Stansberry, J. A., Stapelfeldt, K. R., et al. 2007, ApJ, 658, 1289

Tajiri, T., Kawahara, H., Aizawa, M., et al. 2020, ApJS, 251, 18

Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141

Vida, K., Bódi, A., Szklenár, T., & Seli, B. 2021, A&A, 652, A107

Wahhaj, Z., Koerner, D. W., Ressler, M. E., et al. 2003, ApJ, 584, L27

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Wyatt, M. C., Smith, R., Greaves, J. S., et al. 2007, ApJ, 658, 569

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

<footnote n="1">¹ Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, 1121 Budapest, Konkoly Thege Miklós út 15–17, Hungary e-mail: boldog.adam@csfk.org</footnote>
² CSFK, MTA Centre of Excellence, Budapest, Konkoly Thege Miklós út 15–17, 1121, Hungary

³ MTA-ELTE Exoplanet Research Group, 9700 Szombathely, Szent Imre h. u. 112, Hungary

⁴ ELTE Eötvös Loránd University, Gothard Astrophysical Observatory, 9700 Szombathely, Szent Imre h. u. 112, Hungary

⁵ Department of Astronomy, Stockholm University, AlbaNova University Center, 10691 Stockholm, Sweden

⁶ Sorbonne Université, CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France

⁷ LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, 5 place Jules Janssen, 92195 Meudon, France

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

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

¹⁰ Division Technique INSU, CS20330, 83507 La Seyne sur Mer Cedex, France

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

¹² Dipartimento di Fisica, Universita degli Studi di Torino, via Pietro Giuria 1, 10125, Torino, Italy

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

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

¹⁵ Institut d'astrophysique de Paris, UMR7095 CNRS, Université Pierre & Marie Curie, 98bis blvd. Arago, 75014 Paris, 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

²⁰ Physikalisches Institut, University of Bern, Bern, Switzerland

²¹ 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, 38000 Grenoble, France

²⁷ Center for Space and Habitability, Gesellschaftsstrasse 6, 3012 Bern, Switzerland

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

²⁹ 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

³¹ ESTEC, European Space Agency, 2201AZ Noordwijk, The Netherlands

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

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

³⁴ Space sciences, Technologies and Astrophysics Research (STAR) Institute, Université de Liège, Allée du 6 Août 19C, 4000 Liège, Belgium

³⁵ INAF, Osservatorio Astrofisico di Catania, Via S. Sofia 78, 95123 Catania, Italy

³⁶ 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

⁴¹ 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, Keele, Staffordshire ST5 5BG, UK

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

⁴⁶ 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

⁴⁹ 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

⁵¹ Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, UK
**Appendix A: Periodogram of**
**spot-model-subtracted data**

We applied the Lomb-Scargle periodogram to the spot-model-
subtracted data from both visits before applying the fourth-order
polynomial to check for remaining periods. Figure A.1 shows
the resulting periodogram, where the strongest peak was seen
at 1.79 days. Another peak was present at around 100 min due
to the orbital period of CHEOPS. We concluded that no higher
harmonics of the stellar rotation of 10.39 days were present in
the spot-model-subtracted data.

<figure><img src="image_5.png" /><figcaption>**Fig. A.1.** Lomb-Scargle periodogram of the spot-model-subtracted data. The peak with the maximum power belongs to the period of 1.79 days. We note that the effect of the 98.77 minute rotational period of CHEOPS also appears as a peak around 100 minute.</figcaption></figure>