Active red giants: close binaries versus single rapid rotators

The objective of this work is to determine what fraction of red-giant (RG) stars shows photometric rotational modulation, and understand its origin. One of the underlying questions is the role of close binarity in this population, standing upon the fact that RGs in short-period binary systems (<150 days or so) have been observed to display strong rotational modulation. We select a sample of about 4500 relatively bright RGs observed by Kepler, and show that 370 of them (8%) display rotational modulation. Almost all have oscillation amplitudes below the median of the sample, while 30 of them are not oscillating at all. Of the 85 of these RGs with rotational modulation chosen for follow-up radial-velocity observation and analysis, 34 show clear evidence of spectroscopic binarity. Surprisingly, 26 of the 30 non-oscillators are in this group of binaries. To the contrary, about 85% of the active RGs with detectable oscillations are not part of close binaries. With the help of stellar masses and evolutionary states computed from the oscillation properties, it appears that low-mass red-giant branch stars tend to be magnetically inactive, while intermediate-mass ones tend to be highly active. The opposite trends are true for helium-core burning (red clump) stars, whereby the lower-mass clump stars are comparatively more active and the higher-mass ones less so. In other words, we find that low-mass red-giant branch stars gain angular momentum as they evolve to clump stars, while higher-mass ones lose angular momentum. The trend observed with low-mass stars leads to possible scenarios of planet engulfment or other merging events during the shell-burning phase. Regarding intermediate-mass stars, the rotation periods are long with respect to theoretical expectations reported in the literature, which reinforces the existence of an unidentified sink of angular momentum after the main sequence.


