Planets around evolved intermediate-mass stars III. Planet candidates and long-term activity signals in six open clusters

[abridged]The aim of this work is to search for planets around evolved stars, with a special focus on stars more massive than 2\,M$_\odot$ in light of previous findings that show a drop in planet occurrence around stars above this mass. We used \texttt{kima} to find the Keplerian orbits most capable of explaining the periodic signals observed in RV data. We also studied the variation of stellar activity indicators and photometry in order to discard stellar signals mimicking the presence of planets. We present a planet candidate in the open cluster NGC3680 that orbits the 1.64\,M$_\odot$ star No. 41. The planet has a minimum mass of 5.13M\,$_{J}$ and a period of 1155 days. We also present periodic and large-amplitude RV signals of probable stellar origin in two more massive stars (5.84 and 3.05\,M$_\odot$ in the clusters NGC2345 and NGC3532). Finally, using new data, we revise the RV signals of the three stars analysed in our previous paper. We confirm the stellar origin of the signals observed in NGC2423 No. 3 and NGC4349 No. 127. On the other hand, the new data collected for IC4651 No. 9122 (1.79\,M$_\odot$) seem to support the presence of a bona fide planet of 6.22M\,$_{J}$ at a period of 744 days, although more data will be needed to discard a possible correlation with the CCF-FWHM. The targets presented in this work showcase the difficulties in interpreting RV data for evolved massive stars. The use of several activity indicators (CCF-FWHM, CCF-BIS, \ha), photometry, and long-term observations (covering several orbital and stellar rotational periods) is required to discern the true nature of the signals. However, in some cases, all this information is insufficient, and the inclusion of additional data -- such as the determination of magnetic field variability or RV points in the near-infrared -- will be necessary to identify the nature of the discovered signals.


Introduction
In the last 30 years, more than 5000 planets have been discovered, mainly around main sequence (MS) solar-type stars 1 , that is, FGK dwarf stars and M dwarfs.However, the number of plan- ⋆⋆ RV tables are only available in electronic form at the CDS via anonymous ftp to cdsarc.cds.unistra.fr(130.79.128.5) or via https://cdsarc.cds.unistra.fr/cgi-bin/qcat?J/A+A/ 1 exoplanets.euets detected orbiting intermediate-mass stars is still low despite the increasing frequency of giant planets with stellar mass (e.g.Johnson et al. 2010).Moreover, several studies point to a sharp decrease in the planet occurrence around stars more massive than ∼2 M ⊙ (Reffert et al. 2015;Wolthoff et al. 2022).The search for planets around these more massive stars is crucial to our understanding of the limits of planet formation and survival, but we are still limited by the applicability of different detection techniques to this kind of star.For example, the most successful planet-detection techniques (based on photometric transits and radial velocity (RV)) are not optimal for large, massive, and fastrotating stars, such as early-type MS stars.Nevertheless, both methods have been applied to late-A and early-F stars, revealing a handful of substellar companions (e.g.Desort et al. 2008;Borgniet et al. 2019;Grandjean et al. 2023;Collier Cameron et al. 2010;Sebastian et al. 2022;Vowell et al. 2023).The largest fraction of planets around intermediate-mass stars were detected by successfully applying the abovementioned techniques to K giants, the evolved counterparts of early-type MS stars.Photometric data from the Kepler and TESS missions led to the discovery of a good number of planets around subgiants and red giant branch (RGB) stars (e.g.Lillo-Box et al. 2014;Grunblatt et al. 2022).On the other hand, several longterm RV surveys are being carried out by different teams, such as the Lick survey (Frink et al. 2001), the Okayama Planet Search Program (Sato 2005), the Tautenburg Observatory Planet Search (Hatzes et al. 2005), and the PTPS-TAPAS program (Niedzielski et al. 2015).We refer the reader to Table 1 in Ottoni et al. (2022) for a exhaustive compilation of the current RV surveys around giant stars.Nevertheless, the occurrence rates found by some of these surveys can also be affected by selection effects.Because redder stars tend to show larger RV jitter (e.g.Frink et al. 2001), it is common to apply a cutoff for a maximum B − V.As a consequence, the most metal-rich stars with low log g values are left out of such surveys (Mortier et al. 2013).
Alternatively, the direct imaging technique can be applied to early-type MS stars, especially those with larger masses and young ages.However, due to the inherent biases affecting the direct imaging technique, only substellar companions at large distances can be detected.The largest surveys to date, SHINE (with SPHERE@VLT) and GPIES (with GPI@Gemini South), only found a few substellar companions around A stars (Vigan et al. 2021;Nielsen et al. 2019, respectively).Indeed, most of the above-mentioned surveys (either around dwarf or evolved stars) only found detected planets around stars below ∼2.5 M ⊙ (i.e.spectral type later than A0-A1).Some theoretical studies predict that planet formation for stars more massive than 3 M ⊙ is very difficult because in such stars irradiation overcomes accretion and the snowlines lie at greater distances, hindering the formation of cores (Kennedy & Kenyon 2008).However, the very recent results of the BEAST survey, wherein the direct imaging technique is applied to B-type stars (M ≳2.5 M ⊙ ), show that substellar companions (most likely brown dwarfs) can form around this kind of star, which is contradictory to the low occurrence rates found in RV surveys of evolved stars.A brown dwarf was reported in an orbit of 290 AU around the 9 M ⊙ star µ2 Sco (Squicciarini et al. 2022) and a 11 M J planet candidate was detected in a 556 AU orbit around the 6-10 M ⊙ binary b Cen AB (Janson et al. 2021).
Given the importance of correctly determining the planetoccurrence rates for intermediate-mass stars and the difficulty in obtaining accurate masses for evolved stars, we began an RV survey around giant stars in open clusters.The advantage of open clusters is that the ages and masses of their stars can be much better constrained, meaning the planetary characterisation will be much more reliable.In the first paper of the survey, Lovis & Mayor (2007, hereafter Paper I) presented the discovery of a planet and a brown dwarf candidate orbiting a 2.3 M ⊙ red giant in the open cluster NGC2423 and a 3.8 M ⊙ red giant in NGC4349, respectively.In a subsequent work, Delgado Mena et al. (2018, hereafter Paper II) presented a planet candidate around a 2.1 M ⊙ giant in IC4651 that was found to show a suspicious signal in one activity indicator with a slightly shorter period than that of the RV variations.In addition, the analysis of new data obtained after the publication of Paper I pointed to a probable stellar origin of the RV signals found in NGC2423 and NGC4349.
These results made evident the complexity in analysing these noisy stars, which have long rotational periods that are compatible in many cases with the orbital periods of the candidate substellar companions.Interestingly, the large-amplitude RV sig- nals found in those three objects had periods of close to 700 days, as is also the true for other suspicious cases in the literature such as Aldebaran (Hatzes et al. 2015;Reichert et al. 2019) or γ Draconis (Hatzes et al. 2018).Indeed, the planetoccurrence rate from three combined large RV surveys analysed by Wolthoff et al. (2022) shows a maximum at 720 days, which might be caused by the accumulation of planets with orbital periods around 600 days found orbiting the more massive stars in the surveys (those with M >1.4 M ⊙ ).One of the explanations for this accumulation might be contamination by false positives, which only are starting to be revealed after years of observations.The aim of this work is to present new results of our RV survey for three stars that show periodic RV variations and may host substellar companions.The outline of the paper is as follows: in Sect. 2 we present the data and the stellar parameters.A planet candidate in NGC3680 whose signal might be of stellar origin is presented in Sect. 4. In Sects.5 and 6, we show the cases of two massive stars 2 in NGC2345 and NGC3532, which show largeamplitude RV signals that are likely caused by modulation of stellar magnetic activity.However, these stars also present secondary RV signals that might be caused by a brown dwarf or a planet.In Sect.7, we present additional data for three clusters stars with RV variations mimicking substellar bodies already discussed in Paper II.In Sect.8, we discuss the possible origin of all the signals presented in the previous sections.Finally, in Sect.9, we present our conclusions.
Table 1.Stellar characteristics of the analysed planet-host candidates (the first three stars are presented here for the first time).Stellar parameters (above the horizontal line) were derived by Tsantaki et al. (2023).Distances and magnitudes are extracted from WEBDA and SIMBAD, respectively.

