A new formula for predicting solar cycles based on the current theoretical understanding of the solar cycle from flux transport dynamo is presented. Two important processes---fluctuations in the Babcock-Leighton mechanism and variations in the meridional circulation, which are believed to be responsible for irregularities of the solar cycle---are constrained by using observational data. We take the polar field near minima of the cycle as a measure of the randomness in the Babcock-Leighton process, and the decay rate near the minima as a consequence of the change in meridional circulation. We couple these two observationally derived quantities into a single formula to predict the amplitude of the future solar cycle. Our new formula suggests that the cycle 25 would be a moderate cycle. Whether this formula for predicting the future solar cycle can be justified theoretically is also discussed using simulations from the flux transport dynamo model.

Although the Saha ionization equation is a standard topic in advanced statistical physics and most professional physicists would have some knowledge of it, the exact nature of Meghnad Saha's contributions in this subject is not widely known. Based on an analysis of Saha's original papers and other relevant papers by contemporary scientists, as well as other source materials such as letters exchanged between scientists, I discuss how Saha used the theory of thermal ionization to solve some of the most important astrophysics problems of that era.

The polar magnetic field of the Sun is a manifestation of certain aspects of the dynamo process and is a good precursor for predicting a sunspot cycle before its onset. Although actual synoptic measurements of this field exist only from the mid-1970s, it has now been possible to determine its evolution from the beginning of the twentieth century with the help of various proxies. The recently developed 3D kinematic dynamo model can study the build-up of the Sun's polar magnetic field more realistically than the earlier surface flux transport model.

The meridional circulation of the Sun is observationally found to vary with the solar cycle, becoming slower during the solar maxima. We explain this by constructing a theoretical model in which the equation of the meridional circulation (the $\phi$ component of the vorticity equation) is coupled with the equations of the flux transport dynamo model. We find that the Lorentz force of the dynamo-generated magnetic fields can slow down the \MC\ during the solar maxima in broad conformity with the observations.

Several phenomena connected with the magnetic field of the Sun (the cool sunspots, the hot corona, solar flares, the solar wind) are collectively known as solar activity. This paper discusses how one uses the MHD equations to understand how the magnetic field of the Sun is produced by the dynamo process and then gives rise to these diverse activities, making the Sun the best laboratory for plasma physics in the limit of high magnetic Reynolds number (defined at the end of the Introduction).

This is the talk I gave at the workshop "Solar-Stellar Magnetism: Past, Present and Future" held in Jaipur on 18 February 2018 on the occasion of my 60th year.

The flux transport dynamo, in which the poloidal magnetic field is generated by the Babcock--Leighton mechanism and the meridional circulation plays a crucial role, has emerged as an attractive model for the solar cycle. Based on theoretical calculations done with this model, we argue that the fluctuations in the Babcock--Leighton mechanism and the fluctuations in the meridional circulation are the most likely causes of the irregularities of the solar cycle. With our increased theoretical understanding of how these irregularities arise, it can be possible to predict a future solar cycle by feeding the appropriate observational data in a theoretical dynamo model.

We begin with a review of the predictions for cycle~24 before its onset. After summarizing the basics of the flux transport dynamo model, we discuss how this model had been used to make a successful prediction of cycle~24, on the assumption that the irregularities of the solar cycle arise due to the fluctuations in the Babcock--Leighton mechanism. We point out that fluctuations in the meridional circulation can be another cause of irregularities in the cycle.

Observations of the meridional circulation of the Sun, which plays a key role in the operation of the solar dynamo, indicate that its speed varies with the solar cycle, becoming faster during the solar minima and slower during the solar maxima. To explain this variation of the meridional circulation with the solar cycle, we construct a theoretical model by coupling the equation of the meridional circulation (the $\phi$ component of the vorticity equation within the solar convection zone) with the equations of the flux transport dynamo model. We consider the back reaction due to the Lorentz force of the dynamo-generated magnetic fields and study the perturbations produced in the meridional circulation due to it. This enables us to model the variations of the meridional circulation without developing a full theory of the meridional circulation itself. We obtain results which reproduce the observational data of solar cycle variations of the meridional circulation reasonably well. We get the best results on assuming the turbulent viscosity acting on the velocity field to be comparable to the magnetic diffusivity (i.e. on assuming the magnetic Prandtl number to be close to unity). We have to assume an appropriate bottom boundary condition to ensure that the Lorentz force cannot drive a flow in the subadiabatic layers below the bottom of the tachocline. Our results are sensitive to this bottom boundary condition. We also suggest a hypothesis how the observed inward flow towards the active regions may be produced.

