Aims. Although the time of the Maunder minimum (1645-1715) is widely known as a period of extremely low solar activity, it is still being debated whether solar activity during that period might have been moderate or even higher than the current solar cycle #24. We have revisited all existing evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods. We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of the sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum.Results. The level of solar activity during the Maunder minimum is reassessed on the basis of all available datasets. Conclusions. We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely lower than during the Dalton minimum of around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected with a high confidence level.

Context. Late-type stars rotate differentially owing to anisotropic turbulence in their outer convection zones. The rotation is called solar-like (SL) when the equator rotates fastest and anti-solar (AS) otherwise. Hydrodynamic simulations show a transition from SL to AS rotation as the influence of rotation on convection is reduced, but the opposite transition occurs at a different point in the parameter space. The system is bistable, i.e., SL and AS rotation profiles can both be stable. Aims. We study the effect of a dynamo-generated magnetic field on the large-scale flows, particularly on the possibility of bistable behaviour of differential rotation. Methods. We solve the hydromagnetic equations numerically in a rotating spherical shell that typically covers ±75• latitude (wedge geometry) for a set of different radiative conductivities controlling the relative importance of convection. We analyse the resulting differential rotation, meridional circulation, and magnetic field and compare the corresponding modifications of the Reynolds and Maxwell stresses. Results. In agreement with earlier findings, our models display SL rotation profiles when the rotational influence on convection is strong and a transition to AS when the rotational influence decreases. We find that dynamo-generated magnetic fields help to produce SL differential rotation compared to the hydrodynamic simulations. We do not observe any bistable states of differential rotation. In the AS cases we find coherent single-cell meridional circulation, whereas in SL cases we find multi-cellular patterns. In both cases, we obtain poleward circulation near the surface with a magnitude close to that observed in the Sun. In the slowly rotating cases, we find activity cycles, but no clear polarity reversals, whereas in the more rapidly rotating cases irregular variations are obtained. Moreover, both differential rotation and meridional circulation have significant temporal variations that are similar in strength to those of the Sun. Conclusions. Purely hydrodynamic simulations of differential rotation and meridional circulation are shown to be of limited relevance as magnetic fields, self-consistently generated by dynamo action, significantly affect the flows.

Context. Both dynamo theory and observations of stellar large-scale magnetic fields suggest a change from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation. Aims. We seek to understand this transition using numerical simulations. Methods. We use three-dimensional simulations of turbulent magnetohydrodynamic convection in spherical shell wedges and considered rotation rates between 1 and 31 times the solar value. Results. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition coincides with a change in the rotation profile from antisolar- to solar-like differential rotation with a faster equator and slow poles. In the solar-like rotation regime, the field configuration consists of an axisymmetric oscillatory field accompanied by an m = 1 azimuthal mode (two active longitudes), which also shows temporal variability. At slow (rapid) rotation, the axisymmetric (nonaxisymmetric) mode dominates. The axisymmetric mode produces latitudinal dynamo waves with polarity reversals, while the nonaxisymmetric mode often exhibits a slow drift in the rotating reference frame and the strength of the active longitudes changes cyclically over time between the different hemispheres. In the majority of cases we find retrograde waves, while prograde waves are more often found from observations. Most of the obtained dynamo solutions exhibit cyclic variability either caused by latitudinal or azimuthal dynamo waves. In an activity-period diagram, the cycle lengths normalized by the rotation period form two different populations as a function of rotation rate or magnetic activity level. The slowly rotating axisymmetric population lies close to what in observations is called the inactive branch, where the stars are believed to have solar-like differential rotation, while the rapidly rotating models are close to the superactive branch with a declining cycle to rotation frequency ratio and an increasing rotation rate. Conclusions. We can successfully reproduce the transition from axi- to nonaxisymmetric dynamo solutions for high rotation rates, but high-resolution simulations are required to limit the effect of rotational quenching of convection at rotation rates above 20 times the solar value.

Context. Solar-like differential rotation is characterized by a rapidly rotating equator and slower poles. However, theoretical models and numerical simulations can also result in a slower equator and faster poles when the overall rotation is slow. Aims. We study the critical rotational influence under which differential rotation flips from solar-like (fast equator, slow poles) to an anti-solar one (slow equator, fast poles). We also estimate the non-diffusive (Λ effect) and diffusive (turbulent viscosity) contributions to the Reynolds stress. Methods. We present the results of three-dimensional numerical simulations of mildly turbulent convection in spherical wedge geometry. Here we apply a fully compressible setup which would suffer from a prohibitive time step constraint if the real solar luminosity was used. To avoid this problem while still representing the same rotational influence on the flow as in the Sun, we increase the luminosity by a factor of roughly 10 6 and the rotation rate by a factor of 10 2 . We regulate the convective velocities by varying the amount of heat transported by thermal conduction, turbulent diffusion, and resolved convection. Results. Increasing the efficiency of resolved convection leads to a reduction of the rotational influence on the flow and a sharp transition from solar-like to anti-solar differential rotation for Coriolis numbers around 1.3. We confirm the recent finding of a large-scale flow bistability: contrasted with running the models from an initial condition with unprescribed differential rotation, the initialization of the model with certain kind of rotation profile sustains the solution over a wider parameter range. The anti-solar profiles are found to be more stable against perturbations in the level of convective turbulent velocity than the solar-type solutions. Conclusions. Our results may have implications for real stars that start their lives as rapid rotators implying solar-like rotation in the early main-sequence evolution. As they slow down, they might be able to retain solar-like rotation for lower Coriolis numbers, and thus longer in time, before switching to anti-solar rotation. This could partially explain the puzzling findings of anti-solar rotation profiles for models in the solar parameter regime.

