Rotation periods inferred from the magnetic variability of some Ap stars are incredibly long, exceeding ten years in some cases. An explanation for such slow rotation is lacking. This paper attempts to provide an explanation of the super-slow rotation of the magnetic and thermal patterns of Ap stars in terms of the longitudinal drift of the unstable disturbances of the kink-type (Tayler) instability of their internal magnetic field. The rates of drift and growth were computed for eigenmodes of Tayler instability using stellar parameters estimated from a structure model of an A star. The computations refer to the toroidal background magnetic field of varied strength. The non-axisymmetric unstable disturbances drift in a counter-rotational direction in the co-rotating reference frame. The drift rate increases with the strength of the background field. For a field strength exceeding the (equipartition) value of equal Alfven and rotational velocities, the drift rate approaches the proper rotation rate of a star. The eigenmodes in an inertial frame show very slow rotation in this case. Patterns of magnetic and thermal disturbances of the slowly rotating eigenmodes are also computed. The counter-rotational drift of Tayler instability eigenmodes is a possible explanation for the observed phenomenon of super-slowly rotating Ap stars.

The rise of flux-tubes with intense magnetic fields from the base of the convection zone to the solar surface has been substantiated as a probable mechanism for sunspot formation. The origin of flux-tubes of sufficient strength (~10^5 G) is however uncertain. This paper considers the instability of a large-scale toroidal magnetic field caused by the magnetic suppression of convective heat transport as a candidate for the flux tube forming mechanism. The consideration employs the analytical dependence of the eddy diffusion on the magnetic field supplied by mean-field magnetohydrodynamics. The instability tends to produce regions of increased field strength with spatial scales of an order of 100 Mm at the base of the convection zone. Characteristic growth times of the instability are short compared to the 11-year cycle. The threshold field strength for the onset of the instability increases from several hundred Gauss in the vicinity of the equator to some kilo-Gauss at middle latitudes. Growth rates of unstable disturbances decrease with latitude. These latitudinal trends can be the reason for the observed confinement of sunspot activity to a near-equatorial belt.

Asteroseismology has revealed that cores of red giants rotate about one order of magnitude faster than their convective envelopes. This paper attempts an explanation for this rotational state in terms of the theory of angular momentum transport in stellar convection zones. A differential rotation model based on the theory is applied to a sequence of evolutionary states of a red giant of one solar mass. The model computations show a rotation of about ten times faster in the cores compared to the stellar surface. This rotational state is caused by the non-diffusive downward convective transport of angular momentum. The contrast in rotational rates between core and envelope increases with the radius (age) of the star. Seismologically detected scaling for the spindown of the giants' cores is also reproduced.

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.

Conservation of magnetic helicity by the Hall drift does not prevent Hall instability of helical fields. This conclusion follows from stability analysis of a force-free spatially-periodic Hall equilibrium. The growth rates of the instability scale as $\sigma \propto B^{3/4}\eta^{1/4}$ with the field strength $B$ and magnetic diffusivity $\eta$ and can be large compared to the rate of resistive decay of the background field. The instability deviates the magnetic field from the force-free configuration. The unstable eigenmodes include a fine spatial structure which evolves into current sheets at the nonlinear stage of the instability. The instability catalyses the resistive release of magnetic energy. The energy is released in a sequence of spikes, every spike emits several percent of the total energy. A numerically defined scaling for the energy released in a single spike permits an extrapolation to astrophysically relevant values of the Hall number. The instability can be relevant to magnetic energy release in a neutron star crust and, possibly, in stellar coronae.

The hemispheric asymmetry of sunspot activity observed possesses a regular component varying on a time scale of several solar cycles whose origin and properties are currently debated. This paper addresses the question of whether the long-term hemispheric asymmetry can result from random variations of solar dynamo parameters in time and latitude. Scatter in the observed tilt angles of sunspot groups is estimated to infer constraints on fluctuations in the dynamo mechanism for poloidal field regeneration. A dynamo model with fluctuations in the Babcock-Leighton type $\alpha$-effect is designed in accordance with these constraints and then used to compute a large number of magnetic cycles for statistical analyses of their hemispheric asymmetry. Hemispheric asymmetry in the simulated dynamo results from the presence of an equator-symmetric part in the oscillating magnetic field. The subdominant quadrupolar oscillations are stochastically forced by dominant dipolar oscillations via the equator-symmetric part of the fluctuating $\alpha$-effect. The amplitude and sense of the asymmetry of individual cycles varies on a time scale of the order of four dynamo-cycle periods. The variations are irregular, i.e. not periodic. The model suggests that asymmetry in the polar magnetic fields in the solar minima can be used as a precursor for asymmetry of sunspot activity in the following solar cycle.