We develop a three-dimensional kinematic self-sustaining model of the solar dynamo in which the poloidal field generation is from tilted bipolar sunspot pairs placed on the solar surface above regions of strong toroidal field by using the SpotMaker algorithm, and then the transport of this poloidal field to the tachocline is primarily caused by turbulent diffusion. We obtain a dipolar solution within a certain range of parameters. We use this model to study the build-up of the polar magnetic field and show that some insights obtained from surface flux transport (SFT) models have to be revised. We present results obtained by putting a single bipolar sunspot pair in a hemisphere and two symmetrical sunspot pairs in two hemispheres.We find that the polar fields produced by them disappear due to the upward advection of poloidal flux at low latitudes, which emerges as oppositely-signed radial flux and which is then advected poleward by the meridional flow. We also study the effect that a large sunspot pair, violating Hale's polarity law would have on the polar field. We find that there would be some effect---especially if the anti-Hale pair appears at high latitudes in the mid-phase of the cycle---though the effect is not very dramatic.

In this review, we discuss whether the present solar dynamo models can be extrapolated to explain various aspects of stellar activity. We begin with a summary of the following kinds of data for solar-like stars: (i) data pertaining to stellar cycles from Ca H/K emission over many years; (ii) X-ray data indicating hot coronal activity; (iii) starspot data (especially about giant polar spots); and (iv) data pertaining to stellar superflares. Then we describe the current status of solar dynamo modelling---giving an introduction to the flux transport dynamo model, the currently favoured model for the solar cycle. While an extrapolation of this model to solar-like stars can explain some aspects of observational data, some other aspects of the data still remain to be theoretically explained. It is not clear right now whether we need a different kind of dynamo mechanism for stars having giant starspots or producing very strong superflares.

One important ingredient of flux transport dynamo models is the rise of the toroidal magnetic field through the convection zone due to magnetic buoyancy to produce bipolar sunspots and then the generation of the poloidal magnetic field from these bipolar sunspots due to the Babcock-Leighton mechanism. Over the years, two methods of treating magnetic buoyancy, a local method and a non-local method have been used widely by different groups in constructing 2D kinematic models of the flux transport dynamo. We review both these methods and conclude that neither of them is fully satisfactory, presumably because magnetic buoyancy is an inherently 3D process. We also point out so far we do not have proper understanding of why sunspot emergence is restricted to rather low latitudes.

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.

We first discuss how the flux transport dynamo with reasonably high diffusion can explain both the regular and the irregular features of the solar cycle quite well. Then we critically examine the inadequacies of the model and the challenge posed by some recent observational data about meridional circulation, arriving at the conclusion that this model can still work within the bounds of observational data

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.

The 11-year sunspot cycle has many irregularities, the most promi- nent amongst them being the grand minima when sunspots may not be seen for several cycles. After summarizing the relevant observational data about the irregularities, we introduce the flux transport dynamo model, the currently most successful theoretical model for explaining the 11-year sunspot cycle. Then we analyze the respective roles of nonlinearities and random fluctuations in creating the irregularities. We also discuss how it has recently been realized that the fluctuations in meridional circula- tion also can be a source of irregularities. We end by pointing out that fluctuations in the poloidal field generation and fluctuations in meridional circulation together can explain the occurrences of grand minima.

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.

Parameters of a special kind of \alpha-effect known in dynamo theory as the Babcock-Leighton mechanism are estimated using the data of sunspot catalogs. The estimates evidence the presence of the Babcock-Leighton \alpha-effect on the Sun. Fluctuations of the \alpha-effect are also estimated. The fluctuation amplitude appreciably exceeds the mean value, and the characteristic time for the fluctuations is comparable to the period of the solar rotation. Fluctuations with the parameters found are included in a numerical model for the solar dynamo. Computations show irregular changes in the amplitudes of the magnetic cycles on time scales of centuries and millennia. The calculated statistical characteristics of the grand solar minima and maxima agree with the data on solar activity over the Holocene.

