The possibility of explaining shear flow dynamos in terms of a magnetic shear-current effect (SC) is examined. A competing explanation is the incoherent $\alpha$-shear dynamo effect. Our primary diagnostics is the turbulent magnetic diffusivity tensor, in particular the off-diagonal diffusivity tensor component $\eta_{yx}$, a systematically negative sign of which would imply coherent dynamo action through SC. To be able to measure turbulent transport coefficients from systems with strong magnetic fluctuations, we present an extension of the test-field method that is capable of such measurements in the case where the pressure gradient term is dropped from the MHD equations, which is a nonlinear TFM (NLTFM). The hydrodynamic equation is related to Burger's equation and the resulting flows are referred to as magnetized burgulence. We use both kinetic and magnetic forcings, to mimic cases without and with simultaneous small-scale dynamo action. When the simplified MHD flow is forced kinetically, negative $\eta_{yx}$ values are obtained with exponential growth in both the radial and azimuthal magnetic field components. Using stochastic monochromatic magnetic forcing, the exponential growth is no longer seen, and NLTFM yields positive values of $\eta_{yx}$. In an attempt to recover the exponential growth seen in the kinetically forced case, we also employ an alternative forcing function with lowest $k$ vectors being removed - in this case the exponential growth is recovered, but the NLTFM results do not change dramatically. Our analysis of the dynamo numbers for the coherent SC and incoherent $\alpha$ and SC effects show that the incoherent effects are the main drivers of the dynamo instability in majority of cases. The coherent SC effect is mildly enhancing the dynamo action in the kinetically forced cases, while we find no evidence for magnetic SC in our simulations.

(abridged) Context: Turbulent diffusion of large-scale flows and magnetic fields play major roles in many astrophysical systems. Aims: Our goal is to compute turbulent viscosity and magnetic diffusivity, relevant for diffusing large-scale flows and magnetic fields, respectively, and their ratio, the turbulent magnetic Prandtl number, ${\rm Pm}_{\rm t}$, for isotropically forced homogeneous turbulence. Methods: We use simulations of forced turbulence in fully periodic cubes composed of isothermal gas with an imposed large-scale sinusoidal shear flow. Turbulent viscosity is computed either from the resulting Reynolds stress or from the decay rate of the large-scale flow. Turbulent magnetic diffusivity is computed using the test-field method. The scale dependence of the coefficients is studied by varying the wavenumber of the imposed sinusoidal shear and test fields. Results: We find that turbulent viscosity and magnetic diffusivity are in general of the same order of magnitude. Furthermore, the turbulent viscosity depends on the fluid Reynolds number (${\rm Re}$) and scale separation ratio of turbulence. The scale dependence of the turbulent viscosity is found to be well approximated by a Lorentzian. The results for the turbulent transport coefficients appear to converge at sufficiently high values of ${\rm Re}$ and the scale separation ratio. However, a weak decreasing trend is found even at the largest values of ${\rm Re}$. The turbulent magnetic Prandtl number converges to a value that is slightly below unity for large ${\rm Re}$ whereas for small ${\rm Re}$, we find values between 0.5 and 0.6. Conclusions: The turbulent magnetic diffusivity is in general consistently higher than the turbulent viscosity. The actual value of ${\rm Pm}_{\rm t}$ found from the simulations ($\approx0.9\ldots0.95$) at large ${\rm Re}$ and scale separation ratio is higher than any of the analytic predictions.

Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions have been obtained so far. In this paper we present a semi-global dynamo simulation in the transition region, exhibiting two coexisting dynamo modes, a cyclic and a stationary one, both being dynamically significant. We seek to understand how such a dynamo is driven by analyzing the large-scale flow properties (DR and meridional circulation) together with the turbulent transport coefficients obtained with the test-field method. Neither an $\alpha\Omega$ dynamo wave nor an advection-dominated dynamo are able to explain the cycle period and the propagation direction of the mean magnetic field. Furthermore, we find that the $\alpha$ effect is comparable or even larger than the $\Omega$ effect in generating the toroidal magnetic field, and therefore, the dynamo seems to be $\alpha^2\Omega$ or $\alpha^2$ type. We further find that the effective large-scale flows are significantly altered by turbulent pumping.

