The effect of density stratification on the resonant absorption of MHD waves in coronal loops

The standing quasi-modes of the ideal magnetohydrodynamics (MHD) in a zero-β cylindrical magnetic flux tube that undergoes a longitudinal density stratification and radial density structuring are considered. The radial structuring is assumed to be a linearly varying density profile. Using the relevant connection formulae of the resonant absorption, the dispersion relation for the fast MHD body waves is derived and solved numerically to obtain both the frequencies and damping rates of the fundamental and first-overtone,   k = 1, 2  , modes of both the kink  ( m = 1)  and fluting  ( m = 2)  waves, where k and m are the longitudinal and azimuthal mode numbers, respectively.     

The effects of density stratification on standing fast body oscillations in coronal loops and dissipation

The standing magnetohydrodynamic (MHD) quasi-linear modes in a zero-β cylindrical magnetic flux tube that undergoes a longitudinal density stratification and radial density structuring are considered. The radial structuring is assumed to have a step-like density profile. The dispersion relation for the fast MHD body waves is derived and solved numerically to obtain the frequencies of the fundamental, first-overtone and second-overtone   k = 1, 2, 3  modes of both kink  ( m = 1)  and fluting  ( m = 2)  waves, where k and m are the longitudinal and azimuthal mode numbers, respectively. Damping rates due to both viscous and resistive dissipations in the presence of the density stratification are derived and solved numerically for the first three modes of both kink and fluting waves.     

Resonantly damped oscillations of elliptically shaped stratified emerging coronal loops

The effects of both elliptical shape and stage of emergence of the coronal loop on the resonant absorption of standing kink oscillations are studied. To do so, a typical coronal loop is modeled as a zero-beta longitudinally stratified cylindrical magnetic flux tube. We developed the connection formulae for the resonant absorption of standing transversal oscillations of a coronal loop with an elliptical shape, at various stages of its emergence. Using the connection formulae, the dispersion relation is derived and solved numerically to obtain the frequencies and damping rates of the fundamental and first-overtone kink modes. Our numerical results show that both the elliptical shape and stage of emergence of the loop alter the frequencies and damping rates of the tube as well as the ratio of frequencies of the fundamental and its first-overtone modes. However, the ratio of the oscillation frequency to the damping rate is not affected by the tube shape and stage of its emergence and also is independent of the density stratification parameter.     

Torsional Alfvén waves in stratified and expanding magnetic flux tubes

The effects of both density stratification and magnetic field expansion on torsional Alfvén waves in magnetic flux tubes are studied. The frequencies, the period ratio P ₁/P ₂ of the fundamental and its first-overtone, and eigenfunctions of torsional Alfvén modes are obtained. Our numerical results show that the density stratification and magnetic field expansion have opposite effects on the oscillating properties of torsional Alfvén waves.     

Magneto-Acoustic Waves in Compressible Magnetically Twisted Flux Tubes

The oscillatory modes of a magnetically twisted compressible flux tube embedded in a compressible magnetic environment are investigated in cylindrical geometry. Solutions to the governing equations to linear wave perturbations are derived in terms of Whittaker’s functions. A general dispersion equation is obtained in terms of Kummer’s functions for the approximation of weak and uniform internal twist, which is a good initial working model for flux tubes in solar applications. The sausage, kink and fluting modes are examined by means of the derived exact dispersion equation. The solutions of this general dispersion equation are found numerically under plasma conditions representative of the solar photosphere and corona. Solutions for the phase speed of the allowed eigenmodes are obtained for a range of wavenumbers and varying magnetic twist. Our results generalise previous classical and widely applied studies of MHD waves and oscillations in magnetic loops without a magnetic twist. Potential applications to solar magneto-seismology are discussed.     

MHD surface waves in a complex (longitudinal + sheared) magnetic field

We study the dispersion characteristics of fast MHD surface waves travelling on a plasma slab immersed in a complex magnetic field consisting of a large longitudinal B _0z component and a small sheared B _0y component. The analysis shows that for typical coronal conditions both the sausage and kink waves are generally pseudo-surface waves. The tangential magnetic field, B _0y , modifies the dispersion curves, and for sufficiently large sheared fields there is a transition from pseudo-surface to pure-surface fast kink waves.On leave from Faculty of Physics, Sofia University, BG-1126 Sofia, Bulgaria.     

Linear MHD Wave Propagation in Time-Dependent Flux Tube