Synoptic magnetograms and relevant proxy data were analyzed to study the evolution of the Sun's polar magnetic fields. Time-latitude analysis of large-scale magnetic fields demonstrates cyclic changes in their zonal structure and the polar-field buildup. The time-latitude distributions of the emergent and remnant magnetic flux enable us to examine individual features of recent cycles. The poleward transport of predominantly following polarities contributed much of the polar flux and led to polar-field reversals. Multiple reversals of dominant polarities at the Sun's poles were identified in Cycles 20 and 21. Three-fold reversals were caused by remnant flux surges of following and leading polarities. Time-latitude analysis of solar magnetic fields in Cycles 20--24 revealed zones which are characterized by a predominance of negative (non-Joy's) tilts and appearance of active regions which violate Hale's polarity law. The decay of non-Joy's and anti-Hale's active regions result in remnant flux surges which disturb the usual order in magnetic flux transport and sometimes lead to multiple reversals of polar fields. The analysis of local and large-scale magnetic fields in their causal relation improved our understanding of the Sun's polar field weakening.

Magnetic quenching of turbulent thermal diffusivity leads to instability of the large-scale field with the production of spatially isolated regions of enhanced field. This conclusion follows from a linear stability analysis in the framework of mean-field magnetohydrodynamics that allows for thermal diffusivity dependence on the magnetic field. The characteristic growth time of the instability is short compared to the 11-year period of solar activity. The characteristic scale of the increased field regions measures in tens of mega-meters. The instability can produce magnetic inhomogeneities whose buoyant rise to the solar surface forms the solar active regions. The magnetic energy of the field concentrations coincides in order of magnitude with the energy of the active regions.

The duration of activity growths in solar cycles is on average shorter than the duration of its declines. This asymmetry can result from fluctuations in dynamo parameters. A solar dynamo model with fluctuations in the $\alpha$-effect shows the statistical asymmetry which increases with both fluctuation amplitude and coherence time. An interpretation for the asymmetry origin is suggested, which predicts a correlation between the asymmetry measure and delay of the polar field reversals relative to the activity maxima. Data on the twelve latest solar cycles confirm such a correlation.

The Kepler mission identified huge flares on various stars including some of solar type. These events are substantially more energetic than solar flares, and so they are referred to as superflares. Even a small probability of such a superflare on the Sun would be a menace to modern society. A flare comparable in energy with that of superflares was observed on 24th and 25th September on the binary HK Lac. Unlike the Kepler stars, there are observations of differential rotation for HK Lac. This differential rotation appears to be anti-solar. For anti-solar differential rotation, dynamo models can give magnetic activity waves of dipole symmetry as well as quasi-stationary magnetic configurations with quadrupole symmetry. The magnetic energy of such stationary configurations is usually about two orders of magnitude higher than that associated with activity waves. We believe that this mechanism could provide sufficient energy to produce superflares on late type stars, and present some simple models in support of this idea.

Context. Solar activity cycles vary in amplitude and duration. The variations can be at least partly explained by fluctuations in dynamo parameters. Aims. We want to restrict uncertainty in fluctuating dynamo parameters and find out which properties of the fluctuations control the amplitudes of the magnetic field and energy in variable dynamo cycles. Methods. A flux-transport model for the solar dynamo with fluctuations of the Babcock-Leighton type $\alpha$-effect was applied to generate statistics of magnetic cycles for our purposes. The statistics were compared with data on solar cycle periods to restrict the correlation time of dynamo fluctuations. Results. A characteristic time of fluctuations in the $\alpha$-effect is estimated to be close to the solar rotation period. The fluctuations produce asymmetry between the times of rise and descent of dynamo cycles, the rise time being on average shorter. The affect of the fluctuations on cycle amplitudes depends on the phase of the cycle in which the fluctuations occur. Negative fluctuations (decrease in $\alpha$) in the rise phase delay decay of poloidal field and increase the cycle amplitude in toroidal field and magnetic energy. Negative fluctuation in the decline phase reduces the polar field at the end of a cycle and the amplitude of the next cycle. The low amplitude of the 24th solar cycle compared to the preceding 23rd cycle can be explained by this effect. Positive fluctuations in the descent phase enhance the magnetic energy of the next cycle by increasing the seed poloidal field for the next cycle. The statistics of the computed energies of the cycles suggest that superflares of $\ge 10^{34}$ erg are not possible on the Sun.

Asteroseismological determinations of stellar ages have shown that old main-sequence dwarfs do not obey gyrochronology. Their rotation is slow compared to young stars but faster than gyrochronology predicts. This can be explained by the presence of a maximum rotation period beyond which the large-scale dynamo switches off and stops providing global magnetic fields necessary for stellar spindown. Assuming this explanation, the excess of stellar dynamo parameters over their marginal values can be estimated for given spectral type and rotation rate. The estimation gives the dynamo number for the Sun about 10% above its critical value. The corresponding dynamo model provides - though with some further tuning - reasonable results for the Sun. Following the same approach, the differential rotation and marginal dynamo modes are computed for stars between 0.7 and 1.2 solar masses. With increasing stellar mass, the differential rotation and the ratio of toroidal-to-poloidal field are predicted to increase while the field topology changes from dipolar to mixed quadrupolar-dipolar parity.

A pictorial explanation for shear-Hall instability is suggested and shows that the shear flow is not necessary for the instability because its role can be played by the Hall effect of an inhomogeneous background magnetic field. Linear stability analysis for a simple model of magnetic field varying periodically in space confirms such a 'double Hall' instability. Numerical computations show a considerable increase in Ohmic dissipation rate at the nonlinear stage of instability development. Field dissipation has a spiky character associated with magnetic reconnection in current sheets and X-points. Double Hall instability can be significant for magnetic field dissipation in neutron star crusts and, possibly, in the solar corona.

Recent data from the Kepler mission has revealed the occurrence of superflares in sun-like stars which exceed by far any observed solar flares in release of energy. A natural idea is that the dynamo mechanism in superflaring stars differs in some respect from that in the Sun. We search for a difference in the dynamo-related parameters between superflaring stars and the Sun to suggest a dynamo-mechanism as close as possible to the conventional solar/stellar dynamo but capable of providing much higher magnetic energy. Dynamo based on joint action of differential rotation and mirror asymmetric motions can in principle result in excitation of two types of magnetic fields. First of all, it is well-known in solar physics dynamo waves. The point is that another magnetic configuration with initial growth and further stabilisation is also possible for excitation. For comparable conditions, magnetic field strength of second configuration is much larger rather of the first one just because dynamo do not spend its efforts for periodic magnetic field inversions but use its for magnetic field growth. We analysed available data from the Kepler mission concerning the superflaring stars in order to find tracers of anomalous magnetic activity. Starting from the recent paper, we find that anti-solar differential rotation or anti-solar sign of the mirror-asymmetry of stellar convection can provide the desired strong magnetic field in dynamo models. We confirm this concept by numerical models of stellar dynamos with corresponding governing parameters. We conclude that the proposed mechanism can plausibly explain the superflaring events at least for some cool stars, including binaries, subgiants and, possibly, low-mass stars and young rapid rotators.

A model for the solar dynamo, consistent in global flow and numerical method employed with the differential rotation model, is developed. The magnetic turbulent diffusivity is expressed in terms of the entropy gradient, which is controlled by the model equations. The magnetic Prandtl number and latitudinal profile of the alpha-effect are specified by fitting the computed period of the activity cycle and the equatorial symmetry of magnetic fields to observations. Then, the instants of polar field reversals and time-latitude diagrams of the fields also come into agreement with observations. The poloidal field has a maximum amplitude of about 10 Gs in the polar regions. The toroidal field of several thousand Gauss concentrates near the base of the convection zone and is transported towards the equator by the meridional flow. The model predicts a value of about 10^37 erg for the total magnetic energy of large-scale fields in the solar convection zone.

Mean-field hydrodynamics advanced to clear explanations for the origin and properties of the global meridional flow in stellar convection zones. Qualitative arguments and analysis of basic equations both show that the meridional circulation is driven by non-conservative centrifugal and buoyancy forces and results from a slight disbalance between these two drivers. The deviations from the thermal wind balance are relatively large near the boundaries of convection zones. Accordingly, the meridional flow attains its largest velocities in the boundary layers and decreases inside the convection zone. This picture, however, is neither supported nor dismissed by the conflicting results of recent helioseismic soundings or 3D numerical experiments. The relevant physics of the differential temperature and its possible relation to the solar oblateness are briefly discussed.

Recent data on superflares on sun-like stars and radiocarbon data on solar activity in the past are both indicative of transient epochs of unusually high magnetic activity. We propose an explanation for the grand activity maxima in the framework of a solar dynamo model with fluctuating parameters. Solar-type dynamos are oscillatory because of the combination of the solar-type differential rotation with positive (in the Northern hemisphere) alpha-effect. An artificial reversal of the sign in the alpha-effect changes the dynamo to a steady regime with hundreds of times larger magnetic energy compared to the amplitude of the cyclic dynamo. Sufficiently large and durable fluctuations reversing the sign of the alpha-effect during the growth phase of a magnetic cycle can, therefore, cause a transient change to a steady dynamo with considerably increased magnetic energy. This qualitative scenario for grand activity maxima is supported by computations of the dynamo model with a fluctuating alpha-effect. The computed statistics of several thousand magnetic cycles gives examples of cycles with very high magnetic energy. Our preliminary estimations, however, suggest that the probability of solar superflares is extremely low.

Solar dynamo models require some mechanism for magnetic field concentration near the base of the convection zone in order to generate super-kilogauss toroidal fields with sufficiently large (~10^{24} Mx) magnetic flux. We consider the downward diamagnetic pumping near the base of the convection zone as a possible concentration mechanism and derive the pumping velocities with allowance for the effect of rotation. Transport velocities for poloidal and toroidal fields differ in rotating fluid. The toroidal field is transported downward along the radius only but the pumping velocity for the poloidal field has an equatorward meridional component also. Previous results for cases of slow and rapid rotation are reproduced and the diamagnetic pumping expressions adapted for use in dynamo models are presented.

Helioseismology revealed an increase in the rotation rate with depth just beneath the solar surface. The relative magnitude of the radial shear is almost constant with latitude. This rotational state can be interpreted as a consequence of two conditions characteristic of the near-surface convection: the smallness of convective turnover time in comparison with the rotation period and absence of a horizontal preferred direction of convection anisotropy. The latter condition is violated in the presence of a magnetic field. This raises the question of whether the subphotospheric fields can be probed with measurements of near-surface rotational shear. The shear is shown to be weakly sensitive to magnetic fields but can serve as a probe for sufficiently strong fields of the order of one kilogauss. It is suggested that the radial differential rotation in extended convective envelopes of red giants is of the same origin as the near-surface rotational shear of the Sun.

The magnetic activity of solar-type stars generally increases with stellar rotation rate. The increase, however, saturates for fast rotation. The Babcock-Leighton mechanism of stellar dynamos saturates as well when the mean tilt-angle of active regions approaches ninety degrees. Saturation of magnetic activity may be a consequence of this property of the Babcock-Leighton mechanism. Stellar dynamo models with a tilt-angle proportional to the rotation rate are constructed to probe this idea. Two versions of the model - treating the tilt-angles globally and using Joy's law for its latitude dependence - are considered. Both models show a saturation of dynamo-generated magnetic flux at high rotation rates. The model with latitude-dependent tilt-angles also shows a change in dynamo regime in the saturation region. The new regime combines a cyclic dynamo at low latitudes with an (almost) steady polar dynamo.

Long-term variations of solar activity, including the Grand minima, are believed to result from temporal variations of dynamo parameters. The simplest approximation of dynamo waves is applied to show that cyclic variations of the parameters can lead to an exponential growth or decay of magnetic oscillations depending on the variations frequency. There is no parametric resonance in a dynamo, however: the selective sensitivity to distinct frequencies, characteristic of resonant phenomena, is absent. A qualitative explanation for this finding is suggested. Nonlinear analysis of dynamo-waves reveals the hysteresis phenomenon found earlier in more advanced models. However, the simplified model allows a computation of a sufficiently large number of dynamo-cycles for constructing the distribution function of their amplitudes to reproduce qualitatively two modes of solar activity inferred recently from cosmogenic isotope content in natural archives.

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

With application to inner stellar radiative zones, a linear theory is used to analyze the instability of a dipole-parity toroidal background field, in the presence of density stratification, differential rotation, and realistically small Prandtl numbers. The physical parameters are the normalized latitudinal shear $a$ and the normalized field amplitude $b$. Only the solutions for the wavelengths with the maximal growth rates are considered. If these scales are combined to the radial values of velocity, one finds that the (very small) radial velocity only depends slightly on $a$ and $b$, so that it can be used as the free parameter of the eigenvalue system. The resulting instability-generated tensors of magnetic diffusivity and eddy viscosity are highly anisotropic. The eddy diffusivity in latitudinal direction exceeds the eddy diffusivity in radial direction by orders of magnitude. Its latitudinal profile shows a strong concentration toward the poles which is also true for the effective viscosity which has been calculated via the angular momentum transport of the instability pattern. The resulting effective magnetic Prandtl number reaches values of $O(10^2)$, so that the differential rotation decays much faster than the toroidal background field, which is {the} necessary condition to explain the observed slow rotation of the early red-giant and sub-giant cores by means of magnetic instabilities.

Recent progress in observational studies of magnetic activity in M dwarfs urgently requires support from ideas of stellar dynamo theory. We propose a strategy to connect observational and theoretical studies. In particular, we suggest four magnetic configurations that appear relevant to dwarfs from the viewpoint of the most conservative version of dynamo theory, and discuss observational tests to identify the configurations observationally. As expected, any such identification contains substantial uncertainties. However the situation in general looks less pessimistic than might be expected. Several identifications between the phenomenology of individual stars and dynamo models are suggested. Remarkably, all models discussed predict substantial surface magnetic activity at rather high stellar latitudes. This prediction looks unexpected from the viewpoint of our experience observing the Sun (which of course differs in some fundamental ways from these late-type dwarfs). We stress that a fuller understanding of the topic requires a long-term (at least 15 years) monitoring of M dwarfs by Zeeman-Doppler imaging.

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.

An explanation is suggested for the north-south asymmetry of the polar magnetic field reversal in the current cycle of solar activity. The contribution of the Babcock-Leighton mechanism to the poloidal field generation is estimated using sunspot data for the current activity cycle. Estimations are performed separately for the northern and southern hemispheres. The contribution of the northern hemisphere exceeded considerably that of the southern hemisphere during the initial stage of the cycle. This is the probable reason for the earlier reversal of the northern polar field. The estimated contributions of the Babcock-Leighton mechanism are considerably smaller than similar estimations for the previous activity cycles. A relatively weak (<1G) large-scale polar field can be expected for the next activity minimum.

It can be shown on observational grounds that two basic effects of dynamo theory for solar activity - production of the toroidal field from the poloidal one by differential rotation and reverse conversion of the toroidal field to the poloidal configuration by helical motions - are operating in the Sun. These two effects, however, do not suffice for constructing a realistic model for the solar dynamo. Only when a non-local version of the alpha-effect is applied, is downward diamagnetic pumping included and field advection by the equatorward meridional flow near the base of the convection zone allowed for, can the observed activity cycles be closely reproduced. Fluctuations in the alpha-effect can be estimated from sunspot data. Dynamo models with fluctuating parameters reproduce irregularities of solar cycles including the grand activity minima. The physics of parametric excitation of irregularities remains, however, to be understood.

Surfaces of constant pressure and constant density do not coincide in differentially rotating stars. Stellar radiation zones with baroclinic stratification can be unstable. Instabilities in radiation zones are of crucial importance for angular momentum transport, mixing of chemical species and, possibly, for magnetic field generation. This paper performs linear analysis of baroclinic instability in differentially rotating stars. Linear stability equations are formulated for differential rotation of arbitrary shape and then solved numerically for rotation non-uniform in radius. As the differential rotation increases, r- and g-modes of initially stable global oscillations transform smoothly into growing modes of baroclinic instability. The instability can therefore be interpreted as stability loss to r- and g-modes excitation. Regions of stellar parameters where r- or g-modes are preferentially excited are defined. Baroclinic instability onsets at a very small differential rotation of below 1%. The characteristic time of instability growth is about one thousand rotation periods. Growing disturbances possess kinetic helicity. Magnetic field generation by the turbulence resulting from baroclinic instability in differentially rotating radiation zones is, therefore, possible.

M-dwarfs demonstrate two types of activity: 1) strong (kilogauss) almost axisymmetric poloidal magnetic fields; and 2) considerably weaker nonaxisymmetric fields, sometimes including a substantial toroidal component. Dynamo bistability has been proposed as an explanation. However it is not straightforward to obtain such a bistability in dynamo models. On the other hand, the solar magnetic dipole at times of magnetic field inversion becomes transverse to the rotation axis, while the magnetic field becomes weaker at times far from that of inversion. Thus the Sun resembles a star with the second type of activity. We suggest that M-dwarfs can have magnetic cycles, and that M-dwarfs with the second type of activity can just be stars observed at times of magnetic field inversion. Then the relative number of M-dwarfs with the second type of activity can be used in the framework of this model to determine parameters of stellar convection near the surface.

