# Non-transiting exoplanets as a means of understanding star-planet interactions in close-in systems

C. Gourvès, S. N. Breton, A. Dyrek, A. F. Lanza, R. A. García, S. Mathur, À.

R. G. Santos, A. Strugarek

► To cite this version:

C. Gourvès, S. N. Breton, A. Dyrek, A. F. Lanza, R. A. García, et al. Non-transiting exoplanets as a means of understanding star-planet interactions in close-in systems. Astronomy & Astrophysics - A&A, 2025, 699, pp.A226. 10.1051/0004-6361/202553791. cea-05159459

HAL Id: cea-05159459

https://cea.hal.science/cea-05159459v1

Submitted on 12 Jul 2025

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
# Non-transiting exoplanets as a means of understanding star-planet interactions in close-in systems

C. Gourvès<fnref n="1" />,<fnref n="2" />,*, S. N. Breton<fnref n="2" />, A. Dyrek<fnref n="3" />,<fnref n="1" />, A. F. Lanza<fnref n="2" />,
R. A. García<fnref n="1" />, S. Mathur<fnref n="4" />,<fnref n="5" />, Â. R. G. Santos<fnref n="6" />, and A. Strugarek<fnref n="1" />

Received 17 January 2025 / Accepted 17 May 2025

## ABSTRACT

Previous studies showed evidence of a dearth of close-in exoplanets around fast rotators, which can be explained by the combined action of intense tidal and magnetic interactions between planets and their host star. Detecting more exoplanets experiencing such interactions, with orbits evolving on short timescales, is therefore crucial to improve our understanding of the underlying physical mechanisms. For this purpose, we performed a new search for close-in non-transiting substellar companions in the *Kepler* data, focusing on orbital periods below 2.3 days. We focused on main-sequence solar-type stars and subgiant stars for which a surface rotation period was measured. For each star, we looked for an excess in the power spectral density of the light curve, which could correspond to the signature of a close-in non-transiting companion. We compared our candidates with existing catalogues to eliminate potential contaminants in our sample, and we visually inspected the phase-folded light curve and its wavelet decomposition. We identify 88 stars exhibiting a signature consistent with the presence of a close non-transiting substellar companion. We show that the objects in our sample are located mostly within the dearth zone, emphasising the importance of performing follow-up of such systems in order to gather observational evidence of star-planet interactions.

**Key words.** methods: data analysis – planets and satellites: detection – planet–star interactions – stars: low-mass – stars: solar-type

## 1. Introduction

Since the first observation of a planet orbiting a solar-type star (Mayor & Queloz 1995), more than 5500 exoplanets have been confirmed, the majority detected by the *Kepler* satellite with the transit method (Borucki et al. 2010). In particular, this has opened the possibility of studying star-planet interactions shaping orbital architecture both at the scale of individual star-planet systems (e.g. Patra et al. 2017; Barker et al. 2024) and for larger samples (e.g. McQuillan et al. 2013; Messias et al. 2022; García et al. 2023).

Observing close-in planets, which exhibit orbital periods shorter than 1 day and experience strong tidal (e.g. Mathis 2018) and magnetic (e.g. Strugarek 2018) interactions with their host star, is crucial to calibrate dissipation processes related to these phenomena. Unfortunately, ultra-close-in planets, depending on the planetary mass, its orbital period, and the rotational period of the star, may be pushed away or experience orbital decay and disruption on very short timescales (e.g. Ahuir et al. 2021). This means that we are more likely to observe planets that are less dynamically perturbed and located further away from their host. Among all the planets detected using the transit method, approximately 97% have orbital periods exceeding one day<fnref n="1" />.

This highlights the need for alternative detection methods to study ultra close-in planets. Statistically, non-transiting systems are likely far more numerous than transiting ones, suggesting that observing them could increase our chances of detecting planets very close to their star.

Having rather large modulations in their orbital phase curve, close-in planets can be detected with photometric observations, even if they are not transiting. Until now, only a few non-transiting companions have been found (Millholland & Laughlin 2017; Lillo-Box et al. 2021; Cullen & Bayliss 2024). These modulations are attributed to three main physical processes (Shporer 2017). The first is atmospheric, encompassing stellar light reflected by the companion's atmosphere and thermal emission, due to heating of the companion's atmosphere. The second is the tidal ellipsoidal distortion, induced by the tidal forces exerted by the companion on its host star. The third process is the Doppler boosting effect, caused by variations in the radial velocity between the object and the observer, leading to a modulation in the observed flux of the star. Fig. 1 shows an example of the three modulation components for a non-transiting planet, as well as the composite phase-curve, generated using Eq. (4), (7), (9), and (10) from Shporer (2017) and a black body emission law computed at the planetary equilibrium temperature, as a function of the star-planet-observer phase angle (Seager 2010). A realistic set of stellar and planetary parameters was chosen, analogous to the ones expected for a simplified bookcase-like system. A hot Jupiter with a radius of 1 $R_{jup}$, a mass of 1 $M_{jup}$,

<footnote n="1">¹ Université Paris-Saclay, Université Paris Cité, CEA, CNRS, AIM, 91191 Gif-sur-Yvette, France</footnote>

<footnote n="2">² INAF – Osservatorio Astrofisico di Catania, Via S. Sofia, 78, 95123 Catania, Italy</footnote>

<footnote n="3">³ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA</footnote>

<footnote n="4">⁴ Instituto de Astrofísica de Canarias, La Laguna, Tenerife, Spain</footnote>

<footnote n="5">⁵ Departamento de Astrofísica, Universidad de La Laguna, La Laguna, Tenerife, Spain</footnote>

<footnote n="6">⁶ Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal</footnote>

<footnote>* Corresponding author: clemence.gourves@cea.fr</footnote>

<footnote n="1">¹ According to the NASA Exoplanet Archive, as of April 6, 2025. The database can be accessed at https://exoplanetarchive.ipac.caltech.edu/index.html</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Simulated contributions of the three modulations to the phase-curve amplitude according to the phase angle. The modulation linked to atmospheric processes is shown in black and includes reflection and emission modulations. The ellipsoidal distortion is shown in purple, while the Doppler boosting is in orange. The overall contribution is shown by the composite phase-curve in yellow. Here, an eclipse would occur at the inferior conjunction, while the transit would happen at the superior conjunction with the observer.</figcaption></figure>

and an orbital period of 2 days was considered. The eccentricity was set to zero, while the inclination of the orbit on the line of sight was 60°, the maximum value taken for non-transiting systems (Winn 2014). We opted for a Bond albedo of 0.5 (e.g. Seager 2010), which is typical for a hot Jupiter-type planet, and picked an atmospheric recirculation rate of 50% (Fortney et al. 2021). We chose 1 $M_⊕$, 1 $R_⊕$, and $T_{\text{eff}}$ = 6000 K as stellar parameters. In this configuration, the atmospheric processes outweigh the other two phenomena.

While the atmospheric processes are mostly driven by stellar irradiation, Doppler boosting and tidal ellipsoidal distortion are induced by the bodies' gravitational interaction, and hence they depend on the stellar and planet masses. The greater the companion's orbital velocity and mass, the stronger the Doppler boosting. However, depending on the mass ratio between the star and the planet, for short orbital periods, the ellipsoidal distortion effect might outweigh the other two phenomena in amplitude.

The layout of the paper is as follows. In Sect. 2, we present the sample of *Kepler* light curves we consider, and we explain how we look for a non-transiting companion signature. In Sect. 3, we perform a fit of the photometric modulations to extract additional constraints on the candidates' global parameters. In Sect. 4, we discuss the perspectives that the characterisation of this sample may open.

## 2. Observational data and search methodology

In mono-band photometry and when blindly searched for in a large sample, non-transiting companion signatures may easily be confused with modulations arising from stable active longitudes of fast-rotating young stars. Therefore, prior knowledge of the rotation of the stellar host is necessary to undertake our analysis. To avoid confusing exoplanetary signals with the rotation period (e.g. García et al. 2014), we considered the catalogue by Santos et al. (2019, 2021), which consists of 55 232 main-sequence and subgiant FGKM stars.

To illustrate the procedure we describe below, we selected KIC 5697777, an F-type star with a rotation period of 3.90 days. The results obtained for this star are presented in Fig. 2 and are discussed in this section.

### 2.1. Unravelling signatures of non-transiting objects

We used *Kepler* light curves calibrated with the KEPSEISMIC method<fnref n="2" /> (García et al. 2011; Pires et al. 2015), whereby signals with a period longer than 20 days were filtered out (see Santos et al. 2019) in order to be less affected by the instrumental noise and because we are interested in close-in objects. The light curve of KIC 5697777 is presented in panel a of Fig. 2.

The power spectral density (PSD) of a stellar light curve (Fig. 2, panel b) allows for the decomposition of the signal into its various frequency components. Identifying a high-amplitude peak, typically from a few tens to a few thousands of parts per million, is then akin to detecting a peak that corresponds to the passage of a non-transiting object. As is seen in Fig. 1, if tidal interactions dominate the flux modulation, the second harmonic of the orbital signal is expected to be the largest peak of the harmonic pattern (when still excluding the stellar rotation period signal). Nevertheless, it is the frequency of the first harmonic that corresponds to the orbital period of the non-transiting companion object.

To identify a non-transiting object signal, we computed the PSD of each of the 55 232 stars' light curves. To remove false positives due to low-inclination stellar binary systems, we focused on systems with a signal at a period distinct from the rotation period of the targeted star. In a close binary stellar system, tidal effects typically tend to synchronise the rotational periods of the two stars (e.g. Zahn 1975). In contrast, a planetary-mass companion generally lacks sufficient mass to induce such synchronisation. The aim here is to identify significant high-amplitude peaks in the PSD, different from the peak associated with the stellar rotation period. Given that only close-in objects should be detectable with our method, we focused on objects with short orbital periods and we searched for peaks with a frequency higher than 5 µHz (2.3 days). Below 5 µHz, the main signal in the PSD of a main-sequence solar-like star usually comes from the surface rotational modulation, granulation, and activity. In addition, because the object would be too distant from the star, the phase curve modulation amplitude would be too low to be detectable in visible photometry (e.g. Shporer 2017). A high-amplitude peak with a frequency higher than 5 µHz can then be associated with the presence of a non-transiting object. If we know the stellar rotation period and its harmonics, we do not expect any additional excess power of this nature associated with the dynamics of the star alone, when considering a main-sequence solar-type star (e.g. García & Ballot 2019).