Recent numerical work in helioseismology has shown that a periodically varying subsurface magnetic field leads to a fanning of the $f$-mode, which emerges from the density jump at the surface. In an attempt to model a more realistic situation, we now modulate this periodic variation with an envelope, giving thus more emphasis on localized bipolar magnetic structures in the middle of the domain. Some notable findings are: (i) compared to the purely hydrodynamic case, the strength of the $f$-mode is significantly larger at high horizontal wavenumbers $k$, but the fanning is weaker for the localized subsurface magnetic field concentrations investigated here than the periodic ones studied earlier; (ii) when the strength of the magnetic field is enhanced at a fixed depth below the surface, the fanning of the $f$-mode in the $k\omega$ diagram increases proportionally in such a way that the normalized $f$-mode strengths remain nearly the same in different such cases; (iii) the unstable Bloch modes reported previously in case of harmonically varying magnetic fields are now completely absent when more realistic localized magnetic field concentrations are imposed beneath the surface, thus suggesting that the Bloch modes are unlikely to be supported during most phases of the solar cycle; (iv) the $f$-mode strength appears to depend also on the depth of magnetic field concentrations such that it shows a relative decrement when the maximum of the magnetic field is moved to a deeper layer. We argue that detections of $f$-mode perturbations such as those being explored here could be effective tracers of solar magnetic fields below the photosphere before these are directly detectable as visible manifestations in terms of active regions or sunspots.

We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($\phi$) averaged fields to determine the 27 turbulent transport coefficients of the electromotive force, of which six are related to the $\alpha$ tensor. This method has previously been used either in simulations in Cartesian coordinates or in the geodynamo context and is applied here for the first time to fully compressible simulations of solar-like dynamos. We find that the $\phi\phi$-component of the $\alpha$ tensor does not follow the profile expected from that of kinetic helicity. The turbulent pumping velocities significantly alter the effective mean flows acting on the magnetic field and therefore challenge the flux transport dynamo concept. All coefficients are significantly affected by dynamically important magnetic fields. Quenching as well as enhancement are being observed. This leads to a modulation of the coefficients with the activity cycle. The temporal variations are found to be comparable to the time-averaged values and seem to be responsible for a nonlinear feedback on the magnetic field generation. Furthermore, we quantify the validity of the Parker-Yoshimura rule for the equatorward propagation of the mean magnetic field in the present case.

We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kramers-like opacity prescription as a function of temperature and density or a static profile of a similar shape. We show that the subadiabatic layer arises due to nonlocal energy transport by buoyantly driven downflows in the upper parts of the CZ. Analysis of the force balance of the upflows and downflows confirms that convection is driven by cooling at the surface. We find that the commonly used prescription for the convective enthalpy flux being proportional to the negative entropy gradient does not hold in the stably stratified layers where the flux is positive. We demonstrate the existence of a non-gradient contribution to the enthalpy flux, which is estimated to be important throughout the convective layer. A quantitative analysis of downflows indicates a transition from a tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a height-independent number of strong downdrafts persist. This change of flow topology occurs when a substantial subadiabatic layer is present in the lower part of the CZ.

With the ultimate aim of using the fundamental or $f$-mode to study helioseismic aspects of turbulence-generated magnetic flux concentrations, we use randomly forced hydromagnetic simulations of a piecewise isothermal layer in two dimensions with reflecting boundaries at top and bottom. We compute numerically diagnostic wavenumber-frequency diagrams of the vertical velocity at the interface between the denser gas below and the less dense gas above. For an Alfv\'en-to-sound speed ratio of about 0.1, a 5% frequency increase of the $f$-mode can be measured when $k_xH_{\rm p}=3$-$4$, where $k_x$ is the horizontal wavenumber and $H_{\rm p}$ is the pressure scale height at the surface. Since the solar radius is about 2000 times larger than $H_{\rm p}$, the corresponding spherical harmonic degree would be 6000-8000. For weaker fields, a $k_x$-dependent frequency decrease by the turbulent motions becomes dominant. For vertical magnetic fields, the frequency is enhanced for $k_xH_{\rm p}\approx4$, but decreased relative to its nonmagnetic value for $k_xH_{\rm p}\approx9$.

