Damped transverse oscillations of interacting coronal loops
^{1}
Departament de FísicaUniversitat de les Illes Balears,
07122
Palma de Mallorca,
Spain
email: roberto.soler@uib.es, mluna@iac.es
^{2}
Institut d’Aplicacions Computacionals de Codi Comunitari
(IAC 3), Universitat de les Illes Balears, 07122
Palma de Mallorca,
Spain
^{3}
Instituto de Astrofísica de Canarias, 38200, La
Laguna, Tenerife, Spain
^{4}
Departamento de Astrofísica, Universidad de La
Laguna, 38206,
La Laguna, Tenerife,
Spain
Received:
8
July
2015
Accepted:
4
September
2015
Damped transverse oscillations of magnetic loops are routinely observed in the solar corona. This phenomenon is interpreted as standing kink magnetohydrodynamic waves, which are damped by resonant absorption owing to plasma inhomogeneity across the magnetic field. The periods and damping times of these oscillations can be used to probe the physical conditions of the coronal medium. Some observations suggest that interaction between neighboring oscillating loops in an active region may be important and can modify the properties of the oscillations. Here we theoretically investigate resonantly damped transverse oscillations of interacting nonuniform coronal loops. We provide a semianalytic method, based on the Tmatrix theory of scattering, to compute the frequencies and damping rates of collective oscillations of an arbitrary configuration of parallel cylindrical loops. The effect of resonant damping is included in the Tmatrix scheme in the thin boundary approximation. Analytic and numerical results in the specific case of two interacting loops are given as an application.
Key words: magnetohydrodynamics (MHD) / Sun: atmosphere / Sun: corona / Sun: oscillations / waves
© ESO, 2015
1. Introduction
Transverse oscillations of magnetic loops in the solar corona have been under intense research since the first observational reports by the Transition Region And Coronal Explorer (TRACE) (see, e.g., Nakariakov et al. 1999; Aschwanden et al. 1999). Largeamplitude coronal loop oscillations are usually excited after energetic events such as solar flares, coronal mass ejections, or low coronal eruptions (see Zimovets & Nakariakov 2015). Based on magnetohydrodynamic (MHD) wave theory (e.g., Nakariakov & Verwichte 2005), transverse loop oscillations have been interpreted as standing kink MHD waves. Kink MHD modes are nearly incompressible waves, mainly driven by magnetic tension, and responsible for global transverse motions of the flux tube (see, e.g., Edwin & Roberts 1983; Goossens et al. 2009, 2012). A relevant feature of largeamplitude loop oscillations is that they are strongly damped. It has been shown that resonant absorption, caused by plasma inhomogeneity in the direction perpendicular to the magnetic field, is an efficient damping mechanism of kink MHD waves in coronal loops (see, e.g., Ruderman & Roberts 2002; Goossens et al. 2002). Because of resonant absorption, the energy from the global kink motion of the flux tube is transferred to smallscale, unresolved rotational motions around the nonuniform boundary of the tube (see, e.g., Terradas et al. 2006; Goossens et al. 2014; Soler & Terradas 2015). As a result of this process, the global kink oscillation of the coronal loop is quickly damped in time. The interested reader is referred to Goossens et al. (2011), where the theory and applications of resonant waves in the solar atmosphere are reviewed.
Observations often show that neighboring oscillating loops in an active region interact with each other and exhibit collective behavior (e.g., Schrijver & Brown 2000; Verwichte et al. 2004; Schrijver et al. 2002; White et al. 2013). Interaction between loops can modify the properties of their transverse oscillations with respect to the properties of the classic kink mode of an isolated loop. Therefore, advanced models describing coronal loop oscillations should take into account interactions within loop systems. A number of works have studied collective transverse oscillations in Cartesian geometry (e.g., Díaz et al. 2005; Díaz & Roberts 2006; Luna et al. 2006; Arregui et al. 2007, 2008). In cylindrical geometry, Luna et al. (2008) numerically investigated transverse oscillations of two cylindrical loops, and Ofman (2009) performed numerical simulations in the case of four interacting loops. Concerning analytical works in cylindrical geometry, Luna et al. (2009, 2010) used the Tmatrix theory of scattering (see, e.g., Twersky 1952; Waterman 1969; Bogdan 1987; Keppens et al. 1994) to investigate transverse oscillations of two and three loops (Luna et al. 2009) and of bundles of many loops (Luna et al. 2010) in the β = 0 approximation, where β refers to the ratio of the gas pressure to the magnetic pressure. Soler et al. (2009) later extended the method of Luna et al. (2009, 2010) by incorporating gas pressure and longitudinal flows and studied collective oscillations of flowing prominence threads. These works showed that loop interaction affects the properties of their oscillations. Luna et al. (2008, 2009) found that a system of two loops of arbitrary radii supports four kinklike collective modes. The shift of the collective mode frequencies with respect to the individual kink frequencies of the loops is significant when the distance between loops is small (approximately the size of the loop radius) and when loops have similar densities. Conversely, the oscillating loops show little interaction when they are far from each other and when their densities are substantially different. On the other hand, Van Doorsselaere et al. (2008) and Robertson et al. (2010) used a different method based on bycilindrical coordinates to study transverse oscillations of two pressureless loops in the thin tube (TT) approximation. Of the four kinklike collective modes obtained by Luna et al. (2009) in the Tmatrix theory, only two different modes remain in the TT approximation considered by Van Doorsselaere et al. (2008) and Robertson et al. (2010).
Concerning the damping of the oscillations, Arregui et al. (2007, 2008) investigated resonantly damped oscillations of two slabs, while Terradas et al. (2008) numerically studied the resonant damping of transverse oscillations of a multistranded loop. Those works showed that the process of resonant damping is not compromised by the irregular geometry of a realistic loop model and still produces the efficient attenuation of global transverse oscillations. The damping of transverse oscillations of two cylindrical loops was analytically investigated by Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014), who considered the TT approximation and used bycilindrical coordinates. Results obtained with bycilindrical coordinates should be treated with caution when the distance between loops is small. Geometrical effects intrinsically associated with the bycilindrical coordinates may produce unphysical results. Our purpose is to use the Tmatrix method of Luna et al. (2009, 2010) to investigate resonantly damped oscillations of bundles of loops. The present paper is partially based on unpublished results included in Soler (2010)^{1}. The effect of resonant absortion in the Alfvén continuum is incorporated to the Tmatrix scheme by using the method that combines the jump conditions of the perturbations at the resonance position with the socalled thin boundary (TB) approximation (see, e.g., Sakurai et al. 1991; Goossens et al. 1992). A similar method was previously used by Keppens (1995) to investigate absorption of acoustic waves. We provide a general analytic theory, which is valid for bundles of many transversely nonuniform parallel loops of arbitrary radii. Specific results in the case of two loops are obtained and compared to those given in Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014). Our results are also compared to those of Arregui et al. (2007, 2008) obtained in Cartesian geometry.
This paper is organized as follows. Section 2 contains the description of the equilibrium configuration and the basic governing equations. The general analytic Tmatrix theory of scattering to compute the frequencies and damping rates of collective loop oscillations is given in Sect. 3. Later, the specific case of damped oscillations of two loops is discussed both analytically and numerically in Sect. 4. Finally, some concluding remarks are given in Sect. 5.
2. Model and governing equations
Our equilibrium configuration is composed of N straight and parallel magnetic cylinders of length L embedded in a uniform coronal plasma. The ends of the magnetic tubes are fixed at two rigid walls representing the solar photosphere. We set the zdirection to be along the axes of the tubes. The magnetic field is straight along the zdirection, namely B = Bê_{z}, where B is a constant everywhere. We use the subscripts “i” and “e” to refer to, in general, the internal region of the tubes and the external plasma, respectively. The subscript or superscript “j” is used to refer to a particular loop. We denote by R_{j} the radius of the jth tube. The distance between the centers of the jth and j’th loops is d_{jj′}. We denote by ρ_{j} the internal density of the jth tube, while ρ_{e} denotes the external density, i.e., the density of the coronal environment. In our model ρ_{j} and ρ_{e} are constants. There is a transversely nonuniform transitional layer surrounding each magnetic tube in which the density continuously varies from the internal density, ρ_{j}, to the external density, ρ_{e}. The thicknesses of the nonuniform boundary layer of the jth cylinter is l_{j}. A sketch of the equilibrium configuration in the case of two magnetic tubes (N = 2) is given in Fig. 1.
Fig. 1 Sketch of the equilibrium configuration in the specific case of two transversely nonuniform coronal loops. 

