How Rotating Solar Atmospheric Jets Become Kelvin--Helmholtz Unstable

Recent observations support the propagation of a number of magnetohydrodynamic (MHD) modes which, under some conditions, can become unstable and the developing instability is the Kelvin--Helmholtz instability (KHI). In its nonlinear stage the KHI can trigger the occurrence of wave turbulence which is considered as a candidate mechanism for coronal heating. We review the modeling of tornado-like phenomena in the solar chromosphere and corona as moving weakly twisted and spinning cylindrical flux tubes, showing that the KHI rises at the excitation of high-mode MHD waves. The instability occurs within a wavenumber range whose width depends on the MHD mode number \emph{m}, the plasma density contrast between the rotating jet and its environment, and also on the twists of the internal magnetic field and the jet velocity. We have studied KHI in two twisted spinning solar polar coronal hole jets, in a twisted rotating jet emerging from a filament eruption, and in a rotating macrospicule. The theoretically calculated KHI development times of a few minutes for wavelengths comparable to the half-widths of the jets are in good agreement with the observationally determined growth times only for high order (10 $\mathrm{\leqslant}$ \emph{m} $\mathrm{\leqslant}$ 65) MHD modes. Therefore, we expect that the observed KHI in these cases is due to unstable high-order MHD modes.


INTRODUCTION
Solar jets are ubiquitous in the solar atmosphere and recent observations have revealed that they are related to small scale filament eruptions. They are continuously observed by the Extremeultraviolet Imaging Spectrometer (EIS) (Culhane et al., 2007) on board Hinode (Kosugi et al., 2007) satellite, Atmospheric Imaging Assembly (AIA) (Lemen et al., 2012), on board the Solar Dynamics Observatory (SDO) (Pesnell et al., 2012), as well as from the Interface Region Imaging Spectrograph (IRIS) (De Pontieu et al., 2014) alongside the Earth-based solar telescopes. The physical parameters of various kinds of solar jets have been reported in a series of articles, see for instance, Schmieder et al. (2013), Sterling et al. (2015), Panesar, Sterling, and Moore (2016a), Chandra et al. (2017), Joshi et al. (2017), and references cited in. It was established that more of the solar jets possess rotational motion. Such tornado-like jets, termed macrospicules, were firstly detected in the transition region by Pike and Mason (1998) using observations by the Solar and Heliospheric Observatory (SOHO) (Domingo, Fleck, and Poland, 1995). Rotational motion in macrospicules was also explored by Kamio et al. (2010), Curdt and Tian (2011), Bennett and Erdélyi (2015), Kiss, Gyenge, and Erdélyi (2017), and Kiss, Gyenge, and Erdélyi (2018). Type II spicules, according to De Pontieu et al. (2012) and Martínez-Sykora et al. (2013), along with the coronal hole EUV jets (Nisticò et al., 2009;Liu et al., 2009;Nisticò et al., 2010;Shen et al., 2011;Chen, Zhang, and Ma, 2012;Hong et al., 2013;Young and Muglach, 2014a,b;Moore, Sterling, and Falconer, 2015), and X-ray jets (Moore et al., 2013), can rotate, too. Rotating EUV jet emerging from a swirling flare (Zhang and Ji, 2014) or formed during a confined filament eruption (Filippov et al., 2015) confirm once again the circumstance that the rotational motion is a common property of many kinds of jets in the solar atmosphere.
The first scenario for the numerical modeling of hot X-ray jets was reported by Heyvaerts, Priest, and Rust (1997) and the basic idea was that a bipolar magnetic structure emerges into a unipolar pre-existing magnetic field and reconnects to form hot and fast jets that are emitted from the interface between the fields into contact. Later on, by examining many X-ray jets in Hinode/X-Ray Telescope coronal X-ray movies of the polar coronal holes, Moore et al. (2010) found that there is a dichotomy of polar X-ray jets, namely "standard" and "blowout" jets exist. Fang, Fan, and McIntosh (2014) studied the formation of rotating coronal jets through numerical simulation of the emergence of a twisted magnetic flux rope into a pre-existing open magnetic field. Another scenario for the nature of solar jets was suggested by Sterling et al. (2015), according to which the X-ray jets are due to flux cancellation and/or "minieruptions" rather than emergence. An alternative model for solar polar jets due to an explosive release of energy via reconnection was reported by Pariat, Antiochos, and DeVore (2009). Using three-dimensional MHD simulations, the authors demonstrated that this mechanism does produce massive, high-speed jets. In subsequent two articles (Pariat et al., 2015(Pariat et al., , 2016, Pariat and co-authors, presented several parametric studies of a three-dimensional numerical MHD model for straight and helical solar jets. On the other side, Panesar et al. (2016b) have shown that the magnetic flux cancellation can trigger the solar quiet-region coronal jets and they claim that the coronal jets are driven by the eruption of a small-scale filament, called a "minifilament." The small-scale chromospheric jets, like microspicules, were first numerically modeled by Murawski, Srivastava, and Zaqarashvili (2011). Using the FLASH code, they solved the two-dimensional ideal MHD equations to model a macrospicule, whose physical parameters match those of a solar spicule observed. Another mechanism for the origin of macrospicules was proposed by Kayshap et al. (2013), who numerically modeled the triggering of a macrospicule and a jet.
It is natural to expect, that solar jets, being magnetically structured entities, should support the propagation of different type of MHD waves: fast and slow magnetoacoustic waves and torsional Alfvén waves. All these waves are usually considered as normal MHD modes traveling along the jet. Owing to the presence of a velocity shear near the jet-surrounding plasma interface, every jet can become unstable and the most universal instability which emerges is the Kelvin-Helmholtz (KH) one. The simplest configuration at which one can observe the KHI is the two semi-infinite incompressible magnetized plasmas flowing with different velocities provided that the thin velocity shear at the interface exceeds some critical value (Chandrasekhar, 1961). Recently, Cheremnykh et al. (2018a) theoretically established that shear plasma flows at the boundary of plasma media can generate eight MHD modes, of which only one can be unstable due to the development of the KHI. Ismayilli et al. (2018) investigated a shear instability of the KH type in a plasma with temperature anisotropy under the MHD approximation. The KHI of the magnetoacousic waves propagating in a steady asymmetric slab, and more specifically the effect of varying density ratios was explored by Barbulescu and Erdélyi (2018). A very good review on the KHI in the solar atmosphere, solar wind, and geomagnetosphere in the framework of ideal MHD the reader can find in Mishin and Tomozov, 2016. In cylindrical geometry, being typical for the solar jets, the KHI exhibits itself as a vortex sheet running on the jet-environment boundary, which like in the flat geometry, is growing in time if the axial velocity of the jet in a frame of reference attached to the surrounding plasma exceeds a threshold value (Ryu, Jones, and Frank, 2000). In its nonlinear stage, the KHI trigger the wave turbulence which is considered as one of the main heating mechanisms of the solar corona (Cranmer et al, 2015). The development of the KHI in various cylindrical jet-environment configurations has been studied in photospheric jets (Zhelyazkov and Zaqarashvili, 2012), in solar spicules (Zhelyazkov, 2012;Ajabshirizadeh et al., 2015;Ebadi, 2016), in high-temperature and cool solar surges (Zhelyazkov et al., 2015a,b), in magnetic tubes of partially ionized compressible plasma (Soler et al., 2015), in EUV chromospheric jets (Zhelyazkov, Chandra, and Srivastava, 2016;Bogdanova et al., 2018), in soft X-ray jets (Vasheghani Farahani et al., 2009;Zhelyazkov, Chandra, and Srivastava, 2017), and in the twisted solar wind flows (Zaqarashvili et al., 2014). A review on KHI in the solar atmosphere, including some earlier studies, the reader can find in Zhelyazkov, 2015. The first modeling of the KHI in a rotating cylindrical magnetized plasma jet was done by Bondenson, Iacono, and Bhattacharjee (1987). Later on, Bodo et al. (1989Bodo et al. ( , 1996 carried out a study of the stability of flowing cylindrical jet immersed in constant magnetic field B 0 . The authors used the standard procedure for exploring the MHD wave propagation in cylindrical flows considering that all the perturbations of the plasma pressure p, fluid velocity v, and magnetic field B, are ∝ exp[i(−ωt+kz+mθ)].
Here, ω is the angular wave frequency, k the propagating wavenumber, and m the azimuthal mode number. Using the basic equations of ideal magnetohydrodynamics, Bodo et al. (1989Bodo et al. ( , 1996 derived a Bessel equation for the pressure perturbation and an expression for the radial component of the fluid velocity perturbation. The found solutions in both media (the jet and its environment) are merged at the perturbed tube boundary through the conditions for continuity of the total (thermal plus magnetic) pressure and the Lagrangian displacement. The latter is defined as the ratio of radial velocity perturbation component and the angular frequency in the corresponding medium. The obtained dispersion relation is used for examining the stability conditions of both axisymmetric, m = 0 (Bodo et al., 1989), and nonaxisymmetric, |m| 1 modes (Bodo et al., 1996). In a recent article, Bodo et al. (2016) performed a linear stability analysis of magnetized rotating cylindrical jet flows in the approximation of zero thermal pressure. They focused their analysis on the effect of rotation on the current driven mode and on the unstable modes introduced by rotation. In particular, they found that rotation has a stabilizing effect on the current driven mode only for rotation velocities of the order on the Alfvén speed. The more general case, when both the magnetic field and jet flow velocity are twisted, was studied by Zaqarashvili, Zhelyazkov, and Ofman (2015) and Cheremnykh et al. (2018b), whose dispersion equations for modes with m 2, represented in different ways, yield practically identical results.
The main goal of this review article is to suggest a way of using the wave dispersion relation derived in Zaqarashvili, Zhelyazkov, and Ofman, 2015 to study the possibility for the rising and development of KHI in rotating twisted solar jets. Among the enormous large number of observational studies of rotating jets with different origin or nature, we chose those which provide the magnitudes of axial and rotational speeds, jet width and height alongside the typical plasma parameters like electron number densities and electron temperature of the spinning structure and its environment. Thus, the targets of our exploration are: (i) the spinning coronal hole jet of 2010 August 21 (Chen, Zhang, and Ma, 2012); (ii) the rotating coronal hole jet of 2011 February 8 (Young and Muglach, 2014a), (iii) the twisted rotating jet emerging from a filament eruption on 2013 April 10-11 (Filippov et al., 2015), and (iv) the rotating macrospicule observed by Pike and Mason (1998) on 1997 March 8.
The paper is organized as follows: in the next section, we discuss the geometry of the problem, equilibrium magnetic field configuration and basic physical parameters of the explored jets. Section 3 is devoted to a short, concise, derivation of the wave dispersion relation. Section 4 deals with numerical results for each of the four jets and contain the available observational data. In the last Section 5, we summarize the main findings in our research and outlook the further improvement of the used modeling.

THE GEOMETRY, MAGNETIC FIELD, AND PHYSICAL PARAMETERS IN A JET MODEL
We model whichever jet as an axisymmetric cylindrical magnetic flux tube with radius a and electron number density n i (or equivalently, homogeneous plasma density ρ i ) moving with velocity U . We consider that the jet environment is a rest plasma with homogeneous density ρ e immersed in a homogeneous background magnetic field B e . This field, in cylindrical coordinates (r, φ, z), possesses only an axial component, i.e., B e = (0, 0, B e ). (Note that the label 'i' is abbreviation for interior, and the label 'e' denotes exterior.) The magnetic field inside the tube, B i , and the jet velocity, U , we assume, are uniformly twisted and are given by the vectors respectively. We note, that B iz and U z , are constant. Concerning the azimuthal magnetic and flow velocity components, we suppose that they are linear functions of the radial position r and evaluated at r = a they correspondingly are equal to B iφ (a) ≡ B φ = Aa and U φ = Ωa, where A and Ω are constants. Here, Ω is the jet angular speed, deduced from the observations. Hence, in equilibrium, the rigidly rotating plasma column, that models the jet, must satisfy the following force-balance equation (see, e.g., Chandrasekhar, 1961;Goossens, Hollweg, and Sakurai, 1992) where µ is the plasma permeability and p t = p i + B 2 i /2µ with B 2 i = B 2 iφ (r) + B 2 iz is the total (thermal plus magnetic) pressure. According to equation (2), the radial gradient of the total pressure should balance the centrifugal force and the force owing to the magnetic tension. After integrating equation (2) from 0 to a, taking into account the linear dependence of U φ and B iφ on r, we obtain that where p t (0) = p 1 (0) + B 2 iz /2µ. (Integrating equation (2) from 0 to any r one can find the radial profile of p t inside the tube. Such an expression of p t (r), obtained, however, from an integration of the momentum equation for the equilibrium variables, have been obtained in Zhelyazkov et al., 2018a-see Eq. (2) there.) It is clear from a physical point of view that the internal total pressure (evaluated at r = a) must be balanced by the total pressure of the surrounding plasma which implies that This equation can be presented in the form where p 1 (0) is the thermal pressure at the magnetic tube axis, and p e denotes the thermal pressure in the environment. In the pressure balance equation (3), the number ε 1 ≡ B φ /B iz = Aa/B iz represents the magnetic field twist parameter. Similarly, we define ε 2 ≡ U φ /U z as a characteristics of the jet velocity twist. We would like to underline that the choice of plasma and environment parameters must be such that the total pressure balance equation (3) is satisfied. In our case, the value of ε 2 is fixed by observationally measured rotational and axial velocities while the magnetic field twist, ε 1 , has to be specified when using equation (3). We have to note that equation (3) is a corrected version of the pressure balance equation used in Zhelyazkov et al., 2018a;Zhelyazkov and Chandra, 2018b.
From measurements of n and T for similar coronal hole EUV jets (Nisticò et al., 2009(Nisticò et al., , 2010, we take n inside the jet to be n i = 1.0 × 10 9 cm −3 , and assume that the electron temperature is T i = 1.6 MK. The same quantities in the environment are respectively n e = 0.9 × 10 9 cm −3 and T e = 1.0 MK. Note that the electron number density of the blowout jet observed by Young and Muglach (2014a) is in one order lower. The same applies for its environment. We consider that the background magnetic field for both hole coronal jets is B e = 3 G. The values of n and T of the rotating jet emerging from a filament eruption, observed by Filippov et al. (2015), were evaluated by us and they are n i = 4.65 × 10 9 cm −3 and T i = 2.0 MK, respectively. From the same data set, we have obtained n e = 4.02 × 10 9 cm −3 and T e = 2.14 MK. The background magnetic field, B e , with which the pressure balance equation (3) is satisfied, is equal to 6 G. For the rotating macrospicule we assume that n i = 1.0 × 10 10 cm −3 and n e = 1.0 × 10 9 cm −3 to have at least one order denser jet with respect to the surrounding plasma. Our choice for macrospicule temperature is T i = 5.0 × 10 5 K, while that of its environment is supposed to be T e = 1.0 × 10 6 K. The external magnetic field, B e , was taken as 5 G. All aforementioned physical parameters of the jets are summarized in Table 1. The plasma beta was calculated using (6/5)c 2 s /v 2 A , where c s = (γk B T /m ion ) 1/2 is the sound speed (in which γ = 5/3, k B is the Boltzmann's constant, T the electron temperature, and m ion the ion or proton mass), and v A = B/(µn ion m ion ) 1/2 is the Alfvén speed, in which expression B is the full magnetic field =(B 2 φ + B 2 z ) 1/2 , and n ion is the ion or proton number density.

WAVE DISPERSION RELATION
A dispersion relation for the propagation of high-mode (m 2) MHD waves in a magnetized axially moving and rotating twisted jet was derived by Zaqarashvili, Zhelyazkov, and Ofman (2015) and Cheremnykh et al. (2018b). That equation was obtained, however, under the assumption that both media (the jet and its environment) are incompressible plasmas. As seen from the last column in Table 1, plasma beta is greater than 1 in the first, third, and forth jets which implies that the plasma of each of the aforementioned jets can be considered as a nearly incompressible fluid (Zank and Matthaeus, 1993). It is seen from the same table that the plasma beta of the second jet is less than one as is in each of the jet environments and that is why it is reasonable to treat them as cool media. Thus, the wave dispersion relation, derived, for instance, in Zaqarashvili, Zhelyazkov, and Ofman, 2015, has to be modified. In fact, we need two modified versions: one for the incompressible jet-cool environment configuration, and other for the cool jet-cool environment configuration. We are not going to present in details the derivation of the modified dispersion equations on the basis of the governing MHD equations, but will only sketch the essential steps in that procedure. The main philosophy in deriving the wave dispersion equation is to find solutions for the total pressure perturbation, p tot , and for the radial component, ξ r of the Lagrangian displacement, ξ, and merge them at the tube perturbed boundary through the boundary conditions for their (p tot and ξ r ) continuity (Chandrasekhar, 1961). In the case of the first configuration, we start with the linearized ideal MHD equations, governing the incompressible dynamics of the perturbations in the spinning jet where are the perturbations of fluid velocity and magnetic field, respectively, and p tot is the perturbation of the total pressure, p t = p i + B 2 i /2µ. The Lagrangian displacement, ξ, can be found from the fluid velocity perturbation, v, using the relation (Chandrasekhar, 1961) Further on, assuming that all perturbations are ∝ exp [i (−ωt + mφ + k z z)] and considering that the rotation and the magnetic field twists in the jet are uniform, that is, where Ω and A are constants, from the above set of equations, (4)-(8), we obtain the following dispersion equation of the MHD wave with mode number m (for details see Zaqarashvili, Zhelyazkov, and Ofman, 2015): where and P m (κ e a) = κ e aK ′ m (κ e a) K m (κ e a) .
In above expressions, the prime means differentiation of the Bessel functions with respect to their arguments, are the squared wave amplitude attenuation coefficients in the jet and its environment, in which are the local Alfvén frequencies in both media, and is the Doppler-shifted angular wave frequency in the jet. We note that in the case of incompressible coronal plasma (Zaqarashvili, Zhelyazkov, and Ofman, 2015), κ e = k z , because at an incompressible environment the argument of the modified Bessel function of second kind, K m , and its derivative, The basic MHD equations for an ideal cool plasma are, generally, the same as the set of equations (4)-(8) with equation (6) replaced by the continuity equation Recall that for cold plasmas the total pressure reduces to the magnetic pressure only, that is p The above equation, which defines the density perturbation, is not used in the derivation of the wave dispersion relation because we are studying the propagation and stability of Alfvén-wave-like perturbations of the fluid velocity and magnetic field. Following the standard scenario for deriving the MHD wave dispersion relation (Zhelyazkov and Chandra, 2018b), we finally arrive at: where .
Here, the wave attenuation coefficient in the internal medium has the form while that in the environment, with Ω = 0 and A = 0, is given by Note that (i) both dispersion relations, (10) and (11), have similar forms-the difference is in the expressions for the wave attenuation coefficient inside the jet, namely κ c i = κ i 1 − σ 2 /ω 2 Ai 1/2 ; and (ii) the wave attenuation coefficients in the environments are not surprisingly the same, that is, κ c e = κ e ≡ 1 − ω 2 /ω 2 Ae 1/2 .

NUMERICAL SOLUTIONS, WAVE DISPERSION, AND GROWTH RATE DIAGRAMS
In studying at which conditions the high (m 2) MHD modes in a jet-coronal plasma system become unstable, that is, all the perturbations to grow exponentially in time, we have to consider the wave angular frequency, ω, as a complex quantity: ω ≡ Re(ω) + i Im(ω) in contrast to the wave mode number, m, and propagating wavenumber, k z , which are real quantities. The Re(ω) is responsible for the wave dispersion while the Im(ω) yields the wave growth rate. In the numerical task for finding the complex solutions to the wave dispersion relation (10) or (11), it is convenient to normalize all velocities with respect to the Alfvén speed inside the jet, defined as v Ai = B iz / √ µρ i , and the lengths with respect to a. Thus, we have to search the real and imaginary parts of the nondimensional wave phase velocity, v ph = ω/k z , that is, Re(v ph /v Ai ) and Im(v ph /v Ai ) as functions of the normalized wavenumber k z a. The normalization of the other quantities like the local Alfvén and Doppler-shifted frequencies alongside the Alfvén speed in the environment, v Ae = B e / √ µρ e , requires the usage of both twist parameters, ε 1 and ε 2 , and also of the magnetic fields ratio, b = B e /B iz . The nondimensional form of the jet axial velocity, U z , is given by the Alfvén Mach number M A = U z /v Ai . Another important nondimensional parameter is the density contrast between the jet and its surrounding medium, η = ρ e /ρ i . Hence, the input parameters in the numerical task of finding the solutions to the transcendental equation (10) or (11) Zaqarashvili, Zhelyazkov, and Ofman (2015) have established that KHI in an untwisted (A = 0) rotating flux tube with negligible longitudinal velocity can occur if This inequality says that each MHD wave with mode number m 2, propagating in a rotating jet can become unstable. This instability condition can be used also in the cases of slightly twisted spinning jets, provided that the magnetic field twist parameter, ε 1 , is a number lying in the range of 0.001-0.005, simply because the numerical solutions, for example, to equation (10) show that practically there is no difference between the instability ranges at ε 1 = 0, and at 0.001 or 0.005. An important step in our study is the supposition that the deduced from observations jet axial velocity, U z , is the threshold speed for the KHI occurrence. Then, for fixed values of m, η, U φ = Ωa, v Ai , and b, the inequality (12) can be rearranged to define the upper limit of the instability range on the k z a-axis According to the above inequality, the KHI can occur for nondimensional wavenumbers k z a less than (k z a) rhs . On the other hand, one can talk for instability if the unstable wavelength, λ KH = 2π/k z , is shorter than the height of the jet, H, which means that the lower limit of the instability region is given by: where ∆ℓ is the jet width. Hence, the instability range in the k z a-space is (k z a) lhs < k z a < (k z a) rhs . Note that the lower limit, (k z a) lhs , is fixed by the width and height of the jet, while the upper limit, (k z a) rhs , depends on several jet-environment parameters. At fixed U φ , v Ai , η, and b, the (k z a) rhs is determined by the MHD wave mode number, |m|. As seen from inequality (13), with increasing the m, that limit shifts to the right, that is, the instability range becomes wider. The numerical solutions to the wave dispersion relation (10) confirm this and for given m one can obtain a series of unstable wavelengths, λ KH = π ∆ℓ/k z a, as the shortest one takes place at k z a ≈ (k z a) rhs . For relatively small mode numbers, when even the shortest unstable wavelengths turn out to be a few tens megameters, that could hardly be associated with the observed KH ones. As observations show, the KHI vortex-like structures running at the boundary of the jet, have the size of the width or radius of the flux tube (see, for instance, Fig. 1 in Zhelyazkov et al., 2018a). Therefore, we have to look for such an m, whose instability range would accommodate the expected unstable wavelength presented by its nondimensional wavenumber, k z a = π∆ℓ/λ KH . An estimation of the required mode number for an ε 1 = 0.005-rotating flux tube can be obtained by presenting the instability criterion (12) in the form We will use this inequality for obtaining the optimal m for each of the studied jets by specifying the value of that k z a (along with the other aforementioned input parameters) which corresponds to the expected unstable wavelength λ KH .

Kelvin-Helmholtz Instability in a Standard Polar Coronal Hole Jet
Chen, Zhang, and Ma (2012) observationally studied the jet event of 2010 August 21, which occurred in the coronal hole region, close to the north pole of the Sun. Figure 1 presents the jet's evolution in AIA 304 Å. The jet started around 06:07 UT, reached its maximum height around 06:40 UT. During the evolution of the jet between 06:32 and 06:38 UT, small scale moving blobs appeared on the right boundary. We interpret these blobs, shown by arrows in Figure 1, as evidence of KHI. By tracking six identified moving features in the jet, Chen, Zhang, and Ma (2012) found that the plasma moved at an approximately constant speed along the jet's axis. Inferred from linear and trigonometric fittings to the axial and transverse heights of the six tracks, the authors have found that the mean values of the axial velocity, U z , transfer/rotational velocity, U φ , angular speed, Ω, rotation period, T , and rotation radius, a, are 114 km s −1 , 136 km s −1 , 0.81 • s −1 (or 14.1 × 10 −3 rad s −1 ), 452 s, and 9.8 × 10 3 km, respectively. The height of the jet is evaluated as H = 179 Mm.
It seems reasonable the shortest unstable wavelength, λ KH , to be equal to 10 Mm (approximately half of the jet width, ∆ℓ = 19.6 Mm), which implies that its position on the k z a one-dimensional space is k z a = 6.158. The input parameters, necessary to find out that MHD wave mode number, whose instability range will contain the nondimensional wavenumber of 6.158, using inequality (15), are accordingly (see Table 1) η = 0.9, b = 1.834, v Ai = 112.75 km s −1 , and U φ = 136 km s −1 . (We note, that the values of b and v Ai were obtained with the help of equation (3) assuming that the magnetic field twist is ε 1 = 0.005.) With these entry data, from inequality (15) one obtains that |m| > 15 should provide the required instability region or window. The numerical solutions to equation (10) show that this value of m is overestimatedan m = 11 turns out to be perfect for the case. The discrepancy between the predicted and computed value of |m| is not surprising because inequality (15) yields only an indicative value. The input parameters for finding the solutions to the dispersion equation (10) are as follows: m = 11, η = 0.9, ε 1 = 0.005, ε 2 = 1.2, b = 1.834, and M A = 1.01 (=114/112.75). The results of computations are graphically presented in Figure 2. From that figure, one can obtain the normalized wave phase velocity, Re(v ph /v Ai ), and the normalized growth rate, Im(v ph /v Ai ), of the unstable λ KH = 10 Mm wave, both read at the purple cross points. From the same plot, one can find the instability characteristics at another wavelength, precisely λ KH = 12 Mm, whose position on the k z a-axis is fixed at k z a = 5.131. The values of nondimensional wave phase velocity and growth rate can be read from the green cross points. The KHI wave growth rate, γ KH , growth time, τ KH = 2π/γ KH , and wave velocity, v ph , in absolute units, estimated from the plots in Figure 2, for the two wavelengths, are γ KH ∼ = 23.09 × 10 −3 s −1 , τ KH ∼ = 4.5 min, v ph ∼ = 178 km s −1 , for λ KH = 10 Mm, and γ KH ∼ = 50.65 × 10 −3 s −1 , τ KH ∼ = 2.1 min, v ph ∼ = 202 km s −1 , for λ KH = 12 Mm.
Let us recall that the value of the Alfvén speed used in the normalization is v Ai = 112.75 km s −1 . We see that the two wave phase velocities are slightly super-Alfvénic and when moving along the k z aaxis to the left, the normalized wave velocity becomes higher. If we fix a k z a-position near the lower limit of the unstable region, (k z a) lhs = 0.344, say, at k z a = 0.513, which means λ KH = 120 Mm, the KHI characteristics obtained from the numerical solutions to equation (10) are τ KH = 1.4 min and v ph = 1 473 km s −1 , respectively. As we have discussed in Zhelyazkov et al., 2018a, "the KHI growth time could be estimated from the temporal evolution of the blobs in their initial stage and it was found to be about 2-4 minutes", so the instability developing times of 2.1 and 4.5 min obtained from our plots are in good agreement with the observations. A specific property of the instability k z a-ranges is that for a fixed mode number, m, their widths depend upon ε 1 and with increasing the value of ε 1 , the instability window becomes narrower and at some critical ε 1 its width equals zero. In our case that happens with ε cr 1 = 0.653577 at (k z a) lhs = 0.344. In Figure 3, curves of dimensionless v ph and γ KH have been plotted for several ε 1 values. Note that each larger value of ε 1 implies an increase in B iφ . But that increase in B iφ requires an increase in B iz too, in order the total pressure balance equation (3) to be satisfied under the condition that the hydrodynamic pressure term and the environment total pressure are fixed. The increase in B iz (and in the full magnetic field B i ) implies a decrease both in the magnetic field ratio, b, and in the Alfvén Mach number, M A . Thus, gradually increasing the magnetic field twist ε 1 from 0.005 to 0.653577, we get a series of dispersion and growth rate curves with progressively diminishing parameters b and M A . The red growth rate curve in the right panel of Figure 3 has been obtained for ε cr 1 = 0.653577 with M A = 0.7652 and it visually fixes the lower limit of all other instability windows. The azimuthal magnetic field B cr iφ that stops the KHI, computed at B i = 2.58 G, is equal to 1.4 G.

Kelvin-Helmholtz Instability in a Blowout Polar Coronal Hole Jet
Young and Muglach (2014a) observed a small blowout jet at the boundary of the south polar coronal hole on 2011 February 8 at around 21:00 UT. The evolution of jet observed by the AIA is displayed in Figure 4. The jet activity was between 20:50 to 21:15 UT. This coronal hole is centered around x = −400 arcsec, y = −400 arcsec. The jet has very broad and faint structure and is ejected in the south direction. We could see the evolution of jet in the AIA 193 Å clearly. However, in AIA 304 Å the whole jet is not visible. Moreover, we observe the eastern boundary of the jet in AIA 304 Å. During its evolution in 304 Å we found the blob structures at the jet boundary. These blobs could be due to the KHI as reported in previous observations (see, for example Zhelyazkov et al., 2018a). At the jet initiation/base site, we observed the coronal hole bright points. These bright points are the results of coronal low-laying loops reconnection (Madjarska, 2019).
According to Young and Muglach (2014a) estimations, the jet is extended for H = 30 Mm with a width of ∆ℓ = 15 Mm. The jet duration is 25 min and the bright point is not significantly disrupted by the jet occurrence. The jet n is n i = 1.7 × 10 8 cm −3 , while that of the surrounding coronal plasma we assume to be n e = 1.5 × 10 8 cm −3 . The jet temperature is T i = 1.7 MK and the environment one is T e = 2.0 MK. The jet axial velocity is U z = 250 km s −1 and the rotational one is U φ = 90 km s −1 . Assuming a magnetic field twist ε 1 = 0.025 and B e = 3 G, from equation (3), we obtain η = 0.882, v Ai = 494.7 km s −1 (Alfvén speed in the environment is v Ae = 534.0 km s −1 ), and b = 1.014. We note, that while in the derivation of equation (11) we have neglected the thermal pressures, here, in using equation (3), we kept them. If we anticipate that the shortest unstable wavelength is equal to 7.5 Mm (with k z a = 2π), the mode number m whose instability range would accommodate the aforementioned wavelength, according to inequality (15) must be at least |m| = 71. The numerics show that the suitable m is |m| = 65. Thus, the input parameters for obtaining the numerical solutions to equation (11)  It is seen from the left panel that the unstable m = 65 MHD waves are generally super-Alfvénic. Since the instability developing times of the m = 65 mode are relatively short, that is, much less than the jet lifetime of 25 min, we can conclude that the KHI in this jet is relatively fast.
With the increase in the parameter ε 1 , the instability region, as seen from the right panel of Figure 6, becomes narrower and at the lower limit (k z a) lhs = π/2 with ε cr 1 = 0.10682 and M A = 0.5026 its width is equal to zero. In other words, there is no longer instability. Therefore, the critical azimuthal magnetic field that suppresses the KHI is B cr iφ ∼ = 0.3 G-obviously a relatively small value. Filippov et al. (2015) observationally studied three jets events originated from the active region NOAA 11715 (located on the west limb) on 2013 April 10-11. These authors claim that the jets originated from the emergence of a filament having a null-point (inverted Y) topology. We have considered the second event described in that paper for a detailed study. The jet electron number density, n i , and electron temperature, T i , both listed in Table 1, have been calculated by us using the techniques elaborated by Aschwanden et al. (2013). This technique requires the data from the six 94, 131,171,193,211, and 335 Å AIA/SDO EUV channels. In addition to the electron number densities and electron temperatures in the jet and surrounding plasma, we have also estimated the jet width as ∆ℓ≈30 Mm, its height as H = 180 Mm, and have found the jet lifetime to be 30 min. The two important parameters, axial and azimuthal velocities, according to the observations, are U z = 100 and U φ = 180 km s −1 , respectively. The time evolution of the jet in AIA 304 Å is shown in Figure 7 and we have observed vortex type structures in the eastern side of the jet, which are indicated by arrows. These structures implicitly indicate for occurrence of KHI.

Kelvin-Helmholtz Instability in a Jet Emerging from a Filament Eruption
With the typical n and T e (see Table 1), rotating velocity U φ = 180 km s −1 , assumed B e = 6 G, and ε 1 = 0.1, equation (3) yields η = 0.864, b = 4.36, and v Ai = 44.00 km s −1 (for comparison, the Alfvén speed in the environment is v Ae = 206.3 km s −1 ). We note that the choice of ε 1 was made taking into account the fact that the inclination of the treads of the jet in the event on 2013 April 10, detected by SDO/AIA, yields a relationship between B iφ and B iz , which was evaluated as ε 1 ≈ 0.1. If we assume that the shortest unstable wavelength is λ KH = 12 Mm, which is located at k z a = 2.5π on the k z aaxis, from inequality (15) we find that a MHD wave with |m| = 12 would provide an instability region, accommodating the nondimensional k z a = 2.5π. It turns out that a suitable mode number is m = 10. The wave dispersion and growth rate diagrams are shown in Figure 8. In that instability range one can also find the instability characteristics at k z a = 2π, which corresponds to λ KH = 15 Mm. The input parameters for finding the solutions to equation (10)  The KHI developing or growth times seem reasonable and the wave phase velocities are super-Alfvénic ones.
It is intriguing to see how the width of the instability range will shorten as the magnetic field twist ε 1 is increased. Our numerical computations indicate that for a noticeable contraction of the instability window one should change the magnitude of ε 1 with relatively large steps. The results of such computations are illustrated in Figure 9. It is necessary to underline that at values of ε 1 close to 1, (i) β i becomes <1 and the jet has to be treated as a cool medium, which implies a new wave dispersion relation and probably a higher wave mode number, m; (ii) one cannot use ε 1 > 1, because in that case the instability is of another kind, namely kink instability (Lundquist, 1951;Hood and Priest, 1979;Zaqarashvili et al., 2014). At this 'pathological' case, one cannot reach the lower limit of the instability range, (k z a) lhs = 0.524, and consequently we are unable to evaluate that azimuthal magnetic field, B cr φ , which will stop the KHI onset!

Kelvin-Helmholtz Instability in a Spinning Macrospicule
As we have mentioned in Section 1, Pike and Mason (1998) did a statistical study of the dynamics of solar transition region features, like macrospicules. These features were observed on the solar disk and also on the solar limb by using data from the Coronal Diagnostic Spectrometer (CDS) onboard SOHO. In addition, in their article, Pike and Mason (1998) discussed the unique CDS observations of a macrospicule first reported by Pike and Harrison (1997) along with their own (Pike and Mason) observations from the Normal Incidence Spectrometer (NIS). This spectrometer covers the wavelength range from 307 to 379 Å and that from 513 to 633 Å using a microchannel plate and CCD combination detector. The details of macrospicule events observed near the limb are given in Table I in Pike and Mason, 1998, while those of macrospicule events observed on the disk are presented in Table II. The main finding in the study of Pike and Mason (1998) was the rotation in these features based on the red and blue shifted emission on either side of the macrospicule axes. According to the authors, the detected rotation assuredly plays an important role in the dynamics of the transition region. Using the basic observational parameters obtained by Pike and Mason (1998), Zhelyazkov and Chandra (2019) examined the conditions for KHI rising in the macrospicule. Let us discuss that study as follows.
One observes that at both unstable wavelengths the corresponding phase velocities are supper-Alfvénic. Moreover, the two growth times of 2.2 and ∼0.6 min seem reasonable bearing in mind the fact that the macrospicule lifetime is about 15 minutes, which implies that the KHI at the selected wavelengths is rather fast. The B cr iφ that suppresses the KHI onset, equals 0.57 G and was calculated with ε cr 1 = 0.202085 and M A = 1.2119. Our study (Zhelyazkov and Chandra, 2019) shows that a decrease in the background magnetic field to B e = 4.8 G would require the excitation of MHD wave with mode number m = 48, at which the KHI characteristics at the wavelengths of 3 and 5 Mm are very close to those obtained with m = 52.

SUMMARY AND OUTLOOK
In this article, we have studied the emerging of KHI in four different spinning solar jets (standard and blowout coronal hole jets, jet emerging from a filament eruption, and rotating macrospicule) due to the excitation of high-mode (m 2) MHD waves traveling along the jets. First and foremost, we model each jet as a vertically moving with velocity U cylindrical twisted magnetic flux tube with radius a. There are four basic steps in the modeling as follows: • Topology of jet-environment magnetic and velocity fields For simplicity, we assume that the plasma densities of the jet and its environment, ρ i and ρ e , respectively, are homogeneous. Generally they are different and the density contrast is characterized by the ratio ρ e /ρ i = η. The twisted internal magnetic and velocity fields are supposed to be uniform, that is, represented in cylindrical coordinates, (r, φ, z), by the vectors B i = (0, B iφ (r), B iz ) and U = (0, U φ (r), U z ), where their azimuthal components are considered to be linear functions of the radial position r, viz. B iφ (r) = Ar and U φ (r) = Ωr, where A and Ω (the azimuthal jet velocity) are constants. We note that B iz and U z are also constants. It is convenient the twists of the magnetic field and the flow velocity of the jet to be characterized by the two numbers ε 1 = B iφ (a)/B iz ≡ Aa/B iz and ε 2 = U φ (a)/U z ≡ Ωa/U z , respectively. Note that Ωa is the jet rotational speed U z . The surrounding coronal or chromospheric plasma is assumed to be immobile and embedded in a constant magnetic field B e = (0, 0, B e ). In our study, the density contrast, η, varies from 0.1 to 0.9, the magnetic field twist, ε 1 , can have a wide range of magnitudes from 0.005 to 0.95 (it has to be less than 1 in order to avoid the rising of the kink instability), while the velocity twist parameter, ε 2 , is fixed by the observationally measured rotational and axial speeds.
• Listing of the basic physical parameters and determination of plasmas betas In general, at a fixed density contrast, the plasma beta is controlled by the magnetic field (inside or outside the magnetic flux tube) and the electron temperatures of the jet and surrounding plasma. The values of these physical parameters should satisfy the total pressure balance equation (3) at all levels (equilibrium and perturbational). Our practice is to fix B e , and using equation (3) to determine the internal Alfvén speed defined as v Ai = B iz / √ µρ i . It is worth underlying that the usage of equation (3) requires the specification of ε 1 . In our four cases, for finding the KHI characteristics, we took ε 1 to be equal to 0.005, 0.025, or 0.1. The electron temperatures in the jets are from 500 000 K for the macrospicule to 2.0 MK in the jet emerging from a filament eruption. The electron temperatures of surrounding plasmas are 2.14 MK in the active solar region (Filippov et al., 2015), 2.0 MK in the blowout coronal hole jet (Young and Muglach, 2014a), and 1.0 MK in the environments of the standard coronal hole jet (Chen, Zhang, and Ma, 2012) and the macrospicule. With background magnetic fields of 3 to 6 G, rotating velocities of 40 to 180 km s −1 , and ε 1 = 0.005, the total pressure balance equation (3) yields plasma betas of the first, third, and forth jets greater than 1 and those of the environments and the internal medium of the second jet much less than 1 (see Table 1). With these plasma beta values one can consider the media of the standard coronal hole jet, the rotating jet emerging from a filament eruption, and microspicule as nearly incompressible plasmas, while the internal medium of the blowout coronal hole jet and the surrounding magnetized plasma in the four cases can be treated as a cool medium (Zank and Matthaeus, 1993).
• Solving of the wave dispersion relation and finding the KHI characteristics For finding the solutions to the MHD wave dispersion equations (10) or (11), which are a slight modification of the 'basic' dispersion relation derived in Zaqarashvili, Zhelyazkov, and Ofman, 2015, it is necessary to specify the following input data: the wave mode number, m, the density contrast, η, the two twist parameters ε 1 and ε 2 , the magnetic fields ratio, b (obtainable from equation (3)), and the Alfvén Mach number The roots of the dispersion equations (10) or (11) are the normalized wave phase velocity and instability growth rate as functions of the nondimensional wavenumber k z a. From plots which graphically represent the found solutions, one obtains at the anticipated wavelengths (given by their k z a-values on the horizontal k z a-axes) the corresponding Re(v ph /v Ai ) and Im(v ph /v Ai ) values. From them, one can find in absolute units the KHI growth rate, γ KH , the instability developing or growth time, τ KH = 2π/γ KH , and the corresponding wave phase velocity, v ph . The MHD wave mode numbers at which we were able to calculate the instability characteristics at wavelengths comparable to the radius or width of the jet are between 10 and 65 and the KHI growth times at those wavelengths are of the order on a few minutes, generally in good agreement with the observations. It is curious to note that in searching KHI growth times of the order on few seconds, when studying the dynamics and stability of small-scale rapid redshifted and blueshifted excursions, appearing as high-speed jets in the wings of the Hα line, Kuridze et al. (2016) had to assume the excitation of MHD waves with mode numbers up to 100. A typical property of the instability developing times, owing to the shape of the plotted dispersion curves, is that with increasing the examined wavelength the growth times become shorter-for instance, at λ KH = 10 Mm the KHI developing time in the coronal hole jet (Chen, Zhang, and Ma, 2012) is around 4.5 min, while at λ KH = 12 Mm it is equal to ∼ =2.1 min. A change in B e can influence the MHD mode number m, which would yield an instability region similar or identical to that seen in Figure 2. It is necessary to mention that the width of the instability range except by changing the MHD wave mode number, m, can be regulated by increasing or diminishing the parameter ε 1 .
• Finding the critical azimuthal magnetic field which suppresses the emergence of KHI It was numerically established, that any increase in ε 1 yields to the shortening of the instability range. This observation implies that there should exist some critical ε 1 , at which the upper limit of the instability range coincides with the lower one-in that case the width of the instability window is zero, which means that there is no longer any instability. With such an ε cr 1 , one can calculate the B cr iφ , which stops the KHI appearance. For the rotating blowout coronal hole jet this B cr iφ is relatively small-it is equal to ≈0.3 G, while for the standard coronal hole jet it is 1.4 G. It is worth noticing that due to the specific parameters of the jet emerging from a filament eruption, we were unable to find a B cr iφ , which would stop the KHI onset because with values of ε 1 close to 1 our dispersion relation becomes inappropriate (the internal medium being nearly incompressible becomes a cool one) and we cannot calculate that ε 1 at which the lower limit, (k z a) lhs = 0.524, can be reached (see Figure 9).
In this article, we also corrected the total pressure balance equation used in Zhelyazkov et al., 2018a andChandra, 2018b, which turns out to be erroneous. The true total balance equation is given by equation (3). In fact, the corrected pressure balance equation, used here, changes the mode numbers at which the KHI occurs, namely from m = 12 to m = 11 for the coronal hole jet and from m = 18 to m = 10 for the rotating jet emerging from a filament eruption. The computed KHI developing or growth times in the aforementioned articles, nonetheless, are not changed noticeably-they are of the same order with these computed in this paper. In addition, there is another improvement issue, scilicet when studying how the increasing ε 1 shortens the instability region, one has to apply equation (3)  Our approach in investigating the KHI in rotating twisted solar jets can be improved in the following directions: (i) to assume some radial profile of the plasma density of the jet, which immediately will require additional study on the occurrence of continuous spectra and resonant wave absorption (Goedbloed and Poedts, 2004), alongside to see to what extent these phenomena will influence the instability growth times; (ii) to investigate the impact of the nonlinear azimuthal magnetic and velocity fields radial profiles on the emergence of KHI; and (iii) to derive a MHD wave dispersion relation without any simplifications like considering the jet and its environment as incompressible or cool plasmas-this will show how the compressibility will change the picture. We should also not forget that the nonlinearity, as (Miura and Pritchett, 1982;Miura, 1984) claim, can lead to the saturation of the KHI growth, and to formation of nonlinear waves. Nevertheless, even in its relatively simple form, our way of investigating the conditions under which the KHI develops is flexible enough to explore that event in any rotating solar jet in case that the basic physical and geometry parameters of the jet are provided by observations.

CONFLICT OF INTEREST STATEMENT
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

AUTHOR CONTRIBUTIONS
IZh wrote the substantial parts of the manuscript. RJ wrote the Introduction section and prepared three figures associated with the observations. RC contributed by writing the parts devoted to observations, as well as in the careful proofreading of the text.

FUNDING
The work of IZh and RC was supported by the Bulgarian Science Fund contract DNTS/INDIA 01/7. RJ was funded by the Department of Science and Technology, New Delhi, India as an INSPIRE fellow. RC was also supported from the SERB-DST project no. SERB/F/7455/2017-17.  Growth rates of the unstable m = 11 mode for the same input parameters. The azimuthal magnetic field that corresponds to ε cr 1 = 0.653577 (the instability window with zero width) and stops the KHI onset is equal to 1.4 G.   Im(v ph /v Ai ) k z a ε 2 = 1.8 m = 10 η = 0.864 Figure 9. Growth rate curves of the unstable m = 10 MHD mode propagating along a twisted incompressible emerging from a filament eruption jet at η = 0.864, ε 2 = 1.8, and the following values of ε 1 (from right to left): 0.1 (orange), 0.6 (green), 0.8 (purple), and 0.9 (blue) .