Convective stability in magnetic stars

Dynamical stability of a static axisymmetrical magnetic star with respect to high-order modes of oscillation is investigated by means of the energy method, neglecting the Eulerian perturbation of gravity. The magnetic field is assumed to be continuous across the surface of the star and its first-order spatial derivatives, but it may have both toroidal and poloidal components.The second variation of the potential energy is written in a way which, in the case of apurely toroidal field, and for axisymmetrical and non-axisymmetrical modes, yields Tayler's local stability criteria which are necessary and sufficient conditions for convective stability, and in the case of ageneral field yields a single local stability criterion, which is a sufficient condition for convective stability.     

Simple tests for the ideal MHD stability of line-tied coronal magnetic fields

Methods for investigating the stability of line-tied, cylindrically-symmetric magnetic fields are presented. The energy method is used and the perturbed potential energy integral is manipulated to produce simple tests that predict either stability to general coronal disturbances or instability to localized modes, both satisfying photospheric line-tying. Using these tests the maximum amount of magnetic energy, that can be stored in the coronal magnetic field prior to an instability, can be estimated. The tests are applied to four different classes of equilibria and results are obtained for both arcade and loop geometries.     

Waves and instabilities in a differentially rotating disc containing a poloidal magnetic field

The theory of waves and instabilities in a differentially rotating disc containing a poloidal magnetic field is developed within the framework of ideal magnetohydrodynamics. A continuous spectrum, for which the eigenfunctions are localized on individual magnetic surfaces, is identified but is found not to contain any instabilities associated with differential rotation. The normal modes of a weakly magnetized thin disc are studied by extending the asymptotic methods used previously to describe the equilibria. Waves propagate radially in the disc according to a dispersion relation which is determined by solving an eigenvalue problem at each radius. The dispersion relation for a hydrodynamic disc is re-examined and the modes are classified according to their behaviour in the limit of large wavenumber. The addition of a magnetic field introduces new, potentially unstable, modes and also breaks up the dispersion diagram by causing avoided crossings. The stability boundary to the magnetorotational instability in the parameter space of polytropic equilibria is located by solving directly for marginally stable equilibria. For a given vertical magnetic field in the disc, bending of the field lines has a stabilizing effect and it is shown that stable equilibria exist which are capable of launching a predominantly centrifugally driven wind.     

Mean-field dynamo with cubic non-linearity

The turbulent mean-field dynamo of αω-type with a mean helicity quadratically dependent on the magnetic field is investigated. A nonlinear system of ordinary differential equation is derived for the amplitudes of the magnetic field expansion over the eigenvectors of the linear problem. In a one-mode approximation the non-linear supercritical solution is stable when dγ/d D > 0, where γ is the growth rate of the linear solution and D is the dynamo number. Non-linear interation between two modes of dipole and quadrupole symmetry is considered. The conditions are found for the synchronization and beats of these modes under the assumption that the quadrupole mode is weaker than the dipole one.     

Stability of magnetic equilibria in radio bubbles

Current-carrying flows, in the laboratory and in astrophysical jets, can form remarkably stable magnetic structures. Decades of experience show that such flows often build equilibria that reverse field directions, evolving to a magnetohydrodynamic (MHD) Taylor state, which has remarkable stability properties. We model jets and the magnetic bubbles they build as reversed-field pinch equilibria by assuming the driver current to be stiff in the MHD sense. Taking the jet current as rigid and a fixed function of position, we prove a theorem: that the same, simple MHD stability conditions guarantee stability, even after the jet turns off. This means that magnetic structures harbouring a massive inventory of magnetic energy can persist long after the building jet current has died away. These may be the relic radio 'fossils', 'ghost bubbles' or 'magnetic balloons' found in clusters. These equilibria, which are under magnetic tension, will evolve, retaining the stability properties from that state. The remaining fossil is not a disordered ball of magnetic fields, but a stable structure under tension, able to respond to the slings and arrows of outside forces. Typically their Alfvén speeds greatly exceed the cluster sound speed, and so they can keep out hot cluster plasma, leading to X-ray ghosts. Passing shocks cannot easily destroy them, but can energize and light them up anew at radio frequencies. Bubbles can rise in the hot cluster plasma, perhaps detaching from the parent radio galaxy but stable against Rayleigh–Taylor and other modes.     

