Modeling Kelvin – Helmholtz Instability in Soft X-Ray Solar Jets

Development of Kelvin–Helmholtz (KH) instability in solar coronal jets can trigger the wave turbulence considered as one of the main mechanisms of coronal heating. In this review, we have investigated the propagation of normal MHD modes running on three X-ray jets modeling them as untwisted and slightly twisted moving cylindrical flux tubes. The basic physical parameters of the jets are temperatures in the range of 5.2–8.2MK, particle number densities of the order of 109 cm, and speeds of 385, 437, and 532 km s, respectively. For small density contrast between the environment and a given jet, as well as at ambient coronal temperature of 2.0MK and magnetic field around 7G, we have obtained that the kink (m = 1) mode propagating on moving untwisted flux tubes can become unstable in the first and second jets at flow speeds of ≅348 and 429 km s, respectively. The KH instability onset in the third jet requires a speed of ≅826 km s, higher than the observed one. The same mode, propagating in weakly twisted flux tubes, becomes unstable at flow speeds of ≅361 km s for the first and of 443 km s for the second jet. Except the kink mode, the twisted moving flux tube supports the propagation of higher (m > 1) MHDmodes that can become unstable at accessible jets’ speeds.


Introduction
Jets are considered to be ubiquitous confined plasma ejecta in the solar atmosphere.They have been extensively observed in the solar atmosphere in various wavebands, such as H [1,2], Ca ii H [3][4][5], EUV [6], and soft X-ray [7] in order to understand their multitemperature characteristics.X-ray jets were discovered by the Soft X-Ray Telescope (SXT) on board Yohkoh [8], as transient X-ray energy release and enhancement with apparent collimated ballistic motions of the plasma associated with the flares in X-ray bright points, emerging flux regions, or active regions (for details, see Shibata et al. [7]).As it has been pointed out by Shimojo et al. [9], jets from X-ray bright points in active regions most likely appear at the western edge of preceding sunspots and exhibit a recurrent plasma propulsion in the solar atmosphere.These X-ray jets are confined plasma dynamics with typical morphological properties, for example, (1-40) × 10 4 km length, and the width of 5 × 10 3 -10 5 km.Such jets possess apparent velocities of 10-1000 km s −1 and lifetime of 100-16,000 s [9].The electron densities of the X-ray jets are of the orders of (0.7-4) × 10 9 cm −3 .Their temperatures lie in the range of 3-8 MK with an average temperature of 5.6 MK [10].
In terms of spatial location jets can be classified as polar jets [11] and active region jets [12].A study of polar jet parameters based on Hinode XRT observations was carried out by Savcheva et al. [13] who showed that jets preferably occur inside the polar coronal holes.Culhane et al. [11] have found from Hinode's Extreme-ultraviolet Imaging Spectrometer (EIS) 40  slot observations of a polar coronal hole that jet temperature ranges from 0.4 to 5.0 MK.The jet velocities had typical values that are mostly less than the Sun's escape velocity (618 km s −1 ); therefore, in consequence most of the jets fall back in the lower solar atmosphere after their triggering.Using the XTR on Hinode, Cirtain et al. [14] conclude that X-ray jets in polar coronal holes have two distinct velocities: one near the Alfvén speed (∼800 km s −1 ) and another near the sound speed (200 km s −1 ).Moreover, they were the first to give an evidence for the propagation of Alfvén waves in solar X-ray jets.Kim et al. [15] presented 2 Advances in Astronomy the morphological and kinematic characteristics of three small-scale X-ray/EUV jets simultaneously observed by the Hinode XRT and the Transition Region and Coronal Explorer (TRACE).While observing the coronal jets, for two different wavelength bands, they obtain matching characteristics for their projected speed (90-310 km s −1 ), lifetime (100-2000 s), and size (1.1-5 × 10 5 km).Chifor at al. [12] have reported 2007 January 15/16 observations of a recurring jet situated on the west side of NOAA active region 10938.A strong blue-shifted component and an indication of a weak redshifted component at the base of the jet were observed around  e = 1.6 MK in these jets.The upflow velocities were observed exceeding up to 150 km s −1 .These jets were seen over a range of temperatures between 0.25 and 2.5 MK, while their estimated electron densities lie above 10 11 cm −3 for the high-velocity upflow components.Yang et al. [16] presented simultaneous observations of three recurring jets in EUV and soft X-ray (SXR), which occurred in an active region on 2007 June 5.On comparing their morphological and kinematic properties, the authors have found that EUV and SXR jets had similar onset locations, directions, size, and terminal velocities.The three observed jets were having maximum Doppler velocities ranging from 25 to 121 km s −1 in the Fe xii 195 line and from 115 to 232 km s −1 in the He ii 256 line.Extensive multi-instrument observations obtained simultaneously with the SUMER spectrometer on board the Solar and Heliospheric Observatory (SoHO), with EIS and XRT on board Hinode, and with the Extreme-ultraviolet imagers (EUVI) of the Sun-Earth Connection Coronal and Heliospheric Investigation (SECCHI) instrument suite on board the Ahead and Behind STEREO spacecrafts were performed by Madjarska [17].The dynamic process of X-ray jet formation and evolution has been derived in great detail.In particular, for the first time there was found spectroscopically a temperature of 12 MK (Fe xxiii 263.76 Å) and density of 4 × 10 10 cm −3 in the quiet Sun.The author has clearly identified two types of upflows in which the first one was the collimated upflow along the open magnetic fields, and the second was the formation of a plasma cloud from the expelled bright point small-scale loops.Chandrashekhar et al. [18] studied the dynamics of two jets seen in a polar coronal hole with a combination of EIS and XRT/Hinode data.They found no evidence of helical motions in these events but detected a significant shift of the jet position in a direction normal to the jet axis, with a drift velocity of about 27 and 7 km s −1 , respectively.
The launch of the Solar Dynamics Observatory (SDO) [19] with the Atmospheric Imaging Assembly (AIA) [20,21] opens a new page in observing the solar jets.Moschou et al. [22] have reported high cadence observations of solar coronal jets observed in the extreme-ultraviolet (EUV) 304 Å using Atmospheric Imaging Assembly (AIA) instrument on board SDO.They registered, in fact, coronal hole jets, with speeds of 94 to 760 km s −1 and lifetimes of the order of several tens of minutes.A detailed description of the dynamical behavior of a jet in an on-disk coronal hole observed with AIA/SDO was presented by Chandrashekhar et al. [23].Their study reveals new evidence of plasma flows prior to the jet's initiation along the small-scale loops at the base of the jet.The authors have also found further evidence that flows along the jet consisting of multiple, quasi-periodic smallscale plasma ejections.In addition, spectroscopic analysis estimates temperature as Log 5.89 ± 0.08 K and electron densities as Log 8.75 ± 0.05 cm −3 in the observed jet.Measured properties of the registered transverse wave have provided evidence that strong damping of the wave occurred as it propagates along the jet with speeds of ∼110 km s −1 .Using the magnetoseismological inversion, observed plasma, and wave parameters, the jet's magnetic field is estimated as  = 1.21 ± 0.2 G. Recently, Sterling et al. [24] have reported high-resolution X-ray and extreme-ultraviolet observations of 20 randomly selected X-ray jets that form in coronal holes at the solar polar caps.In each jet, converse to the widely accepted emerging magnetic flux model, a miniature version of the filament eruptions that initiated coronal mass ejections drove the jet producing reconnection process.Formation of a rotating jet during the filament eruption on 2013 April 10-11 on the base of multiwavelength and multiviewpoint observations with STEREO/SECCHI/EUVI and SDO/AIA was reported by Filippov et al. [25].The confined eruption of the filament within a null-point topology, which is also known as an Eiffel tower magnetic field configuration, forms a twisted jet after magnetic reconnection near the null point.The sign of the helicity in the jet is observed the same as that of the sign of the helicity in the filament.It is noteworthy that the untwisting motion of the reconnected magnetic field lines gives rise to the accelerating plasma along the jet.
It is well established that the magnetic reconnection at different heights in the solar atmosphere plays a key role in triggering jet-like events.The magnetic reconnection between open and closed fields (standard reconnection scenario) is one of the well-known processes of the jet's occurrence [15,26].The jets emerging by this kind of mechanism are known as standard jets [3,7,27,28].Some other observational and simulation studies showed that the reconnection at the magnetic null in a fan-spine magnetic topology can also trigger jet-like events [29][30][31][32].Eruptions of small arches, filaments, and flux ropes from within this type of magnetic field configurations can be responsible for reconnection and jets' occurring.These types of jets are known as blowout jets [28,33,34].Moore et al. [35] used the full-disk He ii 304 Å movies from the Atmospheric Imaging Assembly on SDO to study the cool ( ∼ 10 5 K) component of X-ray jets observed in polar coronal holes by XRT.The AIA 304 Å movies revealed that most polar X-ray jets spin as they erupt.The authors examined 54 X-ray jets that were found in polar coronal holes in XRT movies sporadically taken during the first year of continuous operation of AIA (2010 May through 2011 April).These 54 jets were big and bright enough in the XRT images to be categorized as a standard jet or as a blowout jet.From the X-ray movies, 19 of the 54 jets appeared like standard jets, 32 appeared as blowout jets, and three were ambiguous and were not falling in any category.Moore et al. [36] have studied 14 large-scale solar coronal jets observed in Sun's pole.In EUV movies from the SDO/AIA, each jet was very similar to most X-ray and EUV jets erupting in coronal holes.However, each was exceptional in that it went higher than most of the standard coronal jets.They were detected in the outer corona beyond 2.2  ⊙ in images as observed from the Solar and Heliospheric Observatory/Large Angle Spectroscopic Coronagraph (LASCO/C2 coronagraph [37]).
Schmieder et al. [38] proposed a new model between the standard and blowout models, where magnetic reconnection occurs in the bald patches around some twisted field lines, one of whose foot points is open.Pariat et al. [39] included the magnetic field inclination and photospheric field distribution and performed another 3D numerical MHD model for the two different types of jet events: standard and blowout jets.We note also that stereoscopic studies (multiple points of view) have been carried out by the EUVI/SECCHI imagers on board the twin STEREO spacecraft to estimate the expected speed, motion, and morphology of polar coronal jets [40].To sum up, the flux emergence [41,42] and flux cancellation [43,44] are the two main triggering processes that are known to be responsible for jets' occurrence.In few observational and primarily numerical studies the wave-induced reconnection has also been suggested as a cause for the onset of jet-like events [42,[45][46][47].
We consider magnetically structured X-ray solar jets as moving cylindrical magnetic flux tubes that support the excitation/propagation of various kind of magnetohydrodynamic (MHD) oscillations and waves.While in static solar atmospheric plasma the propagating MHD modes are stable, the axial motion of the flux tubes engenders a velocity jump at the tube surface which can trigger a Kelvin-Helmholtz (KH) instability.The KH instability arises at the interface of two fluid layers that move with different speeds (see, e.g., Chandrasekhar [48])-then a strong velocity shear arises near the interface between these two fluids forming a vortex sheet.This vortex sheet becomes unstable to the spirallike perturbations at small spatial scales [49].In cylindrical geometry, when a magnetic flux tube is axially moving, such a vortex sheet is evolved near tube's boundary and it may become unstable against KH instability provided that the tube axial velocity exceeds a critical value [50].Further on, in the nonlinear stage of KH instability, this vortex sheet causes the conversion of the directed flow energy into turbulent energy making an energy cascade at smaller spatial scales [51].
The first modeling of KH instability in X-ray jets was carried out by Vasheghani Farahani et al. [75] who explored transverse wave propagation along the detected by Cirtain et al. [14] coronal hole soft X-ray jets.Vasheghani Farahani et al. analyzed analytically, in the limit of thin magnetic flux tube, the dispersion relation of the kink MHD mode and have obtained that this mode is unstable against the KH instability when the critical jet velocity is equal to 4.47V A = 3576 km s −1 (V A = 800 km s −1 is the Alfvén speed inside the jet).Numerical solving of the same dispersion relation when considering the jet and its environment as cold magnetized plasmas, carried out by Zhelyazkov [76,77], yielded a little bit lower critical flow speed for the instability onset; namely, 4.31V A = 3448 km s −1 .The lowest critical jet speed of 4.025V A = 3220 km s −1 was derived by numerically solving the wave dispersion relation without any approximations, that is, treating both media as compressible plasmas.But even the latter critical jet speed is still too high for the KH instability to be detected/observed in coronal hole soft X-ray jets.The reason for obtaining such high critical speeds is the circumstance that Vasheghani Farahani et al. [75] assumed an electron number density of the order of 10 8 cm −3 and magnetic field strength of 10 G.
Here, we study the propagation of kink and higher MHD modes in standard active region soft X-ray jets, notably jets #8, #11, and #16 of Shimojo and Shibata's set of sixteen observed flares and jets [10], and have shown that with one order higher electron densities, ∼10 9 cm −3 , and moderate magnetic field, ∼7 G, MHD modes in high-speed jets, like these ones, can become unstable against the KH instability at accessible jets speeds, except for jet #16 which requires a higher flow velocity.In the next section we list the basic physical parameters of jet #11 and the topology of magnetic fields inside and outside the moving flux tube modeling jet and also derive the wave dispersion equation for both untwisted and twisted tubes.The physical parameters of other two jets (#8 and #16) will be provided en route in the next section.The numerical solving MHD wave dispersion relations and the discussion of the conditions under which the KH instability can develop in such moving structures are presented in Section 3. The last section summarizes the main results obtained in this article and outlooks our future studies of KH instability in more complex (rotating) solar atmosphere jets.

