Seismic Signatures of Stellar Magnetic Activity—What Can We Expect From TESS?

1. Introduction

The primary seismic signature of stellar activity and activity cycles is a systematic change in mode frequencies over timescales associated with activity cycles, which are typically of the order of months to decades (Baliunas et al., 1995; Saar and Brandenburg, 1999; Vida et al., 2014). Measuring this phenomenon for stars is primarily limited by the baseline length of the observations. For the Sun, however, we have Sun-as-a-star helioseismic observations lasting for decades, covering several 11-year solar activity cycles. From these data it is found that the frequency of low-degree modes close to the frequency of maximum oscillation power ν_max change by approximately 0.4 μHz between activity minimum and maximum (e.g., Jiménez-Reyes et al., 1998; Broomhall, 2017). Furthermore, this shift is strongly correlated with the level of activity along the cycle (Jiménez-Reyes et al., 1998; Broomhall and Nakariakov, 2015).

The first asteroseismic evidence of variations in stellar activity were found in an F5 V star observed by CoRoT (Baglin et al., 2006; Auvergne et al., 2009), HD 49933. For this star García et al. (2010) found that not only did the frequencies vary systematically with time but, as is also observed for the Sun (Pallé et al., 1990; Komm et al., 2000; Broomhall and Nakariakov, 2015), the amplitudes of the modes varied in anti-phase with the frequencies. In a more detailed analysis, Salabert et al. (2011) found that the observed frequency shifts are larger than seen in the Sun but the same dependence on mode frequency is observed, i.e., the magnitude of the frequency shift increases with mode frequency (e.g., Libbrecht and Woodard, 1990; Salabert et al., 2004).

Since then examples of systematic variations in the frequencies of stars observed by Kepler (Borucki et al., 2010; Koch et al., 2010) have been found (Régulo et al., 2016; Salabert et al., 2016b, 2018; Kiefer et al., 2017; Santos et al., 2018). Indeed, Kiefer et al. (2017) found that 23 out of the 24 stars in their sample showed significant and systematic frequency shifts in time, while Santos et al. (2018) found quasi-periodic variations on 60% of their sample of 87 solar-type stars. The sample of potential candidates for observing asteroseismic activity cycles is limited by the fact that activity suppresses the amplitude of acoustic oscillations (Chaplin et al., 2011). Since activity also increases with decreasing rotation period (e.g., Böhm-Vitense, 2007; Mamajek and Hillenbrand, 2008; Suárez Mascareño et al., 2016; Metcalfe and van Saders, 2017) these stars tend to be more active and consequently are less likely to have high signal-to-noise oscillations. The dependence of the length of stellar activity cycles on surface rotation rate or on Rossby number is complex and is actively researched (e.g., Strugarek et al., 2017, 2018; Warnecke, 2018). While it has proven to be difficult to observe complete activity cycles with Kepler, it is still possible to search for systematic variations with time in the properties of p modes. A concise review on the capabilities and prospects of asteroseismology regarding the inference on stellar cycles and activity is also given by Chaplin and Basu (2014).

Stellar magnetic activity also manifests in variability of stellar luminosity (see e.g., Fabbian et al., 2017, and references therein). Over the solar cycle the total irradiance is found to change by about 0.1% in temporal correlation with solar activity (Fröhlich, 2013). This slight increase in luminosity is associated with faculae (Fröhlich, 2013). In contrast to faculae, spots decrease solar and stellar luminosities (Fröhlich, 2013; Shapiro et al., 2017). For the Sun, the increase in luminosity caused by faculae dominates over the decrease in luminosity brought about by sunspots, hence the overall increase in total irradiance between cycle minimum and maximum. However, for other stars it is the decrease in luminosity caused by the presence of spots that is the dominant source of variability (Shapiro et al., 2014). Magnetic features have typical life times which are of the order of the rotation period of a star. This has to be taken into account when a meaningful proxy of stellar magnetic activity is to be constructed from the variability of the stellar luminosity. Indeed, Mathur et al. (2014) defined such a photometric proxy for stellar activity with the S_ph index. It is given by taking the mean of the standard deviations of segments of the time series each having a length 5 times the rotation period of the star (see also Salabert et al., 2016a, 2017). As long cadence photometric time series (CoRoT: 512 s, Kepler: 29.4 min, TESS: 30 min) are sufficient for such studies, the number of available light curves for the investigation of stellar magnetic activity through photometry is rather large (CoRoT: ~163,000 targets, Kepler: ~197,000 targets, TESS: > 10⁷ targets at end of nominal mission). Studies in which the variability of light curves is explored regarding signatures of activity or activity cycles include, but are not limited to Vida et al. (2014), Ferreira Lopes et al. (2015), Salabert et al. (2016a), Montet et al. (2017), and Karoff et al. (2019). In this article, we do not consider the variability of light curves, as activity-related changes are not included in the synthetic time series we use (see Ball et al., 2018 and section 2).

With the launch of NASA's Transiting Exoplanet Survey Satellite (TESS, Ricker et al., 2015), the opportunity of studying seismic signatures of many more stars is at hand. So far, p-mode frequency shifts have only been detected for a few dozen Kepler targets and one CoRoT target, and as of yet, the dependence of the amplitude of the activity-related frequency shifts (which corresponds to the amplitude of stellar activity cycles) on fundamental stellar parameters is not known. With the large number of stars that will be observed by TESS, asteroseismology can become an important additional tool to learn about stellar cycles. As we will discuss in section 6, a mission extension of TESS to 4 or 6 years is however necessary for this. This paper is aimed at a prediction of how many detections of activity-related frequency shifts we can expect in an extended TESS mission.

2. Synthetic TESS Sample

To simulate the yield of seismic signatures of activity that can be expected from TESS, we use the synthetic catalog of Ball et al. (2018), which consists of light curves that realistically mimic the properties of the TESS targets in short-cadence data (120 s)¹. The sample consists of 12,731 stars, covering main-sequence, subgiant, and low-luminosity red giant stars, which were selected from a synthetic Milky Way population simulated with TRILEGAL (Girardi et al., 2005). The selection criteria are the same as for the real Asteroseismic Target List (ATL; see Schofield et al., 2019) of the TESS Asteroseismic Science Consortium (TASC), which are chosen such as to maximize the yield of detected oscillations and optimize the coverage of different stages of stellar evolution. Ball et al. (2018) give a detailed description of their methods to obtain oscillation mode parameters (frequencies, mode damping widths, amplitudes), noise levels, granulation background, and other aspects of the simulated time series. We note that the time series of Ball et al. (2018) are provided without noise, which has to be added as ${N}(0,{\sigma }_{\mathrm{\text{noise}}})·\sqrt{30}$ to the time series. Here, σ_noise is the noise amplitude as given in the catalog of Ball et al. (2018) and the factor $\sqrt{30}$ is accounting for the 2-min sampling of the time series.

Figure 1 shows a Kiel diagram of the synthetic sample with stellar ages given by the color of the dots. We separated the sample into main-sequence stars (in the following called region I) and post-main-sequence stars (region II) along the Terminal Age Main-Sequence (TAMS), which we defined at a core hydrogen abundance content of 10⁻⁵. In Supplementary Material, Figure S1 shows the Kiel diagram as in Figure 1 but color coded for the two regions. As can be seen in these figures, the number of stars on the lower main-sequence is rather small. This is due to the low intrinsic mode amplitudes of solar-like oscillations in cool dwarf stars (see e.g., Kjeldsen and Bedding, 1995). Thus, these stars have a low detection probability for their oscillation modes with TESS data (Campante et al., 2016). A large number of targets are young, early subgiants and stars toward the end of their main-sequence life with 6, 000 K ≲ T_eff ≲ 6, 900 K. Subgiants ascending to the red giant branch and low-luminosity red giants have large mode amplitudes and high detection probabilities, hence, many of the best asteroseismic targets are in this region of the Kiel diagram.

FIGURE 1

Figure 1. Kiel diagram of the synthetic TESS sample of Ball et al. (2018). The age of the stars is given by the color of the dots. The Sun is indicated by its usual symbol. Evolutionary tracks for stars with masses from 0.8 to 2.0 M_⊙ in steps of 0.2 M_⊙ are overlaid as solid black lines.