A solar-type dynamo model in a spherical shell is developed with allowance for random dependence of the poloidal field generation mechanism on time and latitude. The model shows repeatable epochs of strongly decreased amplitude of magnetic cycles similar to the Maunder minimum of solar activity. Random dependence of dynamo-parameters on latitude brakes the equatorial symmetry of generated fields. The model shows the correlation of the occurrence of grand minima with deviations in the dynamo-field from dipolar parity. An increased north-south asymmetry of magnetic activity can, therefore, be an indicator of transitions to grand minima. Qualitative interpretation of this correlation is suggested. Statistics of grand minima in the model are close to the Poisson random process indicating that the onset of a grand minimum is statistically independent of preceding minima.

A linear analysis of baroclinic instability in a stellar radiation zone with radial differential rotation is performed. The instability onsets at a very small rotation inhomogeneity. There are two families of unstable disturbances corresponding to Rossby waves and internal gravity waves. The instability is dynamical: its growth time of several thousand rotation periods is short compared to the stellar evolution time. A decrease in thermal conductivity amplifies the instability. Unstable disturbances possess kinetic helicity thus indicating the possibility of magnetic field generation by the turbulence resulting from the instability.

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.

During a solar magnetic field reversal the magnetic dipole moment does not vanish, but migrates between poles, in contradiction to the predictions of mean-field dynamo theory. We try to explain this as a consequence of magnetic fluctuations. We exploit the statistics of fluctuations to estimate observable signatures. Simple statistical estimates, taken with results from mean-field dynamo theory, suggest that a non-zero dipole moment may persist through a global field reversal. Fluctuations in the solar magnetic field may play a key role in explaining reversals of the dipolar component of the field.

