The rate of magnetic field diffusion plays an essential role in several astrophysical plasma processes. It has been demonstrated that the omnipresent turbulence in astrophysical media induces fast magnetic reconnection, which consequently leads to large-scale magnetic flux diffusion at a rate independent of the plasma microphysics. This process is called "reconnection diffusion" (RD) and allows for the diffusion of fields which are dynamically important. The current theory describing RD is based on incompressible magnetohydrodynamic (MHD) turbulence. In this work, we have tested quantitatively the predictions of the RD theory when magnetic forces are dominant in the turbulence dynamics (Alfv\'{e}nic Mach number $M_A < 1$). We employed the \textsc{Pencil Code} to perform numerical simulations of forced MHD turbulence, extracting the values of the diffusion coefficient $\eta_{RD}$ using the Test-Field method. Our results are consistent with the RD theory ($\eta_{RD} \sim M_A^{3}$ for $M_A < 1$) when turbulence approaches the incompressible limit (sonic Mach number $M_S \lesssim 0.02$), while for larger $M_S$ the diffusion is faster ($\eta_{RD} \sim M_A^{2}$).This work shows for the first time simulations of compressible MHD turbulence with the suppression of the cascade in the direction parallel to the mean magnetic field, which is consistent with incompressible weak turbulence theory. We also verified that in our simulations the energy cascading time does not follow the scaling with $M_A$ predicted for the weak regime, in contradiction with the RD theory assumption. Our results generally support and expand the RD theory predictions.

One of the most enigmatic phenomena related to solar activity is the so-called Maunder minimum phase. It consists of the lowest sunspot's counting ever registered for the Sun, and never so far confirmed for other stars. Since the spectroscopic observations of stellar activity at the Mount Wilson Observatory, the solar analog HD43587 has shown a very low and apparently invariant activity level, which makes it a Maunder minimum candidate. We analyze the chromospheric activity evolution of HD43587 and its evolutive status, intending to unravel the reasons for this low and flat activity. We use an activity measurements dataset available in the literature, and compute the activity S-index from HARPS and NARVAL spectra, to infer a cycle period. Besides, we analyze the CoRoT light-curve of HD43587, and apply gyrochronology and activity calibrations, to determine its rotation period. Finally, based on an evolutionary model and the inferred rotation period, we use the EULAG-MHD code to perform global MHD simulations of HD43587 with the aim of getting some insights about its dynamo process. We confirm the almost flat activity profile, with a cycle period $P_{\mathrm{cyc}} = 10.44\pm3.03$ yrs deduced from the S-index time series, and a long-term trend that might be a period larger than 50 yrs. It was impossible to define a rotation period from the light-curve, however gyrochronology and activity calibrations allow to infer an indirect estimate of $\overline{P}_{\mathrm{rot}} = 22.6 \pm 1.9$ d. Furthermore, the MHD simulations confirm an oscillatory dynamo with a cycle period in good agreement with the observations and a low level of surface magnetic activity. We conclude that this object might be experiencing a "natural" decrease in magnetic activity as a consequence of its age. Nevertheless, the possibility that HD43587 is in a Maunder minimum phase cannot be ruled out.

The dynamo mechanism, responsible for the solar magnetic activity, is still an open problem in astrophysics. Different theories proposed to explain such phenomena have failed in reproducing the observational properties of the solar magnetism. Thus, ab-initio computational modeling of the convective dynamo in a spherical shell turns out as the best alternative to tackle this problem. In this work we review the efforts performed in global simulations over the past decades. Regarding the development and sustain of mean-flows, as well as mean magnetic field, we discuss the points of agreement and divergence between the different modeling strategies. Special attention is given to the implicit large-eddy simulations performed with the EULAG-MHD code.

Various models of solar subsurface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequence of increasing surface convection rates appears to be a significant alteration of the distribution of angular momentum, especially below the tachocline where the rotational frequency predominantly increases at higher latitudes. These hydrodynamic changes correspond to large shifts in the development of the current helicity in this stable layer ($r<0.72R_{\odot}$), significantly altering its impact on the generation of poloidal and toroidal fields at the tachocline and below, acting as a major contributor towards transitions in the dynamo cycle. The enhanced near-surface flow speed manifests in a global shift of the toroidal field ($B_{\phi}$) in the butterfly diagram - from a North-South symmetric pattern to a staggered anti-symmetric emergence.

Various models of solar sub-surface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequence of increasing surface convection rates appears to be a significant alteration of the distribution of angular momentum, especially below the tachocline where the rotational frequency predominantly increases at higher latitudes. These hydrodynamic changes correspond to large shifts in the development of the current helicity in this stable layer ($r<0.72R_{\odot}$), significantly altering its impact on the generation of poloidal and toroidal fields at the tachocline and below, acting as a major contributor towards transitions in the dynamo cycle. The enhanced near-surface flow speed manifests in a global shift of the toroidal field ($B_{\phi}$) in the butterfly diagram -- from a North-South symmetric pattern to a staggered anti-symmetric emergence.

Unveiling the evolution of toroidal field instability, known as Tayler instability, is essential to understand the strength and topology of the magnetic fields observed in early-type stars, in the core of the red giants, or in any stellar radiative zone. We want to study the non-linear evolution of the instability of a toroidal field stored in a stably stratified layer, in spherical symmetry and in the absence of rotation. In particular, we intend to quantify the suppression of the instability as a function of the Brunt-V\"ais\"ala ($\omega_{\rm BV}$) and the Alfv\'en ($\omega_{\rm A}$) frequencies. We use the MHD equations as implemented in the anelastic approximation in the EULAG-MHD code and perform a large series of numerical simulations of the instability exploring the parameter space for the $\omega_{\rm BV}$ and $\omega_{\rm A}$. We show that beyond a critical value gravity strongly suppress the instability, in agreement with the linear analysis. The intensity of the initial field also plays an important role: weaker fields show much slower growth rates. Moreover, in the case of very low gravity, the fastest growing modes have a large characteristic radial scale, at variance with the case of strong gravity, where the instability is characterized by horizontal displacements. Our results illustrate that the anelastic approximation can efficiently describe the evolution of toroidal field instability in stellar interiors. The suppression of the instability as a consequence of increasing values of $\omega_{\rm BV}$ might play a role to explain the magnetic desert in Ap/Bp stars since weak fields are only marginally unstable in the case of strong gravity.

Two fundamental properties of stellar magnetic fields have been determined by observations for solar-like stars with different Rossby numbers (Ro), namely, the magnetic field strength and the magnetic cycle period. The field strength exhibits two regimes: 1) for fast rotation it is independent of Ro, 2) for slow rotation it decays with Ro following a power law. For the magnetic cycle period two regimes of activity, the active and inactive branches, also have been identified. For both of them, the longer the rotation period, the longer the activity cycle. Using global dynamo simulations of solar like stars with Rossby numbers between ~0.4 and ~2, this paper explores the relevance of rotational shear layers in determining these observational properties. Our results, consistent with non-linear alpha^2-Omega dynamos, show that the total magnetic field strength is independent of the rotation period. Yet at surface levels, the origin of the magnetic field is determined by Ro. While for Ro<1 it is generated in the convection zone, for Ro>1 strong toroidal fields are generated at the tachocline and rapidly emerge towards the surface. In agreement with the observations, the magnetic cycle period increases with the rotational period. However, a bifurcation is observed for Ro~1, separating a regime where oscillatory dynamos operate mainly in the convection zone, from the regime where the tachocline has a predominant role. In the latter the cycles are believed to result from the periodic energy exchange between the dynamo and the magneto-shear instabilities developing in the tachocline and the radiative interior.