The effect of a poloidal magnetic field on the linear and adiabatic oscillations of a polytropen=3

The frequencies of the linear and adiabatic oscillations of a gaseous polytrope with a poloidal magnetic field are determined with the aid of a perturbation method. The influence of the poloidal magnetic field on the different types of spheroidal oscillation modes is discussed. The poloidal magnetic field generally strengthens the stability of the oscillation modes and this effect is the largest in the case of the non-radialp-modes.     

Bounds on the stability of 3D magnetic equilibria in the solar corona

A necessary and a sufficient condition are derived for the ideal magnetohydrodynamic stability of any 3D magnetohydrostatic equilibrium using the energy method and incorporating photospheric line-tying. The theory is demonstrated by application to a simple class of theoretical 3D equilibria. The main thrust of the method is the formulation of the stability conditions as two sets of ordinary differential equations together with appropriate boundary conditions which may be numerically integrated along tied field lines one at a time. In the case of the shearless fields with non-negligible plasma pressure treated here the conditions for stability arenecessary and sufficient. The method employs as a trial function a destabilizing ballooning mode, of large wave number vector perpendicular to the equilibrium field lines. These modes may not be picked up in a solution of the full partial differential equations which arise from a direct treatment of the problem.     

Stability of the classical type of relative equilibria of a rigid body in the J 2 problem

The motion of a point mass in the J ₂ problem is generalized to that of a rigid body in a J ₂ gravity field. The linear and nonlinear stability of the classical type of relative equilibria of the rigid body, which have been obtained in our previous paper, are studied in the framework of geometric mechanics with the second-order gravitational potential. Non-canonical Hamiltonian structure of the problem, i.e., Poisson tensor, Casimir functions and equations of motion, are obtained through a Poisson reduction process by means of the symmetry of the problem. The linear system matrix at the relative equilibria is given through the multiplication of the Poisson tensor and Hessian matrix of the variational Lagrangian. Based on the characteristic equation of the linear system matrix, the conditions of linear stability of the relative equilibria are obtained. The conditions of nonlinear stability of the relative equilibria are derived with the energy-Casimir method through the projected Hessian matrix of the variational Lagrangian. With the stability conditions obtained, both the linear and nonlinear stability of the relative equilibria are investigated in details in a wide range of the parameters of the gravity field and the rigid body. We find that both the zonal harmonic J ₂ and the characteristic dimension of the rigid body have significant effects on the linear and nonlinear stability. Similar to the classical attitude stability in a central gravity field, the linear stability region is also consisted of two regions that are analogues of the Lagrange region and the DeBra-Delp region respectively. The nonlinear stability region is the subset of the linear stability region in the first quadrant that is the analogue of the Lagrange region. Our results are very useful for the studies on the motion of natural satellites in our solar system.     

Magnetic Energy Storage and Current Density Distributions for Different Force-Free Models