A threshold limit was applied in the PSDs, to distinguish sharp high-amplitude peaks associated with non-transiting companions from background noise. The statistical distribution of the PSD function is a $\chi^2$ distribution with two degrees of freedom (Woodard 1984). Following Appourchaux et al. (2000), we calculated the power level at which a peak due to noise has a probability, $p_{\text{det}}$, of appearing over a frequency range, $\Delta_{\text{det}}$. Since the number of considered bins in the bandwidth will increase with the length of the time series, this threshold depends on the observation time, $T$, which is about 4 years for the *Kepler* observations. The ratio of the statistical $\chi^2$ threshold

<footnote n="2">² The light curves can be downloaded at the following address: https://archive.stsci.edu/prepds/kepseismic/</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Flux analysis of KIC 5697777. (a) Light curve calibrated according to the KEPSEISMIC method and filtered at 20 days. (b) Power spectral density of the light curve. The rotation period of the star is taken from Santos et al. (2021) and is located within the orange area. The green area shows the high-amplitude signal peak at around 0.5 days. The second high-amplitude signal peak located at around 0.25 days is the second harmonic of the orbital signal. (c) Phase-folding of the light curve with an additional filter to remove the stellar rotational contribution. The red curve is a smoothing of the black curve, obtained with a triangular filter. (d) Wavelet power spectrum of the light curve. A high-amplitude periodicity is highlighted by orange and red, whereas yellow and blue denote a low-power amplitude. The yellow bands with a constant period throughout time correspond to the high-amplitude peak and its second harmonic, identified in panel b. The reddish band corresponds to the rotation period and its harmonics, highlighted in orange in panel b. The shaded area represents the cone of influence of the reliable periodicities.</figcaption></figure>

to the average signal amplitude, $\bar{s}$, is given by $s_{det}/\bar{s} \approx \ln(T) + \ln(\Delta_{det}) - \ln(p_{det})$, with $p_{det}$ the detection probability threshold, $T$ in µs and $\Delta_{det}$ in microhertz. Here, $\Delta_{det}$ is within the range of 5 to 284 µHz, the Nyquist frequency of *Kepler* observations. After some empirical tests, we decided to consider a $s_{det}/\bar{s}$ value of approximately 260, which corresponds to an extremely low detection probability threshold, $p_{det}$, of $10^{-100}$ (this corresponds to a sigma of 259). This choice allowed us to filter out most of the signal related to the harmonics of stellar rotation, while preserving the high-amplitude peaks we were looking for.

We computed $s_{det}/\bar{s}$ by dividing the PSD by a factor of $\alpha = (9/8)^3 \times \text{median (PSD)}$, corresponding to the corrected mean of the PSD. Nevertheless, the background signal of the PSD is non-uniform, predominantly influenced by the convective signal and the magnetic activity of the star, which results in higher amplitudes at low frequencies compared to high frequencies (García & Ballot 2019). Hence, adopting such a low probability threshold, $p_{det}$, reduces the likelihood of obtaining false positives and improves the accuracy of our signal analysis. With this method, we preselected 4788 stars.

Subsequently, we conducted a visual inspection of the PSDs displaying prominent peaks distinct from the stellar rotation period. This is illustrated in panel b of Fig. 2, in which the PSD exhibits a high-amplitude peak around 0.52 days (green area), while the stellar rotation period of KIC 5697777 corresponds to a period of 3.9 days (orange area). For our remaining stars, the PSDs presenting peaks that are part of the harmonic chain of the rotation period are excluded from the catalogue. However, we decided to keep in our sample targets for which the detected peak lies within the harmonic pattern of surface rotation but exhibits a clearly distinct morphology with respect to it; that is, for which we observe a sharp peak amidst the wide rotation harmonics. This analysis step reduced the considered sample to 362 stars.

## 2.2. Comparison with existing catalogues

To eliminate potential contaminants in our sample, we conducted a thorough cross-check with existing catalogues. The aim was to remove any stars previously identified as hosting binaries, brown dwarfs, transiting exoplanets, or identified false positives from the *Kepler* Object of Interest catalogue (KOI, Akeson et al. 2013). The first cross-analysis was performed using data from the National Aeronautics and Space Administration (NASA) Exoplanet Archive<fnref n="3" />, resulting in the removal of four stars confirmed to harbour transiting planets with an orbital period matching the ones detected in our signals (see Appendix A.1 for an example). We also removed 59 stars identified as false positives of transiting planets. A second comparison was carried out with the catalogue from Berger et al. (2018), allowing for the removal of 15 stars associated with binary systems. The detected signal in the PSD is indeed most likely attributable to the binary star system. Lastly, we compared our results to the Gaia DR3 non-single stars catalogue (Pourbaix et al. 2022), which allowed us to exclude another star from our sample, as it is part of a binary system. A total of 79 stars were removed from our sample following these cross-checks and upon completion of this step, our refined sample consists of 283 stars possibly hosting non-transiting exoplanets.

<footnote n="3">³ See Footnote 1 for the url.</footnote>
## 2.3. Phase-folding, wavelet analysis, and visual inspection

To further refine our study and pinpoint potential indicators of yet undetected transits, false positives, or binary systems, we phase-folded the light curves of the 283 candidates on the orbital period that we obtained from the threshold search conducted in Sect. 2.1. To remove the stellar contribution at lower frequencies, we applied a finite impulse response filter<fnref n="4" /> before phase-folding the signal. The phase-folding analysis for KIC 5697777 is presented in panel c of Fig. 2. The complex morphology of the phase curve suggests that different effects (see Fig. 1) combine to create a possible non-transiting modulation. A subsequent visual inspection was performed to exclude any remaining stars for which the phase-folding of the star signal fails to match the expected pattern of a non-transiting exoplanet passage (as in Fig. 2, panel c). As is shown in Appendix A.2, if the trend of the phase-folded light curve differs from a stable and coherent signal over time, the candidate was discarded. This step allowed for the refinement of our overall sample to 245 candidates.

We also computed the wavelet power spectrum (WPS, Torrence & Compo 1998; Liu et al. 2007) of the light curve, and we only kept stars that present stable modulations over time, as is shown in panel d of Fig. 2 for KIC 5697777. The aim is to determine whether the signal attributed to the non-transiting object exhibits a stable period throughout the different quarters of Kepler's observations. If the signal disappears during at least one quarter, this suggests that the signal may result from contamination by a nearby star. An example of a wavelet power spectrum that does not match a non-transiting signal is provided in Appendix A.3. This step removed 92 stars from our sample.

To ensure that the Kepler signal of every star studied in this work was not contaminated by the light of a stellar neighbour with comparable brightness, we performed a search for potential contaminants using the Gaia catalogue (Gaia Collaboration 2023). We looked for Gaia objects within 40" of each target that have an apparent magnitude in the Gaia G band smaller than $G_{candidate} + 1$, where $G_{candidate}$ is the apparent magnitude of our Kepler candidate in the Gaia G band. This results in the identification of 70 Gaia objects near 38 of our targets (see Table A.1 in Appendix A.4). For 11 of these 38 targets, the potential contaminants are located within a few arcseconds, which could result in the contamination of the Kepler light curve. Consequently, we chose to remove these 11 stars from our sample. For the remaining 27 targets, the potential contaminants are at least 10" away, making significant contamination issues unlikely. Similarly to what was done in Noraz et al. (2022), we also performed a visual inspection on the 142 stars remaining at this stage. We individually inspected the targets using the Two Micron All Sky Survey (2MASS) images provided on the Kepler Asteroseismic Science Operations Center (KASOC) database<fnref n="5" /> (Kjeldsen et al. 2010). If the star is within the halo of a brighter star, it was excluded from our sample because the short-period signal might come not from a non-transiting object but rather from the nearby star in question. An example of a star whose signal is likely contaminated by a neighbouring brighter star is shown in Appendix A.4. After this final step of selection, we obtained a list of 88 non-transiting exoplanet host candidates. At this point, it is important

to underline that none of our candidates are in the list of 60 targets presented by Millholland & Laughlin (2017). We recall that they used an automated algorithm in the search for non-transiting short-period giant planets in Kepler, looking for targets with an orbital period between 1 and 5 days. While they modelled the non-transiting signal and stellar variability together, we used our prior knowledge of stellar rotation to filter out the stellar contribution in the light curve. More importantly, we extended the search towards shorter orbital periods.

## 3. Fit of the photometric data

As a final step, we attempted to extract additional information on the global parameters of our candidates by fitting their signature in the light curve. Two stars were excluded from this analysis, as their signatures could not be reliably fitted in the light curve. The signal of KIC 7379583 seems contaminated by residual instrumental drifts, leading to peaks in the phase-folded light curve that prevent the fit from converging. For KIC 11350907, the signal amplitude is 84 ppm, while the random noise is around 292 ppm, making it difficult for the fit to converge.

When only photometric data were available, the resulting signal presented in Fig. 1 depends on several free parameters, leading to degeneracies between different effects. To address this, we distinguish two cases. On the one hand, when the phase-folding of the light curve exhibits a modulation that is approximately sinusoidal, we assumed that the non-transiting signal is dominated by emission and reflection. We then fitted only this component in the data, using the following equation (see Shporer 2017; Lillo-Box et al. 2021):

$$F(\phi) = A_0 - A_{\text{atm}} \cos(2\pi[\phi + \delta]), \quad (1)$$

where $A_0$ is a flux offset parameter, and the phase, $\phi$, is given by

$$\phi = (t - T_0)/P_{\text{orb}}, \quad (2)$$

where $T_0$ is taken to correspond to the time of inferior conjunction, and $\delta$ is a phase offset with respect to this conjunction.

The amplitude, $A_{\text{atm}}$, is given by

$$A_{\text{atm}} = 57 \alpha_{\text{refl}} \sin i \left(\frac{M_{\star}}{M_{\odot}}\right)^{-2/3} \left(\frac{P}{\text{day}}\right)^{-4/3} \left(\frac{R_{\text{p}}}{R_{\text{J}}}\right)^{2} \text{ ppm}, \quad (3)$$

