Guided by the working hypothesis that the Schwabe cycle of solar activity is synchronized by the 11.07 years alignment cycle of the tidally dominant planets Venus, Earth and Jupiter, we reconsider the phase diagrams of sediment accumulation rates in Lake Holzmaar, and of methanesulfonate (MSA) data in the Greenland ice core GISP2, which are available for the period 10000-9000 cal. BP. Since some half-cycle phase jumps appearing in the output signals are, very likely, artifacts of applying a biologically substantiated transfer function, the underlying solar input signal with a dominant 11.04 years periodicity can be considered as mainly phase-coherent over the 1000 years period in the early Holocene. For more recent times, we show that the re-introduction of a hypothesized "lost cycle" at the beginning of the Dalton minimum would lead to a real phase jump. Similarly, by analyzing various series of $^{14}$C and $^{10}$Be data and comparing them with Schove's historical cycle maxima, we support the existence of another "lost cycle" around 1565, also connected with a real phase jump. Viewed synoptically, our results lend greater plausibility to the starting hypothesis of a tidally synchronized solar cycle, which at times can undergo phase jumps, although the competing explanation in terms of a non-linear solar dynamo with increased coherence cannot be completely ruled out.

We argue that the most prominent temporal features of the solar dynamo, in particular the Hale cycle, the Suess-de Vries cycle (associated with variations of the Gnevyshev-Ohl rule), Gleissberg-type cycles, and grand minima can be self-consistently explained by double synchronization with the 11.07-years periodic tidal forcing of the Venus-Earth-Jupiter system and the (mainly) 19.86-years periodic motion of the Sun around the barycenter of the solar system. In our numerical simulation, grand minima, and clusters thereof, emerge as intermittent and non-periodic events on millennial time scales, very similar to the series of Bond events which were observed throughout the Holocene and the last glacial period. If confirmed, such an intermittent transition to chaos would prevent any long-term prediction of solar activity, notwithstanding the fact that the shorter-term Hale and Suess-de Vries cycles are clocked by planetary motion.

Aiming at a consistent planetary synchronization model of both short-term and long-term solar cycles, we start with an analysis of Schove's historical data of cycle maxima. Their deviations (residuals) from the average cycle duration of 11.07 years show a high degree of regularity, comprising a dominant 200-year period (Suess-de Vries cycle), and a few periods around 100 years (Gleissberg cycle). Encouraged by their robustness, we support previous forecasts of an upcoming grand minimum in the 21st century. To explain the long-term cycles, we enhance our tidally synchronized solar dynamo model by a modulation of the field storage capacity of the tachocline with the orbital angular momentum of the Sun, which is dominated by the 19.86-year periodicity of the Jupiter-Saturn synodes. This modulation of the 22.14 years Hale cycle leads to a 193-year beat period of dynamo activity which is indeed close to the Suess-de Vries cycle. For stronger dynamo modulation, the model produces additional peaks at typical Gleissberg frequencies, which seem to be explainable by the non-linearities of the basic beat process, leading to a bi-modality of the Schwabe cycle. However, a complementary role of beat periods between the Schwabe cycle and the Jupiter-Uranus/Neptune synodic cycles cannot be completely excluded.

We report the results of a liquid metal Taylor-Couette experiment in the Rayleigh-stable regime under the influence of an azimuthal magnetic field. We observe that the resulting azimuthal magnetorotational instability (AMRI) from our experimental setup is significantly influenced by the thermal boundary conditions. Even a minimal radial heat flux leads to a symmetry breaking, which results in the AMRI waves traveling either upwards or downwards. We identify the thermal radiation by the central axial current as the heat source responsible for vertical convection in the liquid. Preliminary numerical investigations point towards an interaction between AMRI and thermal convection, which supports our experimental findings.

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.

We discuss a solar dynamo model of Tayler-Spruit type whose Omega-effect is conventionally produced by a solar-like differential rotation but whose alpha-effect is assumed to be periodically modulated by planetary tidal forcing. This resonance-like effect has its rationale in the tendency of the current-driven Tayler instability to undergo intrinsic helicity oscillations which, in turn, can be synchronized by periodic tidal perturbations. Specifically, we focus on the 11.07 years alignment periodicity of the tidally dominant planets Venus, Earth, and Jupiter, whose persistent synchronization with the solar dynamo is briefly touched upon. The typically emerging dynamo modes are dipolar fields, oscillating with a 22.14 years period or pulsating with a 11.07 years period, but also quadrupolar fields with corresponding periodicities. In the absence of any constant part of alpha, we prove the subcritical nature of this Tayler-Spruit type dynamo. The resulting amplitude of the alpha oscillation that is required for dynamo action turns out to lie in the order of 1 m/s, which seems not implausible for the sun. When starting with a more classical, non-periodic part of alpha, even less of the oscillatory alpha part is needed to synchronize the entire dynamo. Typically, the dipole solutions show butterfly diagrams, although their shapes are not convincing yet. Phase coherent transitions between dipoles and quadrupoles, which are reminiscent of the observed behaviour during the Maunder minimum, can be easily triggered by long-term variations of dynamo parameters, but may also occur spontaneously even for fixed parameters. Further interesting features of the model are the typical second intensity peak and the intermittent appearance of reversed helicities in both hemispheres.

The project DRESDYN (DREsden Sodium facility for DYNamo and thermohydraulic studies) conducted at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) provides a new platform for a variety of liquid sodium experiments devoted to problems of geo- and astrophysical magnetohydrodynamics. The most ambitious experiment within this project is a precession driven dynamo experiment that currently is under construction. It consists of a cylinder filled with liquid sodium that simultaneously rotates around two axes. The experiment is motivated by the idea of a precession-driven flow as a complementary energy source for the geodynamo or the ancient lunar dynamo. In the present study we address numerical and experimental examinations in order to identify parameter regions where the onset of magnetic field excitation will be most probable. Both approaches show that in the strongly nonlinear regime the flow is essentially composed of the directly forced primary Kelvin mode and higher modes in terms of standing inertial waves that arise from nonlinear self-interactions. A peculiarity is the resonance-like emergence of an axisymmetric mode that represents a double roll structure in the meridional plane, which, however, only occurs in a very limited range of the precession ratio. This axisymmetric mode turns out to be beneficial for dynamo action, and kinematic simulations of the magnetic field evolution induced by the time-averaged flow exhibit magnetic field excitation at critical magnetic Reynolds numbers around ${\rm{Rm}}^{\rm{c}}\approx 430$, which is well within the range of the planned liquid sodium experiment.

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$.