For a study of the temporal evolution of seismic parameters, it is necessary that the targets are observed for either a long, uninterrupted period of time or that they are observed again after some time. During its nominal 2 year mission, TESS will observe 13 sectors per hemisphere with an area of four times 24° × 24° for a duration of about 27.4 days per sector. The sectors partly overlap, chiefly in the region of the ecliptic poles. Thus, some stars will have up to 13 sectors' worth of data. A detailed description of the TESS observing strategy can be found in Ricker et al. (2015). The time series of Ball et al. (2018) reproduce this observing strategy, i.e., their length depends on the position of the respective star in the sky. Here, we shall assume that the TESS mission is extended after its 2 year nominal duration. We further assume that the observation strategy of the nominal mission will be repeated in possible 2 and 4 year mission extensions, and thus that the same targets are observed again for the same number of sectors after 2 and 4 years. However, the scaling relations for p-mode frequency shifts described below are not TESS-specific.

3. Reconsidering the Scaling Relation for P-mode Frequency Shifts

There are two models for activity related frequency shifts in the literature: Chaplin et al. (2007a) simply use the amplitude of the activity cycle $\Delta R_{\mathrm{\text{HK}}}^{\prime }$ (see Equation 1) as a proxy for the frequency shift and then scale the values such that the solar frequency shift equals the observed 0.4 μHz. Saar and Brandenburg (2002) found that $\Delta R_{\mathrm{\text{HK}}}^{\prime }$ scales with the average logarithmic fraction of the stellar luminosity that is emitted in the Calcium II H and K spectral lines $R_{\mathrm{\text{HK}}}^{\prime }$:

$\begin{array}{ll}\Delta R_{\text{HK}}^{\prime }\propto {R_{\text{HK}}^{\prime }}^{0.77}.& (1)\end{array}$

For the second model, Metcalfe et al. (2007) find that the frequency shifts δν scale as

$\begin{array}{ll}\delta \nu \propto \frac{D}{I}\Delta R_{\text{HK}}^{\prime },& (2)\end{array}$

where I is the inertia of the modes and D is the depth of the perturbation that causes the mode frequencies to shift. They further assume that the depth D scales as the pressure scale height H_p, which is proportional to L^0.25R^1.5M⁻¹ at the stellar photosphere, and that mode inertia scales as MR⁻¹, where R, M, and L are stellar radius, mass, and luminosity, respectively. Thus, Metcalfe et al. (2007) find

$\begin{array}{ll}\delta \nu \propto \frac{R^{2.5}{L}^{0.25}}{M^2}\Delta R_{\text{HK}}^{\prime }& (3)\end{array}$

for the frequency shift amplitude over a full activity cycle.

Here we take a new look at the sensitivity of acoustic mode frequencies to activity. As a starting point we use Equation (1) of Metcalfe et al. (2007) to estimate the frequency shift δν:

$\begin{array}{ll}\delta \nu =\frac{\int dVKS}{2I\nu },& (4)\end{array}$

where K is a kernel for the sensitivity of p modes to perturbations, S is a source function, and ν is the frequency of the mode. Following Metcalfe et al. (2007), we use

$\begin{array}{ll}K\propto |\nabla ·\overrightarrow{\xi }{|}^2,& (5)\end{array}$

$\begin{array}{ll}S=A\delta (D-D_c),& (6)\end{array}$

where $\overrightarrow{\xi }$ is the eigenfunction of the considered mode, A is the amplitude of the source function, and the δ function determines the location of the source D_c. For now, we set A = 1 and reintroduce it later via the $R_{\mathrm{\text{HK}}}^{\prime }$ index.

In the following, we concentrate on radial modes, which simplifies the calculations considerably. The eigenfunctions for radial modes can be described by

$\begin{array}{ll}\overrightarrow{\xi }={\xi }_r(r){\overrightarrow{e}}_r,& (7)\end{array}$

where ξ_r(r) is the radial displacement eigenfunction and ${\overrightarrow{e}}_r$ is the unit vector in the radial direction. We do not introduce the spherical harmonic $Y_0^0$ into the eigenfunction here, as it would only appear as a constant in the following equations.

With this, the kernel function can be written as

$\begin{array}{ll}K=|\frac{2}{r}{\xi }_r(r)+\frac{\partial {\xi }_r(r)}{\partial r}{|}^2.& (8)\end{array}$

Approximating the derivative $\frac{\partial {\xi }_r(r)}{\partial r}\approx \frac{{\xi }_r(r)}{r}$ and assuming that the radial position r which we are considering is close enough to the surface so that the amplitude of the eigenfunction is similar to that at the photospheric stellar radius R we get

$\begin{array}{ll}K\approx \frac{9|{\xi }_r(R){|}^2}{R^2}.& (9)\end{array}$

Thus we have fixed the source to be at D_c = R. The mode inertia I can be approximated as follows:

$\begin{array}{ll}I={\displaystyle \int }dV\rho |\overrightarrow{\xi }{|}^2\approx R^3\overline{\rho }|{\xi }_r(R){|}^2,& (10)\end{array}$

where we used again the eigenfunction of a radial mode, replaced the radial integral ∫ dr by simply multiplying with the photospheric radius R, and replaced the density ρ with the mean density $\overline{\rho }$.

With Equations (8) and (10) we can rewrite Equation (4):

$\begin{array}{l}\delta \nu \approx \frac{\frac{9|{\xi }_r(R){|}^2}{R^2}}{2R^3\overline{\rho }|{\xi }_r(R){|}^2\nu }{\displaystyle \int }dV\\ =\text{const}·\frac{R}{M\nu }.& (11)\end{array}$

Throughout this, we assumed that the source is close enough to the photospheric radius so that D_c = R. All constant factors are absorbed into const. Equation (11) estimates the sensitivity of radial modes to a near-surface perturbation.

The strength of the perturbation and hence the magnitude of the frequency shift is imparted by the strength of the magnetic activity on the star for which a proxy is used. From the synthetic TESS sample, the $R_{\mathrm{\text{HK}}}^{\prime }$ index can be estimated with a relation found by Noyes et al. (1984) for main-sequence stars:

$\begin{array}{ll}log({\tau }_c/P_{\text{rot}})=-\left(0.324-0.400y+0.283y^2-1.325y^3\right),& (12)\end{array}$

where $y=log(R_{\mathrm{\text{HK}}}^{\prime }\times 10^5)$, τ_c is the convective turnover time, and P_rot is the star's rotation period. Both these values are taken from the models of Ball et al. (2018). The (local) convective turnover time τ_c is calculated half a mixing length αH_p/2 above the base of the outer convection zone of each star following (Gilliland, 1986; Landin et al., 2010):

$\begin{array}{ll}{\tau }_c=\frac{\alpha H_p}{v},& (13)\end{array}$

where the solar-calibrated mixing length parameter α = 1.957 is used, H_p is the pressure scale height, and v is the convective velocity.

To take into account that a large fraction of the stars in the synthetic TESS sample are not on the main-sequence (in contrast to the stars considered in Saar and Brandenburg, 2002 and Noyes et al., 1984), we modify the values of $R_{\mathrm{\text{HK}}}^{\prime }$ with three factors: First, we include a factor R⁻² to account for the increasing surface area over which magnetic activity can be distributed, thus affecting the oscillations less. Second, starting at a Rossby number of $\mathrm{\text{Ro}}=\frac{P_{\mathrm{\text{rot}}}}{{\tau }_c}=2$, we divide by $\sqrt{\mathrm{\text{Ro}}-1}$. This is to accommodate the results of Metcalfe and van Saders (2017) who found that, as stars reach a critical Rossby number of ≈ 2, a transition in the magnetic activity appears to happen: Magnetic braking weakens, rotation periods increase more slowly, and cycle periods grow longer and eventually disappear. Finally, we explicitly include the weakening of magnetic activity with the inverse square root of stellar age t_Age [Gyrs] found by Skumanich (1972):

$\begin{array}{ll}{R}_{\text{HK, mod}}^{\prime }=\frac{R_{\text{HK}}^{\prime }}{R^{2}max(\sqrt{\text{Ro}-1},1){\displaystyle t_{\text{Age}}^{0.5}}}.& (14)\end{array}$

With this and Equation (11), the expected frequency shift over a complete stellar activity cycle at the frequency of maximum power ν_max scales as:

$\begin{array}{ll}\delta \nu \propto \frac{R}{M{\nu }_{\text{max}}}\Delta R_{\text{HK},\text{ mod}}^{\prime }.& (15)\end{array}$

