A&A 469, 1135-1143 (2007)
DOI: 10.1051/0004-6361:20077404
J. Terradas^{} - J. Andries^{} - M. Goossens
Centrum voor Plasma-Astrofysica, Katholieke Universiteit Leuven, Celestijnenlaan 200B, 3001 Leuven, Belgium
Received 5 March 2007 / Accepted 7 April 2007
Abstract
Context. Flares and eruptions in the solar corona generate oscillations of loops which have been interpreted as eigenmodes (mainly the fundamental kink mode, although other modes can also be excited). From the theoretical point of view the excitation of the tube eigenmodes due to an initial disturbance has not been studied in much detail.
Aims. The main aim of this work is to calculate for a given initial disturbance the amount of energy that is deposited in the trapped fast mode oscillation, how it depends on the initial perturbation and how it is distributed among the different eigenmodes (kink and fluting and also the longitudinal harmonics).
Methods. We calculate, using analytical expressions, the amplitude and the energy of the oscillation of the magnetic tube for different kinds of initial excitations.
Results. We find that external excitations deposit a small amount of energy in the tube. We show that fluting modes have quite small energies in comparison with the energy of the kink mode (around three orders of magnitude for the first fluting mode). On the contrary, the longitudinal fundamental mode and the longitudinal harmonics have energies of the same order of magnitude. In addition, we find that the loop length and density contrast can be important factors that determine the amount of energy that is trapped by the loop.
Conclusions. The energy deposited in loops is typically six orders of magnitude smaller than the energy of the initial disturbance (for external excitations). However, it strongly depends on the distance of the initial perturbation and also on the loop properties (length and density). Fluting modes in coronal loops are very difficult to excite. Longitudinal harmonics are in principle more easily excited.
Key words: magnetohydrodynamics (MHD) - Sun: corona - Sun: magnetic fields - waves
Several kinds of wave motions have been observed in coronal loops using the TRACE and SOHO satellites and other ground-based instruments. The observed loop oscillations (see for example Aschwanden et al. 2002; Schrijver et al. 2002) have been compared with the information of the theoretical models of coronal loop oscillations (see Roberts 2004,2002). A possible interpretation of these oscillations is that they are the natural motions of the tube. With the use of the interpretation in terms of magnetohydrodynamic (MHD) eigenmodes, coronal seismology, suggested by e.g. Uchida (1970) and Roberts et al. (1984), has been performed by Nakariakov & Ofman (2001) to determine the magnetic field and by Goossens et al. (2002); Aschwanden et al. (2003) and Arregui et al. (2007) to determine the density contrast and the density inhomogeneity length scale.
Although most of the TRACE observations suggest that the kink mode (m=1) is excited since the loop axis is displaced laterally, other modes are probably also excited, for example the sausage and fluting modes. The sausage or pulsating modes are not trapped for thin and long loops but they can exist in fat and dense flare or post-flare loops. Nakariakov et al. (2003) and Aschwanden et al. (2004) have suggested that these modes, having short periods of the order of seconds, have been detected in the radio band. Fluting modes are trapped even in thin and long loops with a frequency quite similar to the frequency of the kink mode. However, up to now there is no observational evidence of such modes.
Additionally, observations indicate that the most excited mode is the fundamental mode in the longitudinal direction, i.e. the loop has a node at the footpoints and a maximum at half the loop length. There is only one example of excitation of the first longitudinal harmonic (Verwichte et al. 2004). This mode has a maximum at a quarter of the loop length and a node at half the loop length (in addition to the nodes at the footpoints). The period of the mode is basically half the period of the fundamental mode, however longitudinal density stratification can modify this ratio. The difference between the real ratio and 2 has been used by Andries et al. (2005) as a seismological tool to obtain information about the density scale height in loops.
Although the main features of the rich spectrum of eigenmodes in a coronal magnetic tube are well know, the problem of the excitation of the eigenmodes for a given initial perturbation has been only partially addressed (see Terradas et al. 2006; Ruderman & Roberts 2006). In the present paper, this problem is investigated in detail and given an initial perturbation, the amount of energy that is distributed among the modes is calculated. We study how the energy of the trapped modes depends on the characteristics of the initial perturbation, namely the shape and location, and also on the loop properties. We focus on the energy that is trapped in the loop and not on the energy that is radiated away to wave leakage. This last issue has been investigated using straight loop models (see for example Terradas et al. 2005; Cally 1986; Pascoe et al. 2007; Terradas et al. 2007; Cally 2003) and also considering curved loops (Selwa et al. 2006; Díaz et al. 2006; Selwa et al. 2007; Brady et al. 2006; Verwichte et al. 2006a; Brady & Arber 2005-c).
This paper is organised as follows. In Sect. 2 a brief description of the tube eigenmodes is given and the energy associated with the different modes is introduced. In Sect. 3 we consider the initial value problem and we describe the method for calculating the energy trapped by the loop given an initial perturbation. The results for simple initial excitations are shown in Sect. 4. More complex perturbations are briefly studied in Sect. 5. In Sect. 6 the main conclusions are drawn.
The main properties of the loop eigenmodes using a straight cylindrical tube model have been studied by several authors (see for example Spruit 1981; Cally 1986,2003; Edwin & Roberts 1983). The dispersion diagram under coronal conditions and in the zero- limit indicates that only fast body waves are allowed. In the long-wavelength limit (i.e., when the tube radius is much smaller than the wavelength of the disturbance) the kink mode (m=1) and the fluting modes ( m=2, 3, 4,...) basically oscillate with the kink frequency. In contrast, the sausage mode (m=0) is leaky at this limit and is only trapped for wavenumbers larger than a cut-off value (see Roberts et al. 1984; Edwin & Roberts 1983).
The frequency of the modes, ,
is determined by solving the dispersion relation which is derived by imposing continuity of the total pressure perturbation, P, and radial velocity, v_{r}, at the loop boundary (r=R). The eigenfunction is written in terms of Bessel
and Hankel functions (J_{m} and H_{m}^{(1)}). In particular, the total pressure
perturbation is
(2) |
In order to compare different initial perturbations,
all perturbations must represent the same amount of
energy. The waveenergy density in the zero-
approximation is (see for example Bray & Loughhead 1974)
Figure 1: Energy density as a function of the radial position for the modes m=1, m=2 and m=3. In this plot the three modes have the same total energy E (energy density integrated over space). The loop length is L=50R and the density contrast is . The sausage mode is not represented since for these loop parameters this mode is not trapped. | |
Open with DEXTER |
From the observations we can estimate the amplitude
and the energy of the oscillating loops. According to the eigenfunction of the kink mode, the
radial displacement at the loop radius R is
We now turn to the theoretical study of how the amplitude of oscillation (i.e. the amplitude introduced in the previous section) depends on the initial disturbance. There are several ways to study the initial value problem. The first method is to solve the time-dependent problem numerically and to derive from the results the amount of energy that is deposited in the normal modes of the structure. The advantage here is that the behaviour of the tube is known at any time. Another method is to use analytical approximations of the solution. This has been done recently by Ruderman & Roberts (2006) who solved the initial value problem using the Laplace transform to determine the motion of the tube. As expected, these authors have found that the asymptotic behaviour of the loop for is given by the fundamental normal mode of the structure. Since we are mostly interested in the energy that is trapped in the loop and not in the transients before the loop settles in the normal mode, we will use the second approach.
Ruderman & Roberts (2006) concentrated on the kink oscillation (m=1). Here we
have extended the results to .
These authors studied initial perturbations
in the total pressure given by the function F
The initial disturbance excites fast waves that propagate towards and away from
the loop. Part of the energy is trapped by the density enhancement while most of
the energy does not reach the loop or is eventually radiated from the tube
(through the leaky modes). Ruderman & Roberts (2006) showed that for
the tube is oscillating with the eigenmode, representing the trapped
energy. Since the tube oscillates with the eigenmode the magnetic pressure is
simply given by Eq. (1). An interesting result of the
work of Ruderman & Roberts (2006) is that an expression for the amplitude is explicitly derived (they use
as the amplitude, in our notation
). The coefficient
depends on the particular form of
the initial perturbation and it is given by
(9) |
In the previous analysis the amplitude
and consequently also the energy deposited in a mode with fixed wavenumbers k and m has been calculated. However, in general a real perturbation is not restricted to a single k and m pair, but it requires a sum over several
different wavenumbers. If the form of the initial perturbation in the azimuthal
and vertical direction is
then we can write
(12) |
From the practical point of view, for a given form of the initial perturbation, P_{0}, we determine the corresponding decomposition as a sum of azimuthal numbers m and longitudinal wavenumbers . By evaluating Eq. (7) the amplitude ( ) associated with each pair of wavenumbers is calculated. First, we solve the dispersion relation and compute the eigenfrequency . Second, the derivative of the dispersion function is calculated and it is evaluated at . Finally, we calculate numerically the integral in the coefficient . Additionally, the energy of the mode is computed. This process is repeated for each set of values m and n.
The aim of this section is to investigate how the energy deposited in the eigenmodes depends on the properties of the initial disturbance and the properties of the loop. We first study in detail simple perturbations. More general excitations are addressed in the following section.
The initial perturbations that we study in this section are simple, in the
sense that they can be represented by a single k and m. Let us consider first
the kink mode (m=1) and the first fluting mode (m=2) and concentrate on the
fundamental longitudinal wavenumber .
For simplicity, we take the
perturbation (in the magnetic pressure) to have the following form in the radial
direction:
Figure 2: a) Example of kink perturbation and b) fluting (m=2) perturbation in the perpendicular plane to the loop axis. For this initial disturbance r_{0}=15R and a=2R. In the vertical direction the dependence is sin(kz). White colour corresponds to positive magnetic pressure perturbations while black colour represents negative values. The circular white region in the centre represents the loop. | |
Open with DEXTER |
In order to compare results for perturbations located at different
positions we need to choose the normalisation constant,
,
in such a way that the energy of the initial perturbations is the same.
For the perturbation given by Eq. (13), since
,
and
,
the total energy is
E_{0} | = | (14) | |
= | (15) |
We have focused on perturbations located in the coronal medium (
)
and we have studied how the displacement at the loop boundary and the total
energy change with r_{0} and a. In Fig. 3 these two quantities are
plotted for the m=1 mode. Both quantities decrease quite rapidly with r_{0} (specially the energy) while they smoothly increase with the width of the
initial perturbation a. This is a consequence of the convolution of the
initial perturbation and the eigenfunction. From Eq. (7) we find that,
since the contribution of the internal part of the loop (see Eq. (8)) is
negligible due to the assumption
,
the dependence with r_{0} and a is
Note that due to the strong dependence on r_{0} in general only a small amount of the energy of the initial perturbation is trapped in the loop. For example, for a perturbation located at a distance similar to the loop length r_{0}=40R (L=50R and a=R), the trapped energy is 9.4 10^{-7}E_{0}.
Figure 3: a) Displacement and b) trapped energy as a function of r_{0} and a for the m=1 mode. For these plots L=50R and . The behaviour of the displacement and energy with r_{0} and a is given by Eqs. (18) and (19) respectively in the limit . | |
Open with DEXTER |
We have done the same analysis for the first fluting mode (m=2) and we have found that the behaviour with r_{0} and a is qualitatively similar to the results for the kink mode, specially in the limit of large r_{0}. Since we have normalised the energy of the initial pulse we can compare the trapped energy of these two modes. The difference in energy between the modes is typically of three orders of magnitude. Therefore, fluting modes will be hardly excited in comparison with the kink modes. This is in accordance with Fig. 1.
We now study the influence of the loop parameters on the amplitude of the displacement and the trapped energy. In particular, we investigate the dependence on the ratio between the length and the loop radius, L/R, and the density contrast between the tube and the coronal environment, .
To study the effect of an initial perturbation on loops with different lengths we again normalise the initial perturbation using Eq. (16). Although the length of the system is now changed, the total energy of the initial perturbation is constant using this normalisation. In Fig. 4 the displacement and energy are plotted as a function of the loop length for different positions of the initial perturbation. There are several aspects of these curves that deserve comment. The displacement and energy have a maximum. This means that there are loops, with a preferred loop length ( ), that are able to trap energy more efficiently than others. The position of this maximum depends on the location of the initial perturbation, and we see from Fig. 4 that increases with r_{0} (see the curves with different line styles). Note that the displacement and energy do not change much with the loop length for long loops ( ). By contrast, short loops () are quite sensitive to the location of the initial perturbation.
Figure 4: a) Displacement and b) energy as a function of the loop length. For all the cases the width of the perturbation is constant a=R. Solid line corresponds to r_{0}=10R, dashed to r_{0}=20R and dotted-dashed corresponds to r_{0}=40R. For these plots . | |
Open with DEXTER |
The behaviour of the displacement as a function of L in two different limits can be determined analytically. Using the thin tube approximation () and the asymptotic behaviour of small arguments in the Bessel and Hankel functions (valid when ) we find that . This is the dependence present in Fig. 4a for . On the contrary, for short loops in comparison with the loop radius the thin tube approximation is no longer valid, but using the approximation for large arguments it is possible to show that , being a constant. Since the displacement tends to zero in the two limits, the curve must show a maximum in between (see Fig. 4).
Figure 5: a) Displacement at r=R and b) energy as a function of the density contrast. The perturbation is located at r_{0}=40R, a=3R. The solid line corresponds to a loop with L=50R while the dashed line corresponds to L=100R. | |
Open with DEXTER |
We have also studied how the results depend on the density contrast. Figure 5 shows the behaviour of the displacement and energy as a function of this parameter. Both the displacement and the energy decrease as is increased, which is the expected behaviour based on physical grounds. For the same initial perturbation high density contrast loops, i.e. basically heavy loops, are less excited than low density loops. We can see from Fig. 5a that, for example, loops with a density contrast of three oscillate, according to our model, with an amplitude that is three times larger than the amplitude of oscillation of loops with a density contrast of eight. Thus, the density contrast is another parameter that can have a strong influence on the amplitude of oscillation of loops. Note that the behaviour of the trapped energy with the density contrast is not in contradiction to the result that heavy loops radiate their energy more slowly than low density contrast loops (see for example Selwa et al. 2007; Verwichte et al. 2006b) due to the emission of leaky waves. This simply means that high density contrast loops are good waveguides because they have the eigenfunction very localised around the loop. But for the same reason, they are not very efficient in trapping energy from an external perturbation, since the trapped energy is proportional to the convolution between the initial perturbation and the eigenfunction.
In this section we investigate the extension of the method described before to more general and realistic perturbations. In particular, azimuthally and longitudinally localised perturbations are studied.
Let us assume that instead of single modes we have an external perturbation that
has the same radial dependence as in Eq. (13) but now it has
a Gaussian dependence on the coordinate, i.e. we consider that the
perturbation is
Figure 6: Example of an azimuthally localised perturbation. The circular white region represents the loop. For this initial disturbance r_{0}=15R, a=2R and . This disturbance is decomposed in Fourier sinus series using Eq. (10). | |
Open with DEXTER |
It is simple to show that the energy of this initial disturbance is
Figure 7: a) Energy of the Fourier components of the initial perturbation calculated using Eq. (21). For this plot r_{0}=40R, a=4R and (n=1, L=50R). b) Trapped energy of each eigenmode. Note that the vertical axis is in logarithmic scale. The difference in energy between consecutive modes increases with m. For example, the trapped energy for the m=1 is 4 10^{2} larger than the energy of the m=2. | |
Open with DEXTER |
We have studied a perturbation with r_{0}=40R, a=4R and . In Fig. 7a the energy of the modes in the angular decomposition is shown. For this particular perturbation the energy distribution in the initial perturbation decreases smoothly with the order of the harmonic m (although the largest coefficient, a_{m}, in the Fourier decomposition of the pressure perturbation corresponds to m=3). The results of trapped energy associated with each m are displayed In Fig. 7b. It is clear that the energy decreases with m and that the kink mode has the largest energy. This behaviour is a direct consequence of the result that the energy trapped by the kink mode is between two and three orders of magnitude larger than the energy of the first fluting mode (see Sect. 4.1). Even for m=1 the trapped energy is quite small (1.04 10^{-6}E_{0}) and, as expected, is lower than that for an excitation with a single m (we found for the same parameters that the energy was 3.77 10^{-6}E_{0}) since in that case all the initial energy was in one single mode. Qualitatively similar results are found when the width of the perturbation in the angular direction, , is changed. Obviously, the energy of the modes in the initial perturbation is different, but again the trapped energy distribution is dominated by the kink mode.
We can also study the energy distribution for a perturbation that has a Gaussian
dependence with the z coordinate,
Figure 8: Example of longitudinally localised perturbation. The white area represents the loop. For this initial disturbance r_{0}=15R, a=2R, and z_{0}=L/4. This disturbance is decomposed in sinus series using Eq. (10) with m fixed (m=1). | |
Open with DEXTER |
Now the energy of the disturbance is simply
We have focused first on a perturbation with m=1, r_{0}=40R, a=4R, and located at half the loop length, z_{0}=L/2. The energy of the modes in the initial perturbation is represented in Fig. 9a. Since the function is even in z there is only a contribution from the odd n ( ). Although the Fourier decomposition of this particular pressure perturbation is dominated by the mode with n=2, the most energetic is the n=5 (due to the dependence on n in the previous equation). Figure 9b shows that the fundamental mode is the most energetic and that the trapped energy decreases with the order of the harmonic (the difference in energy between consecutive modes is typically of the order 10^{2}).
We have also considered the case z_{0}=L/4, i.e. a perturbation located at a quarter of the loop length. The energy decomposition of the initial perturbation is represented in Fig. 10a. Again the kink mode has the largest energy as can be appreciated in Fig. 10b. Nevertheless, the differences in energy between the longitudinal harmonics are not as large as for the z_{0}=L/2 (now for example the fundamental and the first harmonic have energies of the same order of magnitude). This result also suggests that longitudinal harmonics are in principle more easily excited than azimuthal harmonics (with differences in the energy of the order 10^{3}).
Figure 9: a) Energy of the Fourier components (calculated using Eq. (23) with m fixed) of the initial perturbation given by Eq. (22). The m=5 has the largest energy in the decomposition. Only the terms with m odd contribute to the sum. b) Trapped energy of each eigenmode (the vertical axis is in logarithmic scale). For this case r_{0}=40R, a=4R (L=50R), z_{0}=L/2 and . | |
Open with DEXTER |
Figure 10: a) Energy of the Fourier components (calculated using Eq. (11) with m fixed) of the initial perturbation given by Eq. (22). The m=6 has the largest coefficient in the decomposition. b) Trapped energy of each eigenmode (the vertical axis is in logarithmic scale). Note that the fundamental mode has the largest energy (although the first harmonic has the maximum weight in the decomposition). For this case r_{0}=40R, a=4R (L=50R), z_{0}=L/4 and . | |
Open with DEXTER |
The energy that an initial disturbance stores in the eigenmodes of coronal loops has been examined. We have focused on external perturbations and we have found that when they are located far from the loop the trapped energy decreases quickly with the distance of the perturbation (as ). Even for perturbations near the loop the trapped energy is several orders of magnitude smaller that the energy of the initial perturbation. This is mainly because first, not all of the energy of the initial perturbation reaches the tube and second, part of the energy does not get trapped in the loop (and it is probably radiated away through the excitation of the leaky modes). Therefore, most of the energy is in fact in the so-called "leaky continuum'' (see Andries & Goossens 2007) rather than in the eigenmodes representing the trapped energy. However, the estimate of the amount of energy in a particular mode can be a useful parameter that can give information about the possible damping mechanisms acting in the solar corona. Note that in the case of resonant absorption, this quantity is basically the energy that it is eventually dissipated and converted into heat.
According to our simple model, loops with different lengths and densities can respond quite differently to an initial perturbation. The loops with an optimal length are able to trap more efficiently the energy of the perturbation. On the other hand, short loops are less favourable for oscillations. Additionally, their density also plays a significant role in the wave trapping problem. We have found that for the same initial perturbation the trapped energy decreases quickly with the density contrast. These results suggest that these two quantities, i.e. the loop length and density contrast are important factors for the amplitude of oscillation in loops and can help to explain why, as reported from TRACE observations (Aschwanden et al. 2002; Schrijver et al. 2002), some loops oscillate clearly while other loops in the same active region, and even closer to the location of the perturbation, do not show any signature of oscillation.
For different kinds of disturbances we have found that the kink mode always traps the largest amount of energy in comparison with other modes. For example, the energy of the kink mode is typically 10^{3} larger than the energy of the first flute mode suggesting that these modes will rarely be excited for perturbations located far from the loop. This explains why it is justified to concentrate on the kink mode of oscillation since the contribution of the flute modes is rather small. Additionally, this could explain the fact that up to now no flute modes have been detected. However, other explanations are possible, for example that the motions of the fluting modes, displacing the loop boundary but not the loop axis, are not spatially resolved with the actual instruments. In contrast, the longitudinal harmonics are more easily excited. We have shown that for a perturbation localised close to one of the footpoints, the fundamental and first harmonic have energies of the same order of magnitude. Thus, from the observational point of view these modes are likely to be detected. In fact, there is an example of the detection of such a mode (see Verwichte et al. 2004; Andries et al. 2005).
There are estimations of the various types of energy released in real flares. However, it is not clear what amount of energy is associated with MHD waves. From the information of our model and the information of the coronal loop oscillations, this energy can be estimated. We have found that a typical oscillating loop has an energy of . If the perturbation is located at the loop length, then the energy of the initial disturbance should be 10^{6} times larger, i.e. of the order . This quantity is, in order of magnitude, similar to the estimations of the other types of energy. Thus, the energy of the oscillating loops can be potentially used as a seismological tool to determine the amount of energy released in real flares (related to fast MHD waves).
The model used here has several limitations. First, all the magnetic field lines have the same length, which is not realistic if we compare with the three-dimensional structure of bipolar regions. Second, we have studied the linear regime, however in a real flare a shock wave is produced and nonlinear effects are important and also couplings between the different harmonics are possible. An additional restriction of our study is that we have analysed only perturbations in the magnetic pressure, but in a real explosive event the initial disturbance can be produced for example by a density enhancement. In this regard, more general perturbations, especially those that are able to also excite the sausage modes need to be considered.
Acknowledgements
We thank the referee for comments which helped to improve this paper. J. Terradas is grateful to the K.U. Leuven (Department of Mathematics) for a fellowship. He also thanks the Spanish Ministry of Education and Science for the funding provided under the Juan de la Cierva program and for funding provided under grant AYA2006-07637.