The tilt of bipolar magnetic region (BMR) is crucial in the Babcock--Leighton process for the generation of the poloidal magnetic field in Sun. Based on the thin flux tube model of the BMR formation, the tilt is believed to be caused by the Coriolis force acting on the rising flux tube of the strong toroidal magnetic field from the base of the convection zone (BCZ). We analyze the magnetic field dependence of BMR tilts using the magnetograms of Michelson Doppler Imager (MDI) (1996-2011) and Helioseismic and Magnetic Imager (HMI) (2010-2018). We observe that the distribution of the maximum magnetic field ($B_{\rm max}$) of BMRs is bimodal. Its first peak at the low field corresponds to BMRs which do not have sunspots as counterparts in the white light images, whereas the second peak corresponds to sunspots as recorded in both types of images. We find that the slope of Joy's law ($\gamma_0$) initially increases slowly with the increase of $B_{\rm max}$. However, when $B_{\rm max} \gtrsim 2$ kG, $\gamma_0$ decreases. Scatter of BMR tilt around Joy's law systematically decreases with the increase of $B_{\rm max}$. The decrease of observed $\gamma_0$ with $B_{\rm max}$ provides a hint to a nonlinear tilt quenching in the Babcock--Leighton process. We finally discuss how our results may be used to make a connection with the thin flux tube model.

Simulations of magnetohydrodynamics convection in slowly rotating stars predict anti-solar differential rotation (DR) in which the equator rotates slower than poles. This anti-solar DR in the usual $\alpha \Omega$ dynamo model does not produce polarity reversal. Thus, the features of large-scale magnetic fields in slowly rotating stars are expected to be different than stars having solar-like DR. In this study, we perform mean-field kinematic dynamo modelling of different stars at different rotation periods. We consider anti-solar DR for the stars having rotation period larger than 30~days and solar-like DR otherwise. We show that with particular $\alpha$ profiles, the dynamo model produces magnetic cycles with polarity reversals even with the anti-solar DR provided, the DR is quenched when the toroidal field grows considerably high and there is a sufficiently strong $\alpha$ for the generation of toroidal field. Due to the anti-solar DR, the model produces an abrupt increase of magnetic field exactly when the DR profile is changed from solar-like to anti-solar. This enhancement of magnetic field is in good agreement with the stellar observational data as well as some global convection simulations. In the solar-like DR branch, with the decreasing rotation period, we find the magnetic field strength increases while the cycle period shortens. Both of these trends are in general agreement with observations. Our study provides additional support for the possible existence of anti-solar DR in slowly rotating stars and the presence of unusually enhanced magnetic fields and possibly cycles which are prone to production of superflare.

One of the most robust features of the solar magnetic cycle is that the stronger cycles rise faster than the weaker ones. This is popularly known as the Waldmeier Effect, which is known for more than 80 years. This fundamental feature of the solar cycle has not only practical implications, e,g., in predicting the solar cycle, but also implications in understanding the solar dynamo. Here we ask the question whether the Waldmeier Effect exists in other Sun-like stars. To answer this question, we analyze the Ca \II{} H \& K S-index from Mount Wilson Observatory for 21 Sun-like G--K stars. We specifically check two aspects of Waldmeier Effect, namely, WE1: the anti-correlation between the rise times and the peaks and WE2: the positive correlation between rise rates and amplitudes. We show that except HD~16160, HD~81809, HD~155886 and HD~161239, all stars considered in the analysis show WE2. While WE1 is found to be present only in some of the stars studied. Further, the WE1 correlation is weaker than the WE2. Both WE1 and WE2 exist in the solar S-index as well. Similar to the solar cycles, the magnetic cycles of many stars are asymmetric about their maxima. The existence of the Waldmeier Effect and asymmetric cycles in Sun-like stars suggests that the dynamo mechanism which operates in the Sun is also operating in other stars.

Observations of chromospheric and coronal emissions from various solar-type stars show that the stellar magnetic activity varies with the rotation rates of the stars. The faster the star rotates, its magnetic activity gets stronger but activity cycle period does not show a straightforward variation with the rotation rate. For slowly rotating stars, the cycle period decreases with the increase of rotation rate, while for the fast rotators dependency of cycle period on rotation is presently quite complicated. We aim to provide an explanation of these observational trends of stellar magnetic activity using a dynamo model. We construct a theoretical dynamo model for stars of mass 1 $M_\odot$ based on the kinematic flux transport dynamo model including radial pumping near the surface of the stars. The inclusion of this near surface downward radial pumping is found to be necessary to match the observed surface magnetic field in case of the Sun. The main ingredients of our dynamo model, meridional circulation and differential rotation for stars are obtained from a mean-field hydrodynamic model. Our model shows a decrease of cycle period with increasing rotation rate in the slowly rotating regime and a slight increase of cycle period with rotation rate for the rapid rotators. The strength of the magnetic field is found to be increasing as the rotation rate of the star increases. We also find that the parity of the stellar magnetic field changes with rotation. According to our model, the parity flips to quadrupolar from dipolar if the rotation period of the star is less than 17 days.

