Issue |
A&A
Volume 502, Number 1, July IV 2009
|
|
---|---|---|
Page(s) | 315 - 323 | |
Section | The Sun | |
DOI | https://doi.org/10.1051/0004-6361/200811405 | |
Published online | 16 April 2009 |
The effect of elliptic shape on the period ratio P1/P2 of emerging coronal loops
R. J. Morton - R. Erdélyi
Solar Physics and Space Plasma Research Centre (SP2RC), University of Sheffield, Hicks Building, Hounsfield Road, Sheffield S3 7RH, UK
Received 24 November 2008 / Accepted 28 March 2009
Abstract
Aims. We determine the effect of an elliptical shape on the period ratio for the standing transversal oscillations of a longitudinally stratified coronal loop throughout its emergence from the low solar atmosphere into the ubiquitously magnetised corona.
Methods. Under the assumption that elliptical curvature has a negligible effect on eigenfrequencies, the equation that describes the projection of a density profile onto a magnetic flux tube with elliptical shape is obtained in a gravitationally stratified atmosphere. The effect of the elliptical shape on the period ratio of the fundamental mode to the first harmonic (
P1/P2) at various stages of emergence is determined, assuming that the oscillation periods are much shorter than the characteristic time scale of loop emergence.
Results. We find that there are two separate cases of elliptical shape that occur, the minor ellipse and the major ellipse. It is then shown how the period ratio
P1/P2 is dependent upon the ellipticity (), the parameter characterising the stage of emergence (
)
and the density scale height (H). Ellipticity is found to make an important contribution to
P1/P2 for the minor ellipse when compared to its counterpart of standing oscillations of stratified loops with semi-circle or circle-arc shape. The major ellipse was found to have a lesser effect on the period ratio of standing oscillations. We also find the value of
P1/P2 is dependent upon the stage of emergence of the loop, where the greatest contribution from emergence to the ratio of
P1/P2 is when the loop is almost fully emerged. The important implication for magneto-seismological interpretations of the observations of oscillating coronal loops is that measurements of ellipticity and stage of emergence should supplement observations of oscillation periods and should be considered when applying observed frequencies of the fundamental mode and first harmonic to determine the diagnostic properties of these oscillating loops, e.g. the density scale height or strength of magnetic field. Neglecting the determination of ellipticity and stage of emergence may result in a
error in estimating density scale height.
Key words: magnetohydrodynamics (MHD) - plasmas - Sun: corona - waves
1 Introduction
In recent years it has been possible to study the plasma fine structure of the complex and dynamic magnetised solar atmosphere due to the relatively high-resolution imaging and spectroscopic data obtained from the Solar and Heliospheric Observatory (SOHO), the Transitional Region and Coronal Explorer (TRACE) and Hinode. Various magnetohydrodynamic (MHD) waves (e.g. longitudinal, transverse) have been observed and identified within a variety of open (e.g. sunspots, pores, plumes) and closed (e.g. prominences, coronal loops) magnetic structures (see Banerjee et al. 2007; Nakariakov & Verwichte 2005 for the latest reviews of observational results with plenty of further references therein). The study of these magnetic oscillations, in order to derive information on the diagnostic properties of the plasma, is known as solar magneto-seismology and the methodology was first suggested by Uchida (1970) and Roberts et al. (1984) who initially used the term of coronal seismology for specific coronal applications. The field has rapidly developed and the method of magneto-seismology is now an integral acquisition tool for developing our understanding of the large range of dynamical events that are detected in the solar atmosphere, often with applications to the problem of coronal heating (Erdélyi 2006a,2008,2006b; Taroyan 2008).
![]() |
Figure 1:
a) Sketch showing the equilibrium state of an elliptically
curved, magnetic flux tube with plasma density
|
Open with DEXTER |
Measured properties of the MHD oscillations (such as period, spatial distribution of velocity and/or amplitudes, damping rate) may be inverted to retrieve information about otherwise unmeasurable or hardly measurable solar atmospheric parameters (e.g. gravitational or magnetic stratification, structuring, loop expansion, cross-sectional shape, loop geometry, wave leakage, magnetic twist, loop curvature). The basic theoretical model of oscillations in coronal loops, put forward in a seminal paper by Edwin & Roberts (1983), was an ideal and linear MHD theory of a straight and homogenous cylindrical magnetic flux tube. However, SOHO, TRACE and Hinode observations clearly show that the properties of oscillations of various solar magnetic loops can be significantly different from those predicted by the homogenous cylindrical tube model. The source of these differences is thought to be due to various geometric and physical (often second-order) effects, as listed above (e.g. stratification, expansion, structuring, shape, wave leakage, variable background on the time scale of oscillations, cross section geometry, etc. to name a few).
It is clear there are a number of effects to be considered when interpreting solar flux tube wave or oscillation observations using the tools of magneto-seismology. However, there are at least two key measurable quantities that one can use for the magneto-seismology theory: period (or frequency) ratio P1/P2 of the fundamental mode to the first harmonic of oscillations (i.e. temporal seismology) and amplitude distribution of oscillations along the flux tube (i.e. spatio-seismology). For a review on magneto-seismology of kink oscillations see Ruderman & Erdélyi (2009).
Erdélyi & Verth (2007); Verth et al. (2007) found that the spatial position along the loop of the maximum amplitude of the first harmonic of the standing fast kink mode varies when the loop is stratified (referred to as anti-node shift). The predicted anti-node shifts calculated in Verth et al. (2007), are between 600 to 6000 km depending on the density profile used. At present there has been no observation yet of the anti-node shift due to the low spatial resolution of current observational instruments.
The period ratio was first measured by Verwichte et al. (2004), who observed two values in different loops,
P1/P2=1.81
0.25 and 1.64
0.23. These values are in contrast to the theoretically derived value for the period ratio using the uniform tube model, where
.
The systematic deviation between the theoretical and observed values is a serious inspiration to refine modelling efforts, just like such ``discrepancies'' were the
driving impetus for helioseismology (see, e.g. Erdélyi 2006a,b). Andries et al. (2005a,b) and Goossens et al. (2006) suggested that the ratio of the periods of the fundamental mode to the first harmonic,
P1/P2, could be used as a tool to determine the density scale height in the corona. For a specific review on the application of the period ratio
as a tool for coronal seismology see Andries et al. (2009).
To the best of our knowledge, stratification (Andries et al. 2005b; McEwan et al. 2006) and loop expansion due to a non-constant magnetic field (Verth 2007; Verth & Erdélyi 2008; Ruderman et al. 2008) so far seem to be the main causes of the most dominant deviation between the measured and theoretically predicted periods of oscillations. However, magnetic structuring and gravitational stratification are, in a sense, in competition with each other, since their effect on the period ratio of P1/P2 is opposite, i.e. the tube expansion parameter discussed by Verth & Erdélyi (2008) increases the period ratio ( P1/P2>2), while the longitudinal stratification causes the period ratio to decrease ( P1/P2<2). This raises the question: what physical effect(s) or perhaps loop properties are then responsible for the deviation of the measured period ratio from those predicted by magneto-seismology? Could the answer be linked to loop geometry, e.g. curvature, shape, cross-section, stage of emergence, etc.?
There have been a number of other investigations into how the different geometries of coronal loops affect MHD oscillations. Dymova & Ruderman (2006) investigated the dependence of
P1/P2 on
circular shape loop geometry for a longitudinally stratified loop. This idea was underpinned by the first impression of oscillating coronal loops appearing to have circular shape. They modeled the loop as an arc of a circle with curvature radius
with the center at a distance l from the photosphere. It was found that the shape of the loop, determined by the parameter
,
had a marked effect on the ratio
P1/P2. They concluded that in order to determine the scale height from the observed values of
P1/P2, a more precise knowledge of the loop shape is needed.
Investigations of elliptical cross-sectional geometry of loops in cold plasma by Ruderman (2003) and in a finite plasma-
by Erdélyi & Morton (2009), has shown that magnetic loops can support
additional transverse modes compared to the cylindrical model. However, it is shown in Erdélyi & Morton (2009) that the period ratio is practically unaffected by the ellipticity of the loop cross-section. In spite of the high spatial and temporal resolution of the current cohort of instrumentation, this latter theory has yet to be put to the test. We can conclude that curvature and cross-section are less dominant parameters for temporal magneto-seismology in coronal
applications, so we shall ignore their geometric properties in what follows.
In this paper we investigate the effect of loop shape on the period ratio
P1/P2 by examining loops which deviate from a circular shape. Here we analyse the standing
oscillations of loops with an elliptical shape. Further to this, we also investigate how
P1/P2 is influenced by the different consecutive stages of emergence of an elliptical loop. It is assumed that the characteristic time of emergence is much greater than the period of loop oscillations. This assumption is verified observationally as the loop emergence takes place on
day-scale, while transverse oscillations are reported with periods in the sub-hour domain. The projection of an exponential density profile on to a loop with an elliptical shape is calculated and the Sturm-Liouville equation describing a stratified magnetic flux tube (Dymova & Ruderman 2005; Erdélyi & Verth 2007) is solved numerically. The effect on the period ratio
P1/P2 is then
calculated for different values of ellipticity at various stages of emergence. The magnitude of the effect of ellipticity is then compared with the circle-arc model (Dymova & Ruderman 2006). Beside the
theoretical justification given above, in this particular aspect, recent progress by Aschwanden et al. (2008) in 3D reconstruction of loops using the Solar TErrestrial RElations Observatory (STEREO), has provided new and interesting information on the loop geometry. It was found that the majority of the loops studied were in fact non-circular in shape, with a deviation from circularity up to .
2 Modelling of coronal loop geometry
In this paper we consider a thin coronal loop with a circular cross-section which has elliptical shape as shown in Fig. 1a. The effects of circular curvature on the eigenfrequencies of transverse oscillations have been investigated and found to be of the order
,
where R is the radius of the loop cross-section and 2L is the loop length (Terradas et al. 2006; Van Doorsselaere et al. 2004. Note, caution should be exercised as there is a flaw in the derivation in Van Doorsselaere et al. 2004). It was noted by Terradas et al. (2006) that one of the main differences between the straight and curved cylindrical magnetic flux tubes was that curvature introduces preferential directions of oscillations. The introduction of a curved shape on the cylindrical tube changes the problem from 1-dimensional to 2-dimensional. For example, two distinct kink modes can occur in a curved loop, polarised along the vertical (z' axis in Fig. 1a) and along the horizontal axis (x' axis in Fig. 1a). The periods of these two kink modes are similar to each other and to the period of the kink mode found in the model of a straight magnetic tube of Edwin & Roberts (1983). Most importantly, in the context of this study, in the thin tube limit (
)
the difference between these modes is negligible, so a straight tube approximation can be used when studying curved loops. Dymova & Ruderman (2006) applied this ``straightening'' approximation to study the effects of circular shape on the transverse oscillations of longitudinally stratified loops.
A study of oscillations in a magnetic tube with semi-elliptical shape is possible using a variation of the toroidal coordinate system. However, the model appears to have no analytic solution due to the dependence of the equilibrium quantities on all the spatial coordinates. Since Díaz (2006) found that the effects of curvature on oscillations of an elliptically curved magnetic slab are negligible, this may suggest that one would find similar results for a semi-elliptical magnetic tube. In this paper, we progress under the assumption that the effect of elliptical curvature is negligible in the thin tube approximation and leave a study, analytic or numerical, of the semi-elliptical loop till a later date.
Consider a straight magnetic flux tube with a cross-section of radius R and length 2L which is embedded in a magnetic environment (see Fig. 1b). The magnetic field inside
and outside the tube is uniform, with the internal (external) magnetic field of the form
(
). Here,
is the unit vector in the z-direction and we assume
.
The plasma inside and outside the tube is an ideal and cold plasma (i.e.
,
where
is the plasma beta defined by
). The cold plasma approximation is valid in the corona because the characteristic value of
.
Linear perturbations about the equilibrium for a cold plasma are described by the linear, ideal MHD equations

