## Stellar oscillations show differential rotation

The Sun rotates faster at its equator than at its poles. This process is known as differential rotation and is seen in the motion of sunspots. Helioseismology has shown that the effect extends into the Sun's interior. It has not been possible to measure whether other stars also experience equivalent differential rotation. Benomar *et al.* used the Kepler spacecraft to monitor stellar oscillations of a group of Sun-like stars. By decomposing the oscillations into separate frequencies, they searched for signs of differential rotation. Several stars do indeed seem to have equators that spin faster than their poles, and none indicated the opposite pattern.

*Science*, this issue p. 1231

## Abstract

The differentially rotating outer layers of stars are thought to play a role in driving their magnetic activity, but the underlying mechanisms that generate and sustain differential rotation are poorly understood. We report the measurement using asteroseismology of latitudinal differential rotation in the convection zones of 40 Sun-like stars. For the most significant detections, the stars’ equators rotate approximately twice as fast as their midlatitudes. The latitudinal shear inferred from asteroseismology is much larger than predictions from numerical simulations.

Analysis of acoustic oscillations visible on the Sun’s surface by using helioseismology has been critical for constraining its rotation profile. Helioseismology has revealed that the rotation rate of the Sun’s convection zone decreases with latitude (*1*, *2*). This latitudinal differential rotation has a magnitude of 11% of the average rate from equator to midlatitudes (≈45° latitude) and 30% between equator and the poles. At the base of the convection zone, the Sun transitions to solid-body rotation. How such a rotation profile is established and maintained is still poorly understood. However, it is likely that differential rotation plays a role in sustaining the solar magnetic field through a dynamo mechanism (*3*–*5*).

So far, little is known about the latitudinal differential rotation in other stars, and classical methods for investigating it are primarily sensitive to the near-surface layers. Most studies rely on photometric variability from starspots at different latitudes (*6*) or use Doppler imaging to track magnetic features at the surface and their migration in latitude (*7*, *8*). Another approach involves studying the Fourier transform of the spectroscopic line profiles (*9*, *10*).

Asteroseismology provides an opportunity to probe rotation inside stars, including Sun-like pulsators (*11*–*13*), because it studies the resonant frequencies of waves within the body of the star. Among these are acoustic waves that travel at various depths and have a varying sensitivity to rotation with latitude. These acoustic waves can be used to infer the star’s internal rotation both in radius and latitude (*14*–*17*).

The NASA *Kepler* spacecraft has provided high-precision, long-duration photometric time series for many stars, which is necessary for the study of the differential rotation of Sun-like stars with asteroseismology.

Oscillations of Sun-like stars are driven by the stochastic convective motion of material in the outer envelope of the stars. Each mode of oscillation is identified by the overtone number *n* and spherical harmonic functions of angular degree *l* and azimuthal order *m*. Modes with the same *n* and *l* but different *m* appear as multiplets of 2*l* + 1 components. The splitting of modes in a multiplet provides information about the rotation profile and on the physical processes acting within the star (such as flows, internal structure, and tidal forces).

An efficient means for quantifying the impact of rotation on the split frequencies of acoustic waves is through the use of Clebsch-Gordan *a*-coefficients *a*_{1}, *a*_{2}, *a*_{3}, ... (*18*). This decomposition has been used extensively in helioseismology but not in the analysis of other stars. For stars other than the Sun, only low-degree modes (*l* ≤ 3) can be observed, and so only the coefficients *a*_{1}, *a*_{2}, and *a*_{3} can be determined. The coefficient *a*_{1} is an average of the rotation rate. The *a*_{2} coefficient is related to the asphericity (*19*), and *a*_{3} is a measure of the latitudinal differential rotation (*16*). The coefficient *a*_{3} is positive when the pole rotates slower than the equator—that is, if the star exhibits a solarlike rotation profile (fig. S1). Conversely, *a*_{3} is negative for antisolar rotation profiles. Simulations tend to show that fast rotation causes solarlike rotation, whereas slower rotation rates lead to antisolar rotation (*20*–*22*), although simulations found an antisolar rotation at solar rotation rates (*21*, *22*). This is evidently incompatible with solar observations, and constraints from other stars are necessary to resolve this issue.

Pulsations appear as Lorentzian profiles in a power spectrum of the star’s photometric time series (*23*). The analysis of the pulsations is performed by fitting a model to the power spectrum by using a Markov chain Monte Carlo (MCMC) sampling technique (*11*). The contribution of *a*_{1} to mode frequencies is easily measured, but the absolute value of *a*_{3} is approximately two orders of magnitude smaller than *a*_{1}, so it is only possible to measure it by using multiyear observations from space instruments such as *Kepler*.