We compare simple analytical closure models of homogeneous turbulent Boussinesq convection for stellar applications with three-dimensional simulations. We use simple analytical closure models to compute the fluxes of angular momentum and heat as a function of rotation rate measured by the Taylor number. We also investigate cases with varying angles between the angular velocity and gravity vectors, corresponding to locating the computational domain at different latitudes ranging from the pole to the equator of the star. We perform three-dimensional numerical simulations in the same parameter regimes for comparison. The free parameters appearing in the closure models are calibrated by two fitting methods using simulation data. Unique determination of the closure parameters is possible only in the non-rotating case or when the system is placed at the pole. In the other cases the fit procedures yield somewhat differing results. The quality of the closure is tested by substituting the resulting coefficients back into the closure model and comparing with the simulation results. To eliminate the possibilities that the results obtained depend on the aspect ratio of the simulation domain or suffer from too small Rayleigh numbers we performed runs varying these parameters. The simulation data for the Reynolds stress and heat fluxes broadly agree with previous compressible simulations. The closure works fairly well with slow and fast rotation but its quality degrades for intermediate rotation rates. We find that the closure parameters depend not only on rotation rate but also on latitude. The weak dependence on Rayleigh number and the aspect ratio of the domain indicates that our results are generally valid

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

We show that in the presence of a harmonically varying magnetic field the fundamental or $f$-mode in a stratified layer is altered in such a way that it fans out in the diagnostic $k\omega$ diagram, but with mode power also within the fan. In our simulations, the surface is defined by a temperature and density jump in a piecewise isothermal layer. Unlike our previous work (Singh et al. 2014) where a uniform magnetic field was considered, we employ here a nonuniform magnetic field together with hydromagnetic turbulence at length scales much smaller than those of the magnetic fields. The expansion of the $f$-mode is stronger for fields confined to the layer below the surface. In some of those cases, the $k\omega$ diagram also reveals a new class of low frequency vertical stripes at multiples of twice the horizontal wavenumber of the background magnetic field. We argue that the study of the $f$-mode expansion might be a new and sensitive tool to determining subsurface magnetic fields with longitudinal periodicity.

We analyze the nature of dynamo action that produces horizontally averaged magnetic fields in two particular flows that were studied by Roberts (1972, Phil. Trans. R. Soc. A 271, 411), namely his flows II and III. They have zero kinetic helicity either pointwise (flow II), or on average (flow III). Using direct numerical simulations, we determine the onset conditions for dynamo action at moderate values of the magnetic Reynolds number. Using the test-field method, we show that the turbulent magnetic diffusivity is then positive for both flows. However, we demonstrate that for both flows large-scale dynamo action occurs through delayed transport. Mathematically speaking, the magnetic field at earlier times contributes to the electromotive force through the off-diagonal components of the alpha tensor such that a zero mean magnetic field becomes unstable to dynamo action. This represents a qualitatively new mean-field dynamo mechanism not previously described.

