Solar magnetic activity shows both smooth secular changes, such as the Grand Modern Maximum, and quite abrupt drops that are denoted as Grand Minima. Direct numerical simulations (DNS) of convection driven dynamos offer one way of examining the mecha
nisms behind these events. In this work, we analyze a solution of a solar-like DNS that has been evolved for roughly 80 magnetic cycles of 4.9 years, during which epochs of irregular behavior are detected. The emphasis of our analysis is to find physical causes for such behavior. The DNS employed is a semi-global (wedge) magnetoconvection model. For data analysis we use Ensemble Empirical Mode Decomposition (EEMD) and phase dispersion ($D^2$) methods. A special property of the DNS is the existence of multiple dynamo modes at different depths and latitudes. The dominant mode is solar-like. This mode is accompanied by a higher frequency mode near the surface and a low-frequency mode in the bottom of the convection zone. The overall behavior of the dynamo solution is very 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. We interpret the overall irregular behavior to be 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 study LQ Hya photometry for 1982-2014 with the carrier fit (CF) -method and compare our results to earlier photometric analysis and recent Doppler imaging maps. We first utilize different types of statistical methods to estimate various candidates
for the carrier period for the CF method. Secondly, a global fit to the whole data set and local fits to shorter segments are computed with the period that is found to be the optimal one. The harmonic least-squares analysis of all the available data reveals a short period close to 1.6 days as a limiting value for a set of significant frequencies. We interpret this as the rotation period of the spots near the equatorial region. In addition, the distribution of the significant periods is found to be bimodal, hinting of a longer-term modulating period, which we set out to study with a two-harmonic CF model. Weak modulation signal is, indeed retrieved, with a period of roughly 6.9 years. The phase dispersion analysis gives a clear symmetric minimum for coherence times lower than and around 100 days. We interpret this as the mean rotation period of the spots (1.60514 days), and this value is chosen to be used as the carrier period for the CF analysis. With the CF method we seek for any systematic trends in the spot distribution in the global time frame, and locally look for abrupt phase changes earlier reported in rapidly rotating objects. During 2005-2008 the global CF reveals a coherent structure rotating with a period of 1.6037 days, while during most other times the spot distribution appears rather random in phase. The evolution of the spot distribution of the object is found to be very chaotic, with no clear signs of an azimuthal dynamo wave that would persist over longer time scales, although the short-lived coherent structures observed occasionally do not rotate with the same speed as the mean spot distribution.