Recent observations of the magnetic field in pre-main sequence stars suggest that the magnetic field topology changes as a function of age. The presence of a tachocline could be an important factor in the development of magnetic field with higher multipolar modes. In this work we performed MHD simulations using the EULAG-MHD code to study the magnetic field generation and evolution in models that mimic stars at two evolutionary stages. The stratification for both stellar phases was computed by fitting stellar structure profiles obtained with the ATON stellar evolution code. The first stage is at $1.1 \mathrm{Myr}$, when the star is completely convective. The second stage is at $14 \mathrm{Myrs}$, when the star is partly convective, with a radiative core developed up to $30\%$ of the stellar radius. In this proceedings we present a preliminary analysis of the resulting mean-flows and magnetic field. The mean-flow analysis shown that the star rotate almost rigidly on the fully convective phase, whereas at the partially convective phase there is differential rotation with conical contours of iso-rotation. As for the mean magnetic field both simulations show similarities with respect to the field evolution. However, the topology of the magnetic field is different.

Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to the slow down of the axisymmetric poleward horizontal component. The behavior of this large scale flow deep inside the convection zone remains largely unknown. Present helioseismic techniques are not sensitive enough to capture the dynamics of this weak large scale flow. Moreover, the large time of integration needed to map the meridional circulation inside the convection zone, also masks some of the possible dynamics on shorter timescales. In this work we examine the dynamics of the meridional circulation that emerges from a 3D MHD global simulation of the solar convection zone. Our aim is to assess and quantify the behavior of meridional circulation deep inside the convection zone, where the cyclic large-scale magnetic field can reach considerable strength. Our analyses indicate that the meridional circulation morphology and amplitude are both highly influenced by the magnetic field, via the impact of magnetic torques on the global angular momentum distribution. A dynamic feature induced by these magnetic torques is the development of a prominent upward flow at mid latitudes in the lower convection zone that occurs near the equatorward edge of the toroidal bands and that peaks during cycle maximum. Globally, the dynamo-generated large-scale magnetic field drives variations in the meridional flow, in stark contrast to the conventional kinematic flux transport view of the magnetic field being advected passively by the flow.