One peculiar feature of the solar cycle which is yet to be understood properly is the frequent occurrence of double peaks (also known as the Gnevyshev peaks). Not only the double peaks but also multiple peaks and spikes are often observed in any phase of the cycle. We propose that these peaks and spikes are generated due to fluctuations in the Babcock-Leighton process (the poloidal field generation from tilted bipolar magnetic regions). When the polar field develops, large negative fluctuations in the Babcock-Leighton process can reduce the net polar field abruptly. As these fluctuations in the polar field are propagated to the new toroidal field, these can promote double peaks in the next solar cycle. When fluctuations in the polar field occur outside the solar maximum, we observe their effects as spikes or dips in the following sunspot cycle. Using an axisymmetric Babcock-Leighton dynamo model we first demonstrate this idea. Later, we perform a long simulation by including random scatter in the poloidal field generation process and successfully reproduce the double-peaked solar cycles. These results are robust under reasonable changes in the model parameters as long as the diffusivity is not too larger than $10^{12}$ cm$^2$s$^{-1}$. Finally, we analyze the observed polar field data to show a close connection between the short-term fluctuations in the polar field and the double peaks/spikes in the next cycle. Thereby, this supports our theoretical idea that the fluctuations in the Babcock-Leighton process can be responsible for the double peaks/spikes in the observed solar cycle.

In 1844 Schwabe discovered that the number of sunspots increased and decreased over a period of about 11 years, that variation became known as the sunspot cycle. Almost eighty years later, Hale described the nature of the Sun's magnetic field, identifying that it takes about 22 years for the Sun's magnetic polarity to cycle. It was also identified that the latitudinal distribution of sunspots resembles the wings of a butterfly showing migration of sunspots in each hemisphere that abruptly start at mid-latitudes towards the Sun's equator over the next 11 years. These sunspot patterns were shown to be asymmetric across the equator. In intervening years, it was deduced that the Sun (and sun-like stars) possess magnetic activity cycles that are assumed to be the physical manifestation of a dynamo process that results from complex circulatory transport processes in the star's interior. Understanding the Sun's magnetism, its origin and its variation, has become a fundamental scientific objective \-- the distribution of magnetism, and its interaction with convective processes, drives various plasma processes in the outer atmosphere. In the past few decades, a range of diagnostic techniques have been employed to systematically study finer scale magnetized objects, and associated phenomena. The patterns discerned became known as the ``Extended Solar Cycle'' (ESC). The patterns of the ESC appeared to extend the wings of the activity butterfly back in time, nearly a decade before the formation of the sunspot pattern, and to much higher solar latitudes. In this short review, we describe their observational patterns of the ESC and discuss possible connections to the solar dynamo as we depart on a multi-national collaboration to investigate the origins of solar magnetism through a blend of archived and contemporary data analysis with the goal of improving solar dynamo understanding and modeling.

The Sun occasionally goes through Maunder-like extended grand minima when its magnetic activity drops considerably from the normal activity level for several decades. Many possible theories have been proposed to explain the origin of these minima. However, how the Sun managed to recover from such inactive phases every time is even more enigmatic. The Babcock--Leighton type dynamos, which are successful in explaining many features of the solar cycle remarkably well, are not expected to operate during grand minima due to the lack of a sufficient number of sunspots. In this Letter, we explore the question of how the Sun could recover from grand minima through the Babcock--Leighton dynamo. In our three-dimensional dynamo model, grand minima are produced spontaneously as a result of random variations in the tilt angle of emerging active regions. We find that the Babcock-Leighton process can still operate during grand minima with only a minimal number of sunspots and that the model can emerge from such phases without the need for an additional generation mechanism for the poloidal field. The essential ingredient in our model is a downward magnetic pumping which inhibits the diffusion of the magnetic flux across the solar surface.