We have conducted experimental measurements and numerical simulations of a precession driven flow in a cylindrical cavity. The study is dedicated to the precession dynamo experiment currently under construction at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) and aims at the evaluation of the hydrodynamic flow with respect to its ability to drive a dynamo. We focus on the strongly nonlinear regime in which the flow is essentially composed of the directly forced primary Kelvin mode and higher modes in terms of standing inertial waves arising from nonlinear self-interactions. We obtain an excellent agreement between experiment and simulation with regard to both, flow amplitudes and flow geometry. A peculiarity is the resonance-like emergence of an axisymmetric mode that represents a double roll structure in the meridional plane. Kinematic simulations of the magnetic field evolution induced by the time-averaged flow yield dynamo action at critical magnetic Reynolds numbers around ${\rm{Rm}}^{\rm{c}}\approx 430$ which is well within the range of the planned liquid sodium experiment.

It is demonstrated that the azimuthal magnetorotational instability (AMRI) also works with radially increasing rotation rates contrary to the standard magnetorotational instability for axial fields which requires negative shear. The stability against nonaxisymmetric perturbations of a conducting Taylor-Couette flow with positive shear under the influence of a toroidal magnetic field is considered if the background field between the cylinders is current-free. For small magnetic Prandtl number Pm-->0 the curves of neutral stability converge in the Hartmann number/Reynolds number plane approximating the stability curve obtained in the inductionless limit Pm=0. The numerical solutions for Pm=0 indicate the existence of a lower limit of the shear rate. For large Pm the curves scale with the magnetic Reynolds number of the outer cylinder but the flow is always stable for magnetic Prandtl number unity as is typical for double-diffusive instabilities. We are particularly interested to know the minimum Hartmann number for neutral stability. For models with resting or almost resting inner cylinder and with perfect-conducting cylinder material the minimum Hartmann number occurs for a radius ratio of 0.9. The corresponding critical Reynolds numbers are smaller than 10.000.

Magnetorotational instability (MRI) is one of the fundamental processes in astrophysics, driving angular momentum transport and mass accretion in a wide variety of cosmic objects. Despite much theoretical/numerical and experimental efforts over the last decades, its saturation mechanism and amplitude, which sets the angular momentum transport rate, remains not well understood, especially in the limit of high resistivity, or small magnetic Prandtl numbers typical to interiors (dead zones) of protoplanetary disks, liquid cores of planets and liquid metals in laboratory. Using direct numerical simulations, in this paper we investigate the nonlinear development and saturation properties of the helical magnetorotational instability (HMRI) -- a relative of the standard MRI -- in a magnetized Taylor-Couette flow at very low magnetic Prandtl number (correspondingly at low magnetic Reynolds number) relevant to liquid metals. For simplicity, the ratio of azimuthal field to axial field is kept fixed. From the linear theory of HMRI, it is known that the Elsasser number, or interaction parameter determines its growth rate and plays a special role in the dynamics. We show that this parameter is also important in the nonlinear problem. By increasing its value, a sudden transition from weakly nonlinear, where the system is slightly above the linear stability threshold, to strongly nonlinear, or turbulent regime occurs. We calculate the azimuthal and axial energy spectra corresponding to these two regimes and show that they differ qualitatively. Remarkably, the nonlinear state remains in all cases nearly axisymmetric suggesting that HMRI turbulence is quasi two-dimensional in nature. Although the contribution of non-axisymmetric modes increases moderately with the Elsasser number, their total energy remains much smaller than that of the axisymmetric ones.

Helical and azimuthal magnetorotational instabilities operate in rotating magnetized flows with relatively steep negative or extremely steep positive shear. The corresponding lower and upper Liu limits of the shear, which determine the threshold of modal growth of these instabilities, are continuously connected when some axial electrical current is allowed to pass through the rotating fluid. We investigate the nonmodal dynamics of these instabilities arising from the non-normality of shear flow in the local approximation, generalizing the results of the modal approach. It is demonstrated that moderate transient/nonmodal amplification of both types of magnetorotational instability occurs within the Liu limits, where the system is stable according to modal analysis. We show that for the helical magnetorotational instability this magnetohydrodynamic behavior is closely connected with the nonmodal growth of the underlying purely hydrodynamic problem.

The solar cycle appears to be remarkably synchronized with the gravitational torques exerted by the tidally dominant planets Venus, Earth and Jupiter. Recently, a possible synchronization mechanism was proposed that relies on the intrinsic helicity oscillation of the current-driven Tayler instability which can be stoked by tidal-like perturbations with a period of 11.07 years. Inserted into a simple alpha-Omega dynamo model these resonantly excited helicity oscillations led to a 22.14 years dynamo cycle. Here, we assess various alternative mechanisms of synchronization. Specifically we study a simple time-delay model of Babcock-Leighton type dynamos and ask whether periodic changes of either the minimal amplitude for rising toroidal flux tubes or the Omega effect could eventually lead to synchronization. In contrast to the easy and robust synchronizability of Tayler-Spruit dynamo models, our answer for those Babcock-Leighton type models is less propitious.

We examine kinematic dynamo action driven by an axisymmetric large scale flow that is superimposed with an azimuthally propagating non-axisymmetric perturbation with a frequency $\omega$. Although we apply a rather simple large scale velocity field, our simulations exhibit a complex behavior with oscillating and azimuthally drifting eigenmodes as well as stationary regimes. Within these non-oscillating regimes we find parametric resonances characterized by a considerable enhancement of dynamo action and by a locking of the phase of the magnetic field to the pattern of the perturbation. We find an approximate fulfillment of the relationship between the resonant frequency $\omega_{\rm{res}}$ of the disturbed system and the eigenfrequency $\omega_0$ of the undisturbed system given by $\omega_{\rm{res}}=2\omega_0$ which is known from paradigmatic rotating mechanical systems and our prior study [Giesecke et al., Phys. Rev. E, 86, 066303 (2012)]. We find further -- broader -- regimes with weaker enhancement of the growth rates but without phase locking. However, this amplification regime arises only in case of a basic (i.e. unperturbed) state consisting of several different eigenmodes with rather close growth rates. Qualitatively, these observations can be explained in terms of a simple low dimensional model for the magnetic field amplitude that is derived using Floquet theory. The observed phenomena may be of fundamental importance in planetary dynamo models with the base flow being disturbed by periodic external forces like precession or tides and for the realization of dynamo action under laboratory conditions where imposed perturbations with the appropriate frequency might facilitate the occurrence of dynamo action.

