Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between $10^{32}$-$10^{34}$ erg~s$^{-1}$, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range $\sim10^{10}$-$10^{11}$ G. The decay of their weak dipole fields, mediated by the Hall effect and Ohmic dissipation, cannot provide sufficient thermal energy to power their X-ray luminosity, as opposed to magnetars whose X-ray luminosities are comparable. Motivated by the question of producing high X-ray power through magnetic field decay while maintaining a weak dipole field, we explore the evolution of a crustal magnetic field that does not consist of an ordered axisymmetric structure, but rather comprises a tangled configuration. This can be the outcome of a non-self-excited dynamo, buried inside the crust by fallback material following the supernova explosion. We find that such initial conditions lead to the emergence of the magnetic field from the surface of the star and the formation of a dipolar magnetic field component. An internal tangled magnetic field of the order of $10^{14}$ G can provide sufficient Ohmic heating to the crust and power CCOs, while the dipole field it forms is approximately $10^{10}$ G, as observed in CCOs.

Recent observations of Jupiter and Saturn suggest that heavy elements may be diluted in the gaseous envelope, providing a compositional gradient that could stabilise ordinary convection and produce a stably-stratified layer near the core of these planets. This region could consist of semi-convective layers with a staircase-like density profile, which have multiple convective zones separated by thin stably-stratified interfaces, as a result of double-diffusive convection. These layers could have important effects on wave propagation and tidal dissipation that have not been fully explored. We analyse the effects of these layers on the propagation and transmission of internal waves within giant planets, extending prior work in a local Cartesian model. We adopt a simplified global Boussinesq planetary model in which we explore the internal waves in a non-rotating spherical body. We begin by studying the free modes of a region containing semi-convective layers. We then analyse the transmission of internal waves through such a region. The free modes depend strongly on the staircase properties, and consist of modes with both internal and interfacial gravity wave-like behaviour. We determine the frequency shifts of these waves as a function of the number of steps to explore their potential to probe planetary internal structures. We also find that wave transmission is strongly affected by the presence of a staircase. Very large-wavelength waves are transmitted efficiently, but small-scale waves are only transmitted if they are resonant with one of the free modes. The effective size of the core is therefore larger for non-resonant modes.

We reveal and investigate a new type of linear axisymmetric helical magnetorotational instability which is capable of destabilizing viscous and resistive rotational flows with radially increasing angular velocity, or positive shear. This instability is double-diffusive by nature and is different from the more familiar helical magnetorotational instability, operating at positive shear above the Liu limit, in that it works instead for a wide range of the positive shear when ${\rm (i)}$ a combination of axial/poloidal and azimuthal/toroidal magnetic fields is applied and ${\rm (ii)}$ the magnetic Prandtl number is not too close to unity. We study this instability first with radially local WKB analysis and then confirm its existence using a global stability analysis of the magnetized flow between two rotating cylinders with conducting or insulating boundaries. From an experimental point of view, we also demonstrate the presence of the new instability in a magnetized viscous and resistive Taylor-Couette flow with positive shear for such values of the flow parameters, which can be realized in upcoming experiments at the DRESDYN facility. Finally, this instability might have implications for the dynamics of the equatorial parts of the solar tachocline and dynamo action there, since the above two necessary conditions for the instability to take place are satisfied in this region. Our global stability calculations for the tachocline-like configuration, representing a thin rotating cylindrical layer with the appropriate boundary conditions -- conducting inner and insulating outer cylinders -- and the values of the flow parameters, indicate that it can indeed arise in this case with a characteristic growth time comparable to the solar cycle period.

Decades ago S. Lundquist, S. Chandrasekhar, P.H. Roberts and R. J.~Tayler first posed questions about the stability of Taylor-Couette flows of conducting material under the influence of large-scale magnetic fields. These and many new questions can now be answered numerically where the nonlinear simulations even provide the instability-induced values of several transport coefficients. The cylindrical containers are axially unbounded and penetrated by magnetic background fields with axial and/or azimuthal components. The influence of the magnetic Prandtl number $Pm$ on the onset of the instabilities is shown to be substantial. The potential flow subject to axial fields becomes unstable against axisymmetric perturbations for a certain supercritical value of the averaged Reynolds number $\overline{Rm}=\sqrt{Re\cdot Rm}$ (with $Re$ the Reynolds number of rotation, $Rm$ its magnetic Reynolds number). Rotation profiles as flat as the quasi-Keplerian rotation law scale similarly but only for $Pm\gg 1$ while for $Pm\ll 1$ the instability instead sets in for supercritical $Rm$ at an optimal value of the magnetic field. Among the considered instabilities of azimuthal fields, those of the Chandrasekhar-type, where the background field and the background flow have identical radial profiles, are particularly interesting. They are unstable against nonaxisymmetric perturbations if at least one of the diffusivities is non-zero. For $Pm\ll 1$ the onset of the instability scales with $Re$ while it scales with $\overline{Rm}$ for $Pm\gg 1$. - Even superrotation can be destabilized by azimuthal and current-free magnetic fields; this recently discovered nonaxisymmetric instability is of a double-diffusive character, thus excluding $Pm= 1$. It scales with $Re$ for $Pm\to 0$ and with $Rm$ for $Pm\to \infty$.