Meridional flow results from slight deviations from the thermal wind balance. The deviations are relatively large in the boundary layers near the top and bottom of the convection zone. Accordingly, the meridional flow attains its largest velocities at the boundaries and decreases inside the convection zone. The thickness of the boundary layers, where meridional flow is concentrated, decreases with rotation rate, so that an advection-dominated regime of dynamos is not probable in rapidly rotating stars. Angular momentum transport by convection and by the meridional flow produce differential rotation. The convective fluxes of angular momentum point radially inward in the case of slow rotation but change their direction to equatorward and parallel to the rotation axis as the rotation rate increases. The differential rotation of main-sequence dwarfs is predicted to vary mildly with rotation rate but increase strongly with stellar surface temperature. The significance of differential rotation for dynamos has the opposite tendency to increase with spectral type.

Gyrochronology and recent theoretical findings are used to reduce the number of input parameters of differential rotation models. This eventually leads to a theoretical prediction for the surface differential rotation as a function of only two stellar parameters - surface temperature and rotation period - that can be defined observationally. An analytical approximation for this function is suggested. The tendency for the differential rotation to increase with temperature is confirmed. The increase is much steeper for late F-stars compared to G- and K-dwarfs. Slow and fast rotation regimes for internal stellar rotation are identified. A star attains its maximum differential rotation at rotation rates intermediate between these two regimes. The amplitude of the meridional flow increases with surface temperature and rotation rate. The structure of the flow changes considerably between cases of slow and fast rotation. The flow in rapid rotators is concentrated in the boundary layers near the top and bottom of the convection zone with very weak circulation in between.