Figure 2 shows the expected p-mode frequency shifts for the synthetic TESS sample between the activity minima and maxima of a full stellar activity cycle. The values are scaled to a solar frequency shift of 0.4 μHz at ν_max. Stellar age is given by the color of the dots. For the Sun, we assume $R_{\mathrm{\text{HK}}}^{\prime }=-4.901$ (value from Saar and Brandenburg, 1999) and do not estimate it with Equation (12). The spike in shifts at 6,500 K is caused by a transition of how rotation rates are modeled by Ball et al. (2018), cf. their section 2.5 and their Figure 2.

FIGURE 2

Figure 2. Frequency shift amplitudes for full activity cycles for the synthetic TESS sample calculated with Equation (15). The color of the dots indicates stellar age. The Sun is indicated by its usual symbol.

4. Measuring δν

Two main methods have been developed to determine the change in p-mode frequency as a function of time. The first technique uses the cross-correlation (CC) between power spectra of different epochs of the time series (see e.g., Palle et al., 1989; Chaplin et al., 2007b; Régulo et al., 2016; Kiefer et al., 2017). The regions around the frequency of maximum oscillation amplitude ν_max of the power spectra are retained and their cross-correlation function is computed. Typically, only the region in which p-mode peaks can be identified is retained. As the expected frequency shifts between two power spectra of any epoch is only of the order of a few μHz, the complete cross-correlation function does not need to be calculated: restricting the cross-correlation function to a certain lag regime considerably speeds up the computation. The central peak of the cross-correlation function is then fitted to obtain the frequency shift δν_CC between the two power spectra. To estimate the uncertainty σ of the measured frequency shift, we employ the resampling approach described in Régulo et al. (2016).

The second technique uses the individual mode frequencies (peak-bagging: PB) obtained at different epochs. Here, the frequency shift δν_PB, of a particular mode with harmonic degree ℓ and radial order n is given by

$\begin{array}{ll}\delta {\nu }_{\text{PB}}={\nu }_{\ell ,n}(t)-\overline{{\nu }_{\ell ,n}},& (16)\end{array}$

where ν_{ℓ, n}(t) is the mode frequency obtained from a time series observed between t and t + δt, and $\overline{{\nu }_{\ell ,n}}$ is a reference frequency for that mode. Both methods (CC and PB) require the data to be split into time series of length δt, where δt is some fraction of the full baseline of observations available. Ideally, δt will be sufficiently short to resolve changes with time but long enough that the power spectrum has sufficient resolution and signal-to-noise for the modes to be observed clearly. See Santos et al. (2018) for a detailed description of the peak-bagging procedure in determining frequency shifts.

Most of the targets in the synthetic TESS sample have only few sectors' worth of data (one sector: 7,930 stars, two sectors: 2,601 stars, three sectors: 781 stars, four sectors or more: 1,419; see also in Supplementary Material, Tables S1 and S2). We thus calculate the power spectra for the complete synthetic time series of each star without separating them into shorter segments: here, δt is equal to the length of the time series.

In general, the two methods CC and PB are found to produce results in good agreement with each other but the uncertainties associated with the cross-correlation technique tend to be systematically larger than the frequency shifts determined directly from the mode frequencies (see Chaplin et al., 2007b; Santos et al., 2018 for comparisons of the frequency shift uncertainty obtained with the two methods). However, the peak-bagging method requires good resolution and signal-to-noise for the mode frequencies to be determined accurately and precisely and so there are instances where the cross-correlation method is preferable. The PB technique is also much more reliant on input parameters like initial guesses for mode frequencies and widths than the CC method. We do not apply the PB method here, but assume that the frequency shifts measured with both methods are in agreement and that the uncertainties of the PB method are smaller by a factor of two than those obtained with the CC method.

We implemented the CC technique to measure the frequency shifts in the following way: the power spectrum is calculated from the synthetic time series with an oversampling factor of 8. This power spectrum is then shifted in frequency by the amount that can be expected over the fraction of the stellar cycle that is covered in the 2 (4) years that pass between the TESS nominal mission and the second observation of any given star (see section 2). We estimate the length of the activity cycles by

$\begin{array}{ll}{P}_{\text{cyc}}\left[\text{years}\right]=0.5P_{\text{rot}}\left[\text{days}\right]·max\left(\sqrt{\text{Ro}-1},1\right).& (17)\end{array}$

We chose this rather crude estimation as a compromise between the short and long activity cycle branches that can be found in a P_cyc-P_rot diagram like that of Böhm-Vitense (2007) or Metcalfe and van Saders (2017). The lengthening of activity cycles above Ro ≈ 2 found by Metcalfe and van Saders (2017) is again modeled into this with the factor $max(\sqrt{\mathrm{\text{Ro}}-1},1)$. Thus, the input frequency shift is given by

$\begin{array}{ll}\delta {\nu }_{\text{input shift}}=\delta \nu ·min\left(\frac{2n}{P_{\text{cyc}}},1\right),& (18)\end{array}$

where the factor 2 accounts for the fact that the full cycle frequency shift is observed over a half-cycle, between activity minimum and maximum, n = 2 or 4 is the number of years between observations, depending on the length of the TESS mission extension, and δν is calculated with Equation (15). After shifting the power spectrum, we return to only the natural frequencies without oversampling. The oversampling was only introduced to be able to shift by frequencies smaller than the natural frequency resolution, which is ≈ 0.422 μHz for stars with only one sector of data.

The regions around the frequency of maximum oscillation amplitude ν_max of the unshifted and shifted power spectra are then cross-correlated. We retain a width of ±5Δν around ν_max, where Δν is the large frequency separation. Both Δν and ν_max are given in the catalog of Ball et al. (2018). To obtain the frequency shift between the two power spectra, we fit the cross-correlation function with a Lorentzian profile.

5. Influence of Assumptions on the Predicted Shift Amplitudes

There are several assumptions made and approximations applied in the derivation of the scaling relation given by Equation (15). The first assumption is that the magnitude of activity-related frequency shifts should be proportional to the strength of stellar magnetic activity. This is certainly justified by the fact that for the Sun there is a clear positive correlation between the p-mode frequency shifts and proxies for magnetic activity as shown by, e.g., Broomhall and Nakariakov (2015). Also, in a follow-up study to Santos et al. (2018), in which they investigate the dependence of the observed frequency shifts on stellar fundamental parameters, Santos et al. (submitted) find that for a sample of 30 solar-like stars, the amplitude of the frequency shifts decreases with decreasing $R_{\mathrm{\text{HK}}}^{\prime }$. Hence, the $R_{\mathrm{\text{HK}}}^{\prime }$ activity index, which is available for a relatively large number of stars (see e.g., Baliunas et al., 1995 and the database of Egeland, 2018) and is employed in the scaling relation of Chaplin et al. (2007b), is a good starting point for stars similar to the Sun. In principle, spectroscopic activity indices other than $R_{\mathrm{\text{HK}}}^{\prime }$ can be used in any of the frequency shift scaling relations we discuss here; Metcalfe et al. (2007), for example, use the Mg II emission index i_{Mg II}. However, the scaling of the amplitude of the index over a complete activity cycle, as is done for $R_{\mathrm{\text{HK}}}^{\prime }$ in Equation (1), has to be adjusted accordingly.

To account not only for the strength of the activity, as is done in the Chaplin et al. scaling, but also for the physics of the oscillation modes, Metcalfe et al. (2007) derived the sensitivity of mode frequencies to near surface perturbations. We note that the sensitivity of modes we find with Equation (11) is in fact very similar to that found by Metcalfe et al. (2007) and quoted by Karoff et al. (2009) in their Equation (14). Inserting the scaling relation of Kjeldsen and Bedding (1995) ${\nu }_{\mathrm{\text{max}}}\propto M{R}^{-1}{T}_{\mathrm{\text{eff}}}^{-0.5}$ into Equation (11) and using the Stefan-Boltzmann law $\sqrt{T_{\mathrm{\text{eff}}}}\propto {(R^{-2}L)}^{1/8}$, gives $\delta \nu \approx \frac{R^{2.75}{L}^{0.125}}{M^2}$. All quantities are assumed to be normalized to solar values in this discussion. The difference of our Equation (11) to the mode sensitivity factor of Metcalfe et al. (2007) is thus only a factor $R^{-0.25}{L}^{0.125}\propto \sqrt{T_{\mathrm{\text{eff}}}}$, which is in the range 0.91–1.14 for the sample of stars investigated here, i.e., is of order unity. Thus, the approximations we made in the derivation of Equation (11)—approximation of the radial integral by simply multiplying with the stellar radius R, approximation of the derivatives, localizing the perturbation to the surface—yield essentially the same result as the calculations of Metcalfe et al. (2007).