We describe the preliminary design of a magnetograph and visible-light imager instrument to study the solar dynamo processes through observations of the solar surface magnetic field distribution. The instrument will provide measurements of the vector magnetic field and of the line-of-sight velocity in the solar photosphere. As the magnetic field anchored at the solar surface produces most of the structures and energetic events in the upper solar atmosphere and significantly influences the heliosphere, the development of this instrument plays an important role in reaching the scientific goals of The Atmospheric and Space Science Coordination (CEA) at the Brazilian National Institute for Space Research (INPE). In particular, the CEA's space weather program will benefit most from the development of this technology. We expect that this project will be the starting point to establish a strong research program on Solar Physics in Brazil. Our main aim is acquiring progressively the know-how to build state-of-art solar vector magnetograph and visible-light imagers for space-based platforms to contribute to the efforts of the solar-terrestrial physics community to address the main unanswered questions on how our nearby Star works.

The phenomenon of solar "torsional oscillations" (TO) represents migratory zonal flows associated with the solar cycle. These flows are observed on the solar surface and, according to helioseismology, extend through the convection zone. We study the origin of the TO using results from a global MHD simulation of the solar interior that reproduces several of the observed characteristics of the mean-flows and magnetic fields. Our results indicate that the magnetic tension (MT) in the tachocline region is a key factor for the periodic changes in the angular momentum transport that causes the TO. The torque induced by the MT at the base of the convection zone is positive at the poles and negative at the equator. A rising MT torque at higher latitudes causes the poles to speed-up, whereas a declining negative MT torque at the lower latitudes causes the equator to slow-down. These changes in the zonal flows propagate through the convection zone up to the surface. Additionally, our results suggest that it is the magnetic field at the tachocline that modulates the amplitude of the surface meridional flow rather than the opposite as assumed by flux-transport dynamo models of the solar cycle.

