We report in this paper phase-resolved spectropolarimetric observations of the rapidly-rotating fully-convective M4 dwarf V374 Peg, on which a strong, mainly axisymmetric, large-scale poloidal magnetic field was recently detected. In addition to the original data set secured in 2005 August, we present here new data collected in 2005 September and 2006 August. From the rotational modulation of unpolarised line profiles, we conclude that starspots are present at the surface of the star, but their contrast and fractional coverage are much lower than those of non-fully convective active stars with similar rotation rate. Applying tomographic imaging on each set of circularly polarised profiles separately, we find that the large-scale magnetic topology is remarkably stable on a timescale of 1 yr; repeating the analysis on the complete data set suggests that the magnetic configuration is sheared by very weak differential rotation (about 1/10th of the solar surface shear) and only slightly distorted by intrinsic variability. This result is at odds with various theoretical predictions, suggesting that dynamo fields of fully-convective stars should be mostly non-axisymmetric unless they succeed at triggering significant differential rotation.

Stellar magnetic fields can be investigated using several, very complementary approaches. While conventional spectroscopy is capable of estimating the average magnetic strength of potentially complex field configurations thanks to its low sensitivity to the vector properties of the field, spectropolarimetry can be used to map the medium- and large-scale structure of magnetic topologies. In particular, the latter approach allows one to retrieve information about the poloidal and toroidal components of the large-scale dynamo fields in low-mass stars, and thus to investigate the physical processes that produce them. Similarly, this technique can be used to explore how magnetic fields couple young stars to their massive accretion disc and thus to estimate how much mass and angular momentum are transfered to the newly-born low-mass star. We present here the latest results in this field obtained with spectropolarimetry, with special emphasis on the surprising discoveries obtained on very-low mass fully-convective stars and classical T Tauri stars thanks to the ESPaDOnS spectropolarimeter recently installed on the 3.6m Canada-France-Hawaii Telescope.

Magnetohydrodynamic simulations of fully convective, rotating spheres with volume heating near the center and cooling at the surface are presented. The dynamo-generated magnetic field saturates at equipartition field strength near the surface. In the interior, the field is dominated by small-scale structures, but outside the sphere by the global scale. Azimuthal averages of the field reveal a large-scale field of smaller amplitude also inside the star. The internal angular velocity shows some tendency to be constant along cylinders and is ``anti-solar'' (fastest at the poles and slowest at the equator).

An efficient algorithm for calculating radiative transfer on massively parallel computers using domain decomposition is presented. The integral formulation of the transfer equation is used to divide the problem into a local but compute-intensive part for calculating the intensity and optical depth integrals, and a nonlocal part for communicating the intensity between adjacent processors. The waiting time of idle processors during the nonlocal communication part does not have a severe impact on the scaling. The wall clock time thus scales nearly linearly with the inverse number of processors.

In a previous paper (Ouyed et al. 2004) we presented a new model for soft gamma-ray repeaters (SGR), based on the onset of colour superconductivity in quark stars. In this model, the bursts result from the reorganization of the exterior magnetic field following the formation of vortices that confine the internal magnetic field (the \mz effect). Here we extend the model by presenting full 3-dimensional simulations of the evolution of the inclined exterior magnetic field immediately following vortex formation. The simulations capture the violent reconnection events in the entangled surface magnetic field as it evolves into a smooth, more stable, configuration which consists of a dipole field aligned with the star's rotation axis. The total magnetic energy dissipated in this process is found to be of the order of $10^{44} \mathrm{erg}$ and, if it is emitted as synchrotron radiation, peaks typically at $280 \mathrm{keV}$. The intensity decays temporally in a way resembling SGRs and AXPs (anomalous X-ray pulsars), with a tail lasting from a few to a few hundred times the rotation period of the star, depending on the initial inclination between the rotation and dipole axis. One of the obvious consequences of our model's final state (aligned rotator) is the suppression of radio-emission in SGRs and AXPs following their bursting era. We suggest that magnetar-like magnetic field strength alone cannot be responsible for the properties of SGRs and AXPs, while a quark star entering the ``Meissner phase'' is compatible with the observational facts. We compare our model to observations and highlight our predictions.

