Oscillations of coronal loops and second pulsations of solar radio emission

The dispersion properties of the sausage eigenmodes of oscillations in a thin magnetic flux tube are numerically analyzed in terms of ideal magnetohydrodynamics (MHD). The period of the modes accompanied by the emission of MHD waves into the surrounding medium, which leads to acoustic damping of oscillations, is determined by the radius of the tube, not by its length. The dissipation of the sausage oscillations in comparatively high (?0.7R _⊙) and tenuous (?6 × 10⁸ cm^?3) coronal loops is considered. Their Q factor has bound found to be determined by the acoustic damping mechanism. The ratio of the plasma densities outside and inside the loop and the characteristic height of the emission source have been estimated by assuming the quasi-periodic pulsations of meter-wavelength radio emission to be related to the sausage oscillations.     

Coronal Magnetic Flux Rope Equilibria and Magnetic Helicity

1 INTRODUCTIONObservations show that the magnetic helicity of solar magnetic structures has a predominantsign in each hemisphere of the Sun, positive in the southern hemisphere and negative in thenorthern, regardless of the solar cycle (Rust, 1994). The magnetic helicity is strictly conservedin the frame of ideal MHD (WOltjer, 1958), and approximately conserved in the presence ofresistive dissipation and magnetic reconnection in a highly conductive plajsma (Taylor, 1974;Berger, 1984; H…     

Transient Coronal Sigmoids and Rotating Erupting Flux Ropes

To determine the relationship between transient coronal (soft X-ray or EUV) sigmoids and erupting flux ropes, we analyse four events in which a transient sigmoid could be associated with a filament whose apex rotates upon eruption and two further events in which the two phenomena were spatially but not temporally coincident. We find the helicity sign of the erupting field and the direction of filament rotation to be consistent with the conversion of twist into writhe under the ideal MHD constraint of helicity conservation, thus supporting our assumption of flux rope topology for the rising filament. For positive (negative) helicity the filament apex rotates clockwise (counterclockwise), consistent with the flux rope taking on a reverse (forward) S shape, which is opposite to that observed for the sigmoid. This result is incompatible with two models for sigmoid formation: one identifying sigmoids with upward arching kink-unstable flux ropes and one identifying sigmoids with a current layer between two oppositely sheared arcades. We find instead that the observations agree well with the model by Titov and Démoulin (Astron. Astrophys. 351, 707, 1999), which identifies transient sigmoids with steepened current layers below rising flux ropes.     

Dissipation of Longitudinal Oscillations in Stratified Nonisothermal Hot Coronal Loops

We investigate the damping of longitudinal (i.e., slow or acoustic) waves in nonisothermal, hot (T≥ 5.0 MK), gravitationally stratified coronal loops. Motivated by SOHO/SUMER and Yohkoh/SXT observations, and by taking into account a range of dissipative mechanisms such as thermal conduction, compressive viscosity, radiative cooling, and heating, the nonlinear governing equations of one-dimensional hydrodynamics are solved numerically for standing-wave oscillations along a magnetic field line. A semicircular shape is chosen to represent the geometry of the coronal loop. It was found that the decay time of standing waves decreases with the increase of the initial temperature, and the periods of oscillations are affected by the different initial footpoint temperatures and loop lengths studied by the numerical experiments. In general, the period of oscillation of standing waves increases and the damping time decreases when the parameter that characterises the temperature at the apex of the loop increases for a fixed footpoint temperature and loop length. A relatively simple second-order scaling polynomial between the damping time and the parameter determining the apex temperature is found. This scaling relation is proposed to be tested observationally. Because of the lack of a larger, statistically relevant number of observational studies of the damping of longitudinal (slow) standing oscillations, it can only be concluded that the numerically predicted decay times are well within the range of values inferred from Doppler shifts observed by SUMER in hot coronal loops.     

On the Excitation of Leaky Modes in Cylindrical Loops

The role of leaky waves in the coronal loop oscillations observed by TRACE is not yet clearly understood. In this work, the excitation of fast waves in solar coronal loops modelled as dense plasma cylindrical tubes in a uniform straight magnetic field is investigated. We study the trapped and especially leaky modes (whose energy escapes from the tube) that result from an initial disturbance by solving the time-dependent problem numerically. We find that the stationary state of the tube motion is given by the trapped normal modes. By contrast, the transient behaviour between the initial and the stationary phase is dominated by wave leakage. The so-called trig leaky modes are clearly identified since the transient behaviour shows periods and damping times that are in agreement with the values calculated from the normal-mode analysis. Consequently, these radiating modes have physical significance. However, we have not found any evidence for the excitation of other types of modes, such as the principal leaky kink mode. J. Andries is postdoctoral Fellow of the National Fund for Scientific Research – Flanders (Belgium) (F.W.O.-Vlaanderen).     

Strong and Weak Damping of Slow MHD Standing Waves in Hot Coronal Loops

Slow magnetohydrodynamic (MHD) standing wave oscillations in hot coronal loops for both strong (i.e. τ_d/P∼ 1) and weak (i.e. τ_d/P≥ 2) damping are investigated taking account of viscosity, thermal conductivity and optically thin radiation. The individual effect of the dissipative terms is not sufficient to explain the observed damping. However, the combined effect of these dissipative terms is sufficient to explain the observed strong damping, as well as weak damping seen by SUMER. We find that, the ratio of decay time (τ_d) and period (P) of wave, i.e., τ_d/P (which defines the modes of damping, whether it is strong or weak) is density dependent. By varying density from 10⁸ to 10¹⁰ cm⁻³ at a fixed temperature in the temperature range 6 – 10 MK, observed by SUMER, we get two sets of damping: one for which τ _d/P∼ 1 corresponds to strong damping that occurs at lower density and another that occurs at higher density for which τ_d/P ≥ 2 corresponds to weak damping. Contrary to strong-damped oscillations, the effect of optically thin radiation provides some additional dissipation apart from thermal conductivity and viscosity in weak-damped oscillations. We have, therefore, derived a resultant dispersion relation including the effect of optically thin radiation. Solutions of this dispersion relation illustrate how damping time varies with physical parameters of loops in both strong and weak damping cases.     

A New Look at Mode Conversion in a Stratified Isothermal Atmosphere

Recent numerical investigations of wave propagation near coronal magnetic null points (McLaughlin and Hood: Astron. Astrophys. 459, 641, 2006) have indicated how a fast MHD wave partially converts into a slow MHD wave as the disturbance passes from a low-β plasma to a high-β plasma. This is a complex process and a clear understanding of the conversion mechanism requires the detailed investigation of a simpler model. An investigation of mode conversion in a stratified, isothermal atmosphere with a uniform, vertical magnetic field is carried out, both numerically and analytically. In contrast to previous investigations of upward-propagating waves (Zhugzhda and Dzhalilov: Astron. Astrophys. 112, 16, 1982a; Cally: Astrophys. J. 548, 473, 2001), this paper studies the downward propagation of waves from a low-β to high-β environment. A simple expression for the amplitude of the transmitted wave is compared with the numerical solution.     

A Study of Surges: II. On the Relationship between Chromospheric Surges and Coronal Mass Ejections

Liu et al. (Astrophys. J. 628, 1056, 2005a) described one surge – coronal mass ejection (CME) event showing a close relationship between solar chromospheric surge ejection and CME that had not been noted before. In this work, large Hα surges (>72 Mm, or 100 arcsec) are studied. Eight of these were associated with CMEs. According to their distinct morphological features, Hα surges can be classified into three types: jetlike, diffuse, and closed loop. It was found that all of the jetlike surges were associated with jetlike CMEs (with angular widths ≤30 degrees); the diffuse surges were all associated with wide-angle CMEs (e.g., halo); the closed-loop surges were not associated with CMEs. The exclusive relation between Hα surges and CMEs indicates difference in magnetic field configurations. The jetlike surges and related narrow CMEs propagate along coronal fields that are originally open. The unusual transverse mass motions in the diffuse surges are suggested to be due to magnetic reconnections in the corona that produce wide-angle CMEs. For the closed-loop surges, their paths are just outlining stable closed loops close to the solar surface. Thus no CMEs are associated with them.     

The heating of the solar corona

The heating of the solar corona has been a fundamental astrophysical issue for over sixty years. Over the last decade in particular, space-based solar observatories (Yohkoh, SOHO and TRACE) have revealed the complex and often subtle magnetic-field and plasma interactions throughout the solar atmosphere in unprecedented detail. It is now established that any energy release mechanism is magnetic in origin - the challenge posed is to determine what specific heat input is dominating in a given coronal feature throughout the solar cycle. This review outlines a range of possible magnetohydrodynamic (MHD) coronal heating theories, including MHD wave dissipation and MHD reconnection as well as the accumulating observational evidence for quasi-periodic oscillations and small-scale energy bursts occurring in the corona. Also, we describe current attempts to interpret plasma temperature, density and velocity diagnostics in the light of specific localised energy release. The progress in these investigations expected from future solar missions (Solar-B, STEREO, SDO and Solar Orbiter) is also assessed.Received: 6 February 2003, Published online: 14 November 2003 Correspondence to: R. W. Walsh     

Dissipation of Standing Slow Magnetoacoustic Waves in Hot Coronal Loops

The damping of standing slow waves in hot (T>6 MK) coronal loops of semicircular shape is revisited in both the linear and nonlinear regimes. Dissipation by thermal conduction, compressive viscosity, radiative cooling, and heating are examined for nonstratified and stratified loops. We find that for typical conditions of hot SUMER loops, thermal conduction increases the period of damped oscillations over the sound-crossing time, whereas the decay times are mostly shaped by compressive viscosity. Damping from optically thin radiation is negligible. We also find that thermal conduction alone results in slower damping of the density and velocity waves compared to the observations. Only when compressive viscosity is added do these waves damp out at the same rate as the observed rapidly decaying modes of hot SUMER loop oscillations, in contrast to most current work, which has pointed to thermal conduction as the dominant mechanism. We compare the linear predictions with numerical hydrodynamic calculations. Under the effects of gravity, nonlinear viscous dissipation leads to a reduction of the decay time compared to the homogeneous case. In contrast, the linear results predict that the damping rates are barely affected by gravity.     

Global Coronal Seismology

Following the observation and analysis of large-scale coronal-wave-like disturbances, we discuss the theoretical progress made in the field of global coronal seismology. Using simple mathematical techniques we determine average values for the magnetic field together with a magnetic map of the quiet Sun. The interaction between global coronal waves and coronal loops allows us to study loop oscillations in a much wider context, i.e. we connect global and local coronal oscillations.     

How to Reach Superequipartition Field Strengths in Solar Magnetic Flux Tubes

A number of independent arguments indicate that the toroidal flux system responsible for the sunspot cycle is stored at the base of the convection zone in the form of flux tubes with field strength close to 10⁵ G. Although the evidence for such strong fields is quite compelling, how such field strength can be reached is still a topic of debate. Flux expulsion by convection should lead to about the equipartition field strength, but the magnetic energy density of a 10⁵-G field is two orders of magnitude larger than the mean kinetic energy density of convective motions. Line stretching by differential rotation (i.e., the “Ω effect” in the classical mean-field dynamo approach) probably plays an important role, but arguments based on energy considerations show that it does not seem feasible that a 10⁵-G field can be produced in this way. An alternative scenario for the intensification of the toroidal flux system in the overshoot layer is related to the explosion of rising, buoyantly unstable magnetic flux tubes, which opens a complementary mechanism for magnetic-field intensification. A parallelism is pointed out with the mechanism of “convective collapse” for the intensification of photospheric magnetic flux tubes up to field strengths well above equipartition; both mechanisms, which are fundamentally thermal processes, are reviewed.     

Numerical Experiments Base on the Catastrophe Model of Solar Eruptionstwo

On the basis of the catastrophe model developed by Isenberg et al., we have used the NIRVANA code to perform the magnetohydrodynamics (MHD) numerical experiments to look into the various behaviors of the coronal magnetic configuration that includes a current-carrying flux rope for modelling the prominence levitation in the corona. These behaviors include the evolution of the equilibrium height of magnetic flux rope with the background magnetic field, the corresponding interior equilibrium of magnetic flux rope, the dynamic properties of magnetic flux rope after the system loses equilibrium, as well as the impact of the reference radius on the equilibrium height of magnetic flux rope. In our calculations, an empirical model of the coronal density distribution given by Sittler & Guhathakurta is used, and the physical dissipation is included. Our experiments show that a deviation between the simulated equilibrium height of magnetic flux rope and the theoretical result of Isenberg et al. exists, but it is not apparent, and the evolutionary features of the two results are similar. If the magnetic flux rope is initially located at the stable branch of the theoretical equilibrium curve, the magnetic flux rope will quickly reach the equilibrium position after several rounds of oscillations as a result of the self-adjustment of the system; when the system is located at the critical point it will quickly lose equilibrium and evolve to the eruptive state; the impact of the variation of reference radius on the equilibrium height of magnetic flux rope is consistent with the prediction of the theory; in the eruptive state, the kinetic properties of magnetic flux rope are consistent with the results given by the Lin-Forbes model and observation, and the fast-mode shock in front of the magnetic flux rope is observed in our experiments; furthermore, because that the dissipation is included in our numerical experiments, the energy conversion from the magnetic energy to other forms of energy is very apparent in the eruptive process.     

An Estimate of <Emphasis Type="Italic">P</Emphasis>-Mode Damping by Wave Leakage

High-cadence TRACE observations show that outward-propagating intensity disturbances are a common feature in large, quiescent coronal loops. Analysis of the frequency distribution of these modes shows peaks at both three- and five-minute periods, indicating that they may be driven by the solar surface oscillations (p modes). The energy flux contained within the coronal intensity disturbances is of the order of (1.1±0.4)×10³ ergs cm⁻² s⁻¹. A simple order-of-magnitude estimate of the damping rate of the relevant p modes allows us to put an observational constraint on the damping of p modes and shows that leakage into the overlying coronal atmosphere might be able to account for a significant fraction of p-mode damping.     

MHD equilibrium and stability properties of a bipolar current loop

A study is made of equilibrium and stability properties of a semi-toroidal current loop imbedded in a high temperature plasma. The loop carries a toroidal current density J _t and poloidal current density J _p. By explicity including the global curvature of the loop, the net Lorentz and pressure forces acting along the major radius are calculated. Requirement of equilibrium force-balance gives rise to conditions that must be satisfied by the physical parameters and geometry. On the basis of these conditions, we deduce a class of equilibrium semi-toroidal current loops satisfying c ^#X2212;1 J × B ? ▽p = 0. It is found that the averge pressure inside the loop is less than the ambient coronal pressure in equilibrium. Furthermore, this class of equilibria is shown to be stable to a number of destructive MHD modes. The theoretical results are discussed in the context of solar bipolar current loops.     

Oblique propagation and dissipation of Alfvén waves in coronal holes

We investigate the effect of viscosity and magnetic diffusivity on the oblique propagation and dissipation of Alfvén waves with respect to the normal outward direction, making use of MHD equations, density, temperature and magnetic field structure in coronal holes and underlying magnetic funnels. We find reduction in the damping length scale, group velocity and energy flux density as the propagation angle of Alfvén waves increases inside the coronal holes. For any propagation angle, the energy flux density and damping length scale also show a decrement in the source region of the solar wind (< 1.05 R_⊙) where these may be one of the primary energy sources, which can convert the inflow of the solar wind into the outflow. In the outer region (> 1.21 R_⊙), for any propagation angle, the energy flux density peaks match with the peaks of MgX 609.78 Å and 624.78 Å linewidths observed from the Coronal Diagnostic Spectrometer (CDS) on SOHO and the non-thermal velocity derived from these observations, justify the observed spectroscopic signature of the Alfvén wave dissipation.     

Coronal Flux Rope Catastrophe in Octapole Magnetic Fields

As demonstrated by many previous studies, a system consisting of an isolated coronal flux rope and a surrounding background magnetic field exhibits a catastrophic behavior. In particular, if the magnetic field of the system is force-free and axisymmetric in spherical geometry, the magnetic energy at the catastrophic point, referred to as the catastrophic energy threshold, is found to be larger than the corresponding partly or fully open field energy. This paper takes an octapole field as the background and introduces a flux rope within the central arcade of the octapole field. A relaxation method based on time-dependent ideal magnetohydrodynamic (MHD) simulations is used to find axisymmetric force-free field solutions in spherical geometry associated with the flux rope system. With respect to an increase of either the annular flux Φ_p or the axial flux Φ_ϕ of the rope, the system exhibits a catastrophic behavior as expected, and the catastrophic energy threshold is larger than that of the corresponding partly open field, in which the central arcade is opened up, but the remainder remains closed. For a given octapole field, the energy threshold depends on either Φ_p or Φ_ϕ at the catastrophic point, and it increases with increasing Φ_p or decreasing Φ_ϕ. On the other hand, the extent to which the central bipolar component of the octapole field is open also affects the energy threshold. These results differ from those for the bipolar background field case, in which the catastrophic energy threshold is almost independent of the magnetic properties of the flux rope at the catastrophic points and the extent to which the background field is open. The reason for such a difference is briefly discussed.     

Deflection of Coronal Rays by Remote CMEs: Shock Wave or Magnetic Pressure?

We analyze five events of the interaction of coronal mass ejections (CMEs) with the remote coronal rays located up to 90° away from the CME as observed by the SOHO/LASCO C2 coronagraph. Using sequences of SOHO/LASCO C2 images, we estimate the kink propagation in the coronal rays during their interaction with the corresponding CMEs ranging from 180 to 920 km s⁻¹ within the interval of radial distances from 3 R _⊙ to 6 R _⊙. We conclude that all studied events do not correspond to the expected pattern of shock wave propagation in the corona. Coronal ray deflection can be interpreted as the influence of the magnetic field of a moving flux rope within the CME. The motion of a large-scale flux rope away from the Sun creates changes in the structure of surrounding field lines, which are similar to the kink propagation along coronal rays. The retardation of the potential should be taken into account since the flux rope moves at a high speed, comparable with the Alfvén speed.     

Waves and Oscillations in the Corona – (Invited Review)

It has long been suggested on theoretical grounds that MHD waves must occur in the solar corona, and have important implications for coronal physics. An unequivocal identification of such waves has however proved elusive, though a number of events were consistent with an interpretation in terms of MHD waves. Recent detailed observations of waves in events observed by SOHO and TRACE removes that uncertainty, and raises the importance of MHD waves in the corona to a higher level. Here we review theoretical aspects of how MHD waves and oscillations may occur in a coronal medium. Detailed observations of waves and oscillations in coronal loops, plumes and prominences make feasible the development of coronal seismology, whereby parameters of the coronal plasma (notably the Alfvén speed and through this the magnetic field strength) may be determined from properties of the oscillations. MHD fast waves are refracted by regions of low Alfvén speed and slow waves are closely field-guided, making regions of dense coronal plasma (such as coronal loops and plumes) natural wave guides for MHD waves. There are analogies with sound waves in ocean layers and with elastic waves in the Earth's crust. Recent observations also indicate that coronal oscillations are damped. We consider the various ways this may be brought about, and its implications for coronal heating.     

Acoustic damping of fast kink oscillations of coronal loops

The damping of fast kink oscillations of solar coronal loops attributable to the radiation of MHD waves into the surroundings is considered in the thin-tube approximation. The oscillation damping decrement is calculated both by using a new energy method and by solving the dispersion equation for magnetic-tube eigenmodes. The two approaches are in good agreement under appropriate assumptions. The damping is negligible if MHD waves are radiated perpendicular to the magnetic field. The low Q factor of the loop oscillations in active regions found with the TRACE space telescope is associated with the generation of running waves that propagate along magnetic field lines.