Here









It was shown in Dymova & Ruderman (2005) and Erdélyi & Verth (2007), for the first order approximation with respect to the parameter
where
(thin tube approximation), that the solutions to
Eqs. (1) for the squared frequencies of the kink tube oscillations are the eigenvalues of the Sturm-Liouville problem
where
![]() |
(3) |
Here A(z) is the magnetic pressure perturbation and





A comment by Dymova & Ruderman (2006), concerning the radial dependence of density, should also be repeated. From equilibrium conditions, it can be shown that the dependence of density perpendicular to the magnetic field can be chosen arbitrarily. However, we need to consider the density variation in the radial direction as a discontinuous jump between definite magnetic surfaces. A continuous density profile would cause the oscillations of the coronal loop to experience resonant damping (see, e.g. Ruderman & Roberts 2002), however this is well beyond the scope of the current study.
3 Elliptically shaped loops with plasma density exponentially decreasing with height
Now consider a coronal loop with a plasma density that is exponentially decreasing with height, i.e.
where



In order to allow us to make a direct prediction of observables, it is necessary that the values of
used reflect the possible range of values available in the corona. A relationship between temperature and hydrostatic scale height for plasma in hydrostatic equilibrium is
given by (e.g. Aschwanden 2004),
where MK = 106 K. For hot coronal loops with









