A new catalog of magnetically active solar-like oscillators

We present a new catalog of stars for which detected solar-like oscillations and magnetic activity measurements are both available from chromospheric spectroscopic observations. Our results were obtained by exploiting NASA TESS mission light curves for active stars observed within the Mount Wilson Observatory HK project and the HK survey of the Hamburg Robotic Telescope TIGRE. We analyzed the light curves for a total of 191 stars by adopting recent techniques based on Bayesian analysis and model comparison to assess the detection of a power excess originating from solar-like oscillations. We characterized the oscillations in a total of 34 targets, for which we provide estimates for the global asteroseismic parameters of $\nu_\mathrm{max}$ (the frequency of maximum oscillation power), $\Delta\nu$ (the large frequency separation), and for the amplitude of the solar-like oscillation envelope $A_\mathrm{max}$. We provide strong statistical evidence for the detection of solar-like oscillations in 15 stars of our sample, identify six further stars where a detection is likely, and 13 stars for which oscillations cannot be ruled out. The key parameters extracted in this work will be exploited for a detailed stellar modeling of the targets and to calibrate relations that connect the level of the measured magnetic activity to the suppression induced on the global oscillation amplitudes. This opens the possibility of shedding light on the interplay between magnetic fields and oscillations. Because of their relatively high brightness, the targets may also be of interest for future dedicated follow-up observations using both photometry and spectropolarimetry.


Introduction
Photometric observations of brightness variations induced by stellar oscillations have proven to be a powerful way for studying stars in detail (Aerts et al. 2010).Space missions such as CoRoT 1 (Baglin et al. 2006), NASA Kepler and K2 (Borucki et al. 2010;Howell et al. 2014), and TESS 2 (Ricker et al. 2014) have provided and are still providing a wealth of data that can be used to probe fundamental stellar properties such as mass and radius.This enables subsequent modeling aimed at resolving the internal structure and dynamics of stars and at deriving their age.Solar-like oscillations in particular have been observed in tens of thousands of stars with low to intermediate mass because, similarly to the Sun, these oscillations are excited by near-surface convection.This has allowed asteroseismology to consolidate its role over the past years as a key science in fostering our understanding of stars and their evolution.In this context, the study of magnetic fields is of crucial importance because they operate throughout the stellar evolution, starting at the onset of star formation (e.g.McKee & Ostriker 2007).As originally found from solar observations Chaplin et al. (2000), magnetic fields impact on stellar oscillations in different ways (e.g.García et al. 2010;Chaplin et al. 2011a), but their interplay is not yet fully under-Tables 1 and 2 are also available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr(130.79.128.5) or via https:// cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/683/A161 1 COnvection ROtation and planetary Transits. 2Transiting Exoplanet Survey Satellite.
stood (e.g.Bonanno et al. 2014Bonanno et al. , 2019)).The reason mainly is that it is difficult to find usable targets other than the Sun that are suitable for this type of study, because magnetic activity in stars has the main effect of suppressing the oscillation amplitude (e.g.Di Mauro et al. 2022).The vast majority of targets with detected solar-like oscillations currently are evolved stars with low levels of magnetic activity, which are clearly not suited for studies of their magnetism.In this respect, increasing the sample of stars with detected oscillations and available measurements of magnetic activity is a mandatory task that needs to be accomplished.
The TESS mission, although not comparable to Kepler in terms of photometric precision and duration of the observations, is observing a large number of relatively nearby and bright stars with its large sky coverage.These stars include some for which the magnetic activity index obtained from chromospheric spectroscopy observations of the Ca II H&K emission lines, for instance, is available.Many targets with measurements of chromspheric activity are part of the historical Mt.Wilson Observatory HK Project (e.g.Wilson 1978), to which more were added by means of the HK survey of the Hamburg Robotic Telescope (TIGRE; Mittag et al. 2011).A high fraction of these targets does not fall within the Kepler field of view of the nominal mission, and therefore, it was not possible to detect potential oscillations through photometric observations before the advent of TESS.
We exploit the known catalogs of stars with measured magnetic activity from chromospheric spectroscopic observations (Mittag et al. 2011;Marsden et al. 2014;Boro Saikia et al. 2018;Olspert et al. 2018) to analyze the corresponding TESS light curves acquired during mission Cycles 1, 2, and 3 in order to search for potential solar-like oscillations.We exploit a Bayesian framework (Corsaro & De Ridder 2014;Müllner et al. 2021) to characterize the stellar power spectra and to assess the presence of an oscillation power excess.For the stars with a potential detection, we extract the global asteroseismic properties of ν max , ∆ν, and A max using well-consolidated modeling techniques of the stellar power spectra (Kallinger et al. 2014;Corsaro et al. 2017).We compile a new catalog of stars in this way that will be used to perform subsequent stellar modeling and to investigate the relation between the magnetic activity level and suppression of the oscillation amplitude for a wider range of fundamental stellar properties than was possible before.