Nonhelical hydromagnetic forced turbulence is investigated using large scale simulations on up to 256 processors and $1024^3$ meshpoints. The magnetic Prandtl number is varied between 1/8 and 30, although in most cases it is unity. When the magnetic Reynolds number is based on the inverse forcing wavenumber, the critical value for dynamo action is shown to be around 35 for magnetic Prandtl number of unity. For small magnetic Prandtl numbers we find the critical magnetic Reynolds number to increase with decreasing magnetic Prandtl number. The Kazantsev $k^{3/2}$ spectrum for magnetic energy is confirmed for the kinematic regime, i.e. when nonlinear effects are still unimportant and when the magnetic Prandtl number is unity. In the nonlinear regime, the energy budget converges for large Reynolds numbers (around 1000) such that for our parameters about 70% is in kinetic energy and about 30% is in magnetic energy. The energy dissipation rates are converged to 30% viscous dissipation and 70% resistive dissipation. Second order structure functions of the Elsasser variables give evidence for a $k^{-5/3}$ spectrum. Nevertheless, the three-dimensional spectrum is close to $k^{-3/2}$, but we argue that this is due to the bottleneck effect. The bottleneck effect is shown to be equally strong both for magnetic and nonmagnetic turbulence, but it is far weaker in one-dimensional spectra that are normally studied in laboratory turbulence. Structure function exponents for other orders are well described by the She-Leveque formula, but the velocity field is significantly less intermittent and the magnetic field is more intermittent than the Elsasser variables.

Nonhelical hydromagnetic turbulence without an imposed magnetic field is considered in the case where the magnetic Prandtl number is unity. The magnetic field is entirely due to dynamo action. The magnetic energy spectrum peaks at a wavenumber of about 5 times the minimum wavenumber in the domain, and not at the resistive scale, as has previously been argued. Throughout the inertial range the spectral magnetic energy exceeds the kinetic energy by a factor of about 2.5, and both spectra are approximately parallel. At first glance, the total energy spectrum seems to be close to k^{-3/2}, but there is a strong bottleneck effect and it is suggested that the asymptotic spectrum is k^{-5/3}. This is supported by the value of the second order structure function exponent that is found to be \zeta_2=0.70, suggesting a k^{-1.70} spectrum.

The kinematic dynamo problem is investigated for the flow of a conducting fluid in a cylindrical, periodic tube with conducting walls. The methods used are an eigenvalue analysis of the steady regime, and the three-dimensional solution of the time-dependent induction equation. The configuration and parameters considered here are close to those of a dynamo experiment planned in Perm, which will use a torus-shaped channel. We find growth of an initial magnetic field by more than 3 orders of magnitude. Marked field growth can be obtained if the braking time is less than 0.2 s and only one diverter is used in the channel. The structure of the seed field has a strong impact on the field amplification factor. The generation properties can be improved by adding ferromagnetic particles to the fluid in order to increase its relative permeability,but this will not be necessary for the success of the dynamo experiment. For higher magnetic Reynolds numbers, the nontrivial evolution of different magnetic modes limits the value of simple `optimistic' and `pessimistic' estimates.