First results of a new spherical Couette experiment are presented. The liquid metal flow in a spherical shell is exposed to a homogeneous axial magnetic field. For a Reynolds number Re=1000, we study the effect of increasing Hartmann number Ha. The resulting flow structures are inspected by ultrasound Doppler velocimetry. With a weak applied magnetic field, we observe an equatorially anti-symmetric jet instability with azimuthal wave number m=3. As the magnetic field strength increases, this instability vanishes. When the field is increased further, an equatorially symmetric return flow instability arises. Our observations are shown to be in good agreement with linear stability analysis and non-linear flow simulations.

Magnetic fields of planets, stars and galaxies are generated by self-excitation in moving electrically conducting fluids. Once produced, magnetic fields can play an active role in cosmic structure formation by destabilizing rotational flows that would be otherwise hydrodynamically stable. For a long time, both hydromagnetic dynamo action as well as magnetically triggered flow instabilities had been the subject of purely theoretical research. Meanwhile, however, the dynamo effect has been observed in large-scale liquid sodium experiments in Riga, Karlsruhe and Cadarache. In this paper, we summarize the results of some smaller liquid metal experiments devoted to various magnetic instabilities such as the helical and the azimuthal magnetorotational instability, the Tayler instability, and the different instabilities that appear in a magnetized spherical Couette flow. We conclude with an outlook on a large scale Tayler-Couette experiment using liquid sodium, and on the prospects to observe magnetically triggered instabilities of flows with positive shear.

The current-driven, kink-type Tayler instability (TI) is a key ingredient of the Tayler-Spruit dynamo model for the generation of stellar magnetic fields, but is also discussed as a mechanism that might hamper the up-scaling of liquid metal batteries. Under some circumstances, the TI involves a helical flow pattern which goes along with some alpha effect. Here we focus on the chiral symmetry breaking and the related impact on the alpha effect that would be needed to close the dynamo loop in the Tayler-Spruit model. For low magnetic Prandtl numbers, we observe intrinsic oscillations of the alpha effect. These oscillations serve then as the basis for a synchronized Tayler-Spruit dynamo model, which could possibly link the periodic tidal forces of planets with the oscillation periods of stellar dynamos.

Recent years have seen an increased interest in the question of whether the gravitational action of planets could have an influence on the solar dynamo. Without discussing the observational validity of the claimed correlations, we ask for a possible physical mechanism that might link the weak planetary forces with solar dynamo action. We focus on the helicity oscillations that were recently found in simulations of the current-driven, kink-type Tayler instability, which is characterized by an m=1 azimuthal dependence. We show how these helicity oscillations can be resonantly excited by some m=2 perturbation that reflects a tidal oscillation. Specifically, we speculate that the 11.07 years tidal oscillation induced by the Venus--Earth--Jupiter system may lead to a 1:1 resonant excitation of the oscillation of the alpha-effect. Finally, in the framework of a reduced, zero-dimensional alpha--Omega dynamo model we recover a 22.14-year cycle of the solar dynamo.

The helical magnetorotational instability is known to work for resistive rotational flows with comparably steep negative or extremely steep positive shear. The corresponding lower and upper Liu limits of the shear are continuously connected when some axial electrical current is allowed to flow through the rotating fluid. Using a local approximation we demonstrate that the magnetohydrodynamic behavior of this dissipation-induced instability is intimately connected with the nonmodal growth and the pseudospectrum of the underlying purely hydrodynamic problem.

It is known that in a hydrodynamic Taylor-Couette system uniform rotation or a rotation law with positive shear ('super-rotation') are linearly stable. It is also known that a conducting fluid under the presence of a sufficiently strong axial electric-current becomes unstable against nonaxisymmetric disturbances. It is thus suggestive that a cylindric pinch formed by a homogeneous axial electric-current is stabilized by rotation laws with $d\Omega/dR \geq 0$. However, for magnetic Prandtl numbers Pm$\neq 1$ and for slow rotation also rigid rotation and super-rotation support the instability by lowering their critical Hartmann numbers. For super-rotation in narrow gaps and for modest rotation rates this double-diffusive instability even exists for toroidal magnetic fields with rather arbitrary radial profiles, the current-free profile $B_\phi\propto 1/R$ included. For rigid rotation and for super-rotation the sign of the azimuthal drift of the nonaxisymmetric hydromagnetic instability pattern strongly depends on the magnetic Prandtl number. The pattern counterrotates with the flow for Pm$\ll 1$ and it corotates for Pm$\gg 1$ while for rotation laws with negative shear the instability pattern migrates in the direction of the basic rotation for all Pm. An axial electric-current of minimal 3.6 kAmp flowing inside or outside the inner cylinder suffices to realize the double-diffusive instability for super-rotation in experiments using liquid sodium as the conducting fluid between the rotating cylinders. The limit is 11 kAmp if a gallium alloy is used.

According to Rayleigh's criterion, rotating flows are linearly stable when their specific angular momentum increases radially outward. The celebrated magnetorotational instability opens a way to destabilize those flows, as long as the angular velocity is decreasing outward. Using a short-wavelength approximation we demonstrate that even flows with very steep positive shear can be destabilized by azimuthal magnetic fields which are current-free within the fluid. We illustrate the transition of this instability to a rotationally enhanced kink-type instability in case of a homogeneous current in the fluid, and discuss the prospects for observing it in a magnetized Taylor-Couette flow.