The three scaling relations are most easily distinguished for stars reaching the lower red giant branch, which is where the disagreement is largest. In Supplementary Material, Figures S2 and S3 we show the full-cycle frequency shifts of the synthetic TESS sample for the Chaplin et al. and Metcalfe et al. scalings, respectively. The Chaplin et al. (2007a) relation, simply using the $R_{\mathrm{\text{HK}}}^{\prime }$ estimation of Noyes et al. (1984) for all stars from main-sequence to the lower red giant branch, gives maximal full-cycle shift values of ≈ 3 μHz. With estimated cycle lengths of at least 10 years, no shifts larger than ≈ 1.2 μHz should be found in TESS data with a 2 year mission extension. This value increases to ≈ 5 μHz for the new scaling relation presented in this article and to even larger values for the scaling relation of Metcalfe et al. (2007).

As can be seen in Supplementary Material, Figure S3, the full-cycle frequency shifts obtained by employing the scaling relation of Metcalfe et al. (2007) would lead to very large values of over 50 μHz for stars reaching the lower red giant branch. This is why we looked for reasonable adjustments. As the sample of stars on which Noyes et al. (1984) based their results is exclusively comprised of main-sequence stars, we modified the $R_{\mathrm{\text{HK}}}^{\prime }$ index with the factors discussed above. This aims to capture some of the evolution that can naively be expected in stellar magnetic activity. Here, the decrease of activity with $t_{\mathrm{\text{Age}}}^{-0.5}$ found by Skumanich (1972) is again extrapolated from the main-sequence up to the lower red giant branch. The inclusion of this factor must be understood as an attempt to circumvent our ignorance of how exactly activity evolves after the main-sequence. It is completely possible that activity does indeed decrease faster than just proportional to $t_{\mathrm{\text{Age}}}^{-0.5}$ or even stops completely at a certain evolutionary stage. The same applies for the factor $max(\sqrt{\mathrm{\text{Ro}}-1},1)$.

Contemporaneous measurements of the $R_{\mathrm{\text{HK}}}^{\prime }$ index and photometric time series for asteroseismology would facilitate a better differentiation between these scaling relations. If indeed $\delta \nu \propto A·{(R_{\mathrm{\text{HK}}}^{\prime })}^B$, with A, B factors dependent on stellar fundamental parameters regulating mode sensitivity to perturbations (A) and the amplitude of activity over a full cycle (B), then measurements of frequency shifts and activity level can be used to learn about mode physics as well as stellar dynamos. We note that measurements of stellar rotation periods are essential for the prediction of the $R_{\mathrm{\text{HK}}}^{\prime }$ index with Equation (12) or with the modified expression (14) should direct spectroscopic measurements of stellar activity not be available.

The estimation of the cycle period as described by Equation (17) is obviously very crude. If cycle periods are underestimated with this, then the resulting shifts we can expect to observe decrease and our predictions for δν are an upper limit in this case. More, reliably measured cycle periods are needed to replace this estimation with a more rigorous prediction. We remark however, that Equation (17) does include the latest findings concerning the evolution of cycle periods by Metcalfe and van Saders (2017) through the lengthening of cycles with the factor $max(\sqrt{\mathrm{\text{Ro}}-1},1)$.

6. Potential of TESS

We evaluated the frequency shifts obtained from four settings: the TESS mission is either extended by 2 years or by 4 years and the frequency shifts are determined by the cross-correlation method (CC) or the peak-bagging method (PB), where the PB shift results are simply taken as those of the CC method but with half the uncertainty as described in section 4. The predicted number of detections of significant frequency shifts for different subsets of the cohort of TESS ATL short-cadence stars is given in Table 1 for a 4 year mission extension and in Supplementary Material, Table S3 for a 2 year mission extension. We define significant frequency shifts as those which are at least 1σ > 0, and within 3σ of the input frequency shift δν_{input shift}, where σ is the uncertainty of the measured frequency shift. Tightening the second condition to 2σ or 1σ would eliminate some of the more precisely but less accurately measured frequency shifts.

TABLE 1

Table 1. Predicted number of stars with significant frequency shifts for different spectral types and evolutionary stages.

We separated the spectral types according to their effective temperature with boundaries at T_eff = 5, 200 K between G and K and at T_eff = 6, 000 K between F and G. There are no K-type main-sequence stars in the sample. The separation into region I and II is discussed in section 2. Table 1 gives the number of stars for each subset.

The suppression of p-mode amplitudes by activity has not been taken into account in our simulations. However, only stars whose oscillations have a low detection probability and small intrinsic mode amplitudes are likely to be affected by this. These are mainly late-type dwarf stars whose number in the mock sample is rather small, thus the overall number of expected detections is not strongly affected by this. As can be seen in Supplementary Material, Tables S1 and S2 the fraction of stars with detectable shifts increases with the length of the time series as frequency resolution increases. However, even for stars which are observed for only one sector, it should be possible to detect shifts at a significant level, especially for stars with long-lived modes.

With the CC method we were able to measure significant frequency shifts in 171 stars in region I and in 2,275 stars in region II for a 4 year TESS mission extension. The large difference between the two regions reflects the increased amplitude and longer mode life times of solar-like oscillations for stars in region II. Also, the much smaller mode widths of cooler stars allow one to measure even smaller shifts, especially with the CC method. Figure 3 shows a Kiel diagram of the stars with significant shifts measured with the CC method for the 4 year extension. By comparison of Figure 1 with Figure 3 it can be seen that stars with detectable frequency shifts are largely found at effective temperatures below ≲6,500 K. Note the different color tables used in these two figures. For region I, the fraction of F-type stars for which we predict significant shifts is only 2.1% (141 out of 6,676 stars) and only 1.5% (24 out of 1,642) for region II F stars. This reflects the predicted smaller full cycle frequency shifts and large mode widths of F stars. The fraction of stars with predicted detectable frequency shifts increases for G stars: 55.6% (30 out of 54) for region I and 47.0% (765 out of 1,628) for region II. For the G dwarfs of region I this increase is mainly due to the smaller mode line widths compared to F dwarfs. Furthermore, for region II G-type stars the detectability of shifts is enhanced by the larger predicted full cycle shifts and the larger mode amplitudes. A large fraction of K-type region II stars are predicted to have detectable shifts with 54.4% (1,486 out of 2,731). With the smaller uncertainties which we assumed for the PB method, these numbers increase for all subsets of stars, adding up to a total of 3,804 out of 12,731 stars (29.9%) for the complete cohort. Table S3 gives the predicted number of detections for a 2 year TESS extension.

FIGURE 3

Figure 3. As Figure 1 but only stars with significant frequency shifts are shown. Shifts are measured with the cross-correlation technique and a 4 year TESS mission extension is assumed. Note that the color bar is different from Figure 1.

We also tested for the number of false positives caused by the realization noise. For this we set the input shift to zero for all stars and proceeded as described above. We found no false positives detections with the PB method for either a 2 or a 4 year extension even when ommiting the condition that the measured shift be within 3σ of the input shift value.

In Figures 4, 5 the significant frequency shifts obtained with the cross-correlation method for a 4 year TESS extension are plotted for three subsets of stars (rows from top to bottom: all stars, region I, region II) as a function of effective temperature (left column of panels) and rotation period (right column of panels). The color of the dots gives stellar age in Figure 4 and stellar mass in Figure 5. We capped the ordinate in these plots at 5 μHz. Higher frequency shifts tend to have large uncertainties and including them would only distort the presentation.

FIGURE 4

Figure 4. Significant frequency shifts of the synthetic TESS sample for: (A) all stars as a function of effective temperature; (B) all stars as a function of rotation period; (C) stars in region I as a function of effective temperature; (D) stars in region I as a function of rotation period; (E) stars in region II as a function of effective temperature; (F) stars in region II as a function of rotation period. The frequency shifts are measured with the cross-correlation method. Here, a 4 year TESS mission extension is assumed. The color of the dots indicates stellar age as given in the color bar on the right. Note that the range of the abscissa differs from panel to panel.

FIGURE 5

Figure 5. Same as Figure 4 but color coded with stellar mass.

