Incompressible magnetohydrodynamic modes in the thin magnetically twisted flux tube
^{1} Department of Space Plasma, Space Research Institute, 03056 Kiev, Ukraine
email: oleg.cheremnykh@gmail.com
^{2} Department of Automatic Control and Systems Engineering, University of Sheffield, Sheffield, S1 3JD, UK
email: v.fedun@sheffield.ac.uk
^{3} School of Mathematics and Statistics, University of Sheffield, Sheffield, S3 7RH, UK
email: g.verth@sheffield.ac.uk
Received: 7 October 2016
Accepted: 13 April 2017
Context. Observations have shown that twisted magnetic fields naturally occur, and indeed are omnipresent in the Sun’s atmosphere. It is therefore of great theoretical interest in solar atmospheric waves research to investigate the types of magnetohydrodynamic (MHD) wave modes that can propagate along twisted magnetic flux tubes.
Aims. Within the framework of ideal MHD, the main aim of this work is to investigate small amplitude incompressible wave modes of twisted magnetic flux tubes with m ≥ 1. The axial magnetic field strength inside and outside the tube will be allowed to vary, to ensure the results will not be restricted to only cold plasma equilibria conditions.
Methods. The dispersion equation for these incompressible linear MHD wave modes was derived analytically by implementing the long wavelength approximation.
Results. It is shown, in the long wavelength limit, that both the frequency and radial velocity profile of the m = 1 kink mode are completely unaffected by the choice of internal background magnetic twist. However, fluting modes with m ≥ 2 are sensitive to the particular radial profile of magnetic twist chosen. Furthermore, due to background twist, a low frequency cutoff is introduced for fluting modes that is not present for kink modes. From an observational point of view, although magnetic twist does not affect the propagation of long wavelength kink modes, for fluting modes it will either work for or against the propagation, depending on the direction of wave travel relative to the sign of the background twist.
Key words: Sun: atmosphere / Sun: general / Sun: oscillations
© ESO, 2017
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0),
which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
An axially symmetric, vertical and magnetically twisted flux tube is a convenient model for analytical studies of various magnetohydrodynamic (MHD) perturbations. For a long time this approximation was the focus of investigations of MHD wave propagation in solar and space plasmas (see e.g. Priest & Hood 1991; Bennett et al. 1999; Erdélyi & Fedun 2007; Ruderman 2007; LadikovRoev et al. 2013; Cheremnykh et al. 2014) and hightemperature (see e.g. Suydam 1958; Bateman 1978; Galeev & Sudan 1989; Cheremnykh 1989; Burdo et al. 1994) plasmas. This geometry is also a useful approximation in solving fundamental problems of plasma physics (see for example Trubnikov 1966; Filippov 2007; Cheremnykh et al. 1994; Zagorodny & Cheremnykh 2014), to name but a few. In spite of many previous theoretical studies of wave propagation in solar magnetic flux tubes many questions still remain open. There are at least two contradictory opinions on how the radial dependence of the equilibrium azimuthal component of the magnetic field outside of the flux tube should be modelled. Filippov (2007) and Vršnak (2008) assume that external magnetic field decreases with distance from the tube boundary inversely proportional to the radius, that is, as a function of 1 /r. This approximation was previously used by, for example, Erdélyi & Fedun (2006), Ruderman (2015), Giagkiozis et al. (2015). Recently, Giagkiozis et al. (2015) has shown that the wave solution for a background constant twist outside the tube is actually very close to the solution when the twist is proportional to 1 /r. From another point of view, by taking into account plasma conductivity, the external magnetic field does not penetrate significantly through the tube boundary and, therefore, can be neglected (see e.g. Parker 2007; Solov’ev 2011, 2012). In these papers authors applied simple electromagnetic assumptions confirming the absence of azimuthal components of the magnetic field outside the tube.
In the present work, we will examine MHD wave propagation in a magnetic flux tube with an internal twist only. To go beyond cold plasma equilibria conditions the axially aligned magnetic field inside and outside the flux tube are allowed to be different. A similar background model has been used previously by Bennett et al. (1999) and thereafter in a number of other papers, for example Erdélyi & Fedun (2007, 2010). In the framework of ideal MHD, we will assume incompressible linear perturbations and implement the thin tube approximation. Also, we will focus on the analytical solutions related to modes with only m ≥ 1. The analytical dispersion relation and expression for eigenfunctions will be obtained by assuming the small parameter ε = k_{z}a ≪ 1, where a is the radius of the magnetic flux tube and k_{z} is the longitudinal wavenumber.
2. Derivation of linear radial component wave equation
We proceed from the linearised ideal MHD equations for the displacement vector ξ of a finite volume element. By assuming the time dependence of all perturbed physical quantities as exp(−iωt), these equations can be written as (see e.g. Priest 2014; Kadomtsev 1966; Zagorodny & Cheremnykh 2014) (1)where
(2)Here, the symbol δ corresponds to the perturbed quantities, ρ is the equilibrium plasma density, p is the equilibrium plasma pressure, ξ = ξ_{r}e_{r} + ξ_{ϕ}e_{ϕ} + ξ_{z}e_{z} is the displacement vector (where e_{r}, e_{ϕ} and e_{z} are unit vectors of cylindrical coordinates r, ϕ and z), γ is the adiabatic index, ω is the angular frequency, B is the equilibrium magnetic field, τ is the normalised magnetic field, δp_{1} is the perturbation of total plasma pressure (plasma plus magnetic), and κ is the vector of curvature of magnetic field lines. The derivation of the expression for δp_{1} is shown in Appendix A. ξ_{⊥} = ξ−ξ_{∥}τ in Eq. (2) corresponds to the perpendicular component of ξ to the equilibrium magnetic field. All physical quantities depend on the radial coordinate r. For simplicity, the magnetic field is normalised as . We use the cylindrical coordinate system (r,ϕ,z), hence the magnetic surfaces are nested cylinders of radius r and the unperturbed tube axis is parallel to the zaxis (see Fig. 1). Both inside and outside the tube the equilibrium magnetic field is given by (3)which must satisfy the magnetohydrostatic equilibrium: (4)Using the relation and the rules for differentiating unit vectors, ∂e_{r}/∂ϕ = e_{ϕ} and ∂e_{ϕ}/∂ϕ = −e_{r}, Eqs. (1) and (2) reduce to:
To obtain Eqs. (5)−(7) we assumed that all equilibrium quantities depend on r alone. Therefore, we can Fourier decompose the solution as where m is the azimuthal wave number. In this study we shall not deal with the axisymmetric case when m = 0, which describes the sausage mode. Instead, we focus on all the modes with m ≥ 1, corresponding to nonaxially symmetric oscillations which are the kink m = 1 and surface m > 1 modes. A cartoon of the model geometry is shown in Fig. 1.
Fig. 1 Magnetohydrostatic equilibrium of the twisted magnetic flux tube. 

