No Arabic abstract
The Sun rotates differentially with a fast equator and slow pole. Convection in the solar interior is thought to maintain the differential rotation. However, although many numerical simulations have been conducted to reproduce the solar differential rotation, previous high-resolution calculations with solar parameters fall into the anti-solar (fast pole) differential rotation regime. Consequently, we still do not know the true reason why the Sun has a fast-rotating equator. While the construction of the fast equator requires a strong rotational influence on the convection, the previous calculations have not been able to achieve the situation without any manipulations. The problem is called convective conundrum. The convection and the differential rotation in numerical simulations were different from the observations. Here, we show that a high-resolution calculation succeeds in reproducing the solar-like differential rotation. Our calculations indicate that the strong magnetic field generated by a small-scale dynamo has a significant impact on thermal convection. The successful reproduction of the differential rotation, convection, and magnetic field achieved in our calculation is an essential step to understanding the cause of the most basic nature of solar activity, specifically, the 11-year cycle of sunspot activity.
Solar-cycle related variation of differential rotation is investigated through analyzing the rotation rates of magnetic fields, distributed along latitudes and varying with time at the time interval of August 1976 to April 2008. More pronounced differentiation of rotation rates is found to appear at the ascending part of a Schwabe cycle than at the descending part on an average. The coefficient $B$ in the standard form of differential rotation, which represents the latitudinal gradient of rotation, may be divided into three parts within a Schwabe cycle. Part one spans from the start to the $4^{th}$ year of a Schwabe cycle, within which the absolute $B$ is approximately a constant or slightly fluctuates. Part two spans from the $4^{th}$ to the $7^{th}$ year, within which the absolute $B$ decreases. Part three spans from the $7^{th}$ year to the end, within which the absolute $B$ increases. Strong magnetic fields repress differentiation of rotation rates, so that rotation rates show less pronounced differentiation, but weak magnetic fields seem to just reflect differentiation of rotation rates. The solar-cycle related variation of solar differential rotation is inferred to the result of both the latitudinal migration of the surface torsional pattern and the repression of strong magnetic activity to differentiation of rotation rates.
The latitudinal distributions of the yearly mean rotation rates measured respectively by Suzuki in 1998 and 2012 and Pulkkinen $&$ Tuominen in 1998 are utilized to investigate internal-cycle variation of solar differential rotation. The rotation rate at the solar Equator seems to decrease since cycle 10 onwards. The coefficient $B$ of solar differential rotation, which represents the latitudinal gradient of rotation, is found smaller in the several years after the minimum of a solar cycle than in the several years after the maximum time of the cycle, and it peaks several years after the maximum time of the solar cycle. The internal-cycle variation of the solar rotation rates looks similar in profile to that of the coefficient $B$. A new explanation is proposed to address such a solar-cycle related variation of the solar rotation rates. Weak magnetic fields may more effectively reflect differentiation at low latitudes with high rotation rates than at high latitudes with low rotation rates, and strong magnetic fields may more effectively repress differentiation at relatively low latitudes than at high latitudes. The internal-cycle variation is inferred to the result of both the latitudinal migration of the surface torsional pattern and the repression of strong magnetic activity to differentiation.
To explore the physics of large-scale flows in solar-like stars, we perform 3D anelastic simulations of rotating convection for global models with stratification resembling the solar interior. The numerical method is based on an implicit large-eddy simulation approach designed to capture effects from non-resolved small scales. We obtain two regimes of differential rotation, with equatorial zonal flows accelerated either in the direction of rotation (solar-like) or in the opposite direction (anti-solar). While the models with the solar-like differential rotation tend to produce multiple cells of meridional circulation, the models with anti-solar differential rotation result in only one or two meridional cells. Our simulations indicate that the rotation and large-scale flow patterns critically depend on the ratio between buoyancy and Coriolis forces. By including a subadiabatic layer at the bottom of the domain, corresponding to the stratification of a radiative zone, we reproduce a layer of strong radial shear similar to the solar tachocline. Similarly, enhanced superadiabaticity at the top results in a near-surface shear layer located mainly at lower latitudes. The models reveal a latitudinal entropy gradient localized at the base of the convection zone and in the stable region, which however does not propagate across the convection zone. In consequence, baroclinicity effects remain small and the rotation iso-contours align in cylinders along the rotation axis. Our results confirm the alignment of large convective cells along the rotation axis in the deep convection zone, and suggest that such banana-cell pattern can be hidden beneath the supergranulation layer.
We study differential rotation in late-stage shell convection in a 3D hydrodynamic simulation of a rapidly rotating $16M_odot$ helium star with a particular focus on the convective oxygen shell. We find that the oxygen shell develops a quasi-stationary pattern of differential rotation that is neither described by uniform angular velocity as assumed in current stellar evolution models of supernova progenitors, nor by uniform specific angular momentum. Instead, the oxygen shell develops a positive angular velocity gradient with faster rotation at the equator than at the pole by tens of percent. We show that the angular momentum transport inside the convection zone is not adequately captured by a diffusive mixing-length flux proportional to the angular velocity or angular momentum gradient. Zonal flow averages reveal stable large-scale meridional flow and an entropy deficit near the equator that mirrors the patterns in the angular velocity. The structure of the flow is reminiscent of simulations of stellar surface convection zones and the differential rotation of the Sun, suggesting that similar effects are involved; future simulations will need to address in more detail how the interplay of buoyancy, inertial forces, and turbulent stresses shapes differential rotation during late-stage convection in massive stars. Our findings may have implications for neutron star birth spin periods and supernova explosion scenarios that involve rapid core rotation. If convective regions develop positive angular velocity gradients, angular momentum could be shuffled out of the core region more efficiently, potentially making the formation of millisecond magnetars more difficult.
Measuring surface differential rotation (DR) on different types of stars is important when characterizing the underlying stellar dynamo. It has been suggested that anti-solar DR laws can occur when strong meridional flows exist. We aim to investigate the differential surface rotation on the primary star of the RS CVn binary HU Vir by tracking its starspot distribution as a function of time. We also aim to recompute and update the values for several system parameters of the triple system HU Vir (close and wide orbits). Time-series high-resolution spectroscopy for four continuous months was obtained with the 1.2-m robotic STELLA telescope. Nine consecutive Doppler images were reconstructed from these data, using our line-profile inversion code iMap. An image cross-correlation method was applied to derive the surface differential-rotation law for HU Vir. New orbital elements for the close and the wide orbits were computed using our new STELLA radial velocities (RVs) combined with the RV data available in the literature. Photometric observations were performed with the Amadeus Automatic Photoelectric Telescope (APT), providing contemporaneous Johnson-Cousins $V$ and $I$ data for approximately 20 years. This data was used to determine the stellar rotation period and the active longitudes. We confirm anti-solar DR with a surface shear parameter $alpha$ of -0.029 $pm$ 0.005 and -0.026 $pm$ 0.009, using single-term and double-term differential rotation laws, respectively. The best fit is achieved assuming a solar-like double-term law with a lap time of $approx$ 400 d. Our orbital solutions result in a period of 10.387678 $pm$ 0.000003 days for the close orbit and 2726 $pm$ 7 d ($approx$ 7.5 yr) for the wide orbit. A Lomb-Scarge (L-S) periodogram of the pre-whitened $V$-band data reveals a strong single peak providing a rotation period of 10.391 $pm$ 0.008 d.