We present results from three-dimensional non-linear hydrodynamic simulations of a precession driven flow in cylindrical geometry. The simulations are motivated by a dynamo experiment currently under development at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) in which the possibility of generating a magnetohydrodynamic dynamo will be investigated in a cylinder filled with liquid sodium and simultaneously rotating around two axes. In this study, we focus on the emergence of non-axisymmetric time-dependent flow structures in terms of inertial waves which - in cylindrical geometry - form so-called Kelvin modes. For a precession ratio ${\rm{Po}}=\Omega_p/\Omega_c=0.014$ the amplitude of the forced Kelvin mode reaches up to one fourth of the rotation velocity of the cylindrical container confirming that precession provides a rather efficient flow driving mechanism even at moderate values of ${\rm{Po}}$. More relevant for dynamo action might be free Kelvin modes with higher azimuthal wave number. These free Kelvin modes are triggered by non-linear interactions and may constitute a triadic resonance with the fundamental forced mode when the height of the container matches their axial wave lengths. Our simulations reveal triadic resonances at aspect ratios close to those predicted by the linear theory except around the primary resonance of the forced mode. In that regime we still identify various free Kelvin modes, however, all of them exhibit a retrograde drift around the symmetry axis of the cylinder and none of them can be assigned to a triadic resonance. The amplitudes of the free Kelvin modes always remain below the forced mode but may reach up to 6% of the of the container's angular velocity. The properties of the free Kelvin modes will be used in future simulations of the magnetic induction equation to investigate their ability to provide for dynamo action.

The Tayler instability is a kink-type, current driven instability that plays an important role in plasma physics but might also be relevant in liquid metal applications with high electrical currents. In the framework of the Tayler-Spruit dynamo model of stellar magnetic field generation, the question of spontaneous helical (chiral) symmetry breaking during the saturation of the Tayler instability has received considerable interest. Focusing on fluids with low magnetic Prandtl numbers, for which the quasistatic approximation can be applied, we utilize an integro-differential equation approach in order to investigate the saturation mechanism of the Tayler instability. Both the exponential growth phase and the saturated phase are analyzed in terms of the action of the alpha and beta effects of mean-field magnetohydrodynamics. In the exponential growth phase we always find a spontaneous chiral symmetry breaking which, however, disappears in the saturated phase. For higher degrees of supercriticality, we observe helicity oscillations in the saturated regime. For Lundquist numbers in the order of one we also obtain chiral symmetry breaking of the saturated magnetic field.

The next generation dynamo experiment currently under development at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) will consist of a precessing cylindrical container filled with liquid sodium. We perform numerical simulations of kinematic dynamo action applying a velocity field that is obtained from hydrodynamic models of a precession driven flow. So far, the resulting magnetic field growth-rates remain below the dynamo threshold for magnetic Reynolds numbers up to Rm=2000.

We perform a local stability analysis of rotational flows in the presence of a constant vertical magnetic field and an azimuthal magnetic field with a general radial dependence. Employing the short-wavelength approximation we develop a unified framework for the investigation of the standard, the helical, and the azimuthal version of the magnetorotational instability, as well as of current-driven kink-type instabilities. Considering the viscous and resistive setup, our main focus is on the case of small magnetic Prandtl numbers which applies, e.g., to liquid metal experiments but also to the colder parts of accretion disks. We show that the inductionless versions of MRI that were previously thought to be restricted to comparably steep rotation profiles extend well to the Keplerian case if only the azimuthal field slightly deviates from its current-free (in the fluid) profile. We find an explicit criterion separating the pure azimuthal inductionless magnetorotational instability from the regime where this instability is mixed with the Tayler instability. We further demonstrate that for particular parameter configurations the azimuthal MRI originates as a result of a dissipation-induced instability of the Chandrasekhar's equipartition solution of ideal magnetohydrodynamics.

The magnetohydrodynamic stability of axially unbounded cylindrical flows is considered which contain a toroidal magnetic background field with the same radial profile as the linear azimuthal velocity. Chandrasekhar (1956) has shown for ideal fluids the stability of this configuration if the Alfven velocity of the field equals the velocity of the background flow. It is demonstrated for magnetized Taylor-Couette flows at the Rayleigh line, however, that for finite diffusivity such flows become unstable against nonaxisymmetric perturbations where the critical magnetic Reynolds number of the rotation rate does not depend on the magnetic Prandtl number Pm if Pm much << 1. In order to study this new diffusive azimuthal magnetorotational instability, flows and fields with the same radial profile but with different amplitudes are considered. For Pm << 1 the instability domain with the weakest fields and the slowest rotation rates lies below the Chandrasekhar line of equal amplitudes for Alfven velocity and rotation velocity. We find that then the lines of marginal instability scale with the Reynolds number and the Hartmann number. The minimum values of the field strength and the rotation rate which are needed for the instability (slightly) grow for more and more flat rotation. Finally, the corresponding electric current of the background field becomes so strong that the Tayler instability (which even exists without rotation) also appears in the bifurcation map at small Hartmann numbers displacing after all the azimuthal magnetorotational instability.

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.

Using the WKB approximation we perform a linear stability analysis for a rotational flow of a viscous and electrically conducting fluid in an external azimuthal magnetic field that has an arbitrary radial profile B_{phi}(R). In the inductionless approximation, we find the growth rate of the three-dimensional perturbation in a closed form and demonstrate in particular that it can be positive when the velocity profile is Keplerian and the magnetic field profile is slightly shallower than R^{-1}.

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.

The magnetorotational instability (MRI) can destabilize hydrodynamically stable rotational flows, thereby allowing angular momentum transport in accretion disks. A notorious problem for MRI is its questionable applicability in regions with low magnetic Prandtl number, as they are typical for protoplanetary disks and the outer parts of accretion disks around black holes. Using the WKB method, we extend the range of applicability of MRI by showing that the inductionless versions of MRI, such as the helical MRI and the azimuthal MRI, can easily destabilize Keplerian profiles ~ 1/r^(3/2) if the radial profile of the azimuthal magnetic field is only slightly modified from the current-free profile ~ 1/r. This way we further show how the formerly known lower Liu limit of the critical Rossby number, Ro=-0.828, connects naturally with the upper Liu limit, Ro=+4.828.