At numerical resolutions around $512^3$ and above, three-dimensional energy spectra from turbulence simulations begin to show noticeably shallower spectra than $k^{-5/3}$ near the dissipation wavenumber (`bottleneck effect'). This effect is shown to be significantly weaker in one-dimensional spectra such as those obtained in wind tunnel turbulence. The difference can be understood in terms of the transformation between one-dimensional and three-dimensional energy spectra under the assumption that the turbulent velocity field is isotropic. Transversal and longitudinal energy spectra are similar and can both accurately be computed from the full three-dimensional spectra. Second-order structure functions are less susceptible to the bottleneck effect and may be better suited for inferring the scaling exponent from numerical simulation data.

An axisymmetric model of a cool, dynamo-active accretion disc is applied to protostellar discs. Thermally and magnetically driven outflows develop that are not collimated within 0.1 AU. In the presence of a central magnetic field from the protostar, accretion onto the protostar is highly episodic, which is in agreement with earlier work.

According to the kinematic theory of nonhelical dynamo action the magnetic energy spectrum increases with wavenumber and peaks at the resistive cutoff wavenumber. It has previously been argued that even in the dynamical case the magnetic energy peaks at the resistive scale. Using high resolution simulations (up to 1024^3 meshpoints) with no large scale imposed field we show that the magnetic energy peaks at a wavenumber that is independent of the magnetic Reynolds number and about 5 times larger than the forcing wavenumber. Throughout the inertial range the spectral magnetic energy exceeds the kinetic energy by a factor of 2 to 3. Both spectra are approximately parallel. The total energy spectrum seems to be close to k^{-3/2}, but there is a strong bottleneck effect and it is suggested that the asymptotic spectrum is instead k^{-5/3}. This is supported by the value of the second order structure function exponent that is found to be zeta_2=0.70, suggesting a k^{-1.70} spectrum. The third order structure function scaling exponent is very close to unity, in agreement with Goldreich-Sridhar theory. Adding an imposed field tends to suppress the small scale magnetic field. We find that at large scales the magnetic power spectrum follows then a k^{-1} slope. When the strength of the imposed field is of the same order as the dynamo generated field, we find almost equipartition between the magnetic and kinetic power spectra.

Isotropic homogeneous hydromagnetic turbulence is studied using numerical simulations at resolutions of up to 1024^3 meshpoints. It is argued that, in contrast to the kinematic regime, the nonlinear regime is characterized by a spectral magnetic power that is decreasing with increasing wavenumber, regardless of whether or not the turbulence has helicity. This means that the root-mean-square field strength converges to a limiting value at large magnetic Reynolds numbers. The total (magnetic and kinetic) energy spectrum tends to be somewhat shallower than k^{-5/3}, in agreement with the findings of other groups. In the presence of helicity, an inverse cascade develops, provided the scale separation between the size of the computational box and the scale of the energy carrying eddies exceeds a ratio of at least two. Finally, the constraints imposed by magnetic helicity conservation on mean-field theory are reviewed and new results of simulations are presented.

We present an axisymmetric numerical model of a dynamo active accretion disc. If the dynamo-generated magnetic field in the disc is sufficiently strong (close to equipartition with thermal energy), a fast magneto-centrifugally driven outflow develops within a conical shell near the rotation axis, together with a slower pressure driven outflow from the outer parts of the disc as well as around the axis. Our results show that a dynamo active accretion disc can contribute to driving an outflow even without any external magnetic field. The fast outflow in the conical shell is confined by the azimuthal field which is produced by the dynamo in the disc and advected to the disc corona. This part of the outflow has high angular momentum and is cooler and less dense than its surroundings. The conical shell's half-opening angle is typically about 30 degrees near the disc and decreases slightly with height. The slow outflow is hotter and denser.

Recent progress in the theory of solar and stellar dynamos is reviewed. Particular emphasis is placed on the mean-field theory which tries to describe the collective behavior of the magnetic field. In order to understand solar and stellar activity, a quantitatively reliable theory is necessary. Much of the new developments center around magnetic helicity conservation which is seen to be important in numerical simulations. Only a dynamical, explicitly time dependent theory of alpha-quenching is able to describe this behavior correctly.

The theory of large scale dynamos is reviewed with particular emphasis on the magnetic helicity constraint in the presence of closed and open boundaries. In the presence of closed or periodic boundaries, helical dynamos respond to the helicity constraint by developing small scale separation in the kinematic regime, and by showing long time scales in the nonlinear regime where the scale separation has grown to the maximum possible value. A resistively limited evolution towards saturation is also found at intermediate scales before the largest scale of the system is reached. Larger aspect ratios can give rise to different structures of the mean field which are obtained at early times, but the final saturation field strength is still decreasing with decreasing resistivity. In the presence of shear, cyclic magnetic fields are found whose period is increasing with decreasing resistivity, but the saturation energy of the mean field is in strong super-equipartition with the turbulent energy. It is shown that artificially induced losses of small scale field of opposite sign of magnetic helicity as the large scale field can, at least in principle, accelerate the production of large scale (poloidal) field. Based on mean field models with an outer potential field boundary condition in spherical geometry, we verify that the sign of the magnetic helicity flux from the large scale field agrees with the sign of alpha. For solar parameters, typical magnetic helicity fluxes lie around 10^{47} Mx^2 per cycle.

The self-excitation of magnetic field by a spiral Couette flow between two coaxial cylinders is considered. We solve numerically the fully nonlinear, three-dimensional MHD equations for magnetic Prandtl numbers Pm (ratio of kinematic viscosity to magnetic diffusivity) between 0.14 and 10 and kinematic and magnetic Reynolds numbers up to about 2000. In the initial stage of exponential field growth (kinematic dynamo regime), we find that the dynamo switches from one distinct regime to another as the radial width \delta r_B$ of the magnetic field distribution becomes smaller than the separation of the field maximum from the flow boundary. The saturation of magnetic field growth is due to a reduction in the velocity shear resulting mainly from the azimuthally and axially averaged part of the Lorentz force, which agrees with an asymptotic result for the limit of Pm<<1. In the parameter regime considered, the magnetic energy decreases with kinematic Reynolds number as Re^{-0.84}, which is approximately as predicted by the nonlinear asymptotic theory (\sim Re^{-1}). However, when the velocity field is maintained by a volume force (rather than by viscous stress) the dependence of magnetic energy on the kinematic Reynolds number is much weaker.

In ordinary turbulence research it has been a long standing tradition to solve the equations in spectral space giving the best possible accuracy. This is indeed a natural choice for incompressible problems with periodic boundaries, but it is no longer optimal in many astrophysical circumstances. It is argued that lower order spatial derivatives schemes are unacceptable in view of their low overall accuracy, even when mass, momentum, and energy are conserved to machine accuracy. High order finite difference schemes are therefore found to be quite efficient and physically appropriate. They are also easily and efficiently implemented on massively parallel computers. High order schemes also yield sufficient overall accuracy. Our code uses centered finite differences which make the adaptation to other problems simple. Since the code is not written in conservative form, conservation of mass, energy and momentum can be used to monitor to quality of the solution. A third order Runge-Kutta scheme with 2N-storage is used for calculating the time advance.

We consider the effect of vertical outflows on the mean-field dynamo in a thin disk. These outflows could be due to winds or magnetic buoyancy. We analyse both two-dimensional finite-difference numerical solutions of the axisymmetric dynamo equations and a free-decay mode expansion using the thin-disk approximation. Contrary to expectations, a vertical velocity can enhance dynamo action, provided it is not too strong. In the nonlinear regime this can lead to super-exponential growth of the magnetic field.

Dynamo action is investigated in simulations of locally isotropic and homogeneous turbulence in a slab between open boundaries. It is found that a `pseudo-vacuum' boundary condition (where the field is vertical) leads to strong helicity fluxes which significantly reduce the amplitude of the resulting large scale field. On the other hand, if there is a conducting halo outside the dynamo-active region the large scale field amplitude can reach larger values, but the time scale after which this field is reached increases linearly with the magnetic Reynolds number. In both cases most of the helicity flux is found to occur on large scales. From the variety of models considered we conclude that open boundaries tend to lower the saturation field strength compared to the case with periodic boundaries. The rate at which this lower saturation field strength is attained is roughly independent of the strength of the turbulence and of the boundary conditions. For dynamos with less helicity, however, significant field strengths could be reached in a shorter time.

We study the effects of vertical outflows on mean-field dynamos in disks. These outflows could be due to thermal winds or magnetic buoyancy. We analyse numerical solutions of the nonlinear mean-field dynamo equations using a two-dimensional finite-difference model. Contrary to expectations, a modest vertical velocity can enhance dynamo action. This can lead to super-exponential growth of the magnetic field and to higher magnetic energies at saturation in the nonlinear regime.