We explore the cause of the solar cycle variabilities using a novel 3D Babcock--Leighton dynamo model. In this model, based on the toroidal flux at the base of the convection zone, bipolar magnetic regions (BMRs) are produced with statistical properties obtained from observed distributions. We find that a little quenching in BMR tilt is sufficient to stabilize the dynamo growth. The randomness and nonlinearity in the BMR emergences make the poloidal field unequal and cause some variability in the solar cycle. However, when observed scatter of BMR tilts around Joy's law with a standard deviation of $15^\circ$, is considered, our model produces a variation in the solar cycle, including north-south asymmetry comparable to the observations. The morphology of magnetic fields closely resembles observations, in particular the surface radial field possesses a more mixed polarity field. Observed scatter also produces grand minima. In 11,650 years of simulation, 17 grand minima are detected and $11\%$ of its time the model remained in these grand minima. When we double the tilt scatter, the model produces correct statistics of grand minima. Importantly, the dynamo continues even during grand minima with only a few BMRs, without requiring any additional alpha effect. The reason for this is the downward magnetic pumping which suppresses the diffusion of the magnetic flux across the surface. The magnetic pumping also helps to achieve 11-year magnetic cycle using the observed BMR flux distribution, even at the high diffusivity.

Observations suggest that the large-scale convective velocities obtained by solar convection simulations might be over-estimated (convective conundrum). One plausible solution to this could be the small-scale dynamo which cannot be fully resolved by global simulations. The small-scale Lorentz force suppresses the convective motions and also the turbulent mixing of entropy between upflows and downflows, leading to a large effective Prandtl number ($\Pr$). We explore this idea in three-dimensional global rotating convection simulations at different thermal conductivity ($\kappa$), i.e., at different $\Pr$. In agreement with previous non-rotating simulations, the convective velocity is reduced with the increase of $\Pr$ as long as the thermal conductive flux is negligible. A subadiabatic layer is formed near the base of the convection zone due to continuous deposition of low entropy plumes in low-$\kappa$ simulations. The most interesting result of our low-$\kappa$ simulations is that the convective motions are accompanied by a change in the convection structure that is increasingly influenced by small-scale plumes. These plumes tend to transport angular momentum radially inward and thus establish an anti-solar differential rotation, in striking contrast to the solar rotation profile. If such low diffusive plumes, driven by the radiative-surface cooling, are present in the Sun, then our results cast doubt on the idea that a high effective $\Pr$ may be a viable solution to the solar convective conundrum. Our study also emphasizes that any resolution of the conundrum that relies on the downward plumes must take into account angular momentum transport as well as heat transport.

In this paper, we explore the evolution of sunspot latitude distribution and explore its relations with the cycle strength. With the progress of the solar cycle, the distributions in two hemispheres from mid-latitudes propagate toward the equator and then (before the usual solar minimum) these two distributions touch each other. By visualizing the evolution of the distributions in two hemispheres, we separate the solar cycles by excluding this hemispheric overlap. From these isolated solar cycles in two hemispheres, we generate latitude distributions for each cycle, starting from cycle 8 to cycle 23. We find that the parameters of these distributions, namely, the central latitude ($C$), width ($\delta$) and height ($H$) evolve with the cycle number and they show some hemispheric asymmetries. Although the asymmetries in these parameters persist for a few successive cycles, they get corrected within a few cycles and the new asymmetries appear again. In agreement with the previous study, we find that distribution parameters are correlated with the strengths of the cycles, although these correlations are significantly different in two hemispheres. The general trend that stronger cycles begin sunspot eruptions at relatively higher latitudes and have wider bands of sunspot emergence latitudes are confirmed when combining the data from two hemispheres. We explore these features using a flux transport dynamo model with stochastic fluctuations. We find that these features are correctly reproduced in this model. The solar cycle evolution of the distribution center is also in good agreement with observations. Possible explanations of the observed features based on this dynamo model are presented.

We present results from a three-dimensional Babcock--Leighton dynamo model that is sustained by the explicit emergence and dispersal of bipolar magnetic regions (BMRs). On average, each BMR has a systematic tilt given by Joy's law. Randomness and nonlinearity in the BMR emergence of our model produce variable magnetic cycles. However, when we allow for a random scatter in the tilt angle to mimic the observed departures from Joy's law, we find more variability in the magnetic cycles. We find that the observed standard deviation in Joy's law of $\sigma_\delta = 15^\circ$ produces a variability comparable to observed solar cycle variability of $\sim $ 32%, as quantified by the sunspot number maxima between 1755--2008. We also find that tilt angle scatter can promote grand minima and grand maxima. The time spent in grand minima for $\sigma_\delta = 15^\circ$ is somewhat less than that inferred for the Sun from cosmogenic isotopes (about 9% compared to 17%). However, when we double the tilt scatter to $\sigma_\delta = 30^\circ$, the simulation statistics are comparable to the Sun ($\sim $18% of the time in grand minima and $\sim 10$% in grand maxima). Though the Babcock--Leighton mechanism is the only source of poloidal field, we find that our simulations always maintain magnetic cycles even at large fluctuations in the tilt angle. We also demonstrate that tilt quenching is a viable and efficient mechanism for dynamo saturation; a suppression of the tilt by only 1-2$^\circ$ is sufficient to limit the dynamo growth. Thus, any potential observational signatures of tilt quenching in the Sun may be subtle.