In the last decades, force-free-field modelling has been used extensively to describe the coronal magnetic field and to better understand the physics of solar eruptions at different scales. Especially the evolution of active regions has been studied by successive equilibria in which each computed magnetic configuration is subject to an evolving photospheric distribution of magnetic field and/or electric-current density. This technique of successive equilibria has been successful in describing the rate of change of the energetics for observed active regions. Nevertheless the change in magnetic configuration due to the increase/decrease of electric current for different force-free models (potential, linear and nonlinear force-free fields) has never been studied in detail before. Here we focus especially on the evolution of the free magnetic energy, the location of the excess of energy, and the distribution of electric currents in the corona. For this purpose, we use an idealised active region characterised by four main polarities and a satellite polarity, allowing us to specify a complex topology and sheared arcades to the coronal magnetic field but no twisted flux bundles. We investigate the changes in the geometry and connectivity of field lines, the magnetic energy and current-density content as well as the evolution of null points. Increasing the photospheric current density in the magnetic configuration does not dramatically change the energy-storage processes within the active region even if the magnetic topology is slightly modified. We conclude that for reasonable values of the photospheric current density (the force-free parameter α<0.25 Mm⁻¹), the magnetic configurations studied do change but not dramatically: i) the original null point stays nearly at the same location, ii) the field-line geometry and connectivity are slightly modified, iii) even if the free magnetic energy is significantly increased, the energy storage happens at the same location. This extensive study of different force-free models for a simple magnetic configuration shows that some topological elements of an observed active region, such as null points, can be reproduced with confidence only by considering the potential-field approximation. This study is a preliminary work aiming at understanding the effects of electric currents generated by characteristic photospheric motions on the structure and evolution of the coronal magnetic field.     

The effects on the MHD stability of field line tying to the end faces of a cylindrical magnetic loop

We examine the magnetohydrodynamic (MHD) stability of a magnetic loop, taking into account field line tying at its foot points. We use the ideal MHD energy equation to derive a stability equation for a specific class of perturbations.We found that for a loop with large aspect ratio (10) the field line tying effect is negligible to the m = 1 kink mode but important to the localized modes. The stability criterion for high m localized modes is derived and compared with the Suydam criterion. The result shows that for the perturbation of the class studied, there are two effects of field line tying; one is a field line bending effect which is always stabilizing and the other is a shear effect which is stabilizing or destabilizing depending on the sign of the gradient of potential magnetic field. The net effect of field line tying is determined by the sum of these two effects.The result of this work is contrary to the result of Hood and Priest, in which they found that the field line tying effect is significant to the m = 1 mode. We believe that the contradiction comes from their incomplete minimization of the energy equation.     

Axisymmetric magnetic fields in stars: relative strengths of poloidal and toroidal components

In this third paper in a series on stable magnetic equilibria in stars, I look at the stability of axisymmetric field configurations and, in particular, the relative strengths of the toroidal and poloidal components. Both toroidal and poloidal fields are unstable on their own, and stability is achieved by adding the two together in some ratio. I use Tayler's stability conditions for toroidal fields and other analytic tools to predict the range of stable ratios and then check these predictions by running numerical simulations. If the energy in the poloidal component as a fraction of the total magnetic energy is written as E_p / E , it is found that the stability condition is a ( E / U ) < E_p / E ≲ 0.8 where E /U is the ratio of magnetic to gravitational energy in the star and a is some dimensionless factor whose value is of order 10 in a main-sequence star and of order 10³ in a neutron star. In other words, whilst the poloidal component cannot be significantly stronger than the toroidal, the toroidal field can be very much stronger than the poloidal–given that in realistic stars we expect E / U < 10⁻⁶. The implications of this result are discussed in various contexts such as the emission of gravitational waves by neutron stars, free precession and a 'hidden' energy source for magnetars.     

Axisymmetry vs. nonaxisymmetry of a Taylor‐Couette flow with azimuthal magnetic fields

The instability of a supercritical Taylor‐Couette flow of a conducting fluid with resting outer cylinder under the influence of a uniform axial electric current is investigated for magnetic Prandtl number Pm = 1. In the linear theory the critical Reynolds number for axisymmetric perturbations is not influenced by the current‐induced axisymmetric magnetic field but all axisymmetric magnetic perturbations decay. The nonaxisymmetric perturbations with m = 1 are excited even without rotation for large enough Hartmann numbers (“Tayler instability”). For slow rotation their growth rates scale with the Alfvén frequency of the magnetic field but for fast rotation they scale with the rotation rate of the inner cylinder. In the nonlinear regime the ratio of the energy of the magnetic m = 1 modes and the toroidal background field is very low for the non‐rotating Tayler instability but it strongly grows if differential rotation is present. For super‐Alfv´enic rotation the energies in the m = 1 modes of flow and field do not depend on the molecular viscosity, they are almost in equipartition and contain only 1.5 % of the centrifugal energy of the inner cylinder. The geometry of the excited magnetic field pattern is strictly nonaxisymmetric for slow rotation but it is of the mixed‐mode type for fast rotation – contrary to the situation which has been observed at the surface of Ap stars. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The visco-resistive stability of arcades in the corona of the Sun