Rotational shear layers at the boundary between radiative and convective zones, tachoclines, play a key role in the process of magnetic field generation in solar-like stars. We present two sets of global simulations of rotating turbulent convection and dynamo. The first set considers a stellar convective envelope only; the second one, aiming at the formation of a tachocline, considers also the upper part of the radiative zone. Our results indicate that the resulting mean-flows and dynamo properties like the growth rate, saturation energy and mode depend on the Rossby (Ro) number. For the first set of models either oscillatory (with ~2 yr period) or steady dynamo solutions are obtained. The models in the second set naturally develop a tachocline which, in turn, leads to the generation of strong mean magnetic field. Since the field is also deposited into the stable deeper layer, its evolutionary time-scale is much longer than in the models without a tachocline. Surprisingly, the magnetic field in the upper turbulent convection zone evolves in the same time scale as the deep field. These models result in either an oscillatory dynamo with ~30 yr period or in a steady dynamo depending on Ro. In terms of the mean-field dynamo coefficients computed using FOSA, the field evolution in the oscillatory models without a tachocline seems to be consistent with dynamo waves propagating according to the Parker-Yoshimura sign rule. In the models with tachoclines the dynamics is more complex involving other transport mechanisms as well as tachocline instabilities.

To explore the physics of large-scale flows in solar-like stars, we perform 3D anelastic simulations of rotating convection for global models with stratification resembling the solar interior. The numerical method is based on an implicit large-eddy simulation approach designed to capture effects from non-resolved small scales. We obtain two regimes of differential rotation, with equatorial zonal flows accelerated either in the direction of rotation (solar-like) or in the opposite direction (anti-solar). While the models with the solar-like differential rotation tend to produce multiple cells of meridional circulation, the models with anti-solar differential rotation result in only one or two meridional cells. Our simulations indicate that the rotation and large-scale flow patterns critically depend on the ratio between buoyancy and Coriolis forces. By including a subadiabatic layer at the bottom of the domain, corresponding to the stratification of a radiative zone, we reproduce a layer of strong radial shear similar to the solar tachocline. Similarly, enhanced superadiabaticity at the top results in a near-surface shear layer located mainly at lower latitudes. The models reveal a latitudinal entropy gradient localized at the base of the convection zone and in the stable region, which however does not propagate across the convection zone. In consequence, baroclinicity effects remain small and the rotation iso-contours align in cylinders along the rotation axis. Our results confirm the alignment of large convective cells along the rotation axis in the deep convection zone, and suggest that such "banana-cell" pattern can be hidden beneath the supergranulation layer.

Convective turbulent motions in the solar interior, as well as the mean flows resulting from them, determine the evolution of the solar magnetic field. With the aim to get a better understanding of these flows we study anelastic rotating convection in a spherical shell whose stratification resembles that of the solar interior. This study is done through numerical simulations performed with the EULAG code. Due to the numerical formulation, these simulations are known as implicit large eddy simulations (ILES), since they intrinsically capture the contribution of, non-resolved, small scales at the same time maximizing the effective Reynolds number. We reproduce some previous results and find a transition between buoyancy and rotation dominated regimes which results in anti-solar or solar like rotation patterns. Even thought the rotation profiles are dominated by Taylor-Proudman columnar rotation, we are able to reproduce the tachocline and a low latitude near-surface shear layer. We find that simulations results depend on the grid resolution as a consequence of a different sub-grid scale contribution.

Many astrophysical bodies harbor magnetic fields that are thought to be sustained by a dynamo process. However, it has been argued that the production of large-scale magnetic fields by mean-field dynamo action is strongly suppressed at large magnetic Reynolds numbers owing to the conservation of magnetic helicity. This phenomenon is known as {\it catastrophic quenching}. Advection of magnetic fields by stellar and galactic winds toward the outer boundaries and away from the dynamo is expected to alleviate such quenching. Here we explore the relative roles played by advective and turbulent--diffusive fluxes of magnetic helicity in the dynamo. In particular, we study how the dynamo is affected by advection. We do this by performing direct numerical simulations of a turbulent dynamo of $\alpha^2$ type driven by forced turbulence in a Cartesian domain in the presence of a flow away from the equator where helicity changes sign. Our results indicate that in the presence of advection, the dynamo, otherwise stationary, becomes oscillatory. We confirm an earlier result for turbulent--diffusive magnetic helicity fluxes that for small magnetic Reynolds numbers ($\Rm\lesssim 100...200$, based on the wavenumber of the energy-carrying eddies) the magnetic helicity flux scales less strongly with magnetic Reynolds number ($\Rm^{-1/2}$) than the term describing magnetic helicity destruction by resistivity ($\Rm^{-1}$). Our new results now suggest that for larger $\Rm$ the former becomes approximately independent of $\Rm$, while the latter falls off more slowly. We show for the first time that both for weak and stronger winds, the magnetic helicity flux term becomes comparable to the resistive term for $\Rm\gtrsim 1000$, which is necessary for alleviating catastrophic quenching.