where $\alpha_{\text{refl}}$ depends on the albedo and heat recirculation efficiency. We adopted this approach for the 77 targets listed in Table B.1.

On the other hand, when the modulation in the phase-folding was more complex (see panel c of Fig. 2), we attempted to perform the fit considering the three components of the signal,

$$F(\phi) = A_0 - A_{\text{atm}} \cos(2\pi[\phi + \delta]) - A_{\text{ellip}} \cos(4\pi\phi) + A_{\text{boost}} \sin(2\pi\phi), \quad (4)$$

where $A_{\text{ellip}}$ is the amplitude of the tidal deformation,

$$A_{\text{ellip}} = 13\alpha_{\text{ellip}} \sin i \left(\frac{R_{\star}}{R_{\odot}}\right)^3 \times \left(\frac{M_{\star}}{M_{\odot}}\right)^{-2} \left(\frac{P}{\text{day}}\right)^{-2} \left(\frac{M_{\text{p}} \sin i}{M_{\text{J}}}\right) \text{ ppm}, \quad (5)$$

<footnote n="4">⁴ The filter is implemented through a combination of scipy.signal.firls and scipy.signal.filtfit. See pre-process function of the star-privateer module (Breton et al. 2024): https://star-privateer.readthedocs.io/en/latest/usage/api/analysis_pipeine.html#preprocessing</footnote>

<footnote n="5">⁵ The KASOC database can be accessed at https://kasoc.phys.au.dk.</footnote>
with $\alpha_{ellip}$ a coefficient that depends on the limb and gravity darkening of the host star (Esteves et al. 2013). $A_{boost}$ is the amplitude of the Doppler boosting contribution,

$$A_{boost} = 2.7 \alpha_{boost} \left(\frac{P}{\text{day}}\right)^{-1/3} \left(\frac{M_{\star}}{M_{\odot}}\right)^{-2/3} \left(\frac{M_{\text{p}} \sin i}{M_{\text{J}}}\right) \text{ppm}, \quad (6)$$

where $\alpha_{boost}$ is the photon-weighted bandpass-integrated boosting factor. After Millholland & Laughlin (2017), we used a parameterisation in terms of the effective temperature, $T_{\text{eff}}$, for $\alpha_{ellip}$ and $\alpha_{boost}$,

$$\begin{align}
\alpha_{\text{ellip}} &= -(2.2 \times 10^{-4} \text{ K}^{-1}) T_{\text{eff}} + 2.6, \nonumber \\
\alpha_{\text{boost}} &= -(6 \times 10^{-4} \text{ K}^{-1}) T_{\text{eff}} + 7.2. \tag{7}
\end{align}$$

We adopted this approach for nine targets presented in Table B.2 in Appendix B.

In practice, Eq. (1) and (4) were fitted to the data using Bayesian inference. The log-likelihood function is

$$\ln \mathcal{L}_{\phi}(S_{\phi}, \theta) = -\frac{1}{2} \times \sum_{i=1}^{n} \frac{(S_{\phi}(t_i) - \mu_{\phi}(\phi_i, \theta))^2}{\sigma_{\phi}(t_i)^2}, \quad (8)$$

where $S_{\phi}(t)$ is the light curve data and $\sigma_{\phi}(t)$ is the normally distributed noise, considering independent errors. $\mu_{\phi}(\phi, \theta)$ is the model depending on the parameter, $\theta$, as a function of the orbital phase, $\phi$, as expressed in Eq. (1) when following the first approach, and Eq. (4) when following the second approach. We used the UltraNest sampler described in Buchner (2021), which is based on a nested sampling algorithm. To avoid degeneracies between $T_0$ and $\delta$, we set $\delta$ to 0 and fitted for $T_0$, $A_{\text{atm}}$, and $A_0$ for the first approach, and $T_0$, $A_{\text{atm}}$, $A_0$, $A_{\text{boost}}$, and $A_{\text{ellip}}$ for the second approach. We fixed the period, $P_{\text{orb}}$, at the value inferred from the PSD of the light curve. The priors are $\mathcal{U}(-10, 10)$ for $A_0$ and $\mathcal{U}(-10^5, 10^5)$ for $A_{\text{atm}}$, $A_{\text{ellip}}$, and $A_{\text{boost}}$, where $\mathcal{U}(X, Y)$ denotes a uniform probability distribution between the bounds X and Y. We also considered uniform priors for $T_0$ and inferred them with the timing corresponding to the minimum amplitude of the phase-folded light curves. The posterior distributions of two datasets, chosen as examples of each approach, are shown in Fig. B.1 and B.2 of Appendix B.

At this stage, it is necessary to emphasise that the model used to perform the fit relies on a number of simplifying assumptions: the eccentricity is assumed null, the atmospheric contribution from the planet is modelled in the simplest possible way, and the stellar tidal response is considered to arise solely from equilibrium tides, with any potential contribution from dynamical tides being neglected. We therefore advise the reader to keep all these caveats in mind when they consider the results of the fit we present below. Some of the limitations of this approach are discussed in Shporer (2017).

An example of the fit result when the phase-folding of the light curve exhibits a modulation that is approximately sinusoidal is illustrated by Fig. 3. It shows the normalised light curve as a function of orbital phase for KIC 4930913, showing the observational data in grey, along with the best-fit model in orange. The model closely follows the general trend of the data despite significant scatter. The bottom panel displays the residuals, defined as the difference between the data and the model. The residuals are centred around zero and show no evidence of large-scale systematic trends, indicating a satisfactory fit.

An example of the fit result, based on the three components of the signal described by Eq. (4), is illustrated in Fig. 4 for KIC

<figure><img src="image_4.png" /><figcaption>Fig. 3. Phase-folded light curve of KIC 4930913, together with the best-fit model in orange. The red points correspond to data binned in 1-day intervals. The residuals show the difference between the model and the data points and are free of large systematic trends.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Fig. 4. Phase-folded light curve of KIC 4264634, together with the best-fit model in orange. The red points correspond to data binned in 1-day intervals. The scatter of the residuals is higher due to a larger dispersion of the data points but remains free of large systematic trends as well.</figcaption></figure>

4264634. The phase-folding of the light curve of this star shows a more complex pattern (in grey). As is seen in Fig. 3, the best-fit model (in orange) follows the general trend of the data and the residuals are again centred around zero.

The full results of the fits are displayed in Tables B.1 and B.2. For all 86 datasets, the posterior distribution of each parameter as well as the residuals of the fit follow a normal distribution. This indicates that the sampling successfully converged.

Using Eq. (3) and the fitted $A_{\text{atm}}$ parameter, we derived the $\alpha_{\text{ref}} \sin i (R_{\text{p}}/R_{\text{J}})^2$ parameter, obtaining a minimum value of 0.08 and a maximum value of 21.01 (see Table B.1). Provided that both $\alpha_{\text{ref}}$ and $i$ are known, one can then infer the radius of
the object considered. Through Eq. (5) and (6), and the fitted parameters, $A_{\text{ellip}}$ and $A_{\text{boost}}$, we retrieved the $(M_p \sin^2 i/M_J)$ and $(M_p \sin i/M_J)$ parameters, respectively. For the 77 targets analysed in the first approach, the orbital inclination is unknown and cannot be derived from the fits. However, an upper limit can be estimated by considering the geometric condition required for a transit to occur. For a circular orbit, a transit event occurs for $i \ge i_{\text{max}}$ with $i \in [0; \pi/2]$.

$$i_{\text{max}} = \cos^{-1}\left(\frac{R_{\star}}{a}\right), \qquad (9)$$

where a is the semi-major axis of the orbit. Assuming that all of our non-transiting companions follow a circular orbit (Mazeh 2008), they would then have an inclination inferior to $i_{\text{max}}$. To ensure that the inclination value is in compliance with the fit, we inferred an inclination value using the ratio of $(M_p \sin^2 i/M_J)$ and $(M_p \sin i/M_J)$, for the nine targets where $A_{\text{ellip}}$ and $A_{\text{boost}}$ are accessible. For three hosts (KIC 11668095, KIC 5440365, and KIC 8844761), the uncertainties of the parameters $(M_p \sin^2 i/M_J)$ and $(M_p \sin i/M_J)$ are too large to infer $\sin i$ values. Regarding KIC 4373708, KIC 11702835, and KIC 5622796, the $\sin i$ values are inconsistent, suggesting that these stars are probably out of the range of validity of the model and would require more complex models (e.g. Penoyre & Sandford 2019). Nevertheless, we were able to infer the inclination angle for three systems. For KIC 11550689, we obtained an inclination of $6^\circ$, for KIC 5697777, $0.4^\circ$, and for KIC 4264634, $38^\circ$. All of these values are lower than the maximum inclination angles for these systems (see Table B.2). The corresponding masses that we can therefore derive in these three cases are 32, 1412, and 21 M$_{\text{J}}$, respectively. In the cases of KIC 11550689 and KIC 4264634, if the numerous assumptions made in this modelling approach are correct, this would suggest that the non-transiting objects might be brown dwarfs. In the case of KIC 5697777, the inferred mass is approximately 1.3 M$_{\odot}$, which is obviously inconsistent with the hypothesis of a non-transiting substellar companion. Again, this might be because the model we apply is not valid and the inclination we obtain from this inference is not correct.

## 4. Discussion and conclusion

Among our 88 candidates, we have 13 F-type, 23 G-type, 37 K-type, and 15 M-type stars exhibiting a signature consistent with the presence of a non-transiting planetary companion in a close orbit. The properties of each candidate are presented in Table C.1. Among these candidates, we detected three systems for which the second harmonic is the largest, indicating that tidal ellipsoidal distortion outweighs the two other phenomena (see Appendix D).

About 62% of our sample (55 stars) exhibit orbital periods of less than a day, confirming the ability of this method to detect very short-period close-in substellar companions. In comparison, we recall again that approximately only 3% of planets detected with the transit method exhibit orbital periods shorter than 1 day.

Within our set of 88 candidates, two stars have been previously confirmed to harbour at least one planet, with an orbital period different from the one detected in this study (they would have been removed from our sample otherwise; see Sect. 2.2). Details of these two candidates are provided in the Appendix E. The method presented in this work could help increase the