Stellar magnetism plays an important role in stellar evolution theory. Approximatively 10% of observed main sequence (MS) and pre-main-sequence (PMS) radiative stars exhibit surface magnetic fields above the detection limit, raising the question of their origin. These stars host outer radiative envelopes, which are stably stratified. Therefore, they are assumed to be motionless in standard models of stellar structure and evolution. We focus on rapidly rotating, radiative stars which may be prone to the tidal instability, due to an orbital companion. Using direct numerical simulations in a sphere, we study the interplay between a stable stratification and the tidal instability, and assess its dynamo capability. We show that the tidal instability is triggered regardless of the strength of the stratification (Brunt-V{\"a}is{\"a}l{\"a} frequency). Furthermore, the tidal instability can lead to both mixing and self-induced magnetic fields in stably stratified layers (provided that the Brunt-V{\"a}is{\"a}l{\"a} frequency does not exceed the stellar spin rate in the simulations too much). The application to stars suggests that the resulting magnetic fields could be observable at the stellar surfaces. Indeed, we expect magnetic field strengths up to several Gauss.Consequently, tidally driven dynamos should be considered as a (complementary) dynamo mechanism, possibly operating in radiative MS and PMS stars hosting orbital companions. In particular, tidally driven dynamos may explain the observed magnetism of tidally deformed and rapidly rotating Vega-like stars.

We explore magnetic field configurations that lead to the formation of magnetic spots on the surface of neutron stars, and to the displacement of the magnetic dipole axis. We find that a toroidally dominated magnetic field is essential for the generation of a single spot with a strong magnetic field. Once a spot forms, it survives for several million years, even after the total magnetic field has decayed significantly. We find that the dipole axis is not stationary with respect to the neutron star's surface and does not in general coincide with the location of the magnetic spot. This is due to non-axisymmetric instabilities of the toroidal field that displace the poloidal dipole axis at rates that may reach $0.4^{\circ}$ per century. A misaligned poloidal dipole axis with the toroidal field leads to more significant displacement of the dipole axis than the fully aligned case. Finally we discuss the evolution of neutron stars with such magnetic fields on the $P-\dot{P}$ diagram and the observational implications. We find that neutron stars spend a very short time before they cross the Death-Line of the $P-\dot{P}$ diagram, compared to their characteristic ages. Moreover, the maximum intensity of their surface magnetic field is substantially higher than the dipole component of the field. We argue that SGR 0418+5729 could be an example of this type of behaviour, having a weak dipole field, yet hosting a magnetic spot responsible for its magnetar behaviour. The evolution on the pulse profile and braking index of the Crab pulsar, which are attributed to an increase of its obliquity, are compatible with the anticipated drift of the magnetic axis.

The magnetorotational instability (MRI) is thought to be a powerful source of turbulence in Keplerian accretion disks. Motivated by recent laboratory experiments, we study the MRI driven by an azimuthal magnetic field in an electrically conducting fluid sheared between two concentric rotating cylinders. By adjusting the rotation rates of the cylinders, we approximate angular velocity profiles $\omega \propto r^{q}$. We perform direct numerical simulations of a steep profile close to the Rayleigh line $q \gtrsim -2 $ and a quasi-Keplerian profile $q \approx -3/2$ and cover wide ranges of Reynolds ($Re\le 4\cdot10^4$) and magnetic Prandtl numbers ($0\le Pm \le 1$). In the quasi-Keplerian case, the onset of instability depends on the magnetic Reynolds number, with $Rm_c \approx 50$, and angular momentum transport scales as $\sqrt{Pm} Re^2$ in the turbulent regime. The ratio of Maxwell to Reynolds stresses is set by $Rm$. At the onset of instability both stresses have similar magnitude, whereas the Reynolds stress vanishes or becomes even negative as $Rm$ increases. For the profile close to the Rayleigh line, the instability shares these properties as long as $Pm\gtrsim0.1$, but exhibits a markedly different character if $Pm\rightarrow 0$, where the onset of instability is governed by the Reynolds number, with $Re_c \approx 1250$, transport is via Reynolds stresses and scales as $Re^2$. At intermediate $Pm=0.01$ we observe a continuous transition from one regime to the other, with a crossover at $Rm=\mathcal{O}(100)$. Our results give a comprehensive picture of angular momentum transport of the MRI with an imposed azimuthal field.

Since their formation, stars slow down their rotation rates by the removal of angular momentum from their surfaces, e.g. via stellar winds. Despite the complexity of the processes involved, a traditional model, where the removal of angular momentum loss by magnetic fields is prescribed, has provided a useful framework to understand observational relations between stellar rotation and age and magnetic field strength. Here, a spindown model is proposed where loss of angular momentum by magnetic fields is evolved dynamically, instead of being kinematically prescribed. To this end, we evolve the stellar rotation and magnetic field simultaneously over stellar evolution time by extending our previous work on a dynamo model which incorporates the nonlinear feedback mechanisms on rotation and magnetic fields. Our extended model reproduces key observations and explains the presence of the two branches of (fast and slow rotating) stars which have different relations between rotation rate $\Omega$ vs. time (age), magnetic field strength $|B|$ vs. rotation rate, and frequency of magnetic field $\omega_{cyc}$ vs. rotation rate. For fast rotating stars we find: (i) an exponential spindown $\Omega \propto e^{-1.35t}$, with $t$ measured in Gyrs, (ii) $|B|$ saturates for higher rotation rate, (iii) $\omega_{cyc} \propto \Omega^{0.85}$. For slow rotating stars we obtain: (i) a power law spindown $\Omega \propto t^{-0.52}$, (ii) $|B|$ scales almost linearly with rotation rate, (iii) $\omega_{cyc} \propto \Omega^{1.16}$. The results obtained are in good agreement with observations. The Vaughan-Preston gap is consistently explained in our model by the shortest spindown timescale in this transition from fast to slow rotators. Our results highlight the importance of self-regulation of magnetic fields and rotation by direct and indirect interactions involving nonlinear feedback in stellar evolution.

We numerically solve the magnetic induction equation in a spherical shell geometry, with a kinematically prescribed axisymmetric flow that consists of a superposition of a small-scale helical flow and a large-scale shear flow. The small-scale flow is chosen to be a local analog of the classical Roberts cells, consisting of strongly helical vortex rolls. The large-scale flow is a shearing motion in either the radial or the latitudinal directions. In the absence of large-scale shear, the small-scale flow is an efficient dynamo, in agreement with previous results. Adding increasingly large shear flows strongly suppresses the dynamo efficiency, indicating that shear is not always a favourable ingredient in dynamo action.

We study a resistive tearing instability developing in a system evolving through the combined effect of Hall drift in the Electron-MHD limit and Ohmic dissipation. We explore first the exponential growth of the instability in the linear case and we find the fastest growing mode, the corresponding eigenvalues and dispersion relation. The instability growth rate scales as $\gamma \propto B^{2/3} \sigma^{-1/3}$ where $B$ is the magnetic field and $\sigma$ the electrical conductivity. We confirm the development of the tearing resistive instability in the fully non-linear case, in a plane parallel configuration where the magnetic field polarity reverses, through simulations of systems initiating in Hall equilibrium with some superimposed perturbation. Following a transient phase, during which there is some minor rearrangement of the magnetic field, the perturbation grows exponentially. Once the instability is fully developed the magnetic field forms the characteristic islands and X-type reconnection points, where Ohmic decay is enhanced. We discuss the implications of this instability for the local magnetic field evolution in neutron stars' crusts, proposing that it can contribute to heating near the surface of the star, as suggested by models of magnetar post-burst cooling. In particular, we find that a current sheet a few meters thick, covering as little as $1\%$ of the total surface can provide $10^{42}~$erg in thermal energy within a few days. We briefly discuss applications of this instability in other systems where the Hall effect operates such as protoplanetary discs and space plasmas.

Current models of magnetars require extremely strong magnetic fields to explain their observed quiescent and bursting emission, implying that the field strength within the star's outer crust is orders of magnitude larger than the dipole component inferred from spin-down measurements. This presents a serious challenge to theories of magnetic field generation in a proto-neutron star. Here, we present detailed modelling of the evolution of the magnetic field in the crust of a neutron star through 3-D simulations. We find that, in the plausible scenario of equipartition of energy between global-scale poloidal and toroidal magnetic components, magnetic instabilities transfer energy to non-axisymmetric, kilometre-sized magnetic features, in which the local field strength can greatly exceed that of the global-scale field. These intense small-scale magnetic features can induce high energy bursts through local crust yielding, and the localised enhancement of Ohmic heating can power the star's persistent emission. Thus, the observed diversity in magnetar behaviour can be explained with mixed poloidal-toroidal fields of comparable energies.

We consider axially periodic Taylor-Couette geometry with insulating boundary conditions. The imposed basic states are so-called Chandrasekhar states, where the azimuthal flow $U_\phi$ and magnetic field $B_\phi$ have the same radial profiles. Mainly three particular profiles are considered: the Rayleigh limit, quasi-Keplerian, and solid-body rotation. In each case we begin by computing linear instability curves and their dependence on the magnetic Prandtl number Pm. For the azimuthal wavenumber m=1 modes, the instability curves always scale with the Reynolds number and the Hartmann number. For sufficiently small Pm these modes therefore only become unstable for magnetic Mach numbers less than unity, and are thus not relevant for most astrophysical applications. However, modes with m>10 can behave very differently. For sufficiently flat profiles, they scale with the magnetic Reynolds number and the Lundquist number, thereby allowing instability also for the large magnetic Mach numbers of astrophysical objects. We further compute fully nonlinear, three-dimensional equilibration of these instabilities, and investigate how the energy is distributed among the azimuthal (m) and axial (k) wavenumbers. In comparison spectra become steeper for large m, reflecting the smoothing action of shear. On the other hand kinetic and magnetic energy spectra exhibit similar behavior: if several azimuthal modes are already linearly unstable they are relatively flat, but for the rigidly rotating case where m=1 is the only unstable mode they are so steep that neither Kolmogorov nor Iroshnikov-Kraichnan spectra fit the results. The total magnetic energy exceeds the kinetic energy only for large magnetic Reynolds numbers Rm>100.

Motivated by recent advances in Direct Statistical Simulation (DSS) of astrophysical phenomena such as out of equilibrium jets, we perform a Direct Numerical Simulation (DNS) of the helical magnetorotational instability (HMRI) under the generalised quasilinear approximation (GQL). This approximation generalises the quasilinear approximation (QL) to include the self-consistent interaction of large-scale modes, interpolating between fully nonlinear DNS and QL DNS whilst still remaining formally linear in the small scales. In this paper we address whether GQL can more accurately describe low-order statistics of axisymmetric HMRI when compared with QL by performing DNS under various degrees of GQL approximation. We utilise various diagnostics, such as energy spectra in addition to first and second cumulants, for calculations performed for a range of Reynolds and Hartmann numbers (describing rotation and imposed magnetic field strength respectively). We find that GQL performs significantly better than QL in describing the statistics of the HMRI even when relatively few large-scale modes are kept in the formalism. We conclude that DSS based on GQL (GCE2) will be significantly more accurate than that based on QL (CE2).

We investigate the density-shear instability in Hall-MHD via numerical simulation of the full non-linear problem, in the context of magnetar activity. We confirm the development of the instability of a plane-parallel magnetic field with an appropriate intensity and electron density profile, in accordance with analytic theory. We find that the instability also appears for a monotonically decreasing electron number density and magnetic field, a plane-parallel analogue of an azimuthal or meridional magnetic field in the crust of a magnetar. The growth rate of the instability depends on the Hall properties of the field (magnetic field intensity, electron number density and the corresponding scale-heights), while being insensitive to weak resistivity. Since the Hall effect is the driving process for the evolution of the crustal magnetic field of magnetars, we argue that this instability is critical for systems containing strong meridional or azimuthal fields. We find that this process mediates the formation of localised structures with much stronger magnetic field than the average, which can lead to magnetar activity and accelerate the dissipation of the field and consequently the production of Ohmic heating. Assuming a $5\times10^{14}$G magnetic field at the base of crust, we anticipate that magnetic field as strong as $10^{15}$G will easily develop in regions of typical size of a few $10^{2}$ meters, containing magnetic energy of $10^{43}$erg, sufficient to power magnetar bursts. These active regions are more likely to appear in the magnetic equator where the tangential magnetic field is stronger.

We present the first fully self-consistent three-dimensional model of a neutron star's magnetic field, generated by electric currents in the star's crust via the Hall effect. We find that the global-scale field converges to a Hall-attractor state, as seen in recent axisymmetric models, but that small-scale features in the magnetic field survive even on much longer timescales. These small-scale features propagate toward the dipole equator, where the crustal electric currents organize themselves into a strong equatorial jet. By calculating the distribution of magnetic stresses in the crust, we predict that neutron stars with fields stronger than $10^{14}$G can still be subject to starquakes more than $10^5$yr after their formation.

We discuss a novel instability in inertia-less electron magneto-hydrodynamics (EMHD), which arises from a combination of electron velocity shear and electron density gradients. The unstable modes have a lengthscale longer than the transverse density scale, and a growth-rate of the order of the inverse Hall timescale. We suggest that this density-shear instability may be of importance in magnetic reconnection regions on scales smaller than the ion skin depth, and in neutron star crusts. We demonstrate that the so-called Hall drift instability, previously argued to be relevant in neutron star crusts, is a resistive tearing instability rather than an instability of the Hall term itself. We argue that the density-shear instability is of greater significance in neutron stars than the tearing instability, because it generally has a faster growth-rate and is less sensitive to geometry and boundary conditions. We prove that, for uniform electron density, EMHD is "at least as stable" as regular, incompressible MHD, in the sense that any field configuration that is stable in MHD is also stable in EMHD. We present a connection between the density-shear instability in EMHD and the magneto-buoyancy instability in anelastic MHD.

The azimuthal version of the magnetorotational instability (MRI) is a non-axisymmetric instability of a hydrodynamically stable differentially rotating flow under the influence of a purely or predominantly azimuthal magnetic field. It may be of considerable importance for destabilizing accretion disks, and plays a central role in the concept of the MRI dynamo. We report the results of a liquid metal Taylor-Couette experiment that shows the occurrence of an azimuthal MRI in the expected range of Hartmann numbers.

Large-scale planetary or stellar magnetic fields generated by a dynamo effect are mostly attributed to flows forced by buoyancy forces in electrically conducting fluid layers. However, these large-scale fields may also be controlled by tides, as previously suggested for the star $\tau$-boo, Mars or the Early Moon. By simulating a small local patch of a rotating fluid, \cite{Barker2014} have recently shown that tides can drive small-scale dynamos by exciting a hydrodynamic instability, the so-called elliptical (or tidal) instability. By performing global magnetohydrodynamic simulations of a rotating spherical fluid body, we investigate if this instability can also drive the observed large-scale magnetic fields. We are thus interested by the dynamo threshold and the generated magnetic field in order to test if such a mechanism is relevant for planets and stars. Rather than solving the problem in a geometry deformed by tides, we consider a spherical fluid body and add a body force to mimic the tidal deformation in the bulk of the fluid. This allows us to use an efficient spectral code to solve the magnetohydrodynamic problem. We first compare the hydrodynamic results with theoretical asymptotic results, and numerical results obtained in a truely deformed ellipsoid, which confirms the presence of the elliptical instability. We then perform magnetohydrodynamic simulations, and investigate the dynamo capability of the flow. Kinematic and self-consistent dynamos are finally simulated, showing that the elliptical instability is capable of generating dipole dominated large-scale magnetic field in global simulations of a fluid rotating sphere.

So far, numerical studies of double-diffusive layering in turbulent convective flows have neglected the effects of rotation. We undertake a first step into that direction by investigating how Coriolis forces affect a double-diffusive layer inside a rotating spherical shell. For this purpose we have run simulations in a parameter regime where these layers are expected to form and successively increased the rate of rotation with the result that fast rotation is found to have a similar stabilising effect on the overall convective flux as an increase of the stability ratio $R_{\rho}$ has in a non-rotating setup. We have also studied to what extent the regimes of rotational constraints suggested by King, Stellmach, and Buffett (2013) for rotation in the case of Rayleigh-B\'enard convection are influenced by double-diffusive convection: their classification could also be applicable to the case of double-diffusive convection in a spherical shell if it is extended to be also a function of the stability ratio $R_{\rho}$. Furthermore, we examined the ratio of saline and thermal Nusselt numbers and compared our results with models of Spruit (2013), Rosenblum et al. (2011) and Wood, Garaud, and Stellmach (2013). We find our data to be fitted best by Spruit's model. Our result that fast rotation further decreases the convective transport, which is already lowered by double-diffusive convection, could play a major role for e.g. the modeling of the interior of some rapidly rotating giant planets, as gaseous giant planets have recently been proposed to be influenced by double-diffusive convection.

The interaction of differential rotation and toroidal fields that are current-free in the gap between two corotating axially unbounded cylinders is considered. It is shown that nonaxisymmetric perturbations are unstable if the rotation rate and Alfv\'en frequency of the field are of the same order, almost independent of the magnetic Prandtl number Pm. For the very steep rotation law \Omega\propto R^{-2} (the Rayleigh limit) and for small Pm the threshold values of rotation and field for this Azimuthal MagnetoRotational Instability (AMRI) scale with the ordinary Reynolds number and the Hartmann number, resp. A laboratory experiment with liquid metals like sodium or gallium in a Taylor-Couette container has been designed on the basis of this finding. For fluids with more flat rotation laws the Reynolds number and the Hartmann number are no longer typical quantities for the instability. For the weakly nonlinear system the numerical values of the kinetic energy and the magnetic energy are derived for magnetic Prandtl numbers \leq 1. We find that the magnetic energy grows monotonically with the magnetic Reynolds number Rm, while the kinetic energy grows with Rm/\sqrt{Pm}. The resulting turbulent Schmidt number, as the ratio of the `eddy' viscosity and the diffusion coefficient of a passive scalar (such as lithium) is of order 20 for Pm=1, but for small Pm it drops to order unity. Hence, in a stellar core with fossil fields and steep rotation law the transport of angular momentum by AMRI is always accompanied by an intense mixing of the plasma, until the rotation becomes rigid.

We investigate numerically kinematic dynamos driven by flow of electrically conducting fluid in the shell between two concentric differentially rotating spheres, a configuration normally referred to as spherical Couette flow. We compare between axisymmetric (2D) and fully three dimensional flows, between low and high global rotation rates, between prograde and retrograde differential rotations, between weak and strong nonlinear inertial forces, between insulating and conducting boundaries, and between two aspect ratios. The main results are as follows. Azimuthally drifting Rossby waves arising from the destabilisation of the Stewartson shear layer are crucial to dynamo action. Differential rotation and helical Rossby waves combine to contribute to the spherical Couette dynamo. At a slow global rotation rate, the direction of differential rotation plays an important role in the dynamo because of different patterns of Rossby waves in prograde and retrograde flows. At a rapid global rotation rate, stronger flow supercriticality (namely the difference between the differential rotation rate of the flow and its critical value for the onset of nonaxisymmetric instability) facilitates the onset of dynamo action. A conducting magnetic boundary condition and a larger aspect ratio both favour dynamo action.

The onset of convection in a rapidly rotating layer in which a thermal wind is present is studied. Diffusive effects are included. The main motivation is from convection in planetary interiors, where thermal winds are expected due to temperature variations on the core-mantle boundary. The system admits both convective instability and baroclinic instability. We find a smooth transition between the two types of modes, and investigate where the transition region between the two types of instability occurs in parameter space. The thermal wind helps to destabilise the convective modes. Baroclinic instability can occur when the applied vertical temperature gradient is stable, and the critical Rayleigh number is then negative. Long wavelength modes are the first to become unstable. Asymptotic analysis is possible for the transition region and also for long wavelength instabilities, and the results agree well with our numerical solutions. We also investigate how the instabilities in this system relate to the classical baroclinic instability in the Eady problem. We conclude by noting that baroclinic instabilities in the Earth's core arising from heterogeneity in the lower mantle could possibly drive a dynamo even if the Earth's core were stably stratified and so not convecting.

Aims. In the first few minutes of a newborn neutron star's life the Hall effect and differential rotation may both be important. We demonstrate that these two ingredients are sufficient for generating a 'shear-Hall instability' and for studying its excitation conditions, growth rates, and characteristic magnetic field patterns. Methods. We numerically solve the induction equation in a spherical shell, with a kinematically prescribed differential rotation profile {\Omega}(s), where s is the cylindrical radius. The Hall term is linearized about an imposed uniform axial field. The linear stability of individual azimuthal modes, both axisymmetric and non-axisymmetric, is then investigated. Results. For the shear-Hall instability to occur, the axial field must be parallel to the rotation axis if {\Omega}(s) decreases outward, whereas if {\Omega}(s) increases outward it must be anti-parallel. The instability draws its energy from the differential rotation, and occurs on the short rotational timescale rather than on the much longer Hall timescale. It operates most efficiently if the Hall time is comparable to the diffusion time. Depending on the precise field strengths B0, either axisymmetric or non-axisymmetric modes may be the most unstable. Conclusions. Even if the differential rotation in newborn neutron stars is quenched within minutes, the shear-Hall instability may nevertheless amplify any seed magnetic fields by many orders of magnitude.

Helical magnetic background fields with adjustable pitch angle are imposed on a conducting fluid in a differentially rotating cylindrical container. The small-scale kinetic and current helicities are calculated for various field geometries, and shown to have the opposite sign as the helicity of the large-scale field. These helicities and also the corresponding $\alpha$-effect scale with the current helicity of the background field. The $\alpha$-tensor is highly anisotropic as the components $\alpha_{\phi\phi}$ and $\alpha_{zz}$ have opposite signs. The amplitudes of the azimuthal $\alpha$-effect computed with the cylindrical 3D MHD code are so small that the operation of an $\alpha\Omega$ dynamo on the basis of the current-driven, kink-type instabilities of toroidal fields is highly questionable. In any case the low value of the $\alpha$-effect would lead to very long growth times of a dynamo in the radiation zone of the Sun and early-type stars of the order of mega-years.

We explore numerically the flow induced in a spherical shell by differentially rotating the inner and outer spheres. The fluid is also taken to be electrically conducting (in the low magnetic Reynolds number limit), and a magnetic field is imposed parallel to the axis of rotation. If the outer sphere is stationary, the magnetic field induces a Shercliffe layer on the tangent cylinder, the cylinder just touching the inner sphere and parallel to the field. If the magnetic field is absent, but a strong overall rotation is present, Coriolis effects induce a Stewartson layer on the tangent cylinder. The non-axisymmetric instabilities of both types of layer separately have been studied before; here we consider the two cases side by side, as well as the mixed case, and investigate how magnetic and rotational effects interact. We find that if the differential rotation and the overall rotation are in the same direction, the overall rotation may have a destabilizing influence, whereas if the differential rotation and the overall rotation are in the opposite direction, the overall rotation always has a stabilizing influence.

We present axisymmetric numerical calculations of the fluid flow induced in a spherical shell with inner sphere rotating and outer sphere stationary. A magnetic field is also imposed, consisting of particular linear combinations of axial and dipolar fields, chosen to make $B_r=0$ at either the outer sphere, or the inner, or in between. This leads to the formation of Shercliff shear layers at these particular locations. We then consider the effect of increasingly large inertial effects, and show that an outer Shercliff layer is eventually de-stabilized, an inner Shercliff layer appears to remain stable, and an in-between Shercliff layer is almost completely disrupted even before the onset of time-dependence, which does eventually occur though.

The linear stability of MHD Taylor-Couette flows in axially unbounded cylinders is considered, for magnetic Prandtl number unity. Magnetic fields varying from purely axial to purely azimuthal are imposed, with a general helical field parameterized by \beta=B_\phi/B_z. We map out the transition from the standard MRI for \beta=0 to the nonaxisymmetric Azimuthal MagnetoRotational Instability (AMRI) for \beta\to \infty. For finite \beta, positive and negative wave numbers m, corresponding to right and left spirals, are no longer identical. The transition from \beta=0 to \beta\to\infty includes all the possible forms of MRI with axisymmetric and nonaxisymmetric modes. For the nonaxisymmetric modes, the most unstable mode spirals in the opposite direction to the background field. The standard (\beta=0) MRI is axisymmetric for weak fields (including the instability with the lowest Reynolds number) but is nonaxisymmetric for stronger fields. If the azimuthal field is due in part to an axial current flowing through the fluid itself (and not just along the central axis), then it is also unstable to the nonaxisymmetric Tayler instability, which is most effective without rotation. For large \beta this instability has wavenumber m=1, whereas for \beta\simeq 1 m=2 is most unstable. The most unstable mode spirals in the same direction as the background field.

We numerically study the evolution of magnetic fields and fluid flows in a thin spherical shell. We take the initial field to be a latitudinally confined, predominantly toroidal flux tube. For purely toroidal, untwisted flux tubes, we recover previously known radial-shredding instabilities, and show further that in the nonlinear regime these instabilities can very effectively destroy the original field. For twisted flux tubes, including also a poloidal component, there are several possibilities, including the suppression of the radial-shredding instability, but also a more directly induced evolution, brought about because twisted flux tubes in general are not equilibrium solutions of the governing equations.

The Hall effect is an important nonlinear mechanism affecting the evolution of magnetic fields in neutron stars. Studies of the governing equation, both theoretical and numerical, have shown that the Hall effect proceeds in a turbulent cascade of energy from large to small scales. We investigate the small-scale Hall instability conjectured to exist from the linear stability analysis of Rheinhardt and Geppert. Identical linear stability analyses are performed to find a suitable background field to model Rheinhardt and Geppert's ideas. The nonlinear evolution of this field is then modelled using a three-dimensional pseudospectral numerical MHD code. Combined with the background field, energy was injected at the ten specific eigenmodes with the greatest positive eigenvalues as inferred by the linear stability analysis. Energy is transferred to different scales in the system, but not into small scales to any extent that could be interpreted as a Hall instability. Any instabilities are overwhelmed by a late-onset turbulent Hall cascade, initially avoided by the choice of background field, but soon generated by nonlinear interactions between the growing eigenmodes. The Hall cascade is shown here, and by several authors elsewhere, to be the dominant mechanism in this system.

The magnetorotational instability (MRI) is thought to play a key role in the formation of stars and black holes by sustaining the turbulence in hydrodynamically stable Keplerian accretion discs. In previous experiments the MRI was observed in a liquid metal Taylor-Couette flow at moderate Reynolds numbers by applying a helical magnetic field. The observation of this helical MRI (HMRI) was interfered with a significant Ekman pumping driven by solid end-caps that confined the instability only to a part of the Taylor-Couette cell. This paper describes the observation of the HMRI in an improved Taylor-Couette setup with the Ekman pumping significantly reduced by using split end-caps. The HMRI, which now spreads over the whole height of the cell, appears much sharper and in better agreement with numerical predictions. By analyzing various parameter dependencies we conclude that the observed HMRI represents a self-sustained global instability rather than a noise-sustained convective one.

The magnetorotational instability (MRI) plays an essential role in the formation of stars and black holes. By destabilizing hydrodynamically stable Keplerian flows, the MRI triggers turbulence and enables outward transport of angular momentum in accretion discs. We present the results of a liquid metal Taylor-Couette experiment under the influence of helical magnetic fields that show typical features of MRI at Reynolds numbers of the order 1000 and Hartmann numbers of the order 10. Particular focus is laid on an improved experiment in which split end caps are used to minimize the Ekman pumping.

Liu et al. [Phys. Rev. E 74, 056302 (2006)] have presented a WKB analysis of the helical magnetorotational instability (HMRI), and claim that it does not exist for Keplerian rotation profiles. We show that if radial boundary conditions are included, the HMRI can exist even for rotation profiles as flat as Keplerian, provided only that at least one of the boundaries is sufficiently conducting.

We consider the effect of toroidal magnetic fields on hydrodynamically stable Taylor-Couette differential rotation flows. For current-free magnetic fields a nonaxisymmetric m=1 magnetorotational instability arises when the magnetic Reynolds number exceeds O (100). We then consider how this `azimuthal magnetorotational instability' (AMRI) is modified if the magnetic field is not current-free, but also has an associated electric current throughout the fluid. This gives rise to current-driven Tayler instabilities (TI) that exist even without any differential rotation at all. The interaction of the AMRI and the TI is then considered when both electric currents and differential rotation are present simultaneously. The magnetic Prandtl number Pm turns out to be crucial in this case. Large Pm have a destabilizing influence, and lead to a smooth transition between the AMRI and the TI. In contrast, small Pm have a stabilizing influence, with a broad stable zone separating the AMRI and the TI. In this region the differential rotation is acting to stabilize the Tayler instabilities, with possible astrophysical applications (Ap stars). The growth rates of both the AMRI and the TI are largely independent of Pm, with the TI acting on the timescale of a single rotation period, and the AMRI slightly slower, but still on the basic rotational timescale. The azimuthal drift timescale is ~20 rotations, and may thus be a (flip-flop) timescale of stellar activity between the rotation period and the diffusion time.

The magnetorotational instability (MRI) plays a key role in the formation of stars and black holes, by enabling outward angular momentum transport in accretion disks. The use of combined axial and azimuthal magnetic fields allows the investigation of this effect in liquid metal flows at moderate Reynolds and Hartmann numbers. A variety of experimental results is presented showing evidence for the occurrence of the MRI in a Taylor-Couette flow using the liquid metal alloy GaInSn.

We consider the linear stability of dissipative MHD Taylor-Couette flow with imposed toroidal magnetic fields. The inner and outer cylinders can be either insulating or conducting; the inner one rotates, the outer one is stationary. The magnetic Prandtl number can be as small as 10-5, approaching realistic liquid-metal values. The magnetic field destabilizes the flow, except for radial profiles of B$_\phi$(R) close to the current-free solution. The profile with B$_{in}$=B$_{out}$ (the most uniform field) is considered in detail. For weak fields the TC-flow is stabilized, until for moderately strong fields the m=1 azimuthal mode dramatically destabilizes the flow again. There is thus a maximum value for the critical Reynolds number. For sufficiently strong fields (as measured by the Hartmann number) the toroidal field is always unstable, even for Re=0. The electric currents needed to generate the required toroidal fields in laboratory experiments are a few kA if liquid sodium is used, somewhat more if gallium is used. Weaker currents are needed for wider gaps, so a wide-gap apparatus could succeed even with gallium. The critical Reynolds numbers are only somewhat larger than the nonmagnetic values, so such an experiment would require only modest rotation rates.

We consider the flow of an electrically conducting fluid between differentially rotating cylinders, in the presence of an externally imposed toroidal field B_0 (r_i/r) e_phi. It is known that the classical, axisymmetric magnetorotational instability does not exist for such a purely toroidal imposed field. We show here that a non-axisymmetric magnetorotational instability does exist, having properties very similar to the axisymmetric magnetorotational instability in the presence of an axial field.

We study experimentally the flow of a liquid metal confined between differentially rotating cylinders, in the presence of externally imposed axial and azimuthal magnetic fields. For increasingly large azimuthal fields a wave-like disturbance arises, traveling along the axis of the cylinders. The wavelengths and speeds of these structures, as well as the field strengths and rotation rates at which they arise, are broadly consistent with theoretical predictions of such a traveling wave magnetorotational instability.

A recent paper [R. Hollerbach and G. Rudiger, Phys. Rev. Lett. 95, 124501 (2005)] has shown that the threshold for the onset of the magnetorotational instability (MRI) in a Taylor-Couette flow is dramatically reduced if both axial and azimuthal magnetic fields are imposed. In agreement with this prediction, we present results of a Taylor-Couette experiment with the liquid metal alloy GaInSn, showing evidence for the existence of the MRI at Reynolds numbers of order 1000 and Hartmann numbers of order 10.

We study the stability of cylindrical Taylor-Couette flow in the presence of combined axial and azimuthal magnetic fields, and show that adding an azimuthal field profoundly alters the previous results for purely axial fields. For small magnetic Prandtl numbers Pm, the critical Reynolds number Re_c for the onset of the magneto-rotational instability becomes independent of Pm, whereas for purely axial fields it scales as Pm^{-1}. For typical liquid metals, Re_c is then reduced by several orders of magnitude, enough that this new design should succeed in realizing this instability in the laboratory.

We study the magnetorotational instability in cylindrical Taylor-Couette flow, with the (vertically unbounded) cylinders taken to be perfect conductors, and with externally imposed spiral magnetic fields. The azimuthal component of this field is generated by an axial current inside the inner cylinder, and may be slightly stronger than the axial field. We obtain an instability beyond the Rayleigh line, for Reynolds numbers of order 1000 and Hartmann numbers of order 10, and independent of the (small) magnetic Prandtl number. For experiments with R_out = 2 R_in = 10 cm and Omega_out = 0.27 Omega_in, the instability appears for liquid sodium for axial fields of ~20 Gauss and axial currents of \~1200 A. For gallium the numbers are ~50 Gauss and ~3200 A. The vertical cell size is about twice the cell size known for nonmagnetic experiments.

We present numerical simulations of the flow in a rapidly rotating cylindrical annulus. We show that at the rotation rates relevant to the magneto-rotational instability, the flow is strongly constrained by the Taylor-Proudman theorem. As a result, it is controlled almost entirely by the end-plates. We then consider two possible options for minimizing these end-effects, namely (i) simply taking a very long cylinder, and (ii) splitting the end-plates into a series of differentially rotating rings. Regarding option (i), we show that the cylinder would have to be hundreds of times as long as it is wide before end-effects become unimportant in the interior. Since this is clearly not feasible, we turn to option (ii), and show that in order to obtain a smooth angular velocity profile, the end-plates would have to be split into around ten rings. If the end-plates are split into fewer rings, perhaps 3-5, the angular velocity profile will not be smooth, but will instead consist of a series of Stewartson layers at the boundaries from one ring to the next. We suggest therefore that the instabilities one obtains in this system will be the familiar Kelvin-Helmholtz instabilities of these Stewartson layers, rather than the magneto-rotational instability. At best, one might hope to obtain the MRI superimposed on these Kelvin-Helmholtz modes. Any subsequent interpretation of results is thus likely to be quite complicated.

We consider the evolution of magnetic fields under the influence of Hall drift and Ohmic decay. The governing equation is solved numerically, in a spherical shell with r_i/r_o = 0.75. Starting with simple free decay modes as initial conditions, we then consider the subsequent evolution. The Hall effect induces so-called helicoidal oscillations, in which energy is redistributed among the different modes. We find that the amplitude of these oscillations can be quite substantial, with some of the higher harmonics becoming comparable to the original field. Nevertheless, this transfer of energy to the higher harmonics is not sufficient to significantly accelerate the decay of the original field, at least not at the R_B = O(100) parameter values accessible to us, where this Hall parameter R_B measures the ratio of the Ohmic timescale to the Hall timescale. We do find clear evidence though of increasingly fine structures developing for increasingly large R_B, suggesting that perhaps this Hall-induced cascade to ever shorter lengthscales is eventually sufficiently vigorous to enhance the decay of the original field. Finally, the implications for the evolution of neutron star magnetic fields are discussed.

Stars of spectral classes A and late B are almost entirely radiative. CP stars are a slowly rotating subgroup of these stars. It is possible that they possessed long-lived accretion disks in their T Tauri phase. Magnetic coupling of disk and star leads to rotational braking at the surface of the star. Microscopic viscosities are extremely small and will not be able to reduce the rotation rate of the core of the star. We investigate the question whether magneto-rotational instability can provide turbulent angular momentum transport. We illuminate the question whether or not differential rotation is present in CP stars. Numerical MHD simulations of thick stellar shells are performed. An initial differential rotation law is subject to the influence of a magnetic field. The configuration gives indeed rise to magneto-rotational instability. The emerging flows and magnetic fields transport efficiently angular momentum outwards. Weak dependence on the magnetic Prandtl number (~0.01 in stars) is found from the simulations. Since the estimated time-scale of decay of differential rotation is 10^7-10^8 yr and comparable to the life-time of A stars, we find the braking of the core to be an ongoing process in many CP stars. The evolution of the surface rotation of CP stars with age will be an observational challenge and of much value for verifying the simulations.