The role of MHD turbulence in astrophysical environments is still highly debated. An important question that permeates this debate is the transport of magnetic flux. This is particularly important, for instance, in the context of star formation. When clouds collapse gravitationally to form stars, there must be some magnetic flux transport. otherwise the new born stars would have magnetic fields several orders of magnitude larger than the observed ones. Also, the magnetic flux that is dragged in the late stages of the formation of a star can remove all the rotational support from the accretion disk that grows around the protostar. The efficiency of the mechanism which is often invoked to allow the transport of magnetic fields in the different stages of star formation, namely, the ambipolar diffusion, has been lately put in check. We here discuss an alternative mechanism for magnetic flux transport which is based on turbulent fast magnetic reconnection. We review recent results obtained from 3D MHD numerical simulations that indicate that this mechanism is very efficient for decoupling and transport magnetic flux from the inner denser regions to the outskirts of collapsing clouds in the different stages of star formation. We also discuss this mechanism in the context of dynamo processes and speculate that it can play a role both in the solar dynamo and in accretion disk dynamo processes.

Recent observations revealed that small magnetic elements abundant at the solar surface move poleward with a velocity which seems to be lower than the plasma velocity $U_{\theta}$. Guerrero et al. (2011) explained this discrepancy as a consequence of diffusive spreading of the magnetic elements due to a positive radial gradient of $|U_{\theta}|$. As the gradient's sign (inferred by local helioseismology) is still unclear, cases with a negative gradient are studied in this paper. Under this condition, the velocity of the magnetic tracers turns out to be larger than the plasma velocity, in disagreement with the observations. Alternative mechanisms for explaining them independently are proposed. For the turbulent magnetic pumping it is shown that it has to be unrealistically strong to reconcile the model with the observations.

Some selected concepts for the solar activity cycle are briefly reviewed. Cycle modulations through a stochastic alpha effect are being identified with limited scale separation ratios. Three-dimensional turbulence simulations with helicity and shear are compared at two different scale separation ratios. In both cases the level of fluctuations shows relatively little variation with the dynamo cycle. Prospects for a shallow origin of sunspots are discussed in terms of the negative effective magnetic pressure instability. Tilt angles of bipolar active regions are discussed as a consequence of shear rather than the Coriolis force.

We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Sun's surface and following their evolution for a prescribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensional surface-transport model the simulated MFT speed at the surface is always the same as the meridional flow-speed used as input to the model, but is different in a two-dimensional transport model in the meridional (r,theta) plane. The difference depends on the value of the magnetic diffusivity and on the radial gradient of the latitudinal velocity. We have confirmed our results with two different codes in spherical and Cartesian coordinates.