<figure><img src="image_6.png" /><figcaption>Fig. 5. $P_{\text{rot}}$ as a function of $P_{\text{orb}}$ for the 88 non-transiting-exoplanet candidate systems. The size of the dots is a function of the Kepler apparent magnitude of the host star, sourced from Brown et al. (2011), whereas the temperature $T_{\text{eff}}$, sourced from Berger et al. (2020), is colour-coded. The grey sample consists of 576 confirmed single-planet-host main-sequence solar-like stars with a rotation period obtained by García et al. (2023). KIC 5697777 is represented by a diamond. The dashed black line corresponds to the fit of the lower envelope of points obtained by McQuillan et al. (2013).</figcaption></figure>

number of multiplanetary systems detected, which in turn can enable the scientific community to gain greater insights into the formation and evolution of exoplanetary systems.

McQuillan et al. (2013, hereafter M13) identified a depletion zone for planets close to fast rotators. García et al. (2023, hereafter G23) demonstrate that a physically motivated model that includes tidal and magnetic effects is able to reproduce the depleted region. Indeed, the amplitude of tidal effects increases dramatically as the semi-major axis gets smaller, meaning that planets in very short orbits decay very fast towards their host star when in an unstable configuration (Hut 1980). In particular, the main component of the tidal potential scales as $1/a^3$ (e.g. Mathis 2018), and under some conditions the orbital decay timescale is proportional to $a^{13/2}$ (Collier Cameron & Jardine 2018).

Fig. 5 displays the stellar rotation period $P_{\text{rot}}$ as a function of the orbital period $P_{\text{orb}}$ for our sample, linked to the findings of G23 in grey. Our 88 candidates are shown as a function of the effective temperature, $T_{\text{eff}}$, sourced from Berger et al. (2020). The dashed black line corresponds to the fit of the lower envelope of points obtained by M13. In this context, we demonstrate that our candidates, if confirmed, would be located well below this limit, thereby populating a previously deserted area in a region where we expect exoplanets to experience orbital decay on relatively short timescales. We remind the reader that the method used in this work strongly favours the detection of planets on very short orbits. In this study, we have focused on very short orbital periods below 2.3 days. If we had been able to detect non-transiting substellar companions with frequencies exceeding our limit of 2.3 days, the region above the dashed black line would contain a higher density of candidates.

Since the orbital evolution of close-in systems takes place over short timescales, studying non-transiting exoplanets would provide us with additional systems to refine our understanding of star-planet interactions. Given the scarcity of available data on planets in short orbits around fast rotators, such a discovery would be a major asset for improving the calibration of models computing the evolution of the architecture of star-planet systems (e.g. Benbakoura et al. 2019; Ahuir et al. 2021). In
addition, by fitting phase curve models to the light curves, we
have been able to extract additional constraints on the mass and
radius for 86 of the detected objects. Subsequent spectroscopic
ground-based follow-up is necessary to further characterise the
sample identified in this work, in order to confirm the plane-
tary nature of the candidates identified in our sample. Radial
velocity observation campaigns could, for example, provide esti-
mates of the masses of the non-transiting companions detected
in this work. Such detections would be very useful to refine the
set of assumptions that are used to model the signal of non-
transiting planets in photometric light curves. In the future, a
similar search for non-transiting objects could be done with the
K2 (Howell et al. 2014), Transiting Exoplanet Survey Satellite
(TESS, Ricker et al. 2015), and PLAnetary Transits and Oscil-
lations of Stars (PLATO, Rauer et al. 2024) missions. Notably,
the fact that PLATO will focus on bright target observations
(Montalto et al. 2021) will greatly expand the number of sys-
tems accessible to ground-based follow-up, allowing for the
detailed characterisation of a much larger sample of exoplanetary
systems.

Data availability

Full Tables B.1 and C.1 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/699/A226

Acknowledgements. C.G., S.N.B., and A.F.L. would like to thank Dr. I. Pagano for enabling C.G.'s research internship at INAF-OACT during Spring 2023. The authors acknowledge the critical reading and the constructive comments from an anonymous referee that significantly allowed improving the original version of this paper. A.S. and C.G. acknowledge funding from the European Research Council project ExoMagnets (grant agreement no. 101125367). A.S. acknowledges funding from the Origins project DynamEarths. A.S., R.G. and S.M. acknowledge funding from the Programme National de Planétologie (INSU/PNP). C.G., S.N.B., and A.F.L acknowledge support from PLATO ASI-INAF agreement no. 2022-28-HH.0 "PLATO Fase D". R.A.G. acknowledges support from PLATO and SoHO/GOLF CNES grants. A.R.G.S acknowledges the support from the FCT through national funds and FEDER through COMPETE2020 (UIDB/04434/2020, UIDP/04434/2020 & 2022.03993.PTDC) and the support from the FCT through the work contract No. 2020.02480.CEECIND/CP1631/CT0001. S.M. acknowledges support by the Spanish Ministry of Science and Innovation with the grant no. PID2019-107061GB-C66 and through AEI under the Severo Ochoa Centres of Excellence Programme 2020-2023 (CEX2019-000920-S). This paper includes data collected by the *Kepler* mission, and obtained from the MAST data archive at the Space Telescope Science Institute (STScI). Funding for the *Kepler* mission is provided by the NASA Science Mission Directorate. STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. This research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. Software: `matplotlib` (Hunter 2007), `numpy` (Harris et al. 2020), `scipy` (Virtanen et al. 2020), `pandas` (The `pandas` development team 2020), `pymc` (Wiecki et al. 2023), `star-privateer` (Breton et al. 2024). `ultranest` (Buchner 2021), `corner` (Foreman-Mackey 2016), `astropy` (Astropy Collaboration 2022), `apollinaire` (Breton et al. 2022).

References

Ahuir, J., Strugarek, A., Brun, A. S., & Mathis, S. 2021, A&A, 650, A126
Akeson, R. L., Chen, X., Ciardi, D., et al. 2013, PASP, 125, 989
Appourchaux, T., Fröhlich, C., Andersen, B., et al. 2000, ApJ, 538, 401
Astropy Collaboration (Price-Whelan, A. M., et al.) 2022, ApJ, 935, 167

Barker, A. J., Efroimsky, M., Makarov, V. V., & Veras, D. 2024, MNRAS, 527, 5131

Benbakoura, M., Réville, V., Brun, A. S., Le Poncin-Lafitte, C., & Mathis, S. 2019, A&A, 621, A124

Berger, T. A., Huber, D., Gaidos, E., & van Saders, J. L. 2018, ApJ, 866, 99

Berger, T. A., Huber, D., van Saders, J. L., et al. 2020, AJ, 159, 280

Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977

Breton, S. N., García, R. A., Ballot, J., Delsanti, V., & Salabert, D. 2022, A&A, 663, A118

Breton, S. N., Lanza, A. F., Messina, S., et al. 2024, A&A, 689, A229

Brown, T. M., Latham, D. W., Everett, M. E., & Esmerdo, G. A. 2011, AJ, 142, 112

Buchner, J. 2021, J. Open Source Softw., 6, 3001

Collier Cameron, A., & Jardine, M. 2018, MNRAS, 476, 2542

Cullen, C. J., & Bayliss, D. 2024, MNRAS, 531, 1133

Esteves, L. J., De Mooij, E. J. W., & Jayawardhana, R. 2013, ApJ, 772, 51

Foreman-Mackey, D. 2016, J. Open Source Softw., 1, 24

Fortney, J. J., Dawson, R. I., & Komacek, T. D. 2021, J. Geophys. Res.: Planets, 126, e2020JE006629

Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1

García, R. A., & Ballot, J. 2019, Liv. Rev. Sol. Phys., 16, 4

García, R. A., Hekker, S., Stello, D., et al. 2011, MNRAS, 414, L6

García, R. A., Mathur, S., Pires, S., et al. 2014, A&A, 568, A10

García, R. A., Gourvès, C., Santos, Â. R. G., et al. 2023, A&A, 679, L12

Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357

Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398

Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90

Hut, P. 1980, A&A, 92, 167

Kjeldsen, H., Christensen-Dalsgaard, J., Handberg, R., et al. 2010, Astron. Nachr., 331, 966

Lillo-Box, J., Millholland, S., & Laughlin, G. 2021, A&A, 654, A9

Liu, Y., San Liang, X., & Weisberg, R. H. 2007, J. Atmos. Ocean. Technol., 24, 2093

Mathis, S. 2018, in Handbook of Exoplanets, eds. H. J. Deeg, & J. A. Belmonte, 24

Mathur, S., Huber, D., Batalha, N. M., et al. 2017, ApJS, 229, 30

Mayor, M., & Queloz, D. 1995, Nature, 378, 355

Mazeh, T. 2008, in EAS Publications Series, 29, eds. M. J. Goupil, & J. P. Zahn, 1

McQuillan, A., Mazeh, T., & Aigrain, S. 2013, ApJ, 775, L11

Messias, Y. S., de Oliveira, L. L. A., Gomes, R. L., et al. 2022, ApJ, 930, L23

Millholland, S., & Laughlin, G. 2017, AJ, 154, 83

Montalto, M., Piotto, G., Marrese, P. M., et al. 2021, A&A, 653, A98

Morton, T. D., Bryson, S. T., Coughlin, J. L., et al. 2016, ApJ, 822, 86

Noraz, Q., Breton, S. N., Brun, A. S., et al. 2022, A&A, 667, A50

Patra, K. C., Winn, J. N., Holman, M. J., et al. 2017, AJ, 154, 4

Penoyre, Z., & Sandford, E. 2019, MNRAS, 488, 4181

Pires, S., Mathur, S., García, R. A., et al. 2015, A&A, 574, A18

Pourbaix, D., Arenou, F., Gavras, P., et al. 2022, Gaia DR3 Documentation, Chapter 7: Non-single stars, https://gea.esac.esa.int/archive/documentation/GDR3/index.html

Rauer, H., Aerts, C., Cabrera, J., et al. 2024, arXiv e-prints, [arXiv:2406.05447]

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telesc. Instrum. Syst., 1, 014003

Santos, Â. R. G., García, R. A., Mathur, S., et al. 2019, ApJS, 244, 21

Santos, Â. R. G., Breton, S. N., Mathur, S., & García, R. A. 2021, ApJS, 255, 17

