This is our reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars" by T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gomez. We hope that both the comment and our reply lead to fruitful discussions which of the two presented scenarios is more likely.

The current and upcoming high precision photometric surveys such as TESS, CHEOPS, and PLATO will provide the community with thousands of new exoplanet candidates. As a consequence, the presence of such a correlation is crucial in selecting the targets with the lowest RV jitter for efficient RV follow-up of exoplanetary candidates. Studies of this type are also crucial to design optimized observational strategies to mitigate RV jitter when searching for Earth-mass exoplanets. Our goal is to assess the correlation between high-precision photometric variability measurements and high-precision RV jitter over different time scales. We analyze 171 G, K, and M stars with available TESS high precision photometric time-series and HARPS precise RVs. We derived the stellar parameters for the stars in our sample and measured the RV jitter and photometric variability. We also estimated chromospheric Ca II H $\&$ K activity indicator $log(R' _{HK})$, $\textit{v sin i}$, and the stellar rotational period. Finally, we evaluate how different stellar parameters and a RV sampling subset can have an impact on the potential correlations. We find a varying correlation between the photometric variability and RV jitter as function of time intervals between the TESS photometric observation and HARPS RV. As the time intervals of the observations considered for the analysis increases, the correlation value and significance becomes smaller and weaker, to the point that it becomes negligible. We also find that for stars with a photometric variability above 6.5 ppt the correlation is significantly stronger. We show that such a result can be due to the transition between the spot-dominated and the faculae-dominated regime. We quantified the correlations and updated the relationship between chromospheric Ca II H $\&$ K activity indicator $log(R' _{HK})$ and RV jitter.

Young and active stars generally have regular, almost sinusoidal, patterns of variability attributed to their rotation, while the majority of older and less active stars, including the Sun, have more complex and non-regular light-curves which do not have clear rotational-modulation signals. Consequently, the rotation periods have been successfully determined only for a small fraction of the Sun-like stars observed by transit-based planet-hunting missions, such as CoRoT, Kepler, and TESS. This suggests that only a small fraction of such systems have been properly identified as solar-like analogs. We apply a new method for determining rotation periods of low-activity stars, like the Sun. The method is based on calculating the gradient of the power spectrum (GPS) of stellar brightness variations and identifying a tell-tale inflection point in the spectrum. The rotation frequency is then proportional to the frequency of that inflection point. In this paper test this GPS method against available photometric records of the Sun. We apply GPS, autocorrelation functions, Lomb-Scargle periodograms, and wavelet analyses to the total solar irradiance (TSI) time series obtained from the Total Irradiance Monitor (TIM) on the Solar Radiation and Climate Experiment (SORCE) and the Variability of solar IRradiance and Gravity Oscillations (VIRGO) experiment on the SOlar and Heliospheric Observatory (SoHO) missions. We analyse the performance of all methods at various levels of solar activity. We show that the GPS method returns accurate values of solar rotation independently of the level of solar activity. In particular, it performs well during periods of high solar activity, when TSI variability displays an irregular pattern and other methods fail. Our results suggest that the GPS method can successfully determine the rotational periods of stars with both regular and non-regular light-curves.

Twenty-four years after the discoveries of the first exoplanets, the radial-velocity (RV) method is still one of the most productive techniques to detect and confirm exoplanets. But stellar magnetic activity can induce RV variations large enough to make it difficult to disentangle planet signals from the stellar noise. In this context, HD41248 is an interesting planet-host candidate, with RV observations plagued by activity-induced signals. We report on ESPRESSO observations of HD41248 and analyse them together with previous observations from HARPS with the goal of evaluating the presence of orbiting planets. Using different noise models within a general Bayesian framework designed for planet detection in RV data, we test the significance of the various signals present in the HD41248 dataset. We use Gaussian processes as well as a first-order moving average component to try to correct for activity-induced signals. At the same time, we analyse photometry from the TESS mission, searching for transits and rotational modulation in the light curve. The number of significantly detected Keplerian signals depends on the noise model employed, which can range from 0 with the Gaussian process model to 3 with a white noise model. We find that the Gaussian process alone can explain the RV data while allowing for the stellar rotation period and active region evolution timescale to be constrained. The rotation period estimated from the RVs agrees with the value determined from the TESS light curve. Based on the data that is currently available, we conclude that the RV variations of HD41248 can be explained by stellar activity (using the Gaussian process model) in line with the evidence from activity indicators and the TESS photometry.