![]() |
Figure 2:
Sketch showing a partially emerged loop with minor elliptical shape. The geometry of the shape is described by the arc of an ellipse with a minor half-axis of length b and a major half-axis of length a, the height of the loop apex the photosphere |
Open with DEXTER |
Table 1: Approximate values of ellipticity derived for the STEREO loops observed by Aschwanden et al. (2008).
![]() |
Figure 3:
a)- c) The ratio of the periods of the fundamental mode to the first harmonic,
P1/P2, as a function of density stratification parameter,
|
Open with DEXTER |
Recent observations using STEREO (Aschwanden et al. 2008), have provided an insight to the extent of the ellipticity of loops in the corona. Seven loops were reconstructed in three-dimensions and a circle, with radius
,
was fitted through the loop footpoints and the loop top above the
midpoint. The maximum (
)
and minimum radii (
)
of the loops, measured from the center of the fitted circle, were determined (see Table 1 for observed values). The ratio of these radii (
)
to the radius of the fitted circle was taken to define the circularity of loops, i.e.
.
The largest value obtained was
(loop No. 7), a
departure from circularity. Although this is not a strict measure of ellipticity, it clearly shows the departure of
the loops from the circular shape. Taking e.g.
and
from loop 7, we obtain a rough estimate for the ellipticity of
.
The other estimated values of ellipticity range between
0.5-0.76 for the loops 1-7 and the values are shown in Table 1. Caution should be taken with these values as the fitting procedure for the circle does not take into account the inclination of the loop apex from the vertical axis. This means
and
,
in some cases, are larger and smaller, respectively, than the actual maximum and minimum radius. The estimate of the ellipticity is then larger than the actual value. Nonetheless, taking these estimates from the observations into account, we conjecture that loops with large ellipticities would probably be rare within the corona. A typical value of ellipticity that is considered here for the minor elliptical loop is
.
For a loop of total length
approximately 120 Mm, the length of the minor half-axis would be
Mm and major half-axis
Mm for this value of ellipticity.
3.1 Minor elliptical loop
Consider a loop with elliptical shape with a given state of emergence as shown in Fig. 2. There are a number of parameters that characterise the problem and should be central in a
quantitative description of standing oscillations in the prescribed geometry. One parameter describes the state of emergence of the loop, denoted by
![]() |
(6) |
Here l is the distance of the loop center from the photosphere, b is the length of the minor half-axis of the ellipse and