(abridged) Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Here we perform a systematic comparison between Cartesian and spherical geometries as a function of the rotation rate. Aims. We extend the earlier studies by using spherical wedges to obtain turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime. In particular, we want to clarify whether the sharp equatorial profile of the horizontal Reynolds stress found in earlier Cartesian models is reproduced in spherical models. Methods. We employ direct numerical simulations of turbulent convection. In order to reach as high spatial resolution as possible in the spherical runs, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong `banana cells'. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor--Proudman balance.

A hypothesis for sunspot formation is the buoyant emergence of magnetic flux tubes created by the strong radial shear at the tachocline. In this scenario, the magnetic field has to exceed a threshold value before it becomes buoyant and emerges through the whole convection zone. We follow the evolution of a random seed magnetic field with the aim of study under what conditions it is possible to excite the dynamo instability and whether the dynamo generated magnetic field becomes buoyantly unstable and emerges to the surface as expected in the flux-tube context. We perform numerical simulations of compressible turbulent convection that include a vertical shear layer. Like the solar tachocline, the shear is located at the interface between convective and stable layers. We find that shear and convection are able to amplify the initial magnetic field and form large-scale elongated magnetic structures. The magnetic field strength depends on several parameters such as the shear amplitude, the thickness and location of the shear layer, and the magnetic Reynolds number ($\Rm$). Whenever the toroidal magnetic field reaches amplitudes greater a threshold value which is close to the equipartition value, it becomes buoyant and rises into the convection zone where it expands and forms mushroom shape structures. Some events of emergence, i.e. those with the largest amplitudes of the initial field, are able to reach the very uppermost layers of the domain. These episodes are able to modify the convective pattern forming either broader convection cells or convective eddies elongated in the direction of the field. However, in none of these events the field preserves its initial structure.

The small-scale magnetic helicity produced as a by-product of the large-scale dynamo is believed to play a major role in dynamo saturation. In a mean-field model the generation of small-scale magnetic helicity can be modelled by using the dynamical quenching formalism. Catastrophic quenching refers to a decrease of the saturation field strength with increasing Reynolds number. It has been suggested that catastrophic quenching only affects the region of non-zero helical turbulence (i.e. where the kinematic alpha operates) and that it is possible to alleviate catastrophic quenching by separating the region of strong shear from the alpha layer. We perform a systematic study of a simple axisymmetric two-layer alpha-omega dynamo in a spherical shell for Reynolds numbers in the range 1 < Rm < 10^5. In the framework of dynamical quenching we show that this may not be the case, suggesting that magnetic helicity fluxes would be necessary.

Dynamos in the Sun and other bodies tend to produce magnetic fields that possess magnetic helicity of opposite sign at large and small scales, respectively. The build-up of magnetic helicity at small scales provides an important saturation mechanism. In order to understand the nature of the solar dynamo we need to understand the details of the saturation mechanism in spherical geometry. In particular, we want to understand the effects of magnetic helicity fluxes from turbulence and meridional circulation. We consider a model with just radial shear confined to a thin layer (tachocline) at the bottom of the convection zone. The kinetic alpha owing to helical turbulence is assumed to be localized in a region above the convection zone. The dynamical quenching formalism is used to describe the build-up of mean magnetic helicity in the model, which results in a magnetic alpha effect that feeds back on the kinetic alpha effect. In some cases we compare with results obtained using a simple algebraic alpha quenching formula. In agreement with earlier findings, the magnetic alpha effect in the dynamical alpha quenching formalism has the opposite sign compared with the kinetic alpha effect and leads to a catastrophic decrease of the saturation field strength with increasing magnetic Reynolds numbers. However, at high latitudes this quenching effect can lead to secondary dynamo waves that propagate poleward due to the opposite sign of alpha. Magnetic helicity fluxes both from turbulent mixing and from meridional circulation alleviate catastrophic quenching.

We present nonlinear mean-field alpha-Omega dynamo simulations in spherical geometry with simplified profiles of kinematic alpha effect and shear. We take magnetic helicity evolution into account by solving a dynamical equation for the magnetic alpha effect. This gives a consistent description of the quenching mechanism in mean-field dynamo models. The main goal of this work is to explore the effects of this quenching mechanism in solar-like geometry, and in particular to investigate the role of magnetic helicity fluxes, specifically diffusive and Vishniac-Cho (VC) fluxes, at large magnetic Reynolds numbers (Rm). For models with negative radial shear or positive latitudinal shear, the magnetic alpha effect has predominantly negative (positive) sign in the northern (southern) hemisphere. In the absence of fluxes, we find that the magnetic energy follows an Rm^-1 dependence, as found in previous works. This catastrophic quenching is alleviated in models with diffusive magnetic helicity fluxes resulting in magnetic fields comparable to the equipartition value even for Rm=10^7. On the other hand, models with a shear-driven Vishniac-Cho flux show an increase of the amplitude of the magnetic field with respect to models without fluxes, but only for Rm<10^4. This is mainly a consequence of assuming a vacuum outside the Sun which cannot support a significant VC flux across the boundary. However, in contrast with the diffusive flux, the VC flux modifies the distribution of the magnetic field. In addition, if an ill-determined scaling factor in the expression for the VC flux is large enough, subcritical dynamo action is possible that is driven by the action of shear and the divergence of current helicity flux.

In the non-linear phase of a dynamo process, the back-reaction of the magnetic field upon the turbulent motion results in a decrease of the turbulence level and therefore in a suppression of both the magnetic field amplification (the alpha-quenching effect) and the turbulent magnetic diffusivity (the eta-quenching effect). While the former has been widely explored, the effects of eta-quenching in the magnetic field evolution have rarely been considered. In this work we investigate the role of the suppression of diffusivity in a flux-transport solar dynamo model that also includes a non-linear alpha-quenching term. Our results indicate that, although for alpha-quenching the dependence of the magnetic field amplification with the quenching factor is nearly linear, the magnetic field response to eta-quenching is non-linear and spatially non-uniform. We have found that the magnetic field can be locally amplified in this case, forming long-lived structures whose maximum amplitude can be up to ~2.5 times larger at the tachocline and up to ~2 times larger at the center of the convection zone than in models without quenching. However, this amplification leads to unobservable effects and to a worse distribution of the magnetic field in the butterfly diagram. Since the dynamo cycle period increases when the efficiency of the quenching increases, we have also explored whether the eta-quenching can cause a diffusion-dominated model to drift into an advection-dominated regime. We have found that models undergoing a large suppression of the diffusivity produce a strong segregation of magnetic fields that may lead to unsteady dynamo-oscillations. On the other hand, an initially diffusion-dominated model undergoing a small suppression in eta remains in the diffusion-dominated regime.

The sunspot solar cycle has been usually explained as the result of a dynamo process operating in the sun. This is a classical problem in Astrophysics that until the present is not fully solved. Here we discuss current problems and limitations with the solar dynamo modeling and their possible solutions using a kinematic dynamo model within the Babcock-Leighton approximation as a tool. In particular, we discuss the importance of the turbulent magnetic pumping versus the meridional flow circulation in the dynamo operation.

The turbulent pumping effect corresponds to the transport of magnetic flux due to the presence of density and turbulence gradients in convectively unstable layers. In the induction equation it appears as an advective term and for this reason it is expected to be important in the solar and stellar dynamo processes. In this work, we have explored the effects of the turbulent pumping in a flux-dominated Babcock-Leighton solar dynamo model with a solar-like rotation law. The results reveal the importance of the pumping mechanism for solving current limitations in mean field dynamo modeling such as the storage of the magnetic flux and the latitudinal distribution of the sunspots. In the case that a meridional flow is assumed to be present only in the upper part of the convective zone, it is the full turbulent pumping that regulates both the period of the solar cycle and the latitudinal distribution of the sunspots activity.

Flux-dominated solar dynamo models have demonstrated to reproduce the main features of the large scale solar magnetic cycle, however the use of a solar like differential rotation profile implies in the the formation of strong toroidal magnetic fields at high latitudes where they are not observed. In this work, we invoke the hypothesis of a thin-width tachocline in order to confine the high-latitude toroidal magnetic fields to a small area below the overshoot layer, thus avoiding its influence on a Babcock-Leighton type dynamo process. Our results favor a dynamo operating inside the convection zone with a tachocline that essentially works as a storage region when it coincides with the overshoot layer.

Flux-dominated solar dynamo models, which have demonstrated to be quite successful in reproducing most of the observed features of the large scale solar magnetic cycle, generally produce an inappropriate latitudinal distribution of the toroidal magnetic fields, showing fields of large magnitude in polar regions where the radial shear has a maximum amplitude. Employing a kinematic solar dynamo model, we here explore the contribution of both the radial and the latitudinal shear in the generation of the toroidal magnetic fields by varying the shape and the thickness of the solar tachocline. We also explore the effects of the diffusivity profile of the convective zone. We find that the latitudinal component is always dominant over the radial component at producing toroidal field amplification. These results are very sensitive to the adopted diffusivity profile, specially in the inner convection zone. A diagram of the toroidal field at a latitude of 60 degrees versus the diffusivity at the convection layer for different values of the tachocline width has revealed that these fields are mainly eliminated for tachoclines with width d_1>0.08 Ro, for a restric range of diffusivites; or for d_1<0.02 and almost any value of eta_c in the appropriate solar range. For intermediate values of d_1=0.04Ro-0.06Ro, strong toroidal fields should survive at high latitudes in the butterfly diagram and those values are therefore not suitable.

Circulation-dominated solar dynamo models, which employ a helioseismic rotation profile and a fixed meridional flow, give a good approximation to the large scale solar magnetic phenomena, such as the 11-year cycle or the so called Hale's law of polarities. Nevertheless, the larger amplitude of the radial shear at the high latitudes makes the dynamo to produce a strong toroidal magnetic field at high latitudes, in contradiction with the observations of the sunspots (Sporer`s Law). A possible solution was proposed by Nandy and Choudhur (2002) in which a deep meridional flow can conduct the magnetic field inside of a stable layer (the radiative core) and then allow that it erupts just at lower latitudes. We have recently explored this hypothesis in a different approximation, using the magnetic buoyancy mechanism proposed by Dikpaty and Charbonneau (1999) and found that a deep meridional flow pushes the maximum of the toroidal magnetic field towards the solar equator, but, in contrast to Nandy and Choudhuri (2002), a second zone of maximum fields remains at the poles (Guerrero and Munoz (2004)). We here review these results and also discuss a new possible scenario where the tachocline has an ellipsoidal shape, following early helioseismologic observations Charbonneau et. al., (1999), and find that the modification of the geometry of the tachocline can lead to results which are in good agreement with observations and opens the possibility to explore in more detail, through the dynamo model, the place where the magnetic field could be really stored.