The key elements of the Babcock-Leighton dynamos are the generation of poloidal field through decay and dispersal of tilted bipolar active regions and the generation of toroidal field through the observed differential rotation. These models are traditionally known as flux transport dynamo models as the equatorward propagations of the butterfly wings in these models are produced due to an equatorward flow at the bottom of the convection zone. Here we investigate the role of downward magnetic pumping near the surface using a kinematic Babcock-Leighton model. We find that the pumping causes the poloidal field to become predominately radial in the near-surface shear layer, which allows the negative radial shear to effectively act on the radial field to produce a toroidal field. We observe a clear equatorward migration of the toroidal field at low latitudes as a consequence of the dynamo wave even when there is no meridional flow in the deep convection zone. Both the dynamo wave and the flux transport type solutions are thus able to reproduce some of the observed features of the solar cycle, including the 11-year periodicity. The main difference between the two types of solutions is the strength of the Babcock-Leighton source required to produce the dynamo action. A second consequence of the magnetic pumping is that it suppresses the diffusion of fields through the surface, which helps to allow an 11-year cycle at (moderately) larger values of magnetic diffusivity than have previously been used.

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 mechanisms 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.

The small-scale magnetic field is ubiquitous at the solar surface---even at high latitudes. From observations we know that this field is uncorrelated (or perhaps even weakly anticorrelated) with the global sunspot cycle. Our aim is to explore the origin, and particularly the cycle dependence of such a phenomenon using three-dimensional dynamo simulations. We adopt a simple model of a turbulent dynamo in a shearing box driven by helically forced turbulence. Depending on the dynamo parameters, large-scale (global) and small-scale (local) dynamos can be excited independently in this model. Based on simulations in different parameter regimes, we find that, when only the large-scale dynamo is operating in the system, the small-scale magnetic field generated through shredding and tangling of the large-scale magnetic field is positively correlated with the global magnetic cycle. However, when both dynamos are operating, the small-scale field is produced from both the small-scale dynamo and the tangling of the large-scale field. In this situation, when the large-scale field is weaker than the equipartition value of the turbulence, the small-scale field is almost uncorrelated with the large-scale magnetic cycle. On the other hand, when the large-scale field is stronger than the equipartition value, we observe an anticorrelation between the small-scale field and the large-scale magnetic cycle. This anticorrelation can be interpreted as a suppression of the small-scale dynamo. Based on our studies we conclude that the observed small-scale magnetic field in the Sun is generated by the combined mechanisms of small-scale dynamo and tangling of the large-scale field.

Using different proxies of solar activity, we have studied the following features of solar cycle. (i) A linear correlation between the amplitude of cycle and its decay rate, (ii) a linear correlation between the amplitude of cycle $n$ and the decay rate of cycle $(n - 1)$ and (iii) an anti-correlation between the amplitude of cycle $n$ and the period of cycle $(n - 1)$. Features (ii) and (iii) are very useful because they provide precursors for future cycles. We have reproduced these features using a flux transport dynamo model with stochastic fluctuations in the Babcock-Leighton $\alpha$ effect and in the meridional circulation. Only when we introduce fluctuations in meridional circulation, we are able to reproduce different observed features of solar cycle. We discuss the possible reasons for these correlations.