MHD waves and oscillations in sharply structured magnetic plasmas have been studied for static and steady systems in the thin tube approximation over many years. This work will generalize these studies by introducing a slowly varying background density in time, in order to determine the changes to the wave parameters introduced by this temporally varying equilibrium, i.e. to investigate the amplitude, frequency, and wavenumber for the kink and higher order propagating fast magnetohydrodynamic wave in the leading order approximation to the WKB approach in a zero-β plasma representing the upper solar atmosphere. To progress, the thin tube and over-dense loop approximations are used, restricting the results found here to the duration of a number of multiples of the characteristic density change timescale. Using such approximations it is shown that the amplitude of the kink wave is enhanced in a manner proportional to the square of the Alfvén speed, $V_{\mathrm{A}}^{2}$ . The frequency of the wave solution tends to the driving frequency of the system as time progresses; however, the wavenumber approaches zero after a large multiple of the characteristic density change timescale, indicating an ever increasing wavelength. For the higher order fluting modes the changes in amplitude are dependent upon the wave mode; for the m=2 mode the wave is amplified to a constant level; however, for all m≥3 the fast MHD wave is damped within a relatively small multiple of the characteristic density change timescale. Understanding MHD wave behavior in time-dependent plasmas is an important step towards a more complete model of the solar atmosphere and has a key role to play in solar magneto-seismological applications.     

The Effect of a Twisted Magnetic Field on the Nature of Kink MHD Waves

We consider a pressureless plasma in a thin magnetic-flux tube with a twisted magnetic field. We study the effect of twisted magnetic field on the nature of propagating kink waves. To do this, the restoring forces of oscillations in the linear ideal magnetohydrodynamics (MHD) were obtained. In the presence of a twisted magnetic field, the ratio of the magnetic-tension force to the gradient of the magnetic pressure increases for the mode with negative azimuthal wave number, but it decreases for the mode with positive azimuthal wave number. For the kink mode with positive azimuthal mode number, the ratio of the forces is more affected by the twisted magnetic field in dense loops. For the kink mode with negative azimuthal mode number, the perturbed magnetic pressure is negligible under some conditions. The magnetic twist increases (diminishes) the damping of the kink waves with positive (negative) azimuthal mode number due to resonant absorption. Our conclusion is that introducing a twisted magnetic field breaks the symmetry between the nature of the kink waves with positive and negative azimuthal wave number, and the wave can be a purely Alfvénic wave in the entire loop.     

Linear MHD Sausage Waves in Compressible Magnetically Twisted Flux Tubes

Oscillations of magnetic flux tubes are of great importance as they contain information about the geometry and fine structure of the flux tubes. Here we derive and analytically solve in terms of Kummer’s functions the linear governing equations of wave propagation for sausage surface and body modes (m=0) of a magnetically twisted compressible flux tube embedded in a compressible uniformly magnetized plasma environment in cylindrical geometry. A general dispersion relation is obtained for such flux tubes. Numerical solutions for the phase velocity are obtained for a wide range of wavenumbers and for varying magnetic twist. The effect of magnetic twist on the period of oscillations of sausage surface modes for different values of the wavenumber and vertical magnetic field strength is calculated for representative photospheric and coronal conditions. These results generalize and extend previous studies of MHD waves obtained for incompressible or for compressible but nontwisted flux tubes. It is found that magnetic twist may change the period of sausage surface waves of the order of a few percent when compared to counterparts in straight nontwisted flux tubes. This information will be most relevant when high-resolution observations are used for diagnostic exploration of MHD wave guides in analogy to solar-interior studies by means of global eigenoscillations in helioseismology.     

Numerical simulations of driven MHD waves in coronal loops

The solar corona, modeled by a low-, resistive plasma slab, sustains MHD wave propagations due to footpoint motions in the photosphere. Simple test cases are undertaken to verify the code. Uniform, smooth and steep density, magnetic profile and driver are considered. The numerical simulations presented here focus on the evolution and properties of the Alfvén, fast and slow waves in coronal loops. The plasma responds to the footpoint motion by kink or sausage waves depending on the amount of shear in the magnetic field. The larger twist in the magnetic field of the loop introduces more fast-wave trapping and destroys initially developed sausage-like wave modes. The transition from sausage to kink waves does not depend much on the steep or smooth profile. The slow waves develop more complex fine structures, thus accounting for several local extrema in the perturbed velocity profiles in the loop. Appearance of the remnants of the ideal singularities characteristic of ideal plasma is the prominent feature of this study. The Alfvén wave which produces remnants of the ideal x ^–1 singularity, reminiscent of Alfvén resonance at the loop edges, becomes less pronounced for larger twist. Larger shear in the magnetic field makes the development of pseudo-singularity less prominent in case of a steep profile than that in case of a smooth profile. The twist also causes heating at the edges, associated with the resonance and the phase mixing of the Alfvén and slow waves, to slowly shift to layers inside the slab corresponding to peaks in the magnetic field strength. In addition, increasing the twist leads to a higher heating rate of the loop. Remnants of the ideal log ¦x¦ singularity are observed for fast waves for larger twist. For slow waves they are absent when the plasma experiences large twist in a short time. The steep profiles do not favour the creation of pseudo-singularities as easily as in the smooth case.     

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.     

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.     