We develop two different solar dynamo models to verify the hypothesis that a deep meridional flow can restrict the apperance of sunspots below 45 degrees, proposed by Nandy & Choudhuri (2002). In the first one, a single polytropic approximation for the density profile was taken, for both radiative and convective zones. In the second one, two polytropes were used to distinguish between both zones Pinzon & Calvo-Mozo (2001). The magnetic buoyancy mechanism proposed by Dikpati & Charbonneau (1999) was chosen in both models. We, actually, have obtained that a deep meridional flow pushes the maxima of toroidal magnetic field toward the solar equator, but in contrast to Nandy & Choudhuri (2002) a second zone of maximal fields remains at the poles. The second model, although closely resembling the solar standard model of Bahcall, Pinsonneault & Wasserbug (1995); Bahcall, Pinsonneault & Basu (2001), gives solar cyles three times longer than observed.

We present the results of optical observations of the GRB990510 field carried out at different epochs from European Southern Observatory (ESO) telescopes. Deep observations, down to a limiting magnitude of about 27 and 24 in the Bessel-R and Gunn-I band, respectively, were obtained between May 16 and 18 from the ESO NTT-SUSI2 telescope and on May 20 from the ESO 3.6 m (EFOSC2) telescope. These observations, together with other published photometric data, allowed to monitor the faint tail of the decaying Optical Transient (OT) associated to the GRB990510. We discuss the light curves in the different filters (V, R and I) in the light of the recently proposed decay laws. No obvious host associated to the GRB990510 optical afterglow was found in the R and I band images. By comparing the light curves with respect to the theoretical colors of galaxies with different morphology we derived a lower limit of R about 26.6 for the host galaxy.