By comparing panel (A) of Figure 4 with Figure 2 it can be seen that for older stars, as stellar cycles grow longer with stellar age, only a fraction of the predicted shift over a full activity cycle can be detected with data that are separated by only 4 years. As expected, younger and faster rotating stars tend to show larger detectable shifts. The older the stars become (green and blue data points), the smaller the shifts that can be expected to be observed after 2 years become. This can also be seen in panel (B), where the measured shifts are plotted as a function of rotation period: as stars age, their rotation period becomes longer, the level of activity decreases, and the frequency shifts become smaller. Obviously, this by design for our model for the full-cycle frequency shifts (Equation 15) with the factors $t_{\mathrm{\text{age}}}^{-0.5}$ and the decreasing rotation period with age as it is modeled in the synthetic TESS sample (cf. section 2.5 and Figure 2 of Ball et al., 2018). It is noteworthy that for more evolved stars, which have smaller mode damping widths, shifts are detected more precisely and smaller shifts are detectable at a significant level. The lack of stars with a large frequency shift at short rotation periods in panel (B) of Figures 4, 5 is due to the fact that most of these young, faster-rotating stars in the synthetic TESS sample are hot stars (T_eff ≳ 6, 500 K) with broader p-mode peaks for which shifts are harder to measure precisely. Also, the full cycle frequency shifts of these stars are rather small as can be seen from Figure 2.

The discrete structure in the measured frequency shifts in Figures 4, 5 is due to the large number of stars with only one sector of observation (7,930 stars). The frequency resolution of one sector's worth of observation is 0.422 μHz. As the cross-crorrelation function is computed from periodograms with this resolution, the fitted shift value is likely to be a multiple of this value.

For the stars in region I (panels C and D) there is no clear tendency in frequency shift with effective temperature or rotation period. This is partly due to the way that rotation periods are modeled in the synthetic sample. For stars with T_eff > 6, 500 K, most of which are in region I, the rotation period is drawn from a normal distribution with mean 5 days and a standard deviation of 2.1 days, i.e., there is no dependency on age or effective temperature. Overall, only a small percentage of targets with T_eff > 6, 500 K have predicted detections. The same explanation for this as given above applies. Also, stellar mass introduces a selection bias here: due to their shorter evolutionary time scales, the more massive stars in this sample do not have the time to brake their rotation and thus stars with M ≳ 1.8 M_⊙ for which detection of frequency shifts are predicted, are all at rotation periods P_rot ≲ 20 d for stars in region I (panel D) and at P_rot ≲ 40 d for stars in region II (panel F). Similarly, only lower mass stars evolve slowly enough for rotation period to increase and the activity to decrease: note the accumulation of old, low-mass stars at T_eff ≈ 5, 000 K in panels (E), which corresponds to the base of the red giant branch seen in Figure 1.

6.1. Comparison to Real Observations

To test the agreement between our predictions and real observations we concentrate on main-sequence stars as there are currently no measured frequency shifts of more evolved stars. We use the results of Santos et al. (2018) available online at this link. This is the largest available record on activity related frequency shifts to date. Santos et al. (2018) employ a PB approach on segments of Kepler time series with a length of 90 days which overlap by 45 days. We calculate the frequency shift for each star in their sample by taking the difference between the minimum and maximum value of the mean frequency shift, which combines their results for modes with harmonic degree l = 0 and 1.

The left panel of Figure 6 shows the shifts of Santos et al. (2018) in blue and the shifts of the region I stars from this work in red as a function of effective temperature. As can be seen, our predictions cover essentially the same range of shifts as the measurements. The baseline of the Kepler data used by Santos et al. (2018) is typically around 3 years, which is shorter than the baseline we assumed for our simulations. Thus, when assuming that the Kepler data do not cover a full activity cycle, extrapolating onto the same baseline would give somewhat larger shift values for the measured frequency shifts compared to our predictions. However, our sample is dominated by stars closer to the TAMS as can be seen in the right hand panel of Figure 6, which shows a Kiel diagram of the stars of Santos et al. (2018) in blue and our region I stars in red. This biases the predicted shifts more toward lower values, as these stars tend to have slowed down their rotation compared to their younger counterparts of similar mass and effective temperature. Restricting our sample further would however leave us with only a very small number of stars. We will have to wait for data of a hopefully extended TESS mission to be able to truly test our scaling relation Equation (15) and those of Chaplin et al. (2007a) and Metcalfe et al. (2007) on a larger sample of stars.

FIGURE 6

Figure 6. (Left): Significant frequency shifts from region I of the synthetic TESS data as function of effective temperature in red. A 4 year extension of the TESS mission is assumed and the error bars of the CC method are halved to approximate the PB method. Significant frequency shifts found by Santos et al. (2018) are shown in blue. (Right): Kiel diagram of the stars which are included in the left panel.

7. Overestimation of Global Stellar Parameters

The full-cycle frequency shifts as shown in Figure 2 have typical values of 6 μHz for a low-luminosity red giant. For such a star, let ν_max = 324 μHz, Δν = 23.3 μHz, T_eff = 5, 058 K, M = 1.064M_⊙ (target 00197 in the sample of Ball et al., 2018). Using the scaling relation (see e.g., Kjeldsen and Bedding, 1995)

$\begin{array}{ll}M=M_{\odot }{\left(\frac{{\nu }_{\text{max}}}{{\nu }_{\text{max,}\odot }}\right)}^3{\left(\frac{\Delta {\nu }_{\odot }}{\Delta \nu }\right)}^4{\left(\frac{T_{\text{eff}}}{T_{\text{eff,}\odot }}\right)}^{\frac{3}{2}},& (19)\end{array}$

the mass value given above is recovered. The solar reference values are those used in Campante et al. (2016), which are ν_{max, ⊙} = 3, 090 μHz, Δν_⊙ = 135.1 μHz, T_{eff, ⊙} = 5, 777 K. If ν_max is shifted by 6 μHz to higher frequencies, a mass of M_shift = 1.12M_⊙ is obtained from the scaling relation Thus, if ν_max were measured during times of higher activity, stellar mass would be overestimated by 0.06 M_⊙ with this scaling relation. The statistical uncertainties of asteroseismic stellar masses obtained from the standard scaling relations are typically of the order of 10% (Chaplin et al., 2014; Guggenberger et al., 2016). Thus, an activity-related shift to higher frequencies would introduce a systematic error and lead to an overestimation of mass and radius when using the standard scaling relations. Indeed, Gaulme et al. (2016) found in their study of 10 eclipsing binary systems with at least one oscillating red giant that the dynamical masses were consistently smaller than those obtained with the asteroseismic scaling relations. The same arguments apply for stellar radii. We note that for more evolved red giants, a shift to higher frequencies of only 1 μHz would lead to an even larger overestimation: Considering the values of Gaulme et al. (2016) for the star KIC 5786154 (ν_max = 29.75 μHz, Δν = 3.523 μHz, T_eff = 4, 747 K) and using the solar reference values and scaling relation as given above, we find a mass for this star of M = 1.44M_⊙. With ν_max = 30.75 μHz, which incorporates the frequency shift predicted here, the scaling relation gives M = 1.58M_⊙. This example shows that it might be necessary to calibrate the global seismic parameters ν_max and Δν for the level of activity of the star when the scaling relations are to be used for an estimation of mass and radius of more evolved stars. A paper on activity-related frequency shifts in evolved Kepler stars is currently in preparation. As we will show there, frequency shifts of a few tenths of μHz are not uncommon on the red giant branch.

8. Discussion and Conclusion

To estimate the yield of detections of activity-related p-mode frequency shifts in an extended TESS mission, we first derived a scaling relation for full-cycle shifts (Equation 15) for stars from the main-sequence up to low luminosity red giants. We used a catalog of simulated light curves (Ball et al., 2018) and input frequency shifts which can be expected to occur within the 2 (4) years which would pass between a first observation of a given star and a repeated observation by TESS. We then measured these frequency shifts with a cross-correlation technique and found that we can expect to find a couple hundred main-sequence and early subgiant stars and a few thousand more evolved subgiant and low-luminosity red giant stars with activity-related p-mode frequency shifts in an extended TESS mission according to our scaling relation.

Our scaling relation does not account for the varying sensitivity of modes of different harmonic degree: for the Sun it is observed that modes of higher harmonic degree experience larger shifts over the solar cycle (e.g., Broomhall, 2017). For several Kepler stars, larger shifts have been measured for l = 1 modes than for l = 0 modes by Santos et al. (2018). A dedicated study of how the frequency shifts of different harmonic degrees can influence the estimation of fundamental stellar parameters is presented in this collection by Pérez Hernández et al. (2019).