In the last decades, the experimental study of dynamo action has made great progress. However, after the dynamo experiments in Karlsruhe and Riga, the von-Karman-Sodium (VKS) dynamo is only the third facility that has been able to demonstrate fluid flow driven self-generation of magnetic fields in a laboratory experiment. Further progress in the experimental examination of dynamo action is expected from the planned precession driven dynamo experiment that will be designed in the framework of the liquid sodium facility DRESDYN (DREsden Sodium facility for DYNamo and thermohydraulic studies). In this paper, we briefly present numerical models of the VKS dynamo that demonstrate the close relation between the axisymmetric field observed in that experiment and the soft iron material used for the flow driving impellers. We further show recent results of preparatory water experiments and design studies related to the precession dynamo and delineate the scientific prospects for the final set-up.

The electrical current through an incompressible, viscous and resistive liquid conductor produces an azimuthal magnetic field that becomes unstable when the corresponding Hartmann number exceeds a critical value in the order of 20. This Tayler instability, which is not only discussed as a key ingredient of a non-linear stellar dynamo model (Tayler-Spruit dynamo), but also as a limiting factor for the maximum size of large liquid metal batteries, was recently observed experimentally in a column of a liquid metal (Seilmayer et al., Phys. Rev. Lett. 108, 244501, 2012}. On the basis of an integro-differential equation approach, we have developed a fully three-dimensional numerical code, and have utilized it for the simulation of the Tayler instability at typical viscosities and resistivities of liquid metals. The resulting growth rates are in good agreement with the experimental data. We illustrate the capabilities of the code for the detailed simulation of liquid metal battery problems in realistic geometries.

We have performed numerical simulations of the kinematic induction equation in order to examine the dynamo efficiency of an axisymmetric von-K\'arm\'an-like flow subject to time-dependent non-axisymmetric velocity perturbations. The numerical model is based on the setup of the French Von-K\'arm\'an-Sodium dynamo (VKS) and on the flow measurements from a model water experiment conducted at the University of Navarra in Pamplona, Spain. Our simulations show that the interactions of azimuthally drifting flow perturbations with the fundamental drift of the magnetic eigenmode (caused by the inevitable equatorial symmetry breaking of the basic flow) essentially determine the temporal behavior of the dynamo state. We find two distinct regimes of dynamo action that depend on the (prescribed) drift frequency of an ($m=2$) vortex-like flow perturbation. For comparatively slowly drifting vortices we observe a narrow window with enhanced growth-rates and a drift of the magnetic eigenmode that is synchronized with the perturbation drift. The resonance-like enhancement of the growth-rates takes place when the vortex drift frequency roughly equals the drift frequency of the magnetic eigenmode in the unperturbed system. Outside of this small window, the field generation is hampered compared to the unperturbed case, and the field amplitude of the magnetic eigenmode is modulated with approximately twice the vortex drift frequency. The abrupt transition between the resonant regime and the modulated regime is identified as an spectral exceptional point where eigenvalues (growth-rates and frequencies) and eigenfunctions of two previously independent modes collapse.

Many astrophysical phenomena (such as the slow rotation of neutron stars or the rigid rotation of the solar core) can be explained by the action of the Tayler instability of toroidal magnetic fields in the radiative zones of stars. In order to place the theory of this instability on a safe fundament it has been realized in a laboratory experiment measuring the critical field strength, the growth rates as well as the shape of the supercritical modes. A strong electrical current flows through a liquid-metal confined in a resting columnar container with an insulating outer cylinder. As the very small magnetic Prandtl number of the gallium-indium-tin alloy does not influence the critical Hartmann number of the field amplitudes, the electric currents for marginal instability can also be computed with direct numerical simulations. The results of this theoretical concept are confirmed by the experiment. Also the predicted growth rates of the order of minutes for the nonaxisymmetric perturbations are certified by the measurements. That they do not directly depend on the size of the experiment is shown as a consequence of the weakness of the applied fields and the absence of rotation.

The magnetorotational instability (MRI) plays a key role for cosmic structure formation by triggering turbulence in the rotating flows of accretion disks that would be otherwise hydrodynamically stable. In the limit of small magnetic Prandtl number, the helical and the azimuthal version of MRI are known to be governed by a quite different scaling behaviour than the standard MRI with a vertical applied magnetic field. Using the short-wavelength approximation for an incompressible, resistive, and viscous rotating fluid we present a unified description of helical and azimuthal MRI, and we identify the universal character of the Liu limit Ro= -0.8284 for the critical Rossby number. From this universal behaviour we are also lead to the prediction of higher azimuthal wavenumber for rather small ratios of azimuthal to axial applied fields.

The eigenvalues and eigenfunctions of a linear {\alpha}^{2}-dynamo have been computed for different spatial distributions of an isotropic \alpha-effect. Oscillatory solutions are obtained when \alpha exhibits a sign change in the radial direction. The time-dependent solutions arise at so called exceptional points where two stationary modes merge and continue as an oscillatory eigenfunction with conjugate complex eigenvalues. The close proximity of oscillatory and non-oscillatory solutions may serve as the basic ingredient for reversal models that describe abrupt polarity switches of a dipole induced by noise. Whereas the presence of an inner core with different magnetic diffusivity has remarkable little impact on the character of the dominating dynamo eigenmodes, the introduction of equatorial symmetry breaking considerably changes the geometric character of the solutions. Around the dynamo threshold the leading modes correspond to hemispherical dynamos even when the symmetry breaking is small. This behavior can be explained by the approximate dipole-quadrupole degeneration for the unperturbed problem. More complicated scenarios may occur in case of more realistic anisotropies of \alpha- and \beta-effect or through non-linearities caused by the back-reaction of the magnetic field (magnetic quenching).

In the current-driven, kink-type Tayler instability (TI) a sufficiently strong azimuthal magnetic field becomes unstable against non-axisymmetric perturbations. The TI has been discussed as a possible ingredient of the solar dynamo mechanism and a source of the helical structures in cosmic jets. It is also considered as a size limiting factor for liquid metal batteries. We report on a liquid metal TI experiment using a cylindrical column of the eutectic alloy GaInSn to which electrical currents of up to 8 kA are applied. We present results of external magnetic field measurements that indicate the occurrence of the TI in good agreement with numerical predictions. The interference of TI with the competing large scale convection, resulting from Joule heating, is also discussed.