Seager, S. 2010, Exoplanet Atmospheres: Physical Processes, Princeton Series in Astrophysics (Princeton, N.J: Princeton University Press), oCLC: ocn466343197

Shporer, A. 2017, PASP, 129, 072001

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

Strugarek, A. 2018, in Handbook of Exoplanets, eds. H. J. Deeg, & J. A. Belmonte, 25

The pandas development team 2020, https://doi.org/10.5281/zenodo.3509134

Torrence, C., & Compo, G. P. 1998, Bull. Am. Meteorol. Soc., 79, 61

Torres, G., Kane, S. R., Rowe, J. F., et al. 2017, AJ, 154, 264

Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nat. Methods, 17, 261

Wiecki, T., Salvatier, J., Vieira, R., et al. 2023, https://doi.org/10.5281/zenodo.7961775

Winn, J. N. 2014, arXiv e-prints [arXiv:1001.2010]

Woodard, M. F. 1984, PhD thesis, University of California, San Diego, USA

Zahn, J. P. 1975, A&A, 41, 329
# Appendix A: Contaminants

## A.1. *Transiting exoplanet*

In Sect. 2.2, we were able to remove 4 stars confirmed to harbour transiting exoplanets. For example, Morton et al. (2016) discovered Kepler-685b, a gas giant exoplanet orbiting around KIC 3351888 with an orbital period of 1.62 days.

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

**Fig. A.1.** (a) PSD of KIC 3351888. The $P_{\text{rot}}$, shown in the orange area, was measured by Santos et al. (2021). The green area shows the high-amplitude peak signal detected in this work, at around 1.62 days. (b) Phase-folding analysis on the light curve of KIC 3351888. The red curve is a smoothing of the black curve, obtained with a triangular filter.

Panel (a) of Fig. A.1 shows the PSD of KIC 3351888 where a high-amplitude peak signal at 1.62 days can be observed. In panel b of Fig A.1, the phase-folding curve demonstrates a flux variation analogous to that observed during a transit. The signal detected in this work for KIC 3351888 therefore corresponds to the planet Kepler-685b. Because it has already been detected, this candidate is not part of our final sample.

## A.2. *Phase-folding*

By phase-folding the light curve, we can examine the temporal stability of a signal corresponding to a given period. To refine our sample, we performed a visual inspection of the 283 stars obtained after Sects. 2.1 and 2.2. In these sections, we searched for high-amplitude peaks in the PSDs and compared our results with existing catalogues. In Sect. 2.3, we performed a phase-folding analysis followed by a visual inspection.

Fig. A.2 shows the phase-folded light curve of KIC 11154093, which exhibits a complex pattern with multiple distinct flux variations across the folding period. Visually, the observed pattern deviates significantly from the expected behaviour of a signal (see Fig. 1 and panel c of Fig. 2). This complex behaviour might be linked to the presence of multiple stars close to this target, such as seen when looking at the 2MASS

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

**Fig. A.2.** Phase-folding analysis on the light curve of KIC 11154093. The red curve is a smoothing of the black curve, obtained with a triangular filter.

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

**Fig. A.3.** Image from the KASOC website of KIC 11154093 using 2MASS (Skrutskie et al. 2006). The purple cross marks KIC 11154093.

image (see Appendix A.4) of the target in Fig. A.3. As a result, this target was excluded from our sample.

## A.3. *Wavelet analysis*

In Sect. 2.3, we computed the WPS of each of the 245 stellar light curves, and we only kept stars that present stable modulations along time, as shown in panel d of Fig. 2.

Figure A.4 presents the rotational analysis of KIC 2852669, a K-type star included in our sample of 245 stars. The top panel displays the PSD while the bottom panel shows the normalised time-averaged projection of the wavelet transform.

The different analyses all point to an orbital period of around 0.6 days (panel a, Fig. A.4). However, the signal shown on panel b of Fig. A.4 does not show stable modulations along time. We recall that the focal plane of *Kepler* is rotated after each quarter, so the photometric apertures used to monitor a given target have varying levels of background contamination. The periodic signal found in this study is therefore unlikely to be caused by a non-transiting object, and KIC 3233612 was excluded from our sample.
<figure><img src="image_10.png" /><figcaption>**Fig. A.4.** Rotational analysis of KIC 3233612. (a) PSD of the star. The stellar rotation period is shown in the orange area, while the green area shows the high-amplitude peak at around 0.6 days. (b) Period-time analysis of the stellar light curve. High-amplitude periodicity is highlighted by orange and red, whereas yellow and blue denote low-power amplitude. The yellow discontinuous structure appearing around 0.6 days results from background contamination. The shaded area represents the cone of influence of the reliable periodicities.</figcaption></figure>

### A.4. Field contaminants

In Sect. 2.3, targets with a nearby bright object in their aperture are removed from the sample. We first looked for potential contaminants using the Gaia catalogue (Gaia Collaboration 2023) and found 70 Gaia objects near 38 of our targets. We excluded 11 targets for which the potential contaminants are located within a few arcseconds. The 27 remaining stars are located at least 10" away from the potential contaminant, making significant contaminations problems unlikely. The 10 closest contaminants, associated with these 27 retained targets, are listed in Table A.1.

<figure><img src="image_11.png" /><figcaption>**Fig. A.5.** Image from the KASOC website of KIC 12505054 using 2MASS (Skrutskie et al. 2006). The purple cross marks KIC 12505054.</figcaption></figure>

<table><caption>**Table A.1.** Candidates with a contaminant located at least 10" away and less than 40" away from the considered star.</caption><thead><tr><th>KIC</th><th>Gaia ID</th><th>Contaminant Gaia ID</th><th>distance</th></tr></thead><tbody><tr><td>3222610</td><td>2099534981005352576</td><td>2099535015365093120</td><td>17.88</td></tr><tr><td>5731359</td><td>2073820977455843200</td><td>2073821183614272256</td><td>16.51</td></tr><tr><td>5938531</td><td>2104466423796556160</td><td>2104466423796548992</td><td>19.02</td></tr><tr><td>6370174</td><td>2077612883814892160</td><td>2077614361283648128</td><td>20.54</td></tr><tr><td>7511885</td><td>2105555322564080384</td><td>2105555322564081664</td><td>18.85</td></tr><tr><td>7882472</td><td>2103046030870052224</td><td>2103046030870053760</td><td>18.29</td></tr><tr><td>10027247</td><td>2086316205543022336</td><td>2086316205543023104</td><td>14.29</td></tr><tr><td>10291936</td><td>2086348091380882688</td><td>2086348160100359808</td><td>13.29</td></tr><tr><td>11408355</td><td>2134998079291667968</td><td>2135001034229167232</td><td>11.04</td></tr><tr><td>11605209</td><td>2132673883508778880</td><td>2132673883508775168</td><td>18.45</td></tr><tr><td>...</td><td>...</td><td>...</td><td>...</td></tr></tbody></table>

**Notes.** Only the 10 closest contaminants are shown in this table. Gaia ID is the corresponding Gaia ID of the Kepler star. The distance, measured as the angular separation between the star and the contaminant, is expressed in arcseconds.

We then carried out a visual inspection using the 2MASS images (Skrutskie et al. 2006), which can be accessed through the KASOC database (Kjeldsen et al. 2010). Similarly to what was done e.g. in Noraz et al. (2022), we excluded candidates if the star is within the halo of a brighter object.

Figure A.5 shows KIC 12505054, a star located within the halo of another bright object not listed in the KASOC database. As a result, this star was excluded from our final sample. Indeed, a nearby bright object is likely to contaminate the signal obtained for this star.
**Appendix B: Summary of photometric fit results**

<figure><img src="image_12.png" /><figcaption>**Fig. B.1.** Corner plot for the fit of KIC 4930913 showing the posterior distributions of the fitted parameters $A_0$, $A_{\text{atm}}$ and $T_0$. The values and error bars represent the median and $1\sigma$ uncertainty based on the posterior distributions.</figcaption></figure>

In Sect. 3 we computed the fit of the photometric modulation data and retrieved the theoretical values of $A_{\text{atm}}$. By doing so, we were able to deduce the parameter $\alpha_{\text{refl}} \sin i(R_p/R_j)^2$, thanks to Eq. (3). These resulting parameters are presented in Table B.1 for the 15 first targets. Fig. B.1 shows the corner plot of the posterior probability distribution resulting from the fit of KIC 4930913, a star introduced in Sect. 3 to illustrate the fitting of light curves exhibiting approximately sinusoidal modulation. The diagonal panels display the marginalised one dimensional distribution for each parameter: the baseline amplitude $A_0$, the atmospherical amplitude $A_{\text{atm}}$ and the reference time $T_0$. The off-diagonal panels present the two-dimensional joint posterior distributions, with contour levels corresponding to the $1\sigma$, $2\sigma$, and $3\sigma$ credible regions. The results indicate that all parameters are well constrained. In particular, $T_0$ shows a narrow distribution, implying high precision in the timing estimate. The distribution of $A_0$ is broader and slightly asymmetric, reflecting larger uncertainties. No significant correlations are observed between the parameters, suggesting they are independently constrained by the data. We also calculated with Eq. (9) the maximum inclination possible for each system, assuming they all follow a circular orbit. For two stars (KIC 3124224 and KIC 10396339), we do not provide $i_{\text{max}}$, because, using $R_\star$ and $M_\star$ from the reference catalogues adopted, we obtain $a < R_\star$, which is obviously unphysical. Nevertheless, considering that $R_\star$ and $M_\star$ might be incorrect, we decided to keep these two stars in our sample.

When the modulation of the phase-folded light curve was more complex, we performed a fit considering the three components of the signal. The corner plot in Fig. B.2 summarises the posterior probability distributions resulting for KIC 4264634, the example star
<figure><img src="image_13.png" /><figcaption>**Fig. B.2.** Corner plot for the fit of KIC 4264634 showing the posterior distributions of the fitted parameters $A_0$, $A_{\text{atm}}$, $A_{\text{ellip}}$, $A_{\text{boost}}$ and $T_0$. The values and error bars represent the median and $1\sigma$ uncertainty based on the posterior distributions.</figcaption></figure>

