Corrugation Instability of a Coronal Arcade

We analyse the behaviour of linear magnetohydrodynamic perturbations of a coronal arcade modelled by a half-cylinder with an azimuthal magnetic field and nonuniform radial profiles of the plasma pressure, temperature, and the field. Attention is paid to the perturbations with short longitudinal (in the direction along the arcade) wavelengths. The radial structure of the perturbations, either oscillatory or evanescent, is prescribed by the radial profiles of the equilibrium quantities. Conditions for the corrugation instability of the arcade are determined. It is established that the instability growth rate increases with decreases in the longitudinal wavelength and the radial wave number. In the unstable mode, the radial perturbations of the magnetic field are stronger than the longitudinal perturbations, creating an almost circularly corrugated rippling of the arcade in the longitudinal direction. For coronal conditions, the growth time of the instability is shorter than one minute, decreasing with an increase in the temperature. Implications of the developed theory for the dynamics of coronal active regions are discussed.


Introduction
Arcades of plasma loops are a typical feature of solar coronal active regions. Examples of coronal arcades are sets of postflare loops in two-ribbon flares and closed magnetic configurations overlying eruptive magnetic ropes. Arcades are key ingredients of the standard D.Y. Klimushkin et al. model of a flare and an inverse-polarity prominence equilibrium. Dynamical processes in arcade structures are important for understanding several key basic processes in the corona, such as energy balance and release, and prominence formation and evolution. Often, these processes are observed to evolve quasi-periodically. For example, light curves of solar and stellar flares often have quasi-periodic pulsations that could be associated with oscillatory processes (e.g. Nakariakov and Melnikov, 2009) in the flaring arcades. A transverse oscillatory motion with a period of about 1000 s in an arcade situated above an emerging magnetic rope was detected by Kim, Nakariakov, and Cho (2014). Consecutive ignition of quasi-periodic hard X-ray emission with a period of about a minute, propagating along the magnetic polarity inversion line in a flaring arcade, was observed in a solar flare by Grigis and Benz (2005) and Bogachev et al. (2005). A similar effect was detected in the microwave emission of a solar flare (e.g. Kim et al., 2013). A recent comprehensive statistical analysis demonstrated that phenomenon is associated with successive triggering of the energy release process in different loops of the flaring arcade (Kuznetsov et al., 2016). Oscillatory processes in the magnetic arcades supporting quiescent prominences might also be responsible for prominence oscillations (Arregui, Oliver, and Ballester, 2012).
A likely mechanism for the oscillatory processes in arcades is connected with magnetohydrodynamic (MHD) oscillations that have been intensively modelled for more than two decades (see, e.g. Goossens, Erdélyi, and Ruderman, 2011, for a comprehensive review). Several studies have addressed perturbations in the plane perpendicular to the arcade axis. In particular, the propagation of fast magnetoacoustic waves in coronal potential arcades, excited by photospheric footpoint motions, was studied in Cadez, Oliver, and Ballester (1996). The excitation of resonant shear Alfvén oscillations in coronal arcades driven by toroidal footpoint motions was described by Ruderman et al. (1997). Oliver, Hood, and Priest (1996) considered a general two-dimensional arcade equilibrium structure with no longitudinal magnetic field component, and derived governing equations for the coupled slow and fast modes, while the Alfvén mode was decoupled from the magnetoacoustic modes. The propagation of impulsively generated Alfvén waves in an arcade was numerically modelled by Chmielewski et al. (2014). Magnetohydrodynamic waves in sheared coronal arcades, including propagation along the axis of the arcade, were studied by e.g. Arregui, Oliver, and Ballester (2004b,a) in the zero-β limit, and the properties of coupled fast and Alfvén modes were established. Slow waves propagating along the axis of an arcade, i.e. along the neutral line, have been proposed as a candidate process to explain the progression of flaring energy releases along the neutral line that are typical for two-ribbon flares (Nakariakov and Zimovets, 2011). Fast waves guided by a plasma arcade have been considered as a possible mechanism for kink oscillations in coronal active regions (Hindman and Jain, 2015).
Much less attention has been received by MHD instabilities in plasma arcades, which may also be responsible for the progression of flaring energy releases along the neutral line in two-ribbon flares. Another interesting topic that requires understanding regarding MHD instabilities in arcades is the abrupt destabilisation that could trigger a coronal mass ejection. Likewise, the source of the filamentation of an arcade, i.e. the appearance of corrugated fine structure that is often detected in post-flare plasma arcades, is of interest. These processes require consideration of MHD perturbations that have a characteristic scale along the arcade axis (i.e., across the field) much shorter than the lengths of the loops that form the arcade. These perturbations, with a spatial scale much shorter than the parallel field-aligned scale when the field is azimuthal, could be referred to as corrugated perturbations, and the associated instability could be called a corrugation instability. On the other hand, as these perturbations have a small scale across the field, they could also be referred to as "ballooning", but this term must be used with caution in the coronal context, as this term is usually applied in solar physics to perturbations with a high azimuthal wave number m.
In solar physics, only a few studies have been dedicated to ballooning modes so far. These works were mainly dedicated to the study of slightly twisted magnetic ropes, where the ballooning perturbations were associated with high values of m. Ruderman (2007) carried out a theoretical analysis of long-wavelength non-axisymmetric modes, including ballooning modes, of a compressible zero-beta weakly twisted magnetic rope surrounded by a straight and homogeneous magnetic field. The main attention was paid to the kink mode. Karami and Bahari (2012) generalised this study and addressed the oscillation period ratios of different harmonics. González (2006, 2008) developed an energy-principlebased formalism for the estimation of the non-axisymmetric mode parameters, i.e. m > 0.  showed that in a thin magnetic tube with a curved magnetic field, the stable m = 1 mode is dominant when the azimuthal magnetic field outside the tube is absent. However, in the presence of an external magnetic field, this mode becomes unstable (Cheremnykh, Kryshtal, and Tkachenko, 2017). Stepanov, Kopylova, and Tsap (2008) and Tsap et al. (2008) demonstrated that the ballooning instability of arch-shaped twisted magnetic ropes might cause the formation of helmet-shaped structures in the lower solar corona, and might also produce a quasi-periodic modulation of gyrosynchrotron emission in flares.
On the other hand, the theory of the ballooning oscillations and instability of highly twisted magnetic ropes, relevant to the study of the corrugation instability of an arcade, has been intensively elaborated in the context of the terrestrial magnetosphere. Gold (1959) suggested very early on that the hot plasma MHD instability (later coined the ballooning instability) can explain various eruptive phenomena in the terrestrial magnetotail. Walker (1987) showed that the MHD modes in the ballooning limit are represented by the coupled Alfvén and slow magnetosonic branches, while the fast magnetosonic mode is evanescent. Burdo, Cheremnykh, and Verkhoglyadova (2000) showed that the coupling of the Alfvén and slow magnetoacoustic modes gives rise to a ballooning instability in the dipole field. These authors also formulated for the first time the instability criterion outside the local approximation. Cheremnykh and Parnowski (2004), Agapitov, Cheremnykh, and Parnowski (2008) and Leonovich and Kozlov (2013b) considered the instability in a plasma with an inhomogeneity along the field lines. Ma, Hirose, and Liu (2014) studied the instability in the Hall-MHD framework. Klimushkin, Mager, and Pilipenko (2012) and Mager and Klimushkin (2017) generalised the theory of ballooning perturbations accounting for kinetic effects. The exploitation of the coronal-magnetospheric analogy opens up interesting opportunities for knowledge transfer (Nakariakov et al., 2016).
The aim of this paper is to further develop the linear theory of corrugated perturbations in a plasma arcade with an arbitrary radial profile of the equilibrium quantities. The shear of the magnetic field is neglected. Mathematically, this work is based on the set of equations for the coupled Alfvén and slow magnetoacoustic modes derived in Klimushkin (1997Klimushkin ( , 1998 and Cheremnykh and Danilova (2011), and the formalism developed for the study of the spatial structure of the unstable modes in the terrestrial magnetosphere (Klimushkin, 2006;Klimushkin and Chen, 2006;Mazur, Fedorov, and Pilipenko, 2013;Cheremnykh, Klimushkin, and Kostarev, 2014;Cheremnykh, Klimushkin, and Mager, 2016). In Section 2 we review the derivation of the governing equations that describe the radial structure of the perturbations. Section 4 addresses the derivation and study of the dispersion relation for the coupled Alfvén and slow modes with large longitudinal wave numbers perpendicular to the field, which determines the wave vector radial component in the WKB approximation. In Section 5 we study the possibility of mode trapping in the radial direction across the magnetic shells, derive the eigenvalues of the perturbations, consider the conditions under which the modes are unstable, and estimate the growth rates. In Section 6 we summarise the results and discuss their possible observational manifestations. D.Y. Klimushkin et al.

Figure 1
Model representing a twisted magnetic rope filled with a plasma, as a cylinder. The coordinates y, ϕ, and r are directed along the cylinder axis, along the magnetic field line, and in the radial direction, respectively. The horizontal line denotes the photosphere with the line-tying boundary conditions.

Governing Equations
Consider a coronal arcade as half of a plasma cylinder with an azimuthal magnetic field. The magnetic field component along the axis of the cylinder is neglected. The plasma is considered to be of finite β. Thus, in the magnetic geometry in our model corresponds to a z-pinch, with the magnetic surfaces represented by nested coaxial cylinders ( Figure 1). In the solar coronal context, a similar model with a purely poloidal magnetic field has recently been used by Kaneko et al. (2015). The azimuthal angle ϕ shows the direction along the field lines. At ϕ = 0, π/2, which corresponds to the photosphere, the line-tying boundary conditions are applied. The longitudinal y coordinate is directed along the cylinder axis, i.e. the magnetic neutral line. All parameters of the equilibrium are taken to depend only on the radial coordinate r, that is, the field line curvature radius. The coordinates ϕ, y, and r are locally orthogonal to each other. In the following, we focus on ideal MHD processes.
The plasma equilibrium condition has the form where P , J , and B are the equilibrium values of the plasma pressure, electric current density, and magnetic field, respectively. It is convenient to introduce the parameters κ = P −1 dP /dr, that is, the reciprocal spatial scale of the radial non-uniformity of the plasma pressure, and β = 8πP /B 2 , that is, the ratio of the plasma pressure to the magnetic pressure. With this notation, the equilibrium condition can be written as In the following, we consider β, κ, and B as smooth functions of the radial coordinate r.
Within the approximation of infinite plasma conductivity, the component of the wave electric field E along the ambient magnetic field is zero. In the ballooning approximation, a two-dimensional field E is non-rotational (e.g. Klimushkin, 1997), that is, it can be represented as a two-dimensional gradient of some scalar ("a potential"): As the other field variable, the quantity can be used, where ξ is the plasma displacement vector, and γ is the adiabatic index.
Owing to the assumed plasma uniformity with respect to the y direction, we can perform a Fourier transform of any perturbed quantity F as Here, the quantity k y is the component of the wave number along the axis of the arcade. The dependence of F on the radial coordinate r is determined by the specific equilibrium radial profiles of P , J , and B, and the boundary conditions at r = 0, ∞. In the following, we focus on the ballooning limit that is based on the assumption that the characteristic spatial scale of the perturbations across the field, i.e. in the y and r directions, are small in comparison with the spatial scales of the equilibrium quantities. A necessary condition for this is a high k y value: k y κ.
Considering that the perturbations with spatial scale in the radial direction are much smaller than the local radius of the field curvature, we can use the WKB approximation with respect to the radial coordinate r. That is, the r dependence of the perturbations can be represented in the form exp i k r (r, ω)dr, where k r is the radial component of the wave vector. The perpendicular component of the wave vector k ⊥ is k 2 ⊥ = (k 2 r + k 2 y ) 1/2 . A set of the MHD equations under these assumptions was obtained for an arbitrary geometry by Klimushkin (1997Klimushkin ( , 1998 and later, but independently, by Cheremnykh and Danilova (2011) and Mazur, Fedorov, and Pilipenko (2012). For the cylinder geometry, it can be written as Here V A , V S , and V C are the Alfvén, sound, and tube (cusp) speeds, defined as where ρ is the undisturbed plasma density. Equations 6 and 7 determine the structure of the Alfvén and slow modes coupled with each other through the field line curvature; we recall that in the limit k y κ, the fast mode is evanescent (e.g. Walker, 1987). The equations must be supplemented with boundary conditions. On the azimuthal coordinate, the line-tying boundary conditions must be applied (Hood, 1986), which means that the plasma displacement ξ must vanish on the photosphere, where ϕ = 0 and π : As the transverse displacement is proportional to the wave electric field, the condition ξ ⊥ = 0 results in The parallel along the field displacement is related to the function as Klimushkin et al. (e.g. Klimushkin, 1998). Then, the equality ξ = 0 imposes the boundary condition for the function: Equations 6 and 7 supplemented with boundary conditions 9 and 10 are the governing equations used for the analysis presented in the paper.

Alfvén and Slow-Mode Eigenfunctions
Equations 6 and 7 describe the coupling of Alfvén and slow modes. We consider two auxiliary problems connected with uncoupled modes.
We designate N a solution of the boundary problem stated by the Alfvén mode equation with the boundary conditions where N is the parallel (field-aligned) integer harmonic number. The function N is normalised as where ... = π 0 (...)dϕ. Eigenfunctions of problems 11 and 12 are and the corresponding eigenfrequencies are Likewise, we designate N as a solution of the boundary problem stated by the slowmode equation with the boundary conditions normalised as where N and N are two integers. Under these conditions, the eigenfunction is and the corresponding eigenfrequency is Note that the cases ω = AN and ω = CN correspond to well-known Alfvén and slowmode resonances, respectively (e.g. Cheng et al., 1993;Klimushkin, 1997Klimushkin, , 1998Leonovich and Kozlov, 2013a).

Dispersion Relations for the Coupled MHD Modes
We seek the solution of problems 6 and 7 with boundary conditions 9 and 10 as an expansion by harmonics N and N in the form where a N and b N are the coefficients that are yet to be determined. After we insert Equation 22 into Equation 7 and integrate over the angle, we have Then we substitute Equation 22 into Equation 6, which yields The symmetric and antisymmetric (with respect to the arcade apex (ϕ = π/2)) functions must be considered separately. In this paper, the convolution product N N is considered as non-negligible only for the lowest relevant N and N numbers.

Alfvén Mode Eigenfunctions Symmetric with Respect to ϕ = π/2
In this case, the only harmonics that considerably contribute to the dispersion relation are 0 , 2 , and 1 . Then Equation 25 is reduced to the form After we multiply this expression by 1 and integrate over the angle, we have the dispersion relation where expressions 15 and 20 were used with the corresponding values of the integer N . The case k r → ∞ corresponds to the Alfvén and slow-mode resonances that occur at the cylindrical surfaces where the wave frequencies are equal to A1 and C2 , respectively. The opposite case, k r → 0, corresponds to the cut-offs (e.g. Cheremnykh, Klimushkin, and Kostarev, 2014). The cut-off surfaces are determined by the equation where ± is a solution of biquadratic Equation 27 with k 2 ⊥ /k 2 y = 1: with and The value of D is always positive. An instability ( 2 − < 0) can develop on the lower frequency oscillation branch, when the root of the discriminant D is larger than the first term in brackets in Equation 28: which gives the instability threshold: Thus the arcade is unstable when The instability growth rate is The instability criterion 32 qualitatively agrees with the results of Burdo, Cheremnykh, and Verkhoglyadova (2000) that were obtained for the dipole field. Equation 32 should be considered an approximate threshold of the instability. The exact threshold can be found from Equation 25 with ω = 0 and k 2 ⊥ /k 2 y = 1 and boundary conditions 9, where 2m−1 are odd eigenfunctions N given by Equation 14. In this case, we find the instability condition The exact solution of this equation corresponding to the minimum value of −βκr is compared with approximate solution 31 in Figure 2. The approximate threshold of the instability, given by Equation 31, differs only slightly from the exact threshold, especially for V 2 C /V 2 A 1, i.e. for the small β typical for coronal arcades.

Alfvén Mode Eigenfunctions Antisymmetric with Respect to ϕ = π/2
The principal, antisymmetric with respect to the apex harmonic of the Alfvén mode, is 2 . Thus, the convolution N 2 is not zero only for the odd slow-mode harmonics N (with N being an odd integer). We only take into account the principal odd harmonic, 1 . Then Equation 25 becomes (37) D.Y. Klimushkin et al. Multiplying this expression by 2 , and integrating over the angle, we obtain the dispersion relation The Alfvén and slow-mode resonances correspond to the case k r → ∞ and occur at the frequencies A2 and C1 , respectively. The cut-off frequencies, determined by the condition k r → 0, are in this case determined as with and The instability can develop on the lower-frequency branch given by the solution of Equation 39 denoted as "−" when the condition is satisfied. The corresponding instability threshold is in complete agreement with the results of Burdo, Cheremnykh, and Verkhoglyadova (2000). The instability growth rate is As in the symmetric case, the exact instability threshold can be found from Equation 34, while now = 2m , where m = 1, 2, . . . and 2m are even eigenfunctions given by Equation 14. Thus the instability condition is Here −βκr is minimum when m = 1, so in the antisymmetric case, the exact instability threshold coincides with the approximate threshold 42.

General Properties of the Solutions for Symmetric and Antisymmetric Alfvén Modes
In both symmetric and antisymmetric cases, the discriminant D is a positive value, thus the squares of the cut-off frequencies have real values. This agrees with the well-known MHD theorem that is a consequence of the energetic principle, stating that the wave frequency squared of linear perturbations in ideal MHD is always real (Frieman and Rotenberg, 1960;Kadomtsev, 1966). Hence, no over-stabilities can occur in the considered case. With the decrease in β, the lower cut-off frequency − in both cases tends to the slowmode resonance frequency ( C2 or C1 ), while the higher cut-off frequency + tends to the Alfvénic resonant frequency ( A1 or A2 ). Thus, we can conclude that the frequencies − and C characterise the slow-mode oscillation branch, and + and A characterise the Alfvénic oscillation branch. The instability can develop only on the slow-mode branch. This is in accordance with the results of earlier papers (Mazur, Fedorov, and Pilipenko, 2012;Kozlov et al., 2014;Cheremnykh, Klimushkin, and Kostarev, 2014). However, the instability threshold depends on the symmetry properties of the Alfvén mode. Clearly, for low and moderate β values, the instability threshold is lower for the symmetric Alfvén mode than for the antisymmetric mode.
Then, Equations 28 and 39 easily show that the inequalities − < C < + are satisfied. However, the frequency + can be either higher or lower than the Alfvénic resonant frequency A , depending on the profiles β(r), κ(r), and B(r). Following Cheremnykh, Klimushkin, and Kostarev (2014), for the arbitrary ratio k r /k y , the solution of Equations 27 and 38 can be expressed in terms of two pairs of the characteristic frequencies as In the symmetric case, C means C2 , and A means A1 . In the antisymmetric case, C means C1 , and A means A2 . The dependence of the equilibrium parameters on the radial coordinate shows that there are regions of different signs of k 2 r in the radial direction. The regions where k 2 r (r) > 0 can be called the mode localisation regions. In these regions, MHD perturbations have an oscillatory radial structure, cf. the body modes, see Edwin and Roberts (1983), or the leaky modes (e.g. Verwichte, Foullon, and Nakariakov, 2006). Outside these regions, the modes have exponential radial dependences. The cut-off frequencies ± prescribe the cylindrical surfaces where the behaviour of the perturbation changes from oscillatory to evanescent in the radial direction. Equation 45 shows that there are mode localisation regions of two kinds. One of them is bounded by the Alfvénic resonance surface where ω = A (r) and the higher cut-off shell where ω = + (r). In the following, this region is referred to as the Alfvénic mode localisation region. The other region is bounded by the surface of the cusp resonance where ω = C (r) and the lower cut-off shell where ω = − (r). This is the slow-mode localisation region. Because of the inequality − < A < + , these regions do not overlap with each other.
The spatial structure of the perturbation in the mode localisation regions in the considered geometry was studied by Klimushkin (1997) and Cheremnykh, Klimushkin, and Mager (2016). In both Alfvénic and slow-mode localisation regions, the wave was found to have an oscillatory structure across the magnetic shells in the radial direction in the annulus between the cut-off shell and the resonant shell. In such an annulus, the local radial component of the wave vector changes from zero to infinity. On the resonant magnetic surface, the mode is absorbed through mode coupling and subsequent dissipation. Outside these regions, the mode is evanescent. However, this structure corresponds to a stable plasma configuration when all the cut-off and resonant frequencies have real values. However, the lower cut-off frequency D.Y. Klimushkin et al.

Figure 3
Behaviour of the function 2 − (r) in the vicinity of its minimum.

Radial Structure of the Unstable Modes
The instability takes place when ω 2 < 0 for a given k 2 r value. At the resonances, k 2 r → ∞, and the modes of the considered plasma configuration are always stable since the resonant frequencies squared are always positive. At the cut-off frequencies, the situation is different. Equations 28 and 39 show that the Alfvénic (the higher) cut-off frequency squared 2 + is always positive, while the slow-mode (the lower) cut-off frequency squared 2 − can be negative, provided that either condition 32 or 42 is fulfilled. This is the instability condition. Only the slow-mode oscillation branch can be unstable, in agreement with Mazur, Fedorov, and Pilipenko (2012).
In the low-β case, the characteristic Alfvénic frequencies A , + are much higher than the slow-mode characteristic frequencies − , C . Then, in the low-β limit, in the slowmode localisation region Equation 45 can be simplified to where K = ( + / A )k y . As the values of + and A are close to each other, the value of K is of the same order of magnitude as k y . As the instability can develop at low (moreover, negative) values of 2 − (r), we assume that the function 2 − (r) has a minimum at some value of r = r min (Figure 3). Near the magnetic surface r = r min , the function 2 − (r) can be approximated with the use of the Taylor expansion as where x = r − r min , 2 = − 2 − (r min ), and l is a constant linked with the characteristic scale L of the radial non-uniformity (see Figure 3). The region where 2 − < 0 is referred to as an unstable region. Equation 47 shows that the width of this region is approximately 2l.
The value is the characteristic growth rate of the instability. It can be estimated from Equations 28 and 39 as where it was taken that the field line curvature radius at r = r min is r min . As in the unstable case the wave frequency is imaginary, so we can put ω 2 = −γ 2 , where γ 2 > 0. Then Equation 46 can be rewritten as We introduce 2 = l 2 2 − γ 2 2 (50) Figure 4 Dependence of the normalised growth rate squared on the parameter ζ = (2n + 1)/(Kl) for the values 2 C / 2 = 0.5, 1, and 2.
(note that is a function of the growth rate γ ). Then Equation 49 becomes This expression shows that is the radial coordinate of the cut-off magnetic shell. That is, the mode is trapped in the radial direction in the region bounded between the magnetic surfaces with r ± = r min ± (cut-off points). As we are interested mainly in the order of magnitude estimations, the variation of the C across the trapping region can be neglected.
In the cavity (i.e. the wave-trapping region) formed by the radial non-uniformity of the equilibrium parameters, the radial component of the wave vector must satisfy the Bohr-Sommerfeld quantisation condition r+ r− k r (r, γ )dr = π n + 1 2 , where n is an integer corresponding to the radial harmonic number. Substituting k r from Equation 51 into Equation 52 and evaluating the integral, we obtain after some rearrangement where γ n is the growth rate of the n-th radial harmonic. The solution to this equation has the form where ζ = 2n + 1 Kl , and D 1 = 1 + 1 2 The dependence of the growth rate on the parameter ζ is illustrated in Figure 4. Note that the discriminant D 1 is always positive. The figure shows that the perturbations with short wavelengths in the y-direction, along the axis of the arcade, are most unstable, while the perturbations with long wavelengths along the axis of the cylinder are stable. D.Y. Klimushkin et al.

Figure 5
Sketch of an arcade experiencing a corrugation instability. The arrows show the unstable magnetic flux tubesthe "tongues".
In the unstable regime, for corresponding to the perturbations with very short wavelengths along the axis of the arcade and small n, the solution to Equation 53 can be approximated as Thus, for low radial harmonics (n ≈ 1) and Kl 1, the quantity γ 2 n has the same (positive) sign as 2 , and the instability occurs. Therefore, the set of inequalities 32 or 42 and 55 constitutes the corrugation instability condition. In this instability, the wavelength of the unstable perturbation along the arcade is much shorter than its radius, as the unstable perturbations have a "corrugated" spatial structure, see Figure 5. However, there is yet another condition for the instability.
Different radial harmonics have different growth rates. By the order of magnitude, the difference between the increments of two neighbouring harmonics with small n is where R is the radius of the arcade. For larger radial harmonic numbers n, however, the value γ 2 n becomes negative and the perturbations become stable. Thus, the instability conditions 32 or 42 and 55 must be complemented by the inequality n < n cr , where n cr is the critical value of the radial harmonic number, which can be estimated from Equation 54 as This stability condition has a rather simple meaning.

First, Equation 50
shows that in the unstable case γ 2 n > 0, the inequality < l holds. As follows from Equation 47, this regime occurs when the cut-off shells are situated in the region where the slow-mode cut-off frequency squared is negative, 2 − < 0. Vice versa, the stable regime, with γ 2 n < 0, occurs when the cut-off shells are situated outside the unstable region with 2 − < 0. That is, in the radial direction, only the central part of the resonator is situated inside the unstable region: > l.
Second, the characteristic value of the radial component of the wave vector for the critical radial number n cr can be estimated as k r ≈ n cr /l. From Equation 58, we can estimate the ratio of the radial and longitudinal components, Since both ≈ C and + ≈ A , for the stable modes with n n cr we have k r k y . This allows us to estimate the ratio of transverse components of the magnetic field perturbations in this regime, |b r /b y | ≈ k y /k r . Then, in the stable case, the field lines oscillate mainly in the direction along the axis of the arcade, b r b y , while in the unstable case, they are displaced mainly in the radial direction, b r b y . This agrees with results obtained in earlier publications (Agapitov, Cheremnykh, and Parnowski, 2008;Mazur, Fedorov, and Pilipenko, 2012). Thus, the corrugation instability of an arcade distorts its boundary in the radial, i.e. vertical direction. The unstable region has a small size in the direction of the arcade axis, while it includes the whole magnetic flux tube between the footpoints, see Figure 5. A perturbation of this kind could be called a "tongue".

Discussion and Conclusions
Waves guided by various plasma non-uniformities are often considered as a superposition of eigenmodes of certain spatially localised resonators where the perturbations have an oscillatory structure between some reflecting or refractive walls and are evanescent outside (e.g. Nakariakov et al., 2016;Cheremnykh, Klimushkin, and Mager, 2016). In this paper, this approach was applied to the study of unstable small-scale corrugation (or ballooning) perturbations of a plasma arcade filled in with a radially non-uniform plasma of finite temperature. This allowed us to understand the influence of the radial wave structure on the instability of an arcade.
We considered an arcade as a half-cylinder with line-tying boundary conditions at the bottom of the corona. The effects of the magnetic shear were neglected, while the dependence of the perturbation along the axis of the arcade was taken into account. We focused on perturbations with wavelengths along the arcade much smaller than the radius of the arcade. We found that this plasma configuration may be unstable to the corrugated perturbations of the slow magnetoacoustic nature, which is mainly characterised by the radial displacement of the plasma. The arcade is most unstable to the perturbations with low radial harmonic numbers, n ≈ 1, and very short longitudinal (along the cylinder) wavelengths, in comparison with the radius of the rope. In the azimuthal direction, along the magnetic field lines, the wavelength is of the order of the field line length. Thus, the unstable perturbations cause displacements of the plasma in the radial direction, and have an almost circularly corrugated structure in the direction along the neutral line. The discussed instability is most effective in hot arcades with non-negligible values of the plasma parameter β.
The characteristic time of the instability was determined approximately by the ratio of the arcade radius and the sound speed (see Equation 48), i.e. the acoustic travel time across the arcade. For example, in an arcade of radius 50 Mm and filled with the plasma of the temperature of 1 MK, the characteristic time of the instability is estimated as about one minute. For hotter arcades, this time would be shorter, e.g. about twenty seconds, if the temperature were 10 MK. We should stress that the term "corrugated" highlights the ringlike localised (in the longitudinal direction) perturbations of the arcade. In an arcade, these perturbations are of an interchange nature, and hence could be considered to be similar to the fluting, or ballooning, perturbations of a plasma cylinder with the magnetic field along its axis (e.g. considered in Edwin and Roberts, 1983).
With the increase in the radial harmonic number, the ratio of the radial and longitudinal (along the neutral line) components of the magnetic field decreases, as does the instability growth rate. When this ratio approaches unity, the instability ceases. The full set of the instability conditions is determined by inequalities 55 and 57, which prescribe the fieldaligned (azimuthal), longitudinal, and radial structure of the unstable corrugated mode.
As in the unstable modes the magnetic field and plasma are displaced mainly in the radial direction, the instability may lead to the formation of tongues that are highly localised in the direction along the arcade, displacing practically a whole line-tied magnetic flux tube in the radial direction. The magnetic field frozen into a plasma tongue pushed by the instability in the radial direction could reconnect with the overlying magnetic field, leading to the heating of the unstable flux tube. The instability can also act as a trigger to a flare or eruption if the host active region is in a subcritical equilibrium state. If the radial profiles of the equilibrium plasma parameters vary slowly along the arcade, the conditions for the corrugation instability could be fulfilled in only certain segments of the arcade, leading to the formation of a non-regular corrugated pattern. This resembles the pronounced irregular variation of the brightness of individual loops, i.e. the circularly corrugated rippling of the arcade in the longitudinal direction that is typical of post-flare arcades.
A sequence of the tongue emergences occurring at different locations along the arcade could induce a sequence of secondary flares in the host active region. In particular, this scenario could lead to the chaotic displacements of hard X-ray sources observed by Kuznetsov et al. (2016). We would like to stress that as the theory developed in this work considers perturbations with the small scales along the arcade, it is applicable to the short (in the longitudinal direction) arcades as well. Thus, the corrugation instability of plasma arcades may be responsible for several dynamical phenomena that are observed in the solar corona, which justifies the need for a better understanding of this effect.