Observations and data
The stars used in this work have been observed by NASA TESS over the past years during Cycles 1, 2, and 3 for a typical duration of about 27 days (one sector) with a cadence of 2 min.By cross-matching the catalogs of active stars with the MAST3 data archive, we collected usable light curves for a total of 191 stars, 118 of which were taken from Mittag et al. (2011), 20 from BCool Marsden et al. (2014), 24 from Boro Saikia et al. (2018), and 31 from Olspert et al. (2018).
Using the standard PDCSAP4 flux as a starting dataset (e.g.García et al. 2011), we processed each light curve in order to optimize it for the detection of solar-like oscillations.This was done by first applying a 3σ clipping over a shifted light-curve flux with zero mean.Then, a boxcar smoothing was applied, with a varying width adapted to the individual star, and typically ranging from 400 to 800 bins, depending on the level of signal modulation at low frequency.The smoothed light curve was subtracted from the σ clipped mean-shifted light curve, and the mean flux of the smoothed light curve was added to the final flux.In addition, the obtained flux was normalized to its median value, shifted to zero, and multiplied by a factor 10 6 to match the parts-per-million (ppm) units.The procedure involving the 3σ clipping and the subtraction of a smoothed light curve was repeated iteratively three more times to improve the removal of long-trend variations that may hamper the detection of oscillations at intermediate to high frequencies due to the leakage of power caused by high-amplitude signal at low frequencies.An example of this processing is shown in Fig. 1 for the star TIC 95340781 (HD 88737), where an oscillation envelope is detected at ν max ≈ 670 µHz.The first is the most siginficant iteration (first row), allowing the light curve to be flattened out by removing the strongest modulation signal and some prominent outliers (shown in the left panel before the processing was applied, and in central panel after the processing was applied).The second to fourth iterations further increased the flattening of the light curve at shorter timescales, which could be useful to improve the signal-to-noise ratio (S/N) of the oscillation envelope.After the fourth iteration, no further evident improvements were observed, and the procedure was therefore no longer applied beyond this point.The corresponding power spectral density (PSD) is shown in the right panels of the plot, allowing us to keep under control the impact of the filtering on the data product that was subsequently analyzed.We note that the variability shown by the residual light curve on timescales of about a day or shorter typically originates in a real astrophysical signal, such as stellar granulation and oscillations, and therefore, it was intentionally not removed.For reference, the processing of the raw light curve that we described here follows a standard method that is adopted to prepare asteroseismically optimized light curves acquired by NASA Kepler (García et al. 2011).
To facilitate and guide the subsequent analysis of the PSDs, we collected the available spectroscopic information of each target by exploiting the public PASTEL (Soubiran et al. 2010) and Gaia DR3 catalogs (Gaia Collaboration 2016;Creevey et al. 2023).We considered log g and T eff of the star because they can be used to predict the frequency position of a potential oscillation envelope through the well-known relation that connects ν max to the acoustic cutoff frequency (ν ac ) (Brown et al. 1991), and by scaling this relation to the solar values ν max = 3150 µHz, T eff = 5777 K, and log g = 4.43775 dex (Chaplin et al. 2011b).
For stars with multiple measurements of log g and T eff , we simply considered an average.The predicted values for ν max have the main purpose to represent a guess, or starting point, for the analysis presented in Sect.3, meaning that we cannot expect these values to necessarily coincide with the final measurements obtained from this work.Because most of the datasets we adopted show high levels of background noise (in particular, photon noise), the reference values for ν max are also useful to guarantee that our potential detections of solar-like oscillations are compliant with the spectroscopic parameters of the star.This helps us to reduce spurious detections (e.g., caused by the contaminating signal from a background star) and to increase the robustness of our results.