Geometry, Magnetic Field Topology, and MHD Wave Dispersion Relations
We consider the soft X-ray jet as a straight cylinder of radius  and density  i embedded in a uniform field environment with density  e .We study the propagation of MHD waves in two magnetic configurations, notably in untwisted and twisted flux tubes.For an untwisted tube magnetic fields in both media are homogeneous and directed along the -axis of our cylindrical coordinate system (, , ): B i = (0, 0,  i ) and B e = (0, 0,  e ), respectively (see Figure 1).The magnetic field inside the twisted tube is helicoid, B i = (0,  i (),  i ()), while outside the tube the magnetic field is uniform and directed along the tube axis, B e = (0, 0,  e ).Note that we assume a magnetic fields' equilibrium with uniform twist for which the magnetic field inside the tube is B i = (0, ,  i ), where  and  i are constant.The parameter that characterizes the uniform magnetic field twist is the ratio  i ()/ i ≡ , that is,  = / i .Our frame of reference for studying the wave propagation in the jet is attached to the surrounding magnetoplasma-thus k 0 = (0, 0, V 0 ) represents the relative jet velocity, if there is any flow in the environment.The jump of the tangential velocity at the tube boundary then initiates the magnetic KH instability onset when the jump exceeds a critical value.
Before dealing with governing MHD equations, it is necessary to specify what kind of plasma each medium is (the moving tube and its environment).As seen from the Events List in [10], the electron density in jet #11 is  jet ≡  i = 2.9 × 10 9 cm −3 , the temperature is  jet ≡  i = 5.5-6.4MK, and jet speed (in their notation) is  jet = 437 km s −1 .Our choice for environment magnetic field, electron density, and temperature is  e = 6.7 G,  e = 2.6 × 10 9 cm −3 , and  e = 2.0 MK, respectively.With  i = 5.5 MK, the total pressure balance equation (equality of the sum of thermal and magnetic pressures in both media), that is, where  is the magnetic permeability of vacuum, yields the following basic sound and Alfvén speeds in the jet and its environment:  si = 275 km s −1 and V Ai = 111 km s −1 and  se = 166 km s −1 and V Ae = 286 km s −1 , respectively.Accordingly, the plasma betas in both media are  i = 7.369 and  e = 0.403-this implies that, in principle, one can treat the jet as incompressible plasma and its surrounding medium as cool magnetized plasma.When studying the MHD wave propagation in the untwisted magnetic flux tube we shall use two approaches, namely, of compressible plasmas in both media and the simplified limit of incompressible and cool plasmas-a similarity of dispersion curves patterns, obtained from corresponding dispersion relations, will eventually justify the usage of the second approach in particular in the case of twisted tube.There are two important input parameters in the modeling KH instability in moving magnetic flux tubes, which are the density contrast,  e / i ≡  = 0.896, and the ratio of the axial external and internal magnetic fields,  e / i ≡  = 2.44.For the twisted tube the second parameter has the form  twist =  e / i .
In a system of cylindrical coordinates, the equilibrium physical variables (density, fluid velocity, and pressure) are functions of the radial coordinate  only.Then, their perturbations can be Fourier-analyzed putting them proportional to exp[i(− +  +   )], where  is the angular wave frequency (that, in general, can be a complex quantity),  is the mode number (a positive or negative integer), and   is the axial wave number.We can eliminate all except two of the perturbations (perturbation  tot of the total (thermal + magnetic) pressure and the radial component   of the Lagrangian displacement ) to get the following governing equations [78]: The coefficients ,  1 ,  2 , and  3 are functions of the equilibrium variables  0 , B 0 , and k 0 and of the Dopplershifted frequency Ω = −k ⋅ k 0 and have the following forms: where Here  A is the Alfvén frequency and  c is the cusp frequency; the other notation is standard.
Eliminating   from (2), one obtains the well-known second-order ordinary differential equation [79][80][81] By means of the solutions to (5) in both media, one can find the corresponding expressions for   , and after merging these solutions, together with those for  tot , through appropriate boundary conditions at the interface  = , one can derive the dispersion relation of the normal modes propagating in the moving magnetic flux tube.