Open with DEXTER 
We adopt the β = 0 approximation, where β refers to the ratio of the thermal pressure to the magnetic pressure. This is an appropriate approximation to investigate transverse waves in the solar corona. In the β = 0 approximation, the ideal MHD equations governing linear perturbations superimposed on the static equilibrium state are where ξ is the plasma Lagrangian displacement, b is the magnetic field Eulerian perturbation, ρ is the density, and μ is the magnetic permittivity.
We assume the temporal dependence of perturbations as exp(−iωt), where ω is the oscillation frequency. In the case of transversely nonuniform tubes, the global transverse oscillations are quasimodes whose frequency is complex owing to damping by resonant absorption, i.e., ω = ω_{R} + iω_{I}, where ω_{R} and ω_{I} are the real and imaginary parts of the frequency, respectively. The real part of ω is related to the period and the imaginary part corresponds to the damping rate of the oscillations. We consider that the oscillating flux tubes are linetied at the photosphere, which acts as a perfectly reflecting wall in this model because its density is much higher than the coronal density. Hence, we assume perturbations to be proportional to exp(ik_{z}z), with k_{z} the longitudinal wavenumber. For standing oscillations, k_{z} given by (3)We shall restict ourselves to the fundamental mode of oscillation, so we take n = 1.
In the regions with constant density, Equations (1) and (2) can be reduced to following equation, (4)where P′ = B·b/μ is the total pressure Eulerian pertubation and the subscript ⊥ refers to the direction perpendicular to the magnetic field. Thus, denotes the perpendicular part of the ∇^{2} operator. In turn, the quantity k_{⊥} plays the role of the perpendicular wavenumber and is defined as (5)where is the square of the Alfvén frequency and is the square of the Alfvén velocity. We stress that Eq. (4) is only valid in the regions with constant density, so it does not apply within the nonuniform boundaries of the loops.
3. Tmatrix theory of scattering
Equation (4) is the twodimensional Helmholtz equation. To solve Eq. (4) we use the scattering theory in its Tmatrix formalism (see, e.g., Waterman 1969). In the solar context, the Tmatrix theory has previously been used to investigate the scattering and absorption properties of bundles of magnetic flux tubes (e.g., Bogdan & Zweibel 1985; Bogdan & Cattaneo 1989; Keppens et al. 1994; Keppens 1995, among others). Luna et al. (2009, 2010) used of this technique to compute the eigenmodes of systems of magnetic tubes. Because of the inhomogeneity of the tubes in the transverse direction, the modes with frequencies between the internal Alfvén frequencies of the loops and the external Alfvén frequency are resonant in the Alfvén continuum. As a result, the oscillations are damped by resonant absorption. The effect of resonant absorption was not considered by Luna et al. (2009, 2010). Here we extend their theory to consider resonant damping.
3.1. Solutions in the internal and external plasmas
We use local polar coordinates associated with the jth loop. We denote by r_{j} and ϕ_{j} the radial and azimuthal coordinates, respectively, of the coordinate system whose origin is located at the center of the jth tube. We can define an equivalent coordinate system in each tube. In this local coordinate system, the solution to Eq. (4) in the internal region of the jth tube can be expressed as (6)where m is the azimuthal wavenumber, J_{m} is the usual Bessel function of the first kind of order m, and are constants. Unlike the case of isolated tubes (see, e.g., Edwin & Roberts 1983), the solution is not entirely described by a single value of m. Because of interaction between tubes, the values of m are coupled. For transverse, kinklike oscillations the dominant terms in the expansion are those with m = ± 1, but the contribution from other m’s is not negligible unless the tubes are far from each other (Luna et al. 2009).
The solution to Eq. (4) in the external region is written using the principle of superposition, which is applicable to linear waves. The total external solution is computed by adding the net contributions of all flux tubes, namely(7)where is the net contribution of the jth tube to the external solution. The central idea behind the scattering theory is that the solution of Eq. (4) in the external plasma can be decomposed into several fields with different physical meanings, namely the total, exciting, and scattered fields. Here we give an overview of the method. Interested readers are referred to Luna et al. (2009, 2010) for extensive explanations.
In the external plasma, the total field associated with the jth cylinder can be expressed as(8)where and are the usual Hankel functions of the first and second kind, respectively, and and constants. The first term on the righthand side of Eq. (8) represents outgoing waves from the jth tube, and the second term represents incoming waves toward the jth tube. Importantly, we note that(9)The reason for this inequality is that the outgoing wave associated with a particular tube contributes as an incoming wave for all the other tubes. In other words, is not the net external solution. To overcome this problem, the total field associated with the jth cylinder, , is decomposed into a scattered field, , and a exciting field, . Conceptually, the scattered field represents the actual contribution of the various tubes to the net external solution, whereas the exciting field can be understood as the crosstalk mechanism responsible for interaction between flux tubes (see details in, e.g., Bogdan & Cattaneo 1989).
The full net external solution, , is defined so that it corresponds to the sum of the scattered fields associated with all tubes, namely (10)Conversely, the exciting field associated with the jth tube is defined as the difference between the full net contribution and the scattered field of the jth tube, namely (11)Waterman (1969) introduced the Tmatrix operator, T^{j}, which linearly relates the scattered and exciting fields as (12)Bogdan (1987) showed that for cylindrical scatterers the Tmatrix is diagonal, and Keppens et al. (1994) gave an expression of its elements, namely (13)where and are the same constants the appear in Eq. (8). We can use Eq. (13) to eliminate and write all the expressions in terms of alone. With the help of these last formulae, and after some algebraic manipulations using wellknown properties of the Bessel functions, we can rewrite Eq. (8) as(14)from where it is straightforward to identify both exciting and scattered fields, namely Finally, we use the expression of into Eq. (10) to arrive at the total net solution in the external plasma, namely (17)Equations (6) and (17) formally describe the total pressure perturbation in the interior and in the exterior of the tubes, respectively. However, we note that these expressions do not apply in the nonuniform boundary layers. The Tmatrix elements, , contain the information about how the solutions are connected across the nonuniform boundaries of the tubes.
3.2. Tmatrix elements in the thin boundary approximation
At this stage we incorporate the effect of the nonuniform boundary layers. In the nonuniform layer of the jth tube the global wave modes are resonant in the Alfvén continuum at the resonant position, r_{j} = r_{A,j}, where the global oscillation frequency matches the local Alfvén frequency. The resonant position, r_{A,j}, is defined through the resonant condition . We use the TB approximation and restrict ourselves to l_{j}/R_{j} ≪ 1. The TB approximation assumes that the jump of the perturbations across the resonant layer is the same as their jump across the whole nonuniform layer. Thus, the connection formulae of the wave perturbations across the resonance are used as jump conditions for the total pressure and the Lagrangian displacement at the boundaries of the tubes. This method and its applications have been reviewed by Goossens et al. (2011). The TB approximation within the formalism of the Tmatrix theory has previously been used by Keppens et al. (1994) and Keppens (1995).
General expressions of the connection formulae for the perturbations across the resonant layer can be found in, e.g., Sakurai et al. (1991). In local coordinates, the connection formulae for the total pressure Eulerian pertubation, P′, and the radial component of the Lagrangian displacement, ξ_{r}, at r_{j} = r_{A,j} are (18)where [X] = X_{e}−X_{j} denotes the jump of the quantity X across the resonant layer and ρΔ_{A}_{j} is defined as (19)where we have used the resonant condition . After imposing the jump conditions given in Eq. (18), we obtain both the Tmatrix elements and the dispersion relation of the collective modes.
On the one hand, the Tmatrix elements are (20)where the prime ′ denotes the derivative of the Bessel or Hankel function with respect to its argument. In the absence of resonant damping, Eq. (20) consistently reverts to Eq. (17) of Luna et al. (2009).
On the other hand, the constants satisfy an algebraic system of equations, namely (21)for −∞ <m< ∞, where ϕ_{jj′} is the angle formed by the vector positions of the centers of the two tubes with respect to the global reference frame. Equation (21) is a system with an infinite number of algebraic equations. The condition that there is a nontrivial solution, i.e., the determinant formed by the coefficients set equal to zero, provides us with the dispersion relation (see details in Luna et al. 2009, 2010). The solution of the dispersion relation is the complex frequency of the damped collective quasimode. The imaginary part of the frequency is the resonant damping rate. In the case of Luna et al. (2009, 2010), the solution frequency was real because of the absence of resonant damping.
The dispersion relation is a complicated expression and has to be solved by numerical methods. For practical computational purposes, in Eq. (21) the indices m and m′ must be truncated into a finite number. The truncation term must be large enough to avoid inaccuracy of the solutions. Then, the dispersion relation can be solved using standard numerical routines to find the roots of transcendental equations.
4. Application to damped oscillations of two loops
We note that the Tmatrix theory described in the previous Section is valid for an arbitrary number of loops. Here we apply that method to investigate damped transverse oscillations in a loop system composed of two magnetic tubes alone. First we derive approximate expressions of the period and damping time in the case of two identical thin tubes. Later, we consider two tubes with different properties and arbitrary radii and perform a numerical study.
4.1. Approximate solutions for two identical thin tubes
We consider the paradigmatic case of two tubes with identical properties. We denote by ρ_{i}, R, and l the internal density, the radius, and the thickness of the boundary layer of both tubes, respectively. The external density remains ρ_{e}. The distance between the centers of the two tubes is d.
To obtain an approximate dispersion relation for kinklike modes we assume that the contribution from the azimuthal wavenumbers m = ± 1 to the collective oscillation is much more important than the contributions from other values of m. Thus, we take m,m′ = ± 1 in Eq. (21) and roughly neglect the other terms. As pointed out by Luna et al. (2009), the coupling between the different values of m gets stronger as the separation between the cylinders decreases. In terms of the parameters of the model, this means that our approximation applies to the case of large separations, i.e., R/d ≪ 1.
The system defined in Eq. (21) now becomes an algebraic system of four equations for the unknowns , , , and . We obtain the dispersion relation from the condition that there is a nontrivial solution of the system, namely (22)where we used the property that T_{11} = T_{− 1−1} according to the symmetry relations of the Bessel and Hankel functions (see Abramowitz & Stegun 1972). Equation (22) is an equation for ω, which is enclosed in the definitions of T_{11} and k_{⊥ e}. To go further analytically, we perform a series expansion of the Hankel functions for small arguments and retain the first term alone. This would be approximately valid for thin tubes. After some algebraic manipulations, Eq. (22) can be recast as (23)Now, we look for a simplified expression of T_{11}. We use the TT approximation, i.e., R/L ≪ 1. In Eq. (20) we perform an asymptotic expansion of the Bessel and Hankel functions for small arguments. Equation (20) becomes (24)where we also assumed r_{A} ≈ R for simplicity. We take m = 1, substitute Eq. (24) into Eq. (23), and arrive at the following expression, (25)where the parameter δ is defined as (26)Equation (25) is the dispersion relation for resonantly damped kink modes of two identical tubes in the TT and TB approximations and for large separations. In the limit that the tubes are far from each other, R/d → 0 and so δ → 1. Then, Eq. (25) consistently reverts to the dispersion relation of resonant kink modes in an isolated thin tube (e.g., Goossens et al. 1992).
4.1.1. Frequency of the oscillations
First we neglect the presence of the nonuniform boundary layers and study undamped oscillations. We drop the last term on the lefthand side of Eq. (25) and the dispersion relation simplifies to(27)The exact solution is (28)where ζ = ρ_{i}/ρ_{e} is the density contrast and ω_{k} is the frequency of the kink mode in an individual thin tube (Edwin & Roberts 1983), namely (29)Equation (28) can be compared to Eq. (51) of Van Doorsselaere et al. (2008). Both expressions agree if the parameter E of Van Doorsselaere et al. (2008) is approximated as , which is valid for R/d ≪ 1. In agreement with Van Doorsselaere et al. (2008), in the TT approximation we obtain two kinklike solutions corresponding to the + and − signs in Eq. (28). As explained by Luna et al. (2009), there are actually four kinklike modes beyond the TT approximation. For arbitrary radii, the lowfrequency solution (that with the + sign in Eq. (28)) splits in the S_{x} and A_{y} modes of Luna et al. (2009), while the highfrequency solution (that with the − sign in Eq. (28)) becomes their S_{y} and A_{x} solutions.
From Eq. (28) we compute the period of the oscillations, P = 2π/ω, as (30)with P_{k} = 2π/ω_{k} the period of the kink mode of an isolated tube. For large separations the tubes feel little interaction. Then, the periods of the two solutions consistently tend to that of the kink mode of an isolated tube.
4.1.2. Damping rate
Here we incorporate the effect of resonant damping and consider the full expression of the dispersion relation (Eq. (25)). To obtain an approximate expression of the damping rate, we substitute ω = ω_{R} + iω_{I} in Eq. (25) and assume weak damping, i.e.,  ω_{I}  ≪  ω_{R} . Then, we are allowed to neglect terms of order and higher orders. After some algebraic maniputations (not given here for the sake of simplicity) we arrive at an expression for the ratio ω_{I}/ω_{R}, namely (31)To simplify Eq. (31) we consider a smooth monotonic profile for the density in the nonuniform boundaries so that we can express the derivative of the density profile as (32)with ℱ a factor that depends on the form of the profile. For instance, ℱ = 4 /π^{2} for a linear profile and ℱ = 2 /π for a sinusoidal profile. In addition, we approximate ω_{R} by the value of the frequency in the undamped case (Eq. (28)) and substitute the expression of δ (Eq. (26)) to arrive at (33)For a linear density profile, Eq. (33) becomes Eq. (86) of Robertson & Ruderman (2011) if the approximation valid for R/d ≪ 1 is again performed in their expression.
Now, we define the damping time as τ_{D} = 1/  ω_{I}  and use Eq. (33) to obtain the expression for the ratio τ_{D}/P. By keeping terms up to (R/d)^{2}, the expression for the damping ratio is(34)where (35)is the damping ratio of the individual kink mode (see, e.g., Ruderman & Roberts 2002; Goossens et al. 2002). As for individual kink modes, the damping ratio of collective modes is inversely proportional to the thickness of the nonuniform layer. Equation (34) shows that the solution with the + sign is more efficiently damped by resonant absorption than the solution with the − sign. This difference in the damping rates of the two modes qualitatively agrees with the results of Robertson & Ruderman (2011).
4.2. Numerical study
Here we perform a numerical study beyond the approximate limits studied before. To do so, we consider the full dispersion relation obtained from Eq. (21) by using the general expression of the Tmatrix elements (Eq. (20)). The truncated dispersion relation is then numerically solved. Convergence tests of the results have been performed to make sure that the truncation value of the azimuthal series is large enough for the error on the solutions to be negligible. In short, we found that for kinklike modes the truncation value, namely m_{t}, has no important effect unless m_{t}< 5 and the two loops are next to each other (d/R ≈ 2). When m_{t}> 5 and/or the loops are far from each other, we essentially obtain that the results are independent of m_{t}. We have used m_{t} = 30 in all computations given here, which reduces the error and ensures the excellent converge of the solutions even when d/R ≈ 2.
The numerical method follows a twostep procedure. First, we solve the dispersion relation in the absence of nonuniform boundary layers. In that case, the solution is a real frequency that corresponds to an approximation to ω_{R}. This approximate value of ω_{R} is used to compute the resonance positions and the derivative of the density profile at the resonances. We assume sinusoidal density profiles in the nonuniform boundary layers. Then, we use these parameters to solve the complete dispersion relation, which now includes the effect of resonant damping. The frequency obtained from the second run is complex, so that it provides us with a more accurate value of ω_{R} and also gives us the value of ω_{I}.
4.2.1. Identical loops
We initially study the case of two identical tubes. We set the Cartesian coordinates system so that the xyplane is perpendicular to the axes of the loops. In that plane, the centers of the two loops are located on the xaxis. We use the notation introduced by Luna et al. (2008) to denote the four kinklike modes present in a twoloop configuration. The modes are labeled as S_{x}, A_{x}, S_{y}, and A_{y}, where S and A denote symmetric or antisymmetric motions of the two loops, respectively, and the subscripts x and y indicate the main direction of polarization of the oscillations in the coordinates system defined above. The eigenfunctions of these four modes in the case of loops without nonuniform boundary layers can be found in Fig. 2 of Luna et al. (2008).
Fig. 2 Numerical results in the case of two identical coronal loops. a) Dependence of ω_{R}/ω_{A,e} on d/R, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on d/R. The meaning of the various lines is indicated within the panels. The dashed lines correspond to the analytic approximations in the limit d/R ≫ 1 (Eqs. (28) and (33)). We have used ζ = 5, L/R = 100, and l/R = 0.2. Panels c) and d) show the same results as panels a) and b), respectively, but as functions of log _{10}^{(}d/R−2^{)}. 