Small levels of turbulence can be present in stellar radiative interiors due to, e.g., instability of rotational shear. In this paper we estimate turbulent transport coefficients for stably stratified rotating stellar radiation zones. Stable stratification induces strong anisotropy with a very small ratio of radial-to-horizontal turbulence intensities. Angular momentum is transported mainly due to the correlation between azimuthal and radial turbulent motions induced by the Coriolis force. This non-diffusive transport known as the Lambda-effect has outward direction in radius and is much more efficient compared to the effect of radial eddy viscosity. Chemical species are transported by small radial diffusion only. This result is confirmed using direct numerical simulations combined with the test-scalar method. As a consequence of the non-diffusive transport of angular momentum, the estimated characteristic time of rotational coupling (< 100 Myr) between radiative core and convective envelope in young solar-type stars is much shorter compared to the time-scale of Lithium depletion (~1 Gyr).

The components of the total stress tensor (Reynolds stress plus Maxwell stress) are computed within the quasilinear approximation for a driven turbulence influenced by a large-scale magnetic background field. The conducting fluid has an arbitrary magnetic Prandtl number and the turbulence without the background field is assumed as homogeneous and isotropic with a free Strouhal number St. The total large-scale magnetic tension is always reduced by the turbulence with the possibility of a `catastrophic quenching' for large magnetic Reynolds number Rm so that even its sign is reversed. The total magnetic pressure is enhanced by turbulence with short correlation time (`white noise') but it is reduced by turbulence with long correlation time. Also in this case the sign of the total pressure may reverse but only for special turbulences with sufficiently large St> 1. The turbulence-induced terms of the stress tensor are suppressed by strong magnetic fields. For the tension term this quenching grows with the square of the Hartmann number of the magnetic field. For microscopic (i.e. small) diffusivity values the magnetic tension term becomes thus highly quenched even for field amplitudes much smaller than their equipartition value. In the opposite case of large-eddy simulations the magnetic quenching is only mild but then also the turbulence-induced Maxwell tensor components for weak fields remain rather small.

Helioseismology shows that the regions occupied by convection and differential rotation inside the sun almost coincide. This supports the leading theoretical concept for the origin of differential rotation as a result of interaction between convection and rotation. This talk outlines the current state of the differential rotation theory. Numerical models based on the theory reproduce the observed solar rotation quite closely. The models also compute meridional flow and predict that the flow at the bottom of the convection zone is not small compared to the surface. Theoretical predictions for stellar differential rotation as a function of the stellar mass and rotation rate are discussed and compared with observations. The implications of the differential rotation models for solar and stellar dynamos are briefly discussed.

The contribution of the Babcock-Leighton mechanism to the generation of the Sun's poloidal magnetic field is estimated from sunspot data for three solar cycles. Comparison of the derived quantities with the A-index of the large-scale magnetic field suggests a positive answer to the question posed in the title of this paper.

A combination of diamagnetic pumping and a nonlocal alpha-effect of the Babcock-Leighton type in a solar dynamo model helps to reproduce observations of solar magnetic activity. The period of the solar cycle can be reproduced without reducing magnetic diffusivity in the bulk of the convection zone below the standard mixing-length value of $10^{13}$ cm$^2$s$^{-1}$. The simulated global fields are antisymmetric about the equator and the toroidal-to-poloidal field ratio is about a thousand. The time-latitude diagrams of magnetic fields in the model without meridional flow, however, differ from observations. Only when the meridional flow is included and the alpha-effect profile peaking at mid latitudes is applied, can the observational butterfly diagrams be reproduced.

Helicity and \alpha effect driven by the nonaxisymmetric Tayler instability of toroidal magnetic fields in stellar radiation zones are computed. In the linear approximation a purely toroidal field always excites pairs of modes with identical growth rates but with opposite helicity so that the net helicity vanishes. If the magnetic background field has a helical structure by an extra (weak) poloidal component then one of the modes dominates producing a net kinetic helicity anticorrelated to the current helicity of the background field. The mean electromotive force is computed with the result that the \alpha effect by the most rapidly growing mode has the same sign as the current helicity of the background field. The \alpha effect is found as too small to drive an \alpha^{2} dynamo but the excitation conditions for an \alpha\Omega dynamo can be fulfilled for weak poloidal fields. Moreover, if the dynamo produces its own \alpha effect by the magnetic instability then problems with its sign do not arise. For all cases, however, the \alpha effect shows an extremely strong concentration to the poles so that a possible \alpha\Omega dynamo might only work at the polar regions. Hence, the results of our linear theory lead to a new topological problem for the existence of large-scale dynamos in stellar radiation zones on the basis of the current-driven instability of toroidal fields.

The first results of the solar dynamo model that allows for the diamagnetic effect of inhomogeneous turbulence and the nonlocal alpha-effect due to the rise of magnetic loops are discussed. The nonlocal alpha-effect is not subject to the catastrophic quenching related to the conservation of magnetic helicity. Given the diamagnetic pumping, the magnetic fields are concentrated near the base of the convection zone, although the distributed-type model covers the entire thickness of the convection zone. The magnetic cycle period, the equatorial symmetry of the field, its meridional drift, and the polar-to-toroidal field ratio obtained in the model are in agreement with observations. There is also some disagreement with observations pointing the ways of improving the model.