The stability of ballooning modes in coronal arcades is studied using linear visco-resistive MHD. Rigid wall conditions are adopted for modelling the photospheric line-tying of the magnetic field. The full Braginskii viscosity stress tensor is used and particular attention is given to the effect of the viscosity coefficient ₃ which was left out of an earlier investigation by Van der Linden, Goossens, and Hood (1987, 1988). The numerical results for shearless arcades show that the coefficient ₃ has a stabilizing effect. However, for realistic values of the equilibrium quantities the stabilizing effect by ₃ can be neglected in comparison with the strong stabilizing effect of the perpendicular viscosity. The effect of magnetic field strength and mode number on stability are determined. In particular it is found that there exists a critical field strength for every mode number such that the mode is stable for weaker fields and unstable for stronger fields.     

Linear MHD Sausage Waves in Compressible Magnetically Twisted Flux Tubes

Oscillations of magnetic flux tubes are of great importance as they contain information about the geometry and fine structure of the flux tubes. Here we derive and analytically solve in terms of Kummer’s functions the linear governing equations of wave propagation for sausage surface and body modes (m=0) of a magnetically twisted compressible flux tube embedded in a compressible uniformly magnetized plasma environment in cylindrical geometry. A general dispersion relation is obtained for such flux tubes. Numerical solutions for the phase velocity are obtained for a wide range of wavenumbers and for varying magnetic twist. The effect of magnetic twist on the period of oscillations of sausage surface modes for different values of the wavenumber and vertical magnetic field strength is calculated for representative photospheric and coronal conditions. These results generalize and extend previous studies of MHD waves obtained for incompressible or for compressible but nontwisted flux tubes. It is found that magnetic twist may change the period of sausage surface waves of the order of a few percent when compared to counterparts in straight nontwisted flux tubes. This information will be most relevant when high-resolution observations are used for diagnostic exploration of MHD wave guides in analogy to solar-interior studies by means of global eigenoscillations in helioseismology.     

Force-free magnetic arcades relevant to two-ribbon solar flares

Simple analytic models for the passive evolution of arcade-like magnetic fields through a series of force-free equilibria are presented. At the photospheric boundary, the normal magnetic field component is prescribed together with either the longitudinal field component or the photospheric shear. Analytic progress is made by considering either cylindrically symmetric solutions or using the separation of variables technique. Two distinct cylindrically symmetric force-free fields are obtained that possess the same normal field component and photospheric shear. The scond field contains a magnetic bubble. As the shear increases beyond a critical value, so the magnetic energy of the first configuration exceeds that of the second. The possibility is therefore suggested of an eruption of the first field outwards towards the second. Such an eruptive instability is proposed as the origin of a two-ribbon solar flare.A new analytic solution to the force-free field equations, of separable form, is discovered and it is pointed out that the existence of shear in a magnetic field does not preclude it from being potential.Now at AWRE, Aldermaston, Reading, Berkshire.     

Stability of relative equilibria in arbitrary axisymmetric gravitational and magnetic fields

Relative equilibria occur in a wide variety of physical applications, including celestial mechanics, particle accelerators, plasma physics, and atomic physics. We derive sufficient conditions for Lyapunov stability of circular orbits in arbitrary axisymmetric gravitational (electrostatic) and magnetic fields, including the effects of local mass (charge) and current density. Particularly simple stability conditions are derived for source‐free regions, where the gravitational field is harmonic (∇²U = 0) or the magnetic field irrotational (∇ × B = 0). In either case the resulting stability conditions can be expressed geometrically (coordinate‐free) in terms of dimensionless stability indices. Stability bounds are calculated for several examples, including the problem of two fixed centers, the J₂ planetary model, galactic disks, and a toroidal quadrupole magnetic field. This revised version was published online in July 2006 with corrections to the Cover Date.     