The magnetorotational instability (MRI) triggers turbulence and enables outward transport of angular momentum in hydrodynamically stable rotating shear flows, e.g., in accretion disks. What laws of differential rotation are susceptible to the destabilization by axial, azimuthal, or helical magnetic field? The answer to this question, which is vital for astrophysical and experimental applications, inevitably leads to the study of spectral and geometrical singularities on the instability threshold. The singularities provide a connection between seemingly discontinuous stability criteria and thus explain several paradoxes in the theory of MRI that were poorly understood since the 1950s.

The magnetorotational instability (MRI) triggers turbulence and enables outward transport of angular momentum in hydrodynamically stable accretion discs. By using the WKB approximation and methods of singular function theory, we resolve two different paradoxes of MRI that appear in the limits of infinite and vanishing magnetic Prandtl number. For the latter case we derive a new strict limit of the critical Rossby number. This new limit of ${\rm Ro_c}=-0.802$, which appears for a finite Lundquist number of ${\rm Lu}=0.618$, extends the formerly known inductionless Liu limit of ${\rm Ro_c}=-0.828$ valid at ${\rm Lu}=0$.

Starting from the present version of the Riga dynamo experiment with its rotating magnetic eigenfield dominated by a single frequency we ask for those modifications of this set-up that would allow for a non-trivial magnetic field behaviour in the saturation regime. Assuming an increased ratio of azimuthal to axial flow velocity, we obtain energy oscillations with a frequency below the eigenfrequency of the magnetic field. These new oscillations are identified as magneto-inertial waves that result from a slight imbalance of Lorentz and inertial forces. Increasing the azimuthal velocity further, or increasing the total magnetic Reynolds number, we find transitions to a chaotic behaviour of the dynamo.

The magnetorotational instability (MRI) plays a crucial role for cosmic structure formation by enabling turbulence in Keplerian disks which would be otherwise hydrodynamically stable. With particular focus on MRI experiments with liquid metals, which have small magnetic Prandtl numbers, it has been shown that the helical version of this instability (HMRI) has a scaling behaviour that is quite different from that of the standard MRI (SMRI). We discuss the relation of HMRI to SMRI by exploring various parameter dependencies. We identify the mechanism of transfer of instability between modes through a spectral exceptional point that explains both the transition from a stationary instability (SMRI) to an unstable travelling wave (HMRI) and the excitation of HMRI in the inductionless limit. For certain parameter regions we find new islands of the HMRI.

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.

Using a homotopic family of boundary eigenvalue problems for the mean-field $\alpha^2$-dynamo with helical turbulence parameter $\alpha(r)=\alpha_0+\gamma\Delta\alpha(r)$ and homotopy parameter $\beta \in [0,1]$, we show that the underlying network of diabolical points for Dirichlet (idealized, $\beta=0$) boundary conditions substantially determines the choreography of eigenvalues and thus the character of the dynamo instability for Robin (physically realistic, $\beta=1$) boundary conditions. In the $(\alpha_0,\beta,\gamma)-$space the Arnold tongues of oscillatory solutions at $\beta=1$ end up at the diabolical points for $\beta=0$. In the vicinity of the diabolical points the space orientation of the 3D tongues, which are cones in first-order approximation, is determined by the Krein signature of the modes involved in the diabolical crossings at the apexes of the cones. The Krein space induced geometry of the resonance zones explains the subtleties in finding $\alpha$-profiles leading to spectral exceptional points, which are important ingredients in recent theories of polarity reversals of the geomagnetic field.

It is widely known that cosmic magnetic fields, i.e. the fields of planets, stars, and galaxies, are produced by the hydromagnetic dynamo effect in moving electrically conducting fluids. It is less well known that cosmic magnetic fields play also an active role in cosmic structure formation by enabling outward transport of angular momentum in accretion disks via the magnetorotational instability (MRI). Considerable theoretical and computational progress has been made in understanding both processes. In addition to this, the last ten years have seen tremendous efforts in studying both effects in liquid metal experiments. In 1999, magnetic field self-excitation was observed in the large scale liquid sodium facilities in Riga and Karlsruhe. Recently, self-excitation was also obtained in the French "von Karman sodium" (VKS) experiment. An MRI-like mode was found on the background of a turbulent spherical Couette flow at the University of Maryland. Evidence for MRI as the first instability of an hydrodynamically stable flow was obtained in the "Potsdam Rossendorf Magnetic Instability Experiment" (PROMISE). In this review, the history of dynamo and MRI related experiments is delineated, and some directions of future work are discussed.

The asymmetric time dependence and various statistical properties of polarity reversals of the Earth's magnetic field are utilized to infer some of the most essential parameters of the geodynamo, among them the effective (turbulent) magnetic diffusivity, the degree of supercriticality, and the relative strength of the periodic forcing which is believed to result from the Milankovic cycle of the Earth's orbit eccentricity. A time-stepped spherically symmetric alpha^2-dynamo model is used as the kernel of an inverse problem solver in form of a downhill simplex method which converges to solutions that yield a stunning correspondence with paleomagnetic data.

On 11 November 1999, a self-exciting magnetic eigenfield was detected for the first time in the Riga liquid sodium dynamo experiment. We report on the long history leading to this event, and on the subsequent experimental campaigns which provided a wealth of data on the kinematic and the saturated regime of this dynamo. The present state of the theoretical understanding of both regimes is delineated, and some comparisons with other laboratory dynamo experiments are made.

Using a spherical symmetric mean field alpha^2-dynamo model for Earth's magnetic field reversals, we show the coexistence of the noise-induced phenomena coherence resonance and stochastic resonance. Stochastic resonance has been recently invoked to explain the 100 kyr periodicity in the distribution of the residence time between reversals. The comparison of the resulting residence time distribution with the paleomagnetic one allows for some estimate of the effective diffusion time of the Earth's core which may be 100 kyr or slightly below rather than 200 kyr as it would result from the molecular resistivity.