Direct numerical simulations of isotropically forced homogeneous stationary turbulence with an imposed passive scalar concentration gradient are compared with an analytical closure model which provides evolution equations for the mean passive scalar flux and variance. Triple correlations of fluctuations appearing in these equations are described in terms of relaxation terms proportional to the quadratic correlations. Three methods are used to extract the relaxation timescales tau_i from direct numerical simulations. Firstly, we insert the closure ansatz into our equations, assume stationarity, and solve for tau_i. Secondly, we use only the closure ansatz itself and obtain tau_i from the ratio of quadratic and triple correlations. Thirdly we remove the imposed passive scalar gradient and fit an exponential decay law to the solution. We vary the Reynolds (Re) and P\'eclet (Pe) numbers while keeping their ratio at unity and the degree of scale separation and find for large Re fair correspondence between the different methods. The ratio of the turbulent relaxation time of passive scalar flux to the turnover time of turbulent eddies is of the order of three, which is in remarkable agreement with earlier work. Finally we make an effort to extract the relaxation timescales relevant for the viscous and diffusive effects. We find two regimes which are valid for small and large Re, respectively, but the dependence of the parameters on scale separation suggests that they are not universal.

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.

Certain aspects of the mean-field theory of turbulent passive scalar transport and of mean-field electrodynamics are considered with particular emphasis on aspects of compressible fluids. It is demonstrated that the total mean-field diffusivity for passive scalar transport in a compressible flow may well be smaller than the molecular diffusivity. This is in full analogy to an old finding regarding the magnetic mean-field diffusivity in an electrically conducting turbulently moving compressible fluid. These phenomena occur if the irrotational part of the motion dominates the vortical part, the P\`eclet or magnetic Reynolds number is not too large, and, in addition, the variation of the flow pattern is slow. For both the passive scalar and the magnetic cases several further analytical results on mean-field diffusivities and related quantities found within the second-order correlation approximation are presented, as well as numerical results obtained by the test-field method, which applies independently of this approximation. Particular attention is paid to non-local and non-instantaneous connections between the turbulence-caused terms and the mean fields. Two examples of irrotational flows, in which interesting phenomena in the above sense occur, are investigated in detail. In particular, it is demonstrated that the decay of a mean scalar in a compressible fluid under the influence of these flows can be much slower than without any flow, and can be strongly influenced by the so-called memory effect, that is, the fact that the relevant mean-field coefficients depend on the decay rates themselves.

When scale separation in space and time is poor, the alpha effect and turbulent diffusivity have to be replaced by integral kernels. Earlier work in computing these kernels using the test-field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean-field equations are solved for oscillatory alpha-shear dynamos as well as alpha^2 dynamos in which alpha is antisymmetric about the equator, making this dynamo also oscillatory. In both cases, the critical values of the dynamo number is lowered by the fact that the dynamo is oscillatory.

Evidence for a parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of an instability in magnetohydrodynamics appears to be the first of its kind and is similar to chiral symmetry breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are being discussed.

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.

Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, although the {\alpha} - {\Omega} mode is generally expected to dominate because it is the fastest growing. In an {\alpha} - {\Omega} dynamo, most of the field amplification is produced by the shear. As differential rotation is an ubiquitous source of shear in astrophysics, such dynamos are believed to be the source of most astrophysical large-scale magnetic fields. Aims. We study the stability of oscillatory migratory {\alpha} - {\Omega} type dynamos in turbulence simulations. Methods. We use shearing-sheet simulations of hydromagnetic turbulence that is helically forced at a wavenumber that is about three times larger than the lowest wavenumber in the domain so that both {\alpha} - {\Omega} and {\alpha}2 dynamo action is possible. Results. After initial dominance and saturation, the {\alpha} - {\Omega} mode is found to be destroyed by an orthogonal {\alpha}2 mode sustained by the helical turbulence alone. We show that there are at least two processes through which this transition can occur. Conclusions. The fratricide of {\alpha} - {\Omega} dynamos by its {\alpha}2 sibling is discussed in the context of grand minima of solar and stellar activity. However, the genesis of {\alpha} - {\Omega} dynamos from an {\alpha}2 dynamo has not yet been found.

A strong toroidal field can exist in form of a magnetic layer in the overshoot region below the solar convection zone. This motivates a more detailed study of the magnetic buoyancy instability with rotation. We calculate the alpha effect due to helical motions caused by a disintegrating magnetic layer in a rotating density-stratified system with angular velocity Omega making an angle theta with the vertical. We also study the dependence of the alpha effect on theta and the strength of the initial magnetic field. We carry out three-dimensional hydromagnetic simulations in Cartesian geometry. A turbulent EMF due to the correlations of the small scale velocity and magnetic field is generated. We use the test-field method to calculate the transport coefficients of the inhomogeneous turbulence produced by the layer. We show that the growth rate of the instability and the twist of the magnetic field vary monotonically with the ratio of thermal conductivity to magnetic diffusivity. The resulting alpha effect is inhomogeneous and increases with the strength of the initial magnetic field. It is thus an example of an "anti-quenched" alpha effect. The alpha effect is nonlocal, requiring around 8--16 Fourier modes to reconstruct the actual EMF based on the actual mean field.

The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha effect and turbulent diffusion.

Aims: The test-field method for computing turbulent transport coefficients from simulations of hydromagnetic flows is extended to the regime with a magnetohydrodynamic (MHD) background. Methods: A generalized set of test equations is derived using both the induction equation and a modified momentum equation. By employing an additional set of auxiliary equations, we derive linear equations describing the response of the system to a set of prescribed test fields. Purely magnetic and MHD backgrounds are emulated by applying an electromotive force in the induction equation analogously to the ponderomotive force in the momentum equation. Both forces are chosen to have Roberts flow-like geometry. Results: Examples with an MHD background are studied where the previously used quasi-kinematic test-field method breaks down. In cases with homogeneous mean fields it is shown that the generalized test-field method produces the same results as the imposed-field method, where the field-aligned component of the actual electromotive force from the simulation is used. Furthermore, results for the turbulent diffusivity tensor are given, which are inaccessible to the imposed-field method. For MHD backgrounds, new mean-field effects are found that depend on the occurrence of cross-correlations between magnetic and velocity fluctuations. For strong imposed fields, $\alpha$ is found to be quenched proportional to the fourth power of the field strength, regardless of the type of background studied.

The effect of a dynamo-generated mean magnetic field of Beltrami type on the mean electromotive force is studied. In the absence of the mean magnetic field the turbulence is assumed to be homogeneous and isotropic, but it becomes inhomogeneous and anisotropic with this field. Using the testfield method the dependence of the alpha and turbulent diffusivity tensors on the magnetic Reynolds number Rm is determined for magnetic fields that have reached approximate equipartition with the velocity field. The tensor components are characterized by a pseudoscalar alpha and a scalar turbulent magnetic diffusivity etat. Increasing Rm from 2 to 600 reduces etat by a factor ~5, suggesting that the quenching of etat is, in contrast to the 2-dimensional case, only weakly dependent on Rm. Over the same range of Rm, however, alpha is reduced by a factor ~14, which can qualitatively be explained by a corresponding increase of a magnetic contribution to the alpha effect with opposite sign. The level of fluctuations of alpha and etat is only 10% and 20% of the respective kinematic reference values.

The turbulent magnetic diffusivity tensor is determined in the presence of rotation or shear. The question is addressed whether dynamo action from the shear-current effect can explain large-scale magnetic field generation found in simulations with shear. For this purpose a set of evolution equations for the response to imposed test fields is solved with turbulent and mean motions calculated from the momentum and continuity equations. The corresponding results for the electromotive force are used to calculate turbulent transport coefficients. The diagonal components of the turbulent magnetic diffusivity tensor are found to be very close together, but their values increase slightly with increasing shear and decrease with increasing rotation rate. In the presence of shear, the sign of the two off-diagonal components of the turbulent magnetic diffusion tensor is the same and opposite to the sign of the shear. This implies that dynamo action from the shear--current effect is impossible, except perhaps for high magnetic Reynolds numbers. However, even though there is no alpha effect on the average, the components of the alpha tensor display Gaussian fluctuations around zero. These fluctuations are strong enough to drive an incoherent alpha--shear dynamo. The incoherent shear--current effect, on the other hand, is found to be subdominant.

There are various analytical approaches to the mean electromotive force $\cal E =< \vec{u}\times \vec{b}>$ crucial in mean--field electrodynamics, with $\vec{u}$ and $\vec{b}$ being velocity and magnetic field fluctuations. In most cases the traditional approach, restricted to the second--order correlation approximation, has been used. Its validity is only guaranteed for a range of conditions, which is narrow in view of many applications, e.g., in astrophysics. With the intention to have a wider range of applicability other approaches have been proposed which make use of the so--called $\tau$--approximation, reducing correlations of third order in $\vec{u}$ and $\vec{b}$ to such of second order. After explaining some basic features of the traditional approach a critical analysis of the approaches of that kind is given. It is shown that they lead in some cases to results which are in clear conflict with those of the traditional approach. It is argued that this indicates shortcomings of the $\tau$--approaches and poses serious restrictions to their applicability. These shortcomings do not result from the basic assumption of the $\tau$--approximation. Instead, they seem to originate in some simplifications made in order to derive $\cal E$ without really solving the equations governing $\vec{u}$ and $\vec{b}$. A starting point for a new approach is described which avoids the conflict.

A comparison is made between mean-field models and direct numerical simulations of rotating magnetoconvection and the geodynamo. The mean-field coefficients are calculated with the fluid velocity taken from the direct numerical simulations. The magnetic fields resulting from mean-field models are then compared with the mean magnetic field from the direct numerical simulations.

We study the stability/establishment of dipolar magnetostatic equilibrium configurations in new--born neutron stars (NSs) in dependence on the rotational velocity $\Omega$ and on the initial angle $\alpha$ between rotation and magnetic axis. The NS is modeled as a sphere of a highly magnetized ($B \sim 10^{15}$G) incompressible fluid of uniform density which rotates rigidly. For the initial dipolar background magnetic field, which defines the magnetic axis, two different configurations are assumed. We solve the 3D non--linear MHD equations by use of a spectral code. The problem in dimensionless form is completely defined by the initial field strength (for a fixed field geometry), the magnetic Prandtl number $\Pm$, and the normalized rotation rate. The evolution of the magnetic and velocity fields is considered for initial magnetic field strengths characterized by the ratio of ohmic diffusion and initial \Alf{} travel times $\ttOhm/\ttAO \approx 1000$, for $\Pm = 0.1, 1, 10$, and the ratio of rotation period and initial \Alf{} travel time, $P/\ttAO = 0.012, 0.12, 1.2, 12$. We find hints for the existence of a unique stable dipolar magnetostatic configuration for any specific $\alpha$, independent of the initial field geometry. Comparing NSs possessing the same field structure at the end of their proto--NS phase, it turns out that sufficiently fast rotating NSs ($P\la6 $ms) with $\alpha \la 45^0$ retain their magnetar field, while the others lose almost all of their initial magnetic energy by transferring it into magnetic and kinetic energy of relatively small--scaled fields and continue their life as radio pulsars with a dipolar surface field of $10^{12...13}$G.

We study convective motions taken from hydrodynamic simulations of rotating proto--neutron stars (PNSs) with respect to their ability to excite a dynamo instability which may be responsible for the giant neutron star magnetic fields. Since it is impossible to simulate the magnetic field evolution employing the actual magnetic Reynolds numbers ($\Rm$) resulting from the hydrodynamic simulations, (smallest) critical $\Rm$s and the corresponding field geometries are derived on the kinematic level by rescaling the velocity amplitudes. It turns out that the actual values of $\Rm$ are by many orders of magnitude larger than the critical values found. A dynamo might therefore start to act vigorously very soon after the onset of convection. But as in general dynamo growth rates are non--monotonous functions of $\Rm$ the later fate of the magnetic field is uncertain. Hence, no reliable statements on the existence and efficiency of PNS dynamos can be drawn without considering the interplay of magnetic field and convection from the beginning. Likewise, in so far as convection inside the PNS is regarded to be essential in re--launching the supernova explosion, a revision of its role in this respect could turn out to be necessary.

In the paper we comment on (R\"udiger & Shalybkov, Phys. Rev. E. 69, 016303 (2004) (RS)), the instability of the Taylor--Couette flow interacting with a homogeneous background field subject to Hall effect is studied. We correct a falsely generalizing interpretation of results presented there which could be taken to disprove the existence of the Hall--drift induced magnetic instability described in Rheinhardt and Geppert, Phys. Rev. Lett. 88, 101103. It is shown that in contrast to what is suggested by RS, no additional shear flow is necessary to enable such an instability with a non--potential magnetic background field, whereas for a curl--free one it is. In the latter case, the instabilities found in RS in situations where neither a hydrodynamic nor a magneto--rotational instability exists are demonstrated to be most likely magnetic instead of magnetohydrodynamic. Further, some minor inaccuracies are clarified.

In former papers we showed that during the decay of a neutron star's magnetic field under the influence of the Hall--drift, an unstable rise of small--scale field structures at the expense of the large--scale background field may happen. This linear stability analysis was based on the assumption of a uniform density throughout the neutron star crust, whereas in reality the density and all transport coefficients vary by many orders of magnitude. Here, we extend the investigation of the Hall--drift induced instability by considering realistic profiles of density and chemical composition, as well as background fields with more justified radial profiles. Two neutron star models are considered differing primarily in the assumption on the core matter equation of state. For their cooling history and radial profiles of density and composition we use known results to infer the conductivity profiles. These were fed into linear calculations of a dipolar field decay starting from various initial configurations. At different stages of the decay, snapshots of the magnetic fields at the equator were taken to yield background field profiles for the stability analysis. The main result is that the Hall instability may really occur in neutron star crusts. Characteristic growth times are in the order of \lesssim 10^4 ... 10^6 yrs depending on cooling age and background field strength. The influence of the equation of state and of the initial field configuration is discussed.

There is growing evidence, based on both X-ray and radio observations of isolated neutron stars, that besides the large--scale (dipolar) magnetic field, which determines the pulsar spin--down behaviour, small--scale poloidal field components are present, which have surface strengths one to two orders of magnitude larger than the dipolar component. We argue in this paper that the Hall--effect can be an efficient process in producing such small--scale field structures just above the neutron star surface. It is shown that due to a Hall--drift induced instability, poloidal magnetic field structures can be generated from strong subsurface toroidal fields, which are the result of either a dynamo or a thermoelectric instability acting at early times of a neutron star's life. The geometrical structure of these small--scale surface anomalies of the magnetic field resembles that of some types of ``star--spots''. The magnetic field strength and the length--scales are comparable with values that can be derived from various observations.

There exists both theoretical and observational evidence that the magnetic field decay in neutron stars may proceed in a pronounced non--linear way during a certain episode of the neutron star's life. In the presence of a strong magnetic field the Hall--drift dominates the field evolution in the crust and/or the superfluid core of neutron stars. Analysing observations of $P$ and $\dot{P}$ for sufficiently old isolated pulsars we gain strong hints for a significantly non--linear magnetic field decay. Under certain conditions with respect to the geometry and strength of a large--scale magnetic background field an instability is shown to occur which rapidly raises small--scale magnetic field modes. Their growth rates increase with the background field strength and may reach $\sim 10^4$ times the ohmic decay rate. Consequences for the rotational and thermal evolution as well as for the cracking of the crust of neutron stars are discussed.

In the presence of a strong magnetic field and under conditions as realized in the crust and the superfluid core of neutron stars the Hall--drift dominates the field evolution. We show by a linear analysis that for a sufficiently strong large--scale background field depending at least quadratically on position in a plane conducting slab an instability occurs which rapidly generates small--scale fields. Their growth rates depend on the choice of the boundary conditions, increase with the background field strength and may reach $10^3$ times the ohmic decay rate. The effect of that instability on the rotational and thermal evolution of neutron stars is discussed.