Nonlinear mean-field models of the solar dynamo show long-term variability, which may be relevant to different states of activity inferred from long-term radiocarbon data. This paper is aimed to probe the dynamo hysteresis predicted by the recent mean-field models of Kitchatinov \& Olemskoy (2010) with direct numerical simulations. We perform three-dimensional simulations of large-scale dynamos in a shearing box with helically forced turbulence. As initial condition, we either take a weak random magnetic field or we start from a snapshot of an earlier simulation. Two quasi-stable states are found to coexist in a certain range of parameters close to the onset of the large-scale dynamo. The simulations converge to one of these states depending on the initial conditions. When either the fractional helicity or the magnetic Prandtl number is increased between successive runs above the critical value for onset of the dynamo, the field strength jumps to a finite value. However, when the fractional helicity or the magnetic Prandtl number is then decreased again, the field strength stays at a similar value (strong field branch) even below the original onset. We also observe intermittent decaying phases away from the strong field branch close to the point where large-scale dynamo action is just possible. The dynamo hysteresis seen previously in mean-field models is thus reproduced by 3D simulations. Its possible relation to distinct modes of solar activity such as grand minima is discussed.

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. We study the effect of a dynamo-generated magnetic field on the large-scale flows, particularly on the possibility of bistable behavior of differential rotation. We solve the hydromagnetic equations numerically in a rotating spherical shell for a set of different radiative conductivities controlling the relative importance of convection. 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 get coherent single-cell meridional circulation, whereas in SL cases we get multi-cellular patterns. In both cases, we obtain poleward circulation near the surface with a magnitude close to that observed in the Sun. Moreover, both differential rotation and meridional circulation have significant magnetic cycle-related variations that are similar in strength to those of the Sun. 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.

Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the $\alpha$ effect. For stronger fields, these effects are usually suppressed or quenched, and additional anisotropies are introduced. Using different variants of the test-field method, we determine the quenching of the turbulent transport coefficients for the forced Roberts flow, isotropically forced non-helical turbulence, and rotating thermal convection. We see significant quenching only when the mean magnetic field is larger than the equipartition value of the turbulence. Expressing the magnetic field in terms of the equipartition value of the {\it quenched} flows, we obtain for the quenching exponents of the turbulent magnetic diffusivity about 1.3, 1.1, and 1.3 for Roberts flow, forced turbulence, and convection, respectively. However, when the magnetic field is expressed in terms of the equipartition value of the unquenched flows these quenching exponents become about 4, 1.5, and 2.3, respectively. For the $\alpha$ effect, the exponent is about 1.3 for the Roberts flow and 2 for convection in the first case, but 4 and 3, respectively, in the second. In convection, the quenching of turbulent pumping follows the same power law as turbulent diffusion, while for the coefficient describing the $\bf \Omega \times \bf J$ effect nearly the same quenching exponent is obtained as for $\alpha$. For forced turbulence, turbulent diffusion proportional to the second derivative along the mean magnetic field is quenched much less, especially for larger values of the magnetic Reynolds number. However, we find that in corresponding axisymmetric mean-field dynamos with dominant toroidal field the quenched diffusion coefficients are the same for the poloidal and toroidal field constituents.

The Sun has a polar magnetic field which oscillates with the 11 year sunspot cycle. This polar magnetic field is an important component of the dynamo process which is operating in the solar convection zone and produces the sunspot cycle. We have systematic direct measurements of the Sun's polar magnetic field only from about mid 1970s. There are, however, indirect proxies which give us information about this field at earlier times. The Ca K spectroheliograms taken in Kodaikanal Solar Observatory during 1904 - 2007 have now been digitized with the 4k x 4k CCD and have higher resolution (0.86 arcsec) than the other available historical datasets. From these Ca-K spectroheliograms, we have developed a completely new proxy (Polar Network Index, PNI) for the Sun's polar magnetic field. We calculate the PNI from the digitized images using an automated algorithm and calibrate our measured PNI against the polar field as measured by the Wilcox Solar Observatory for the period of 1976 - 1990. This calibration allows us to estimate polar fields for the earlier period up to 1904. The dynamo calculations done with this proxy as input data reproduce the Sun's magnetic behavior for the past century reasonably well.

We attempt to provide a quantitative theoretical explanation for the observations that Ca II H/K emission and X-ray emission from solar-like stars increase with decreasing Rossby number (i.e., with faster rotation). Assuming that these emissions are caused by magnetic cycles similar to the sunspot cycle, we construct flux transport dynamo models of $1M_{\odot}$ stars rotating with different rotation periods. We first compute the differential rotation and the meridional circulation inside these stars from a mean-field hydrodynamics model. Then these are substituted in our dynamo code to produce periodic solutions. We find that the dimensionless amplitude $f_m$ of the toroidal flux through the star increases with decreasing rotation period. The observational data can be matched if we assume the emissions to go as the power 3-4 of $f_m$. Assuming that the Babcock-Leighton mechanism saturates with increasing rotation, we can provide an explanation for the observed saturation of emission at low Rossby numbers. The main failure of our model is that it predicts an increase of magnetic cycle period with increasing rotation rate, which is the opposite of what is found observationally. Much of our calculations are based on the assumption that the magnetic buoyancy makes the magnetic flux tubes to rise radially from the bottom of the convection zone. On taking account of the fact that the Coriolis force diverts the magnetic flux tubes to rise parallel to the rotation axis in rapidly rotating stars, the results do not change qualitatively.