Recent observations of solar type stars with the Kepler satellite by Maehara et al. have revealed the existence of superflares (with energy of 10^33 - 10^35 erg) on Sun-like stars, which are similar to our Sun in their surface temperature (5600 K - 6000 K) and slow rotation (rotational period > 10 days). From the statistical analysis of these superflares, it was found that superflares with energy 10^34 erg occur once in 800 years and superflares with 10^35 erg occur once in 5000 years on Sun-like stars. In this paper, we examine whether superflares with energy of 10^33 - 10^35 erg could occur on the present Sun through the use of simple order-of-magnitude estimates based on current ideas relating to the mechanisms of the solar dynamo.

We point out the difficulties in carrying out direct numerical simulation of the solar dynamo problem and argue that kinematic mean-field models are our best theoretical tools at present for explaining various aspects of the solar cycle in detail. The most promising kinematic mean-field model is the flux transport dynamo model, in which the toroidal field is produced by differential rotation in the tachocline, the poloidal field is produced by the Babcock--Leighton mechanism at the solar surface and the meridional circulation plays a crucial role. Depending on whether the diffusivity is high or low, either the diffusivity or the meridional circulations provides the main transport mechanism for the poloidal field to reach the bottom of the convection zone from the top. We point out that the high-diffusivity flux transport dynamo model is consistent with various aspects of observational data. The irregularities of the solar cycle are primarily produced by fluctuations in the Babcock--Leighton mechanism and in the meridional circulation. We summarize recent work on the fluctuations of meridional circulation in the flux transport dynamo, leading to explanations of such things as the Waldmeier effect.

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.

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.

The periodically varying Lorentz force of the periodic solar magnetic field generated by the solar dynamo can induce two kinds of motions: torsional oscillations and periodic variations in the meridional circulation. Observational evidence now exists for both these kinds of motions. We discuss our ongoing effort in theoretically studying the variations of the meridional circulation. Then we present our theoretical model of torsional oscillations, which addresses the question why these oscillations start before sunspot cycles at latitudes higher than where sunspots are seen.

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.

Since a universally accepted dynamo model of grand minima does not exist at the present time, we concentrate on the physical processes which may be behind the grand minima. After summarizing the relevant observational data, we make the point that, while the usual sources of irregularities of solar cycles may be sufficient to cause a grand minimum, the solar dynamo has to operate somewhat differently from the normal to bring the Sun out of the grand minimum. We then consider three possible sources of irregularities in the solar dynamo: (i) nonlinear effects; (ii) fluctuations in the poloidal field generation process; (iii) fluctuations in the meridional circulation. We conclude that (i) is unlikely to be the cause behind grand minima, but a combination of (ii) and (iii) may cause them. If fluctuations make the poloidal field fall much below the average or make the meridional circulation significantly weaker, then the Sun may be pushed into a grand minimum.

After summarizing the relevant observational data, we discuss how a study of flux tube dynamics in the solar convection zone helps us to understand the formation of sunspots. Then we introduce the flux transport dynamo model and assess its success in modelling both the solar cycle and its departures from strictly periodic behaviour.

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.

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.

The most promising model for explaining the origin of solar magnetism is the flux transport dynamo model, in which the toroidal field is produced by differential rotation in the tachocline, the poloidal field is produced by the Babcock--Leighton mechanism at the solar surface and the meridional circulation plays a crucial role. After discussing how this model explains the regular periodic features of the solar cycle, we come to the questions of what causes irregularities of solar cycles and whether we can predict future cycles. Only if the diffusivity within the convection zone is sufficiently high, the polar field at the sunspot minimum is correlated with strength of the next cycle. This is in conformity with the limited available observational data.

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.

Assuming that the torsional oscillation is driven by the Lorentz force of the magnetic field associated with the sunspot cycle, we use a flux transport dynamo to model it and explain its initiation at a high latitude before the beginning of the sunspot cycle.

Although the Sun's torsional oscillation is believed to be driven by the Lorentz force associated with the sunspot cycle, this oscillation begins 2--3 yr before the sunspot cycle. We provide a theoretical explanation of this with the help of a solar dynamo model having a meridional circulation penetrating slightly below the bottom of the convection zone, because only in such dynamo models the strong toroidal field forms a few years before the sunspot cycle and at a higher latitude.

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.