Dispersion Relation of MHD Modes in an Untwisted Flux
Tube.In an untwisted magnetic flux tube, the coefficient where The cusp frequency,  c , is usually expressed via the so-called tube speed,  T , notably  c =    T , where [82] The solutions for  tot can be written in terms of modified Bessel functions:   ( 0i ) inside the jet and   ( 0e ) in its surrounding plasma.We note that wave attenuation coefficients,  0i and  0e , in both media are calculated from (7) with replacing the sound and Alfvén speeds with the corresponding values for each medium.Recall that, in evaluating  0e , the wave frequency is not Doppler-shiftedit is simply .By expressing the Lagrangian displacements  i and  e in both media via the derivatives of corresponding Bessel functions and by applying the boundary conditions for continuity of the pressure perturbation  tot and   across the interface,  = , one obtains the dispersion relation of normal MHD modes propagating in a flowing compressible jet surrounded by a static compressible plasma [65,83,84] Due to the flowing plasma, the wave frequency is Dopplershifted inside the jet.We recall that for the kink mode ( = 1) one defines the so-called kink speed [82] which, as seen, is independent of sound speeds and characterizes the propagation of transverse perturbations.We will show that notably the kink mode can become unstable against KH instability.When the jet is considered as incompressible plasma and its environment as a cool one, the MHD wave dispersion relation (9) keeps its form, but the two attenuation coefficients  0i,e become much simpler, namely, respectively.