The solar activity cycle is successfully modeled by the flux transport dynamo, in which the meridional circulation of the Sun plays an important role. Most of the kinematic dynamo simulations assume a one-cell structure of the meridional circulation within the convection zone, with the equatorward return flow at its bottom. In view of the recent claims that the return flow occurs at a much shallower depth, we explore whether a meridional circulation with such a shallow return flow can still retain the attractive features of the flux transport dynamo (such as a proper butterfly diagram, the proper phase relation between the toroidal and poloidal fields). We consider additional cells of the meridional circulation below the shallow return flow---both the case of multiple cells radially stacked above one another and the case of more complicated cell patterns. As long as there is an equatorward flow in low latitudes at the bottom of the convection zone, we find that the solar behavior is approximately reproduced. However, if there is either no flow or a poleward flow at the bottom of the convection zone, then we cannot reproduce solar behavior. On making the turbulent diffusivity low, we still find periodic behavior, although the period of the cycle becomes unrealistically large. Also, with a low diffusivity, we do not get the observed correlation between the polar field at the sunspot minimum and the strength of the next cycle, which is reproduced when diffusivity is high. On introducing radially downward pumping, we get a more reasonable period and more solar-like behavior even with low diffusivity.

Having advanced knowledge of solar activity is important because the Sun's magnetic output governs space weather and impacts technologies reliant on space. However, the irregular nature of the solar cycle makes solar activity predictions a challenging task. This is best achieved through appropriately constrained solar dynamo simulations and as such the first step towards predictions is to understand the underlying physics of the solar dynamo mechanism. In Babcock-Leighton type dynamo models, the poloidal field is generated near the solar surface whereas the toroidal field is generated in the solar interior. Therefore a finite time is necessary for the coupling of the spatially segregated source layers of the dynamo. This time delay introduces a memory in the dynamo mechanism which allows forecasting of future solar activity. Here we discuss how this forecasting ability of the solar cycle is affected by downward turbulent pumping of magnetic flux. With significant turbulent pumping the memory of the dynamo is severely degraded and thus long term prediction of the solar cycle is not possible; only a short term prediction of the next cycle peak may be possible based on observational data assimilation at the previous cycle minimum.

The tachocline is important in the solar dynamo for the generation and the storage of the magnetic fields. A most plausible explanation for the confinement of the tachocline is given by the fast tachocline model in which the tachocline is confined by the anisotropic momentum transfer by the Maxwell stress of the dynamo generated magnetic fields. We employ a flux transport dynamo model coupled with the simple feedback formula of this fast tachocline model which basically relates the thickness of the tachocline to the Maxwell stress. We find that this nonlinear coupling not only produces a stable solar-like dynamo solution but also a significant latitudinal variation in the tachocline thickness which is in agreement with the observations.

The occurrence of grand minima like the Maunder minimum is an intriguing aspect of the sunspot cycle. We use the flux transport dynamo model to explain the grand minima, showing that they arise when either the poloidal field or the meridional circulation falls to a sufficiently low value due to fluctuations. Assuming these fluctuations to be Gaussian and determining the various parameters from the data of the last 28 cycles, we carry on a dynamo simulation with both these fluctuations. The results are remarkably close to the observational data.

We propose that grand minima in solar activity are caused by simultaneous fluctuations in the meridional circulation and the Babcock-Leighton mechanism for the poloidal field generation in the flux transport dynamo model. We present the following results: (a) fluctuations in the meridional circulation are more effective in producing grand minima; (b) both sudden and gradual initiations of grand minima are possible; (c) distributions of durations and waiting times between grand minima seem to be exponential; (d) the coherence time of the meridional circulation has an effect on the number and the average duration of grand minima, a coherence time of about 30 years being consistent with observational data. We also study the occurrence of grand maxima and find that the distributions of durations and waiting times between grand maxima are also exponential, like the grand minima. Finally we address the question whether the Babcock-Leighton mechanism can be operative during grand minima when there are no sunspots. We show that an alpha-effect restricted to the upper portions of the convection zone can pull the dynamo out of the grand minima and can match various observational requirements if the amplitude of this alpha-effect is suitably fine-tuned.