Observations and stellar parameters
The RV survey on which this work is based is fully described in Paper I. The objects analysed in this survey were first observed over a time-period of nearly 5 years (from March 2005 to October 2009, ESO periods 75-83, PI: Lovis).In total, 142 stars were monitored within 17 open clusters using the HARPS spectrograph (Mayor et al. 2003) at the ESO 3.6m telescope (R ∼ 115000), and some of them were also observed with CORALIE (1.2 m-Swiss Telescope, La Silla Queloz et al. 2000;Udry et al. 2000) in previous years.For those stars showing large RV variations, we collected more observations between March 2017 and March 2022 (ESO periods 99-108, PI: Delgado Mena).In addition, we obtained a few RV points for a number of stars in the ESO archive (periods 91-94, PI: Alves, Canto-Martins).The RV values are provided in the online tables associated with this publication.The observations were made in objA mode (no simultaneous RV calibration) and the exposures times were estimated in order to have individual spectra with a signal-to-noise ratio (S/N) of at least ∼ 30 at ∼ 6000 Å .This gives a typical RV photon-noise of ∼ 3.5 m s −1 , which is sufficient to detect massive planets around the surveyed stars (with V magnitudes between 7 and 12).
We note that we applied a small negative offset to the RV points taken after May 2015, when an upgrade of the HARPS fibres took place (Lo Curto et al. 2015).The shift was calculated by extrapolating (with a linear fit) the measurements of Table 3 in Lo Curto et al. (2015) to the average full width at half maximum (FWHM) of the cross-correlation function (CCF) for each star (see Appendix A).This offset is larger for stars with larger FWHM.The offset values for the stars3 discussed in this paper are 22.3 m s −1 (NGC3680 No. 41) and 27.8 m s −1 (NGC3532 No. 670).No offset is required for NGC2345 No. 50 because we only have data taken after 2015.In addition, the FWHM needs to be corrected due to an instrumental focus drift affecting the measurements obtained before the fibre change in 2015.The equations for this correction are provided by Gomes da Silva et al. (2012).After the fibre change, the cause of the focus drift was corrected but both the FWHM and the bisector inverse slope (BIS, Queloz et al. 2001) of the CCF 4 present an offset with respect to previous values.As opposed to the RV offset, there is no published data with which to calculate the offset in the FWHM and BIS values.We refer the reader to Appendix A for more details about our attempts to correct this offset.
The stellar parameters of effective temperature (T eff ), surface gravity (log g), metallicity ([Fe/H]), and microturbulence (ξ t ) for most of our target stars in these clusters were presented by Santos et al. (2009Santos et al. ( , 2012) ) and were improved by Delgado Mena et al. (2016), who also derived stellar ages, masses, radii, and Li abundances.With the addition of more observations and some new targets, we presented a new set of more precise parameters for the complete sample of stars (Tsantaki et al. 2023).In this latter analysis, we derived stellar parameters with the spectral synthesis method and made use of MARCS model atmospheres, which are more appropriate for giant stars.In addition, the masses and radii of the stars were revised using Gaia DR2 parallaxes.We refer the reader to Tsantaki et al. (2023) for further information.

Radial velocity, stellar activity, and photometry analysis
The analysis of the RV data was carried out with kima (Faria et al. 2018).This code facilitates a Keplerian fit where the number of planets is a free parameter to be estimated from the data, using the Diffusive Nested Sampling algorithm (Brewer et al. E. Delgado Mena et al.: Planets around evolved intermediate-mass stars 2010) to sample from the posterior distribution of the model parameter.kima also allows the inclusion of Gaussian processes to model stellar activity, although we did not consider this component for the analysis of the stars presented here (as in many cases the signals show coherence).In the present work, we use kima by inputting wide and uninformative priors for the orbital parameters, that is, the orbital period (log-uniform prior from 0.2 to 2000 days), the semi-amplitudes (modified log-uniform prior from 1 m/s to 1 km/s), and the eccentricities (uniform prior between 0 and 1).These priors, and indeed the whole analysis with kima, are not informed by the periodogram of the RVs or the activity indicators.
As a double check, we performed a second analysis of each dataset with the yorbit algorithm (Ségransan et al. 2011) in order to fit the whole dataset with a model composed of a Keplerian function.Yorbit uses a hybrid method based on a fast linear algorithm (Levenberg-Marquardt) and genetic operators (breeding, mutations, crossover), and is optimised to explore the parameter space for Keplerian fitting of RV datasets.In our case, the global search for the best-fitting orbital parameters was made with a genetic algorithm.As we show in the following subsections, the photon noise RV errors of our observations are typically below 3 m s −1 .However, the RV variations are dominated by the stellar jitter, which is much higher than the photon noise in this kind of star (typically of ± 15-20 m s −1 , as discussed in Hekker et al. 2008;Delgado Mena et al. 2018;Lovis & Mayor 2007).Therefore, when using Yorbit, we added in quadrature a 15 m s −1 noise to the photon noise before fitting the data.In the case of kima, the RV jitter is a free parameter (with a uniform prior) that is already considered during the fitting process.
To gain insight into possible activity modulations that could interfere with RV, we used ACTIN5 (Gomes da Silva et al. 2018) to obtain the activity indices from the spectra based on the Hα, Na i D 1 and D 2 , and He i D 3 lines.The line parameters for the Na i and He i indices are described in Gomes da Silva et al. (2011).For the Hα index, we used two versions with different bandpasses: Hα16, using a 1.6 square passband (as used in our previous work) and Hα06 using a 0.6 passband (optimal for FGK dwarfs; see Gomes da Silva et al. 2022).A description of the ACTIN flux and indices determination is available at Gomes da Silva et al. (2021, Appendix A).In general, we could not use the Ca ii H&K for any of the stars due to the very low S/N (< 3) in the spectral orders containing those lines.In addition, the Na i and He i indicators do not show any signal for our stars and therefore we do not show them in the figures.
To further investigate RV variations caused by stellar atmospheric phenomena, we also used the FWHM.These values and their errors are provided by the HARPS pipeline.Moreover, we also analysed the BIS (Queloz et al. 2001), but we note that this diagnostic of line asymmetry loses sensitivity for slow rotators (as some of the stars we study here) (e.g.Saar et al. 1998;Santos et al. 2003;Queloz et al. 2009;Santos et al. 2014).Nevertheless, we have already observed significant correlations between RV and BIS (e.g.see the case of NGC2423 No. 3 in Paper II).We note that we only obtain the periodograms for activity indicators of HARPS data because the number of observations per star with CORALIE is too low for meaningful results.
Most of the stars in our survey were observed within the All Sky Automated Survey (ASAS) from the las Campanas Observatory (Chile) with observations available between the end of 2000 and 2009 (Pojmanski 2002).We downloaded the light curves in V magnitude from The ASAS-3 Photometric V-band