Open with DEXTER 
Since the displacement vector ξ(r) is a function of the radius alone, then the problem becomes onedimensional. For more convenient analysis of the perturbations in the Eqs. (5)−(7) we changed the ϕ and z components of the displacement vector ξ and wave vector k to the components directed along the binormal (subscript b) and along the magnetic field lines (subscript ∥):
(8)In these more convenient variables, the components of the perturbed magnetic field can be obtained from Eq. (2):
(9)and Eqs. (5)−(7) can be written as
where the values of the total pressure δp_{1} and ∇·ξ are also represented via radial (ξ_{r}), binormal (ξ_{b}) and field aligned (ξ_{∥}) displacements:
(13)Equations (10)−(13) are the starting point for further analysis. Let us reduce the Eqs. (10)−(12) to the one single equation for the radial component of the displacement vector (ξ_{r}). From Eqs. (11, 12) one can obtain
Here, By substituting Eqs. (14) and (15) into Eq. (13), we can obtain
(16)where From Eqs. (4), (10), (14) and (16) we obtain the governing wave equation for the linear radial component (ξ_{r}): (17)In Cheremnykh (2015) it was shown that Eq. (17) is equivalent to the well known HainLüst equation (see e.g. Hain & Lüst 1958) and also it was shown that from this equation we can obtain Suydam’s criterion (see e.g. Suydam 1958) and stability criterion for ballooning modes.
3. Incompressible perturbations in the long wavelength approximation
Let us assume that plasma perturbation is incompressible, that is the velocity perturbation δv = ∂ξ/∂t is very small relative to the sound speed, (γ → ∞). It follows from the second equation of (13) and Eqs. (14), (15) that (18)For an incompressible perturbation the expression (16) for δp_{1} becomes (19)By substituting Eq. (19) into Eq. (18), we obtain the condition for the divergence of an incompressible flow, that is, (20)Equation (20) can also be obtained in another way. By substituting δp_{1} from Eq. (13) into Eq. (12), we obtain . By assuming , we also arrive at the same Eq. (20).
For incompressible perturbations Eq. (17) reduces to the equation (21)and coincides with the Eq. (14.36) in Miyamoto (2005). Equation (21) has only one singular point at Alfvén frequency, when . To obtain this equation we assumed that k_{∥}c_{S} → ∞ and, therefore, the longitudinal component wave vector k_{∥} cannot vanish.
For perturbations with a small azimuthal wave number m, Eq. (21) can be simplified further by implementing the long wavelength approximation, that is, for the case when ε = k_{z}a ≪ 1. In this limit, for a homogeneous longitudinal magnetic field, B_{z} = const., Eq. (21) is given by
(22)where To obtain Eq. (22) we also used the relevant long wavelength approximation that Further analysis of the Eq. (22) requires an equilibrium model for the flux tube and the boundary conditions for the perturbed quantities on its surface.
4. Boundary conditions
To obtain the first boundary condition we used the equation of incompressibility given by Eq. (20) and relation Eq. (13). By following the methodology presented in Jackson (1998) and Priest & Hood (1991), let us integrate Eq. (20) over a small interval 2δ near flux tube surface (i.e, within the layer a−δ and a + δ, a is the tube radius, δ ≪ a) results in
By taking into account the continuity of the integrated function, this equation can be represented as Assuming that δ → 0, the first boundary condition can be represented in the form: (23)This condition requires that the radial plasma displacement is continuous at r = a. Now let us obtain the second boundary condition. Equation (22) can be rewritten as (24)By taking into account that the magnetic field and plasma equilibrium parameters are different inside and outside of the magnetic flux tube, Eq. (24) has different solutions for r>a and r<a. The inner and outer solutions should agree for values of r in an intermediate region if δ → 0. By assuming that radial plasma displacement ξ_{r} is continuous, the matching condition at the boundary (see e.g. Soloviev 1975; Priest & Hood 1991) can be obtained by integrating Eq. (24) between a−δ and a + δ: (25)For the function φ inside and outside the flux tube, Eq. (25) represents dispersion the relation for MHD oscillations. According to Eq. (19) the total pressure perturbation is given by This means that Eq. (25) is equivalent to (26)Equation (26) is the dynamic boundary condition that is often applied to study MHD perturbations of magnetic flux tubes (Erdélyi & Fedun 2006; Bennett et al. 1999; Erdélyi & Fedun 2007, to name but a few).
5. General dispersion relation for m ≥ 1 modes
In this section, Eq. (22) together with the boundary conditions Eqs. (23) and (26) will be used to obtain the dispersion equation governing oscillations of the magnetic tube. In the long wavelength approximation in Eq (22) the last term is proportional to ε^{2} and hence a good approximation can be neglected, resulting in (27)This equation has been obtained previously by Wesson (1978) in the study of stability of hightemperature plasma. Changing back to the physical variable ξ_{r}, it is easy to show that Eq. (27) is equivalent to (28)For the specific kink mode value of m = 1, from Eq. (28) we obtain: (29)Assuming zero background magnetic twist outside the tube but an arbitrary twist inside: (30)The physical solution of Eq. (29), for a trapped mode with background magnetic field Eq. (30) bounded at r = 0 and tending to ξ_{r} = 0 as r → ∞ is given as (31)Physically, the radial displacement inside the tube is restricted to be constant with the radius for the m = 1 mode (see Eq. (31)) or else Eq. (29) would give a singularity at r = 0. Mathematically, this singularity can be eliminated only in the case where B_{ϕ} is not finite at the tube axis which is unphysical. The fact that Eq. (31) is the same as for an untwisted tube is a new and interesting analytical result since previously Ruderman (2007), for example, only demonstrated this for the particular internal background magnetic twist of B_{ϕ} ∝ r. Here we have shown that the radial profile of ξ_{r} the kink mode, in the thin tube approximation, is completely independent of any prescribed background internal twist. This is also in agreement with the purely numerical study of Terradas & Goossens (2012) who solved the ideal linearised MHD equations in the zeroβ regime using the PDE2D code (Sewell 2005). Terradas & Goossens found that the kink mode frequency in the long wavelength approximation was not affected by their particular choice of a quadratic radial profile of B_{ϕ} shown in Eq. (8) of their paper, which is consistent with our more general analytical result.
By substituting Eq. (31) into the boundary condition Eq. (25) we obtain: (32)Here, subscripts i and e correspond to the internal and external parts of the magnetic flux tube correspondingly. By taking into account that
(33)from Eq. (32) we obtain the dispersion relation (34)which again, is the same as the kink mode in the thin tube approximation without twist. Therefore, by Eqs. (34) and (31), both the frequency and ξ_{r} eigenfunctions are unaffected by the choice of internal twist for the kink mode.
By inspection of Eq. (28) it can be seen that the case is entirely different for modes with m ≥ 2. In such cases the dispersion relation and ξ_{r} eigenfunction will depend substantially on the radial profile of magnetic twist. However, unlike the m = 1 kink mode, Eq. (28) is less tractable for analysis when the radial profile of B_{ϕ} is arbitrary for the modes with m ≥ 2. Hence, in the next section we shall choose a specific radial twist profile which will enable us to do this.
6. Internal background magnetic twist with B_{ϕ} ∝ r
To make analysis of Eq. (28) more straightforward for fluting modes with m ≥ 2 we will assume that inside the tube the magnetic twist varies linearly and outside it is zero, for example Bennett et al. (1999), Parker (2007), Solov’ev (2011, 2012): (35)where B_{zi} and B_{ze} are constant internal and external magnetic fields. We note that after introduction of this specific magnetic geometry (i.e. Eq. (35)), the quantity loses its dependence on r. From Eqs. (27) and (35) we obtained that inside and outside the flux, tube Eq. (28) can be represented as (36)By applying the boundary conditions given by Eq. (23) the solution of Eq. (36), which describes perturbation of the plasma cylinder border for m ≥ 1, is (37)Here ξ_{a} = ξ_{r}(r = a). By substituting Eq. (37) into Eq. (25) we obtain the dispersion relation in the form:
(38)This dispersion relation is valid for positive azimuthal m and longitudinal k_{z} wave numbers. It is important to note that this relation is invariant under the substitution (m,k_{z}) → (−m,−k_{z}), resulting in:
(39)Hence, Eq. (39) describes wave propagation in both directions along the twisted magnetic tube. If current is absent, meaning that B_{ϕi} = 0, the two last terms vanish and Eq. (39) describes the kink mode (see e.g. Ryutov & Ryutova 1976; Spruit 1982; Edwin & Roberts 1983). Although m in Eq (39) can either be positive or negative, in the following sections we analyse only the case when m> 0. (40)where ξ_{a} = ξ_{r}^{(}r = a^{)} = φ_{a}/r. This function is constant up to the boundary of the flux tube and then decreasing to the infinitely small values as r → ∞ (see Fig. 2).
Fig. 2 Radially dependent eigenfunctions ξ_{r}(r) for axially symmetric magnetic flux tube with a free boundary are shown for azimuthal wave numbers m = 1, 2, 3, and 4. 

