Kepler main-sequence solar-like stars: surface rotation and magnetic-activity evolution
A. R. G. Santos D. Godoy-Rivera A. J. Finley S. Mathur R. A. García S. N. Breton A. -M. Broomhall

1



] 
 
 





Kepler main-sequence solar-like stars: surface rotation and magnetic-activity evolution
    [
    April 29, 2024
=======================================================================================







§ 


While the mission's primary goal was focused on exoplanet detection and characterization,  made and continues to make extraordinary advances in stellar physics. Stellar rotation and magnetic activity are no exceptions.  allowed for these properties to be determined for tens of thousands of stars from the main sequence up to the red giant branch. From photometry, this can be achieved by investigating the brightness fluctuations due to active regions, which cause surface inhomogeneities, or through asteroseismology as oscillation modes are sensitive to rotation and magnetic fields. This review summarizes the rotation and magnetic activity properties of the single main-sequence solar-like stars within the  field. We contextualize the  sample by comparing it to known transitions in the stellar rotation and magnetic-activity evolution, such as the convergence to the rotation sequence (from the saturated to the unsaturated regime of magnetic activity) and the Vaughan-Preston gap. While reviewing the publicly available data, we also uncover one interesting finding related to the intermediate-rotation gap seen in  and other surveys. We find evidence for this rotation gap in previous ground-based data for the X-ray luminosity. Understanding the complex evolution and interplay between rotation and magnetic activity in solar-like stars is crucial, as it sheds light on fundamental processes governing stellar evolution, including the evolution of our own Sun. 

 

§ KEYWORDS:
 Stars: activity, Stars: evolution, Stars: late-type, Stars: low-mass, Stars: magnetic field, Stars: rotation, Stars: solar-type, starspots 




§ INTRODUCTION


Low-mass stars with convective outer layers, also known as solar-like stars, can sustain internal dynamos and have the potential to harbor magnetic activity. Magnetic fields and magnetic cycles are generated by an interaction between differential rotation and convection <cit.>. As the strong magnetic fields emerge at the stellar photosphere, they form active regions, where dark spots appear, usually in pairs or groups of opposite polarity <cit.>. Such active regions can be associated with eruptive events like flares and coronal mass ejections <cit.>. As active regions decay their magnetic fields disperse and concentrate at the edges of the convective cells, forming the bright faculae <cit.>. In the chromosphere, these regions of intermediate magnetic-field strength appear as bright plage. All these phenomena are part of the star's magnetic activity, which varies at different timescales. In particular, the Sun undergoes an 11-year cycle of magnetic activity. Analogously, other stars are known to exhibit magnetic activity cycles <cit.>, with cycle periods ranging from a few years to over 20 years. As the dynamo mechanism is powered by the interplay of convection and rotation <cit.>, the cycle and rotation periods are found to be related, with slow-rotating stars having longer cycles than fast-rotating stars. Stars are often seen to group along two branches: active and inactive <cit.>. However, there is still debate whether these are properly determined or even exist <cit.>.

The level of magnetic activity is also intrinsically linked to rotation <cit.>. At the beginning of their main-sequence lifetime, stars are relatively fast rotators and exhibit high activity levels <cit.>. Stars gradually lose angular momentum due to their magnetized winds <cit.>, in a process known as magnetic braking. The rate at which the stars spin down depends on their rotation, with faster rotators losing angular momentum faster than slower rotators. Eventually, stars will converge into the so-called slow-rotation sequence (see for example Fig. 3 in ), and from that point onwards the rotation rate decays proportionally to the square root of the age, known as the Skumanich spin-down law <cit.>. The spin-down process is also mass-dependent, with lower-mass stars taking longer to converge into the rotation sequence than higher-mass stars, but once reached, lower-mass stars spin down faster than higher-mass stars <cit.>. The magnetic activity also decays with time: as stars evolve and spin down, they gradually become less active <cit.>. Therefore, generally, fast-rotating stars have stronger magnetic activity than slow-rotating stars. This activity-rotation relationship can be represented as a function of the Rossby number, Ro <cit.>. Ro can be defined as the ratio between the star's rotation period and its convective turnover timescale. The latter corresponds to the typical timescale for convective motions in the stars' envelopes and remains mostly constant during the main-sequence lifetime, increasing as the effective temperature decreases <cit.>. Metallicity, however, can complicate this picture. Due to a larger opacity, metal-rich stars have deeper convection zones than their metal-poor counterparts <cit.>. A deeper convection zone leads to a more vigorous dynamo, consequently to higher magnetic activity. Observational evidence for this effect has been found in both large samples <cit.> and in the particular case of the seismic solar-analog HD 173701 (KIC 8006161). This ∼1M_⊙ star exhibits a magnetic cycle with more than twice the amplitude of the solar cycle <cit.>. Strong magnetic activity, in turn, would lead to a more efficient loss of the angular momentum. Thus, metal-rich stars are expected to spin down faster than metal-poor stars <cit.> and, so far, this theoretical expectation has been supported by observations <cit.>.

The Skumanich spin-down law led to the development of gyrochronology <cit.>, enabling estimation of stellar ages from surface rotation. Analogously, magnetochronology and magnetogyrochronology relations have also been established <cit.>. As rotation and magnetic activity measurements are available for large numbers of stars, these techniques are powerful tools for estimating stellar ages. However, the evolution of stellar rotation and magnetic activity is not yet fully understood, as we will discuss below. 

The stars' magnetic fields can be measured through Zeeman Broadening of spectral lines and Zeeman Doppler Imaging <cit.>. These techniques shed light on the evolution of stellar magnetic fields and allow the large-scale magnetic field topology to be recovered. The ZDI technique has been used to follow entire stellar magnetic cycles which, in the case of 61 Cyg A, share many similarities to that of the solar cycle <cit.>. However, due to the challenges of directly measuring magnetic fields, indirect measures are often used. These are magnetic activity proxies and link to phenomena associated with the presence of strong magnetic fields, which can be constrained by spectroscopic and photometric observations. 

In spectroscopy, magnetic activity is often constrained through the analysis of particular spectral lines. Examples of such lines are the Ca ii H & K lines in the near ultra-violet (NUV) and the Ca ii infrared triplet (IRT). In the presence of active regions, these absorption lines show emission at their cores due to chromospheric heating. By measuring such emission and removing the basal and photospheric contributions, these lines provide proxies for magnetic activity in the chromosphere <cit.>. Since active regions go in and out of view as the star rotates, these activity proxies can show short-term quasi-periodic variations, allowing to constrain rotation periods <cit.>. Another activity proxy that can be derived from spectra is the spot filling factor <cit.>, which follows upon the fact that magnetic spots in solar-like stars are cool features in comparison to their surroundings. 
 
In white-light photometry, active regions lead to variations in the stellar brightness, known as rotation modulation. Once again, as active regions co-rotate, in this case with the stellar surface, the periodicity of these brightness variations is related to the surface rotation period at the latitudes where the active regions are formed <cit.>. The amplitude of the brightness variations relates to the surface coverage by active regions <cit.> and is found to be well correlated with chromospheric activity proxies <cit.>. One of the advantages associated with the advent of photometric space missions is that they allow the retrieval of magnetic activity proxies and rotation for stellar samples orders of magnitude larger than those from ground-based surveys.

Stellar flares can also be detected and characterized from white-light photometry <cit.>, providing a proxy for magnetic activity.
Another widely used magnetic activity proxy is X-ray emission from hot plasma confined in coronal loops <cit.>, which is correlated with the stellar wind mass-loss rates <cit.>. Other activity indicators used in the literature are the Hα emission and NUV excess <cit.>. 

Asteroseismology can also be effective in measuring rotation and magnetic activity in solar-like stars, thanks to the acoustic modes (p modes) being sensitive to both properties. In this case, solar-like acoustic oscillations give information mostly from subphotosphere layers <cit.>. Rotation can be measured from the splitting of the modes' azimuthal orders if the stellar inclination angle is not too small <cit.>. In addition to an average rotation <cit.>, asteroseismology also allows us to obtain information about surface latitudinal differential rotation <cit.>. Unfortunately, main-sequence solar-like stars pose a significant challenge in the determination of their radial differential rotation due to uncertainties in observations and stellar models <cit.>. Moreover, the number of visible modes is limited, as well as their sensitivity to greater depths in the stellar interiors. Magnetic activity affects different properties of the acoustic modes <cit.>. Particularly, the mode frequencies are observed to increase with the magnetic activity level, while mode amplitudes decrease. For low-degree modes, those that are possible to observe for stars other than the Sun, modes of different angular degrees are affected differently by stellar activity <cit.>, depending on the latitudes where active regions emerge (active latitudes). This fact reveals another capability of asteroseismology, in this case, to constrain active latitudes in stars by investigating the magnetic signatures in modes of different angular degrees, as it was done for the well-characterized solar-analog HD 173701 by <cit.>. However, the suppression of mode amplitudes by magnetic activity prevents the detection of acoustic modes in stars with strong magnetic activity <cit.>. Therefore, seismic samples are biased towards weakly active slow rotators. 

The  mission <cit.>, launched by NASA (National Aeronautics and Space Administration), provided one of the most significant contributions to the expansion of stars with known surface rotation and measured activity levels <cit.>. Comparatively, Gaia, launched by ESA (European Space Agency), has yielded a much larger number of rotation measurements <cit.>. Nevertheless,  and Gaia yields are complementary <cit.>, with  probing typically slower rotators than Gaia, including stars similar to our Sun (for reference, at 5000 K, the 5th and 95th percentiles of   distribution are 8.2 and 38.8 days, while the analogous limits for Gaia are 0.4 and 12.6 days).  revealed two potential deviations to the Skumanich spin-down law, whose origins are still under debate.

Most of the  main-sequence sample has already converged to the slow-rotation sequence, where the Skumanich spin-down law is generally assumed to be valid. However,  data suggests the existence of a transition within this regime, with the surface rotation distribution being bimodal <cit.>, which is particularly evident at low temperatures, resulting on an intermediate-rotation gap (∼ 15 days at 4500 K). Spin-down stalling, likely associated with this gap, is evident in stellar clusters with ages around 1 Gyr <cit.>. Recently the intermediate-rotation gap was found in K2 <cit.> and ground-based data for partially convective stars <cit.>, but is absent in fully convective stars <cit.>. Different hypotheses to explain the gap were proposed, with the core-envelope coupling theory gaining traction <cit.>. In this scenario, the angular momentum transfer between the stars' fast core and slow envelope would momentarily stall the spin-down. Once the coupling is completed, the Skumanich-like spin-down would resume. 

The second transition  unveiled concerns relatively old main-sequence stars, around the age of the Sun and older. Given their asteroseismic ages, some of these stars spin faster than expected if their spin-down was consistent with the Skumanich law <cit.>. This observation led to the formulation of the weakened magnetic braking (WMB), which would take place around the middle of the main-sequence lifetime, at a given critical Ro <cit.>. Despite the observational support for WMB, its physical cause(s) remains unclear. As the efficiency of angular momentum transport is primarily governed by the stellar magnetic field <cit.>, one possible explanation is that the stellar dynamo becomes less efficient, leading to a weaker or more complex magnetic field configuration <cit.>. This hypothesis has been investigated using spectropolarimetric observations of solar-like stars, which recover the large-scale magnetic field strength and topology <cit.>. Other explanations range from the influence of latitudinal differential rotation <cit.>, to decreases in the stellar wind mass-loss rates due to closed magnetic fields <cit.> or less efficient wind heating/acceleration <cit.>. As the Sun lies around this transition, attempts have been made to measure the present-day solar wind braking torque <cit.>. Both observations and numerical models show the braking torque to be a factor of two to three times smaller than required by the Skumanich relation <cit.>, lending support to the weakened braking hypothesis.

Looking to the magnetic-activity evolution, another transition may exist, the so-called Vaughan-Preston (VP) gap <cit.>. The VP gap is characterized by a lack of stars with intermediate Ca ii H & K emission. However, there is extensive debate in the community surrounding its existence. In particular, when exploring larger and more complete samples, the VP gap attenuates and in some cases almost disappears <cit.>. Although <cit.> did not find a clear gap, the authors also found evidence supporting a phase of rapid evolution, consistent with the parameter space of the VP gap. Therefore, it is not clear yet whether the VP gap is a result of a transition in the magnetic-activity evolution or a result of observational bias. So far, there is no evidence for it in  data. 

These recent discrepancies between the observations and the expected behavior reinforce the need for a better understanding of rotation and magnetic-activity evolution. This review places the  rotational sample in the context of the known transitions during the main sequence and describes them in more detail in the following sections. 



§ KEPLER MAIN-SEQUENCE SOLAR-LIKE ROTATION SAMPLE


The  mission provided exquisite data for stellar physics. In addition to high-precision photometry,  monitored the same stars in a continuous, long-term manner, spanning up to 4 years of observations. Still,  data are not free of systematics and instrumental artifacts. Therefore, it is important to correct and calibrate them <cit.>, while preserving stellar signals at long timescales, such as the rotation modulation of slow rotators. Once processed,  data constrained surface rotation periods and photometric magnetic activity for several tens of thousands of solar-like stars from the main-sequence to the red-giant phase <cit.>. While the focus of this review is the main-sequence (MS) solar-like stars, it is worth noting that magnetic fields and activity are also found in earlier spectral types <cit.>.

We begin with the sample of 55,252 stars with known rotation periods from <cit.>[KEPSEISMIC light curves were adopted in these works and are available on MAST (Mikulski Archive for Space Telescopes): DOI: 10.17909/t9-mrpw-gc07; https://archive.stsci.edu/prepds/kepseismic/https://archive.stsci.edu/prepds/kepseismic/.], which included subgiant stars by design. To select solely the MS stars we adopt the selection criteria based on the color-magnitude diagram (CMD) from Gaia Data Release 3 <cit.>, as detailed in Appendix A of <cit.>. The magnitudes were corrected for extinction and the selection criteria also remove potential binaries and outliers that sit above or below the MS in the CMD, as well as targets with large (≥1.2) renormalized united weighted error <cit.>, Gaia radial velocity variables <cit.>, stars in the Gaia non-single-star sample <cit.>, and eclipsing binaries <cit.>. As we do not yet fully understand the rotational signals from targets flagged as close-in binary candidates in <cit.>, we keep those that pass the criteria (1,311 targets). This leaves us with a reference  sample of 34,898 single MS stars with known rotation rates. The selection criteria are relatively stringent to ensure a clean sample. Applying the same criteria to the sample of <cit.> would reduce it to 21,685 stars. Comparing the respective clean samples, we verify that the latest catalog still pushed the upper edge of the rotation-period distribution towards slower rotators <cit.>.

The top panel of Figure <ref> shows the  rotation MS sample. The dashed lines mark the upper and lower edges of the  distribution, corresponding to the 95th and 5th percentiles, whose origin is discussed in more detail below. In general, hotter stars are fast rotators than cooler stars, which is expected, as for most of the MS, the magnetic braking is more efficient for less massive stars <cit.>. Another feature that can be seen in the  sample is the so-called intermediate- gap, leading to a bimodal  distribution. For cooler solar-like stars (K and M), a lower-density region, in between two populations or sequences of stars, can be found. For G dwarfs, the lower density region disappears but the  distribution is still bimodal <cit.>. The dot-dashed line indicates the intermediate- gap computed for the clean sample as described in <cit.>. 



For the  sample, the magnetic activity is quantified from the rotational modulation in the light curve using the activity proxy , computed in <cit.> as defined by <cit.>. The bottom panels of Fig. <ref> show the activity-rotation diagram for the  MS sample split by spectral type ( boundaries at 3700, 5200, and 6000 K). Generally, fast rotators are more magnetically active than slow rotators. The black solid lines show the respective density contours, as reference for the figures below. GKM stars are color-coded according to their location concerning the intermediate- gap. The color code for the F stars indicates whether their  is below or above the Kraft break <cit.>. While the activity and rotation of cool F stars are correlated, similarly to the case of the GKM stars, hot F stars tend to be fast rotators with low activity levels <cit.>. The solid blue lines show the  distribution for F stars cooler and hotter than the Kraft break, whose median values are indicated by the dashed and dot-dashed lines (287.8 ppm and 185.8 ppm, respectively). This behavior can be explained by the shallow convective envelopes of the latter, which are unable to produce a strong magnetic field yielding inefficient magnetic braking <cit.>.

Below we discuss the transitions and detection biases in the  rotation sample. We discuss other transitions that happen in a regime not probed by  or, if probed, are not observed in  data. In the activity-rotation figures, we opt to show  and not Ro, as  can be measured directly from the observations. In addition, the comparison between Ro numbers from different studies is not straightforward, as Ro depends on the adopted definition for the convective turnover timescale. Nevertheless, to split the samples into different regimes, when possible, we adopt the Ro number and the respective transitions as determined by the authors in the respective original studies.




§ FROM SATURATED TO UNSATURATED: CONVERGENCE TO THE ROTATION SEQUENCE


The members of young open clusters exhibit a wide range of rotation periods <cit.>, consistent with stars transitioning from an ultra-fast rotation (a result of spin-up due to contraction onto the main sequence) to the converged rotation sequence. These two regimes are separated by a lower-density region, i.e. a gap at ultra-fast rotation <cit.>. The signature of this transition can also be identified in the magnetic activity of low-mass stars, particularly in their X-ray emission <cit.>. Two regimes in the coronal magnetic activity are found by <cit.>: a saturated regime where X-ray emission is almost independent on ; and a unsaturated regime where X-ray emission strongly depends on rotation, corresponding to the converged rotation sequence. The ultra-fast- bimodality found by <cit.> marks the transition from saturated to unsaturated magnetic activity. <cit.> proposed that this transition is related to the core-envelope coupling, where the core and envelope would be coupled for stars in the converged rotation sequence (B2003's interface sequence), and decoupled for ultra-fast rotators (B2003's convective sequence). However, fully convective M dwarfs also follow the same activity-rotation relation characterized by the saturated and unsaturated regimes <cit.>. Alternatively, <cit.> attributed the rotation gap to a transition between weak (saturated) to strong (unsaturated) coupling with the stellar wind, which in turn might be related to a transition from complex to simple magnetic field morphology <cit.>. 

Physically motivated models of angular momentum evolution can describe the observed convergence of rotation, particularly reproducing the mass-dependence of the spin-down <cit.>. Lower-mass solar-like stars spend more time in the saturated regime than higher-mass solar-like stars, which converge earlier into the unsaturated regime. Moreover, the efficiency of magnetic braking depends on the rotation period: stars born with fast rotation spin down faster than stars born with slow rotation. Eventually, they converge into the same sequence, the unsaturated regime, where they lose angular momentum following the Skumanich law and, thus, gyrochronology becomes valid. 

Placing the  sample in context, most of the MS stars have already converged into the unsaturated regime. This is not surprising as the  sample corresponds to field stars, with a mix of populations of typical ages of several Gyr, and thus stars have had time to spin down. In fact, the transition between the saturated and unsaturated regimes happens at faster rotation rates than the lower edge of the   distribution. Both saturated regime and transition are very sparsely populated in the  field. Furthermore, this region of the parameter space, particularly P_rot<7  days, is found to be dominated by tidally-synchronized binaries, as determined by <cit.> and <cit.>. Close-in binary candidates identified by <cit.> that survived the selection criteria (Sect. <ref>) were kept in the  MS sample, but they are found to occupy the same region of the parameter space as the tidally-synchronized binaries (with some common targets). It is, thus, unclear whether the surviving  rapidly-rotating targets are young solar-like stars still in the saturated regime or tidally-synchronized binaries, which were not identified as such yet.

Figure <ref> compares the parameter space of rotation and activity for the X-ray emission sample in <cit.> and for the  MS sample. The left panel show the - diagram, while the right panels show the activity- diagram, where F stars were ignored as they are only a few. The color index is used to compute the convective turnover timescale, according to Eq. 10 in <cit.>, which in turn, together with , is used to compute Ro. Ro can vary across different works, depending on the definition of the convective turnover timescale (e.g. the location where it is measured). In this review, we adopt the Ro values as determined by each work, but we focus on the observable . Later in this section, we will compare the different Ro values. In what follows, the subscript of Ro indicates to which work they refer, particularly we use the initial of the first author and the year of the publication: e.g. W2011 for <cit.>.





In Fig. <ref>, we split the X-ray emission sample into saturated and unsaturated regimes according to <cit.>, with the transition Ro_W2011 set at 0.13 (saturated: Ro_W2011≤0.13; unsaturated: Ro_W2011>0.13). The lower edge of the  distribution of the  sample lies above the transition between saturated and unsaturated regimes in X-ray luminosity. This comparison emphasizes the fact that, with , we have access mostly to the stars in the unsaturated regime. Interestingly, it is noticeable that the relation between X-ray emission and rotation (right-hand side) shows a change in slope around the intermediate- gap, i.e. where the shape of the  contours change. Indeed, one can notice that change in the original figures by <cit.>. This is also shown in Fig. <ref>, where the color code indicates the distance to the intermediate- gap (δlog P_rot). The top right panel displays the X-ray emission against Ro_W2011. Stars in the saturated regime (brown hexagons) are partly omitted to better show the unsaturated regime. Around the intermediate- gap (lighter symbols), there is a decrease in the dispersion compared to the remainder of the unsaturated regime. At  longer than the gap (δlog P_rot>0), stars seem to follow a steeper relation than before. This can be seen through the comparison between the model by <cit.> and the smoothed data (solid red line). In the unsaturated regime, at δlog P_rot<0, the smoothed line closely follows the model, while deviating at δlog P_rot>0. Indeed, splitting the unsaturated regime in two and fitting them separately, we find slopes of -1.92 and -2.46 for δlog P_rot<0 and δlog P_rot>0, respectively. The behavior change can also be seen in the bottom panels, where particularly for the K stars it is noticeable that stars with δlog P_rot<0 and δlog P_rot>0 follow different sequences. This transition in the unsaturated regime was naturally not seen by the authors and is noticed now thanks to the knowledge acquired through the  sample, which revealed this intermediate- gap. It is now clear that the unsaturated regime itself presents multiple regimes. 

The transition from saturated to unsaturated regime corresponding to the ultra-fast- bimodality can also be found in other magnetic activity proxies besides X-ray emission. For the Pleiades, with an age of ∼125 Myr <cit.>, the saturated and unsaturated regimes are identified, for example, in the photometric magnetic activity measured from K2 light curves <cit.> and in the spot filling factor measured from APOGEE (Apache Point Observatory for Galactic Evolution Experiment) spectra <cit.>.

<cit.> discovered a new contribution to the stellar brightness variations, the mid-frequency continuum (MFC). The MFC corresponds to an excess of power between ∼20 and ∼300 μHz in comparison with the models for the acoustic background in the power spectrum, which account for photon-shot noise, activity, and two granulation components. Interestingly, the MFC scales with Ro_B2021, suggesting that it is related to stellar magnetism. Given its timescale, the MFC might be related to the supergranular internetwork (see  for a review). However, the MFC and the photometric magnetic activity (σ_H, measured as the amplitude of the rotation harmonics, H, in the power spectrum, which is well correlated with ) do not vary in phase or follow a similar relation with Ro_B2021. <cit.> found that the MFC also shows two regimes, with the MFC transition taking place at smaller Ro_B2021 (shorter ) than that associated with the ultra-fast- bimodality. This is illustrated in the top panels of Fig. <ref>, where the stars saturated in MFC are shown by the turquoise squares. For the photometric magnetic activity, the transition between the saturated and unsaturated regime (brown hexagons and orange triangles) is consistent with the ultra-fast- bimodality, i.e. with the saturated and unsaturated X-ray regimes. In the top panel of Fig. <ref>, we take the MFC transition from <cit.>, while for the photometric activity, we adopt logRo_B2021,H=-0.7, which is slightly shifted from the value indicated by the authors (logRo_B2021,H=-0.5). This change is motivated by the fact that the stars lying around that Ro_B2021 have already transitioned to the unsaturated regime. , , Ro_B2021, harmonic and MFC amplitudes (σ_H and σ_MFC) are adopted from <cit.>.

In comparison, <cit.> used APOGEE spectra of the Pleiades cluster to estimate the spot filling factor (f_spot), based on the temperature contrast between spots and quiet surroundings. As the activity proxies estimated from K2 light curves and APOGEE spectra are both related to spots, they are expected to show the same behavior. Indeed, that is the case, except with the rotation sequence from <cit.>, which is located at slightly shorter  than that from <cit.>. To split the f_spot sample into saturated and unsaturated regimes for Fig. <ref>, we consider the transition at log Ro_C2022,f_spot=-0.677 as determined by <cit.> with their power-law model. The authors identified potential binary or multiple systems, which are neglected for the graphical representation.



The saturated-unsaturated transition associated with the ultra-fast- bimodality can also be seen in the chromospheric activity proxy measured from the emission in the Ca ii IRT (log R'_IRT). <cit.> investigated the rotation and magnetic activity of NGC 3532, whose age is estimated to be around 300 Myr old <cit.>. The rotation and magnetic-activity data of this cluster show saturated and unsaturated regimes (bottom panels in Fig. <ref>), consistent with the ultra-fast- bimodality. We adopt , , log R'_IRT, and Ro_F2021 from <cit.>. According to the authors all the stars with Ro_F2021<0.06 are in the saturated regime, while stars with Ro_F2021>0.11 have converged to the rotation sequence. Stars in between would be transitioning from one to the other regime.



Figure <ref> compares the activity-Ro diagrams for the stellar samples in Figs. <ref> and <ref>. Each panel shows the Ro values as computed in the different works, where the respective best fits are overlaid. Because of the different definitions, the transition from saturated to unsaturated regime happens at different Ro values. However, as seen above, the transitions for the different activity proxies correspond to the same parameter space in terms of the observable . We prefer to focus on  here, but when comparing Ro from different works, it is recommended to place them in a uniform scale (e.g., by normalizing to the solar value according to the respective definition).

<cit.> tracked down the transition from the saturated to unsaturated regime in the average magnetic field strength, ⟨ B⟩, of M dwarfs (and some K) from CARMENES spectra (Fig. <ref>). Similarly, to the magnetic activity proxies above, the transition identified by the authors at their Ro_R2022=1.3 takes place near the lower edge of the  distribution of  M dwarfs. For the  sample, this edge is not well defined at low  due to small sample sizes. Also, the M-dwarf sample of <cit.> is near the location where the gap closes for fully convective stars <cit.>. Some of the ⟨ B ⟩ measurements are an upper limit. These are mostly located in the more dense region with small ⟨ B⟩, constituting 37% of the stars with  longer than the intermediate- gap. We split the ⟨ B ⟩ unsaturated regime into two according to their location in relation to the intermediate- gap of the  sample. In Fig. <ref> we also show log⟨ B⟩ as a function of Ro_R2022, which was computed using the convective turnover timescale and  provided by <cit.>.



Gaia also provides photometric data that allow the measurement of rotation periods. Focusing on the targets with a relatively large number of visits and long temporal coverage, Gaia reports  and the respective photometric activity proxy (A_Gaia, being the amplitude of the rotation signal) for several hundreds of thousands of stars <cit.>. Indeed, the Gaia DR3 sample, with the most reliable  estimates, includes of 474,026 stars <cit.>. To complement the rigorous  vetting and selection by the authors, which removed evolved stars, we neglect potential binaries (52,933 stars) according to the same criteria used for the  sample (Sect. <ref>). Figure <ref> compares the final Gaia (421,093 stars) and  samples (34,898 single MS stars).  and A_Gaia were taken from <cit.>, while  was taken from <cit.>. The parameter spaces probed by Gaia and the  barely overlap as noted by <cit.> and <cit.>. Since Gaia is not very sensitive to long  values due to its scanning pattern, all its M stars are still in the saturated regime (Fig. <ref>). Gaia G and K stars are spread between both the saturated regime and the “tip” of the unsaturated regime. For those in the unsaturated regime, it is possible to recognize the expected trend with  generally decreasing with increasing . Gaia also unveiled the existence of stars with P_rot<1 day and very low activity (ultra-fast rotator branch), in contrast to the stars in the saturated regime. <cit.> hypothesize that stars can either evolve from the saturated regime directly to the unsaturated regime, or to the ultra-fast rotator branch first and from this to the unsaturated regime. This region of the parameter space cannot be explored through the  sample, as its majority has already converged to the unsaturated regime.



The comparisons above suggest that the lower edge of the  distribution reflects the transition to the unsaturated regime <cit.>.



§ INTERMEDIATE-P_ROT GAP: REGIMES WITHIN THE UNSATURATED REGIME


The intermediate- gap and the associated  bimodality were first discovered in the  field of view by <cit.>. One of the original hypotheses was that the  bimodality resulted from two different star formation episodes and it was a feature of the  field <cit.>. However, since then it has also been recovered in the different campaigns of K2 <cit.>, which focused on different fields of view, and in ground-based data (from the Zwicky Transient Facility survey) covering the full northern hemisphere <cit.>. Therefore, these findings suggest that the  bimodality and intermediate- gap are linked to stellar evolution. 

A second hypothesis was postulated by <cit.> and <cit.>. These works found evidence for the stars below the intermediate- gap being spot-dominated, while those above the gap being facula-dominated. This led to the proposition that this gap would be related to the transition from the spot- to facula-dominated and it would result from observational biases due to the canceling between dark spots and bright faculae <cit.>. However, other observational evidence and successful modeling support a third hypothesis.  

The hypothesis proposes that this gap has its origin in the core-envelope coupling <cit.>, which leads to a stalling in the envelope's spin-down followed by a period of quick evolution once the coupling is completed <cit.>. Before and after this transition, stars' envelopes (and thus their surface rotation periods) would follow the Skumanich spin-down law. Starting with a decoupled core-envelope, the surface of the stars below the gap would brake due to magnetized winds. During the coupling between the fast core and the slow envelope, the surface spin-down would stall for a relatively short timescale. 
Once the coupling is completed, the spin-down would resume. In this light, fully convective stars would not face this transition. Indeed, that was what <cit.> found in ground-based photometric data of field stars. The authors retrieve a rotational gap in the partially convective stars, but not for fully convective M dwarfs.

The rotational sequence of the 1 Gyr NGC 6811 open cluster is consistent with a stalling in the spin-down of K dwarfs <cit.>, which fits the core-envelope coupling hypothesis. The stalling can be seen through the overlap between the rotational sequences of Praesepe <cit.> and NGC 6811 in the regime of K-dwarfs <cit.>, which would not be expected according to the Skumanich spin-down law. The theoretical models by <cit.>, which account for a mass-dependent core-envelope coupling timescale, are able to reproduce the observations. In particular, the models match the stalling of the spin-down for K dwarfs around 1 Gyr, consistent with the observations of NGC 6811. The mass-dependent coupling leads to a kink in the rotation sequences at older ages, i.e. the rotation-period sequence is no longer monotonic at these ages. As the coupling timescale increases with decreasing mass <cit.>, this kink moves towards lower masses with age. This feature is seen in the observational data of older clusters <cit.>.



Fig. <ref> compares the rotational sequences of stellar clusters of different ages with the  sample. The rotation period and effective temperature for the clusters' members are gathered from different studies <cit.>. Praesepe members are plotted for reference. Most of Praesepe stars with T_eff≳3500 K have already converged to the rotation sequence, while cooler stars have not. The rotation sequence of Praesepe and younger clusters (Fig. <ref>) shows a general monotonic behavior, with  decreasing with . The remainder of the clusters show evidence of a stalling in the spin-down associated with the intermediate- gap, as discussed above: NGC 6811 partially overlaps with Praesepe; and the older clusters show a kink in the rotation sequence at relatively low temperatures. For the 2.7 Gyr Ruprecht 147,  increases with  between ∼ 4000 and ∼ 4800 K, while elsewhere it decreases with . For the 4 Gyr M67,  the rotation sequence starts by decreasing with  until ∼ 3900 K, where the behavior inverts. The  values at ∼ 3900 K are shorter than those at ∼ 4600 K, suggesting that  increases within this interval, despite the lack of data.

Finally, we remind that the imprint of the transition associated with the intermediate- gap can be seen in the relation between X-ray emission and rotation period (Figs. <ref> and <ref>). Similarly, it can also be identified in the relation between chromospheric emission and rotation (Fig. <ref>, discussed below).



§ VAUGHAN-PRESTON GAP: DISCONTINUITY IN THE MAGNETIC ACTIVITY 


Focusing on stars with longer periods than the intermediate- gap, another gap has been observed, possibly related to a transition at later stages of stellar evolution. Contrary to the other sections, this transition is not seen in the rotation period. Instead, this transition is seen in the chromospheric activity, where there is a scarceness of stars with intermediary Ca ii H & K emission <cit.>. This is known as the Vaughan-Preston (VP) gap. The origin of the VP gap has been contested since its finding: is it of astrophysical origin or an observational bias? Indeed, when investigating larger stellar samples, the gap is not observed anymore <cit.>. Nevertheless, <cit.> concluded that their results are consistent with chromospheric-emission bimodality and associated transitions.

Figure <ref> compares the Mount Wilson Observatory <cit.> sample with the  sample: , Ro_E2017, log  R'_HK, and (B-V) were adopted from <cit.>. For the MWO sample (monitored over two decades), we transform the color index (B-V) into , taking into account the metallicity when available. We retrieve metallicity estimates from <cit.>, <cit.>, and APOGEE <cit.> for 60 stars, with a mean value of -0.07 dex. For the remainder, we opt to assume solar metallicity. For the (B-V)- conversion, we use the publicly available routine in [https://github.com/sczesla/PyAstronomyhttps://github.com/sczesla/PyAstronomy] <cit.> based on <cit.>[https://pyastronomy.readthedocs.io/en/latest/pyaslDoc/aslDoc/aslExt_1Doc/ramirez2005.htmlhttps://pyastronomy.readthedocs.io/en/latest/pyaslDoc/aslDoc/aslExt_1Doc/ramirez2005.html]. One of the targets is ignored because it is outside of the valid parameter space (B-V>1.51). The stars below and above the intermediate- gap (Sect. <ref>) are shown in orange and beige, respectively: the dash-dotted line obtained for  was used to split the stars. The VP gap located at log R'_HK=-4.75 is shown by the blue dotted line. As shown by <cit.>, there is no evidence of the VP gap in the  data, which is several times larger than the ground-based chromospheric emission samples. On the one hand, this could mean that indeed there is an observational bias in the ground-based surveys, and when increasing the sample size, the gap is no longer found. On the other hand, although  observations are relatively long-term, their 4-yr length is still limiting especially when dealing with a variable property such as magnetic activity. Taking solar data,  <cit.> showed how the inferences from 4-yr observations on magnetic activity and its variation can change over time, depending on the phase of the cycle. The  variation strongly depends on the average activity level of the star and on its rotation rate. Consequently, the limiting 4-year monitoring of stars can lead to a smearing of the data points, potentially hiding the possible VP gap. In this case, the large sample size can contribute to a greater concealment.





§ MIDLIFE TRANSITION: WEAKENING OF THE MAGNETIC BRAKING


In recent years, evidence for a possible weakening of the magnetic braking was found <cit.>. The surface rotation of the  MS asteroseismic targets was found to be too fast in comparison to what was expected given their asteroseismic age <cit.>. To explain the observations, <cit.> proposed that the magnetic braking would become less efficient at a critical Ro, the so-called weakening of the magnetic braking (WMB). In this case,  evolves according to the Skumanich relation until reaching the critical Ro where its evolution is no longer driven by magnetic braking.

The origin of this transition is presently uncertain. One theory suggests that stellar differential rotation weakens around the critical Ro <cit.>, disrupting the stellar dynamo process <cit.>. In numerical simulations, this can be associated with a shift from solar-like differential rotation, which has a fast equator and slow poles, to solid body or even anti-solar rotation <cit.>. It has been suggested that this transition could result in a change in magnetic topology from large-scale to small-scale fields <cit.>. However, as the large-scale magnetic field typically governs the efficiency of the magnetic braking <cit.>, and spectropolarimetric observations have shown it does not abruptly disappear at the critical Ro <cit.>, it appears more likely that the overall stellar magnetic field strength weakens <cit.>. This decrease may also reduce the stellar wind mass-loss rate <cit.>, further weakening the wind braking. In other theories, the evolution of the latitudinal differential rotation is sufficient to stall the rotation-evolution of stars around the critical Ro, due to the active latitudes that couple surface rotation to the stellar wind <cit.>.


The magnetic cycles of Sun-like stars are also observed to become longer <cit.>, with some appearing to lose all cyclic variability akin to the Maunder minimum observed for the Sun <cit.>. Interestingly, the Sun is near this critical Ro, which raised the question of whether the Sun could be in transition <cit.> and whether the Maunder Minimum could be a symptom of this. Observations point towards the Sun's magnetic braking being two to three times smaller than required by the Skumanich relation <cit.>. However, the magnetic braking timescale for the Sun is around 10-100 Myrs, so this disagreement could be explained by long-term variation in solar activity that exceeds the available ∼10,000 years of cosmogenic radionuclide records <cit.>.

<cit.> and <cit.> identified an over-density of stars close to the upper edge of the   distribution. This might support the WMB as old stars would remain with a largely unchanged surface rotation, which in turn would result in a pileup of stars with different activity levels near this edge. This upper edge is also recovered in <cit.> who simulated the  field population using angular momentum evolution models accounting for the WMB. However, we must note that the  surface rotation sample does not probe very well this region of the parameter space as it coincides approximately with the current detection limit, which depends on the magnitude <cit.>. In both cases, if the stars follow the Skumanich spin-down law or if they face a WMB, stars beyond this point in the main sequence are weakly active and have small amplitude brightness variations due to active regions. Thus, the   upper edge can be either the result of one of these effects or the combination of both.

Asteroseismology could come to the rescue as it is easier to detect acoustic oscillations in weakly active stars <cit.>. This is also illustrated in Fig. <ref>, as the seismic sample tends to be close to the upper edge of  distribution with small  values. Sadly, the  main-sequence asteroseismic sample with a high enough signal-to-noise ratio to measure rotational splittings is relatively small (94 stars of different mass and metallicity). Nevertheless, <cit.> was able to retrieve seismic rotation periods, P_seismic, for this sample confirming that old stars are rotating faster than expected. Therefore, this discrepancy cannot be solely explained by an observational bias in the  sample. 

Figure <ref> compares the  seismic-rotation and surface-rotation samples. The 94 stars from <cit.> are shown in dark blue, where their rotation period corresponds to that measured through asteroseismology. Cross-matching this sample with <cit.>, we find 56 stars (light blue), out of the 94, with measured surface rotation and . In the right, only the 56 stars are shown (twice, in light and dark blue), as  is not available for the remainder. We opt to show both  and P_seismic to demonstrate that both sets of measurements occupy the same parameter space. However, we note that there are some discrepancies between the seismic surface rotation periods <cit.>. Some differences are expected, arising from differential rotation, as the observed acoustic modes are most sensitive to the subsurface layers <cit.> and might be sensitive to different latitudes in comparison to the active-region latitudes. Additionally, there are uncertainties inherent to both techniques. While seismic detections are easier for weakly active stars, which tend to be slow rotators, these are associated with small effects on the acoustic modes, hampering the estimation of P_seismic. Furthermore, mode linewidths increase with effective temperature <cit.>, which poses a challenge because, depending on their values, they can lead to an overlap between azimuthal components, preventing constraints on P_seismic. There is also a correlation between rotational splitting and inclination angle <cit.>. As our sensitivity to the active latitudes depends on the inclination angle, inclination has an impact on the amplitude of the light curve modulation due to active regions. Nevertheless, the impact on surface rotation estimate should not be significant. Surface rotation measurements can also be biased toward the half of the true  <cit.>, but this corresponds to a small percentage of the targets and efforts were made to avoid such misestimations <cit.>.
Finally, photometric pollution by nearby stars can also influence the  estimate.



While, as described above, there is evidence supporting the WMB, we emphasize that this regime is at the  detection limit of the surface rotation and that the seismic sample is relatively small. Factors such as the spectral type (many are F stars, with thin convective layers), differential rotation, active latitudes, and chemical composition <cit.> can impact rotation and magnetic-activity evolution. Based on ground-based spectroscopic data of solar twins, <cit.> concluded that their results favor a gradual rotation evolution, rather than a broken spin-down law. Future observations are needed to shed light on the evolution of these properties beyond the middle of the star's lifetime.



§ CONCLUSIONS AND PERSPECTIVES


Undoubtedly pivotal for exoplanetary and stellar physics, with the end of its second life –the K2 mission–,  concluded its operations in the fall of 2018. Despite this, to this day,  data continue to be the focal point of numerous studies, consistently contributing to significant discoveries in the field. Considering stellar rotation and magnetic activity, CoRoT <cit.> had already provided an important sneak-peak into the potential of space-based photometry <cit.>. But  allowed us to investigate for the first time the rotation and magnetic activity of tens of thousands of solar-like stars. In particular, despite its observational limitations,  extended the number of old weakly active stars, like the Sun, with measured rotation and magnetic activity.

In this review, concerning , we focused mostly on the information from the rotation modulation and briefly on asteroseismic inferences. However, it is important to acknowledge that flares can also provide information on the magnetic activity of stars and the environment around them <cit.>. We place the  sample in the context of established or potential transitions in the evolution of rotation and magnetic activity of MS solar-like stars.  revealed two unexpected transitions: the intermediate- gap and the WMB. The former is thought to be a transition in the stellar evolution related to the core-envelope coupling, while the physical mechanisms driving the latter are still under debate. One of the main reasons for this debate arises from observational limitations, because of the small number of stars in this regime with measured rotation periods, as these stars are weakly active and below the photometric threshold detection.

Stars are born with strong magnetic activity and rapid rotation rates. As they spin down due to magnetized winds, they eventually converge onto a narrow rotation sequence, transitioning from a saturated regime to an unsaturated regime in magnetic activity.  however did not probe the saturated regime well, as most of its stars are relatively old. The transition between saturated and unsaturated regimes coincides roughly with the lower edge of the   distribution (P_rot∼ 8 days at T_eff=5000 K).

The unsaturated regime itself splits into different regimes separated by the intermediate- gap (P_rot∼ 12 days at T_eff=5000 K). Currently, the most plausible hypothesis for this observed transition is the core-envelope coupling. This coupling would momentarily stall the surface spin-down due to the angular momentum transfer between the fast-rotating core and the slow-rotating envelope. This transition can also be noticed in the X-ray emission sample in <cit.>, which was not established before.

In ground-based data, evidence of a possible gap or transition at intermediate magnetic activity has also been found, corresponding to the Vaughan-Preston (VP) gap (P_rot∼ 30 days at T_eff=5000 K). Although there is no evidence of the VP gap in  data, it cannot be discarded. On one hand,  could support the case for observational bias and the nonexistence of such a gap. On the other hand, the 4-year timespan of  observations may be insufficient to fully cover magnetic cycles, especially of solar analogs, leading to a cumulative dispersion that may be enough to conceal the VP gap. This shows the need for longer-term observations to reach a more complete knowledge of the magnetic activity.

The  rotation distribution is characterized by a well-defined upper edge (P_rot∼ 40 days at T_eff=5000 K), with the Sun lying near this edge. In recent years, signs of a change in stellar evolution around this edge have been found in seismic observations of solar-like stars, leading to the postulation of the weakening of magnetic braking (WMB). However, the WMB is still under debate, due to observational limitations arising from the hard-to-detect small-amplitude signals. In addition, instrumental artifacts hamper the recovery of rotation periods for slow rotators <cit.>. If these artifacts can be better mitigated, it will be possible to extend the number of detections in this regime.

This review also shows that the observed behavior of the different activity proxies is similar. In spite of being sensitive to different layers of the stellar atmosphere and different magnetic features, the activity proxies still relate to each other. As a consequence, the well-established transitions in the stellar evolution are found in all. This also reinforces the validity of photometry to measure and investigate magnetic activity.

Gaia has prompted an extraordinary expansion of the number of stars with known rotation and magnetic activity proxies, allowing these properties to be measured for several hundreds of thousands of stars. However, the parameter spaces of  and Gaia samples do not overlap significantly. Indeed they are complementary: with Gaia it is possible to investigate young ultra-fast rotators, which were mostly absent from  data. Gaia rotation sample yielded the discovery of ultra-fast rotators with unexpected weak magnetic activity. Furthermore, as can be seen above, it is also possible that some of the Gaia stars are crossing the intermediate- gap. Thus, in spite of the small sample size in this region of the parameter space, Gaia might also provide insights into this transition.

Currently, NASA's TESS <cit.> is mostly sensitive to the fast strongly active rotators. Nevertheless, future TESS extended missions can expand the reach of this mission, but the photometric precision still will not allow the detection of small amplitude signals.

Future ESA's PLATO <cit.> mission, planned to be launched in 2026, will observe stars with high precision and continuously for at least two years. PLATO is expected to greatly increase the number of MS solar-like stars and subgiants with seismic detections. The relatively long-term observations will also allow us to measure rotation and magnetic activity both from rotational modulation due to active regions and from asteroseismology. However, as discussed above, one can argue that observations longer than two years are needed to provide a better characterization of magnetic activity. Still, thanks to its high precision, PLATO will be crucial in expanding the observational limit towards more weakly active slower rotators, moving the upper edge of the observed  distribution towards longer , which is currently a major challenge in the field. This will provide more information on the rotation and magnetic-activity evolution beyond the age of the Sun, particularly in the regime of the proposed WMB. PLATO is also expected to increase the overlap between seismic and surface rotation samples, both by increasing seismic detections with high SNR and by pushing the detection limit of surface rotation. This overlap will provide complementary information about the stars, particularly accurate ages and different measures of magnetic activity and rotation. This will greatly improve our understanding of stellar evolution and also differential rotation. Moreover, PLATO will observe bright stars, which is a clear advantage in relation to , meaning that complementary ground-based spectroscopic observations will also be possible. This will allow having independent constraints on rotation and magnetic activity, as well as better characterization of the atmospheric parameters of the stars, namely effective temperature and metallicity. Particularly, better metallicity measurements will help to improve our knowledge of the impact of chemical composition on the magnetic-activity and rotation evolution.

Having observations in different wavelengths will also bring complementary information on the magnetic activity in the different layers of the stellar atmosphere. Furthermore, while dark spots dominate in the passband of the photometric missions mentioned above at the rotation timescale <cit.>, in other wavelengths the bright facula or plage become the most prevalent. As these bright features live longer than spots, they have the potential to produce more stable signals. Therefore, moving away from the optical can increase the rotation yields for stars like the Sun and older <cit.>. Again, this will be particularly important to investigate the WMB and understand the evolution beyond the solar age. Still, long-term observations are required to properly characterize magnetic activity.

PLATO will also be important to investigate the existence of the VP gap, however, the smearing effect described above can be worse in the case of the 2-year observations (in comparison to the  4-year observations). PLATO will also provide more observations for stars around the intermediate- gap and perhaps provide better age constraints for these than what we had before. This can help us to better depict and understand the transition associated with this gap.

One of the most relevant applications of rotation and magnetic activity is age-dating stars. Asteroseismology provides the most precise way to estimate stellar ages, but it is available for a small number of stars. For example, magnetic activity is known to suppress the already small amplitudes of acoustic modes, making seismic detections not possible for relatively high activity levels. In contrast, the detection and characterization of rotation modulation thanks to active regions is easier for such active stars. In principle, this would allow us to provide stellar ages for a large number of stars through gyrochronology, magnetochronology, and/or gyromagnetochronology. However, one needs to better understand the evolution of rotation and magnetic activity. Therefore, unless we understand (or at the very least calibrate) all of these transitions discussed above, the power of age-dating using rotation and magnetic-activity proxies is heavily endangered <cit.>. Furthermore, a better understanding of stellar magnetism is also important for exoplanetary physics. Magnetic activity and rotation can disguise or even mimic planetary signals, hampering the detection and characterization of planets <cit.>. Also, these stellar properties have a significant impact on the habitability of the planets, leading for example to the loss of their atmospheres and shaping the orbital architecture of the systems <cit.>. This also reinforces the need for reliable stellar ages, including those determined from rotation and/or magnetic activity, as they will shed light on the evolution of the planetary systems and on the evolution of exoplanet atmospheres.



§ AUTHOR CONTRIBUTIONS


ARGS: Conceptualization, Formal analysis, Writing – original draft; DGR: Formal analysis, Writing – review & editing; AJF: Writing – original draft; SM: Conceptualization, Writing – review & editing; RAG: Conceptualization, Writing – review & editing; SNB: Writing – review & editing; A-MB: Conceptualization, Writing – review & editing.



§ FUNDING

This work was supported by Fundação para a Ciência e a Tecnologia (FCT) through national funds and by Fundo Europeu de Desenvolvimento Regional (FEDER) through COMPETE2020 - Programa Operacional Competitividade e Internacionalização by these grants: UIDB/04434/2020, UIDP/04434/2020, & 2022.03993.PTDC. 
ARGS acknowledges the support from the FCT through the work contract No. 2020.02480.CEECIND/CP1631/CT0001. DGR and SM acknowledge support from the Spanish Ministry of Science and Innovation (MICINN) grant no. PID2019-107187GB-I00. AJF acknowledges support from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 810218 WHOLESUN). SM acknowledges support from the Spanish Ministry of Science and Innovation (MICINN) with the Ramón y Cajal fellowship no. RYC-2015-17697, PID2019-107061GB-C66 for PLATO, and through AEI under the Severo Ochoa Centres of Excellence Programme 2020–2023 (CEX2019-000920-S). RAG acknowledges the support from PLATO and GOLF CNES grants. SNB acknowledges support from PLATO ASI-INAF agreement n. 2015-019-R.1-2018. A-MB acknowledges the support from STFC consolidated grant ST/T000252/1. 






§ DATA AVAILABILITY STATEMENT



The data used in this review is publically available in the respective original studies as reported.





Frontiers-Harvard 