Magnetohydrodynamic Waves in an Asymmetric Magnetic Slab

Analytical models of solar atmospheric magnetic structures have been crucial for our understanding of magnetohydrodynamic (MHD) wave behaviour and in the development of the field of solar magneto-seismology. Here, an analytical approach is used to derive the dispersion relation for MHD waves in a magnetic slab of homogeneous plasma enclosed on its two sides by non-magnetic, semi-infinite plasma with different densities and temperatures. This generalises the classic magnetic slab model, which is symmetric about the slab. The dispersion relation, unlike that governing a symmetric slab, cannot be decoupled into the well-known sausage and kink modes, i.e. the modes have mixed properties. The eigenmodes of an asymmetric magnetic slab are better labelled as quasi-sausage and quasi-kink modes. Given that the solar atmosphere is highly inhomogeneous, this has implications for MHD mode identification in a range of solar structures. A parametric analysis of how the mode properties (in particular the phase speed, eigenfrequencies, and amplitudes) vary in terms of the introduced asymmetry is conducted. In particular, avoided crossings occur between quasi-sausage and quasi-kink surface modes, allowing modes to adopt different properties for different parameters in the external region.     

Non-linear damping of visco-resistive Alfvén waves in solar spicules

Interaction of Alfvén waves with plasma inhomogeneities generates phase mixing which can lead to dissipate Alfvén waves and to heat the solar plasma. Here we study the dissipation of Alfvén waves by phase mixing due to viscosity and resistivity variations with height. We also consider nonlinear magnetohydrodynamic (MHD) equations in our theoretical model. Non-linear terms of MHD equations include perturbed velocity, magnetic field, and density. To investigate the damping of Alfvén waves in a stratified atmosphere of solar spicules, we solve the non-linear MHD equations in the x–z plane. Our simulations show that the damping is enhanced due to viscosity and resistivity gradients. Moreover, energy variations is influenced due to nonlinear terms in MHD equations.     

Interaction of particles with ion cyclotron waves and magnetosonic waves. Observations from GEOS 1 and GEOS 2

Coordinated observations involving ion composition, thermal plasma, energetic particle, and ULF magnetic field data from GEOS 1 and 2 often reveal the presence of electromagnetic ion cyclotron and magnetosonic waves, which are distinguished by their respective polarization characteristics and frequency spectra. The ion cyclotron waves are identified by a magnetic field perturbation that lies in a plane perpendicular to the Earth's magnetic field B₀ and propagate along B₀. They are associated with the abundance of cold He⁺ in the presence of anisotropic pitch angle distributions of ions having energies E > 20 keV, and were observed at frequencies near the He⁺ gyrofrequency. The magnetosonic waves are characterized by a magnetic field perturbation parallel to B₀ and thus seem to be propagating perpendicular to the Earth's magnetic field. They often occur at harmonics (not always including the fundamental) at the proton gyrofrequency and are associated with phase-space-density distributions that peak at energies E ～ 5–30 keV and at a pitch angle of 90°. Such a ring-like distribution is shown to excite instability in the magnetosonic mode near harmonics of the proton gyrofrequency. Magnetosonic waves are associated in other cases with sharp spatial gradients in energetic ion intensity. Such gradients are encountered in the early afternoon sector (as a consequence of the drift shell distortion caused by the convection electric field) and could likewise constitute a source of free energy for plasma instabilities.     

Helical magnetorotational instability of Taylor‐Couette flows in the Rayleigh limit and for quasi‐Kepler rotation