We performed a combined measurement of the *a*_{1} and *a*_{3} coefficients for 40 stars of mass between 0.9 and 1.5 *M*_{☉} (where *M*_{☉} is the mass of the Sun), observed as part of the *Kepler* LEGACY sample (*24*, *25*). By integrating the probability density functions of *a*_{3} determined by fitting the power spectra, we computed the detection significance for either solar or antisolar latitudinal differential rotation (Fig. 1 and table S1). Applying a detection threshold at a probability of 84% (excluding *a*_{3} = 0 with a detection significance >1σ for a Gaussian), we found that none of the stars unambiguously show antisolar rotation, whereas 32% (13 stars) show significant solarlike rotation. Five stars have a detection probability of more than 97.5% (or a significance >2σ). The excess of solarlike rotators may be due to an observational limitation. Antisolar differential rotation is theoretically expected (*20*, *26*) for slow rotators, in which *a*_{1} and *a*_{3} are difficult to measure. Our sample consists predominantly of fast rotators, so we may only be sensitive to stars with solarlike rotation profiles. Our most statistically significant detections have relatively high rotation rates compared with those below the detection threshold.

The Rossby number is an indicator of the influence of rotation on a fluid compared with convective motion. Current theory (*20*) predicts that solarlike rotators should have a Rossby number less than unity. For the majority of our significant detections, we found a Rossby number less than 0.8 (fig. S2).

We next sought to estimate the rotational shear between the equator and higher latitudes. For this, we selected two representative stars (HD 173701 and HD 187160) in terms of *a*_{1}, *a*_{3}, and internal structure (thickness of convection zone). Their power spectra are shown in figs. S3 and S4 with the best-fitting model. HD 173701 is a cool dwarf of spectral type G8 (the Sun is of spectral type G2), with a mass 0.974 ± 0.029 *M*_{☉} and age 4.69 ± 0.44 billion years (*24*). It is known to be very active (*27*), with irregular variability on time scales that range from 8 to 40 days, likely caused by starspots. The probability density function deduced from the fit indicates that *a*_{1} = 574.6 ± 82.0 nHz and *a*_{3} = 28.4 ± 12.5 nHz (fig. S5). HD 187160 is of spectral type F9 and has a mass and age of 1.09 ± 0.03 *M*_{☉} and 3.28 ± 0.16 billion years. Thus, HD 187160 is similar to a younger version of the Sun. It has a visible magnitude of 7.4 and a faint, cool spectroscopic binary companion of magnitude 8.7 (*28*). The average internal rotation period is ~9 days (*a*_{1} = 1185.0 ± 51.7 nHz). The light curve of HD 187160 shows variability at two distinct time scales, one at ≈9.4 days and another at ≈17.5 days. We attribute this variability to starspots on the surfaces of HD 187160 and its fainter noninteracting companion, respectively. We measured *a*_{3} = 53.7 ± 25.0 nHz (fig. S6).

We have evaluated the potential systematic errors associated with the mode-fitting methodology using simulated power spectra for HD 173701 and HD 187160 (*23*). These systematic errors are small, meaning that *a*_{3} in both cases is positive, with a probability higher than 97.5%. This confirms that these stars have poles that rotate slower than the equator.

To obtain a latitudinal rotation profile, we used an inversion of the *a*-coefficients (*16*). Because we only measured *a*_{1} and *a*_{3}, we were restricted to a two-parameter rotation model for this. We assumed a solarlike rotation profile Ω(θ) = Ω_{0} + 3Ω_{1}(5cos^{2}θ – 1)/2 in the convection zone, where θ is the colatitude (complementary angle of the latitude so that θ = 0 is the pole and θ = 90° is the equator). The average internal rotation rate is given by Ω_{0}, and the term in Ω_{1} ∝ Ω_{pole} – Ω_{eq} measures the contrast in rotation (latitudinal shear) from the equator to the pole in the convection zone. In the equation, Ω_{pole} and Ω_{eq} are the pole and equator rotation rates, respectively. The convection zone extends from the surface down to 0.687 for HD 173701 and to 0.785 for HD 187160, where is the stellar radius. Because of its relatively fast rotation, HD 187160 was investigated for possible variations of *a*_{1} with the overtone *n* (fig. S7). We did not find evidence of significant variation. Studies of radial differential rotation of low-mass main-sequence stars indicate that variation in *a*_{1} should not exceed 30% (*29*), and we therefore do not expect this to have an effect on the inversion.