Up until now there is no confirmed detection of frequency shifts in subgiant and low-luminosity red giant stars. This may either be because their activity cycles are too long to be detectable with the available time series of Kepler and CoRoT, they no longer have activity cycles, or nobody has looked for this phenomenon in these stars yet. We also note that it is entirely possible that the scaling relation we have derived here is not valid for stars that have evolved far off the main-sequence even with the adjustments we have included to account for suppression of activity, see section 5.

A dedicated search for activity-related frequency shifts (or frequency offsets in the case of very long cycles, which could also be attributed to activity see section 7), should therefore be carried out with the existing Kepler data sets. Such results can then also be used to refine the estimates we presented here. Given our scaling relation, Equation (15) and the large number of confirmed oscillating subgiants and red giants (e.g., Huber et al., 2014; Yu et al., 2018) such a study could also help to decide between the scaling relations for the p-mode frequency shifts by Chaplin et al. (2007a), Metcalfe et al. (2007), and the one presented here.

So far, only a few dozen solar-like stars, most of them Kepler targets, have measured p-mode frequency shifts. An extended sample of seismic detections of stellar activity and activity cycles will also be useful to determine the fundamental parameters which govern stellar dynamos and their evolution as stars age. Benomar et al. (2018) have shown that it is possible to measure stellar differential rotation with asteroseismology. This can be connected with photometric measurements of activity and the temporal variation of p-mode frequencies to learn more about the intricate connection of rotation, its change through stellar evolution, and stellar activity cycles.

Carrying out searches for and investigations of stellar magnetic activity with seismic measures (mode frequencies, and amplitudes) are greatly enhanced in their reliability by contemporaneous spectroscopic measurements of magnetic activity. As, e.g., Karoff et al. (2018) have shown, it is feasible to employ seismology together with ground-based observation to study stellar cycles and activity in detail. Stars near the ecliptic poles, which are observed almost continuously during TESS's observation (615 stars in the synthetic sample have 13 sectors of data; a similar number is expected for the real observations) can be used for equivalent studies to fill the P_rot-P_cyc diagram. Long-term ground-based spectroscopic measurements of the level of activity of these targets would be expedient for such studies. The rationale to optimally select targets for such campaigns was laid out by Karoff et al. (2009) and Karoff et al. (2013). Furthermore, the investigation of stellar variability can strengthen the seismic detection of magnetic activity as, e.g., García et al. (2010) and Salabert et al. (2016b) have shown.

Our predictions for frequency shifts measured after 4 years for stars on the main-sequence and early subgiant stars (region I) agree qualitatively with those found by Kiefer et al. (2017) and Santos et al. (submitted) for Kepler stars of similar evolutionary states. The decrease of shifts with decreasing effective temperature seen by Santos et al. (submitted) is less clear in our simulations. This is most likely due to a difference in the underlying sample of stars. We will have to await TESS data from a hopefully extended mission, calculate models of the observed stars, measure the frequency shifts, and then compare predictions from our scaling relation and those of Chaplin et al. (2007a) and Metcalfe et al. (2007) to these observed shifts. As also mentioned above, measurements of a spectroscopic activity proxy like $R_{\mathrm{\text{HK}}}^{\prime }$ and reliable measurements of the stellar rotation periods and age would facilitate this considerably.

We found that it should be possible to reliably measure even small average frequency shifts between two TESS observations each with a length of only 1 month. This should also enable the search for photometric and asteroseismic signatures of short activity cycles in those young, fast rotating stars which have many sectors of data in TESS observations (see Vida et al., 2014; Reinhold et al., 2017 for such studies with Kepler data). Obviously, long and uninterrupted photometric time series with high stability are optimal for asteroseismology in general and specifically for the search for seismic signatures of activity cycles. Even though most of the targets of TESS only have few consecutive sectors of observation, time series provided by TESS, especially by an extended mission, can be used for such searches. In the future, PLATO (Rauer et al., 2014) will pick up the torch, observe two target fields for 2–3 years each, and provide data for asteroseismology. These data will make another major contribution in the seismic sounding of stellar activity cycles.

Author Contributions

RK led the work on this article and wrote most of the text. A-MB contributed to the writing of sections 1, 4, and the overall refinement of the writing. WB generated and provided the dataset which is utilized in this article, computed additional necessary quantities, and helped in the writing of the relevant sections.

Funding

This work was funded by STFC consolidated grant ST/L000733/1. WB acknowledges support from the UK Science and Technology Facilities Council (STFC).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

RK and A-MB acknowledge the support of the STFC consolidated grant ST/L000733/1.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2019.00052/full#supplementary-material

Footnotes

1. ^The datasets used for this study can be found in this repository: https://zenodo.org/record/1470155#.XTpX4jfhXIU.

References

Auvergne, M., Bodin, P., Boisnard, L., Buey, J.-T., Chaintreuil, S., Epstein, G., et al. (2009). The CoRoT satellite in flight: description and performance. Astronomy Astrophy. 506, 411–424. doi: 10.1051/0004-6361/200810860

CrossRef Full Text | Google Scholar

Baglin, A., Auvergne, M., Barge, P., Deleuil, M., Catala, C., Michel, E., et al. (2006). “Scientific Objectives for a Minisat: CoRoT,” in The CoRoT Mission Pre-Launch Status - Stellar Seismology and Planet Finding, Vol 1306, eds M. Fridlund, A. Baglin, J. Lochard, and L. Conroy (Noordwijk: ESA Special Publication), 33.

Baliunas, S. L., Donahue, R. A., Soon, W. H., Horne, J. H., Frazer, J., Woodard-Eklund, L., et al. (1995). Chromospheric variations in main-sequence stars. Astrophy. J. 438, 269–287.

Google Scholar

Ball, W. H., Chaplin, W. J., Schofield, M., Miglio, A., Bossini, D., Davies, G. R., and Girardi, L. (2018). A synthetic sample of short-cadence solar-like oscillators for TESS. Astrophys. J. Suppl. Series 239:34. doi: 10.3847/1538-4365/aaedbc

CrossRef Full Text | Google Scholar

Benomar, O., Bazot, M., Nielsen, M. B., Gizon, L., Sekii, T., Takata, M., et al. (2018). Asteroseismic detection of latitudinal differential rotation in 13 Sun-like stars. Science 361, 1231–1234. doi: 10.1126/science.aao6571

PubMed Abstract | CrossRef Full Text | Google Scholar

Böhm-Vitense, E. (2007). Chromospheric activity in G and K main-sequence stars, and what it tells us about stellar dynamos. Astrophys. J. 657, 486–493. doi: 10.1086/510482

CrossRef Full Text | Google Scholar

Borucki, W. J., Koch, D., Basri, G., Batalha, N., Brown, T., Caldwell, D., et al. (2010). Kepler planet-detection mission: introduction and first results. Science 327, 977–980. doi: 10.1126/science.1185402

PubMed Abstract | CrossRef Full Text | Google Scholar

Broomhall, A.-M. (2017). A helioseismic perspective on the depth of the minimum between solar cycles 23 and 24. Solar Phys. 292:67. doi: 10.1007/s11207-017-1068-5

CrossRef Full Text | Google Scholar

Broomhall, A.-M., and Nakariakov, V. M. (2015). A comparison between global proxies of the Sun's magnetic activity cycle: inferences from helioseismology. Solar Phys. 290, 3095–3111. doi: 10.1007/s11207-015-0728-6

CrossRef Full Text | Google Scholar

Campante, T. L., Schofield, M., Kuszlewicz, J. S., Bouma, L., Chaplin, W. J., Huber, D., et al. (2016). The asteroseismic potential of TESS: exoplanet-host stars. Astrophys. J. 830:138. doi: 10.3847/0004-637X/830/2/138

CrossRef Full Text | Google Scholar

Chaplin, W. J., and Basu, S. (2014). Inferences on stellar activity and stellar cycles from asteroseismology. Space Sci. Rev. 186, 437–456. doi: 10.1007/s11214-014-0090-2

CrossRef Full Text | Google Scholar

Chaplin, W. J., Basu, S., Huber, D., Serenelli, A., Casagrande, L., Silva Aguirre, V., et al. (2014). Asteroseismic fundamental properties of solar-type stars observed by the NASA kepler mission. Astrophys. J. Suppl. Series 210:1. doi: 10.1088/0067-0049/210/1/1

CrossRef Full Text | Google Scholar