Context. Considerable effort has been put into using light curves observed by space telescopes such as CoRoT, Kepler and TESS for determining stellar rotation periods. While rotation periods of active stars can be reliably determined, the light curves of many older and less active stars (e.g. stars similar to the Sun) are quite irregular, which hampers determination of their periods. Aims. We examine the factors causing the irregularities in stellar brightness variations and develop a method for determining rotation periods of low activity stars with irregular light curves. Methods. We extend the Spectral And Total Irradiance Reconstruction (SATIRE) approach for modelling solar brightness variations to Sun-like stars. We calculate the power spectra of stellar brightness variations for various combinations of parameters defining the surface configuration and evolution of stellar magnetic features. Results. The short lifetime of spots in comparison to the stellar rotation period as well as the interplay between spot and facular contributions to brightness variations of stars with near solar activity cause irregularities in their light curves. The power spectra of such stars often lack a peak associated with the rotation period. Nevertheless, the rotation period can still be determined by measuring the period where the concavity of the power spectrum plotted in the log-log scale changes sign, i.e. by identifying the position of the inflection point. Conclusions. The inflection point of the (log-log) power spectrum is found to be a new diagnostic for stellar rotation periods that is shown to work even in cases where the power spectrum shows no peak at the rotation rate.

A former member of the Hyades cluster, $\iota$ Horologii ($\iota$ Hor) is a planet-hosting Sun-like star which displays the shortest coronal activity cycle known to date ($P_{\rm cyc}\sim$ 1.6 yr). With an age of $\sim$$625$ Myr, $\iota$ Hor is also the youngest star with a detected activity cycle. The study of its magnetic properties holds the potential to provide fundamental information to understand the origin of cyclic activity and stellar magnetism in late-type stars. In this series of articles, we present the results of a comprehensive project aimed at studying the evolving magnetic field in this star and how this evolution influences its circumstellar environment. This paper summarizes the first stage of this investigation, with results from a long-term observing campaign of $\iota$ Hor using ground-based high-resolution spectropolarimetry. The analysis includes precise measurements of the magnetic activity and radial velocity of the star, and their multiple time-scales of variability. In combination with values reported in the literature, we show that the long-term chromospheric activity evolution of $\iota$ Hor follows a beating pattern, caused by the superposition of two periodic signals of similar amplitude at $P_{1} \simeq 1.97 \pm 0.02$ yr and $P_{2} \simeq 1.41 \pm 0.01$ yr. Additionally, using the most recent parameters for $\iota$ Hor b in combination with our activity and radial velocity measurements, we find that stellar activity dominates the radial velocity residuals, making the detection of additional planets in this system challenging. Finally, we report here the first measurements of the surface longitudinal magnetic field strength of $\iota$ Hor, which displays varying amplitudes within $\pm4$ G and served to estimate the rotation period of the star ($P_{\rm rot} =~7.70^{+0.18}_{-0.67}$ d).

We report on the determination of the lithium abundance [A(Li)] of 52 solar-like stars. For 41 objects the A(Li) here presented corresponds to the first measurement. We have measured the equivalent widths of the 6708\AA\ lithium feature in high-resolution spectroscopic images ($R \sim 80\,000$), obtained at the Observatorio Astrof\'isico Guillermo Haro (Sonora, Mexico), as part of the first scientific observations of the revitalized Lunar and Planetary Laboratory (LPL) Echelle Spectrograph, now known as the Cananea High-resolution Spectrograph (CanHiS). Lithium abundances were derived with the Fortran code MOOG, using as fundamental input a set of atmospheric parameters recently obtained by our group. With the help of an additional small sample with previous A(Li) determinations, we demonstrate that our lithium abundances are in agreement, to within uncertainties, with other works. Two target objects stand out from the rest of the sample. The star BD+47 3218 ($T_{\rm eff}$ = 6050$\pm$52 K, A(Li) = 1.86$\pm$ 0.07 dex) lies inside the so-called lithium desert in the the A(Li)--$T_{\rm eff}$ plane. The other object, BD+28 4515, has an A(Li) = 3.05$\pm$0.07 dex, which is the highest of our sample and compatible with the expected abundances of relatively young stars.