The fundamental properties of the stellar models assumed for the inversion and the results of the inversion are provided in tables S2 and S3, respectively. The two-dimensional rotation profiles and the probability density of the rotation rate are shown in Fig. 2 at different latitudes for both stars. By construction, the profile is symmetric around the equator and around the rotation axis, so only one quadrant is shown. The rotation rate in the interior matches that of the envelope at a latitude of 26.5° (θ = 63.5°) because of the structure of the two-zone rotation profile. As is shown by the probability density of the rotation profile, the uncertainty on the rotation rate increases substantially toward the poles, beyond a value of ~45° latitude. At high latitudes, the modes become much less sensitive to rotation, thus yielding less information. This limit is imposed by the lack of visible modes with an angular degree higher than *l* = 2 in the spectrum. We found that both HD 173701 and HD 187160 exhibit a latitudinal shear from equator to 45° latitude ΔΩ_{45}/Ω_{eq}, which is approximately five times greater than that of the Sun.

As suggested by numerical simulations, the fast rotation of HD 187160 means that the functional form of its rotation profile may not be solarlike (*23*). We have considered this possibility for HD 187160 (fig. S8) and found that alternative rotation profiles have larger latitudinal shear than the solarlike rotation profile. However, computation of the Bayesian evidence (necessary to determine the goodness of the fit) shows that the solarlike profile is the most likely model.

Because of loss of angular momentum, such as from magnetic braking, a star’s rotation is expected to slow with its age (*1*). This age-rotation relation is apparent in our data: Fig. 3, A and B, shows *a*_{1} as a function of the age given by stellar models (*25*).

As shown in Fig. 3C, the latitudinal shear between the equator and midlatitude for the subset of 13 stars with significant detections is ≈60%, albeit with a large scatter. HD 173701 and HD 187160 are representative of this ensemble, with a differential rotation of ≈50%. The Sun has a significantly lower shear factor than that of the considered ensemble. The difference is of more than 1σ from the dispersion of the ensemble.

This unexpectedly large shear poses a challenge to theoretical models. The balance between angular momentum transport (owing to anisotropy in the turbulent flow) and small-scale flows, which act as enhanced turbulent viscosity, plays a dominant role in regulating latitudinal shear (*30*). The large shear we found indicates a correspondingly large anisotropy in the turbulence, leading to efficient angular momentum transport and suppression of turbulent viscosity. In a typical stellar convection zone, turbulent anisotropy, driven by rotation, substantially affects large-scale flows. Thus, enhanced angular momentum transport and diminished turbulent viscosity amplify large-scale flows and suppress small-scale flows. This could be caused by a very efficient small-scale dynamo, in which small-scale flows are suppressed (*31*). In addition, large-scale magnetic fields tend to reduce shear through the Lorentz force (*32*). However, our results indicate that Lorentz-force feedback is ineffective in the stars we investigated. Thus, the large-scale magnetic field is likely transported efficiently into the deeper regions of the star in which rigid rotation is expected. For this, magnetic pumping is a candidate mechanism (*33*).

## Supplementary Materials

www.sciencemag.org/content/361/6408/1231/suppl/DC1

Materials and Methods

Figs. S1 to S8

Tables S1 to S4

This is an article distributed under the terms of the Science Journals Default License.

## References and Notes

**Acknowledgments:**The authors thank T. R. White and S. Kamiaka for their helpful discussions.

**Funding:**This work is supported by NYUAD institute grant G1502. Research funding from the German Aerospace Center (Grant 50OO1501) and the Max Planck Society (PLATO Science) is acknowledged. Funding for the Stellar Astrophysics Centre is provided by The Danish National Research Foundation (grant DNRF106). The research was supported by the ASTERISK project (ASTERoseismic Investigations with SONG and

*Kepler*) funded by the European Research Council (grant agreement 267864). Funding for the

*Kepler*mission is provided by the NASA Science Mission directorate.

**Author contributions:**O.B., M.B., and L.G. designed the research. O.B. performed seismic analyses and lead in writing the manuscript. M.B. performed seismic inversion. M.B.N. analyzed the light curves variability. All authors discussed the results and contributed to sections of the manuscript.

**Competing interests:**None declared.

**Data and materials availability:**This paper includes data collected by the

*Kepler*mission, provided by the Kepler Asteroseismic Science Operations Center (KASOC) and available at https://doi.org/10.7910/DVN/8SK6OL. Further details on the data used are listed in table S4. The TAMCMC analysis code is available at https://github.com/OthmanB/TAMCMC-C.