We believe the Babcock--Leighton process of poloidal field generation to be the main source of irregularity in the solar cycle. The random nature of this process may make the poloidal field in one hemisphere stronger than that in the other hemisphere at the end of a cycle. We expect this to induce an asymmetry in the next sunspot cycle. We look for evidence of this in the observational data and then model it theoretically with our dynamo code. Since actual polar field measurements exist only from 1970s, we use the polar faculae number data recorded by Sheeley (1991) as a proxy of the polar field and estimate the hemispheric asymmetry of the polar field in different solar minima during the major part of the twentieth century. This asymmetry is found to have a reasonable correlation with the asymmetry of the next cycle. We then run our dynamo code by feeding information about this asymmetry at the successive minima and compare with observational data. We find that the theoretically computed asymmetries of different cycles compare favourably with the observational data, the correlation coefficient being 0.73. Due to the coupling between the two hemispheres, any hemispheric asymmetry tends to get attenuated with time. The hemispheric asymmetry of a cycle either from observational data or from theoretical calculation statistically tends to be less than the asymmetry in the polar field (as inferred from the faculae data) in the preceding minimum. This reduction factor turns out to be 0.38 and 0.60 respectively in observational data and theoretical simulation.

Although systematic measurements of the solar polar magnetic field exist only from mid 1970s, other proxies can be used to infer the polar field at earlier times. The observational data indicate a strong correlation between the polar field at a sunspot minimum and the strength of the next cycle, although the strength of the cycle is not correlated well with the polar field produced at its end. This suggests that the Babcock Leighton mechanism of poloidal field generation from decaying sunspots involves randomness, whereas the other aspects of the dynamo process must be reasonably ordered and deterministic. Only if the magnetic diffusivity within the convection zone is assumed to be high, we can explain the correlation between the polar field at a minimum and the next cycle. We give several independent arguments that the diffusivity must be of this order. In a dynamo model with diffusivity like this, the poloidal field generated at the mid latitudes is advected toward the poles by the meridional circulation and simultaneously diffuses towards the tachocline, where the toroidal field for the next cycle is produced. To model actual solar cycles with a dynamo model having such high diffusivity, we have to feed the observational data of the poloidal field at the minimum into the theoretical model. We develop a method of doing this in a systematic way. Our model predicts that cycle 24 will be a very weak cycle. Hemispheric asymmetry of solar activity is also calculated with our model and compared with observational data.

The mean value of the normalized current helicity in solar active regions is on the order of 1e-8 1/m, negative in the northern hemisphere, positive in the southern hemisphere. Observations indicate that this helicity has a subsurface origin. Possible mechanisms leading to a twist of this amplitude in magnetic flux tubes include the solar dynamo, convective buffeting of rising flux tubes, and the accretion of weak external poloidal flux by a rising toroidal flux tube. After briefly reviewing the observational and theoretical constraints on the origin of helicity, we present a recently developed detailed model for poloidal flux accretion.