Dispersion Relation of MHD Normal Modes in a Twisted
Flux Tube.Inside the tube ( ⩽ ), where  i = , the quantities  B and  Ai take the forms respectively.For incompressible plasma, we redefine (without the loss of generality) the coefficients ,  1 ,  2 , and  3 by dividing them by  4 to obtain Radial displacement   is expressed through the total pressure perturbation as The solution to this equation obviously depends upon the magnetic field and density profile.With aforementioned coefficients ,  1 ,  2 , and  3 , evaluated for the jet's medium, (5) reduces to the modified Bessel equation where The solution to (15) bounded at the tube axis is where   is the modified Bessel function of order  and  i is a constant.Lagrangian displacement  i , by using ( 14) can be written as where the prime sign means a differentiation with respect the Bessel function argument.
For the cool environment with a straight-line magnetic field  e =  e and homogeneous density  e , the  1−3 and  coefficients take the form The total pressure perturbation outside the tube obeys the same Bessel equation as (15), but  2 0i is replaced by which coincides with the attenuation coefficient in the cool environment of an untwisted flux tube.The solution bounded at infinity now is where   is the modified Bessel function of order  and  e is a constant.
In this case, the Lagrangian displacement can be written as and the Alfvén frequency is simplified to Here, V Ae =  e / √  e is the Alfvén speed in the surrounding magnetized plasma.
The boundary conditions which merge the solutions of the Lagrangian displacement and total pressure perturbation inside and outside the twisted magnetic flux tube have the forms [85]  i    = =  e    = , where total pressure perturbations  tot i and  tot e are given by ( 17) and ( 21), respectively.With the help of these boundary conditions we derive the dispersion relation of the normal MHD modes propagating along a twisted magnetic flux tube with axial mass flow k 0 where, as we already mentioned, Ω = −k ⋅k 0 is the Dopplershifted wave frequency in the moving medium, This dispersion equation is similar to the dispersion equation of normal MHD modes in a twisted flux tube surrounded by incompressible plasma [65]-the only difference is that, there in (10),  e ≡  0e =   .

Numerical Solutions and Results
Firstly we shall study the dispersion characteristics of the kink ( = 1) mode in untwisted moving magnetic flux tube (in two approaches, notably (i) considering the jet and its surrounding plasma as compressible media and (ii) treating the jet as an incompressible medium while its environment is assumed to be cool plasma) and, later on, explore the same thing for the kink ( = 1) and higher ( > 1) MHD modes propagating in a twisted moving flux tube.Two dispersion relations ( 9) and ( 25) are transcendent equations in which the wave frequency, , is a complex quantity: Re() + Im(), while the axial wave number,   , is a real variable.The appearance of KH instability (i.e., Im() > 0) is determined primarily by the jet velocity and in searching for a critical/threshold value of it, we will gradually change velocity magnitude, V 0 , from zero to that critical value (and beyond).
Our numerical task is to solve the dispersion relation in complex variables, obtaining the real and imaginary parts of the wave frequency, or as is usually done, of the wave phase velocity V ph = /  , as functions of the axial wave number,   , at various magnitudes of the velocity shear between the soft X-ray jet and its environment, k 0 .
For numerical solving the wave dispersion relations, it is practically to normalize all variables and we do that normalizing the wave length,  = 2/  , to the tube radius, , that implies a dimensionless wave number   .All speeds are normalized with respect to the Alfvén speed inside the jet, V Ai .For normalizing the Alfvén speed in the ambient coronal plasma, V Ae , we need the density contrast, , and the ratio of the magnetic fields  =  e / i , to find V Ae /V Ai = / √ .The normalization of sound speeds in both media requires the specification of the reduced plasma betas, βi,e =  2  si,e /V 2 Ai,e .In the dimensionless computations, the flow speed, V 0 , will be presented by the Alfvén Mach number  A = V 0 /V Ai .