In this paper, we estimate the magnetic Reynolds number of a typical protostar before and after deuterium burning, and claim for the existence of dynamo process in both the phases, because the magnetic Reynolds number of the protostar far exceeds the critical magnetic Reynolds number for dynamo action. Using the equipartition of kinetic and magnetic energies, we estimate the steady-state magnetic field of the protostar to be of the order of kilo-gauss, which is in good agreement with observations.

Prediction of the Sun's magnetic activity is important because of its effect on space environment and climate. However, recent efforts to predict the amplitude of the solar cycle have resulted in diverging forecasts with no consensus. Yeates et al. (2008) have shown that the dynamical memory of the solar dynamo mechanism governs predictability and this memory is different for advection- and diffusion-dominated solar convection zones. By utilizing stochastically forced, kinematic dynamo simulations, we demonstrate that the inclusion of downward turbulent pumping of magnetic flux reduces the memory of both advection- and diffusion-dominated solar dynamos to only one cycle; stronger pumping degrades this memory further. Thus, our results reconcile the diverging dynamo-model-based forecasts for the amplitude of solar cycle 24. We conclude that reliable predictions for the maximum of solar activity can be made only at the preceding minimum--allowing about 5 years of advance planning for space weather. For more accurate predictions, sequential data assimilation would be necessary in forecasting models to account for the Sun's short memory.

One of the most striking aspects of the 11-year sunspot cycle is that there have been times in the past when some cycles went missing, a most well-known example of this being the Maunder minimum during 1645-1715. Analyses of cosmogenic isotopes (C14 and Be10) indicated that there were about 27 grand minima in the last 11,000 yr, implying that about 2.7% of the solar cycles had conditions appropriate for forcing the Sun into grand minima. We address the question how grand minima are produced and specifically calculate the frequency of occurrence of grand minima from a theoretical dynamo model. We assume that fluctuations in the poloidal field generation mechanism and the meridional circulation produce irregularities of sunspot cycles. Taking these fluctuations to be Gaussian and estimating the values of important parameters from the data of last 28 solar cycles, we show from our flux transport dynamo model that about 1-4% of the sunspot cycles may have conditions suitable for inducing grand minima.

The compatibility of the fast tachocline scenario with a flux transport dynamo model is explored. We employ a flux transport dynamo model coupled with simple feedback formulae relating the thickness of the tachocline to the amplitude of the magnetic field or to the Maxwell stress. The dynamo model is found to be robust against the nonlinearity introduced by this simplified fast tachocline mechanism. Solar-like butterfly diagrams are found to persist and, even without any parameter fitting, the overall thickness of the tachocline is well within the range admitted by helioseismic constraints. In the most realistic case of a time and latitude dependent tachocline thickness linked to the value of the Maxwell stress, both the thickness and its latitude dependence are in excellent agreement with seismic results. In the nonparametric models, cycle related temporal variations in tachocline thickness are somewhat larger than admitted by helioseismic constraints; we find, however, that introducing a further parameter into our feedback formula readily allows further fine tuning of the thickness variations.

The sunspot number varies roughly periodically with time. However the individual cycle durations and the amplitudes are found to vary in an irregular manner. It is observed that the stronger cycles are having shorter rise times and vice versa. This leads to an important effect know as the Waldmeier effect. Another important feature of the solar cycle irregularity are the grand minima during which the activity level is strongly reduced. We explore whether these solar cycle irregularities can be studied with the help of the flux transport dynamo model of the solar cycle. We show that with a suitable stochastic fluctuations in a regular dynamo model, we are able to reproduce many irregular features of the solar cycle including the Waldmeier effect and the grand minimum. However, we get all these results only if the value of the turbulent diffusivity in the convection zone is reasonably high.

Guided by the recent observational result that the meridional circulation of the Sun becomes weaker at the time of the sunspot maximum, we have included a parametric quenching of the meridional circulation in solar dynamo models such that the meridional circulation becomes weaker when the magnetic field at the base of the convection zone is stronger. We find that a flux transport solar dynamo tends to become unstable on including this quenching of meridional circulation if the diffusivity in the convection zone is less than about 2 * 10^{11} cm^2/s. The quenching of alpha, however, has a stabilizing effect and it is possible to stabilize a dynamo with low diffusivity with sufficiently strong alpha-quenching. For dynamo models with high diffusivity, the quenching of meridional circulation does not produce a large effect and the dynamo remains stable. We present a solar-like solution from a dynamo model with diffusivity 2.8 * 10^{12} cm^2/s in which the quenching of meridional circulation makes the meridional circulation vary periodically with solar cycle as observed and does not have any other significant effect on the dynamo.