Parent-star characteristics
Our sample contains five giant stars in the open cluster NGC3680 (Age = 1.78 Ga) with an average metallicity of [Fe/H] = -0.15± 0.02 dex (see Tables 1 and 5 of Tsantaki et al. 2023).One of these stars, NGC3680 No. 41 with T eff = 4612 ± 12 K, log g = 2.45 ± 0.04 dex, M * = 1.64 ± 0.05 M ⊙ , and R * = 11.44 ± 0.47 R ⊙ (see Table 1) seems to be on the first ascent of the RGB, and close to the luminosity bump (see Fig. 1).The mean RV of the giants in this cluster is 1.23 ± 0.65 km s −1 while the mean RV of NGC3680 No. 41 is 1.59 km s −1 , and therefore this star is likely a cluster member.This is also supported by its parallax value.In order to obtain an estimation of the maximum rotational period of this star (which is relevant for the interpretation of the observed signals), we can use the projected rotational velocity (v sin i) and the radius.The v sin i of the sample stars were estimated by Tsantaki et al. (2023) using spectral synthesis and by considering a fixed macroturbulence velocity (with the empirical formula for giants by Hekker & Meléndez 2007), which in this case has a value of 5.04 km s −1 .For this star, the v sin i is very low and below the detectability threshold (i.e. the value of the macroturbulence velocity is large and can perfectly account for the broadening of the spectra).Therefore, we cannot provide a reliable estimation of the stellar rotational period but we can safely claim that this star is a slow rotator.

Radial-velocity analysis
For NGC3680 No. 41, we have 18 RV measurements obtained with HARPS over a time period of 9 years (between 2013 and 2022) and 5 measurements with CORALIE obtained over one year (2003)(2004), therefore with a large time gap between the two data series.The average photon noise RV uncertainty is 2.6 and 2.1 m s −1 for HARPS and CORALIE data, respectively.
An analysis using generalised Lomb-Scargle (GLS) periodograms (e.g.Zechmeister & Kürster 2009) was performed for RV (see upper panel of Fig. 2).For each star, the minimum period to search in the periodogram was defined as the minimum cadence of observations when these are done with a separation of more than 10 days.Otherwise, the minimum period is set at 10 days.The false-alarm probability (FAP) was computed by bootstrapping the data and statistically significant peaks were considered for values above FAP = 1%.A significant signal with a period of ∼ 1127 days can be clearly observed in the periodogram of HARPS data only.Using kima, the two RV sets are well fitted by a Keplerian function with P = 1155 days, K = 74.79m s −1 , and e = 0.21 (see Table 5) with a 16.35 m s −1 dispersion of the residuals.The posterior distribution for the fitted extra white noise is centred at 20±5 m s −1 , which is in agreement with the added jitter when using yorbit (as explained in Sect.3).Considering the mass of the star, these values correspond to the expected signal induced by a planet with 5.13 M J and a 2.53 AU semi-major axis.The phase curve of the best-fit solution can be seen in Fig. 3.The results using yorbit are very similar and the data can be best fit with a single Keplerian with a period of 1154 days, K = 74.7 m s −1 , and e = 0.19 (figures from yorbit can be seen in the Appendix, Fig. B.1) The dispersion of the residuals is 16 m s −1 and the reduced χ 2 is 1.37.In the following subsections, we assess whether or not this signal could be of stellar origin by evaluating different stellar activity or variability indicators.

Photometry
We found 664 photometric measurements classified as good quality (these are given the grade A or B, with average errors of 0.045 mag), which show a ∼ 0.3 peak-to-peak variability in V magnitude.In Fig. 4, we can see that the GLS periodogram shows no signals above the FAP level and therefore the RV variations are likely not linked to photometric variability.