The magnetorotational instability (MRI) of differential rotation under the simultaneous presence of axial and azimuthal components of the (current‐free) magnetic field is considered. For rotation with uniform specific angular momentum the MHD equations for axisymmetric perturbations are solved in a local short‐wave approximation. All the solutions are overstable for B_z · B_ϕ ≠ 0 with eigenfrequencies approaching the viscous frequency. For more flat rotation laws the results of the local approximation do not comply with the results of a global calculation of the MHD instability of Taylor‐Couette flows between rotating cylinders. – With B_ϕ and B_z of the same order the traveling‐mode solutions are also prefered for flat rotation laws such as the quasi‐Kepler rotation. For magnetic Prandtl number Pm → 0 they scale with the Reynolds number of rotation rather than with the magnetic Reynolds number (as for standard MRI) so that they can easily be realized in MHD laboratory experiments. – Regarding the nonaxisymmetric modes one finds a remarkable influence of the ratio B_ϕ/B_z only for the extrema. For B_ϕ ≫ B_z and for not too small Pm the nonaxisymmetric modes dominate the traveling axisymmetric modes. For standard MRI with B_z ≫ B_ϕ, however, the critical Reynolds numbers of the nonaxisymmetric modes exceed the values for the axisymmetric modes by many orders so that they are never prefered. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonlinear Generation of Fluting Perturbations by Kink Mode

We study the excitation of fluting perturbations in a magnetic tube by an initially imposed kink mode. We use the ideal magnetohydrodynamic (MHD) equations in the cold-plasma approximation. We also use the thin-tube approximation and scale the dependent and independent variables accordingly. Then we assume that the dimensionless amplitude of the kink mode is small and use it as an expansion parameter in the regular perturbation method. We obtain the expression for the tube boundary perturbation in the second-order approximation. This perturbation is a superposition of sausage and fluting perturbations. The amplitude of the fluting perturbation takes its maximum at the middle of the tube, and it monotonically decreases with the distance from the middle of the tube.     

Nonaxisymmetric Oscillations of Thin Twisted Magnetic Tubes

In this paper we study nonaxisymmetric oscillations of thin twisted magnetic tubes taking the density variation along the tube into account. We use the approximation of the zero-beta plasma. The magnetic field outside the tube is straight and homogeneous; however, it is twisted inside the tube. We assume that the azimuthal component of the magnetic field is proportional to the distance from the tube axis and that the tube is only weakly twisted (i.e., the ratio of the azimuthal and axial components of the magnetic field is small). Using the asymptotic analysis we show that the eigenmodes and eigenfrequencies of the kink and fluting oscillations are described by a classical Sturm – Liouville problem for a second-order ordinary differential equation. The main result is that the twist does not affect the kink mode.     

Viscous Damping of Non-Adiabatic MHD Waves in an Unbounded Solar Coronal Plasma

The damping of MHD waves in solar coronal magnetic field is studied taking into account thermal conduction and compressive viscosity as dissipative mechanisms. We consider viscous homogeneous unbounded solar coronal plasma permeated by a uniform magnetic field. A general fifth-order dispersion relation for MHD waves has been derived and solved numerically for different solar coronal regimes. The dispersion relation results three wave modes: slow, fast, and thermal modes. Damping time and damping per periods for slow- and fast-mode waves determined from dispersion relation show that the slow-mode waves are heavily damped in comparison with fast-mode waves in prominences, prominence–corona transition regions (PCTR), and corona. In PCTRs and coronal active regions, wave instabilities appear for considered heating mechanisms. For same heating mechanisms in different prominences the behavior of damping time and damping per period changes significantly from small to large wavenumbers. In all PCTRs and corona, damping time always decreases linearly with increase in wavenumber indicate sharp damping of slow- and fast-mode waves.     

The effect of non-uniform magnetic field on the slow mode oscillations

In this paper, the slow MHD mode oscillations of the coronal plasma are studied. The aim is to identify the effect of structuring (such as magnetic field, temperature, density, and pressure) on the frequencies of oscillations. We modelled the coronal medium as a low-β plasma with longitudinally density and pressure stratifications and a weakly inhomogeneous magnetic field varied slowly with height and radial directions. The linearized ideal MHD equations reduced to a single Klein–Gordon differential equation for square of oscillatory frequencies. The eigenfunctions and analytical dispersion relations are derived. The dispersion relations were solved numerically. In the case of uniform magnetic field, the previous studies verified. Our numerical results show that, the frequencies and their ratios are very sensitive functions of pressure scale height, and slightly varying functions of inhomogeneity parameter of magnetic field. By changing the magnetic field strength between the apex and footpoints of the loop about 50%, the frequencies ratio are changed about 5%. We concluded that, the pressure scale height and temperature gradient are first order effects and inhomogeneity of magnetic field is a second order effect on the slow mode oscillations.