A series of stellar models of spectral type G is computed to study the rotation laws resulting from mean-field equations. The rotation laws of the slowly rotating Sun, the fast rotating MOST stars epsilon Eri and kappa1 Cet and the rapid rotators R58 and LQ Lup can easily be reproduced. We also find that differences in the depth of the convection zone cause large differences in the surface rotation law and that the extreme surface shear of HD 171488 can only be explained with a artificially shallow convection layer. We also check the thermal wind equilibrium in fast-rotating G dwarfs and find that the polar subrotation (dOmega/dz<0) is due to the barocline effect and that the equatorial superrotation (dOmega/dr>0) is due to the Lambda effect as part of the Reynolds stresses. In the bulk of the convection zones where the meridional flow is slow and smooth the thermal wind equilibrium actually holds between the centrifugal and the pressure forces. It does not hold, however, in the bounding shear layers including the equatorial region where the Reynolds stresses dominate.

The nonlocal alpha-effect of Babcock-Leighton type is not prone to the catastrophic quenching due to conservation of magnetic helicity. This is shown with a dynamo model, which jointly applies the nonlocal alpha-effect, the diamagnetic pumping, and dynamical equation for the magnetic alpha-effect. The same model shows catastrophic quenching when the alpha-effect is changed to its local formulation. The nonlocal model shows the preferred excitation of magnetic fields of dipolar symmetry, which oscillate with a period of about ten years and have a toroidal-to-polar fields ratio of about a thousand.

A new version of a numerical model of stellar differential rotation based on mean-field hydrodynamics is presented and tested by computing the differential rotation of the Sun. The model is then applied to four individual stars including two moderate and two fast rotators to reproduce their observed differential rotation quite closely. A series of models for rapidly rotating ($P_{rot} = 1$ day) stars of different masses and compositions is generated. The effective temperature is found convenient to parameterize the differential rotation: variations with metallicity, that are quite pronounced when the differential rotation is considered as a function of the stellar mass, almost disappear in the dependence of differential rotation on temperature. The differential rotation increases steadily with surface temperature to exceed the largest differential rotation observed to date for the hottest F-stars we considered. This strong differential rotation is, however, found not to be efficient for dynamos when the efficiency is estimated with the standard $C_\Omega$-parameter of dynamo models. On the contrary, the small differential rotation of M-stars is the most dynamo-efficient. The meridional flow near the bottom of the convection zone is not small compared to the flow at the top in all our computations. The flow is distributed over the entire convection zone in slow rotators but retreats to the convection zone boundaries with increasing rotation rate, to consist of two near-boundary jets in rapid rotators. The implications of the change of the flow structure for stellar dynamos are briefly discussed.

In a density-stratified turbulent medium the cross helicity <u'.B'> is considered as a result of the interaction of the velocity fluctuations and a large-scale magnetic field. By means of a quasilinear theory and by numerical simulations we find the cross helicity and the mean vertical magnetic field anti-correlated. In the high-conductivity limit the ratio of the helicity and the mean magnetic field equals the ratio of the magnetic eddy diffusivity and the (known) density scale height. The result can be used to predict that the cross helicity at the solar surface exceeds the value of 1 Gauss km/s. Its sign is anti-correlated with that of the radial mean magnetic field. Alternatively, we can use our result to determine the value of the turbulent magnetic diffusivity from observations of the cross helicity.

The stability of toroidal magnetic fields in rotating radiative stellar zones is studied for realistic values of both the Prandtl numbers. The two considered models for the magnetic geometry represent fields with odd and even symmetry with respect to the equator. In the linear theory in Boussinesq approximation the resulting complex eigenfrequency (including growth rate and drift rate) are calculated for a given radial wavenumber of a nonaxisymmetric perturbation with m=1. The ratio of the Alfven frequency, \Omega_A, to the rate of the basic rotation, \Omega, controls the eigenfrequency of the solution. For strong fields with \Omega_A > \Omega the solutions do not feel the thermal diffusion. The growth rate runs with \Omega_A and the drift rate is close to -\Omega so that the magnetic pattern will rest in the laboratory system. For weaker fields with \Omega_A < \Omega the growth rate strongly depends on the thermal conductivity. For fields with dipolar parity and for typical values of the heat conductivity the resulting very small growth rates are almost identical with those for vanishing gravity. For fields with dipolar symmetry the differential rotation of any stellar radiative zone (like the solar tachocline) is shown as basically stabilizing the instability independent of the sign of the shear. Finally, the current-driven kink-type instability of a toroidal background field is proposed as a model for the magnetism of Ap stars. The recent observation of a lower magnetic field treshold of about 300 Gauss for Ap stars is understood as corresponding to the minimum magnetic field producing the instability.

Deviations from axial symmetry are necessary to maintain self-sustained MRI-turbulence. We define the parameters region where nonaxisymmetric MRI is excited and study dependence of the unstable modes structure and growth rates on the relevant parameters. We solve numerically the linear eigenvalue problem for global axisymmetric and nonaxisymmetric modes of standard-MRI in Keplerian disks with finite diffusion. For small magnetic Prandtl number the microscopic viscosity completely drops out from the analysis so that the stability maps and the growth rates expressed in terms of the magnetic Reynolds number Rm and the Lundquist number S do not depend on the magnetic Prandtl number Pm. The minimum magnetic field for onset of nonaxisymmetric MRI grows with Rm. For given S all nonaxisymmetric modes disappear for sufficiently high Rm. This behavior is a consequence of the radial fine-structure of the nonaxisymmetric modes resulting from the winding effect of differential rotation. It is this fine-structure which presents severe resolution problems for the numerical simulation of MRI at large Rm. For weak supercritical magnetic fields only axisymmetric modes are unstable. Nonaxisymmetric modes need stronger fields and not too fast rotation. If Pm is small its real value does not play any role in MRI.