Stellar activity and line profile analysis
The GLS periodograms of the FWHM and the BIS do not show any significant period (see Fig. 2).However, the RV variations from HARPS seem to be correlated with the BIS (a Pearson's coefficient value of 0.49 and Spearman rank coefficient of 0.54; see Table .3 and Fig. 5), although these correlations are not significant (with a p-value above 0.01).On the other hand, for the CORALIE data, the RV-BIS correlation is much stronger (Pearson's coefficient value of 0.85, though this is based on five points and is not significant, with a p-value of 0.08).In addition, the analysis of the Hα06 index does not show any periodic behaviour but the GLS of the Hα16 index has a peak at the 1% FAP level for a period of 312 days.As we cannot tightly constrain the rotational period of the star, it is difficult to confidently discern whether or not this activity signal is a manifestation of the rotation of the star.As a further test, we detrended the RV points based on a linear fit to the RV versus BIS correlation.When performing a periodogram of the detrended RVs, the peak at 1127 days still shows a strong power of close to 0.8, but below the 1% FAP line.Although the activity indicators do not exhibit any periodic behaviour with the same period as the RV, the possible RV-BIS correlation mentioned above casts doubt on the planetary hypothesis and is further discussed in Sect.8. Furthermore, we note that the correlation between the BIS and the RV residuals of the Keplerian fit becomes weaker (with a Pearson's coefficient value of 0.06), which could point to a stellar origin of the variation.1), is the most evolved star in the cluster and seems to be ascending the early asymptotic giant branch (AGB; see Fig. 1).This star is very interesting because it is the youngest and most massive in the full sample.Despite the fact that the other stars in the cluster have very few observations, we are able determine the mean RV of the cluster as 58.61 ± 0.56 km s −1 while the mean RV of NGC2345 No. 50 is 59.21 km s −1 .Although the scatter is large due to the large RV variations observed in the star, its RV is within ∼ 1σ of the RV dispersion in the cluster, and taking into consideration its parallax, this star seems to be a cluster member.We note that recent studies in this cluster by Holanda et al. (2019) and Alonso-Santiago et al. ( 2019) also catalogue this star as a member and obtain similar stellar parameters.

Large RV variations in
Both mentioned studies place this star in the RGB, as also suggested by its not very low 12 C/ 13 C ratio.We estimated a v sin i of 5.27 km s −1 in Tsantaki et al. ( 2023) for a fixed macroturbulence velocity of 6.06 km s −1 (with the empirical formula for bright giants by Hekker & Meléndez 2007).Using the stellar radius, we can derive the maximum rotational period via P max = 2πR * /vsini.This leads to a maximum rotational period of ∼ 1463 days.

Radial velocity analysis
A total of 24 RV measurements were obtained with HARPS over ∼5 years (between 2017 and 2022).The average photon noise error in RV is 3.3 m s −1 .The GLS periodogram was performed for RV (see upper panel of Fig. 6), which shows a strong signal (well above the 0.1% FAP level) at a period of 1007 days.The analysis with kima reveals that the maximum likelihood is obtained for a model with one Keplerian with a period of 1001 days, K = 488.17m s −1 , and e = 0.03.Considering the mass of the star, these values would correspond to the expected signal induced by a body with 3.5 AU semi-major axis and with 79.34 M J , a value that places this companion at the threshold between brown dwarf and low-mass star.The phase curve of the best-fit solution can be seen in Fig. 7.The posterior distribution for the fitted extra white noise is centred at 55±18 m s −1 .This value is relatively high but is to be expected in very evolved stars, as is the case for this star (Hekker et al. 2008).The residuals of the fit have a dispersion of 60.4 m s −1 , which appears large and might be indicative of an additional signal.
Indeed, the analysis performed by yorbit with a single Keplerian provides very similar output parameters (P=1001 days, K = 488.3m s −1 and e = 0.036), including the large dispersion of the residuals and a large reduced χ 2 of 22.5 (see Fig. B.2).We therefore attempted to make a two-Keplerian fit with yorbit to model both the large-amplitude signal (probably caused by activity, as we show below) and a possible existing substellar body in the system.We obtained a solution for a system with a low-mass star and a brown dwarf, with the period of the 'lowmass star' being 1007 days (close to the period of the single Keplerian fit explained above) and with a slightly larger semiamplitude (K = 554 m s −1 ) corresponding to a mass of 87.36 M J .The brown-dwarf candidate would have a period of 561 days, a mass of 23.27 M J , and an eccentricity of 0.28 (see Fig. C.1).In this case, the dispersion of the residuals is lower, 32 m s −1 , and the reduced χ 2 is 8.1; this model therefore better represents the observed data than the single Keplerian fit.We also performed a GLS periodogram of the residuals, which did not show any significant peak.A complete analysis of activity and photometry signals is presented in the following subsection, with the aim being to better understand if the brown-dwarf candidate is bona fide.

Photometry
We found 501 photometric measurements in the ASAS-3 catalogue classified as good quality (they are given a grade A or B, with average errors of 0.049 mag) that show a ∼ 0.2 peak-to-peak variability in V magnitude (not considering the observations in the first epoch showing a large scatter).In Fig. 8, the GLS periodogram for this dataset shows several signals well above the FAP = 0.1% line with periods of 3000, 1239, 329, 422, and 286 days.Considering the errors associated with the determination of stellar radius and especially with the determination of v sin i (due to the degeneracy with macroturbulence velocity), the signal at 1239 days is compatible with the maximum rotational period estimation of 1463 days.On the other hand, none of these photometric signals are close in period to the large RV variation observed at 1007 days, nor with the signal at 561 days of possible brown dwarf origin.

Stellar activity and line-profile analysis
When studying the RV signals of this star, we first noted a very strong correlation with the FWHM (with a Pearson's coefficient value of -0.72 and a p-value of 7•10 −5 , also strong and significant with the Spearman statistics; see Table 3), which led us to suspect the large RV variability observed.On the other hand, the correlation with the BIS is not significant.We also performed the GLS periodograms for the different activity indicators (see Fig. 6).In this case, the FWHM shows a non-significant peak at 957 days and the BIS has a non-significant peak at 1409 days.It is possible that this signal is related to the 1239 d signal observed in photometry.The most significant signals are seen in the Hα line.The Hα06 index shows a strong signal above the 1% FAP level at 962 days, close to the 1007 d period of the RV.Furthermore, there are two signals at 415 and 586 days, though these are not significant.The situation is much more clear when using the Hα16 index, which shows a strong signal above the 1% FAP level at 1113 days and another signal above the 1% Table 4. Summary of the strongest periods found in the different periodograms.First set of columns: Maximum stellar rotational period and periods of the most significant peak(s) in each GLS periodogram that match the most significant RV period.We highlight in blue the periods that are compatible with the estimated stellar rotation period.Second set of columns: For the stars with a possible two-Keplerian fit, we list the periods of the periodogram peaks that match the period of the hypothetical second companion in the system.The periodogram peaks that lie above the 1% FAP level are marked in boldface.4 for a summary of the different peaks in the periodograms.
In addition, we checked whether the Gaia data reveal the presence of a low-mass star companion through elevated astrometric noise.We consulted the renormalised unit weight error (RUWE 7 ) statistic, which has a value of 1.05 for this star (well below 1.4) and therefore does not show evidence of a companion.This fact further supports the conclusion that the large-7 https://gea.esac.esa.int/archive/documentation/GDR2/Gaia_archive/chap_datamodel/sec_dm_main_tables/ ssec_dm_ruwe.html  1 and 5 of Tsantaki et al. 2023).One of its members, NGC3532 No. 670, with T eff = 4347 ± 11 K, log g = 1.75 ± 0.05 dex, M * = 3.05 ± 0.23 M ⊙ , and R * = 40.95± 2.37 R ⊙ (see Table 1), seems to be at the tip of the RGB but it is not clear whether it is approaching the tip of the RGB (see Fig. 1).The mean RV of the giants in this cluster is 4.03 ± 0.48 km s −1 while the mean RV of NGC3532 No. 670 is 4.20 km s −1 .This similar value, together with its parallax, suggests that this star is likely a cluster member.We estimated a v sin i of 4.64 kms in Tsantaki et al. (2023) for a fixed macroturbulence velocity of 4.52 kms (with the empirical formula by Hekker & Meléndez 2007, for giant stars).This leads to a maximum rotational period of ∼ 447 days.

Radial velocity analysis
For this star, we collected a total of 29 RV measurements with HARPS over ∼9 years (between 2013 and 2022) although there is a gap of 2.5 years without observations.The average photon noise error in RV is 1.3 m s −1 , which is due to the fact that this star is considerably brighter than the previous two described in Sections 4 and 5, respectively.The time series of the RV show a clear variation of large amplitude.This is also demonstrated by the GLS periodogram of the RV (see upper panel of Fig. 9), which presents a strong signal (well above the 1% FAP level) at a period of 842 days.The most probable solution obtained with kima is that with a single brown dwarf of 22.7 M J with a period of 842.4 days and an eccentricity of 0.12 (K = 235.3m s −1 ;

Photometry
We found 549 photometric measurements classified as good quality (they are given a grade A or B, with average errors of 0.046 mag) that show a ∼ 0.3 peak-to-peak variability in V magnitude (excluding the probable outlier points above V = 7.2).The GLS periodogram depicted in Fig. 11 shows four signals above the FAP = 0.1% line with periods of 3333, 412, 1363, and 503 days.None of these periods match the signal at 842 days found in the RV time series or the planet candidate signal in the two-Keplerian fit.However, the period at 412 days is compatible with the estimated rotational period of the star and this would suggest that the RV signal at 842 days is a harmonic of the rotational period.

Stellar activity and line profile analysis
For this star, we observed that the RV variations are not strongly correlated with the BIS, nor with the FWHM (see Table 3   the Spearman and Pearson coefficient and p-values), providing tentative evidence that the large-amplitude RV signal might be of true brown dwarf or planetary origin.However, as we also observed for similar stars in our previous work (e.g.NGC4349 No. 127, Paper II), a lack of correlation can be produced if these activity indicators are not in phase with the RV.This can be assessed with the GLS periodogram of the different activity indicators (shown in Fig. 9).The FWHM shows a signal at 811 days that is probably related to the RV signal at 842 days and a second signal at 619 days also above the 1% FAP level.Therefore, the lower-amplitude signal for the two-Keplerian orbit seems to also be of stellar origin.Indeed, the window function has a significant peak at 363 days that is introducing yearly aliases in the periodograms.The yearly aliases of the 842 d signal fall at 254 and 641 days, whereas the yearly aliases of 625 days fall at 230 and 872 days.Therefore, these two signals are yearly aliases of each other.On the other hand, the BIS does not show any significant period.As was found for NGC2345 No. 50, the strongest signals are observed for the Hα16 indicator, which has a peak at 889 days that is well above the 1% FAP level.Curiously, this star shows no significant signal in the Hα06 indicator.

Revision of previous planetary candidates
In this section, we present a summary of the updated results with the new data collected for the targets of Paper II, where we presented the cases of three stars that show large-amplitude and long-term signals disguised as substellar companions.

IC4651No9122
We have collected 15 additional RV points since 2018 (for a total of 74).In paper II, this star is shown to have a strong RV signal at 747 days but with signals in FWHM and Hα16 at 714 days (above or close to the FAP level), which led us to refrain from claiming the presence of a planet around this star.However, in the GLS of the updated time series (see Fig. 12), none of the activity indicator signals has a significant peak at the same period as the most significant peak in the RV at 744 days.We also checked that the RV signal at 744 days is not an alias of other signals.The window function shows four significant peaks at 30, 365, 1785, and 4466 days.The three longer-period signals introduce significant aliases of the 744 d signal in the RV periodogram at periods of 245, 1275, and 893 days, respectively.The RV analysis with kima provides a two-Keplerian fit as the solution with the largest likelihood.The time series and orbital phase of the two planet candidates with periods of 743 days (K = 99.6 m s −1 and e = 0.13) and 384 days (K = 27.7 m s −1 and e = 0.25) can be seen in Fig. D.1.However, the period for the possible second planet in the system is very close to two significant peaks in the FWHM periodogram at 364 and 393 days.This case again highlights the importance of the FWHM as a reliable stellar indicator.Nevertheless, we note that these signals in FWHM could simply be caused by the window function, which has a peak at 366 days, and is not related to activity.On the other hand, if those are true stellar signals, then it is worth noting that they might be the first harmonic of the signal at 744 days (P/2 = 372 days), which would cause us to doubt the hypothesis that IC4651No9122 is a true planet candidate.With the data in hand, we can at least be sure that the signal fitted with the second keplerian at 384 days is either of stellar origin (reflected in the FWHM if real) or is caused by the window function or a harmonic of the 744 days signal.Therefore, we consider only the signal at 744 days as a planet candidate and repeated the fit with kima, fixing the number of planets to just one (see Fig. 13).In this case, as expected, the residuals of the fit have a larger dispersion (20.39 m s −1 ) and the posterior distribution for the fitted extra white noise is also centred at a larger value (20.5±1.8 m s −1 ).The final planet candidate has a period of 746 days (K = 94.98 m s −1 and e = 0.04) and a minimum mass of 6.22M J (see Table 5 for the detailed parameters).
We note that, in the present work, we are also using the updated stellar parameters from Tsantaki et al. (2023), and now the estimated v sin i is 0.68 kms for a fixed macroturbulence velocity of 4.98 kms (with the empirical formula by Hekker & Meléndez 2007, for giant stars).This leads to a maximum rotational period of ∼ 995 days.This period is close to the period of 985 days seen in the Hα16 periodogram; however, it is not significant and prevents us from relating them.Therefore, we cannot completely rule out the possibility that the 744 d RV signal is related to the star's rotation (if this were shorter in reality).The lack of similar periods in the activity indicator periodograms (as well as the lack of strong correlations; see Table 3) suggests that this is a bona fide planet but the proximity of the two significant FWHM peaks to the first harmonic of the strong RV signal warns us of the possibility of a non-planetary signal.We note that if the recent FWHM values were corrected to take into account the observed offset (see a more detailed discussion in Appendix A), the GLS periodogram of the FWHM would show a significant peak at 687 days.This period seems safely far from the RV period at 744 days.Additionally, the peaks seen at 364 and 393 d would go below the 1% FAP level, which would support the planet-candidate hypothesis.Nevertheless, we can see with this example how the usefulness of the FWHM indicator can be limited if corrections need to be applied.Therefore, additional data will be needed to further assess the nature of the signals observed here.

NGC2423No3
We collected a further 16 RV points for this star (raising the total to 92, with 32 of them from CORALIE).In our previous work, this star had a strong RV signal at 698 days but with a strong RV-BIS correlation and a peak in the BIS periodogram close to the FAP level at a similar period.The GLS of the complete dataset is shown in Fig. 14.The RV periodogram presents the most significant peak at 697 days (with two aliases at 828 and 614 days caused by the long gap in time between the two sets of observations, specifically by the peak at 4506 days in the window function periodogram).In addition, the window function shows several significant signals around one year, which in turn produce aliases around 230 days in the RV periodogram.As opposed to previous findings, with the new data the BIS periodogram shows a significant peak at 697 days.Moreover, the RV-BIS correlation has a significant Pearson coefficient as well as Spearman rank statistics (see Table 3).With the new data, both kima and yorbit would fit the RV time series with a two-Keplerian orbit caused by planetary bodies at 700 and 435 days (with K ∼ 130 and 30 m s −1 , respectively, see Fig D .2).However, the long-period signal appears to be of stellar origin given the correlation with the BIS, and the second signal has a period similar to one of the peaks seen in the periodogram of the photometric data (period of 417 days, see Fig. 8 in Paper II).In addition, both the FWHM and Hα16 periodograms show significant peaks at 380 and 393 days, respectively.
We note that with the updated stellar parameters from Tsantaki et al. (2023), we estimated a v sin i of 2.19 kms (for a fixed macroturbulence velocity of 4.88 kms), which leads to a maximum rotational period of ∼ 407 days.These findings lead us to suspect that all these signals around 400 days are caused by rotational modulation of active regions, which in turn also produces the 435 d RV signal fitted in the second Keplerian orbit.

NGC4349No127
For this target, we previously worked with 46 HARPS measurements up to 2009 (as presented in our previous publication) and now we have 11 additional points (between 2018-2022).This star previously showed a large RV signal at 672 days but with a strong and clear signal in FWHM at 666 days (and also the periodogram of Hα16 had a close to significant peak at 689 days., see Paper II).Those findings indicated that the brown dwarf presented in Paper I was not real.We collected additional data to evaluate the stability of the RV and FWHM signals (as some examples in the literature have shown changes in phase for longterm RV signals; e.g.Aldebaran Reichert et al. 2019) and to learn more about its origin.With the new data, the RV time series can also be fitted with a single Keplerian with the same period, of namely 674 days (K ∼ 226 m s −1 ; see Fig. D.3).However, this long-term signal is also observed in the FWHM periodogram (see Fig. 15), as was previously found to be the case, but not in the Hα periodograms.However, we note that we find a weak but significant correlation between RV and Hα16 (see Table 3).Interestingly, we do not observe a change of phase in the RV (see the orbital phase in Fig. D.3), but if we fit the FWHM with a sinusoidal with a period of 675 days (see Fig. D.4), there seems to be a change in the phase.In addition, the window function shows a strong peak at 4719 days caused by the long gap in between the observations, which in turn produces significant aliases of 588 days (seen in the RV and FWHM periodograms) and 785 days, which can be seen in the FWHM periodogram.
For this star, we also have updated stellar parameters from Tsantaki et al. (2023), for which we estimated a v sin i of 4.81 kms (for a fixed macroturbulence velocity of 4.66 kms), which leads to a maximum rotational period of ∼ 399 days.The periodogram of the photometric data shows two significant peaks at 428 and 342 days, which may be related to the stellar rotation (see Fig. 12 in Paper II).As discussed in this latter work, we cannot discard that the estimated rotational period (<399 days) and the 342 d signal in the photometry are the first harmonics (P/2) of the RV period.Taken together, these findings point to a non-planetary origin of the signal.

The importance of each activity indicator
The detailed analysis presented in this work (and our previous study) demonstrates the difficulty in confirming planetary signals in evolved intermediate-mass stars.Our most important finding is that there is not a single (e.g.universal) activity indicator that we can trust to unveil a spurious planetary signal; all of them have to be analysed in order to understand the origin of the different RV signals, and in some cases this might not be sufficient.We find hints that the CCF-BIS might be more useful for less massive stars (around 2 M ⊙ ), because we find a possible RV-BIS correlation for NGC3680 No. 41 and a clear signal in the BIS periodogram (as well as a significant strong RV-BIS correlation) for NGC2423 No. 3.
On the other hand, the CCF-FWHM and Hα lines seem to be more useful for tracking activity signals in more massive and evolved stars (NGC2345 No. 50, NGC3532 No. 670, and NGC4349 No. 127); see Table 4 for a summary of the periodogram results.However, we notice that conclusions on the origin of a given signal are highly dependent on the way the Hα line is measured.Meanwhile NGC2345 No. 50 has significant peaks at both the Hα16 and Hα06 indicators, NGC3532 No. 670 only shows a significant peak for Hα16, and NGC4349 No.127 shows no significant peak for any of them.Nevertheless, although not significant, the signal of Hα16 for NGC4349 No. 127 is much stronger than that of Hα06.
We acknowledge that the sample analysed here is still not sufficiently large for us to reach clear conclusions, but there are indications that Hα16 is a more reliable indicator than Hα06 for these evolved stars (the opposite is true for FGK dwarfs, for which Hα06 is more appropriate; Gomes da Silva et al. 2022).A thorough analysis of stellar activity indicators in evolved stars is outside the scope of the present paper but will be explored in a future work.On the contrary, measuring the Hα line with a very wide passband can result in spurious periods in the periodogram, because there are contaminating lines in the wings of this strong line.This is the explanation for the recent confirmation of a brown dwarf around the largest star to date (HD 18184; Lee et al. 2023), whose RV variations were previously deemed to be of stellar origin given the similar period in the Hα periodogram (Bang et al. 2018).However, the work by Lee et al. (2023) warns that a ± 1Å passband around the line centre (and not a ± 2Å as used in Bang et al. 2018) is needed to avoid the blending lines.We note that our Hα16 index is measured with a ± 0.8 Å passband around the line centre and is therefore not affected by nearby blending lines (see Figs. E.1 and E.2 in the Appendix).
Finally, we believe that the most critical indicator for this kind of star is the CCF-FWHM, because it has proven useful for detecting the probable stellar origin of the largest-amplitude signals (in the three most massive and evolved stars) but also that of the weaker RV signals (that can be fitted with a second Keplerian) found in stars with a broad range of stellar masses (the cases of NGC3532 No. 670, IC4651 No. 9122, and NGC2423 No. 3).We stress the fact that for NGC4349 No. 127, despite having nearly identical parameters to NGC3532 No. 670, only the FWHM has served to cast doubt on the planetary origin of the RV variability.Finally, the photometric data have been useful in some cases as well, but the long rotational periods of our targets make it difficult to find data of sufficient quantity to track the real period of the stars.(Wolthoff et al. 2022).Nevertheless, the targets presented here are a fraction of the complete sample and a more in depth analysis on planet frequencies will be presented in a future work along with results for the full sample.
In recent years, an increasing number of suspected planet candidates have been found around evolved stars with orbital periods in the range of ∼ 600-800 days.This is the case for NGC4349 No. 127, NGC2423 No. 3, NGC3532 No. 670 discussed here, Pollux (β Gem) (Hatzes et al. 2006;Aurière et al. 2021), Aldebaran (Hatzes et al. 2015;Reichert et al. 2019), and γ Draconis (Hatzes et al. 2018).In addition, the work by Tala Pinto et al. (2020) presents two evolved stars, 3CnC and 44UMa, with periodic RV variations of ∼ 800 days that do not show variations with a similar period in photometry or Hα index.However, given the similarity in stellar parameters and luminosities with the above-mentioned cases and the lack of line-shape-stability indicators such as FWHM (not available in non-stabilised spectrographs, as used in their work), Tala Pinto et al. (2020) classify these objects as planet candidates.The luminous red giant HD 81817 (4.3 M ⊙ , 83.8 R ⊙ ) also shows a secondary RV signal with a period of ∼ 630 days (Bang et al. 2020) but the Hα line shows a significant peak at a similar period as well.
The accumulation of planet candidates in that period range (and more specifically between 300 and 800 days) has also been noted by Döllinger & Hartmann (2021), who in addition points out that this only happens for stars with a radius of greater than 21 R ⊙ (see their Fig. 1).These authors also warn of the possibility that the observed periodic signals are of stellar origin and propose that a large stellar radius and a low stellar metallicity should be seen as a warning signal that RV signals may not be of planetary nature.It is worth noting that the three objects with the clearest RV variations in our sample have large radii (they are the most evolved in each cluster) and metallicities below -0.1 dex.
As in many cases the rotational periods of these large, slow stars are compatible with the RV periods, the simplest explanation for the RV variability would be the rotational modulation of active regions, such as spots.However, long-lived stable spots in evolved single stars have only been reported in the literature for very active stars hosting strong magnetic fields, those that are believed to be the descendants of magnetic Ap-Bp stars.This is the case for EK Eri (Aurière et al. 2011) or β Ceti (Tsvetkova et al. 2013) with P rot of 309 and 215 days, respectively, which also present large photometric variations.
However, the stars presented in this work seem to belong to another group of red giants; these show longer rotational periods, weak large-scale magnetic fields (sub-Gauss level), and higherthan-expected chromospheric activity (Aurière et al. 2015).The best studied case in the literature is Pollux, which shows RV variations with a period of 590 days -that is stable for 25 years-and no apparent S-index variability (Hatzes et al. 2006).Nevertheless, the planet candidate around Pollux was challenged by the discovery of magnetic field variations with a similar period to the RV variations (Aurière et al. 2014).In a more recent work with additional data, Aurière et al. (2021) found a longer period (660 days) for the weak magnetic field variations, but this could still compatible with the RV period due to uncertainties.Therefore, doubt remains as to the true existence of a planet orbiting Pollux, as is also true for the cases presented in the present work.In the cases of red giants with weak magnetic fields such as Pollux, if the detected RV variations are due to the magnetic field, they would not be caused by photometric spots but by magnetic plages or other magnetic structures locally reducing convection (Aurière, priv. comm).Therefore, the lack of photometric variation in red giants (and thus of spots) is not a guarantee that the RVs are of planetary origin.
We do not have sufficient S/N for our sample of stars in the blue region of the spectra to be able to study the S-index variability in order to compare it with the behaviour of Pollux and other giants with weak magnetic fields like Aldebaran.However, many of our stars show Hα variations, which is also a sign of chromo-E.Delgado Mena et al.: Planets around evolved intermediate-mass stars spheric activity.It would be very interesting to further study the magnetic activity of these stars, although it is very challenging to detect weak magnetic fields for relatively dull stars (V ∼ 10) with the currently available instrumentation.
An alternative explanation for the large RV variations in evolved stars is the presence of a poorly understood class of oscillations called oscillatory convective modes (Saio et al. 2015).This class of oscillations is theoretically possible for 2 M ⊙ stars with high luminosity, that is, with log (L/L ⊙ ) ≳ 3 dex; see more detailed discussion in Paper II.The clearest candidate in our sample would therefore be NGC2345 No. 508 , with log (L/L ⊙ ) ∼ 3.9, though NGC4349 No. 127 and NGC3532 No. 670 also have high luminosities, as do other candidates in the literature (see Fig. 15 in Tala Pinto et al. 2020).

Conclusions
In this work, we present some of the results of a long-term survey of RVs in a sample of more than 140 giant stars in 17 open clusters.Data covering such a large time span are crucial in order to find long-period planets, whereas the probability of discovering short-period planets around giant stars diminishes as the stars evolve and their radii expand.The long-term monitoring of giant stars is also essential for distinguishing RV variations caused by the rotational modulation of active regions, because the rotational periods of evolved stars are much longer than those of dwarfs (from a few hundred to more than one thousand days).From the six stars analysed in this work, we are only confident in the planetary origin of the RV signals coming from two of those systems, and the star in both of them has a mass of below 2 M ⊙ .Nevertheless, we believe that additional data are needed to fully confirm the presence of planets around those stars.On the other hand, the more massive stars in the sample show periodic variations similar to the RV in one or several stellar activity indicators (Hα, BIS and FWHM), with the FWHM of the CCF being the most critical in our opinion.We also note that in the cases where a two-Keplerian fit was a good solution to model the RV variability, the Hα16 index and especially the FWHM were decisive in ruling out the presence of a smaller body in these systems.
The cases presented here and others in the literature demonstrate the difficulty in finding even large planets around intermediate-mass evolved stars.It is therefore of utmost importance to obtain a comprehensive picture of the different phenomena causing RV variability in evolved stars, and especially of their magnetic fields and their manifestations on different timescales.Once we understand such phenomena, we will be in a much better position to detect giant planets but also lowermass planets currently hidden in the stellar jitter.Near-infrared (NIR) spectra can be combined with optical spectra in order to discern between stellar activity and planetary signals, as the nonplanetary RV signals are wavelength-dependent, whilst planetary signals are achromatic (e.g.Figueira et al. 2010;Trifonov et al. 2015;Carleo et al. 2020;Carmona et al. 2023).In addition, the study of other stellar-activity indicators in the NIR will help us to better understand the magnetic activity of these stars and its relation with their RV variability.Therefore, as a further step towards a complete understanding of the variability in red giants, we will observe these stars with the new high-resolution NIR spectrograph NIRPS (Bouchy et al. 2017;Wildi et al. 2017) now operational on the 3.6m ESO telescope in La Silla (Chile).

E
. Delgado Mena et al.: Planets around evolved intermediate-mass stars

Fig. 1 .
Fig. 1.Hertzsprung-Russell diagram for the six open clusters analysed in this work.The upper row contains the older clusters (with less massive stars, as shown in the colour scale) and the bottom row contains the younger clusters, with more massive stars.The objects with periodic RV variations are depicted with a star symbol.

Fig. 2 .
Fig. 2. GLS periodogram of the RV, FWHM, BIS, and stellar activity indexes for NGC3680 No. 41 HARPS data.The blue dashed line marks the period of the planet candidate.The period of the strongest peak in each periodogram is marked with a red solid line.The horizontal line marks the 1% FAP level.The last panel shows the periodogram for the window function.

Fig. 3 .
Fig. 3.One Keplerian fit for NGC3680No41 using kima.The top panel shows the phase curve for the signal, with the residuals shown below, highlighting the rms of the residual RVs.The bottom panel shows the RV data from HARPS (red points) and CORALIE (blue points) and representative samples from the posterior distribution.The systemic RV has been subtracted and is shown at the top.
NGC2345 No. 50 mimicking the presence of massive companions 5.1.Parent star characteristics Our sample contains only four giant stars in the open cluster NGC2345 (Age = 0.07 Ga) with an average metallic-

Fig. 4 .
Fig. 4. GLS of V magnitude for NGC3680 No. 41 using ASAS-3 data.The horizontal lines indicate the FAP at levels of 0.1%, 0.5%, and 1%.The vertical red lines show the periods of the three signals with the highest significance.

Fig. 6 .
Fig. 6.GLS periodogram of the RV, FWHM, BIS, and stellar activity indexes for NGC2345 No. 50 HARPS data.The blue dashed line marks the period of the planet candidate.The period of the strongest peak in each periodogram is marked with a red solid line.The horizontal line marks the 1% FAP level.The last panel shows the periodogram for the window function.

Fig. 7 .
Fig. 7.One Keplerian fit for NGC2345No50 using kima.The top panel shows the phase curve for the signals, with the residuals shown below, highlighting the rms of the residual RVs.The bottom panel shows the RV data from HARPS and representative samples from the posterior distribution.The systemic RV has been subtracted and is shown at the top.

Fig. 8 .
Fig. 8. GLS periodogram of V magnitude for NGC2345 No. 50 using ASAS-3 data.The horizontal lines indicate the FAP at the levels of 0.1%, 0.5%, and 1%.The vertical red lines show the periods of the five signals with the highest significance.The time series are fitted by a sinusoidal with the period showing the highest significance.

Fig. 9 .
Fig. 9. GLS periodogram of the RV, FWHM, BIS, and stellar activity indexes for NGC3532 No. 670 HARPS data.The blue dashed line marks the period of the planet candidate.The period of the strongest peak in each periodogram is marked with a solid red line.The horizontal line marks the 1% FAP level.The last panel shows the periodogram for the window function.

Fig. 10 .
Fig. 10.One Keplerian fit for NGC3532 No. 670 using kima.The top panel shows the phase curve for the signals, with the residuals shown below, highlighting the rms of the residual RVs.The bottom panel shows the RV data from HARPS and representative samples from the posterior distribution.The systemic RV has been subtracted and is shown at the top.

Fig. 11 .
Fig. 11.GLS periodogram of V magnitude for NGC3532 No. 670 using ASAS-3 data.The horizontal lines indicate the FAP at levels of 0.1%, 0.5%, and 1%.The vertical red lines show the periods of the four signals with the highest significance.The time series are fitted by a sinusoidal with the period showing the highest significance. for

EFig. 12 .
Fig. 12. GLS periodogram of the RV, FWHM, BIS, and stellar activity indexes for IC4651 No. 9122 HARPS data.The blue dashed line marks the period of the planet candidate.The period of the strongest peak in each periodogram is marked with a red solid line.The horizontal line marks the 1% FAP level.The last panel shows the periodogram for the window function.

Fig. 13 .
Fig. 13.One Keplerian fit for IC4651 No. 9122 using kima.The top panel shows the phase curve for the signal, with the residuals shown below, highlighting the rms of the residual RVs.The bottom panel shows the RV data from HARPS and representative samples from the posterior distribution.The systemic RV has been subtracted and is shown at the top.

EFig. 14 .Fig. 15 .
Fig. 14.GLS periodogram of the RV, FWHM, BIS, and stellar activity indexes for NGC2423 No. 3 HARPS data.The blue dashed line marks the period of the planet candidate.The period of the strongest peak in each periodogram is marked with a red solid line.The horizontal line marks the 1% FAP level.The last panel shows the periodogram for the window function.

Fig
Fig. A.1.RV offset correction after fiber change as a function of FWHM (new observations have larger RV).Triangles are observed standard stars by LoCurto.The 60 sample stars with observations pre and post fiber change are coloured in blue, red or green circles (with their mean FWHM pre fiber change and after correcting from focus drift) and extrapolated using polynomials of 1st, 2nd or 3rd order fitting the standard stars, respectively.

Fig
Fig. A.2. FWHM values before and after the 2015 fiber change (marked with a black dashed line) for the 5 stars in this paper.The red points are the final FWHM values after considering the focus drift correction which only affect the data pre-2015.The horizontal blue dashed lines mark the median values before and after fiber change to calculate the observed FWHM offset.The bottom-right plot shows the FWHM offset (which is positive for most of the stars) for all the stars in our sample with more than 5 observations before and after the fiber change.
Fig. D.3.One-Keplerian fit for NGC4349 No. 127 using kima.The top panel shows the phase curve for the signal, with the residuals shown below, highlighting the rms of the residual RVs.The bottom panel shows the RV data from HARPS and representative samples from the posterior distribution.The systemic RV has been subtracted and is shown at the top.

Fig
Fig. D.4.Time series for the FWHM of NGC4349 No. 127 with a sinusoidal function fit with period 675 days.

Table 2 .
Gaia DR3 identification for each target.

Table 3 .
Pearson's ρ (P) and Spearman rank (S) correlation coefficients, together with their respective p-values, for the correlations between RV and different activity indicators of only the HARPS data.We mark in boldface the strong correlations (those that have a coefficient with an absolute value of higher than 0.5) and the significant correlations, even if weak (those with a p-value below 10 −2 ).

Table 5 .
(Jones et al. 2016)5)parameters for the planet candidates.From the six targets studied in this work, we find suggestions that it is more probable to find bona fide planets around less massive and less evolved stars.NGC3680 No. 41 and IC4651 No. 9122 both have a mass of around 1.7 M ⊙ and a radius of around 12 R ⊙ and seem to be in the first ascent of the RGB (see Fig.1).On the other hand, for stars more massive than 2 M ⊙ , we cannot confirm the planetary origin for any of the RV signals detected in our more massive targets.These four stars are also among the most evolved in their clusters.NGC2345 No. 50 is the most evolved of the full sample, already ascending towards the AGB.NGC 4349 No. 127 and NGC3532 No. 670 seem to be very close to the tip of the RGB and NGC2423 No. 3 is also approaching the RGB tip.This is in agreement with the decrease in planetary occurrence at M ≳2 M ⊙ found by the Lick(Reffert et al. 2015)and EXPRESS surveys(Jones et al. 2016)and the combination of both with the PPPS survey