chosen in Sect. 3 to illustrate the fit which accounts for the three components of the signal. As in Fig. B.1, the diagonal panels show the marginalised one-dimensional distributions for each fitted parameter: the baseline amplitude, $A_0$, the atmospheric amplitude, $A_{\text{atm}}$, the reference time, $T_0$, the ellipsoidal amplitude, $A_{\text{ellip}}$, and the Doppler boosting amplitude, $A_{\text{boost}}$. The off-diagonal panels illustrate the two-dimensional joint posterior distributions, with contours corresponding to the $1\sigma$, $2\sigma$, and $3\sigma$ credible intervals. All parameters are reasonably well constrained, with the amplitude parameters displaying broader distributions, especially $A_0$ and $A_{\text{atm}}$, which exhibit mild asymmetries. The joint distributions suggest weak correlations between some of the amplitude parameters, particularly between $A_{\text{atm}}$ and $A_{\text{ellip}}$. Overall, the parameters remain largely uncorrelated. This supports the robustness of the fit and the ability of the model to independently constrain each physical contribution. The resulting parameters for all of the nine targets studied here are presented in Table B.2. In both Tables B.1 and B.2, the reported values of $A_{\text{atm}}$, $A_{\text{beam}}$ and $A_{\text{ellip}}$ correspond to their absolute values. Thanks to Eq. (5) and (6) we were also able to determine the parameters ($M_p \sin^2 i/M_J$) and ($M_p \sin i/M_J$). With such calculations, we can in principle retrieve the value of inclination of the non-transiting object *i*. This was possible for three of our systems (KIC 11550689, 5697777, and 4264634). For three other stars, the uncertainties on our calculated parameters are too large to infer on the value *i*. For the last three stars, the uncertainties are not sufficiently high enough to explain the obtained inconsistent results. We are therefore led to believe that these stars are out of the validity range of the model used in this article and would probably require more complex models (see Penoyre & Sandford 2019).
<table><caption>**Table B.1.** Photometric fit results of 15 out of the 77 candidates showing quasi-sinusoidal modulations</caption><thead><tr><th>KIC</th><th>A<sub>atm</sub> [ppm]</th><th>α<sub>refl</sub> sin i(R<sub>p</sub>/R<sub>j</sub>)<sup>2</sup></th><th>i<sub>max</sub> [°]</th></tr></thead><tbody><tr><td>2707504</td><td>185.91<sup>+0.63</sup><sub>-0.64</sub></td><td>6.3026<sup>+0.0214</sup><sub>-0.0216</sub></td><td>82</td></tr><tr><td>3124224</td><td>72.51<sup>+0.74</sup><sub>-0.76</sub></td><td>0.1416<sup>+0.0014</sup><sub>-0.0015</sub></td><td>—</td></tr><tr><td>3222610</td><td>133.07<sup>+0.96</sup><sub>-0.96</sub></td><td>1.8074<sup>+0.0130</sup><sub>-0.0131</sub></td><td>80</td></tr><tr><td>3757251</td><td>233.86<sup>+2.48</sup><sub>-2.66</sub></td><td>0.7963<sup>+0.0084</sup><sub>-0.0091</sub></td><td>73</td></tr><tr><td>3954961</td><td>218.78<sup>+1.24</sup><sub>-1.18</sub></td><td>3.2079<sup>+0.0181</sup><sub>-0.0174</sub></td><td>79</td></tr><tr><td>4175216</td><td>279.84<sup>+0.80</sup><sub>-0.78</sub></td><td>2.9049<sup>+0.0083</sup><sub>-0.0081</sub></td><td>78</td></tr><tr><td>4283417</td><td>190.32<sup>+1.83</sup><sub>-1.82</sub></td><td>2.4066<sup>+0.0231</sup><sub>-0.0230</sub></td><td>80</td></tr><tr><td>4641549</td><td>92.29<sup>+2.31</sup><sub>-2.34</sub></td><td>1.9074<sup>+0.0477</sup><sub>-0.0484</sub></td><td>79</td></tr><tr><td>4859400</td><td>40.74<sup>+0.75</sup><sub>-0.69</sub></td><td>0.7340<sup>+0.0135</sup><sub>-0.0125</sub></td><td>67</td></tr><tr><td>4862155</td><td>65.28<sup>+0.73</sup><sub>-0.72</sub></td><td>1.4593<sup>+0.0163</sup><sub>-0.0160</sub></td><td>72</td></tr><tr><td>4916039</td><td>93.12<sup>+1.87</sup><sub>-1.84</sub></td><td>0.4643<sup>+0.0093</sup><sub>-0.0092</sub></td><td>71</td></tr><tr><td>4930913</td><td>98.32<sup>+0.84</sup><sub>-0.84</sub></td><td>1.2077<sup>+0.0103</sup><sub>-0.0103</sub></td><td>73</td></tr><tr><td>4945547</td><td>60.85<sup>+0.97</sup><sub>-1.00</sub></td><td>0.4127<sup>+0.0066</sup><sub>-0.0068</sub></td><td>74</td></tr><tr><td>4996173</td><td>39.76<sup>+0.14</sup><sub>-0.14</sub></td><td>0.3498<sup>+0.0012</sup><sub>-0.0012</sub></td><td>49</td></tr><tr><td>5006208</td><td>52.85<sup>+0.71</sup><sub>-0.62</sub></td><td>2.6370<sup>+0.0356</sup><sub>-0.0312</sub></td><td>75</td></tr><tr><td>...</td><td>...</td><td>...</td><td>...</td></tr></tbody></table>

**Notes.** The full table is available at the CDS.

<table><caption>**Table B.2.** Photometric fit results for the nine candidates showing complex modulations</caption><thead><tr><th>KIC</th><th>A<sub>atm</sub> [ppm]</th><th>α<sub>refl</sub> sin i(R<sub>p</sub>/R<sub>j</sub>)<sup>2</sup></th><th>A<sub>boost</sub> [ppm]</th><th>A<sub>ellip</sub> [ppm]</th><th>(M<sub>p</sub> sin i/M<sub>j</sub>)</th><th>(M<sub>p</sub> sin<sup>2</sup> i/M<sub>j</sub>)</th><th>i<sub>max</sub> [°]</th></tr></thead><tbody><tr><td>4264634</td><td>134.15<sup>+1.41</sup><sub>-1.38</sub></td><td>1.3686<sup>+0.0144</sup><sub>-0.0141</sub></td><td>258.19<sup>+1.34</sup><sub>-1.29</sub></td><td>197.59<sup>+1.16</sup><sub>-1.22</sub></td><td>12.97<sup>+0.07</sup><sub>-0.06</sub></td><td>8.02<sup>+0.05</sup><sub>-0.05</sub></td><td>78</td></tr><tr><td>4373708</td><td>301.70<sup>+1.31</sup><sub>-1.41</sub></td><td>1.4014<sup>+0.0061</sup><sub>-0.0065</sub></td><td>9.12<sup>+1.84</sup><sub>-1.79</sub></td><td>163.97<sup>+1.33</sup><sub>-1.24</sub></td><td>0.45<sup>+0.09</sup><sub>-0.09</sub></td><td>1.44<sup>+0.01</sup><sub>-0.01</sub></td><td>69</td></tr><tr><td>5440365</td><td>0.73<sup>+0.85</sup><sub>-0.86</sub></td><td>0.0156<sup>+0.0181</sup><sub>-0.0184</sub></td><td>0.30<sup>+0.84</sup><sub>-0.83</sub></td><td>31.52<sup>+0.85</sup><sub>-0.85</sub></td><td>0.03<sup>+0.09</sup><sub>-0.09</sub></td><td>1.81<sup>+0.05</sup><sub>-0.05</sub></td><td>76</td></tr><tr><td>5622796</td><td>2.42<sup>+0.79</sup><sub>-0.73</sub></td><td>0.1068<sup>+0.0348</sup><sub>-0.0320</sub></td><td>2.32<sup>+0.79</sup><sub>-0.81</sub></td><td>45.29<sup>+0.82</sup><sub>-0.75</sub></td><td>0.28<sup>+0.09</sup><sub>-0.10</sub></td><td>9.33<sup>+0.17</sup><sub>-0.15</sub></td><td>81</td></tr><tr><td>5697777</td><td>7.42<sup>+0.62</sup><sub>-0.60</sub></td><td>0.0766<sup>+0.0064</sup><sub>-0.0062</sub></td><td>78.80<sup>+0.47</sup><sub>-0.44</sub></td><td>41.29<sup>+0.47</sup><sub>-0.42</sub></td><td>9.89<sup>+0.06</sup><sub>-0.06</sub></td><td>0.07<sup>+0.00</sup><sub>-0.00</sub></td><td>16</td></tr><tr><td>8844761</td><td>0.57<sup>+0.60</sup><sub>-0.55</sub></td><td>0.0049<sup>+0.0052</sup><sub>-0.0048</sub></td><td>0.30<sup>+0.55</sup><sub>-0.57</sub></td><td>25.94<sup>+0.56</sup><sub>-0.57</sub></td><td>0.01<sup>+0.03</sup><sub>-0.03</sub></td><td>0.89<sup>+0.02</sup><sub>-0.02</sub></td><td>78</td></tr><tr><td>11550689</td><td>2.09<sup>+1.44</sup><sub>-1.49</sub></td><td>0.0053<sup>+0.0036</sup><sub>-0.0037</sub></td><td>88.35<sup>+1.21</sup><sub>-1.22</sub></td><td>62.71<sup>+1.27</sup><sub>-1.13</sub></td><td>3.31<sup>+0.05</sup><sub>-0.05</sub></td><td>0.34<sup>+0.01</sup><sub>-0.01</sub></td><td>66</td></tr><tr><td>11668095</td><td>3.17<sup>+2.39</sup><sub>-2.44</sub></td><td>0.0823<sup>+0.0620</sup><sub>-0.0635</sub></td><td>2.26<sup>+2.46</sup><sub>-2.48</sub></td><td>179.92<sup>+2.35</sup><sub>-2.34</sub></td><td>0.14<sup>+0.15</sup><sub>-0.15</sub></td><td>39.18<sup>+0.51</sup><sub>-0.51</sub></td><td>83</td></tr><tr><td>11702835</td><td>6.82<sup>+1.40</sup><sub>-1.51</sub></td><td>0.0355<sup>+0.0073</sup><sub>-0.0079</sub></td><td>35.77<sup>+1.49</sup><sub>-1.50</sub></td><td>156.55<sup>+1.43</sup><sub>-1.56</sub></td><td>1.74<sup>+0.07</sup><sub>-0.07</sub></td><td>2.93<sup>+0.03</sup><sub>-0.03</sub></td><td>74</td></tr></tbody></table>
## Appendix C: Final sample