The question is addressed whether stellar differentially rotating radiative zones (like the solar tachocline) excite nonaxisymmetric r-modes which can be observed. To this end the hydrodynamical stability of latitudinal differential rotation is studied. The amount of rotational shear required for the instability is estimated in dependence of the character of radial stratification and the flow patterns excited by the instability are found. The eigenvalue equations for the nonaxisymmetric disturbances are formulated in 3D and then solved numerically. Radial displacements and entropy disturbances are included. The equations contain the 2D approximation of strictly horizontal displacements as a special limit. The critical magnitude of the latitudinal differential rotation for onset of the instability is considerably reduced in the 3D theory compared to the 2D approximation. The instability requires a subadiabatic stratification. It does not exist in the bulk of convection zone with almost adiabatic stratification but may switch on near its base in the region of penetrative convection. Growth rates and symmetry types of the modes are computed in dependence on the rotation law parameters. The S1 mode with its transequatorial toroidal vortices is predicted as the dominating instability mode. The vortices show longitudinal drift rates retrograde to the basic rotation which are close to that of the observed weak r-mode signatures at the solar surface.

The solar tachocline is shown as hydrodynamically stable against nonaxisymmetric disturbances if it is true that no cos^{4}\theta term exists in its rotation law. We also show that the toroidal field of 200 Gauss amplitude which produces the tachocline in the magnetic theory of Ruediger & Kitchatinov (1997) is stable against nonaxisymmetric MHD disturbances -- but it becomes unstable for rotation periods slightly slower than 25 days. The instability of such weak fields lives from the high thermal diffusivity of stellar radiation zones compared with the magnetic diffusivity. The growth times, however, result as very long (of order of 10\^5 rotation times). With estimations of the chemical mixing we find the maximal possible field amplitude to be ~500 Gauss in order to explain the observed lithium abundance of the Sun. Dynamos with such low field amplitudes should not be relevant for the solar activity cycle. With nonlinear simulations of MHD Taylor-Couette flows it is shown that for the rotation-dominated magnetic instability the resulting eddy viscosity is only of the order of the molecular viscosity. The Schmidt number as the ratio of viscosity and chemical diffusion grows to values of ~20. For the majority of the stellar physics applications, the magnetic-dominated Tayler instability will be quenched by the stellar rotation.

The property of inhomogeneous turbulence in conducting fluids to expel large-scale magnetic fields in the direction of decreasing turbulence intensity is shown as important for the magnetic field dynamics near the base of a stellar convection zone. The downward diamagnetic pumping confines a fossil internal magnetic field in the radiative core so that the field geometry is appropriate for formation of the solar tachocline. For the stars of solar age, the diamagnetic confinement is efficient only if the ratio of turbulent magnetic diffusivity of the convection zone to the (microscopic or turbulent) diffusivity of the radiative interiour is larger than 10^5. Confinement in younger stars require still larger diffusivity ratio. The observation of persistent magnetic structures on young solar-type stars can thus provide evidences for the nonexistence of tachoclines in stellar interiors and on the level of turbulence in radiative cores.

Stability of toroidal magnetic field in a stellar radiation zone is considered for the cases of uniform and differential rotation. In the rigidly rotating radiative core shortly below the tachocline, the critical magnetic field for instability is about 600 G. The unstable disturbances for slightly supercritical fields have short radial scales ~1 Mm. Radial mixing produced by the instability is estimated to conclude that the internal field of the sun can exceed the critical value of 600 G only marginally. Otherwise, the mixing is too strong and not compatible with the observed lithium abundance. Analysis of joint instability of differential rotation and toroidal field leads to the conclusion that axisymmetric models of the laminar solar tachocline are stable to nonaxisymmetric disturbances. The question of whether sun-like stars can posses tachoclines is addressed with positive answer for stars with rotation periods shorter than about two months.

A theoretical sunspot model is provided including magnetic suppression of the diffusivities and also a strong stratification of density and temperature. Heat diffusion alone with given magnetic field and zero mean flow only produces (after a very long relaxation time) dark spots without any bright ring. Models with full dynamics of both field and flow, however, provide rings and also the observed correlation of ring temperature excess and the spot size. The rings are formed as the result of heat transport by the resulting flow system and increased thermal diffusivity due to reduced magnetic quenching around spots.

The questions of how strong magnetic fields can be stored in rotating stellar radiative zones without being subjected to pinch-type instabilities and how much radial mixing is produced if the fields are unstable are addressed. Linear equations are derived for weak disturbances of magnetic and velocity fields which are global in horizontal dimensions but short--scaled in radius. The equations are solved to evaluate the stability of toroidal field patterns with one or two latitudinal belts under the influence of a rigid basic rotation. Hydrodynamic stability of latitudinal differential rotation is also considered. The magnetic instability is essentially three--dimensional. It does not exist in a 2D formulation with strictly horizontal disturbances on decoupled spherical shells. Only stable (magnetically modified) r-modes are found in this case. The instability recovers in 3D. The most rapidly growing modes for the Sun have radial scales smaller than 1 Mm. The finite thermal conductivity makes a strong destabilizing effect. The marginal field strength for the onset of the instability in the upper part of the solar radiative zone is about 600 G. The toroidal field can only slightly exceed this critical value for otherwise the radial mixing produced by the instability would be too strong to be compatible with the observed lithium abundance. Also the threshold for hydrodynamic instability of differential rotation which exists in 2D is lowered in 3D. When radial displacements are included, the value of 28% for critical shear is reduced to 21%.

Rather weak fossil magnetic fields in the radiative core can produce the solar tachocline if the field is poloidal and almost horizontal in the tachocline region, i.e. if the field is confined within the core. This particular field geometry is shown to result from a shallow (<~1Mm) penetration of the meridional flow existing in the convection zone into the radiative core. I.e., two conditions are crucial for a magnetic tachocline theory: (i) the presence of meridional flow of a few meters per second at the base of the convection zone, and (ii) a magnetic diffusivity inside the tachocline smaller than 10^8 cm^2s^-1. Numerical solutions for both confined poloidal fields and the resulting tachocline structures are presented. We find that the tachocline thickness runs as <I>B</I>_p^-1/2 with the poloidal field amplitude falling below 5% of the solar radius for <I>B</I>^p > 5 mG. The resulting toroidal field amplitude inside the tachocline of about 100 G does not depend on the <I>B</I>^p. The hydromagnetic stability of the tachocline is only briefly discussed. For the hydrodynamic stability of latitudinal differential rotation we found that the critical 29% of the 2D theory of Watson (1981) are reduced to only 21% in 3D for marginal modes of about 6 Mm radial scale.

