Here, we study the dichotomy of the escaping atmosphere of the newly discovered close-in exoplanet AU Mic b. On one hand, the high EUV stellar flux is expected to cause a strong atmospheric escape in AU Mic b. On the other hand, the wind of this young star is believed to be very strong, which could reduce or even inhibit the planet's atmospheric escape. AU Mic is thought to have a wind mass-loss rate that is up to $1000$ times larger than the solar wind mass-loss rate ($\dot{M}_\odot$). To investigate this dichotomy, we perform 3D hydrodynamics simulations of the stellar wind--planetary atmosphere interactions in the AU Mic system and predict the synthetic Ly-$\alpha$ transits of AU Mic b. We systematically vary the stellar wind mass-loss rate from a `no wind' scenario to up to a stellar wind with a mass-loss rate of $1000~\dot{M}_\odot$. We find that, as the stellar wind becomes stronger, the planetary evaporation rate decreases from $6.5\times 10^{10}$ g/s to half this value. With a stronger stellar wind, the atmosphere is forced to occupy a smaller volume, affecting transit signatures. Our predicted Ly-$\alpha$ absorption drops from $\sim 20\%$, in the case of `no wind' to barely any Ly-$\alpha$ absorption in the extreme stellar wind scenario. Future Ly-$\alpha$ transits could therefore place constraints not only on the evaporation rate of AU Mic b, but also on the mass-loss rate of its host star.

Stellar high-energy radiation (X-ray and extreme ultraviolet, XUV) drives atmospheric escape in close-in exoplanets. Given that stellar irradiation depends on the stellar magnetism and that stars have magnetic cycles, we investigate how cycles affect the evolution of exoplanetary atmospheric escape. Firstly, we consider a hypothetical HD209458b-like planet orbiting the Sun. For that, we implement the observed solar XUV radiation available over one and a half solar cycles in a 1D hydrodynamic escape model of HD209458b. We find that atmospheric escape rates show a cyclic variation (from 7.6 to 18.5 $\times$ 10$^{10}$ g s$^{-1}$), almost proportional to the incident stellar radiation. To compare this with observations, we compute spectroscopic transits in two hydrogen lines. We find non-detectable cyclic variations in Ly$\alpha$ transits. Given the temperature sensitiveness of the H$\alpha$ line, its equivalent width has an amplitude of 1.9 mA variation over the cycle, which could be detectable in exoplanets such as HD209458b. We demonstrate that the XUV flux is linearly proportional to the magnetic flux during the solar cycle. Secondly, we apply this relation to derive the cyclic evolution of the XUV flux of HD189733 using the available magnetic flux observations of the star from Zeeman Doppler Imaging over nearly a decade. The XUV fluxes are then used to model escape in HD189733b, which shows escape rate varying from 2.8 to 6.5 $\times$ 10$^{10}$ g s$^{-1}$. Like in the HD209458b case, this introduces variations in Ly$\alpha$ and H$\alpha$ transits, with H$\alpha$ variations more likely to be observable. Finally, we show that a strong stellar flare would enhance significantly Ly$\alpha$ and H$\alpha$ transit depths.

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.

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.

The observed convective flows on the photosphere (e.g., supergranulation, granulation) play a key role in the Babcock-Leighton (BL) process to generate large-scale polar fields from sunspots fields. In most surface flux transport (SFT) and BL dynamo models, the dispersal and migration of surface fields are modeled as an effective turbulent diffusion. We present the first kinematic 3D FT/BL model to explicitly incorporate realistic convective flows based on solar observations. The results obtained are generally in good agreement with the observed surface flux evolution and with non-convective models that have a turbulent diffusivity on the order of $3 \times 10^{12}$ cm$^2$ s$^{-1}$ (300 km$^2$ s$^{-1}$). However, we find that the use of a turbulent diffusivity underestimates the dynamo efficiency, producing weaker mean fields and shorter cycle.

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.

In this thesis, various studies are performed leading to better understanding of the 11-year solar cycle and its theoretical modeling with the flux transport dynamo model. Although this is primarily a theoretical thesis, there is a part dealing with the analysis of observational data. The various proxies of solar activity from various observatory including the sunspot area records of Kodaikanal Observatory have been analyzed to study the irregular aspects of solar cycles and an analysis has been carried out on the correlation between the decay rate and the next cycle amplitude. Theoretical analysis starts with explaining how the magnetic buoyancy has been treated in the flux transport dynamo models, and advantages and disadvantages of different treatments. It is found that some of the irregular properties of the solar cycle in the decaying phase can only be well explained using a particular treatment of the magnetic buoyancy. Next, the behavior of the dynamo with different spatial structures of the meridional flow based on recent helioseismology results has been studied. A theoretical model is constructed considering the back reaction due to the Lorentz force on the meridional flows which explains the observed variation of the meridional flow with the solar cycle. Finally, some results with 3D FTD models are presented. This 3D model is developed to handle the Babcock-Leighton mechanism and magnetic buoyancy more realistically than previous 2D models and can capture some important effects connected with the subduction of the magnetic field in polar regions, which are missed in 2D surface flux transport models. This 3D model is further used to study the evolution of the magnetic fields due to a turbulent non-axisymmetric velocity field and to compare the results with the results obtained by using a simple turbulent diffusivity coefficient.

The observed convective flows on the photosphere (e.g., supergranulation, granulation) play a key role in the Babcock-Leighton (BL) process to generate large-scale polar fields from sunspots fields. In most surface flux transport (SFT) and BL dynamo models, the dispersal and migration of surface fields is modeled as an effective turbulent diffusion. Recent SFT models have incorporated explicit, realistic convective flows in order to improve the fidelity of convective transport but, to our knowledge, this has not yet been implemented in previous BL models. Since most Flux-Transport (FT)/BL models are axisymmetric, they do not have the capacity to include such flows. We present the first kinematic 3D FT/BL model to explicitly incorporate realistic convective flows based on solar observations. Though we describe a means to generalize these flows to 3D, we find that the kinematic small-scale dynamo action they produce disrupts the operation of the cyclic dynamo. Cyclic solution is found by limiting the convective flow to act only on the vertical radial component of the magnetic field. The results obtained are generally in good agreement with the observed surface flux evolution and with non-convective models that have a turbulent diffusivity on the order of $3 \times 10^{12}$ cm$^2$ s$^{-1}$ (300 km$^2$ s$^{-1}$). However, we find that the use of a turbulent diffusivity underestimates the dynamo efficiency, producing weaker mean fields and shorter cycle than in the convective models. Also, the convective models exhibit mixed polarity bands in the polar regions that have no counterpart in solar observations. Also, the explicitly computed turbulent electromotive force (emf) bears little resemblance to a diffusive flux. We also find that the poleward migration speed of poloidal flux is determined mainly by the meridional flow and the vertical diffusion.

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.

Long-term sunspot observations are key to understand and predict the solar activities and its effects on the space weather.Consistent observations which are crucial for long-term variations studies,are generally not available due to upgradation/modifications of observatories over the course of time. We present the data for a period of 90 years acquired from persistent observation at the Kodaikanal observatory in India. We use an advanced semi-automated algorithm to detect the sunspots form each calibrated white-light image. Area, longitude and latitude of each of the detected sunspots are derived. Implementation of a semi-automated method is very necessary in such studies as it minimizes the human bias in the detection procedure. Daily, monthly and yearly sunspot area variations obtained from the Kodaikanal, compared well with the Greenwich sunspot area data. We find an exponentially decaying distribution for the individual sunspot area for each of the solar cycles. Analyzing the histograms of the latitudinal distribution of the detected sunspots, we find Gaussian distributions, in both the hemispheres, with the centers at $\sim$15$^{\circ}$ latitude. The height of the Gaussian distributions are different for the two hemispheres for a particular cycle. Using our data, we show clear presence of Waldmeier effect which correlates the rise time with the cycle amplitude. Using the wavelet analysis, we explored different periodicities of different time scales present in the sunspot area times series.

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.

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.