Open with DEXTER 
Fig. 3 Numerical results in the case of two identical coronal loops. a) Dependence of ω_{R}/ω_{A,e} on L/R, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on L/R. The meaning of the various lines is indicated within the figure. We have used ζ = 5, l/R = 0.2, and d/R = 2.5. We note that the horizontal axes of both panels are in logarithmic scale. 

Open with DEXTER 
Fig. 4 Numerical results in the case of two identical coronal loops. a) Dependence of ω_{R}/ω_{A,e} on l/R, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on l/R. The meaning of the various lines is indicated within the figure. We have used ζ = 5, L/R = 100, and d/R = 2.5. 

Open with DEXTER 
Fig. 5 Numerical results in the case of two nonidentical coronal loops with R_{1} = R_{2} = R. a) Dependence of ω_{R}/ω_{A,e} on l_{2}/R, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on l_{2}/R. The meaning of the various lines is indicated within the figure. We have used ζ_{1} = ζ_{2} = 5, L/R = 100, l_{1}/R = 0.2 and d/R = 2.5. 

Open with DEXTER 
Figure 2 shows the dependence of ω_{R} and the ratio ω_{I} /ω_{R} on the separation between loops, d/R, for a particular set of parameters given in the caption of the figure. We note that the curves corresponding to the S_{x} and A_{y} modes, and those of the A_{x} and S_{y} modes, are almost superimposed because we are in the TT regime (we used L/R = 100). Concerning the behavior of ω_{R}, the frequencies of the four solutions tend to the kink frequency of an isolated loop in the limit d/R ≫ 1. Conversely, the smaller the separation between loops, the more important the splitting of the collective frequencies with respect to the kink frequency of an isolated loop. Figure 2a can be compared to Fig. 3 of Luna et al. (2008) and with Fig. 4 of Van Doorsselaere et al. (2008). However, unlike in those previous works we note that in our case the frequencies of the highfrequency solutions (A_{x} and S_{y}) do not tend to the external Alfvén frequency when d/R → 2. The reason for this difference is probably that the A_{x} and S_{y} modes are strongly damped when d/R → 2, and this has some impact on the real part of the frequency as well.
On the other hand, Fig. 2b shows that the highfrequency modes (A_{x} and S_{y}) are more efficiently damped by resonant absorption than the lowfrequency modes (S_{x} and A_{y}). This result agrees with that of Robertson & Ruderman (2011) for large separations. However, the behavior of ω_{I} /ω_{R} obtained here for small separations is dramatically different from that of Robertson & Ruderman (2011). They found that the oscillations become undamped in the limit d/R → 2, while we find that the modes remain damped. Although the analysis of Robertson & Ruderman (2011) is mathematically correct, we find no physical reason for which the oscillations should be undamped when the loops are close to each other. Our computations show that the damping of the lowfrequency modes (S_{x} and A_{y}) is roughly independent on d/R, whereas the damping of the highfrequency modes (A_{x} and S_{y}) gets stronger as the separation between loops is reduced. As pointed out by Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014), the physical significance of the results obtained with bycilindrical coordinates should be treated with caution when the separation between the tubes is small.
Panels c and d of Fig. 2 show the same results as panels a and b, respectively, but now log _{10}^{(}l/R−2^{)} is used in the horizontal axes. These additional graphs are included to show in more detail the behavior of the solutions obtained with the Tmatrix method for small separations between tubes.
We have overplotted in Fig. 2 the analytic approximations of ω_{R} and ω_{I} /ω_{R} given in Eqs. (28) and (33). These approximations were derived in the limit d/R ≫ 1 and reasonably agree with the numerical solutions when d/R ≳ 3. As expected, the approximations do not work well for small separations. The analytic approximations were derived considering the contributions from m = ± 1 alone, but the contribution of high m’s to the full solution is important for small separation between loops (see Luna et al. 2009).
Figure 3 shows the dependence of the solutions on L/R. This figure is included to show that the almost degenerate couples S_{x}–A_{y} and A_{x}–S_{y} split into four different solutions for small values of L/R beyond the TT regime. We point out, however, that the impact of the value of L/R on the solutions is not relevant when realistic values of this parameter are considered. The TT limit used by Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014) is therefore adequate.
Now we plot in Fig. 4 the dependence of ω_{R} and ω_{I} /ω_{R} on the nonuniform layer thickness, l/R. We consider a small separation, namely d/R = 2.5, and the remaining parameters are the same as in Fig. 2. Consistently, when l/R = 0 the modes are undamped. The real part of the frequency of the lowfrequency modes is almost independent of l/R, while their ω_{I} /ω_{R} is roughly linear with l/R. Conversely, the real part of the frequency of the lowfrequency modes decreases when l/R increases, and their ω_{I} /ω_{R} is only linear with l/R for small values of this parameter. As discussed before, the behavior of the lowfrequency modes (S_{x} and A_{y}) is similar to that of the kink mode of an isolated loop. However, the highfrequency modes (A_{x} and S_{y}) seem to be more affected by the interaction between loops and show a somewhat different behavior when l/R increases. We note that because of the TB approximation we are restricted to consider small values of l/R.
It is useful to relate the present results with those of Arregui et al. (2007, 2008), who studied the damping of transverse oscillations of two nonuniform slabs. They found that the ratio ω_{I} /ω_{R} corresponding to the symmetric kink mode of the two slabs is weakly dependent of the separation between the two slabs (see Arregui et al. 2007, their Fig. 4). In turn, Arregui et al. (2008) found that the antisymmetric kink mode of the two slabs is more efficiently damped than the symmetric kink mode (see their Fig. 6). The symmetric and antisymmetric kink modes of two slabs would be equivalent to the S_{x} and A_{x} modes of two cylinders. Thus, our results in cylindrical geometry are consistent with previous findings in Cartesian geometry.
It is also convenient to consider the physical arguments of Arregui et al. (2007, 2008) to explain why the highfrequency modes damp more efficiently than the lowfrequency modes. Arregui et al. (2007, 2008) related the efficiency of the damping with the magnitude of the total pressure perturbation within the resonant layers. According to Andries et al. (2000), the efficiency of the resonant coupling between the global transverse mode and the Alfvén continuum modes is proportional to the total pressure perturbation squared. Soler (2010) plotted the square of the total pressure perturbation corresponding to the S_{x} and A_{x} modes (see his Figs. 9.7 and 9.8) and found that, when the quantities are normalized, the perturbation of the A_{x} mode reaches a larger value in the resonant layers than that of the S_{x} solution. This result qualitatively explains why resonant damping is more efficient for the A_{x} mode than for the S_{x} mode. Equivalently, a similar reasoning help us understand the different attenuation of the S_{y} and A_{y} modes. Nevertheless, a more robust study of the process of resonant absorption in twodimensional configurations would be needed for a complete understanding of the different damping rates (see Russell & Wright 2010).
Fig. 6 Numerical results in the case of two nonidentical coronal loops with R_{1} = R_{2} = R. a) Dependence of ω_{R}/ω_{A,e} on ζ_{1}, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on ζ_{1}. The dotted lines correspond to the Alfvén frequencies of the two loops, and the meaning of the remaining lines is indicated within the figure. We have used ζ_{2} = 3, L/R = 100, l_{1}/R = l_{2}/R = 0.2 and d/R = 3. 