Open with DEXTER 
The square of the frequency of fluting modes (m ≥ 2) given by Eq. (38), in contrast to kink modes (m = 1), has a minimum value when and is equal to (41)In contrast to the kink mode, Eq. (41) shows that a low frequency cutoff is introduced for fluting modes due to background twist. For all m ≥ 2 modes the eigenfunction ξ_{r} has a form of power function Eq. (37) and describes perturbations which are localised at the surface of the twisted magnetic flux tube (see Fig. 2). It is also interesting to note that the sign of the ξ_{r} eigenfunctions ξ_{r} do not vary as r increases.
7. Comparison with previous results
Equation (35) describes the same background magnetic field configuration as Bennett et al. (1999), to allow for direct comparison. In the paper by Bennett et al. the authors obtain the specific dispersion relation which describes sausage (m = 0) waves in an incompressible magnetic flux tube with uniform twist embedded within an untwisted magnetic environment, that is:
(42)Here I_{m}(x) and K_{m}(x) are modified Bessel functions of the first and second kind with an imaginary argument, (43)and (44)The subscript α refers to i and e, which corresponds to the quantities either inside and outside of the tube, respectively. We note that Bennett et al. (1999) have shown numerically (see. Fig. 5 of their paper) that for the kink mode (m = 1), the phase and group velocities are approximately equal, such that: (45)which is in a great agreement with Eq. (34). We now want to check that the dispersion relations given by Eqs. (38) and (42) are equivalent to each other in the long wavelength approximation when m = 0 and furthermore, that the eigenfunctions are also the same. As mentioned in Bennett et al. (1999), it is difficult to analyse analytically the modes with m ≥ 1 and hence this was not attempted in their work. This is partly due to the fact that modes with m ≥ 1 can exist with both and . Equation (42) was obtained for the particular case of . Let us show that these conditions do not actually affect the final form of the dispersion equation in the long wavelength approximation. In the case in which only one change in the left hand side of Eq. (42) will appear due to the replacement (see Appendix B, Eq. (B.8)) (46)Here J_{m}(x) is the Bessel function of the first kind. In the longwave approximation m_{0}a ≈ k_{z}a ≪ 1, therefore the Bessel functions satisfy following relations: (47)By applying relations shown in Eq. (47) we arrive to the conclusion that, independent of the sign of , Eq. (42) can be written in the form: (48)After some algebra it can be shown that Eq. (48) coincides with the previously obtained dispersion relation in Eq. (38).
We now wish to further confirm our results by comparing the eigenfunctions of Eqs. (37) and (40) with those of Bennett et al. (1999). Bennett et al. obtained a solution (see Eq. (B.6)) for the total pressure perturbation δp_{1} when : (49)If , the expression for δp_{1} has the following form: (50)Since in the long wavelength approximation, functions I_{m}(x), J_{m}(x) and K_{m}(x) satisfy relations: (51)δp_{1}, independent of the sign of , is given by: (52)According to Eqs. (B.1), (B.4) and (B.5) the total pressure perturbation δp_{1} is associated with radial displacement ξ_{r} as: (53)From Eqs. (52) and (53) the radial displacement ξ_{r} depends on r as follows: (54)This is consistent with the radial dependence of the eigenfunctions previously derived from dispersion relation (38) and shown in Eqs. (37) and (40).
Now we shall compare our results with those obtained by Ruderman (2007) who studied linear nonaxisymmetric oscillations of a thin magnetic tube in presence of the weak internal magnetic twist, that is B_{ϕ}(r ≤ a) ≪ B_{zi} and B_{ϕ}(r>a) = 0. Ruderman obtained the dispersion relation for incompressible MHD perturbations in the form of (55)It can be seen that this equation is the same as Eq. (39) when B_{ϕ}(a) /a = A and B_{zi} = B_{ze} = B_{0}.
Therefore, the dispersion relation of Eq. (39) and the resulting expressions for the eigenfunctions in Eqs. (37) and (40) are in excellent agreement with the previous results of Bennett et al. (1999), who studied more specialised cases of our more general plasma and magnetic field background configuration.
8. Conclusions
In this work we analysed the incompressible linear MHD modes of a twisted magnetic flux tube. Special attention was given to the problem of finding the eigenvalues and eigenfunctions of modes with m ≥ 1 in the long wavelength limit. The description of these modes leads to significant mathematical difficulties as noted in Bennett et al. (1999). However, in this current work, by using the long wavelength approximation we have made significant analytical progress. Notably, we have shown that the dispersion relation given in Eq. (34) for the m = 1 kink mode is completely unaffected by the radial profile of background internal magnetic twist. It was also found from the ξ_{r} eigenfunction that the kink mode is bodylike in character (see Eq. (40)) and the higher order modes (m ≥ 2) are surfacelike (see Eq. (37)). Hence, these results are the same as known previously in the case of no background magnetic twist for the particular choice of having B_{ϕ} ∝ r inside the tube. However, it was shown in Eq. (41), that a low frequency cutoff was introduced for fluting modes due to the presence of background twist, in contrast to the cutoff free propagation for the kink mode.
It can be seen by the derived dispersion relation shown in Eq. (34) that the phase speed of the kink mode, equal to the group speed in the long wavelength limit, will not be affected at all by the presence of internal background magnetic twist. When there is no internal twist, all surface modes tend to the kink speed in the long wavelength limit. However, this is not the case when twist is present, since it can either work for or against the propagation speed of these modes. For example, when m, k_{z} and B_{ϕ} are all the same sign, the resulting phase speed is increased relative to the kink speed. If the sign of B_{ϕ} is opposite to that of m and k_{z}, then the speed is reduced relative to the kink speed. This is analogous to MHD wave mode propagation along a magnetic flux tube in the presence of fieldaligned flow. That is to say, if the wave is travelling in the same or opposite direction to that of the flow, the speed is increased/decreased relative to the case when no flow is present. Hence, this presents a challenge in interpreting what differences in observed counterpropagating MHD wavemode speeds could be caused by, meaning a magnetic twist, fieldaligned flow, or a combination of both. Futhermore, to our knowledge m ≥ 2 modes have still yet to be identified in solar atmospheric observations of thin twisted magnetic flux tubes, for example, chromospheric fibrils, mottles, and spicules. However, the current work suggests their sensitivity to magnetic twist would make them a very interesting future case study.
Acknowledgments
All authors thank the referee and Prof. M. Ruderman for the constructive suggestions that improved the paper. O.C. and A.K. would like to thank the Ukrainian Scientific and Technical Center, PN 6060; Integrated ScientificProgrammes of the National Academy of Science of Ukraine on Space Research and Plasma Physics. V.F., G.V. thank the STFC and Royal SocietyNewton Mobility Grant for the support received. V.F. also thanks the Newton Fund MAS/CONACyT Mobility Grants Program.
References
 Bateman, G. 1978, MHD instabilities (Cambridge, Mass.: MIT Press) [Google Scholar]
 Bennett, K., Roberts, B., & Narain, U. 1999, Sol. Phys., 185, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Burdo, O. S., Cheremnykh, O. K., Revenchuk, S. M., & Pustovitov, V. D. 1994, Plasma Phys. Control. Fusion, 36, 641 [NASA ADS] [CrossRef] [Google Scholar]
 Cheremnykh, O. K. 1989, Nucl. Fusion, 29, 1899 [CrossRef] [Google Scholar]
 Cheremnykh, O. K. 2015, Kinemat. Phys. Celest. Bodies, 31, 213 [NASA ADS] [CrossRef] [Google Scholar]
 Cheremnykh, O. K., Andrushchenko, Z. M., Edenstrasser, J. W., & Taranov, V. B. 1994, Phys. Plasmas, 1, 2525 [NASA ADS] [CrossRef] [Google Scholar]
 Cheremnykh, O. K., Klimushkin, D. Yu., & Kostarev, D. V. 2014, Kinemat. Phys. Celest. Bodies, 30, 209 [NASA ADS] [CrossRef] [Google Scholar]
 Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Erdélyi, R., & Fedun, V. 2006, Sol. Phys., 238, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Erdélyi, R., & Fedun, V. 2007, Sol. Phys., 246, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Erdélyi, R., & Fedun, V. 2010, Sol. Phys., 263, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Filippov, B. P. 2007, Eruptive process on the Sun (Moscow: Fizmatlit), 216 [in Russian] [Google Scholar]
 Galeev, A. A., & Sudan, R. N. 1989, Basic plasma physics. Handbook of plasma physics (Amsterdam: NorthHolland) [Google Scholar]
 Giagkiozis, I., Fedun, V., Erdélyi, R., & Verth, G. 2015, ApJ, 810, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Hain, K., & Lüst, R. 1958, Z. Naturforsh. A., 13, 936 [NASA ADS] [Google Scholar]
 Jackson, J. D. 1998, Classical Electrodynamics, 3rd edn. (John Wiley & Sons Ltd) [Google Scholar]
 LadikovRoev, Y. P., Cheremnykh, S. O., Yatsenko, V. A. 2013, J. Automat. Info. Sci., 45, 48 [CrossRef] [Google Scholar]
 Kadomtsev, B. B. 1966, Rev. Plasma Phys., 2, 153 [NASA ADS] [Google Scholar]
 Miyamoto, K. 2005, Plasma Physics and Controlled Nuclear Fusion, Springer Series on Atomic, Optical, and Plasma Physics, 38 (Berlin, Heidelberg: SpringerVerlag) [Google Scholar]
 Parker, E. N. 2007, Conversations on Electric and Magnetic Fields in the Cosmos (Princeton: Princeton University Press) [Google Scholar]
 Priest, E. R. 1982, Solar Magnetohydrodynamics (Dordrecht: D. Reidel Pub. Co.) [Google Scholar]
 Priest, E. 2014, Magnetohydrodynamics of the Sun (Cambridge, UK: Cambridge University Press) [Google Scholar]
 Priest, E. R., & Hood, A. W. 1991, Advances in Solar System Magnetohydrodynamics (CUP) [Google Scholar]
 Ruderman, M. S. 2007, Sol. Phys., 246, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. S. 2015, A&A, 575, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ryutov, D. A., & Ryutova, M. P. 1976, Sov. J. Exp. Theor. Phys., 43, 491 [NASA ADS] [Google Scholar]
 Sewell, G. 2005, The Numerical Solution of Ordinary and Partial Differential Equations, 2nd edn. (New Jersey: Wiley) [Google Scholar]
 Solov’ev, A. A. 2011, Astron. Rep., 55, 1025 [NASA ADS] [CrossRef] [Google Scholar]
 Solov’ev, A. A. 2012, Astrophys. Space Sci. Proc., 30, 203 [NASA ADS] [CrossRef] [Google Scholar]
 Soloviev, L. S. 1975, Rev. Plasma Phys., 6, 239 [NASA ADS] [Google Scholar]
 Spruit, H. C. 1982, Sol. Phys., 75, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Suydam, B. R. 1958, in Proc. of the Second UN Int. Conf. on the Peaceful Uses of Atomic Energy, Geneva, 31, 157 [Google Scholar]
 Terradas, J., & Goossens, M. 2012, A&A, 584, 112 [CrossRef] [EDP Sciences] [Google Scholar]
 Trubnikov, B. A. 1966, Plasma theory (Energoatomizdat) [in Russian] [Google Scholar]
 Vršnak, B. 2008, Annales Geophysicae, 26, 3089 [NASA ADS] [CrossRef] [Google Scholar]
 Wesson, J. A. 1978, Nucl. Fusion, 18, 87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zagorodny, A. G., & Cheremnykh, O. K. 2014, Introduction to the plasma physics, 696 (Kiev: Naukova Dumka) [in Russian] [Google Scholar]