Context. For moderate and slow rotation, magnetic activity of solar-like stars is observed to strongly depend on rotation, while for rapid rotation, only a very weak or no dependency is detected. These observations do not yet have a solid explanation in terms of dynamo theory. Aims. To work towards such an explanation, we numerically investigated the rotational dependency of dynamo drivers in solar-like stars, that is, stars that have a convective envelope of similar thickness as in the Sun. Methods. We ran semi-global convection simulations of stars with rotation rates from 0 to 30 times the solar value, corresponding to Coriolis numbers, Co, of 0 to 110. We measured the turbulent transport coefficients describing the magnetic field evolution with the help of the test-field method, and compared with the dynamo effect arising from the differential rotation, self-consistently generated in the models.Results. The trace of the α tensor increases for moderate rotation rates with Co 0.5 and levels off for rapid rotation. This behavior is in agreement with the kinetic α based on the kinetic helicity, if one takes into account the decrease of the convective scale with increasing rotation. The α tensor becomes highly anisotropic for Co 1, α rr dominates for moderate rotation (1

Context. Solar magnetic activity shows both smooth secular changes, such as the modern Grand Maximum, and quite abrupt drops that are denoted as grand minima, such as the Maunder Minimum. Direct numerical simulations (DNS) of convection-driven dynamos offer one way of examining the mechanisms behind these events. Aims.In this work, we analyze a solution of a solar-like DNS that was evolved for roughly 80 magnetic cycles of 4.9 years and where epochs of irregular behavior are detected. The emphasis of our analysis is to find physical causes for such behavior. Methods. The DNS employed is a semi-global (wedge-shaped) magnetoconvection model. For the data analysis tasks we use Ensemble Empirical Mode Decomposition and phase dispersion methods, as they are well suited for analyzing cyclic (non-periodic) signals. Results. A special property of the DNS is the existence of multiple dynamo modes at different depths and latitudes. The dominant mode is solar-like (equatorward migration at low latitudes and poleward at high latitudes). This mode is accompanied by a higher frequency mode near the surface and at low latitudes, showing poleward migration, and a low-frequency mode at the bottom of the convection zone. The low-frequency mode is almost purely antisymmetric with respect to the equator, while the dominant mode has strongly fluctuating mixed parity. The overall behavior of the dynamo solution is extremely complex, exhibiting variable cycle lengths, epochs of disturbed and even ceased surface activity, and strong short-term hemispherical asymmetries. Surprisingly, the most prominent suppressed surface activity epoch is actually a global magnetic energy maximum; during this epoch the bottom toroidal magnetic field obtains a maximum, demonstrating that the interpretation of grand minima-type events is non-trivial. The hemispherical asymmetries are seen only in the magnetic field, while the velocity field exhibits considerably weaker asymmetry. Conclusions. We interpret the overall irregular behavior as being due to the interplay of the different dynamo modes showing different equatorial symmetries, especially the smoother part of the irregular variations being related to the variations of the mode strengths, evolving with different and variable cycle lengths. The abrupt low-activity epoch in the dominant dynamo mode near the surface is related to a strong maximum of the bottom toroidal field strength, which causes abrupt disturbances especially in the differential rotation profile via the suppression of the Reynolds stresses.

We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kramers-like opacity prescription as a function of temperature and density or a static profile of a similar shape. We show that the subadiabatic layer arises due to nonlocal energy transport by buoyantly driven downflows in the upper parts of the CZ. Analysis of the force balance of the upflows and downflows confirms that convection is driven by cooling at the surface. We find that the commonly used prescription for the convective enthalpy flux being proportional to the negative entropy gradient does not hold in the stably stratified layers where the flux is positive. We demonstrate the existence of a non-gradient contribution to the enthalpy flux, which is estimated to be important throughout the convective layer. A quantitative analysis of downflows indicates a transition from a tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a heightindependent number of strong downdrafts persist. This change of flow topology occurs when a substantial subadiabatic layer is present in the lower part of the CZ.

We present results from four convectively-driven stellar dynamo simulations in spherical wedge geometry. All of these simulations produce cyclic and migrating mean magnetic fields. Through detailed comparisons, we show that the migration direction can be explained by an αΩ dynamo wave following the Parker-Yoshimura rule. We conclude that the equatorward migration in this and previous work is due to a positive (negative) α effect in the northern (southern) hemisphere and a negative radial gradient of Ω outside the inner tangent cylinder of these models. This idea is supported by a strong correlation between negative radial shear and toroidal field strength in the region of equatorward propagation.