Open with DEXTER 
4.2.2. Nonidentical loops
Here we consider two loops with different properties and compare our results to those of Gijsen & Van Doorsselaere (2014). We use subscripts 1 and 2 to refer to the two different loops. For simplicity, we take R_{1} = R_{2} = R and L/R = 100 in all following computations.
First, we assume the same density contrast in the two loops, namely ζ_{1} = ζ_{2} = 5, and vary l_{2}/R while l_{1}/R is kept fixed to l_{1}/R = 0.2. These results are shown in Fig. 5 and can be compared to those already presented in Fig. 4 in the case of l_{1} = l_{2}. Importantly, we find that the collective oscillations remain damped when l_{2}/R = 0. Although this is a rather particular situation, the results have interesting implications. When l_{2}/R = 0 resonant absorption only occurs in the boundary layer of loop #1. However, this is enough for the collective oscillations of the two loops to be efficiently damped. We note that in the study by Gijsen & Van Doorsselaere (2014) the thicknesses of the nonuniform layers are linked to the coordinate system.
Now we consider the case of two loops with different density contrasts. We fix ζ_{2} = 3 and compute the solutions as functions of ζ_{1}. These results are given in Fig. 6, where the remaining parameters used in the computations are specified in the caption. For consistency, we keep the same notation as before to denote the various modes according to the ordering of their frequencies, although they do not represent truly collective oscillations if ζ_{1} ≠ ζ_{2} (Luna et al. 2009). Figure 6 can be compared to Fig. 6 of Gijsen & Van Doorsselaere (2014). To make a proper comparison, we note that Gijsen & Van Doorsselaere (2014) plotted the signed damping ratio, while here we plot the absolute value.
For the real part of the frequency (Fig. 6a), we find the same results as Luna et al. (2009). When ζ_{1}<ζ_{2}, the highfrequency modes are associated with loop #1 alone, whereas the lowfrequency modes represent individual oscillations of loop #2. The opposite happens when ζ_{1}>ζ_{2}. Conversely, when ζ_{1} ≈ ζ_{2} the four modes approach and interact in the form of what is known as an avoided crossing. Only in that case do the modes represent truly collective oscillations (Luna et al. 2009).
Figure 6a also shows that the highfrequency modes are always within the Alfvén continua of the two loops, i.e., the frequencies of the A_{x} and S_{y} modes are always larger than ω_{A,1} and ω_{A,2} and smaller than ω_{A,e}. This is true except in the limit ζ_{1} → 1, where the frequencies of the A_{x} and S_{y} modes tend to ω_{A,e}. On the contrary, the lowfrequency S_{x} and A_{y} modes are below the Alfvén continuum of loop #1 when ζ_{1} ≲ 2 and below the Alfvén continuum of loop #2 when ζ_{1} ≳ 5. This result may have implications for the damping by resonant absorption.
Figure 6b shows the damping ratio of the modes as a function of ζ_{1}. The result for the highfrequency modes can be understood as follows. When ζ_{1} → 1 the frequencies of the highfrequency modes tend to ω_{A,e} and, as a consequence, these modes become undamped in that limit. When 1 <ζ_{1}<ζ_{2}, the damping ratio increases when ζ_{1} increases. When 1 <ζ_{1}<ζ_{2} the highfrequency modes represent individual oscillations of loop #1. Then, the damping ratio reaches a maximum when ζ_{1} ≈ ζ_{2}. When ζ_{1}>ζ_{2} the damping ratio saturates to a constant value because the highfrequency modes represent now individual oscillations of loop #2, and the value of ζ_{2} is fixed in the computations. The behavior of the damping of the highfrequency modes agrees with that plotted by Gijsen & Van Doorsselaere (2014) in their Fig. 6.
The overall behavior of the damping ratio of the lowfrequency modes presented in Fig. 6b also agrees with Gijsen & Van Doorsselaere (2014). The damping ratio of the lowfrequency modes is roughly constant when ζ_{1}<ζ_{2} and increases when ζ_{1}>ζ_{2}. The lowfrequency modes are below the Alfvén continuum of loop #1 when ζ_{1} ≲ 2 and below the Alfvén continuum of loop #2 when ζ_{1} ≳ 5, but this has no important impact on the damping. Again, these results can be understood by considering that the lowfrequency modes are associated with loop #2 when ζ_{1}<ζ_{2}, while they are associated with loop #1 when ζ_{1}>ζ_{2}.
The lowfrequency modes computed here do not show the pronounced minimum of the damping rate seen in the solution plotted by Gijsen & Van Doorsselaere (2014) when ζ_{1} ≈ ζ_{2}. There are several effects that may explain this difference. The most obvious one is the different geometry considered in Gijsen & Van Doorsselaere (2014) and here. Another possible explanation is that the density profile in the nonuniform layers used by Gijsen & Van Doorsselaere (2014) is different from that used here. A linear density profile is used in Gijsen & Van Doorsselaere (2014), so that the derivative of density at the resonance position is independent of the frequency of the mode. Here we use a sinusoidal profile and take into account that the position of the resonance and the value of the derivative of density at the resonance position are functions of the mode frequency.
In Fig. 6b, the damping rate of the lowfrequency modes shows a small bump around ζ_{1} ≈ 2, when the real part of the frequency approximately crosses the internal Alfvén frequency of loop #1. The reason for this bump is that the lowfrequency modes intersect with and avoidcross the fluting modes that cluster toward the internal Alfvén frequency. This bump is absent from Fig. 6 of Gijsen & Van Doorsselaere (2014) probably because coupling between kink and fluting modes is not described in the TT approximation.
5. Concluding remarks
In this paper we have extended the analytic Tmatrix theory of scattering of Luna et al. (2009, 2010) to investigate resonantly damped oscillations of an arbitrary configuration of parallel cylindrical coronal loops. After presenting the general theory, we have performed a specific application in the case of two loops. This work is partially based on unpublished results included in Soler (2010), where collective damped oscillations of prominence threads were studied.
We have compared our results to those of the papers by Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014). They investigated the damping of collective oscillations of two loops in the TT approximation and used a method based on bicylindrical coordinates. In general, the results of Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014) are in good agreement with the present results, specially when the separation between loops is large. However, when the separation between the loops is small, i.e., for separations of a few radii, the results of those previous works show important differences with the present findings. For instance, Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014) obtained that by decreasing the distance between loops, the efficiency of resonant damping is reduced. In their computations, both low and highfrequency modes become undamped when the loops are in contact. However, this result lacks of a physical explanation and contradicts previous findings in Cartesian geometry (Arregui et al. 2007, 2008). In our computations, we find that the damping of the highfrequency modes gets stronger by decreasing the separation between loops, while the damping of the lowfrequency modes is roughly independent of the separation. Our solutions do not become undamped when the two tubes are in contact. Thus, the results obtained here are consistent with previous results by Arregui et al. (2007, 2008) of collective oscillations of two slabs.
Although the mathematical analysis of Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014) is flawless, their results may be affected by unavoidable geometrical problems related to the bicylindrical coordinates when the loops are close to each other. In bicylindrical coordinates the shapes of nonuniform boundary layers are not symmetric and change when the separation between tubes decreases. The nonuniform layers get thicker in the outer parts of the tubes and thinner in the inner parts. As already mentioned by Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014), these geometrical limitations may lead to unphysical results for small separations. The Tmatrix method used here is not constrained by the geometrical problems of the bicylindrical coordinates. Therefore, we may conclude that the results given here are more generally applicable than those of Robertson & Ruderman (2011) and Gijsen & Van Doorsselaere (2014) when the loops are close to each other.
Because of the TB approximation we were restricted to consider small values of l/R, i.e., l/R ≪ 1. The effect of thick nonuniform layers could be included into the Tmatrix formalism with the method of Frobenius used by Soler et al. (2013) in the case of an isolated loop. This would substantially increase the mathematical complexity of the problem but, on the other hand, it would provide a more accurate description of the damping of largely nonuniform loops. It has been shown by Soler et al. (2014) that the error on the damping rate associated with the use of the TB approximation can be important when the loops are largely nonuniform. Apart from a numerical factor, the damping rate in the TB approximation is independent of the specific density profile considered within the nonuniform boundary, but the density profile can have a more important impact when the nonuniform layers are thick. In addition, the real part of the frequency depends on l/R beyond the limit l/R ≪ 1. The effect of thick nonuniform boundaries on collective loop oscillations could be explored in the future.
The method given here to compute resonantly damped collective oscillations can have multiple applications in the future. For instance, damped oscillations of a coronal arcade could be studied by modeling the arcade as a long line of parallel loops. Another interesting application is the investigation of oscillations of loops formed by many strands. As shown by Luna et al. (2010), the global oscillation of the whole loop would be determined by the interaction of the oscillations of the individual strands. The resonant absorption process working in the individual strands would affect the damping on the global loop motion (see Terradas et al. 2008). In principle, the presence of multiple resonances in the system may cause the transverse oscillations of a multistranded loop to damp more quickly than the oscillations of an equivalent monolithic loop. This idea could be confirmed using the Tmatrix method. In addition, as shown by Soler et al. (2009), the effects of gas pressure and mass flow along the loops can easily be included in the Tmatrix formalism, thus extending the applicability of the method.
The full text of Soler (2010) is available at http://www.uib.es/depart/dfs/Solar/thesis_robertosoler.pdf
Acknowledgments
We thank Ramón Oliver for useful comments on a draft of this paper. R.S. acknowledges support from MINECO through a “Juan de la Cierva” grant and through projects AYA201122846 and AYA201454485P, from MECD through project CEF110012, from the “Vicerectorat d’Investigació i Postgrau” of the UIB, and from FEDER funds. M.L. acknowledges support by MINECO through projects AYA201124808 and AYA201455078P. M.L. is also grateful to ERC2011StG 277829SPIA.
References
 Abramowitz, M., & Stegun, I. A. 1972, Handbook of Mathematical Functions (National Bureau of Standards – U.S. Government Printing Office) [Google Scholar]
 Andries, J., Tirry, W. J., & Goossens, M. 2000, ApJ, 531, 561 [NASA ADS] [CrossRef] [Google Scholar]
 Arregui, I., Terradas, J., Oliver, R., & Ballester, J. L. 2007, A&A, 466, 1145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Arregui, I., Terradas, J., Oliver, R., & Luis Ballester, J. 2008, ApJ, 674, 1179 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J., Fletcher, L., Schrijver, C. J., & Alexander, D. 1999, ApJ, 520, 880 [NASA ADS] [CrossRef] [Google Scholar]
 Bogdan, T. J. 1987, ApJ, 318, 888 [NASA ADS] [CrossRef] [Google Scholar]
 Bogdan, T. J., & Cattaneo, F. 1989, ApJ, 342, 545 [NASA ADS] [CrossRef] [Google Scholar]
 Bogdan, T. J., & Zweibel, E. G. 1985, ApJ, 298, 867 [NASA ADS] [CrossRef] [Google Scholar]
 Díaz, A. J., & Roberts, B. 2006, Sol. Phys., 236, 111 [NASA ADS] [CrossRef] [Google Scholar]
 Díaz, A. J., Oliver, R., & Ballester, J. L. 2005, A&A, 440, 1167 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Gijsen, S. E., & Van Doorsselaere, T. 2014, A&A, 562, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goossens, M., Hollweg, J. V., & Sakurai, T. 1992, Sol. Phys., 138, 233 [NASA ADS] [CrossRef] [Google Scholar]
 Goossens, M., Andries, J., & Aschwanden, M. J. 2002, A&A, 394, L39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goossens, M., Terradas, J., Andries, J., Arregui, I., & Ballester, J. L. 2009, A&A, 503, 213 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goossens, M., Erdélyi, R., & Ruderman, M. S. 2011, Space Sci. Rev., 158, 289 [NASA ADS] [CrossRef] [Google Scholar]
 Goossens, M., Andries, J., Soler, R., et al. 2012, ApJ, 753, 111 [NASA ADS] [CrossRef] [Google Scholar]
 Goossens, M., Soler, R., Terradas, J., Van Doorsselaere, T., & Verth, G. 2014, ApJ, 788, 9 [NASA ADS] [CrossRef] [Google Scholar]
 Keppens, R. 1995, Sol. Phys., 161, 251 [NASA ADS] [CrossRef] [Google Scholar]
 Keppens, R., Bogdan, T. J., & Goossens, M. 1994, ApJ, 436, 372 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., Terradas, J., Oliver, R., & Ballester, J. L. 2006, A&A, 457, 1071 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Luna, M., Terradas, J., Oliver, R., & Ballester, J. L. 2008, ApJ, 676, 717 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., Terradas, J., Oliver, R., & Ballester, J. L. 2009, ApJ, 692, 1582 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., Terradas, J., Oliver, R., & Ballester, J. L. 2010, ApJ, 716, 1371 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., & Verwichte, E. 2005, Liv. Rev. Sol. Phys., 2, 3 [Google Scholar]
 Nakariakov, V. M., Ofman, L., Deluca, E. E., Roberts, B., & Davila, J. M. 1999, Science, 285, 862 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ofman, L. 2009, ApJ, 694, 502 [NASA ADS] [CrossRef] [Google Scholar]
 Robertson, D., & Ruderman, M. S. 2011, A&A, 525, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Robertson, D., Ruderman, M. S., & Taroyan, Y. 2010, A&A, 515, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ruderman, M. S., & Roberts, B. 2002, ApJ, 577, 475 [NASA ADS] [CrossRef] [Google Scholar]
 Russell, A. J. B., & Wright, A. N. 2010, A&A, 511, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sakurai, T., Goossens, M., & Hollweg, J. V. 1991, Sol. Phys., 133, 227 [NASA ADS] [CrossRef] [Google Scholar]
 Schrijver, C. J., & Brown, D. S. 2000, ApJ, 537, L69 [NASA ADS] [CrossRef] [Google Scholar]
 Schrijver, C. J., Aschwanden, M. J., & Title, A. M. 2002, Sol. Phys., 206, 69 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R. 2010, Ph.D. Thesis, Departament de Fisica, Universitat de les Illes Balears, http://www.uib.es/depart/dfs/Solar/thesis_robertosoler.pdf [Google Scholar]
 Soler, R., & Terradas, J. 2015, ApJ, 803, 43 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Oliver, R., & Ballester, J. L. 2009, ApJ, 693, 1601 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Goossens, M., Terradas, J., & Oliver, R. 2013, ApJ, 777, 158 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Goossens, M., Terradas, J., & Oliver, R. 2014, ApJ, 781, 111 [NASA ADS] [CrossRef] [Google Scholar]
 Terradas, J., Oliver, R., & Ballester, J. L. 2006, ApJ, 642, 533 [NASA ADS] [CrossRef] [Google Scholar]
 Terradas, J., Arregui, I., Oliver, R., et al. 2008, ApJ, 679, 1611 [NASA ADS] [CrossRef] [Google Scholar]
 Twersky, V. 1952, Acoustical Soc. Am. J., 24, 42 [NASA ADS] [CrossRef] [Google Scholar]
 Van Doorsselaere, T., Ruderman, M. S., & Robertson, D. 2008, A&A, 485, 849 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verwichte, E., Nakariakov, V. M., Ofman, L., & Deluca, E. E. 2004, Sol. Phys., 223, 77 [NASA ADS] [CrossRef] [Google Scholar]
 Waterman, P. C. 1969, Acoustical Soc. Am. J., 45, 1417 [NASA ADS] [CrossRef] [Google Scholar]
 White, R. S., Verwichte, E., & Foullon, C. 2013, ApJ, 774, 104 [NASA ADS] [CrossRef] [Google Scholar]
 Zimovets, I. V., & Nakariakov, V. M. 2015, A&A, 577, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Figures