Kelvin-Helmholtz Instability in Untwisted Flux Tubes.
Among the various MHD wave spectra that exist in a static magnetic field flux tube of compressible plasma, surrounded by compressible medium, the most interesting for us is the kink-speed wave whose speed for jet #11, according to (10), is equal to The normal modes propagating in a static homogeneously magnetized flux tube can be pure surface waves, pseudosurface (body) waves, or leaky waves (see Cally [86]).The type of the wave crucially depends on the ordering of the basic speeds in both media (the flux tube and its surrounding plasma), more specifically of sound and Alfvén speeds as well as corresponding tube speeds.In our case the ordering is Since  <  and  <  Te along with  <  Te < , according to Cally's classification, the kink mode in a rest flux tube must be pure surface mode type S − + (for detail, see Cally [86]).For an S − + type wave  should lie between  and :  ⩽  ⩽ , and also  <  Te .In our normalization, the normalized wave velocity should be bracketed between 1 and √ βi = 2.478, that is, 1 ⩽ /(  V Ai ) ⩽ 2.478.The typical wave dimensionless wave velocity /(  V Ai ) is the normalized kink speed  k /V Ai = 1.9166, which, as expected, lies between 1 and 2.478.Concerning the relation between normalized wave phase velocity and external tube speed, in our case we have the opposite inequality; that is,  >  Te , that in Cally's normalization reads as 1.2805 > 0.8649 and in ours as 1.9166 > 1.2937.The reason for that unexpected change in the mutual relations between two velocities is the relatively bigger value of parameter  (=2.44) that yields a larger magnitude of  k , than, for example, in the case when  is close to 1-then the inequality  <  Te is satisfied (see, e.g., Zhelyazkov et al. [87]).
A reasonable question is how the flow will change the dispersion characteristics of the kink ( = 1) mode.Calculations show that the flow shifts upwards the kink-speed dispersion curve and splits it into two separate curves [65].(A similar duplication is observed for the tube-speed V Ti dispersion curves, too.)The evolution of the pair of kinkspeed dispersion curves can be seen in Figure 2(a)-the input parameters for solving the dispersion equation ( 9) of the kink mode ( = 1) traveling in a moving flux tube of compressible plasma surrounded by compressible coronal medium are  = 0.896,  = 2.44, βi = 6.141, and βe = 0.336.During calculations the Alfvén Mach number,  A , was varied from zero (static plasma) to values at which we obtained unstable solutions.Let us first note that at  A = 0 we get the kink-speed dispersion curve which at a very small dimensionless wave number    (=0.005) yields the value of the normalized wave phase velocity equal to 1.9168-very close to the previously estimated magnitude of 1.9166.This observation implies that our code for solving the wave dispersion relation is correct.For small Alfvén Mach numbers the pair of kink-speed modes travel with velocities  A ∓  k /V Ai [65] (in that case, their dispersion curves go practically parallel).At higher  A , however, the behavior of each curve of the pair  A ∓  k /V Ai turns out to be completely different.As seen from Figure 2(a), for  A ⩾ 3.785 both kink curves break and merge forming a family of semiclosed dispersion curves.A further increasing in  A leads to a separation of those curves in opposite directions.This behavior of the kink ( = 1) mode dispersion curves signals us that we are in a range of the Re(V ph /V Ai )-  plane where one can expect the occurrence of KH instability.A prediction of the Alfvén Mach number  A at which the instability will occur can be found from the inequality [88] that gives for the kink ( = 1) mode  A > 3.839.In our case the KH instability starts at that  A , for which the left-hand side semiclosed curve disappears-this happens at Advances in Astronomy  cr A = 3.865.The growth rates of unstable kink modes are plotted in Figure 2(b).The red curves in all diagrams denote the marginal dispersion/growth rate curves: for values of the Alfvén Mach number smaller than  cr A the kink mode is stable; otherwise it becomes unstable and the instability is of the KH type.With  cr A = 3.865 the kink mode will be unstable when the velocity of the moving flux tube is higher than 429 km s −1 -a speed, which is below than the observationally measured jet speed of 437 km s −1 .We would like to underline that all the stable kink modes in the untwisted ( = 0) moving tube modeling jet #11 are pure surface modes while the unstable ones are not-the latter becomes partly surface and partly leaky modes (their external attenuation coefficients,  0e s, are complex quantities with positive imaginary parts).This circumstance means that wave energy is radiated outward in the surrounding medium, which allows us to claim the KH instability plays a dual role: once in its nonlinear stage the instability can trigger wave turbulence and simultaneously the propagating KH-mode is radiating its energy outside.
When numerically solving (9) for the kink ( = 1) mode treating the jet's plasma as incompressible medium and its environment as cool plasma, we obtain dispersion curves' and growth rates' patterns very similar to those shown in Figure 2 (see Figure 3).Computations show, as expected, that the stable kink mode is a nonleaky surface mode, while the unstable one possesses a real internal attenuation coefficient,  0i =   , not changed by the instability, but the external one,  0e , becomes a complex quantity with positive real and negative imaginary parts.Now the threshold Alfvén Mach number is a little bit bigger and yields a critical jet speed for the instability onset of ≅442 km s −1 , being with 5 km s −1 higher than the observationally measured 437 km s −1 .We note that, as a rule, jet's incompressible plasma approximation yields slightly higher threshold Alfvén Mach numbers than the model of compressible media.This ascertainment shows how sensitive to jet's and its environment's media treatment the occurrence of KH instability of the kink mode in our jet is -one can expect instability onset at an accessible jet speed only if both media are considered as compressible magnetized plasmas.
The basic physical parameters of jet #8, namely  i = 4.1 × 10 9 cm −3 and  i = 5.2 MK, for a low density contrast of  = 0.963 assuming as before that the temperature of surrounding plasma is  e = 2.0 MK and the background magnetic field,  e , is equal to 7 G, yield the following sound and Alfvén speeds:  se = 166 km s −1 ,  si = 267.5 km s −1 , V Ae = 242.8km s −1 , and V Ai = 52.22km s −1 .Accordingly, plasma betas in both media are  i = 31.5 and  e = 0.56, while the magnetic fields ratio  = 4.563 defines  i = 1.53 G.The speed ordering, V Ai <  se < V Ae <  si , according to Cally's [86] classification, tells us that in a rest magnetic flux tube the propagating wave must be a pure surface mode of type S − + .The numerical computations confirm this, as well as reproduce the normalized kink speed of 3.3342 within three places after the decimal point.With the flow inclusion, the behavior of the pair kink-speed dispersion curves in the region of the instability onset, as seen from Figure 4(a), is completely different; more specifically, the lower semiclosed kink-speed dispersion curves at  A = 6.625 and 6.65 (the orange curve) correspond to pseudosurface (body) waves, while the higher kink-speed dispersion curves are associated with surface waves.Note that at  A = 6.668 the surface wave dispersion curves brakes in two parts: the right-hand one merges with the lower kink-speed curve at    = 0.2337 (see the blue line) while its left-hand side forms a narrow semiclosed dispersion curve of pure surface mode.The threshold Alfvén Mach number is equal to 6.6704 (the prediction one according to (30) is 6.67) which implies that the critical flow speed for KH instability onset is equal to 348.3 km s −1 -a value far below the jet speed of 385 km s −1 .This circumstance guarantees that even if one considers the jet medium as incompressible plasma and its environment as cool plasma we will get an accessible jet's speed for instability onset.
In a similar way, calculating the sound and Alfvén speeds for jet #16 and its environment (with  i = 1.0 × 10 9 cm −3 ,  i = 7.4 MK,  = 0.95,  e = 2.0 MK, and  e = 6.7 G) we get  se = 166 km s −1 ,  si ≅ 319 km s −1 , V Ae ≅ 474 km s −1 , V Ai = 350 km s −1 , and plasma betas  i ≅ 1.0 and  e = 0.14.The magnetic fields ratio, , is equal to 1.32 giving  i ≅ 5.1 G.In this case, the speed ordering is  se <  si < V Ai < V Ae and the kink ( = 1) mode propagating in a rest magnetic flux tube is a pseudosurface (body) wave of type B − + (see Table I in [86]).The numerical calculations confirm this mode type and yield a normalized value of the kink speed very close to the predicted one of 1.1858.Now the pattern of stable kink-speed dispersion curves is more complicated (see Figure 5(a))-while the regular lower kink-speed curves correspond to pseudosurface (body) waves, the higher kinkspeed dispersion curves have the form of narrow and long semiclosed up to  A = 2.3 loops and are in fact bulk waves: both attenuation coefficients are purely imaginary numbers.As expected, the dispersion curves of unstable kink mode possess complex attenuation coefficients (with positive real and imaginary parts for  0i and positive real and negative imaginary part for  0e ).The threshold Alfvén Mach number being equal to 2.359305 determines a critical flow velocity of 825.8 km s −1 for KH instability onset which is much higher than the jet #16 speed of 532 km s −1 .This consideration shows that even a relatively high-speed observed soft X-ray jet Advances in Astronomy cannot become unstable-the reason for that in this case is the low density of the jet's plasma.