The conventional magnetic induction equation that governs hydromagnetic dynamo action is transformed into an equivalent integral equation system. An advantage of this approach is that the computational domain is restricted to the region occupied by the electrically conducting fluid and to its boundary. This integral equation approach is first employed to simulate kinematic dynamos excited by Beltrami-like flows in a finite cylinder. The impact of externally added layers around the cylinder on the onset of dynamo actions is investigated. Then it is applied to simulate dynamo experiments within cylindrical geometry including the von Karman sodium (VKS) experiment and the Riga dynamo experiment. A modified version of this approach is utilized to investigate magnetic induction effects under the influence of externally applied magnetic fields which is also important to measure the proximity of a given dynamo facility to the self-excitation threshold.

A brief overview is given over recent results on the spectral properties of spherically symmetric MHD alpha2-dynamos. In particular, the spectra of sphere-confined fluid or plasma configurations with physically realistic boundary conditions (BCs) (surrounding vacuum) and with idealized BCs (super-conducting surrounding) are discussed. The subjects comprise third-order branch points of the spectrum, self-adjointness of the dynamo operator in a Krein space as well as the resonant unfolding of diabolical points. It is sketched how certain classes of dynamos with a strongly localized alpha-profile embedded in a conducting surrounding can be mode decoupled by a diagonalization of the dynamo operator matrix. A mapping of the dynamo eigenvalue problem to that of a quantum mechanical Hamiltonian with energy dependent potential is used to obtain qualitative information about the spectral behavior. Links to supersymmetric Quantum Mechanics and to the Dirac equation are indicated.

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.

The asymmetric shape of reversals of the Earth's magnetic field indicates a possible connection with relaxation oscillations as they were early discussed by van der Pol. A simple mean-field dynamo model with a spherically symmetric $\alpha$ coefficient is analysed with view on this similarity, and a comparison of the time series and the phase space trajectories with those of paleomagnetic measurements is carried out. For highly supercritical dynamos a very good agreement with the data is achieved. Deviations of numerical reversal sequences from Poisson statistics are analysed and compared with paleomagnetic data. The role of the inner core is discussed in a spectral theoretical context and arguments and numerical evidence is compiled that the growth of the inner core might be important for the long term changes of the reversal rate and the occurrence of superchrons.

A new class of semi-analytically solvable MHD alpha^2-dynamos is found based on a global diagonalization of the matrix part of the dynamo differential operator. Close parallels to SUSY QM are used to relate these models to the Dirac equation and to extract non-numerical information about the dynamo spectrum.

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.

The intention of the ''von Karman sodium'' (VKS) experiment is to study the hydromagnetic dynamo effect in a highly turbulent and unconstrained flow. Much effort has been devoted to the optimization of the mean flow and the lateral boundary conditions in order to minimize the critical magnetic Reynolds number and hence the necessary motor power. The main focus of this paper lies on the role of ''lid layers'', i.e. layers of liquid sodium between the impellers and the end walls of the cylinder. First, we study an analytical test flow to show that lid layers can have an ambivalent effect on the efficiency of the dynamo. The critical magnetic Reynolds number shows a flat minimum for a small lid layer thickness, but increases for thicker layers. For the actual VKS geometry it is shown that static lid layers yield a moderate increase of the critical magnetic Reynolds number by approximately 12 per cent. A more dramatic increase by 100 until 150 per cent can occur when some rotational flow is taken into account in those layers. Possible solutions of this problem are discussed for the real dynamo facility.

The irregular polarity reversals of the Earth's magnetic field have attracted much interest during the last decades. Despite the fact that recent numerical simulations of the geodynamo have shown nice polarity transitions, the very reason and the basic mechanism of reversals are far from being understood. Using a paradigmatic mean-field dynamo model with a spherically symmetric helical turbulence parameter alpha we attribute the essential features of reversals to the magnetic field dynamics in the vicinity of an exceptional point of the spectrum of the non-selfadjoint dynamo operator. At such exceptional (branch) points of square root type two real eigenvalues coalesce and continue as a complex conjugated pair of eigenvalues. Special focus is laid on the comparison of numerically computed time series with paleomagnetic observations. It is shown that the considered dynamo model with high supercriticality can explain the observed time scale and the asymmetric shape of reversals with a slow decay and a fast field recovery.

In a recent paper (Phys. Rev. Lett. 94 (2005), 184506; physics/0411050) it was shown that a simple mean-field dynamo model with a spherically symmetric helical turbulence parameter alpha can exhibit a number of features which are typical for Earth's magnetic field reversals. In particular, the model produces asymmetric reversals, a positive correlation of field strength and interval length, and a bimodal field distribution. All these features are attributable to the magnetic field dynamics in the vicinity of an exceptional point of the spectrum of the non-selfadjoint dynamo operator. The negative slope of the growth rate curve between the nearby local maximum and the exceptional point makes the system unstable and drives it to the exceptional point and beyond into the oscillatory branch where the sign change happens. A weakness of this reversal model is the apparent necessity to fine-tune the magnetic Reynolds number and/or the radial profile of alpha. In the present paper, it is shown that this fine-tuning is not necessary in the case of higher supercriticality of the dynamo. Numerical examples and physical arguments are compiled to show that, with increasing magnetic Reynolds number, there is strong tendency for the exceptional point and the associated local maximum to move close to the zero growth rate line. Although exemplified again by the spherically symmetric alpha^2 dynamo model, the main idea of this ''self-tuning'' mechanism of saturated dynamos into a reversal-prone state seems well transferable to other dynamos. As a consequence, reversing dynamos might be much more typical and may occur much more frequently in nature than what could be expected from a purely kinematic perspective.

The spectra of self-adjoint operators in Krein spaces are known to possess real sectors as well as sectors of pair-wise complex conjugate eigenvalues. Transitions from one spectral sector to the other are a rather generic feature and they usually occur at exceptional points of square root branching type. For certain parameter configurations two or more such exceptional points may happen to coalesce and to form a higher order branch point. We study the coalescence of two square root branch points semi-analytically for a PT-symmetric 4x4 matrix toy model and illustrate numerically its occurrence in the spectrum of the 2x2 operator matrix of the magnetohydrodynamic alpha^2-dynamo and of an extended version of the hydrodynamic Squire equation.