Introduction
Thanks to the observations of the space photometers CoRoT (Baglin et al. 2009), Kepler (Borucki et al. 2011), and now TESS (Ricker et al. 2015), red giant (RG) stars are known to commonly display solar-like (SL) oscillations (De Ridder et al. 2009). As for main-sequence (MS) SL stars, these global oscillations are stochastically excited by turbulent convection in their outer convective envelopes. From a revised classification of the Kepler targets (Berger et al. 2018) based on the second data release of the ESA GAIA mission (Gaia Collaboration et al. 2018) and using the revised temperatures from the Kepler Stellar Properties Catalog (KSPC DR25, Mathur et al. 2017), Hon et al. (2019) detected oscillations in about 92 % of the 21, 000 RGs observed by the original Kepler mission. We note that this fraction is much larger than the 40 % detection rate reported for MS SL stars (Chaplin et al. 2011;Mathur et al. 2019), whose light curves are dominated by observational noise . Astero-seismology is thus a very powerful tool to analyze the evolution of stars after the MS as it can be performed on almost all targets of vast samples of stars.
Contrarily to MS SL stars (e.g., García et al. 2014a), only a few percent of the RGs show a pseudo-periodic photometric variability likely originating from surface spots (Ceillier et al. 2017), which we hereafter call "rotational modulation". The most commonly adopted explanation is that RGs are slow rotators (e.g., Gray 1981;de Medeiros et al. 1996) where no significant dynamo-driven magnetic fields are generated in their large convective envelope, implying no detectable surface spots. In this respect, it is worth noticing that no direct detection of magnetic field has been reported so far for stars on the red giant branch (RGB) with masses below 1.5 M (Charbonnel et al. 2017). Nevertheless, measuring the rotation rate of RG stars is of prime importance to understand the evolution of stellar rotation after the MS. For this reason, much effort has lately been devoted to measure their rotation rates from their radiative cores Article number, page 1 of 18 arXiv:2004.13792v1 [astro-ph.SR] 28 Apr 2020 (e.g., Beck et al. 2012;Mosser et al. 2012b;Deheuvels et al. 2015;Gehan et al. 2018) to their surfaces (e.g., Carlberg et al. 2011;Costa et al. 2015;Deheuvels et al. 2015;Tayar et al. 2015;Ceillier et al. 2017).
From high-resolution (HR) spectroscopic measurements of projected rotational velocities v sin i by Carlberg et al. (2011), about 2 % of RGs are rapid rotators, using a threshold of v sin i = 10 km s −1 as is commonly adopted in the literature (e.g., Drake et al. 2002;Massarotti et al. 2008;Costa et al. 2015;Tayar et al. 2015). Although the details of stellar dynamo theory are not yet fully understood, the empirical correlation between magnetic activity and rotation (e.g., Noyes et al. 1984) lends support to the expectation that the fraction of stars with detected rotational modulation should be similar to those with rapid rotation (see also Ceillier et al. 2017). Despite this correlation, as of yet there is no theoretical calculation that finds the condition for spot generation to be v sin i ≥ 10 km s −1 . However, Ceillier et al. (2017) reported the detection of starspot variability in about 2 % of the light curves from a systematic search in the Kepler oscillating RGs, which matches the fraction of rapid RG rotators of Carlberg et al. (2011). In the rest of the paper, we consider that RGs with rotational modulation are fast rotators relatively to the bulk of RG stars.
Currently, the observed rotation in both intermediate and low mass stars present puzzling features compared with the standard theoretical predictions. Intermediate mass stars (M ∈ [1.5 − 3] M ) do not have a convective envelope during the MS, which means no large scale magnetic fields, and hence no angular momentum loss during it (e.g., Durney & Latour 1978;Tayar et al. 2015). The rotation periods of A and late B-type stars cover a very broad range during the MS up to v sin i ∼ 300 km s −1 (Zorec & Royer 2012). When the rapid rotators reach the RGB, they develop a convective envelope while still spinning fast, and they keep having a significantly large rotation (v sin i ≥ 10 km s −1 ) during the subsequent core-helium-burning stage, also known as the red clump (RC) in reference to the colormagnitude diagram (e.g., Girardi 2016). This theoretical expectation is actually not observed in practice: Ceillier et al. (2017) find 1.9 % of M ≥ 2M RGs with rotational modulation, and Tayar & Pinsonneault (2018) report a distribution of v sin i peaking between 3 and 5 km s −1 with maximum values at 10 km s −1 . As regards low-mass stars (M ≤ 1.3M ), which have convective envelopes and thus lose angular momentum during the MS, we expect slow rotation rates when they evolve. From stellar evolution models, Tayar et al. (2015) indicate that low-mass stars should show v sin i ≤ 0.3 km s −1 during the RGB, and at maximum an order of magnitude larger on the RC. Surprisingly, measurements of v sin i by Tayar et al. (2015) and photometric modulation by Ceillier et al. (2017) reported that between 7 % and 15 % of RGs with masses M ≤ 1.1M are fast rotators. Such low-mass rapidly rotating RGs are thus suspected to have either been spun up by tidal interactions with a stellar companion, or have merged with or engulfed a stellar or sub-stellar companion.
In parallel to this, it has also been observed that RGs in close eclipsing binary (EB) systems have peculiar photometric properties (Gaulme et al. 2014). Among the 35 Kepler RGs that are confirmed to belong to EBs, 18 display regular solar-like oscillations -with their amplitudes matching the empirical expectations (e.g., Kallinger et al. 2014) -and no rotational modulation. All of the remaining systems display rotational modulation, 7 with partially suppressed oscillations, and 10 with no detectable oscillations (Gaulme et al. 2014. Gaulme et al. (2016) showed that the nondetection of oscillations is not an observational bias. This is ob-served in the closest systems, where most orbits are circularized, and rotation periods are either synchronized or in a spin-orbit resonance. Such a configuration is observed for systems whose orbital periods are shorter than about P orb ∼ 150 days, and the sum of the stellar radii relatively to the system's semi-major axis (R 1 + R 2 )/a 7 % of the semi-major axis. In the following, we consider a binary system to be close when P orb 200 days. Gaulme et al. (2014) suggested that the surface activity and its concomitant mode suppression originate from tidal interactions. During the RGB, close binary systems reach a tidal equilibrium where stars are synchronized and orbits circularized (e.g., Verbunt & Phinney 1995;Beck et al. 2018). Red giant stars are spun up during synchronization, which leads to the development of a dynamo mechanism inside the convective envelope, leading to surface spots. The magnetic field in the envelope likely reduces the turbulent excitation of pressure waves by partially inhibiting convection. Since spots can also absorb acoustic energy, these two effects lead to the suppression of oscillations.
This phenomenon observed with RGs in EBs led us to hypothesize that a large fraction of the RGs with rotational modulation -hereafter "active" RGs -and suppressed oscillations actually belong to short-period (non-eclipsing) binary systems. This connection between fast rotating RGs and binarity has been considered for a long time (e.g., Simon & Drake 1989;Carlberg et al. 2011), but no dedicated observations have been led to test it. Only archived data from large surveys were considered in recent works (Tayar et al. 2015;Ceillier et al. 2017), leading to rather inconclusive results on that matter.
The goal of the present paper is firstly to determine what fraction of active RGs are in close binary systems, and secondly to understand the origin of magnetic activity in the RGs that do not belong to close binaries. We base our work on a selfconsistent approach where we select a sample of about 4500 RGs observed by Kepler that are not known to belong to any binary system, and that are expected to exhibit SL oscillations (Sect. 2). Among them, we show that about 370 (≈ 8 %) display significant rotational modulation, which appears to be correlated with weak oscillations (Sect. 3.1). We then look for radial velocity (RV) variability in a subset of 85 targets from HR spectroscopic measurements (Sect. 3.2), and find a significant fraction of spectroscopic binaries (SBs). In the last section (3.3), we classify and quantify the possible scenarios that could have led to the active RGs that are not in binary systems.

A sample of 4500 RGs observed by Kepler
The first step of our work consists of measuring the fraction of RGs displaying clear rotational modulation, and reassessing the correlation between rotational modulation and suppression of SL oscillations pointed out by e.g., García et al. (2010), Chaplin et al. (2011), Gaulme et al. (2014), or Mathur et al. (2019. Our main concern at this step is to avoid observational biases by including targets that are not suitable for our study. We started from the updated RG catalog of Berger et al. (2018), and selected a subsample of targets whose oscillations should be detectable provided that they are regular RGs. We base our selection criterion on the conclusions of Mosser et al. (2019), who quantified the asteroseismic performance of Kepler as a function of apparent magnitude and observation duration. This led us to consider RGs with Kepler magnitude m Kep ≤ 12.5 and observation duration longer than 13 quarters (≈ 3 years). We also limit ourselves to RGs with radii less than 15R to avoid be-  Berger et al. (2018) in the range [0.2, 100] R and [3000, 7500] K. Radii are from GAIA and temperatures from the Kepler MAST archive. Temperature axis is reverse, that is with higher T eff toward left. Color-coding represents logarithmic number density. The colored part of the diagram indicates the location of our sample. The red line indicates the limit between subgiant and giant stars according to Berger et al. (2018), such as the RGs lay at the top-right corner. According to Berger et al. (2018), the discontinuity in T eff near 4000 K is an artifact due to systematic shifts in T eff scales in the DR25 Kepler Stellar Properties Catalog.
ing biased by targets whose oscillations at maximum amplitude are at frequencies ν max 15 µHz, which could be filtered out during light-curve processing. We also exclude those with radii less than 4R to avoid oscillation modes beyond Kepler's longcadence Nyquist frequency (≈ 283 µHz). We note that 75 % of the Kepler RG sample have radii between 4 and 15 R . In other words we focus on the "mainstream" RGs. A total of 4580 out of the ≈ 21, 000 RGs meet these criteria (See Fig. 1).
As misclassification can happen in large catalogs, we excluded stars that were clearly not single RGs from a visual inspection of the Kepler light curves and their frequency spectra. This downselection filtered out obvious binary systems as double RG oscillators, eclipsing binaries, ellipsoidal binaries, and highly eccentric non-eclipsing binaries ("heartbeat" stars", e.g., Welsh et al. 2011;Beck et al. 2014), but also classical pulsators as δ Scuti or γ Doradus. We found 88 light curves including a binary signal, and 21 with either a δ Scuti or a γ Doradus pulsator. We identified another 6 with a damaged Kepler light curve, that is where photometry is clearly wrong. We thus ended up with a sample of 4465 stars classified as RGs according to Berger et al. (2018) and with no obvious indication of misclassification or peculiarity.

The Kepler light curves
We primarily worked with the Kepler public light curves that are available on the Mikulski Archive for Space Telescopes (MAST) 1 . Two types of time series are available: the Simple Aperture Photometry (SAP) and the Pre-search Data Conditioning Simple Aperture Photometry (PDC-SAP) light curves. The latter consist of time series that were corrected for discontinuities, systematic errors and excess flux due to aperture crowding (Twicken et al. 2010 during the process (e.g., García et al. 2014b;Gaulme et al. 2014).
We thus made use of the SAP data to preserve any possible longterm signal. This choice entails our own detrending and stitching operation on the light curves while ensuring that the rotational modulation is preserved after each interruption of the time series. The methods employed to clean the time series are detailed in Gaulme et al. (2016).
To reinforce our results -especially the flagging of rotational modulation and the measurement of rotational periods -, we repeated the whole light processing and analysis with the public Kepler Light Curves Optimized For Asteroseismology (KEP-SEISMIC, García et al. 2011;Pires et al. 2015) 2 , which are optimized for asteroseismic analysis. KEPSEISMIC light curves were available for 4444 out of the 4465 targets. Except for a few cases, the results obtained with either light curves were in agreement. In the following, our results reflect the final analysis done by working both the SAP and KESPSEISMIC time series.

Surface rotation periods
In a first step, the presence of rotational modulation is automatically searched, before being checked by visual inspection. Rotational modulation is considered whenever the standard deviation of the photometric time series is larger than 0.1 %. We note that this definition is equivalent to the photometric index S ph defined by Mathur et al. (2014) in the case no photometric period is measured. Then from both the power spectrum and the autocorrelation of the times series, the algorithm automatically checks for surface activity, and provides an estimate of its fundamental period in case of positive detection. Rotational modulation is considered for periods ranging from 1 to 180 days. From the power spectrum, peaks are considered to be significant when the statistical null hypothesis test is less than 1 % with respect to the average background noise level. Autocorrelation is used to check the estimate obtained from the power spectrum, by verifying that the highest peak lies in the same period range within a factor 1/2 or 2, as misestimates can happen. Visual inspection was then used to decide the actual rotation period in case of disagreement between power spectrum and autocorrelation periods.
Beyond measuring rotation periods, we need to be careful with possible photometric contamination from nearby stars. Indeed, it can happen that a considered RG does not have any spot but that a neighboring MS stars display spots and that the light curves leak into each other. Ceillier et al. (2017) tracked this possible source of errors by making use of the crowding parameter available on the Kepler database at MAST. The crowding parameter is defined as the fraction of the flux coming from the considered target in the light curve. The lower the crowding parameters, the higher the chance of contamination. In their case, they considered all the 17,377 RGs with detected oscillations at the time, which included many faint targets. Their Fig. 5 shows the comparison between the histogram of rotational periods with and without targets where the crowding factor is less than 0.98. In their case, they found that the targets with low crowding represent 32 % of their active RG sample, and are a significant source of bias as they constitute a clear group of outlying P rot with values less than 30 days. In our case, Fig. 2 shows that the stars with a crowding factor less than 0.98 (13 % of our active RGs) do not constitute a well defined group, so we can legitimately consider that targets with a lower crowding factor are not due to contamination in most cases, and that we are not affected by this problem. It is not surprising since we focus on purpose on the brightest part of the sample to avoid observational biases of that kind.

Oscillation analysis
The detection of oscillations is based on the envelope of the autocorrelation function (EACF) developed by Mosser & Appourchaux (2009), which we apply to the power spectrum of the time series ( Fig. 5 panel i). A first run of the EACF leads to initial estimates of the mean frequency separation ∆ν and ν max , independently from the background noise level. The second step consists of refining ∆ν and ν max by taking into account the background noise. The stellar granulation and accurate values of ν max and mode amplitude H max are estimated by fitting the power spectrum as commonly performed in asteroseismology (Kallinger et al. 2014), and already used in Gaulme et al. (2016). Following Kallinger et al. (2014), the power density spectrum is fitted by the following function: where the N is the function describing the noise, η is a damping factor originating from the data sampling, B is the sum of three "Harvey" functions (super Lorentzian functions centered on 0), and G is the Gaussian function that accounts for the oscillation excess power: The terms ν max and H max are the central frequency and height of the Gaussian function. A second estimate of ∆ν is performed again with EACF from the whitened power spectral (power spectrum divided by background function). Oscillations are considered to be detected when both the maximum value of the EACF is larger than 8 (Mosser & Appourchaux 2009) and peaks are visible in the spectrum. We then compute the stellar masses and radii thanks to the asteroseismic scaling relations that were originally proposed by Kjeldsen & Bedding (1995) for SL MS oscillators, and then successfully applied to RGs (e.g., Mosser et al. 2013). The stellar effective temperatures, also needed in the scaling relations, are those from the Kepler GAIA data release 2 with their corresponding errors. We employ the asteroseismic scaling relations as proposed by Mosser et al. (2013) for RGs, i.e., where ν max, = 3104 µHz, ∆ν = 138.8 µHz, T eff, = 5777 K, and where the observed ∆ν is converted into an asymptotic one, such as ∆ν as = 1.038 ∆ν.
In addition to the properties of the pressure p modes (ν max , ∆ν) we analyze the properties of the mixed modes, which result of the interaction between the p-modes that resonate in the convective envelope and the gravity g modes that resonate in the radiative core. At first order, p modes are evenly spaced by ∆ν as a function of frequency, whereas g modes are evenly spaced by ∆Π as a function of period. From the dipole (l = 1) mixed The dark area corresponds to the correlated signal -i.e., the oscillations -where its abscissa indicates ν max and its ordinate 2/∆ν. Panel (j): échelle diagram associated with the large frequency spacing automatically determined from the EACF plot. The power density spectrum is smoothed by a boxcar over three bins and cut into ∆ν = 11.9-µHz chunks; each is then stacked on top of its lower-frequency neighbor. This representation allows for visual identification of the modes. Darker regions correspond to larger peaks in power density. The x-axis is the frequency modulo the large frequency spacing (i.e., from 0 to ∆ν), and the y-axis is the frequency. modes it is possible to measure the mean period spacing of the dipole g modes ∆Π 1 , which is related to the physical properties of the core. For RGs this information is precious as it allows us to unambiguously distinguish the H-shell (RGB) from the Hecore (RC) burning stars. Since the CoRoT space mission, such a type of analysis is commonly performed with RG stars Mosser et al. 2011;Stello et al. 2013). According to ) when ∆ν ≥ 9.5 µHz, we assume the stars to be on the RGB. When mixed modes are detectable, the evolutionary status is determined according to the methods described in Vrard et al. (2016). Overall, we have the status of 3405 stars out of the 4465: we used 2307 published statuses that were already available in Vrard et al. (2016), and we determined the status of 1098 others. Figure 3 shows the period spacing of the dipole mixed modes ∆Π 1 as a function of the large frequency separation ∆ν for the subsample of active RGs. We note a larger dispersion of the period spacing ∆Π 1 for RGBs with ∆ν ≤ 7 µHz. This issue is well known (Dupret et al. 2009;Grosjean et al. 2014): RGBs with low ∆ν are more evolved than those with larger values, and the coupling between p and g modes is less efficient. Hence, gravity-dominated mixed modes display lower signal-to-noise ratio.
Finally, our analysis of the oscillation spectra also considers the question of the RGs with depleted dipole modes (e.g., Mosser et al. 2012a;García et al. 2014b), which is possibly caused by the presence of fossil magnetic fields in radiative cores of intermediate mass RGs that had convective cores during the MS (Fuller et al. 2015;Stello et al. 2016;Mosser et al. 2017). Despite the suppression of dipole modes is a priori not related to rotational modulation -the former is supposed to be caused by core magnetic fields, the latter by surface ones -, we thought it was interesting to flag it in case of an unlikely correlation. We applied the method presented by Mosser et al. (2012a) to measure the relative visibilities of the modes of degrees l = 0, 1, 2, and 3 to the whole sample, whenever oscillations were detected. Oscillation height H max (ppm 2 µHz −1 ) as a function of the frequency at maximum amplitude ν max (µHz) of the 4465 red giants considered in this work, minus the 30 for which no oscillations were detected, i.e. 4435 total. The parameter H max is the height of the Gaussian function used to fit the oscillation contribution to the power spectral density. The frequency at maximum amplitude ν max is the center of the same Gaussian function. Light gray points indicate all the targets, green are l = 1 depleted, red squares are the systems with significant surface activity, and blue circles indicate the cases of ambiguous detection of surface activity.

A dedicated pipeline to analyze the sample
We developed an automatic pipeline -the Data Inspection Tool, hereafter DIT -to analyze the whole sample, which cleans the light curves and extracts the properties of both rotational modulation and solar-like oscillations. The code is similar to that described in Gaulme & Guzik (2019) to look for stellar pulsators in eclipsing binaries. Beyond the measurements of rotational modulation, background and oscillations, the DIT produces a vetting sheet for each target that includes a series of figures allowing the user to appreciate the quality of the measurements. An example of vetting sheet is displayed in Fig. 5. All vetting sheets were visually inspected. The objective of the visual inspection was to double check the detection of rotational modulation, assess the value of the rotation period, and look for oscillations that had been missed by the DIT. Besides, it allowed us to exclude the obvious binaries and classical pulsators as explained in the previous section. Regarding rotational modulation, we did not strictly stick to the S ph ≥ 0.1 % threshold, as some photometric discontinuities could cause S ph larger than that. Besides, some actual photometric modulation could have lower S ph . The distinction between a spurious photometric modulation, arising from imperfect light curve cleaning, and actual stellar signal was most of the time clear. However, in conditions of low signal-to-noise ratio (S/N), flagging the detection of stellar activity could be fragile. This is why we used two classes of flags for activity: clear and ambiguous detection (See Sect. 3.1).

High-resolution spectroscopy
As indicated in the introduction, we study a subsample of RGs displaying rotational modulation with HR spectroscopy to look for RV variability. Our intend is not to perform a detailed analysis of these spectra such as retrieving chemical composition or v sin i, which will be part of a future work.
From October 2018 to June 2019, we were granted time on the échelle spectrometer of the 3.5-m telescope of the Astrophysical Research Consortium (ARC) at Apache Point observatory (APO), which covers the whole visible domain at an average resolution of 31, 000. This time allowed us to monitor a set of 51 RGs with clear rotational modulation. Even though the ARC échelle spectrometer was not designed for fine RV measurements, it has successfully been used for this purpose in earlier works about RGs (e.g., Rawls et al. 2016;Gaulme et al. 2016). The measurement error reported in these papers is about 0.5 km s −1 for an RG spectrum with S /N between 10 and 20. In practise our spectra have S/N ranging from 10 to 25. With ARCES data, we thus consider a target to be an SB when the dispersion is larger than σ RV = 2 km s −1 . The ARCES Optical spectra were processed and analyzed in the same way as in Gaulme et al. (2016) and we refer to this paper for details.
Beyond the ARC telescope, we used infrared high-resolution spectroscopic data of another 34 targets that were available in the archived data (release 16, Ahumada et al. 2019) of the APO Galactic Evolution Experiment (APOGEE, Majewski et al. 2017), which is part of the Sloan Digital Sky Survey IV (e.g., Blanton et al. 2017). Regarding APOGEE data, we directly use the archived RVs. This spectrometer has an average resolution of 22, 500 and is known to be stable, with a noise level lower than 0.1 km s −1 (Deshpande et al. 2013). Knowing that the RV jitter of RGs with surface gravities 2 ≤ log g ≤ 3 -as we have hereis less than 0.1 km s −1 (Hekker et al. 2008), we consider an RG to be an SB when the RV dispersion is larger than σ RV = 0.5 km s −1 .

Occurrence of red giants with rotational modulation
Among the 4465 RG stars that we have selected, our main conclusions regarding their Kepler light curves are the following: 1. A total of 298 RGs display clear rotational modulation, i.e., in 6.65 % of the cases. For another 73 targets (1.63 %), the detection of surface activity is ambiguous (low signal-tonoise ratio, S /N). If we assume the ambiguous detections as a proxy of our error in judgment, we can consider that 7.5 ± 0.8 % of the RGs display rotational modulation. We thus detect from three to four times more RGs with rotational modulation than the 2 % reported by Ceillier et al. (2017). As indicated in Sect. 2.3, the main difference with their work is that they considered the whole sample of oscillating RGs known at the time, which included faint stars and likely discarded stars with suppressed oscillations. By comparing Fig. 7 with Fig. 4 of Ceillier et al. (2017), we clearly detect more stars in the ∆ν range [6.5, 9.5] µHz, which contains many intermediate-mass RCs. 2. We detect oscillations in 99.32 % of our sample, i.e. in all except 30, which is much more than previously reported. This difference is not surprising because Hon et al. (2019) considered the whole 21,000 RG catalog which contains faint stars, including crowded fields. Besides, their automatic pipeline may have encountered difficulties when dealing with low S /N oscillation spectra, especially in the presence of rotational modulation. We interestingly note that the fraction 8 % of RGs where Hon et al. (2019) do not detect oscillations is equal to the fraction of our RGs displaying activity. 3. All of the 30 targets with no detectable oscillations display a significant rotational modulation. In other words, if an RG does not oscillate, rotational modulation is always detected. 4. Almost all the active RGs display unusually low-amplitude oscillations, as well as most of those with ambiguous status. This can be seen in Fig. 6, which plots a proxy of the oscillation amplitude -the height H max of the Gaussian function employed to fit the oscillation envelope in the power spectral density of the time series -as a function of ν max . We distinguish three groups of oscillators: the regular, the dipoledepleted, and the active ones. The dipole-depleted RGs appears as a distinct group because if l = 1 modes are absent, the oscillation power that we fit with a Gaussian function is reduced 3 . 5. Measurements of dipolar mode visibilities as performed in Mosser et al. (2012a) and Mosser et al. (2017) reveal that 17 % and 10 % of inactive and active RGs display l = 1 depleted modes, respectively. To complement these results, a visual flagging of depleted dipole modes led us to identify 12 % and 20 % in the same two groups. We conclude that the fraction of depleted dipolar-mode oscillators is very similar in both samples of inactive and active RGs (Fig. 6), and cannot be considered as a significant source of biases in our analysis.

Occurrence of close binaries among the active red giants
According to the previous section, 7.5±0.8 % of the RGs display rotational modulation. We now determine how many of them 3 Dipole-depleted RGs would not appear as a distinct group if our proxy of the oscillation amplitude were the amplitude of radial modes.  actually belong to close binary systems thanks to multi-epoch high-resolution spectroscopic measurements. The rate of SBs actually represents a lower limit of the fraction of close binaries because some systems are seen almost pole-on and display very low line-of-sight velocities. However, the RV semi-amplitude of the typical binary system that we expect -synchronized, circularized, period of about 50 days, solar masses -is larger than about 5 km s −1 in 90 % of the cases, by assuming random inclination angles (see Appendix A). In addition, we show that five measurements randomly distributed over a 200-day range with at least seven days between consecutive observations are very likely to detect RV variability. With the ARCES spectrometer, we observed four times most of our targets, and up to seven times some targets whose RV variability was confusing and needed more data. The number of visits per star is reported in Table B.1.
As indicated in Sect. 2.6, we could monitor 51 targets with ARCES on the 3.5-m ARC telescope 4 , and make use of 34 archived RV values from the APOGEE/SDSS IV survey 5 . For ARCES, we selected the brightest targets among the active ones to maximize the number of observable stars, and we made sure that a significant fraction of the targets did not show any oscillations. Overall, our spectroscopic sample is composed of 85 RGs with clear rotational modulation, including 30 with no oscillations and 55 with oscillations. The result appears to be very different depending whether oscillations are detectable (Fig.  8). Among the 30 non-oscillating stars, 26 are clear SBs, i.e. ∼ 85 %. In contrast, among the 55 oscillating ones, only 8 or 9 are SBs, i.e., ∼ 15 %.
To check the significance of the SB fraction, we looked for RV data of a null sample, where no rotational modulation is observed. A total of 212 RGs observed more than four times are available in the APOGEE DR16 archive. Among them two are SBs, that is about 1 %. We therefore consider the enhanced fraction of close binaries among the RGs with rotational modulation (and weak oscillations) to be significant.
We are aware that the incidence of binarity is often reported to be over 50 % (e.g., Eggleton 2006). However, this number includes a broad variety of systems, from contact binaries with orbits as short as 0.2 days to wide ones with periods of thousands of years. Since we focus on binaries with at least one member on the RG phase with orbital periods between about 15 to 200 days, our null sample is not representative of the absolute incidence of binaries. The catalog of approximately 3,000 EBs discovered by Kepler (e.g., Prša et al. 2011;Kirk et al. 2016) helps understanding the order of magnitude of the rate of binaries that we find. The orbital periods of Kepler EBs range from about 0.1 to 1100 days. Only 20 % of them have orbits longer than 15 days and about 50 systems have an RG component (Gaulme & Guzik 2019), that is ∼ 2 % of the total. By assuming a 50 % rate of binaries we thus expect about 1 % of SBs with an RG component and an orbit longer than two weeks.

Rapid rotation without close binarity
Among the RGs with rotational modulation, it thus appears that about 85 % of the active RGs with detectable oscillations are not part of close binaries. With the help of additional information, we shed light on the origin of their activity.

Two groups of single active RGs
We retrieve the masses, radii and evolutionary status of the giants thanks to an asteroseismic analysis described in Sect. 2.4. By evolutionary status, we mean the RGB where stars burn H 4 Actually, we monitored 53 targets with ARCES but one (KIC 11551404) is not part of the sample of 4465 RGs because its GAIA radius is too small (2.7 ± 0.2R ). We picked it before Berger et al. (2018) was published because it displayed strong rotational modulation and was part of earlier Kepler RG catalogs. The RV dispersion clearly indicates it is a close binary system. We also observed KIC 3459199 which actually has an EB signal in its light curve. It is also an SB. We decided to keep these two in Table B.1 for the records. 5 347 RGs displaying rotational modulation were actually available from the APOGEE database, but only 36 were observed more than three times over a time span longer than 10 days, and 4 were already part of our ARCES sample. Two more were observed only twice -both are non-oscillating -but it was enough to prove they are SBs. in shell around an inert He core, versus the RC, where He burning occurs in the core. We distinguish the red clump 1 (RC1) which consists of lower-mass stars (M 2M ) that reached the clump after the He flash, from the red clump 2 (RC2) composed   Table A.1). Figure 9 shows the distributions of mass and radius for both the inactive and the active samples. A striking result is the very different mass distribution between the two samples (radius distributions are more similar). For the inactive RGs, the mass distribution is consistent with previous reports based on CoRoT or Kepler observations (e.g., Miglio et al. 2013). We note that it peaks at 1.3M for RGBs and 1.15M for RC, which could be compatible with a significant mass loss during the late stages of the RGB phase. It also could be explained by the fact that we employed a unique asteroseismic scaling relation, whereas they should not be exactly calibrated in the same way for RGB and RC stars (Miglio et al. 2012;Khan et al. 2019). Regarding the active sample, we first note a general depletion of stars from 1.2 to 2 M . In particular, the RGB sample appears to be depleted of low-mass stars and overpopulated of intermediate mass stars. The RC population is bimodal with a large group of RC1 stars on a narrow range (0.8-1.1 M ), and a large group of (mostly RC2s) on a broad range (1.8-2.8 M ). Hereafter we often distinguish two groups by referring to low-mass RGBs and RCs (LM-RCs and LM-RGBs) with M ≤ 1.5M , and high-mass RGBs and RCs (HM-RGBs and HM-RCs) above.

Fraction of active RGs as a function of mass
If we consider the RGs with unknown evolutionary status to be equally distributed among the RGBs and RCs, we can estimate the absolute fraction of RGs in a given evolutionary status that display surface activity. We note them r RGB,act and r RC,act . For instance, for the RGB group: where the fraction of active RGs with a known evolutionary status r act,status = N act,status /N act = 200/340 = 58.8 %, and the fraction of inactive RGs with a known status is r inac,status = 3205/4077 = 78.6 % (see Table A.1). That way, 6.6 % of RGBs and 8.2 % of RCs are active. However, we observe a strong variation of the occurrence of active RGs as both a function of mass and evolutionary status ( Fig. 10 and Table A.1). Let us first focus on the low mass part of the sample (M ≤ 1.1M ) that was considered by Tayar et al. (2015) and Ceillier et al. (2017). Consistently with these two papers, we estimate the fraction of rapidly rotating RGs to be r act = 10.5 %. Our study brings new information thanks to the evolutionary status: it arises that only 3.4 % of the RGBs are active, against 11.7 % of the RCs. The fact that the fraction of low-mass active RGBs is two to three times smaller than that of active RCs is a strong indicator that the observed rotational modulation originates from an increased rotational rate happening on the RGB. Engulfment or merging of a stellar or substellar companion is a likely hypothesis (Simon & Drake 1989).
As regards stars with M ≥ 2M , we find that about 60 % of the RGBs and 25 % of the RC2s are active. First of all, the smaller fraction of active RCs with respect to RGBs is compatible with loss of angular momentum starting during the RGB. Secondly, such a large rate of active RGs among intermediatemass RGs was not reported by Ceillier et al. (2017). As commented earlier, they likely missed this fraction of the sample because of catalog selection biases. At the same time, we confirm the conclusions of Tayar & Pinsonneault (2018) who report that there is no large fraction of rapid rotators among intermediate mass RGs according to spectroscopic v sin i measurements if we stick to the 10 km s −1 criterion. Indeed, by combining the rotation periods with the asteroseismic radii, we compute the rotation velocities at equator v sin(i = 90 • ). Figure 11 (middle panel) shows v sin(i = 90 • ) for the whole sample and for the individual evolutionary stages. Regarding the HM-RC stars, only ≈ 4 % show v sin(i = 90 • ) ≥ 10 km s −1 .
Interestingly, we notice that the only group that shows a significant fraction (≈ 37 %) of stars with v sin(i = 90 • ) ≥ 10 km s −1 is the LM-RC group, which is suspected to have engulfed a stellar or substellar companion. Finally, by considering all of our RGs with detected rotational modulation and oscillations, the fraction of stars with v sin(i = 90 • ) ≥ 10 km s −1 is about 17 %. By assuming that all of the inactive RGs have v sin i < 10 km s −1 , it means that 1.4 % of our whole sample show v sin(i = 90 • ) ≥ 10 km s −1 , which is a little less than the 2 % reported by Carlberg et al. (2011).

Rotation rate versus activity
Firstly, we observe a strong difference of the photometric activity index S ph between the non-oscillating (almost all close binaries) and the oscillating active RGs (mostly single RGs) from our sample. From Figs. 11 and 12 Photometric index S ph (percent) as a function of rotation period P rot (days) for the RGs with rotational modulation. The photometric index S ph is plotted in log scale, P rot in linear scale. Blue squares indicate the non-oscillating RGs, red disks the RGBs, darker green upwardpointing triangles the low-mass RCs, and lighter downward-pointing triangles high-mass RCs. The markers that display black edges are confirmed spectroscopic binaries. Black dots indicate RV stable RGs. The gray background reflects the distribution of S ph of the inactive stars (the darker, the higher).
and Benbakoura et al. (submitted), the non-oscillating RGs in EBs are almost all both synchronized and circularized. The few non-locked systems are pseudo-synchronized. It means that the rotation rate is not the only factor driving the amplitude of rotational modulation and oscillations: tidal locking causes in a way or another more photometric contrast. Understanding that aspect goes beyond the scope of this paper. Regarding the single active RGs, our results show that rotational modulation appears at rotation rates much smaller than the "fast rotation" threshold that is usually considered. We indeed detect rotational modulation for RGs with v sin i ≥ 2 km s −1 . To put this result in its usual theoretical context (e.g., Noyes et al. 1984;Wright et al. 2011), we consider the generation of magnetic field in the frame of the turbulent dynamo scenario (e.g., Charbonneau 2014). In such a scenario, the efficiency of the magnetic field generation is characterized by the Rossby number Ro = P rot /τ c . The efficient dynamo regime, and hence surface magnetic fields and spots, is associated with Ro ≤ 1. We thus need an estimate of the convective turnover timescale τ c for the stars in our sample.
In Fig. 13, isolevels of τ c are overplotted with the data; evolutionary tracks are also shown for reference. We use the τ c extracted from the stellar evolution tracks from the YaPSI database (Spada et al. 2017). It should be noted that the values of τ c shown in the Figure were calculated using the mixing-length treatment of convection (see Spada et al. 2017 for details), which is adequate to provide an estimate to the leading order.
From Fig. 13, it is apparent that all of our ≈ 4500 RGs could be active if their rotational periods were less than about 250 days. In this T eff -log g diagram (Kiel diagram), active and inactive RGs overlap, which tends to indicate that the inactive stars rotate slower, likely with rotation periods longer than about 250 days. Among the RGB stars (panel a) the active sub-sample is concentrated around the region of the first dredge-up, i.e., near the base of the RGB. For the RC stars (panel b of Figure 13), two distinct sub-populations of active stars exist, roughly corresponding to the low-mass and the intermediate mass regimes. The former have log g ≈ 2.4, while the latter have larger surface gravity and are distributed with approximately the same scatter in log g and T eff . Both characteristics are consistent with the predictions of classic stellar evolution theory for the RC and the secondary RC, respectively (see, e.g., Girardi 2016).
This result also needs be compared with spectro-polarimetric observations of magnetic RGs (e.g., Aurière et al. 2011Aurière et al. , 2015Konstantinova-Antova et al. 2012 and the subsequent the- describes spectropolarimetric observations of a sample of 48 magnetically active RGs, and report clear detection of magnetic fields in 29 stars. Rotation periods are known for most of these stars, and range from about 6 to 600 days. They observe that the magnetic RGs are distributed along a so-called "magnetic strip" in the Hertzsprung-Russell diagram, which is a region where τ c is maximum. Such a condition is met during the first dredge up phase at the bottom of the RGB, as well as for RC stars. This result is fully compatible with ours. The only notable difference with their results is that they identify a correlation between the chromospheric activity and the stellar rotation periods. We do not find such a trend between S ph and P rot (Fig. 12). This could be due to the fact that their sample spans a broader range of periods (most of our single RGs show periods from 40 to 170 days). Moreover, the chromospheric index S index and the Zeeman mea-surements are not sensitive to the stellar inclination contrarily to S ph , which makes trends clearer to appear. Aurière et al. (2015) also identify a few outliers that could be descendent of magnetic Ap stars (e.g., Aurière et al. 2011;Borisova et al. 2016;Tsvetkova et al. 2019). These outliers are characterized by an enhanced chromospheric activity S index with respect to their peers with similar rotation periods. We identify one star among our single RGs with strong activity that could match this description. KIC 5430224 has S ph = 1.1 % but is RV stable with no particularly fast rotation (P rot ≈ 55 days). It shows very weak oscillations. To the contrary the activity of KIC 8160175 (S ph = 2.6 %) seems to be caused by contamination of a nearby star (see Appendix B) and cannot be counted as a possible Ap descendent.

Conclusions
The original objective of this work was to determine what fraction of RGs shows activity (rotational modulation), and understand its origin. One of the underlying questions was the role of close binarity in this population, standing upon the fact that RGs in close binary systems (P orb 200 days) have been observed to display strong rotational modulation.
To avoid as much as possible being influenced by observational biases, we carefully selected a subsample of the RGs observed by the Kepler satellite during its original four-year mission. This sample was picked from the Berger et al. (2018) stellar classification of the Kepler field based on the GAIA DR2. From their sample of RGs, we selected the targets whose light curves should not be limited by the photon noise, meaning that the oscillations of a regular RG should be detectable. This led us to consider the brightest stars (m Kep ≤ 12.5 mag) that were observed the longest (more than 3 years). We added a cut on radii, to make sure that the oscillation range would fall between 15 µHz and the sampling cut-off (Nyquist frequency) at ≈ 283 µHz. The final sample is composed of 4465 RG stars (Fig. 15).
Our first result is the clear detection of SL oscillations in 99.3 % of the sample, which is a much larger fraction than reported in previous studies (e.g., Stello et al. 2013;Hon et al. 2019). We explain this mainly thanks to the fact that our special sample is not affected by photon noise for searching for RG oscillations. We also paid much attention to search for oscillations in low S/N conditions. The second important finding is that between 6.7 % (only clear detection) and 8.3 % (including low S/N detection) of RGs display rotational modulation, which is much larger than reported by previous studies, in particular by Ceillier et al. (2017). We also clearly show that the active RGs present oscillations with lower amplitudes with respect to the inactive sample (Fig. 6). This latter fact -increased activity leads to lower modes -was previously observed for MS stars (e.g., Chaplin et al. 2011;Mathur et al. 2019), but also for RGs in close multiple systems (e.g., Derekas et al. 2011;Gaulme et al. 2014, and Benbakoura et al., submitted.).
From high-resolution spectroscopic observations specially conducted with the échelle spectrometer of the 3.5-m ARC telescope, as well as public SDSS/APOGEE observations, we could determine the role of close binarity in the origin of the active RGs. It appears that among the 30 non-oscillating active RGs, 26 are clear SBs, one lacks in statistically significant RV variations, and another three are possible SBs. This means that almost all of the non-oscillating RGs belong to close binary systems. Among the rest of the sample, that is, the active RGs with partially suppressed but detectable oscillations, the fraction of close binaries is still significantly larger than in among the inactive ones. We count about 15 % of SBs whereas we identify only 1 % in a test sample composed of inactive RGs. This also means that a large majority of active RGs with detectable oscillations do not belong to close SBs.
We shed light on the origin of the active single RGs by analyzing their distributions in terms of physical properties (mass, radius) and evolutionary status (RGB, RC). It arises that most of the active RGs fall into two distict groups: a low-mass (M ≤ 1.5M ) and an intermediate-mass range (1.8 ≤ M ≤ 2.8M ). The low-mass part of the active RGs is mostly composed of stars on the RC: about 3 % of the RGBs are active while about 12 % of the RCs are active. This means that a fraction of the RGs in that mass range gains angular momentum between the RGB and the RC. This observation favors the scenarios of planet engulfment or merging with small stellar objects. Indeed, low-mass stars ignite Helium after the He flash, and present radii up to ∼ 200 R at the tip of the RGB (Fig. 14), and may swallow planets, converting their orbital momentum into spin and friction loss. Another option is a stellar engulfment. Highly eccentric systems (0.3 < e < 0.8) as those composing the heartbeat systems (e.g., Beck et al. 2014;Kuszlewicz et al. 2019) are typically composed of a Sun-mass star and a much smaller companion (0.2 M ). Near the tip of the RGB, the most massive star has a radius larger than the periastron distance and thus merges with the companion.
If this scenario of planet or stellar engulfment is true, it is interesting to notice that it regards a relatively narrow mass range. This can be qualitatively explained by three factors. Firstly, since more massive stars (M ≥ 2M ) do not reach the clump through the He flash, their radii do not exceed ≈ 20 R , so they are less susceptible to engulf planets. Moreover, because they are more massive, they gain relatively less angular momentum by swallowing planets than less massive stars. Lastly, the orbital angular momentum of close-in planets increases as a function of the semi-major axis, and thus is less efficient with a close planet with respect to a planet like Venus, orbiting at 150 R .
As regards intermediate-mass stars, we observe that a large fraction of the RGBs is active (over 50 % when M ≥ 2M ), and a large but reduced fraction of active RCs (about 25 %).
In opposition to the low-mass range, this observation implies that intermediate-mass stars tend to loose angular momentum between the RGB and the clump. Both the fraction of the active RGB in this mass range and the fact that they loose momentum between the RGB and the RC are compatible with the expected scenario for intermediate-mass stars. They do not loose angular momentum during the MS, thus still spin fast on the RGB, and a little less once on the clump. However, the values are still in contrast with recent theoretical studies (e.g., Tayar & Pinsonneault 2018). From the observed distribution of rotation rates of A and B-type MS stars, about half of the RGs are expected to display v sin i ≥ 10km s −1 (e.g., Ceillier et al. 2017). From our asteroseismic radii and photometric rotation periods, we measure that only about 4 % of the active intermediate-mass RGs display such a rapid rotation, i.e. about 2 % of that mass range by including the inactive ones. This reinforces previous reports of an unidentified sink of angular momentum after the MS.
Regarding the connection between the emergence of surface activity, we observe that RGs with v sin i ≥ 2km s −1 show rotational modulation. This is actually not surprising when we compare the rotation periods of the active RGs with the convective turnover time τ c . In a turbulent dynamo scenario, the generation of magnetic fields starts being efficient if the rotation period is lower than τ c . All the active RGs meet these conditions (Fig. 13). A correlated interesting fact is that the active RGs that belong to binary systems and that do not display oscillations show a photometric activity index S ph about an order of magnitude larger than single active RGs. We could explain this by the fact that the non-oscillating active RGs have rotation periods shorter in average than the single active RGs. However, we observe that single RGs with rotation periods between 30 and 80 days still show an S ph about ten times less than the RGs in SBs with the same periods. This implies that tidal locking somehow leads to larger magnetic fields, and also a more efficient oscillation suppression. Investigating that aspect will be part of future works. Last, Mathur et al. (2019) report that MS stars with S ph ≥ 0.2 % never show SL oscillations. We note that this threshold is not the same for RGs, as we detect oscillations up to levels of about S ph ≈ 1 %.
Further developments of this project will cover several aspects. Firstly, complementary high-resolution measurements of the identified SBs will help characterizing better the types of systems that compose the population of close binaries with an RG. In particular, we can infer the mass ratio for the SB2s. We will also study the possibility of leading infrared interferometric measurements to resolve the brightest and closest systems and retrieve the mass and radii of the individual components. As regards the other topics covered in this work, we will perform a detailed analysis of the HR spectra and look for chemical peculiarities of the LM-RC stars, in order to spot possible signatures of merging history. Among possible tracers, we will especially monitor the Lithium absorption line as suggested by e.g., Soares-Furtado et al. (2020), which we have already identified for a handful of targets (see Table B.1). As regards the fast rotating intermediate-mass RGs, we will investigate in more details the analysis of mixed modes to estimate the rotation rate of the core, to measure the evolution of angular momentum between the RGB and RC phases.  Fig. 15. Diagram that summarizes the sample we analyze in this paper. The hatched areas indicate the fraction of spectroscopic binaries among a given population. They represent about 85 % of the non-oscillating active RGs, 15 % of the active RGs with partially suppressed oscillations, and 2 % of the inactive RGs according to our spectroscopic sample. We note that the percentages in the active sample slightly differ from those reported in Table A.1 because they refer to the 370 active stars (with and without detected oscillations) in the present diagram, whereas they refer to the 340 RGs with oscillations in Table A.1.
M. Vrard was supported by FCT -Fundação para a Ciência e a Tecnologia through national funds and by FEDER through COMPETE2020 -Programa Operacional Competitividade e Internacionalização in the context of the grants: It is important to know what to expect in terms of RV measurements. From the studies of Rawls et al. (2016), Gaulme et al. (2016) and Benbakoura et al. (submitted), surface activity of RGs in EBs was observed with pseudo-periods ranging from 20 to 180 days, for systems with orbital periods ranging from about 15 to 120 days. The large majority of them are composed of an RG on the RGB, and an MS star (from M to F types) on loweccentricity orbits (e 0.1). The RG masses M 1 range from 0.9 to 1.6 M with typical mass ratios q = M 2 /M 1 of 0.7 to almost 1. We note the presence of an outlier, KIC 9246715, which is an an eccentric system (e = 0.36) composed of a pair of 2.1M RG stars, likely to be on the horizontal branch (Rawls et al. 2016). Another outlier to be noticed is KIC 8702921 which displays q = 0.16. For our simulations, we therefore consider a range of systems with orbital periods from 10 days to 200 days, individual masses from 1 to 2.5M , mass ratios from and 0.1 to 1, and eccentricities from 0 to 0.4. To complete our simulations we consider the archetypal system to have a M 1 = 1.2M , M 2 = 1.0M , P orb = 40 days, and e = 0. We remind that the RV of the RG component is written as: where γ 1 is the barycentric RV, K 1 the RV semi-amplitude, ω the argument of periastron, V the true anomaly, and e the eccentricity. The semi-amplitude K 1 is a function of the stellar masses and orbital properties: For circular orbits (e = 0), the RV signal is a sine curve with an amplitude K 1 . Figure A.1 shows the effect of the total mass and mass ratio of a given system on the semi-amplitude K 1 of the RV variations of the RG component. In the ranges of systems that we consider, it appears that the typical value of K 1 is ≈ 30 km s −1 , while it can range from about 3 to 90 km s −1 in some more extreme cases. Figure A.2 shows the impact of eccentricity on the values of K 1 . It appears that it has negligible impact when we consider systems with eccentricities lower than 0.4.
The inclination of the orbital plane i determines the amplitude of the actual RV signal. By assuming a uniform distribution of inclination angles, we note that 90 % of the sin i are larger than 0.16. For example, if a system has a RV semi amplitude peaking at 30 km s −1 if observed in its orbital plane, it will display an RV signal peaking more than 4.7 km s −1 in 90 % of the cases, when observed from a random location. It is thus rare than inclination makes an RV signal totally disappear.
The sampling of the measurements is the last factor to determine the likelihood of detecting an SB. We made a simulation of the RV signal of a system with zero eccentricity and a semi amplitude K. We consider the observations happening over a period of 200 days, with at random measurements dates provided that they are separated by at least 7 days from one another. We ran 100, 000 of simulated RV curves with random measurements, by considering either 2, 3, 4, 5 or 10 measurements. As expected, two measurements are statistically not sufficient to detect RV scatter. These simulations led us to observe between 4 and 7 times our targets to assess or not RV variability.

Appendix B: Search for contaminants
We have investigated the presence of nearby stars with the targets we monitored with the 3.5-m telescope at APO to flag possible false positive active RGs. We downloaded the target pixel files from the Kepler MAST and plotted the light curve of each individual pixel for each target to see if the rotational modulation was leaking from a nearby star (see details in Gaulme et al. 2013).
Most of them look fine, being the only target in the field of view. However, we note possible contamination for the following. There are two cases of two stars with similar magnitudes that are nearby in the field of view, but still separable: KIC 3458643 (10.2 and 10.8 magnitude possible RGs), and KIC 12266731.
Besides, three targets are inseparable from a nearby star (1-2 pixel difference in position), but the apparent companions are faint enough to not likely matter. GAIA also lists the brighter star in both cases as a giant. These targets are KIC 5166899, where a 10.8 mag RG overlaps with a 13.8 magnitude nearby star, KIC 8326469 (8.3 mag RG and 12.7 magnitude apparent companion), and KIC 11037219 (10.2 mag RG and 14.1 magnitude apparent Article number, page 15 of 18 A&A proofs: manuscript no. main_V6_revised companion). Another case of stars that lie within 1 or 2 pixels is KIC 5112741: it is composed of a 12.4 and a 12.6 mag stars. However, according to the pixel files, the activity signal appears to come from the target. The most suspicious targets with possible contaminants are: KIC 2833697 (and nearby KIC 2833701) with 9.3 and 9.5 magnitudes respectively, both listed as RGs by GAIA data, and KIC 8160175 with a 10.9 and a 12.2 magnitude star. In the latter case, the photometric signal seem to come from the fainter star. These observations about possible contamination are reported in the "comment" column of Table B.1.  Tayar et al. (2015) and Ceillier et al. (2017). The "known status" column indicates the number and fraction of RGs whose evolutionary status could be determined. The RGB and RC columns report the number and fractions of these types of RGs for each mass range, and the "l = 1 depleted" column indicate the number and fraction of RGs with depleted dipolar modes, according to the Mosser et al. (2017) criterion. The terms "active" and "inactive" refer to RGs with or without detected rotational modulation. The label "frac" indicates the fraction of active RGBs and RCs among all RGBs and RCs of the considered mass range. These fractions take into account the fraction of known evolutionary status. Most of the content of this table (not the l = 1 depleted RGs) can be visualized in Fig. 10 Tayar et al. (2015) and Ceillier et al. (2017) 0.5 < M < 1.   .
Evolutionary status, comments, and literature information are reported in the last three columns. When errors are available, the least significant digits in brackets after the value indicates the uncertainty. Among the comments, the abbreviation "BF x km s −1 " indicates the width of the broadening function (Rucinski 2002), which is used to point out fast rotators. The abbreviation "Li line" or "