Whether the upcoming cycle 24 of solar activity will be strong or not is being hotly debated. The solar cycle is produced by a complex dynamo mechanism. We model the last few solar cycles by `feeding' observational data of the Sun's polar magnetic field into our solar dynamo model. Our results fit the observed sunspot numbers of cycles 21-23 extremely well and predict that cycle~24 will be about 35% weaker than cycle~23.

Following an earlier proposal by Choudhuri (2003) for the origin of twist in the magnetic fields of solar active regions, we model the penetration of a wrapped up background poloidal field into a toroidal magnetic flux tube rising through the solar convective zone. The rise of the straight, cylindrical flux tube is followed by numerically solving the induction equation in a comoving Lagrangian frame, while an external poloidal magnetic field is assumed to be radially advected onto the tube with a speed corresponding to the rise velocity. One prediction of our model is the existence of a ring of reverse current helicity on the periphery of active regions. On the other hand, the amplitude of the resulting twist depends sensitively on the assumed structure (diffuse vs. concentrated/intermittent) of the active region magnetic field right before its emergence, and on the assumed vertical profile of the poloidal field. Nevertheless, in the model with the most plausible choice of assumptions a mean twist comparable to the observations can result rather naturally. Our results indicate that the contribution of this mechanism to the twist can be quite significant, and under favourable circumstances it can potentially account for most of the current helicity observed in active regions.

We present here our response to Dikpati et al.'s criticism of our recent solar dynamo model.

We present calculations of helicity based on our solar dynamo model and show that the results are consistent with observational data.

We explore a two-dimensional kinematic solar dynamo model in a full sphere, based on the helioseismically determined solar rotation profile and with an $\alpha$ effect concentrated near the solar surface, which captures the Babcock--Leighton idea that the poloidal field is created from the decay of tilted bipolar active regions. The meridional circulation, assumed to penetrate slightly below the tachocline, plays an important role. Some doubts have recently been raised regarding the ability of such a model to reproduce solar-like dipolar parity. We specifically address the parity issue and show that the dipolar mode is preferred when certain reasonable conditions are satisfied, the most important condition being the requirement that the poloidal field should diffuse efficiently to get coupled across the equator. Our model is shown to reproduce various aspects of observational data, including the phase relation between sunspots and the weak, diffuse field.

Millisecond pulsars, with magnetic fields weaker by three to four orders compared to those of ordinary pulsars, are presumed to be neutron stars spun up by binary accretion. We expect the magnetic field to get screened by the accreted material. Our simulation of this screening mechanism shows, for the first time, that the field decreases by a purely geometric factor $\sin^{-7/2} \theta_{\rm P,i}$ before freezing to an asymptotic value, where $ \theta_{\rm P,i}$ is the initial angular width of the polar cap. If $ \theta_{\rm P,i}$ lies in the range $5^o$--$10^o$, then the magnetic field diminution factor turns out to be $\sim 10^3$--$10^4$ in conformity with observational data.

In response to the claim by Dziembowski et al. (2001) that the solar radius decreases with magnetic activity at the rate of 1.5 km/yr, we consider the theoretical question whether a radius variation is expected with the solar cycle. If the radius variation is caused by the magnetic pressure of toroidal flux tubes at the bottom of the convection zone, then the dynamo model of Nandy and Choudhuri (2002) would suggest a radius decrease with magnetic activity, in contrast to other dynamo models which would suggest a radius increase. However, the radius decrease is estimated to be only of the order of hundreds of metres.

A possible mechanism for screening of the surface magnetic field of an accreting neutron star, by the accreted material, is investigated. In particular, we investigate the nature of the evolution of the internal field configuration in the case of a) a polar cap accretion and b) a spherical accretion.

Mean field dynamo theory deals with various mean quantities and does not directly throw any light on the question of existence of flux tubes. We can, however, draw important conclusions about flux tubes in the interior of the Sun by combining additional arguments with the insights gained from solar dynamo solutions. The polar magnetic field of the Sun is of order 10 G, whereas the toroidal magnetic field at the bottom of the convection zone has been estimated to be 100,000 G. Simple order-of-magnitude estimates show that the shear in the tachocline is not sufficient to stretch a 10 G mean radial field into a 100,000 G mean toroidal field. We argue that the polar field of the Sun must get concentrated into intermittent flux tubes before it is advected to the tachocline. We estimate the strengths and filling factors of these flux tubes. Stretching by shear in the tachocline is then expected to produce a highly intermittent magnetic configuration at the bottom of the convection zone. The meridional flow at the bottom of the convection zone should be able to carry this intermittent magnetic field equatorward, as suggested recently by Nandy and Choudhuri (2002). When a flux tube from the bottom of the convection zone rises to a region of pre-existing poloidal field at the surface, we point out that it picks up a twist in accordance with the observations of current helicities at the solar surface.

A possible mechanism for screening of the surface magnetic field of an accreting neutron star, by the accreted material, is investigated. We model the material flow in the surface layers of the star by an assumed 2-D velocity field satisfying all the physical requirements. Using this model velocity we find that, in the absence of magnetic buoyancy, the surface field is screened (submergence of the field by advection) within the time scale of material flow of the top layers. On the other hand, if magnetic buoyancy is present, the screening happens over a time scale which is characteristic of the slower flow of the deeper (hence, denser) layers. For accreting neutron stars, this longer time scale turns out to be about $10^5$ years, which is of the similar order as the accretion time scale of most massive X-ray binaries.

We develop a model of the solar dynamo in which, on the one hand, we follow the Babcock-Leighton approach to include surface processes like the production of poloidal field from the decay of active regions, and, on the other hand, we attempt to develop a mean field theory that can be studied in quantitative detail. One of the main challenges in developing such models is to treat the buoyant rise of toroidal field and the production of poloidal field from it near the surface. We build up a dynamo model with two contrasting methods of treating buoyancy. In one method, we incorporate the generation of the poloidal field near the solar surface by Durney's procedure of double ring eruption. In the second method, the poloidal field generation is treated by a positive alpha-effect concentrated near the solar surface, coupled with an algorithm for handling buoyancy. The two methods are found to give qualitatively similar results.