It is shown that the alpha^2-dynamo of Magnetohydrodynamics, the hydrodynamic Squire equation as well as an interpolation model of PT-symmetric Quantum Mechanics are closely related as spectral problems in Krein spaces. For the alpha^2-dynamo and the PT-symmetric model the strong similarities are demonstrated with the help of a 2x2 operator matrix representation, whereas the Squire equation is re-interpreted as a rescaled and Wick-rotated PT-symmetric problem. Based on recent results on the Squire equation the spectrum of the PT-symmetric interpolation model is analyzed in detail and the Herbst limit is described as spectral singularity.

Using a mean-field dynamo model with a spherically symmetric helical turbulence parameter alpha which is dynamically quenched and disturbed by additional noise, the basic features of geomagnetic polarity reversals are shown to be generic consequences of the dynamo action in the vicinity of exceptional points of the spectrum. This simple paradigmatic model yields long periods of constant polarity which are interrupted by self-accelerating field decays leading to asymmetric polarity reversals. It shows the recently discovered bimodal field distribution, and it gives a natural explanation of the correlation between polarity persistence time and field strength. In addition, we find typical features of coherence resonance in the dependence of the persistence time on the noise.

The spectral branching behavior of the 2x2 operator matrix of the magneto-hydrodynamic alpha^2-dynamo is analyzed numerically. Some qualitative aspects of level crossings are briefly discussed with the help of a simple toy model which is based on a Z_2-graded-pseudo-Hermitian 2x2 matrix. The considered issues comprise: the underlying SU(1,1) symmetry and the Krein space structure of the system, exceptional points of branching type and diabolic points, as well as the algebraic and geometric multiplicity of corresponding degenerate eigenvalues.

With only a few exceptions, the numerical simulation of cosmic and laboratory hydromagnetic dynamos has been carried out in the framework of the differential equation method. However, the integral equation method is known to provide robust and accurate tools for the numerical solution of many problems in other fields of physics. The paper is intended to facilitate the use of integral equation solvers in dynamo theory. In concrete, the integral equation method is employed to solve the eigenvalue problem for a hydromagnetic dynamo model with a spherically symmetric, isotropic helical turbulence parameter alpha. Three examples of the function alpha(r) with steady and oscillatory solutions are considered. A convergence rate proportional to the inverse squared of the number of grid points is achieved. Based on this method, a convergence accelerating strategy is developed and the convergence rate is improved remarkably. Typically, quite accurate results can be obtained with a few tens of grid points. In order to demonstrate its suitability for the treatment of dynamos in other than spherical domains, the method is also applied to alpha^2 dynamos in rectangular boxes. The magnetic fields and the electric potentials for the first eigenvalues are visualized.

The homogeneous dynamo effect is at the root of cosmic magnetic field generation. With only a very few exceptions, the numerical treatment of homogeneous dynamos is carried out in the framework of the differential equation approach. The present paper tries to facilitate the use of integral equations in dynamo research. Apart from the pedagogical value to illustrate dynamo action within the well-known picture of Biot-Savart's law, the integral equation approach has a number of practical advantages. The first advantage is its proven numerical robustness and stability. The second and perhaps most important advantage is its applicability to dynamos in arbitrary geometries. The third advantage is its intimate connection to inverse problems relevant not only for dynamos but also for technical applications of magnetohydrodynamics. The paper provides the first general formulation and application of the integral equation approach to time-dependent kinematic dynamos in finite domains. For the spherically symmetric alpha^2 dynamo model it is shown how the general formulation is reduced to a coupled system of two radial integral equations for the defining scalars of the poloidal and the toroidal field components. The integral equation formulation for spherical dynamos with general velocity fields is also derived. Two numerical examples, the alpha^2 dynamo model with radially varying alpha, and the Bullard-Gellman model illustrate the equivalence of the approach with the usual differential equation method.

The isospectrality problem is studied for the operator of the spherical hydromagnetic alpha^2-dynamo. It is shown that this operator is formally pseudo-Hermitian (J-symmetric) and lives in a Krein space. Based on the J-symmetry, an operator intertwining Ansatz with first-order differential intertwining operators is tested for its compatibility with the structure of the alpha^2-dynamo operator matrix. An intrinsic structural inconsistency is obtained in the set of associated matrix Riccati equations. This inconsistency is interpreted as a no-go theorem which forbids the construction of isospectral alpha^2-dynamo operator classes with the help of first-order differential intertwining operators.

Until recently, the existence of oscillatory mean-field dynamos of the alpha^2-type with spherically symmetric and isotropic alpha was an open question. We find such dynamos by means of an evolutionary strategy, and we illustrate the spectral properties of the corresponding dynamo operators.

After the dynamo experiment in November 1999 had shown magnetic field self-excitation in a spiraling liquid metal flow, in a second series of experiments emphasis was placed on the magnetic field saturation regime as the next principal step in the dynamo process. The dependence of the strength of the magnetic field on the rotation rate is studied. Various features of the saturated magnetic field are outlined and possible saturation mechanisms are discussed.

For a simple spherically symmetric mean-field dynamo model we investigate the possibility of determining the radial dependence of the coefficient $\alpha$. Growth rates for different magnetic field modes are assumed to be known by measurement. An evolutionary strategy (ES) is used for the solution of the inverse problem. Numerically, we find quite different $\alpha$-profiles giving nearly the same eigenvalues. The ES is also applied to find functions $\alpha(r)$ yielding zero growth rates for the lowest four magnetic field modes. Additionally, a slight modification of the ES is utilized for an ''energetic'' optimization of $\alpha^2$-dynamos. The consequences of our findings for inverse dynamo theory and for the design of future dynamo experiments are discussed.

The paper deals with the integral equation approach to steady kinematic dynamo models in finite domains based on Biot-Savart's law. The role of the electric potential at the boundary is worked out explicitly. As an example, a modified version of the simple spherical $\alpha$-effect dynamo model proposed by Krause and Steenbeck is considered in which the $\alpha$-coefficient is no longer constant but may vary with the radial coordinate. In particular, the results for the original model are re-derived. Possible applications of this integral equation approach for numerical simulations of dynamos in arbitrary geometry and for an ''inverse dynamo theory'' are sketched.

In an experiment at the Riga sodium dynamo facility, a slowly growing magnetic field eigenmode has been detected over a period of about 15 seconds. For a slightly decreased propeller rotation rate, additional measurements showed a slow decay of this mode. The measured results correspond satisfactory with numerical predictions for the growth rates and frequencies.