Fig. 1 Sketch of the equilibrium configuration in the specific case of two transversely nonuniform coronal loops. 

Open with DEXTER  
In the text 
Fig. 2 Numerical results in the case of two identical coronal loops. a) Dependence of ω_{R}/ω_{A,e} on d/R, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on d/R. The meaning of the various lines is indicated within the panels. The dashed lines correspond to the analytic approximations in the limit d/R ≫ 1 (Eqs. (28) and (33)). We have used ζ = 5, L/R = 100, and l/R = 0.2. Panels c) and d) show the same results as panels a) and b), respectively, but as functions of log _{10}^{(}d/R−2^{)}. 

Open with DEXTER  
In the text 
Fig. 3 Numerical results in the case of two identical coronal loops. a) Dependence of ω_{R}/ω_{A,e} on L/R, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on L/R. The meaning of the various lines is indicated within the figure. We have used ζ = 5, l/R = 0.2, and d/R = 2.5. We note that the horizontal axes of both panels are in logarithmic scale. 

Open with DEXTER  
In the text 
Fig. 4 Numerical results in the case of two identical coronal loops. a) Dependence of ω_{R}/ω_{A,e} on l/R, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on l/R. The meaning of the various lines is indicated within the figure. We have used ζ = 5, L/R = 100, and d/R = 2.5. 

Open with DEXTER  
In the text 
Fig. 5 Numerical results in the case of two nonidentical coronal loops with R_{1} = R_{2} = R. a) Dependence of ω_{R}/ω_{A,e} on l_{2}/R, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on l_{2}/R. The meaning of the various lines is indicated within the figure. We have used ζ_{1} = ζ_{2} = 5, L/R = 100, l_{1}/R = 0.2 and d/R = 2.5. 

Open with DEXTER  
In the text 
Fig. 6 Numerical results in the case of two nonidentical coronal loops with R_{1} = R_{2} = R. a) Dependence of ω_{R}/ω_{A,e} on ζ_{1}, where ω_{A,e} = k_{z}v_{A,e} is the external Alfvén frequency. b) Dependence of ω_{I} /ω_{R} on ζ_{1}. The dotted lines correspond to the Alfvén frequencies of the two loops, and the meaning of the remaining lines is indicated within the figure. We have used ζ_{2} = 3, L/R = 100, l_{1}/R = l_{2}/R = 0.2 and d/R = 3. 

Open with DEXTER  
In the text 