We explore the importance of meridional circulation variations in modelling the irregularities of the solar cycle by using the flux transport dynamo model. We show that a fluctuating meridional circulation can reproduce some features of the solar cycle like the Waldmeier effect and the grand minimum. However, we get all these results only if the value of the turbulent diffusivity in the convection zone is reasonably high.

We discuss one important aspect of Waldmeier effect which says that the stronger cycles rise rapidly than weaker cycles. We studied four different data set of solar activity indices, and find strong linear correlation between rise rates and amplitudes of solar activity. We study this effect theoretically by introducing suitable stochastic fluctuations in our regular solar dynamo model.

We propose that at the beginning of the Maunder minimum the poloidal field or amplitude of meridional circulation or both fell abruptly to low values. With this proposition, a flux transport dynamo model is able to reproduce various important aspects of the historical records of the Maunder minimum remarkably well.

Meridional circulation is an important ingredient in flux transport dynamo models. We have studied its importance on the period, the amplitude of the solar cycle, and also in producing Maunder-like grand minima in these models. First, we model the periods of the last 23 sunspot cycles by varying the meridional circulation speed. If the dynamo is in a diffusion-dominated regime, then we find that most of the cycle amplitudes also get modeled up to some extent when we model the periods. Next, we propose that at the beginning of the Maunder minimum the amplitude of meridional circulation dropped to a low value and then after a few years it increased again. Several independent studies also favor this assumption. With this assumption, a diffusion-dominated dynamo is able to reproduce many important features of the Maunder minimum remarkably well. If the dynamo is in a diffusion-dominated regime, then a slower meridional circulation means that the poloidal field gets more time to diffuse during its transport through the convection zone, making the dynamo weaker. This consequence helps to model both the cycle amplitudes and the Maunder-like minima. We, however, fail to reproduce these results if the dynamo is in an advection-dominated regime.

We discuss two aspects of the Waldmeier Effect, namely (1) the rise times of sunspot cycles are anti-correlated to their strengths (WE1) and (2) the rates of rise of the cycles are correlated to their strengths (WE2). From analysis of four different data sets we conclude that both WE1 and WE2 exist in all the data sets. We study these effects theoretically by introducing suitable stochastic fluctuations in our regular solar dynamo model.

We confirm that the evidence for the Waldmeier effect WE1 (the anti-correlation between rise times of sunspot cycles and their strengths) and the related effect WE2 (the correlation between rise rates of cycles and their strengths) is found in different kinds of sunspot data. We explore whether these effects can be explained theoretically on the basis of the flux transport dynamo models of sunspot cycles. Two sources of irregularities of sunspot cycles are included in our model: fluctuations in the poloidal field generation process and fluctuations in the meridional circulation. We find WE2 to be a robust result which is produced in different kinds of theoretical models for different sources of irregularities. The Waldmeier effect WE1, on the other hand, arises from fluctuations in the meridional circulation and is found only in the theoretical models with reasonably high turbulent diffusivity which ensures that the diffusion time is not more than a few years.

The accretion disk around a compact object is a nonlinear general relativistic system involving magnetohydrodynamics. Naturally the question arises whether such a system is chaotic (deterministic) or stochastic (random) which might be related to the associated transport properties whose origin is still not confirmed. Earlier, the black hole system GRS 1915+105 was shown to be low dimensional chaos in certain temporal classes. However, so far such nonlinear phenomena have not been studied fairly well for neutron stars which are unique for their magnetosphere and kHz quasi-periodic oscillation (QPO). On the other hand, it was argued that the QPO is a result of nonlinear magnetohydrodynamic effects in accretion disks. If a neutron star exhibits chaotic signature, then what is the chaotic/correlation dimension? We analyze RXTE/PCA data of neutron stars Sco X-1 and Cyg X-2, along with the black hole Cyg X-1 and the unknown source Cyg X-3, and show that while Sco X-1 and Cyg X-2 are low dimensional chaotic systems, Cyg X-1 and Cyg X-3 are stochastic sources. Based on our analysis, we argue that Cyg X-3 may be a black hole.

We propose that the poloidal field at the end of the last sunspot cycle before the Maunder minimum fell to a very low value due to fluctuations in the Babcock--Leighton process. With this assumption, a flux transport dynamo model is able to explain various aspects of the historical records of the Maunder minimum remarkably well on choosing the parameters of the model suitably to give the correct growth time.