Asteroseismic data analysis
The analysis of the stellar PSD was carried out in order to estimate ν max , detect potential solar-like oscillations, and when this was verified, extract the additional global asteroseismic parameters ∆ν and A max .In the following sections, we describe the details of this analysis.

Background fitting and ν max
After the PSD was obtained as described in Sect.2, it was modeled so that the level of the background signal could be characterized, namely the signal that is not related to solar-like oscillations.This is essential to be able to quantify the S/N of a potential solar-like oscillation envelope, and therefore, to test for its detection.For this purpose, we followed the representation introduced by Kallinger et al. (2014), which comprises the signal from long-trend variations that occur at low frequency (as the result of the combined effects of stellar activity, supergranulation, and other possible instrumental variations), from granulation activity (modeled as a series of Harvey-like profiles; Harvey 1985), and from photometric white noise.The region containing the solar-like oscillations was instead modeled by a typical Gaussian function.The model we describe here is represented by the following equation: where W is the instrumental noise, and the number of Harveylike components can range from one to three, depending on the complexity of the PSD (a i and b i represent the amplitude A161, page 2 of 9 Fig. 1.Data preparation process for the light curve of TIC 95340781 (HD 88737).Each row represents a different iteration (four in total), as outlined in Sect. 2. The first column shows the light curve before the filtering is applied (in orange), the second column shows the (residual) light curve after the filtering is applied (in gray), and the third column shows the PSD of the residual light curve (in blue).The solid green line shows the smoothing of the light curve at each step, and the yellow line in the third column represents the smoothing of the PSD (i.e., an approximation of the background level; see Sect. and characteristic frequency, respectively, of the ith component), and in turn on the capability of the data to constrain the model.This model clearly permits us to test the hypothesis that the oscillation envelope is present (and therefore, detected) in the given dataset.We followed the method presented by Corsaro et al. (2017) and fit the different background models for every star by means of the public Bayesian inference software Diamonds 5 (Corsaro & De Ridder 2014) and its extension, called Background, which specifically implements the above models (and others as well) for the background fitting in stellar PSDs.This code is particularly suited for this type of application because the fitting parameters are often highly correlated to one another (e.g.Corsaro & De Ridder 2014), which makes it more challenging to find a good solution in the parameter space.The high computational and sampling efficiency is made possible by the implementation of a nested-sampling Monte Carlo algorithm (Skilling 2004), which has the additional key advantage of providing an estimate of the Bayesian evidence (or marginal likelihood) as a direct output for every model that is fit to the data.The Bayesian evidence is crucial to assess the detection, as we discuss below.An example of a background fit to the PSD is shown in Fig. 2 for the stars TIC 408842743 (HD 187691) and TIC 38511251 (HD 23249), where varying background models were adopted.We chose the adequate model in a Bayesian model comparison framework, as presented by Corsaro et al. (2017).In prac-5 High-DImensional And multi-MOdal NesteD Sampling, https:// github.com/EnricoCorsaro/DIAMONDStice, this analysis was conducted by first testing the most complex model (i.e., the model with the largest number of Harveylike profiles) that can be fit by the code to the PSD of a given star.Then, the Bayesian evidence of this model is compared to that of a subsequent model that is slightly simpler (i.e., with one Harvey-like component less than the first model).We therefore selected the most complex model that was tested to determine whether it has the highest Bayesian evidence, and otherwise, we proceeded by including an additional even more simplified model and repeated the model comparison process.We note that when too many Harvey-like profiles were included with respect to what the data were capable of constraining, the code failed to converge to a result.The given model was then discarded because it cannot be reliably included in the model comparison process.In this case, a less complex model was replaced to be tested in addition to the previous model.For most stars, we find that the model accounting for two Harvey-like profiles is favored, followed by the one Harvey-like profile model, and finally, by the three Harvey-like profile model.This is expected because, on the one hand, the instrumental noise level is often so high that the signal in the PSD appears to be nearly flat for most of the observed frequency range (up to the Nyquist frequency), with a power level increasing only toward the low-frequency side.On the other hand, the presence of oscillation envelopes that are suppressed with respect to what we would expect when no magnetic activity is at play further hampers the detection of a second granulation component (the high-frequency component; Corsaro et al. 2017), which makes the three Harvey-like profile model less likely and adequate at reproducing the given datasets.
A161, page 3 of 9 Fig. 2. Example of two different background models fit using Diamonds + Background.The models show two Harvey-like profiles for TIC 408842743 and three for TIC 38511251.This results in a larger number of free parameters that need to be estimated in the latter case with respect to the former.The PSD of the star is shown in gray, and its smoothing by ∆ν is shown in black.The dot-dashed yellow line represents the level of the white noise, and the Harvey-like components are indicated with dot-dashed blue lines.The dot-dashed magenta line shows the Gaussian envelope of the solar-like oscillations.The overall fit of the background without (with) the Gaussian envelope is marked with a solid red (dashed green) line.
The preliminary step before fitting the model to the data was creating priors for each of the free parameters of the model.This was done by exploiting a public Python routine that is available within the Background code package6 .In particular, the routine in charge of building up the prior distributions operates as follows: (i) It evaluates a median signal in a frequency range close to the Nyquist frequency to obtain a guess on the level of the photometric white noise W. (ii) It uses the input guess of ν max found from the PASTEL catalog (by combining log g and T eff as in Eq. ( 1)) to evaluate the priors on amplitude and characteristic frequency of the Harvey-like components of the granulation activity, as shown by Kallinger et al. (2014) and by Corsaro et al. (2017).(iii) It obtains the maximum level of the signal at very low frequency (below the characteristic frequency of any subsequent Harvey-like profile corresponding to granulation) to provide an estimate of the amplitude for the low-frequency component, which is not directly related to the granulation activity of the star.Additionally, the priors on the Gaussian envelope of the solar-like oscillations were generated by using the input ν max to compute an estimate of the standard deviation of the envelope through the empirical relation (e.g.Mosser et al. 2012), and an estimate of ∆ν obtained from the scaling relation (e.g.Stello et al. 2009;Huber et al. 2011), which was used to evaluate a smoothing of the PSD by a window proportional to ∆ν.This smoothed PSD was then adopted to determine the highest power in the range ν max ± 3∆ν, so that a prior on H osc can be set.All prior distributions were made uniform, which is found to provide better computational speed during the nested-sampling Monte Carlo process than using more complex distributions.

Bayesian detection of solar-like oscillations
After this first set of fits was carried out for every star in the sample, we proceeded with the evaluation of the detection of the oscillation power excess.We followed a procedure along the lines presented by Corsaro & De Ridder (2014) and Müllner et al. (2021), who adopted the same fitting code.In addition to the first set of background models described above, we considered a new set of background models that contained the same components as those that were fit at the beginning, except that no Gaussian envelope was included.This new set of models therefore allowed us to test the hypothesis that oscillations are absent, or are at least not detected, in the observed A161, page 4 of 9 stellar PSD.We used the same input prior distributions as we adopted for the first set of fits, but now removed those related to the parameters of the Gaussian envelope.When the fits for the second set of models were also finalized, we performed the Bayesian model comparison by evaluating the Bayes factor as follows: where E 1 is the Bayesian evidence computed by fitting the model incorporating the Gaussian envelope of the solar-like oscillations (i.e., the detection hypothesis), and E 0 is the Bayesian evidence when the Gaussian envelope is not included (i.e., the nondetection hypothesis).We followed the standard Jeffreys empirical scale of strength for the evidence (Jeffreys 1961), which suggests that only when ln B 1,0 ≥ 5.0 can we claim a clear detection of the oscillation signal, while falling into regimes of moderate to weak detection for 1 ≤ ln B 1,0 < 5 (see also Trotta 2008), meaning that there is still a chance that power excess is present, but more data are needed to confirm this.We present our results in Sect. 4

Measuring ∆ν
This part of the analysis aims to measure the global asteroseismic parameter ∆ν.This analysis can be conducted only after a solar-like oscillation envelope is found.In our case, we applied this additional analysis not only to all stars for which a detection provided a strong evidence condition, but also more in general for stars for which an oscillation envelope was either still likely or could not be completely ruled out.
To estimate ∆ν, we relied on a classic method that enabled us to identify a regular structure in the data, namely the autocorrelation function (ACF; e.g.see Mosser & Appourchaux 2009).The middle panels of Fig. 3 show an example of this application for the stars TIC 375621179 and TIC 408842743.The searching range for the estimation of ∆ν was considered to be ±30% of an initial ∆ν guess based on ν max , which in turn was evaluated using the ∆ν -ν max relation provided by Huber et al. (2011).The region of the PSD we used to evaluate the ACF was centered around ν max .In principle, this region should be set to be as large as possible so that it encompasses a larger number of oscillations modes, but at the same time, it should not be too large to avoid including extended low-S/N regions that can hamper the ACF signal.We therefore compared the extension of the range when considering either ν max ± 2σ env or ν max ± 3∆ν, and for each star, we took the larger of the two.This is because when σ env is at least two times ∆ν, the oscillation envelope spans over several radial orders, and it is therefore more likely that the chances to identify a stronger ACF signal are increased.The ACF 2 was then fit in two steps using a nonlinear least-squares fit each time.In the first step, we adopted a combination of a Gaussian and a quadratic function (from four to six terms), and in the second step, we used the results from the first step to set up the starting points for a second fit using just a simple Gaussian function (three terms only, hence without any quadratic component).We took as the best fit the fit whose Gaussian had a centroid that was closest to the position of the maximum of the ACF 2 .The final estimate of ∆ν was therefore the centroid of the best-fit Gaussian function, with a corresponding uncertainty obtained from the same fit.
To provide an additional level of support to the analysis presented here and in Sect.3.1, we adopted another method for extracting the global asteroseismic properties of ν max and ∆ν.As a check, we applied this different method to the stars belonging to the TIGRE catalog, which is the largest catalog of those adopted for the sample selection in Sect. 2. In this case, we used light curves obtained with the lightkurve code package (Lightkurve Collaboration 2018).In order to determine ν max , we applied a 2D autocorrelation function (ACF2D) method (Huber et al. 2009;Viani et al. 2019) to the stellar power spectrum according to varying window widths (10-100 µHz) and spacing parameters.In this method, we gradually applied a shift to the data to find an optimal correlation.To diagnose the ν max values, we checked the highest S/N oscillation frequency lag and also the correlation metric in the ACF2D.We determined the central oscillation frequency of each segment that was related to the strength of the ACF2D in the oscillation frequency lag.As the ACF2D signal increases, the maximum oscillation frequency finds the best value.Then, we determined the peak of the ACF2D and considered the center of the Gaussian function fit to the peak as the final ν max value.In the ν max region, we further probed the ACF2D to possibly estimate ∆ν.Depending on ν max , the full width at half maximum (FWHM) of the mode envelope was shaped.Using well-established empirical relations in Mosser et al. (2010) and Lund et al. (2017), we collapsed the calculated ACF2D region to determine an approximate ∆ν in a manner similar to the method used in ν max .It is well known that the ACF2D is sensitive to regularly shaped peaks.Therefore, we fit a Gaussian function to the peak between the nearest empirical ∆ν to determine the observed ∆ν.For the ν max and ∆ν uncertainties, we applied the method of Huber et al. (2009).For high-S/N levels, the uncertainties of ν max and ∆ν are above 10 µHz and 0.1 µHz, respectively.A comparison of the measurements obtained from the Bayesian background fitting presented in Sect.3.2 (for ν max ) and the ACF analysis presented at the beginning of this section (for ∆ν), with those obtained from the ACF2D method (for both ν max and ∆ν), yields an agreement of ∼4% for ν max and ∼8% for ∆ν on average.This agrees with the uncertainties obtained in this work.

Measuring A max
The last part of the analysis is related to the estimation of the global oscillation amplitude A max .The amplitudes of solarlike oscillations are known to be related to the fundamental properties of stars, such as mass, luminosity, and temperature (e.g.Kjeldsen & Bedding 1995;Corsaro et al. 2013).However, since oscillation amplitudes originate from the interplay between the excitation and damping mechanism, they are subject to stochastic effects.To properly measure the global oscillation amplitudes, we therefore relied on the recipe developed by Kjeldsen et al. (2005Kjeldsen et al. ( , 2008)).This involves the smoothing of the PSD by a Gaussian FWHM of 4∆ν, the subtraction of the background level, and the rescaling by c/∆ν, with c = 3.04 to obtain the power per radial mode.An example of the result of this analysis is shown in the right panel of Fig. 3.After the smoothed amplitude spectrum was obtained, A max was evaluated as the amplitude at the position of ν max .We note that this value may not necessarily coincide with the absolute maximum of the smoothed amplitude curve.As an uncertainty in A max , we considered the largest variation of A max found within the range ν max ± σ env , which is representative of the breadth of the oscillation envelope region.Our rationale stems from the hypothesis that a flatter oscillation envelope, indicative of a lower S/N, might result in a less constrained A max compared to narrower, more peaked envelopes.The extent and flatness of this envelope assume significance because a pronounced suppression effect is induced by the magnetic fields.Our uncertainty therefore reflects A161, page 5 of 9 Fig. 3. Example of measuring ∆ν and A max for the stars TIC 375621179 (HD 188512, top) and TIC 408842743 (HD 187691, bottom).Left panels: Stellar PSD (in gray).The background fit without the Gaussian envelope is overlaid (in red).The fitted Gaussian envelope above the background is marked with a dashed blue line.The vertical blue band corresponds to a 3σ Bayesian credible region around ν max (vertical dashed orange line) as obtained from Diamonds.Middle panels: ACF 2 function to estimate ∆ν as computed on the original PSD (light gray curve) and on the PSD smoothed by ∆ν/10 (thick black curve).The dot-dashed blue line marks the position of ∆ν from the scaling relation.The selected Gaussian fit is shown as the orange curve, and its centroid is indicated by the vertical dashed orange line.The final estimate of ∆ν is represented by the green bullet, and its 1σ uncertainty is overlaid.Right panels: Smoothed amplitude spectrum (in gray).The estimate of A max and its 1σ uncertainty are marked by the magenta bullet at the position of ν max (vertical dotted blue line).The shaded blue region depicts the range over which the uncertainty on A max is computed, ν max ± σ env .
how much A max is subject to change within the central region of the oscillation envelope.Furthermore, our delineation of the uncertainty is inherently linked to ν max , as expected because A max is by definition evaluated at the position of ν max .At the same time, this uncertainty represents a more conservative estimate than is derived from only mapping the Bayesian credible region estimated for the parameter ν max .We note that this is a different definition of the uncertainty than was given by Huber et al. (2011), who relied on Monte Carlo simulations to inject different realizations of the photometric noise.This latter definition inherently underestimates the uncertainty as it overlooks potential effects such as variability induced by the stellar magnetic cycle, stochastic oscillations, and fluctuations in granulation amplitude.
Nonetheless, it is crucial to note that while the S/N and flatness of the oscillation envelope contribute significantly to the level of statistical significance determined through the model comparison (as detailed in Sect.3.2 and Tables 1 and 2), other factors concurrently influence this determination.They include the frequency range of the oscillation region, the shape of the background signal below the oscillation envelope, potential deviations of the oscillation envelope from an ideal Gaussian bell curve, and the number of data bins within the modeled region.Additionally, the as-yet-uncalibrated effect of suppression induced by magnetic activity in this specific TESS target sample (e.g.Chaplin et al. 2011a) further complicates this assessment.This implies that a more precise estimate of A max is not necessarily associated with a stronger detection of the oscillation envelope.A detailed study of the reliability of the adopted uncertainties in A max is envisioned for a subsequent work, where the amplitudes will be analyzed in the context of the suppression effect induced by the magnetic activity.

Results
For the sake of completeness, in addition to the stars for which a detection of the oscillation envelope is highly likely (i.e., ln B 1,0 ≥ 5), we decided to report all the stars resulting in |ln B 1,0 | < 5, meaning that the detection is either a likely scenario or that at least it cannot be totally ruled out.The value of the Bayes factor for each star in the sample therefore provides a statistical weight on the detection, so that the user can decide whether to rely on the result, depending on the needs.All targets for which a power excess is not favored (i.e., ln B 1,0 < 0) still exhibit a potential indication of an oscillation envelope that future observations may eventually confirm or disprove.We list all the results in Table 1 for the stars for which we find a strong detection condition, meaning that they satisfy the Bayesian model selection process because they show a clear and prominent oscillation power excess, and in Table 2, we list all the other stars for which we were unable to draw a firm conclusion.The stars are listed by increasing HD identifier for each catalog, where the asteroseismic parameters of ν max , ∆ν, and A max are obtained as outlined in Sect.3.Moreover, Tables 1 and 2 appear to show that many of the stars analyzed in this work belong to the catalog source of Olspert et al. (2018), for which a magnetic activity cycle period and rotation period measurements are also available.This makes these targets even more suitable and interesting for future dedicated studies.
A161, page 6 of 9 Because of the low S/N that is typically found in the datasets analyzed in this work, the fitting of the background models was a challenging task.In order to make Diamonds + Background converge to a stable solution, it was mandatory to perform the fits multiple times for many of the targets.This was done to improve the choice of the initial prior distributions in a way that allows the likelihood maximum in the parameter space of the solutions to be fully resolved.This was conducted by means of a hierarchical approach, using the former posterior solution as input to refine the prior distributions of a subsequent inference.
Despite the low S/N of the data, the measurement of ∆ν through the computation of the ACF 2 produced a clear signal of a comb-like pattern, that is, reliable results in most stars, well in line with the prediction obtained from the scaling relation (Huber et al. 2011; within just a few percent).The middle panels of Fig. 3 clearly show that the ACF 2 on the smoothed PSD performs significantly better than that applied on the original PSD.This is because when smoothing the PSD by some fraction of ∆ν (in this case, ∆ν/10), most of the stochasticity that characterizes this type of data can easily be removed, resulting in a net improvement in showing the comb-like structure that is typical of the solar-like oscillations.
Our recipe for measuring A max is not absolute; other definitions may be used.However, the way of measuring the global oscillation amplitudes described in Sect. 3 has commonly been adopted in the literature, even for observations based on radial velocity measurements rather than photometric ones.This clearly allows us to rely on a standard reference definition, which can be used to place our new measurements into context with previous ones, such as those published by Huber et al. (2011) for a large number of Kepler targets, in particular, for the stars observed within the Sounding Stellar Cycles with Kepler program by Karoff et al. (2009; see also Karoff et al. 2013;Bonanno et al. 2014), comprising targets with both solarlike oscillations and measured magnetic activity.
We note that for the star TIC 118572803 (HD 22049), which has the highest ν max of the sample, the analysis of ∆ν and A max was more challenging than for other stars.It was not possible to constrain A max to better than a mere 100% precision.Even ∆ν appears to deviate with respect to the scaling relation value (it is found to be below it by more than 10%).We ascribe these discrepancies to two main aspects, namely the fact that ν max of this star is close to the Nyquist frequency of the dataset, which was set to about 4167 µHz, and therefore does not permit fully resolving the extension of the potential power excess, and the very low S/N of the expected asteroseismic signal, which in turn can be attributed to the combined effect of magnetic activity and small oscillation amplitudes that are typical of a star in this regime of high oscillation frequencies (e.g.Corsaro et al. 2013).As expected, the Bayesian model comparison for this star suggests an inconclusive detection process, so that the reported result has to be considered with caution.

Discussion and conclusion
The analysis conducted in this work enabled us to derive global asteroseismic properties of ν max , ∆ν, and A max for a sample of 34 stars exhibiting measurable levels of magnetic activity, which are available from the catalogs of Mittag et al. (2011), Marsden et al. (2014), Boro Saikia et al. (2018), Olspert et al. (2018), comprising also observations from the Mount Wilson HK project (except for HD 118972).Fifteen of these 34 stars show strong evidence for the detection of a solar-like oscillation envelope, while for 6 further targets, our analysis suggests that the presence of oscillations is a likely scenario, but it cannot be confirmed with the current datasets.For an additional 13 stars, we cannot rule out that the oscillation envelope is not detected (ln B 1,0 ≤ 0), and therefore, we decided to include these targets in our compiled list.
Tables 1 and 2 show that the sample we analyzed covers stars from the main sequence up to the red giant branch, that is, the sample spans a considerable range in stellar evolution.We note that despite the 27 days of observation, the highest S/N stars (corresponding to stars with the highest values of ln B 1,0 ) exhibit a PSD that can be further analyzed in terms of detailed oscillation properties, that is, individual mode oscillation frequencies, line widths, and amplitudes, through a dedicated peak bagging analysis (see e.g.Corsaro et al. 2015Corsaro et al. , 2020)).This guarantees an additional level of accuracy in the obtained modeling results and opens the possibility of investigating the impact of magnetic activity on the individual oscillation mode parameters.
The precision attained on our inferred global asteroseismic properties is adequate for a subsequent analysis to investigate how these properties relate to the measured level of magnetic activity.This will be done in a subsequent work currently in preparation through the inclusion of the stellar ages and of the suppression of the oscillation amplitudes.Following the work done by Bonanno et al. (2014), Bonanno & Corsaro (2022), the results provided here set the basis for a dedicated modeling of the targets with the aim to understand how stellar ages (which will be obtained through an asteroseismic-based inference) are connected to the level of oscillation amplitude suppression, rotation, and metallicity, that is, of magnetic activity in stars.
Future missions such as ESA PLATO (Rauer et al. 2014) will enable a significant extension of this type of studies to stars covering a wider range of fundamental stellar properties, and they will likely increase the final outcome by about two orders of magnitude already through the inclusion of stars belonging to the P1 and P2 (primary) samples of bright targets.Solar-like oscillations constitute a key source of information that enable us to infer the internal structure of stars.This work in conjunction with the analysis of the interplay between magnetic activity cycles and the dynamical properties of rotation and convection (e.g.Bonanno & Corsaro 2022) therefore represents a main path to follow for the years to come in order to shed light on the role of magnetic activity in the stellar evolution.Finally, the stars provided here are all relatively bright stars, and therefore, additional follow-up studies involving spectropolarimetry and ground-based photometry observations might be possible in principle.This will enable the inference of the magnetic field topology and will thus contribute to improve our understanding of the behavior and role of magnetic fields in stars with solar-like oscillations.
Fig.1.Data preparation process for the light curve of TIC 95340781(HD 88737).Each row represents a different iteration (four in total), as outlined in Sect. 2. The first column shows the light curve before the filtering is applied (in orange), the second column shows the (residual) light curve after the filtering is applied (in gray), and the third column shows the PSD of the residual light curve (in blue).The solid green line shows the smoothing of the light curve at each step, and the yellow line in the third column represents the smoothing of the PSD (i.e., an approximation of the background level; see Sect.3.1) by a default window of 30 bins.The red arrow marks the position of a power excess typical of solar-like oscillations, here represented by the ν max parameter and the modeled Gaussian function (dot-dashed red line) obtained from the background fit (see Sect. 3).The light-curve units change from counts during the first iteration to ppm for the remainder of the iterations.

Table 1 .
Global asteroseismic parameters ν max , ∆ν, and A max and their corresponding 1σ uncertainties as obtained in this work.