Chaplin, W. J., Bedding, T. R., Bonanno, A., Broomhall, A.-M., García, R. A., Hekker, S., et al. (2011). Evidence for the impact of stellar activity on the detectability of solar-like oscillations observed by kepler. Astrophys. J. Lett. 732:L5. doi: 10.1088/2041-8205/732/1/L5

CrossRef Full Text | Google Scholar

Chaplin, W. J., Elsworth, Y., Houdek, G., and New, R. (2007a). On prospects for sounding activity cycles of Sun-like stars with acoustic modes. Mon. Notices R. Astron. Soc. 377, 17–29. doi: 10.1111/j.1365-2966.2007.11581.x

CrossRef Full Text | Google Scholar

Chaplin, W. J., Elsworth, Y., Miller, B. A., Verner, G. A., and New, R. (2007b). Solar p-mode frequencies over three solar cycles. Astrophys. J. 659, 1749–1760. doi: 10.1086/512543

CrossRef Full Text | Google Scholar

Egeland, R. (2018). Sunstardb: a database for the study of stellar magnetism and the solar-stellar connection. Astrophys. J. Suppl. Series 236:19. doi: 10.3847/1538-4365/aab771

CrossRef Full Text | Google Scholar

Fabbian, D., Simoniello, R., Collet, R., Criscuoli, S., Korhonen, H., Krivova, N. A., et al. (2017). The variability of magnetic activity in solar-type stars. Astronom. Nachricht. 338, 753–772. doi: 10.1002/asna.201713403

CrossRef Full Text | Google Scholar

Ferreira Lopes, C. E., Leão, I. C., de Freitas, D. B., Canto Martins, B. L., Catelan, M., and De Medeiros, J. R. (2015). Stellar cycles from photometric data: CoRoT stars. Astronomy Astrophy. 583:A134. doi: 10.1051/0004-6361/201424900

CrossRef Full Text | Google Scholar

Fröhlich, C. (2013). Total solar irradiance: what have we learned from the last three cycles and the recent minimum? Space Sci. Rev. 176, 237–252. doi: 10.1007/s11214-011-9780-1

CrossRef Full Text | Google Scholar

García, R. A., Mathur, S., Salabert, D., Ballot, J., Régulo, C., Metcalfe, T. S., et al. (2010). CoRoT reveals a magnetic activity cycle in a sun-like star. Science 329:1032. doi: 10.1126/science.1191064

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaulme, P., McKeever, J., Jackiewicz, J., Rawls, M. L., Corsaro, E., Mosser, B., et al. (2016). Testing the asteroseismic scaling relations for red giants with eclipsing binaries observed by kepler. Astrophys. J. 832:121. doi: 10.3847/0004-637X/832/2/121

CrossRef Full Text | Google Scholar

Gilliland, R. L. (1986). Relation of chromospheric activity to convection, rotation, and pre–main-sequence evolution. Astrophys. J. 300:339.

Google Scholar

Girardi, L., Groenewegen, M. A. T., Hatziminaoglou, E., and da Costa, L. (2005). Star counts in the Galaxy. Simulating from very deep to very shallow photometric surveys with the TRILEGAL code. Astronomy Astrophy. 436, 895–915. doi: 10.1051/0004-6361:20042352

CrossRef Full Text | Google Scholar

Guggenberger, E., Hekker, S., Basu, S., and Bellinger, E. (2016). Significantly improving stellar mass and radius estimates: a new reference function for the Δν scaling relation. Mon. Notices R. Astron. Soc. 460, 4277–4281. doi: 10.1093/mnras/stw1326

CrossRef Full Text | Google Scholar

Huber, D., Silva Aguirre, V., Matthews, J. M., Pinsonneault, M. H., Gaidos, E., García, R. A., et al. (2014). Revised stellar properties of kepler targets for the quarter 1-16 transit detection run. Astrophys. J. Suppl. Series 211:2. doi: 10.1088/0067-0049/211/1/2

CrossRef Full Text | Google Scholar

Jiménez-Reyes, S. J., Régulo, C., Pallé, P. L., and Roca Cortés, T. (1998). Solar activity cycle frequency shifts of low-degree p-modes. Astronomy Astrophy. 329, 1119–1124.

Google Scholar

Karoff, C., Metcalfe, T. S., Chaplin, W. J., Elsworth, Y., Kjeldsen, H., Arentoft, T., and Buzasi, D. (2009). Sounding stellar cycles with Kepler - I. Strategy for selecting targets. Mon. Notices R. Astron. Soc. 399, 914–923. doi: 10.1111/j.1365-2966.2009.15323.x

CrossRef Full Text | Google Scholar

Karoff, C., Metcalfe, T. S., Chaplin, W. J., Frandsen, S., Grundahl, F., Kjeldsen, H., et al (2013). Sounding stellar cycles with Kepler - II. Ground-based observations. Mon. Notices R. Astron. Soc. 433, 3227–3238. doi: 10.1093/mnras/stt964

CrossRef Full Text | Google Scholar

Karoff, C., Metcalfe, T. S., Montet, B. T., Jannsen, N. E., Santos, A. R. G., Nielsen, M. B., et al. (2019). Sounding stellar cycles with Kepler - III. Comparative analysis of chromospheric, photometric, and asteroseismic variability. Mon. Notices R. Astron. Soc. 485, 5096–5104. doi: 10.1093/mnras/stz782

CrossRef Full Text | Google Scholar

Karoff, C., Metcalfe, T. S., Santos, Â. R. G., Montet, B. T., Isaacson, H., Witzke, V., et al. (2018). The influence of metallicity on stellar differential rotation and magnetic activity. Astrophys. J. 852:46. doi: 10.3847/1538-4357/aaa026

CrossRef Full Text | Google Scholar

Kiefer, R., Schad, A., Davies, G., and Roth, M. (2017). Stellar magnetic activity and variability of oscillation parameters: An investigation of 24 solar-like stars observed by Kepler. Astronomy Astrophy. 598:A77. doi: 10.1051/0004-6361/201628469

CrossRef Full Text | Google Scholar

Kjeldsen, H., and Bedding, T. R. (1995). Amplitudes of stellar oscillations: the implications for asteroseismology. Astronomy Astrophy. 293, 87–106.

Google Scholar

Koch, D. G., Borucki, W. J., Basri, G., Batalha, N. M., Brown, T. M., Caldwell, D., et al. (2010). Kepler mission design, realized photometric performance, and early science. Astrophys. J. Lett. 713, L79–L86. doi: 10.1088/2041-8205/713/2/L79

CrossRef Full Text | Google Scholar

Komm, R. W., Howe, R., and Hill, F. (2000). Width and energy of solar p-modes observed by global oscillation network group. Astrophys. J. 543, 472–485. doi: 10.1086/317101

CrossRef Full Text | Google Scholar

Landin, N. R., Mendes, L. T. S., and Vaz, L. P. R. (2010). Theoretical values of convective turnover times and Rossby numbers for solar-like, pre-main sequence stars. Astronomy Astrophy. 510:A46. doi: 10.1051/0004-6361/200913015

CrossRef Full Text | Google Scholar

Libbrecht, K. G., and Woodard, M. F. (1990). Solar-cycle effects on solar oscillation frequencies. Nature 345, 779–782. doi: 10.1038/345779a0

CrossRef Full Text | Google Scholar

Mamajek, E. E., and Hillenbrand, L. A. (2008). Improved age estimation for solar-type dwarfs using activity-rotation diagnostics. Astrophys. J. 687, 1264–1293. doi: 10.1086/591785

CrossRef Full Text | Google Scholar

Mathur, S., Salabert, D., García, R. A., and Ceillier, T. (2014). Photometric magnetic-activity metrics tested with the Sun: application to Kepler M dwarfs. J. Space Weather Space Clim. 4:A15. doi: 10.1051/swsc/2014011

CrossRef Full Text | Google Scholar

Metcalfe, T. S., Dziembowski, W. A., Judge, P. G., and Snow, M. (2007). Asteroseismic signatures of stellar magnetic activity cycles. Mon. Notices R. Astron. Soc. 379, L16–L20. doi: 10.1111/j.1745-3933.2007.00325.x

CrossRef Full Text | Google Scholar

Metcalfe, T. S., and van Saders, J. (2017). Magnetic evolution and the disappearance of sun-Like activity cycles. Solar Phys. 292:126. doi: 10.1007/s11207-017-1157-5