![]() |
(7) |
where a is the length of the major half-axis. A value of



The projection of the exponential density profile, Eq. (4), on to the minor elliptical loop gives,
where



where t1, t2 are parametric angles related to the polar angles


The integrals in Eq. (9) are known as incomplete elliptic integrals of the second kind. It is well known that the elliptic integrals in Eq. (9) cannot be expressed in terms of elementary functions, however, various approximations of the integrals do exist (see, e.g. Kennedy 1954). The substitution of Eqs. (8) in to (2) gives
where




![]() |
(12) |
which is the equation found by Dymova & Ruderman (2006) for the circular loop. Further to this, in the limit

![]() |
(13) |
which is the equation for a straight magnetic flux tube with a homogenous density profile, the eigenvalue solution to which is
![]() |
(14) |
where m is the mode number.
![]() |
Figure 4:
Showing the value of
P1/P2 as a function of the loop shape parameter, |
Open with DEXTER |
![]() |
Figure 5:
Shown is the percentage difference between the values of
P1/P2 when comparing the minor elliptical and the circle-arc geometries for a range of
|
Open with DEXTER |
In order to progress in solving Eq. (11), we proceed by solving the equation numerically
using the shooting method (see, e.g. Zwillinger 1989). This involves selecting an initial value for
and then solving Eq. (11).
is then adjusted till the boundary
conditions are satisfied.
In Figs. 3a-c, the ratio of
P1/P2 is plotted against the ratio of the height of the loop apex above the photosphere to the scale height (
). The different panels show the period ratio of a coronal loop at various stages of emergence, where we assume the
characteristic time for emergence is much greater than the period of an oscillation within the loop. Figure 3a shows
P1/P2 for a loop at an early stage of emergence (
), Fig. 3b is a half emerged loop (
), i.e. the semi-ellipse, and Fig. 3c is a loop at a late stage of emergence (
). It is clear from Fig. 3 that the ratio of
P1/P2 for the elliptical loop is dependent upon the level of stratification (
)
and ellipticity (
). Although perhaps not immediately obvious, the ratio
P1/P2 is also dependent on the stage of emergence (
). The effect of the ellipticity, for small values of
,
is to increase the value of
P1/P2 in comparison to the circle arc model. This is due to the loop apex becoming flatter as the ellipticity increases, hence the density along the loop takes on a profile more similar to that of the linear profile suggested by Erdélyi & Verth (2007) and reduces the amount of stratification. Further to this, it can be seen in Fig. 3c that as
increases,
P1/P2 for the elliptical loops tends to the values obtained for the circle-arc model. This also occurs for the
cases, however, it is not so easily seen. Also in Fig. 3c the value of
P1/P2 for the elliptical loops can be seen to decrease below the circle arc model. The reason for this is because when the value of
becomes sufficiently large,
P1/P2 reaches an almost constant value for all further values of
and hence, is almost independent of
.
For the circle arc model this happens at a lower value of
than for the elliptic model. The value of
P1/P2 for the elliptic model then continues to decrease for increasing
,
while
P1/P2 stays almost constant for the circle arc model. We can conclude from this that the effect produced by the shape of the loop on
P1/P2 is limited by the scale height, for small enough values of H there is a maximum effect from loop shape. Further, the effect of the elliptical shape
on
P1/P2 is felt at smaller values of scale heights than the circle-arc model.
The effect of ellipticity is shown explicitly in Fig. 3d, where
P1/P2 is plotted against the ellipticity for different arbitrarily fixed though observationally
typical values of
.
As
,
the value of
P1/P2 tends to 2 for all values of
.
This is expected because as the ellipticity tends to high values, the length of the minor axis is much less than the major axis (i.e.
). The coronal loop could then be approximated by a straight tube lying horizontally along the surface of the photosphere. The tube would then have an almost homogenous density profile hence the ratio of
P1/P2 tends to the value obtained for the homogenous tube, i.e. 2.
The effect of emergence, quantised by ,
on loop oscillations is emphasised in Fig. 4. For
(circular-arc) the effect of
on
P1/P2 is
almost linear, with the value of
P1/P2 decreasing as the loop emerges. Increasing the ellipticity causes the effect of
to have a parabolic profile which has minimum effect on the period ratio around
.
This suggests that the degree of stratification in the loop decreases as the loop starts to emerge (
), and once it is half emerged (
)
the degree of stratification starts to increases as the loop enters the late stage of emergence (
). It is also clear from Fig. 4 that during an early stage of emergence, the circular shaped loop has greater degree of stratification than the elliptic shape. This is expected as the elliptically shaped loop is flatter at the loop apex. However, during the later stages of emergence, the stratification of the elliptically shaped loop becomes greater. This happens at smaller values of
for decreasing values of scale height (increasing values of
).
Let us now study how the modelling of a loop with an elliptic shape compares to modelling a loop with a circular shape in regards to the ratio of the periods of the fundamental and first harmonic
oscillations. The relative percentage difference of the period ratio between the circle-arc and elliptical models is plotted in Fig. 5. It can be seen clearly here that the effect of ellipticity compared to the circle-arc model is not constant for different values of
.
There is a maximum value for the difference between the elliptical and circle-arc model, occurring at different values of
for the various values of
.
The observations of multimode oscillations by, e.g. Van Doorsselaere et al. (2008), have a measurement error of approximately
,
and it is seen in Fig. 5 that the percentage difference between the circle-arc and elliptical models can be greater than this error. The ellipticity of the coronal loop should then be measured and determined when trying to deduce the density scale height from multimode oscillations.
From the observations of transverse oscillations by Van Doorsselaere et al. (2008), we can obtain estimates of the density scale height. Van Doorsselaere et al. (2008) calculated a density scale height of 109 Mm, a value which is almost twice that of the expected hydrostatic value. This estimate was debated when loop expansion was also taken into account e.g. Ruderman et al. (2008) and Verth et al. (2008). Ruderman et al. (2008) estimated the scale height for a range of values of tube expansion ()
and found the estimated the values of scale height were smaller than estimates without expansion, e.g. for
the estimated scale height is
Mm, and for
the scale height is
Mm. We now consider how the elliptical shape of the loop combined with loop expansion effects the estimated value of scale height. The result for the Van Doorsselaere et al. (2008) observations with two different values of
are shown in Fig. 6 (labeled vD). Note that values in Fig. 6 are calculated for
held constant, this means that as
changes so does the length of the loop. The ellipticity can be seen to have a salient effect on the estimation of the scale height, e.g. for
and
the scale height estimate is
Mm, while for
(
)
Mm. This latter example shows there could be up to a factor of 1/3 uncertainty in obtaining the scale height. In Verth et al. (2008), a value of
P1/P2=1.54 was observed and the loop under observation
had a measured
.
Assuming a semi-circular loop, the scale height was calculated as H=7.65 Mm. The estimates of scale height with the inclusion of ellipticity is shown in
Fig. 6 (labeled VEJ). The ellipticity now produces an uncertainty of factor 1/4 in the density scale height, which is about 2 Mm. Either way, whichever measurement is correct, we can conclude that the elliptical shape has a profound importance in determining the scale height from the period ratio of oscillating loops.
![]() |
Figure 6:
The calculated values of density scale height, taking into account
elliptical loop shape, using the observations of Van Doorsselaere et al. (2008) with corrections due to tube expansion ( |
Open with DEXTER |
3.2 Major elliptical loop
Consider a loop with elliptical geometry as shown in Fig. 7a. The stage of emergence is now denoted by
, where
![]() |
(15) |
Here




The function

![]() |
(17) |
where t1 is given in Eq. (10) and
![]() |
(18) |
The substitution of Eqs. (16) into (2) gives
where

As with Eq. (11), (19) has no analytical solution. Solving numerically, the ratio of periods
P1/P2 is plotted against
in Fig. 7b. It can be concluded from Fig. 7b that the major elliptical loop has a much smaller effect on the period ratio of
P1/P2 than the minor elliptical loop. The major elliptical loop also causes the value of
P1/P2 to decrease compared to the circle arc model. This is a qualitatively different result from the model of oscillating minor elliptical loops. For values of
there is a very small deviation from the semi-circular loop. A value of
for a loop of typical length 120 Mm would have
Mm and
Mm, and for
,
Mm and
Mm. As mentioned earlier, we conjecture that such a strong ellipticity, i.e.
,
would probably be unlikely in coronal loops. For the early and late stages of emergence (e.g.
,
-0.5) a similar feature is found as the semi-emerged ellipse, i.e. only at large values of ellipticity do we see appreciable difference from the circle-arc case in the period ratio of standing oscillations.
![]() |
Figure 7:
a) Sketch showing the a loop with major elliptical shape.
The shape is described by the arc of an ellipse with a minor half-axis of length b and a major half-axis of length a, the height of the loop apex above the photosphere, |
Open with DEXTER |
4 Conclusions
In this paper we investigated the period ratio of transverse oscillations of a coronal loop with an elliptical shape, at various stages of its emergence from the sub-photosphere into the solar corona. The effect of curvature of the loop on eigenfrequencies of standing oscillations was neglected and we concentrated on the influence of the loop shape and the state of emergence. This paper generalises earlier studies, e.g. Andries et al. (2005a) and Dymova & Ruderman (2006), who investigated the effect of a circular loop shape on the parameter P1/P2 (the period ratio of the fundamental mode to the first harmonic of standing oscillations).
In order to model an oscillating coronal loop with an elliptical shape embedded in a gravitationally stratified solar plasma, an exponential density profile was projected onto an elliptical shape determining the loop and the parameters found to characterise the geometry were the ellipticity (), the ratio of loop height above the photosphere to the scale height (
)
and the measure of the stage of emergence (
). We considered two types of elliptical loop that can occur, the minor elliptical loop where minor axis of the ellipse is the vertical axis of the loop (see Fig. 2), and the major elliptical loop where major axis of the ellipse is the vertical axis of the loop (see Fig. 7a).
Numerical solutions of the Sturm-Liouville problem, as shown by Eq. (2), were found for both minor and major elliptical loops, and the relationship between
P1/P2,
,
and
was shown in detail.
We found that for the minor ellipse, a situation that occurs most plausibly under coronal conditions, increasing ellipticity increases the value of the ratio
P1/P2 compared to the
circular arc model for all values of .
It was also discovered that as
reaches a large enough value, the value of
P1/P2 becomes almost independent of
.
This implies that loop shape has a limited effect on longitudinal stratification and
P1/P2 under circumstances only when
.
The effect of the elliptical shape on
P1/P2 was found to be felt at larger values of
than the circle arc model (see, e.g
Fig. 3c).
The difference between the period ratio
P1/P2 of standing kink oscillations of loops with elliptical and circular shapes was demonstrated in Fig. 5, and it shows that
the percentage difference can reach %. This means the effect of ellipticity on standing oscillations could be measured even by the current cohort of high-resolution instrumentation.
Further to this, Fig. 4 also indicates that
P1/P2 is strongly effected by the value of
(stage of emergence of the loop), with the strength of the effect increasing as
and the scale height (
)
increases. It can easily be calculated that a loop which is in a late stage of emergence (i.e.
)
can show a deviation of up to 10% in
P1/P2 from that of the semi-emerged loop (
). Not including the stage of emergence could lead to an error of up to
upon estimating the density scale height. For a coronal loop of typical length (100 Mm) embedded in a gravitationally stratified atmosphere this could mean an error of around 6 Mm when estimating the scale height.
The major ellipse, however, is found to have a less pronounced effect on
P1/P2 compared to the circle arc case. Only at high values of ellipticity (
)
is the difference between
the frequency ratio of the major ellipse and circle arc models noticeable. As discussed in Sect. 3.1, there is observational evidence that loops will not have such a high degree
of ellipticity and the greatest value of ellipticity estimated for the seven loops observed by Aschwanden et al. (2008) provides initial evidence for this claim.
The present error bars on measurements of multimode oscillations are approximately 3% (private communication with authors of Verth et al. 2008). The difference between the minor ellipse and circle arc surpassed this for reasonable values of ellipticity, e.g. those derived from STEREO observations. Therefore, the ellipticity of the loop is an important parameter along with the stage of emergence and should be taken into account when using observational data from multimode oscillations to estimate the density scale height and ultimately determining the solar magnetic field by the technique of magneto-seismology.
Acknowledgements
The authors thank M. Ruderman and G. Verth for a number of useful discussions. R.E. acknowledges M. Kéray for patient encouragement. The authors are also grateful to NSF, Hungary (OTKA, Ref. No. K67746) and the Science and Technology Facilities Council (STFC), UK for the financial support they received.
References
- Andries, J., Arregui, I., & Goossens, M. 2005a, ApJ, 624, L57 [NASA ADS] [CrossRef]
- Andries, J., Goossens, M., Hollweg, J. V., Arregui, I., & Van Doorsselaere, T. 2005b, A&A, 430, 1109 [NASA ADS] [CrossRef] [EDP Sciences]
- Andries, J., Van Doorsselaere, T., Roberts, B., et al. 2009, Space Sci. Rev., submitted (26 pages) (In the text)
- Aschwanden, M. J. 2004, Physics of the Solar Corona (Praxis Publishing) (In the text)
- Aschwanden, M. J., Wülser, J.-P., Nitta, N. V., & Lemen, J. R. 2008, ApJ, 679, 827 [NASA ADS] [CrossRef] (In the text)
- Banerjee, D., Erdélyi, R., Oliver, R., & O'Shea, E. 2007, Sol. Phys., 246, 3 [NASA ADS] [CrossRef]
- Díaz, A. J. 2006, A&A, 456, 737 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Dymova, M. V., & Ruderman, M. S. 2005, Sol. Phys., 229, 79 [NASA ADS] [CrossRef] (In the text)
- Dymova, M. V., & Ruderman, M. S. 2006, A&A, 459, 241 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 179 [NASA ADS] [CrossRef] (In the text)
- Erdélyi, R. 2006a, Phil. Trans. Roy. Soc. A, 364, 351 [NASA ADS] [CrossRef]
- Erdélyi, R. 2006b, in Proceedings of SOHO 18/GONG 2006/HELAS I, Beyond the spherical Sun, ed. K. Fletcher, & M. J. Thompson, ESA-SP, 624, 15
- Erdélyi, R. 2008, in Waves and Oscillations in the Solar Atmosphere, Princeton Series in Astrophysics, ed. B. N. Dwivedi, & U. Narain, 61
- Erdélyi, R., & Morton, R. J. 2009, A&A, 494, 295 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Erdélyi, R., & Verth, G. 2007, A&A, 462, 743 [NASA ADS] [CrossRef] [EDP Sciences]
- Goossens, M., Andries, J., & Arregui, I. 2006, Phil. Trans. Roy. Soc. A, 364, 433 [NASA ADS] [CrossRef] (In the text)
- Kennedy, E. C. 1954, The American Mathematical Monthly, 61, 613 [CrossRef] (In the text)
- McEwan, M. P., Donnelly, G. R., Díaz, A. J., & Roberts, B. 2006, A&A, 460, 893 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Nakariakov, V. M., & Verwichte, E. 2005, Liv. Rev. in Sol. Phys., 2, 3 [NASA ADS]
- Roberts, B., Edwin, P. M., & Benz, A. O. 1984, ApJ, 279, 857 [NASA ADS] [CrossRef] (In the text)
- Ruderman, M. S. 2003, A&A, 409, 287 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Ruderman, M. S., & Roberts, B. 2002, ApJ, 577, 475 [NASA ADS] [CrossRef] (In the text)
- Ruderman, M. S., & Erdélyi, R. 2009, Space Sci. Rev., submitted (In the text)
- Ruderman, M. S., Verth, G., & Erdélyi, R. 2008, ApJ, 686, 694 [NASA ADS] [CrossRef] (In the text)
- Taroyan, Y. 2008, in IAU Symp., Waves and Oscillations in the Solar Atmosphere: Heating and Magneto-Seismology, ed. R. Erdélyi, & C. A. Mendoza-Briceño, IAUS, 247, 184 (In the text)
- Terradas, J., Oliver, R., & Ballester, J. L. 2006, ApJ, 650, L91 [NASA ADS] [CrossRef]
- Testa, P., Peres, G., Reale, F., & Orlando, S. 2002, ApJ, 580, 1159 [NASA ADS] [CrossRef] (In the text)
- Uchida, Y. 1970, PASJ, 22, 341 [NASA ADS] (In the text)
- Van Doorsselaere, T., Debosscher, A., Andries, J., & Poedts, S. 2004, A&A, 424, 1065 [NASA ADS] [CrossRef] [EDP Sciences]
- Van Doorsselaere, T., Nakariakov, V. M., & Verwichte, E. 2008, in IAU Symp., Waves and Oscillations in the Solar Atmosphere: Heating and Magneto-Seismology, ed. R. Erdélyi, & C. A. Mendoza-Briceño, IAUS, 247, 140 (In the text)
- Verth, G. 2007, Astron. Nacht., 328, 764 [NASA ADS] [CrossRef] (In the text)
- Verth, G., & Erdélyi, R. 2008, A&A, 486, 1015 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Verth, G., van Doorsselaere, T., Erdélyi, R., & Goossens, M. 2007, A&A, 475, 341 [NASA ADS] [CrossRef] [EDP Sciences]
- Verth, G., Erdélyi, R., & Jess, D. B. 2008, ApJ, 687, L45 [NASA ADS] [CrossRef] (In the text)
- Verwichte, E., Nakariakov, V. M., Ofman, L., & DeLuca, E. E. 2004, Sol. Phys., 223, 77 [NASA ADS] [CrossRef] (In the text)
- Zwillinger, D. 1989, Handbook of differential equations (New York: Academic Press) (In the text)
All Tables
Table 1: Approximate values of ellipticity derived for the STEREO loops observed by Aschwanden et al. (2008).
All Figures
![]() |
Figure 1:
a) Sketch showing the equilibrium state of an elliptically
curved, magnetic flux tube with plasma density
|
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Sketch showing a partially emerged loop with minor elliptical shape. The geometry of the shape is described by the arc of an ellipse with a minor half-axis of length b and a major half-axis of length a, the height of the loop apex the photosphere |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
a)- c) The ratio of the periods of the fundamental mode to the first harmonic,
P1/P2, as a function of density stratification parameter,
|
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Showing the value of
P1/P2 as a function of the loop shape parameter, |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Shown is the percentage difference between the values of
P1/P2 when comparing the minor elliptical and the circle-arc geometries for a range of
|
Open with DEXTER | |
In the text |
![]() |
Figure 6:
The calculated values of density scale height, taking into account
elliptical loop shape, using the observations of Van Doorsselaere et al. (2008) with corrections due to tube expansion ( |
Open with DEXTER | |
In the text |
![]() |
Figure 7:
a) Sketch showing the a loop with major elliptical shape.
The shape is described by the arc of an ellipse with a minor half-axis of length b and a major half-axis of length a, the height of the loop apex above the photosphere, |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.