<table><caption>Table C.1. Global parameters of 16 out of the 88 candidates likely to host non-transiting companions.</caption><thead><tr><th rowspan="2">KIC</th><th rowspan="2">P<sub>orb</sub><br>[days]</th><th colspan="2">semi-major axis</th><th rowspan="2">A<sub>periodogram</sub><br>[ppm]</th><th rowspan="2">P<sub>rot</sub><br>[days]</th><th rowspan="2">T<sub>eff</sub><br>[K]</th><th rowspan="2">Mass<br>[M<sub>☉</sub>]</th><th rowspan="2">Kpmag<br>[mag]</th><th rowspan="2">log g<br>[dex]</th><th rowspan="2">[Fe/H]<br>[dex]</th><th rowspan="2">flag</th></tr><tr><th>[AU]</th><th></th></tr></thead><tbody><tr><td>2707504</td><td>1.83354<sup>+0.00110</sup><sub>-0.00112</sub></td><td>0.02721</td><td><sup>+0.00019</sup><sub>-0.00020</sub></td><td>202<sup>+59</sup><sub>-39</sub></td><td>18.3<sup>+1.9</sup><sub>-1.9</sub></td><td>5135<sup>+179</sup><sub>-179</sub></td><td>0.80<sup>+0.09</sup><sub>-0.07</sub></td><td>13.68</td><td>4.56<sup>+0.05</sup><sub>-0.07</sub></td><td>-0.06<sup>+0.30</sup><sub>-0.30</sub></td><td>0, 0, 0</td></tr><tr><td>3124224</td><td>0.16550<sup>+0.00001</sup><sub>-0.00001</sub></td><td>0.01037</td><td><sup>+0.00001</sup><sub>-0.00001</sub></td><td>73<sup>+22</sup><sub>-13</sub></td><td>74.4<sup>+6.5</sup><sub>-6.5</sub></td><td>5206<sup>+128</sup><sub>-122</sub></td><td>1.36<sup>+0.04</sup><sub>-0.06</sub></td><td>13.91</td><td>3.77<sup>+0.02</sup><sub>-0.03</sub></td><td>0.28<sup>+0.16</sup><sub>-0.15</sub></td><td>0, 0, 0</td></tr><tr><td>3222610</td><td>1.01284<sup>+0.00038</sup><sub>-0.00028</sub></td><td>0.01722</td><td><sup>+0.00009</sup><sub>-0.00007</sub></td><td>144<sup>+90</sup><sub>-43</sub></td><td>4.8<sup>+0.3</sup><sub>-0.3</sub></td><td>4260<sup>+47</sup><sub>-70</sub></td><td>0.66<sup>+0.02</sup><sub>-0.02</sub></td><td>14.42</td><td>4.62<sup>+0.03</sup><sub>-0.02</sub></td><td>0.06<sup>+0.12</sup><sub>-0.12</sub></td><td>0, 0, 0</td></tr><tr><td>3757251</td><td>0.44853<sup>+0.00008</sup><sub>-0.00009</sub></td><td>0.00862</td><td><sup>+0.00003</sup><sub>-0.00003</sub></td><td>228<sup>+46</sup><sub>-36</sub></td><td>15.8<sup>+1.3</sup><sub>-1.3</sub></td><td>3624<sup>+82</sup><sub>-86</sub></td><td>0.42<sup>+0.01</sup><sub>-0.02</sub></td><td>15.77</td><td>4.71<sup>+0.01</sup><sub>-0.01</sub></td><td>0.05<sup>+0.12</sup><sub>-0.13</sub></td><td>0, 0, 0</td></tr><tr><td>3954961</td><td>1.04475<sup>+0.00023</sup><sub>-0.00021</sub></td><td>0.01789</td><td><sup>+0.00007</sup><sub>-0.00006</sub></td><td>238<sup>+74</sup><sub>-47</sub></td><td>13.1<sup>+1.0</sup><sub>-1.0</sub></td><td>4413<sup>+70</sup><sub>-61</sub></td><td>0.70<sup>+0.03</sup><sub>-0.02</sub></td><td>14.84</td><td>4.59<sup>+0.02</sup><sub>-0.01</sub></td><td>0.20<sup>+0.11</sup><sub>-0.11</sub></td><td>0, 0, 0</td></tr><tr><td>4175216</td><td>0.93558<sup>+0.00035</sup><sub>-0.00038</sub></td><td>0.01506</td><td><sup>+0.00008</sup><sub>-0.00008</sub></td><td>272<sup>+41</sup><sub>-30</sub></td><td>13.7<sup>+1.1</sup><sub>-1.1</sub></td><td>4327<sup>+74</sup><sub>-69</sub></td><td>0.52<sup>+0.03</sup><sub>-0.03</sub></td><td>13.82</td><td>4.62<sup>+0.02</sup><sub>-0.02</sub></td><td>0.01<sup>+0.11</sup><sub>-0.10</sub></td><td>0, 0, 0</td></tr><tr><td>4264634</td><td>0.88671<sup>+0.00038</sup><sub>-0.00051</sub></td><td>0.01493</td><td><sup>+0.00008</sup><sub>-0.00010</sub></td><td>316<sup>+134</sup><sub>-71</sub></td><td>12.3<sup>+0.8</sup><sub>-0.8</sub></td><td>3940<sup>+70</sup><sub>-94</sub></td><td>0.56<sup>+0.02</sup><sub>-0.02</sub></td><td>14.82</td><td>4.61<sup>+0.03</sup><sub>-0.02</sub></td><td>0.16<sup>+0.12</sup><sub>-0.13</sub></td><td>0, 0, 0</td></tr><tr><td>4283417</td><td>0.99830<sup>+0.00061</sup><sub>-0.00057</sub></td><td>0.01662</td><td><sup>+0.00012</sup><sub>-0.00011</sub></td><td>185<sup>+41</sup><sub>-32</sub></td><td>13.2<sup>+1.0</sup><sub>-1.0</sub></td><td>4117<sup>+68</sup><sub>-68</sub></td><td>0.61<sup>+0.02</sup><sub>-0.02</sub></td><td>15.28</td><td>4.63<sup>+0.03</sup><sub>-0.03</sub></td><td>0.06<sup>+0.12</sup><sub>-0.11</sub></td><td>0, 0, 0</td></tr><tr><td>4373708</td><td>0.46245<sup>+0.00013</sup><sub>-0.00014</sub></td><td>0.01007</td><td><sup>+0.00004</sup><sub>-0.00004</sub></td><td>293<sup>+43</sup><sub>-33</sub></td><td>9.4<sup>+0.7</sup><sub>-0.7</sub></td><td>4802<sup>+86</sup><sub>-80</sub></td><td>0.64<sup>+0.04</sup><sub>-0.03</sub></td><td>14.83</td><td>4.53<sup>+0.03</sup><sub>-0.03</sub></td><td>0.09<sup>+0.12</sup><sub>-0.12</sub></td><td>0, 0, 0</td></tr><tr><td>4641549</td><td>1.18275<sup>+0.00034</sup><sub>-0.00038</sub></td><td>0.02124</td><td><sup>+0.00009</sup><sub>-0.00010</sub></td><td>101<sup>+63</sup><sub>-26</sub></td><td>21.2<sup>+1.7</sup><sub>-1.7</sub></td><td>5662<sup>+110</sup><sub>-103</sub></td><td>0.91<sup>+0.05</sup><sub>-0.05</sub></td><td>15.76</td><td>4.51<sup>+0.03</sup><sub>-0.04</sub></td><td>-0.14<sup>+0.14</sup><sub>-0.16</sub></td><td>0, 0, 0</td></tr><tr><td>4859400</td><td>0.97366<sup>+0.00034</sup><sub>-0.00027</sub></td><td>0.01984</td><td><sup>+0.00010</sup><sub>-0.00008</sub></td><td>44<sup>+13</sup><sub>-8</sub></td><td>21.5<sup>+2.7</sup><sub>-2.7</sub></td><td>5875<sup>+134</sup><sub>-131</sub></td><td>1.10<sup>+0.12</sup><sub>-0.08</sub></td><td>14.00</td><td>4.02<sup>+0.06</sup><sub>-0.05</sub></td><td>-0.06<sup>+0.15</sup><sub>-0.16</sub></td><td>0, 0, 0</td></tr><tr><td>4862155</td><td>1.04625<sup>+0.00044</sup><sub>-0.00042</sub></td><td>0.02209</td><td><sup>+0.00012</sup><sub>-0.00012</sub></td><td>71<sup>+18</sup><sub>-13</sub></td><td>9.5<sup>+0.6</sup><sub>-0.6</sub></td><td>6574<sup>+150</sup><sub>-141</sub></td><td>1.31<sup>+0.06</sup><sub>-0.06</sub></td><td>13.91</td><td>4.23<sup>+0.04</sup><sub>-0.04</sub></td><td>0.01<sup>+0.14</sup><sub>-0.16</sub></td><td>0, 0, 0</td></tr><tr><td>4916039</td><td>0.49881<sup>+0.00012</sup><sub>-0.00013</sub></td><td>0.01043</td><td><sup>+0.00004</sup><sub>-0.00004</sub></td><td>101<sup>+27</sup><sub>-24</sub></td><td>32.4<sup>+2.7</sup><sub>-2.7</sub></td><td>4809<sup>+97</sup><sub>-91</sub></td><td>0.61<sup>+0.04</sup><sub>-0.04</sub></td><td>15.44</td><td>4.59<sup>+0.03</sup><sub>-0.03</sub></td><td>-0.05<sup>+0.13</sup><sub>-0.13</sub></td><td>0, 0, 0</td></tr><tr><td>4930913</td><td>0.80635<sup>+0.00019</sup><sub>-0.00022</sub></td><td>0.01638</td><td><sup>+0.00006</sup><sub>-0.00007</sub></td><td>95<sup>+134</sup><sub>-53</sub></td><td>16.4<sup>+1.9</sup><sub>-1.9</sub></td><td>5726<sup>+77</sup><sub>-77</sub></td><td>0.90<sup>+0.07</sup><sub>-0.05</sub></td><td>13.96</td><td>4.35<sup>+0.13</sup><sub>-0.11</sub></td><td>-0.12<sup>+0.15</sup><sub>-0.15</sub></td><td>0, 0, 0</td></tr><tr><td>4945547</td><td>0.56461<sup>+0.00016</sup><sub>-0.00015</sub></td><td>0.01217</td><td><sup>+0.00005</sup><sub>-0.00005</sub></td><td>66<sup>+22</sup><sub>-17</sub></td><td>24.9<sup>+2.0</sup><sub>-2.0</sub></td><td>4954<sup>+98</sup><sub>-95</sub></td><td>0.75<sup>+0.03</sup><sub>-0.04</sub></td><td>14.46</td><td>4.59<sup>+0.02</sup><sub>-0.03</sub></td><td>-0.16<sup>+0.12</sup><sub>-0.13</sub></td><td>0, 0, 0</td></tr><tr><td>4996173</td><td>0.47867<sup>+0.00008</sup><sub>-0.00008</sub></td><td>0.01386</td><td><sup>+0.00004</sup><sub>-0.00004</sub></td><td>43<sup>+12</sup><sub>-9</sub></td><td>9.1<sup>+0.6</sup><sub>-0.6</sub></td><td>6712<sup>+142</sup><sub>-140</sub></td><td>1.55<sup>+0.07</sup><sub>-0.07</sub></td><td>10.53</td><td>4.04<sup>+0.03</sup><sub>-0.03</sub></td><td>0.17<sup>+0.13</sup><sub>-0.16</sub></td><td>0, 0, 0</td></tr><tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr></tbody></table>