Appendix A: Obtaining an expression for δp_{1}
From the following equations:
(A.1)By using the equation of equilibrium (A.2)after some algebra from Eq. (A.1) we obtain (A.3)By using the vector relations
(A.4)where τ = B/B we can rewrite Eq. (A.3) as
(A.5)The first term on the RHS of Eq. (A.5) can be transformed as follows:
(A.6)Here κ = ^{(}τ· ▽^{)}τ is the vector of curvature of the magnetic field lines. To obtain Eq. (A.6) we take into account that displacement vector ξ is represented in the form ξ = ξ_{⊥} + ξ_{∥}τ, where subscripts ⊥ and ∥ correspond to the perpendicular and parallel components to the equilibrium magnetic field. Therefore, Eq. (A.5) can be rewritten in the form:
(A.7)From the second and third terms in Eq. (A.7) we can obtain
(A.8)From Eqs. (A.7) and (A.8) it follows that:
(A.9)The first term on the RHS of Eq (A.9) is equal to:
(A.10)Finally, we obtain the second expression in Eq. (2), that is, (A.11)
Appendix B: Exact dispersion relation for the magnetic field described as Eq. (35) for
To obtain the exact dispersion relation for the case (35) we have used the following set of Eqs: (13) and (19), (20) and (21). From Eqs. (19) and (21) we have:
(B.1)Equation (B.1) together with Eqs. (19) and (13) represent a system of equations with respect to ξ_{r} and δp_{1}. From this system it is easy to obtain the equation for δp_{1}:
(B.3)From Eq. (B.3) it follows that, for the magnetic field (35), the following relations are satisfied: (B.4)Therefore, (B.5)By taking into account Eqs. (B.4) and Eq. (B.5), from (B.2) we obtain
(B.6)Equation (B.6) coincides with the Eq. (13) of work Bennett et al. (1999) and satisfies the boundary conditions (23) and (25).
Following Bennett et al. (1999), we introduce (see (43)). We assume , that gives a solution (B.6) in the form (B.7)To obtain Eq. (B.7) we take into account that m_{0} = k_{z} for r>a. From Eqs. (B.1), (23), (25) and (B.7), we obtain desired dispersion relation: (B.8)Used in Eq. (B.8) variables and parameters are explained in Eq. (43).
All Figures
Fig. 1 Magnetohydrostatic equilibrium of the twisted magnetic flux tube. 

Open with DEXTER  
In the text 
Fig. 2 Radially dependent eigenfunctions ξ_{r}(r) for axially symmetric magnetic flux tube with a free boundary are shown for azimuthal wave numbers m = 1, 2, 3, and 4. 

Open with DEXTER  
In the text 