Stability analysis of non-steady mhd-equilibria

Following the work of Bernsteinet al. (1958), Frieman and Rotenberg (1960) and Unno (1968) a formalism is developed which allows to examine the adiabatic stability of a perfectly conducting, rotating and self-gravitating plasma in non-steady equilibrium. Using this method the stability of a plasma in a dynamical phase of its evolution can be predicted. Global stability investigations are carried out which are based on a variation of the total energy of the system and, in general, lead to sufficient conditions for stability. The formalism is applied to the stability of a horizontal magnetic field in a medium stratified by a gravitational field.     

Kink mode <Emphasis Type="Italic">m</Emphasis> = 1 in a thin plasma filament with discontinuous vertical magnetic field

In this paper, we study the conditions of realization and stability of kink modes with azimuthal wave numbers m = ±1 in a cylindrical plasma filament with a twisted magnetic field and a homogeneous current along its axis. We assume that there are vertical constant magnetic fields inside and outside of the filament; the filament is surrounded by current-free plasma; and outside of its boundary, the azimuthal magnetic field decreases inversely in proportion to the distance from the filament’s border. The dispersion equations for stable and unstable modes are obtained in the approximation of “thin” plasma filament. The analysis of the equations for the case of discontinuous vertical magnetic field at the filament’s boundary is provided. The conditions of propagation of the wave modes have been defined. We have obtained that the unstable modes with m = ±1 cannot be realized. The results of this work can be applied to the interpretation of the solar magnetic flux tubes’ behavior using measurements provided by the spacecrafts.     

Non-axisymmetric magnetostatic equilibrium

This paper considers the structural properties of a sunspot-like magnetic flux tube which lacks perfect axisymmetry. The flux tube is taken to be in static equilibrium with an atmosphere in a uniform gravity. Assuming the departure from axisymmetry to be slight, the equations for the first order non-axisymmetric part of the equilibrium are derived in cylindrical coordinates. These first order equations reduce to a linear second order hyperbolic partial differential equation in the r-z plane. Whereas Cauchy type boundary conditions are appropriate for hyperbolic equations, physical considerations dictate the specification of boundary conditions on a closed surve for our problem of interest. The construction of solutions to this boundary value problem is illustrated with three analytically soluble cases, where the zero-order axisymmetric equilibria are chosen to have magnetic field geometry of different complexity. A physical discussion of the results is given.     

A numerical study of the pre-ejection, magnetically-sheared corona as a free boundary problem

A class of magnetostatic equilibria with axial symmetry outside a unit sphere in the presence of plasma pressure and an r ^–2 gravitational field is constructed. The structure contains a localized current-carrying region confined by a background bipolar potential field, and the shape of the region changes subject to the variation of the electric current. The continuity requirement for the magnetic field and plasma pressures at the outer boundary of the cavity defines a free boundary problem, which is solved numerically using a spectral boundary scheme. The model is then used to study the expansion of the current-carrying region, caused by the buildup of magnetic shear, against the background confining field. The magnetic shear in our model is induced by the loading of an azimuthal field, accompanied by a depletion of plasma density.We show that due to the additional effect of confinement by the dense surrounding plasma, the energy of the magnetic field can exceed the energy of its associated open field, presumably a necessary condition for the onset of coronal mass ejections. (However, the plasma beta of the confining fluid is higher than that in the outer boundary of a realistic helmet-streamer structure.) Furthermore, under the assumption that coronal mass ejections are driven by magnetic buoyancy, the result from our model study lends further support to the notion of a suspended magnetic flux rope in the low-density cavity of a helmet-streamer as a promising pre-ejection configuration.The National Center for Atmospheric Research is sponsored by the National Science Foundation.