A plane-shear flow in a fluid with forced turbulence is considered. If the fluid is electrically-conducting then a mean electromotive force (EMF) results even without basic rotation and the magnetic diffusivity becomes a highly anisotropic tensor. It is checked whether in this case self-excitation of a large-scale magnetic field is possible (so-called WxJ-dynamo) and the answer is NO. The calculations reveal the cross-stream components of the EMF perpendicular to the mean current having the wrong signs, at least for small magnetic Prandtl numbers. After our results numerical simulations with magnetic Prandtl number of about unity have only a restricted meaning as the Prandtl number dependence of the diffusivity tensor is rather strong. If, on the other hand, the turbulence field is stratified in the vertical direction then a dynamo active alpha-effect is produced. The critical magnetic Reynolds number for such a self-excitation in a simple shear flow is slightly above 10 like for the other - but much more complicated - flow patterns used in existing dynamo experiments with liquid sodium or gallium.

We show that the MHD theory that explains the solar tachocline by an effect of the magnetic field can work with the decay modes of a fossil field in the solar interior if the meridional flow of the convection zone penetrates slightly the radiative zone beneath. An equatorward flow of about 10 m/s penetrating to a maximum depth of 1000 km below the convection zone is able to generate almost horizontal field lines in the tachocline region so that the internal field is almost totally confined to the radiative zone. The theory of differential solar rotation indeed provides meridional flows of about 10 m/s and a penetration depth of < 1000 km for viscosity values that are characteristic of a stable tachocline.

The influence of the basic rotation on anisotropic and inhomogeneous turbulence is discussed in the context of differential rotation theory. An improved representation for the original turbulence leads to a $\Lam$-effect which complies with the results of 3D numerical simulations. The resulting rotation law and meridional flow agree well with both the surface observations ($\partial\Om/\partial r<0$ and meridional flow towards the poles) and with the findings of helioseismology. The computed equatorward flow at the bottom of convection zone has an amplitude of about 10 m/s and may be significant for the solar dynamo. The depth of the meridional flow penetration into the radiative zone is proportional to $\nu^{0.5}_\mathrm{core}$, where $\nu_\mathrm{core}$ is the viscosity beneath the convection zone. The penetration is very small if the tachocline is laminar.

The differential rotation of anti-solar type detected by observations for several stars may result from a fast meridional flow. The sufficiently intensive meridional circulation may be caused by large-scale thermal inhomogeneities or, perhaps, by tidal forcing from a companion star. First results of simulations of the anti-solar rotation of a giant star with magnetically induced thermal inhomogeneities are presented. Perspectives for observational check of the theory are discussed.

The influence of the Hall effect on the global stability of cool Kepler disks under the influence of an axial magnetic field is considered. For sufficiently large magnetic Reynolds numbers Rm the magnetorotational instability (MRI) exists in a finite interval of magnetic field amplitudes, $B$_{min} < $B$ < $B$_{max}. For Kepler disks the pure MRI needs both rather high Rm (representing the needed electrical conductivity) as well as $B$_{min} of order 0.1 G. The magnetic field pattern resulting from our global and linear calculations is of quadrupolar parity. For magnetic fields $antiparallel$ to the rotation axis the Hall effect reduces the minimum magnetic Reynolds number by about one order of magnitude. The $B$_{min}, however, is even (sightly) increased (see Fig. 6). For magnetic fields $parallel$ to the rotation axis the Hall effect drives its own instability without the action of the Lorentz force. The corresponding critical magnetic Reynolds number proves to be larger with Hall effect (Rm ~ 10) than without Hall effect (Rm ~ 7) so that the Hall effect for parallel fields even disturbs the formation of MHD-instability in cool protoplanetary disks. If the disk is supercritical then the main result of the Hall effect for positive fields is the strong reduction of the minimum magnetic field amplitude which is necessary to start the instability. Observations must show whether in star-forming regions the rotation axis and the magnetic field orientation are correlated or are anticorrelated. If the magnetic fields are high enough then our model predicts the dominance of fields antiparallel to the rotation axis.

The three components of the heat-flux vector $F =\rho C_p <u^\prime T^\prime> are numerically computed for a stratified rotating turbulent convection using the NIRVANA code in a flat box. The latitudinal component $F_\theta$ proves to be negative (positive) in the northern (southern) hemisphere so that the heat always flows towards the poles. As a surprise, the radial heat-flux $F_r$ peaks at the equator rather than at the poles (Taylor numbers O(10^6)). The same behavior is observed for the radial turbulence intensity $<u_r^{\prime 2}>$ which for \emph{free} turbulence is also believed to peak at the poles (see Eq. (19) below). As we can show, however, the consequences of this unexpected result (also obtained by Kaepylae, Korpi and Tuominen 2004) for the theory of differential rotation are small as mainly the $F_\theta$ is responsible to solve the `Taylor number puzzle'. In all our simulations the azimuthal component $F_\phi$ proves to be negative so that the rotating turbulence produces an westwards directed azimuthal heat-flux which should be observable. Fluctuations with higher temperature are expected to be anticorrelated with their own angular velocity fluctuations. We find this rotation-induced result as understandable as the $F_\phi$ is closely related to the radial \Lambda-effect which is known to be also negative in stratified and rapidly rotating convection zones.

A linear but global numerical model for the magnetorotational instability (MRI) in disk geometry is considered to estimate the instability parameters for galaxies. Similarity rules suggested by a local analysis are applied to reveal a universal behavior approached by the results of global calculations for large magnetic Prandtl numbers. The findings are used to estimate the MRI characteristics for galaxies with their very large magnetic Prandtl numbers which cannot be attained in any numerical simulations. The resulting minimum field for the instability, Bmin 10-25 G, is small compared to any seed fields currently discussed. The growth times of MRI are estimated to be on the order of the rotation period of the inner rigidly-rotating core, i.e. ~100 Myr. Global MRI excites preferentially the magnetic field modes of quadrupolar symmetry.