Notes. The full table is available at the CDS.

We present in Table C.1 the overall properties of 16 out of our 88 selected candidates likely to host non-transiting companions. The orbital period and amplitude are derived from this work. The effective temperature $T_{\text{eff}}$, surface gravity log g, metallicity and mass measurements parameters are primarily sourced from Berger et al. (2020) while, for stars not included in this catalogue, the properties are taken from Mathur et al. (2017). The rotation periods are taken from Santos et al. (2019, 2021) and the Kepler apparent magnitudes Kpmag come from Brown et al. (2011). The semi-major axis is calculated with the Kepler's third law, assuming the planetary mass is negligible compared to the stellar mass.

The flag column contains a sequence of three binary values (0 or 1) for each entry. Each digit represents a separate flag. If the first value is set to 1, it indicates that tidal forces outweighs atmospheric processes and Doppler boosting in the orbital signal. The second value, if set to 1, means that the target's orbital period is a harmonic of the stellar rotation period (see Appendix D). Finally, the third value indicates whether an exoplanet has already been found orbiting around the targeted star, but at a different orbital period than the one detected in this work (see Appendix E). In such cases, the value is set to 1. The most common set of flags is 0, 0,0, indicating that the majority of the sample does not exhibit these characteristics.

We estimate the uncertainties on the orbital period $P_{\text{orb}}$ and the periodogram amplitude $A_{\text{periodogram}}$ of the modulation by modelling the main peak attributed in the PSD to the non-transiting planet with a Gaussian profile with central frequency $v_0 = 1/P_{\text{orb}}$, height $H$, and standard deviation $\sigma$. Using the pyMC software (Wiecki et al. 2023), we perform a Hamiltonian Monte-Carlo exploration of the distribution of the set $\theta = \{P_{\text{orb}}, H, \sigma\}$, assuming a $\chi^2$ with two-degrees of freedom as likelihood function for the PSD (Woodard 1984). The periodogram amplitude $A_{\text{periodogram}}$ of the modulation, in ppm, is computed as $A^2_{\text{periodogram}} = 2\sqrt{\pi}\sigma H$. The parameter estimate is taken as the median of the distribution. The lower uncertainty is computed as the difference between the median value of the distribution and the 16th percentile, while the upper uncertainty is computed as the difference between the median value of the distribution and the 84th percentile.
## Appendix D: Systems with an orbital period corresponding to a harmonic of the stellar rotation period

In our 88 stars catalogue, the orbital period of 3 candidates corresponds to a harmonic of the stellar surface rotation period.

<figure><img src="image_14.png" /><figcaption>**Fig. D.1.** (a) PSD of KIC 7801559. The $P_{\text{rot}}$, shown in the orange area, was measured by Santos et al. (2021). The green area shows the high-amplitude signal peak at around 1.34 days. (b) Phase-folding analysis on the light curve of KIC 7801559. The red curve is a smoothing of the black curve, obtained with a triangular filter.</figcaption></figure>

<table><caption>**Table D.1.** Candidates with an orbital period corresponding to a harmonic of the stellar rotation period</caption><thead><tr><th>KIC</th><th><i>P</i><sub>orb</sub> [days]</th><th><i>P</i><sub>rot</sub> [days]</th><th>Harmonic</th></tr></thead><tbody><tr><td>5560359</td><td>1.96325<sup>+0.00081</sup><sub>-0.00073</sub></td><td>3.9<sup>+0.5</sup><sub>-0.5</sub></td><td>2</td></tr><tr><td>7801559</td><td>1.37431<sup>+0.00076</sup><sub>-0.00051</sub></td><td>2.8<sup>+0.2</sup><sub>-0.2</sub></td><td>2</td></tr><tr><td>8193786</td><td>2.18821<sup>+0.00120</sup><sub>-0.00149</sub></td><td>6.4<sup>+0.6</sup><sub>-0.6</sub></td><td>3</td></tr></tbody></table>

**Notes.** The rotation periods are taken from the catalogues by Santos et al. (2019, 2021).

Panel (a) of Fig. D.1 presents the PSD of KIC 7801559, an F-type star with a $P_{\text{rot}}$ of 2.8 days. The companion's orbital period detected in this work corresponds to the second harmonic of the stellar rotation period. On panel (a), the high-amplitude peak corresponding to the orbital signal (in the green area), seems embedded in the rotational harmonic pattern. Nevertheless, the signal shows a very stable periodic behaviour, distinct from the quasi-periodic signature observed in the case of spot modulation, both in the PSD and the light curve phase-folding in panel (b) of Fig. D.1. Therefore, we choose to keep this candidate in our catalogue and suggest that this could be the signature of systems in a spin-orbit resonance. This verification was also done

for KIC 5560359 and KIC 8193786. The companion orbital period and the stellar rotation period of these candidates are shown in Table D.

## Appendix E: Systems with a confirmed exoplanet

In our 88 stars sample, we detected 2 stars with previously confirmed exoplanets, but at an orbital period different from the one detected in this study. The first star is KIC 10027247, an M-type star with a super-Earth exoplanet, named Kepler-1229 b (Morton et al. 2016). Its mass is 2.54 $M_{\oplus}$, and it takes 86.8 days to complete one orbit of its star at a distance of 0.3006 Astronomical Unit (AU). Notably, its orbit is within the habitable zone (Torres et al. 2017). The second star in our catalogue that already hosts a confirmed exoplanet is KIC 9717943, an F-type star. Kepler-1169 b is a terrestrial exoplanet with a mass of 0.779 $M_{\oplus}$ (Morton et al. 2016). Its orbital period is 6.11 days, and it is located at 0.0709 AU from its host star. The properties of these two systems are shown in Table E.

<table><caption>**Table E.1.** Stellar parameters of two candidates with an already confirmed exoplanet.</caption><thead><tr><th>KIC</th><th>9717943</th><th>10027247</th></tr></thead><tbody><tr><td><i>T</i><sub>eff</sub> [K]</td><td>6000<sup>+106</sup><sub>-105</sub></td><td>3576<sup>+75</sup><sub>-73</sub></td></tr><tr><td>Stellar mass [M<sub>☉</sub>]</td><td>1.17<sup>+0.06</sup><sub>-0.07</sub></td><td>0.48<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td>Kpmag [mag]</td><td>12.72</td><td>15.47</td></tr><tr><td>log <i>g</i> [dex]</td><td>4.34<sup>+0.03</sup><sub>-0.03</sub></td><td>4.67<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td>[Fe/H] [dex]</td><td>0.21<sup>+0.14</sup><sub>-0.15</sub></td><td>0.21<sup>+0.12</sup><sub>-0.12</sub></td></tr><tr><td>Amplitude [ppm]</td><td>35<sup>+16</sup><sub>-11</sub></td><td>126<sup>+25</sup><sub>-18</sub></td></tr><tr><td><i>P</i><sub>rot</sub> [days]</td><td>12.5<sup>+1.6</sup><sub>-1.6</sub></td><td>17.5<sup>+1.5</sup><sub>-1.5</sub></td></tr><tr><td>semi-major axis [AU]</td><td>0.03092<sup>+0.00020</sup><sub>-0.00020</sub></td><td>0.01078<sup>+0.00004</sup><sub>-0.00004</sub></td></tr><tr><td><i>P</i><sub>orb</sub> [days]</td><td>1.83809<sup>+0.00097</sup><sub>-0.00097</sub></td><td>0.58875<sup>+0.00013</sup><sub>-0.00013</sub></td></tr><tr><td>Planet name</td><td>Kepler-1169 b</td><td>Kepler-1229 b</td></tr><tr><td>Confirmed planet's <i>P</i><sub>orb</sub> [days]</td><td>6.11009<sup>+0.00003</sup><sub>-0.00003</sub></td><td>86.82899<sup>+0.00107</sup><sub>-0.00107</sub></td></tr></tbody></table>

**Notes.** The $T_{\text{eff}}$, log *g*, metallicity, and mass measurements are taken from Mathur et al. (2017). The rotation periods and *Kepler* apparent magnitudes are taken from Santos et al. (2019, 2021) while the orbital periods of the confirmed planets are from Morton et al. (2016). The semi-major axis is calculated with the Kepler's third law, assuming the planetary mass is negligible compared to the stellar mass. The $P_{\text{orb}}$ value was determined in this work for a newly identified, non-transiting object.