Kelvin-Helmholtz Instability in Twisted Flux
Tubes.A natural question that immediately raises is how the twist of the internal magnetic field, B i , will change the conditions for instability occurrence and the shape of the dispersion curves.When we model the soft X-ray jet as a moving twisted magnetic flux tubes, there are two types of instabilities that can develop in the jet: namely, kink instability due to the twist of the magnetic field and KH instability owing to the tangential discontinuity of plasma velocity at the tube boundary.According to Dungey and Loughhead [89], the kink instability will occur provided  i () > 2 i , or in our notation, as  > 2. Further on, in order to prevent the developing of kink instability, we will consider weakly twisted flux tubes,  < 1, and, thus, only KH instability can take place in such a jet configuration.When numerically solving dispersion equation (25) we need the normalization of the local Alfvén frequencies  Ai,e .The normalization of these frequencies is performed by multiplying each of them by the tube radius, , and dividing by the Alfvén speed where  twist =  e / i =  √ 1 +  2 .For small values of , we will take  twist ≅ .We begin our calculations for the kink ( = 1) mode in the moving twisted tube (modeling jet #11) with  = 0.896,  twist ≅  = 2.44, and  = 0.025.The results of numerical task are presented in Figure 6, from which one sees that the kink mode dispersion and growth rate curves are very similar to those shown in Figure 3.The critical flow velocity for emerging KH instability now is V cr 0 = 442.6 km s −1 , calculated by using the "reduced" Alfvén speed V Ai / √ 1 +  2 = 110.96km s −1 .KH instability of the ( = −1) mode will occur if the jet speed exceeds 442 km s −1 , of which value is beyond the observationally derived jet speed of 437 km s −1 .We note that while the wave attenuation coefficient of the  = 1 mode propagating in untwisted moving magnetic flux tube of incompressible plasma surrounded by a cool medium is not changed by the flow, in the twisted flux tube (in the same approximations) that attenuation coefficient becomes a complex number with positive imaginary part, like the attenuation coefficient in the cool environment.
According to the instability criterion (30), the occurrence of KH instability of higher MHD modes would require lower threshold Alfvén Mach numbers (and accordingly lower critical jet speeds) than those of the kink ( = ±1) mode.Our computations show, however, that, with  = 0.025, the threshold Alfvén Mach numbers are generally still high to initiate the KH instability.Their magnitudes for the  = 2, 3, and 4 modes are 3.965, 3.943, and 3.926, respectively.The corresponding critical velocities for instability onset are accordingly equal to 440.0, 437.5, and 435.6 km s −1 .As seen, only the  = 4 MHD mode can become unstable against KH instability.A clutch of unstable dispersion curves and normalized growth rates are shown in Figure 7.Here we observe a new phenomenon: the KH instability starts at some critical   -number along with the corresponding threshold Alfvén Mach number.That critical dimensionless wave number for the  = 4 MHD mode (look at Figure 7) is equal to 0.528.If we assume that jet #11 has a width Δℓ = 5 × 10 3 km, then the critical wavelength for a KH instability emergence is  =4 cr ≅ 29.8 Mm.A substantial decrease in the threshold Alfvén Mach number in moving twisted magnetic flux tubes can be obtained by increasing the magnetic field twist parameter, say taking, for instance,  = 0.4.The dispersion curves and normalized wave growth rates of the kink ( = ±1) mode are plotted in Figures 8 and 9.One is immediately seeing the rather complicated forms of both the dispersion curves and dimensionless growth rates of the  = −1 kink mode-such a complication was absent at  = 0.025.We note that computations were performed at  twist = / √ 1 +  2 , and the reference "reduced" Alfvén speed inside the jet now is V Ai / √ 1 +  2 = 103 km s −1 .With this Alfvén speed the critical velocities for the instability occurrence are equal to 416 km s −1 (for the  = 1 mode) and ≅408 km s −1 for the  = −1 mode, respectively.In solving dispersion equation (25) for the  = 2 MHD mode, we start with a little bit bigger Alfvén Mach number than that predicted by instability criterion (30) and have obtained three different kinds of dispersion curves and normalized growth rates (see Figure 10), notably curves with distinctive change of the normalized wave phase velocity in a relatively narrow   -interval (the orange and blue curves in Figure 10(a)), two piece-wise dispersion curves (the purple and green curves in the same plot), and one almost linear dispersion curve (the red one which is in fact the marginal dispersion curve obtained for the threshold Alfvén Mach number equal to 4.062).Not less interesting are the dimensionless growth rate curves corresponding to these three kinds of dispersion curves.For the threshold Alfvén Mach number of 4.062 we obtained one instability window and instability starts at the critical dimensionless wave number of 1.944, or equivalently at  =2 cr ≅ 8.0 Mm.For the slightly superthreshold Alfvén Mach numbers of 4.075 and 4.1 we got two different instability windows (see the green and purple curves in Figure 10(b)) and a further increase in  A leads to merging of those two instability windows-see the blue and orange curves in the same plot.The critical flow velocity at which KH instability arises is equal to ≅418 km s −1 .We have skipped the complicated dispersion and growth rate curves pictured in Figure 10 for the  = 3 and  = 4 MHD modes-we have calculated only a few curves for Alfvén Mach numbers close to the corresponding threshold ones (see Figures 11 and 12).One can observe that for the magnetic field twist parameter  = 0.4 the critical dimensionless wave numbers are shifted far on the righttheir values for both modes are equal to 3.619 and 4.512, respectively, that yield  =3 cr ≅ 4.3 Mm and  =4 cr ≅ 3.5 Mm.The critical jet speeds at which KH instability emerges are equal correspondingly to ≅422 and 424 km s −1 .We note that all the threshold Alfvén Mach numbers of higher MHD modes ( ⩾ 2) traveling on the moving twisted magnetic flux tube are larger than the predicted ones, but while for  = 0.025 they (Alfvén Mach numbers) are decreasing with increasing the mode number, at  = 0.4 we have just the opposite ordering.We would like to notice that finding the marginal dispersion and growth rate curves of the higher modes at  = 0.4 turns out to be a laborious computational task.
To finish our survey on KH instability of MHD modes in moving twisted magnetic flux tubes, we will briefly consider how a weak twist of the internal magnetic field,  = 0.025, will change the critical jet speeds for the occurrence of KH instability in jet #8.We are not going to graphically present the dispersion curves and growth rates of unstable modes because simply there is nothing special in their shapein the incompressible plasma jet approximation and cool environment at small magnetic field twist the curves look, more or less, similar.Recall that at density contrast of 0.963 and magnetic fields ratio  = 4.56, the threshold Alfvén Mach number for the kink ( = 1) mode is equal to 6.9047,

Summary and Conclusion
In this paper, we have explored the conditions under which MHD waves with various mode numbers ( = ±1, 2, 3, and 4) propagating along untwisted and twisted magnetic flux tubes (representing three observed X-ray solar jets) can become unstable against the KH instability.We have used two models of a soft X-ray jet, namely, considering the jet and its surrounding coronal medium as compressible magnetized plasmas and another simplified model representing the jet as an incompressible medium and its environment as cool coronal plasma.A comparison of the dispersion curves and normalized wave growth rates of the kink ( = 1) mode calculated on the base of these two models shows that one obtains similar dispersion curves' and growth rates' patterns (see Figures 2 and 3).This circumstance justifies the usage of the second, simplified, model in studying the propagation characteristics of MHD modes in moving twisted magnetic flux tubes.We must, however, immediately underline that the onset of KH instability in incompressible magnetic flux tube always requires slightly greater threshold Alfvén Mach numbers and correspondingly higher critical jet speeds.Our numerical computations show the following: (i)  rest magnetic flux tube, but unstable ones due to high enough jet speed become partly surface and partly leaky waves.That is why it is more appropriate to term them generalized surface waves.
(ii) The KH instability of the kink mode ( = ±1) starts at Alfvén Mach numbers very close to the predicted ones.Numerically found threshold  A s for the higher MHD modes are, in fact, greater than their initial evaluations.Accessible flow velocities for instability onset in jet #11 have been derived for moving untwisted magnetic flux tube of compressible plasma (429 km s −1 ) and for the  = 4 MHD mode propagating in weakly twisted magnetic flux tube ( = 0.025) of incompressible plasma surrounded by cool coronal medium (435.6 km s −1 ).A nonnegligible decrease of the critical jet speed for all considered modes running in twisted moving flux tube is obtained when the magnetic field twist parameter is  = 0.4-then the critical velocities of the  = ±1, 2, 3, and 4 modes are equal to 416, ≅408, 418, ≅422, and 424 km s −1 , respectively.The instability onset of all considered MHD modes in jet #8 occurs at flow velocities, in general, lower than the jet speed of 385 km s −1 -their magnitudes for the kink ( = 1) mode running on untwisted or slightly twisted magnetic flux tube (with magnetic field twist parameter  = 0.025) are equal to 348.3 and 360.6 km s −1 , respectively.The corresponding critical flow velocities for the higher ( = 2, 3, and 4) MHD modes turn out to lie between aforementioned speeds, namely, being equal to ≅359, 35,7 and 355.5 km s −1 .Jet #16 possessing the highest flow speed of 532 km s −1 in Shimojo and Shibata's Events List [10] is, however, stable against the KH instability-the critical flow velocity at which the kink ( = 1) mode would become unstable propagating in untwisted magnetic flux tube is equal to ≅826 km s −1 -much higher than the jet speed.The critical velocities in twisted tubes are even larger, equal to 852.4,848.3, 844.8, and 842.6 km s −1 for the  = 1-4 MHD modes, respectively.
(iii) The instability onset for a given mode  is very sensitive to the main input parameters: density contrast, , magnetic field ratio, , the magnetic field twist , and especially the background magnetic field  e .If we assume that, for jet #11  e = 10 G, Alfvén speeds in the two media have new values, notably V Ai = 320.38km s −1 and V Ae = 427.5 km s −1 .With these new Alfvén speeds the input parameters for solving wave dispersion relations ( 9) and ( 25) have changed and more specifically βi = 0.7376, βe = 0.1506, and  = 1.2636.Now the critical jet velocity for the instability onset dramatically increases-it becomes 751 km s −1 for the kink mode in untwisted tube and 774 km s −1 in weakly twisted ( = 0.025) flux tube.
(iv) If we accept the Paraschiv's et al. [90] suggestion that plasma density in the region crossed by coronal jets has to be the sum  cor +  jet or, on average, two times higher than the plasma in the surrounding corona  cor ; then the total pressure balance equation, again for jet #11, can be satisfied at  e = 10 G and we have the following set of input parameters:  = 0.437, V Ai = 105.7 km s −1 , V Ae = 427.5 km s −1 , βi = 6.7782, βe = 0.1506, and  = 3.781.With these parameters the critical flow velocities for KH instability occurrence reduce to 560 km s −1 for the kink mode in untwisted flux tube and to 575 km s −1 for the same mode in weakly twisted ( = 0.025) flux tube.Note that within this approach one obtains critical jet speeds that are approximately with 200 km s −1 less than those calculated for the same background magnetic field of 10 G.Here one raises the big question: "which approach, the standard one or that of Paraschiv's et al. [90], will be appropriate for studying/modeling the KH instability in solar X-ray (and other) jets?"The perfect answer to this question will be obtained when we will try to model a really observationally detected KH instability in soft X-ray jets like that of Foullon et al. [57].
It is intriguing to see what critical jet speeds we will obtain for a numerically modeled solar X-jet-a case in point is the study of Miyagoshi and Yokoyama [91].These authors have presented MHD numerical simulations of solar X-ray jets based on magnetic reconnection model that includes chromospheric evaporation.Peculiar to their study is that total pressure balance equation excludes the magnetic pressure inside the jet saying that it is very weak; that is (in their notation), where  ≡  e is the background (coronal) magnetic field.Supposing that  cor = 10 9 cm −3 ,  cor = 10 6 K,  = 10 G, and  jet = 4.5 × 10 9 cm −3 , they obtain from their total balance equation the temperature of the jet,  jet = 6.7 × 10 6 K-indeed a reasonable value.However, to satisfy the standard total pressure balance (1) (assuming a very small, but finite value for the magnetic field in the jet,  i ), with the aforementioned values for both the temperature and jet density, we were forced to double the coronal density, that is, to take  cor ≡  e = 2.0 × 10 9 cm −3 .In such a case the sound and Alfvén speeds in the jet and its environment are  si ≅ 304 km s −1 ,  se = 117 km s −1 , V Ai = 47.67 km s −1 , and V Ae = 487.5 km s −1 .The ordering of sound and Alfvén speeds is the same as for the #11 X-ray jet in the Shimojo and Shibata paper [10] that implies the propagation of pure surface stable modes.The magnetic field inside the jet is  i ≅ 1.47 G-hence the magnetic fields ratio is  = 6.817.Thus the input parameters for solving (9) and ( 25) are  = 0.444, βi = 40.5872,and βe = 0.0695 and, as we already found,  = 6.817.Note that this relatively high value of  would suspect a rather big value of the threshold Alfvén Mach number-the instability criterion yields  A > 12.42.Numerically found threshold  A for initiating KH instability of the kink ( = 1) mode in moving untwisted magnetic flux tube of compressible plasma is 12.435 that yields a critical velocity of 607 km s −1 .For the case of weakly twisted flux tube ( = 0.025) in the approximation incompressible jet and cool plasma environment we obtained for the same kink mode a critical velocity of 620 km s −1 .The critical speeds for the higher ( = 2-4) modes are accordingly of 615, 609, and 606 km s −1 .
The three examples of coronal X-jets embedded in a background magnetic field of 10 G show that KH instability of basically the kink ( = 1) mode can occur if jets' speeds lie in the range of 500-700 km s −1 .Such velocities are generally accessible for high-speed jets; lower critical speeds between 400 and 500 km s −1 can be observed/detected at moderate external magnetic fields in the range of 6-7 G.This requirement is in agreement with Pucci et al. [34] evaluation that the magnetic field inside a standard X-ray jet would be equal to 2.8 G, while for a blowout jet that value should be around 4.5 G.To be honest, we can say that the most of soft X-ray solar jets with speeds below 400-450 km s −1 are stable against the KH instability-this instability can develop only in relatively dense high-speed jets.Note, however, that each soft X-ray jet is a unique event-it requires a separate careful exploration of the instability conditions.Arising KH instability in turn might trigger wave turbulence considered as an effective mechanism for coronal heating.
Our simplified model of studying the possibilities for emerging of KH instability in fast soft X-ray jets has to be improved by assuming a radial density gradient or some magnetic field and flow velocity shears.A challenge also is the modeling of the same instability in rotating and oscillatory swaying coronal hole X-ray jets.

Figure 1 :
Figure 1: Equilibrium magnetic fields of a soft X-ray solar jet in an untwisted flux tube (a) and in a weakly twisted flux tube (b).

Figure 2 :Figure 3 :
Figure 2: (a) Dispersion curves of stable and unstable kink ( = 1) MHD mode propagating in a moving untwisted magnetic flux tube of compressible plasma (modeling jet #11) at  = 0.896 and  = 2.44.Unstable dispersion curves located above the middle of the plot have been calculated for four values of the Alfvén Mach number  A = 3.865, 3.89, 3.915, and 3.94.(b) The normalized growth rates of the unstable mode for the same values of  A .Red curves in both plots correspond to the onset of KH instability.

Figure 4 :
Figure 4: (a) Dispersion curves of stable and unstable kink ( = 1) MHD mode propagating in a moving untwisted magnetic flux tube of compressible plasma (modeling jet #8) at  = 0.963 and  = 4.56.Unstable mode's dispersion curves are pictured in the middle of the plot for  A = 6.6704, 6.68, 6.69, and 6.7.(b) The corresponding normalized growth rates of unstable modes for the same values of  A .

Figure 5 :
Figure 5: (a) Dispersion curves of stable and unstable kink ( = 1) MHD mode propagating in a moving untwisted magnetic flux tube of compressible plasma (modeling jet #16) at  = 0.95 and  = 1.32.Dispersion curves of unstable modes located in the middle of the plot have been calculated for  A = 2.359305, 2.37, 2.38, and 2.39.(b) Plots of corresponding normalized growth rates.

Figure 6 :Figure 7 :
Figure 6: (a) Dispersion curves of unstable kink ( = 1) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at  = 0.896 and  = 0.025 and for  A = 3.989, 4.0 (green curve), 4.011 (purple curve), and 4.022.(b) The normalized growth rates for the same values of  A .

Figure 8 :Figure 9 :
Figure 8: (a) Dispersion curves of unstable ( = 1) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at  = 0.896 and  = 0.4 and for  A = 4.04, 4.07 (green curve), 4.1 (purple curve), and 4.13.(b) The normalized growth rates for the same values of  A .

Figure 10 :
Figure 10: (a) Dispersion curves of unstable ( = 2) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at the same parameters as in Figure 8 and for  A = 4.2, 4.15, 4.1, 4.075, and 4.062.(b) The normalized growth rates for the same values of  A .

Figure 11 :
Figure 11: (a) Dispersion curves of unstable ( = 3) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at the same parameters as in Figure 8 and for  A = 4.093, 4.1, 4.109, and 4.125.(b) The normalized growth rates for the same values of  A .

Figure 12 :
Figure 12: (a) Dispersion curves of unstable ( = 4) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at the same parameters as in Figure 8 and for  A = 4.12, 4.13, 4.14, and 4.15.(b) The normalized growth rates for the same values of  A .