CrossRef Full Text | Google Scholar

Montet, B. T., Tovar, G., and Foreman-Mackey, D. (2017). Long-term photometric variability in kepler full-frame images: magnetic cycles of sun-like stars. Astrophys. J. 851:116. doi: 10.3847/1538-4357/aa9e00

CrossRef Full Text | Google Scholar

Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., and Vaughan, A. H. (1984). Rotation, convection, and magnetic activity in lower main-sequence stars. Astrophys. J. 279, 763–777.

Google Scholar

Pallé, P., Régulo, C., and Roca Cortés, T. (1990). “The spectrum of solar p-modes and the solar activity cycle,” in Progress of Seismology of the Sun and Stars, Lecture Notes in Physics, Volume 367, eds Y. Osaki and H. Shibahashi (Berlin; Heidelberg: Springer), 129–134.

Google Scholar

Palle, P. L., Regulo, C., and Roca Cortes, T. (1989). Solar cycle induced variations of the low L solar acoustic spectrum. Astronomy Astrophy. 224, 253–258.

Google Scholar

Pérez Hernández, F., García, R. A., Mathur, S., Santos, A. R. G., and Régulo, C. (2019). Influence of magnetic activity on the determination of stellar parameters through asteroseismology. Frontiers Astronomy Space Sci. 6:41. doi: 10.3389/fspas.2019.00041

CrossRef Full Text | Google Scholar

Rauer, H., Catala, C., Aerts, C., Appourchaux, T., Benz, W., Brandeker, A., et al. (2014). The PLATO 2.0 mission. Exp. Astronomy 38, 249–330. doi: 10.1007/s10686-014-9383-4

CrossRef Full Text | Google Scholar

Régulo, C., García, R. A., and Ballot, J. (2016). Magnetic activity cycles in solar-like stars: the cross-correlation technique of p-mode frequency shifts. Astronomy Astrophy. 589:A103. doi: 10.1051/0004-6361/201425408

CrossRef Full Text | Google Scholar

Reinhold, T., Cameron, R. H., and Gizon, L. (2017). Evidence for photometric activity cycles in 3203 Kepler stars. Astronomy Astrophy. 603:A52. doi: 10.1051/0004-6361/201730599

CrossRef Full Text | Google Scholar

Ricker, G. R., Winn, J. N., Vanderspek, R., Latham, D. W., Bakos, G. Á., Bean, J. L., et al. (2015). Transiting exoplanet survey satellite (TESS). J. Astron. Telescopes Instr. Syst. 1:014003. doi: 10.1117/1.JATIS.1.1.014003

CrossRef Full Text | Google Scholar

Saar, S. H., and Brandenburg, A. (1999). Time evolution of the magnetic activity cycle period. II. Results for an expanded stellar sample. Astrophys. J. 524, 295–310.

Google Scholar

Saar, S. H., and Brandenburg, A. (2002). A new look at dynamo cycle amplitudes. Astronom. Nachricht. 323, 357–360. doi: 10.1002/1521-3994(200208)323:3/4<357::AID-ASNA357>3.0.CO;2-I

CrossRef Full Text | Google Scholar

Salabert, D., Fossat, E., Gelly, B., Kholikov, S., Grec, G., Lazrek, M., et al. (2004). Solar p modes in 10 years of the IRIS network. Astronomy Astrophy. 413, 1135–1142. doi: 10.1051/0004-6361:20031541

CrossRef Full Text | Google Scholar

Salabert, D., García, R. A., Beck, P. G., Egeland, R., Pallé, P. L., Mathur, S., et al. (2016a). Photospheric and chromospheric magnetic activity of seismic solar analogs. Observational inputs on the solar-stellar connection from Kepler and Hermes. Astronomy Astrophy. 596:A31. doi: 10.1051/0004-6361/201628583

CrossRef Full Text | Google Scholar

Salabert, D., García, R. A., Jiménez, A., Bertello, L., Corsaro, E., and Pallé, P. L. (2017). Photospheric activity of the Sun with VIRGO and GOLF. Comparison with standard activity proxies. Astronomy Astrophy. 608:A87. doi: 10.1051/0004-6361/201731560

CrossRef Full Text | Google Scholar

Salabert, D., Régulo, C., Ballot, J., García, R. A., and Mathur, S. (2011). About the p-mode frequency shifts in HD 49933. Astronomy Astrophy. 530:A127. doi: 10.1051/0004-6361/201116633

CrossRef Full Text | Google Scholar

Salabert, D., Régulo, C., García, R. A., Beck, P. G., Ballot, J., Creevey, O. L., et al. (2016b). Magnetic variability in the young solar analog KIC 10644253. Observations from the Kepler satellite and the HERMES spectrograph. Astronomy Astrophy. 589:A118. doi: 10.1051/0004-6361/201527978

CrossRef Full Text | Google Scholar

Salabert, D., Régulo, C., Pérez Hernández, F., and García, R. A. (2018). Frequency dependence of p-mode frequency shifts induced by magnetic activity in Kepler solar-like stars. Astronomy Astrophy. 611:A84. doi: 10.1051/0004-6361/201731714

CrossRef Full Text | Google Scholar

Santos, Â. R. G., Campante, T. L., Chaplin, W. J., Cunha, M. S., Lund, M. N., Kiefer, R., et al. (2018). Signatures of magnetic activity in the seismic data of solar-type stars observed by kepler. Astrophys. J. Suppl. Series 237:17. doi: 10.3847/1538-4365/aac9b6

CrossRef Full Text | Google Scholar

Schofield, M., Chaplin, W. J., Huber, D., Campante, T. L., Davies, G. R., Miglio, A., et al. (2019). The asteroseismic target list for solar-like oscillators observed in 2 minute cadence with the transiting exoplanet survey satellite. Astrophys. J. Suppl. Series 241:12. doi: 10.3847/1538-4365/ab04f5

CrossRef Full Text | Google Scholar

Shapiro, A. I., Solanki, S. K., Krivova, N. A., Cameron, R. H., Yeo, K. L., and Schmutz, W. K. (2017). The nature of solar brightness variations. Nat. Astronomy 1, 612–616. doi: 10.1038/s41550-017-0217-y

CrossRef Full Text | Google Scholar

Shapiro, A. I., Solanki, S. K., Krivova, N. A., Schmutz, W. K., Ball, W. T., Knaack, R., et al. (2014). Variability of Sun-like stars: reproducing observed photometric trends. Astronomy Astrophy. 569:A38. doi: 10.1051/0004-6361/201323086

CrossRef Full Text | Google Scholar

Skumanich, A. (1972). Time scales for CaII emission decay, rotational braking, and lithium depletion. Astrophys. J. 171, 565–567.

Google Scholar

Strugarek, A., Beaudoin, P., Charbonneau, P., and Brun, A. S. (2018). On the sensitivity of magnetic cycles in global simulations of solar-like stars. Astrophys. J. 863:35. doi: 10.3847/1538-4357/aacf9e

CrossRef Full Text | Google Scholar

Strugarek, A., Beaudoin, P., Charbonneau, P., Brun, A. S., and do Nascimento, J. D. (2017). Reconciling solar and stellar magnetic cycles with nonlinear dynamo simulations. Science 357, 185–187. doi: 10.1126/science.aal3999

PubMed Abstract | CrossRef Full Text | Google Scholar

Suárez Mascareño, A., Rebolo, R., and González Hernández, J. I. (2016). Magnetic cycles and rotation periods of late-type stars from photometric time series. Astronomy Astrophy. 595:A12. doi: 10.1051/0004-6361/201628586

CrossRef Full Text | Google Scholar

Vida, K., Oláh, K., and Szabó, R. (2014). Looking for activity cycles in late-type Kepler stars using time-frequency analysis. Mon. Notices R. Astron. Soc. 441, 2744–2753. doi: 10.1093/mnras/stu760

CrossRef Full Text | Google Scholar

Warnecke, J. (2018). Dynamo cycles in global convection simulations of solar-like stars. Astronomy Astrophy. 616:A72. doi: 10.1051/0004-6361/201732413

CrossRef Full Text | Google Scholar

Yu, J., Huber, D., Bedding, T. R., Stello, D., Hon, M., Murphy, S. J., et al. (2018). Asteroseismology of 16,000 kepler red giants: global oscillation parameters, masses, and radii. Astrophys. J. Suppl. Series 236:42. doi: 10.3847/1538-4365/aaaf74

CrossRef Full Text | Google Scholar