Ion cyclotron waves,instabilities and solar wind heating

The effect of alpha particles on the dispersion relation of ion cyclotron waves and its influence on the heating of the solar wind plasma are investigated. The presence of alpha particles can dramatically change the dispersion relation of ion cyclotron waves, and significantly influence the way that ion cyclotron waves heat the solar wind plasma. We find that a spectrum of ion cyclotron waves affects the thermal anisotropy of the solar wind protons and other ions differently in interplanetary space: When alpha particles have a speed u _α>0.5v _A, and both protons and alpha particles have a thermal anisotropy T _⊥/T _∥>1, ion cyclotron waves heat protons in the direction perpendicular to the magnetic field, cool them in the parallel direction, and exert the opposite effect on alpha particles.     

Incompressible MHD Turbulence

The inertial range of incompressible MHD turbulence is most conveniently described in terms of counter propagating waves. Shear Alfvén waves control the cascade dynamics. Slow waves play a passive role and adopt the spectrum set by the shear Alfvén waves. Cascades composed entirely of shear Alfvén waves do not generate a significant measure of slow waves. MHD turbulence is anisotropic with energy cascading more rapidly along k _⊥ than along k _∥. Anisotropy increases with k _⊥ such that the excited modes are confined inside a cone bounded by k _∥∝ k _perp ^2/3. The opening angle of the cone, θ(k _⊥)∝ k _⊥ ^-1/3, defines the scale dependent anisotropy. MHD turbulence is generically strong in the sense that the waves which comprise it are critically damped. Nevertheless, deep inside the inertial range, turbulent fluctuations are small. Their energy density is less than that of the background field by a factor θ²(k _⊥)≪. MHD cascades are best understood geometrically. Wave packets suffer distortions as they move along magnetic field lines perturbed by counter propagating wave packets. Field lines perturbed by unidirectional waves map planes perpendicular to the local field into each other. Shear Alfvén waves are responsible for the mapping's shear and slow waves for its dilatation. The former exceeds the latter by θ^-1(k _⊥)≫ 1 which accounts for dominance of the shear Alfvén waves in controlling the cascade dynamics. This revised version was published online in July 2006 with corrections to the Cover Date.     

Excitation of whistler waves driven by an electron temperature anisotropy

A 3-D particle simulation of excitation of whistler waves driven by an electron temperature anisotropy (T > T ) is presented. Results show that whistler waves can have appreciable growth driven by the anisotropy. The maximum intensity of the excited whistler waves increases as a quadratic function of the anisotropy. Due to the presence of a threshold, one needs a relatively large electron temperature anisotropy above threshold to generate large-amplitude whistler waves. The average amplitude of turbulence in the context of whistler waves is up to as large as about 1% of the ambient magnetic field when T /T . The total energy density of the whistler turbulence is adequate for production of relativistic electrons in solar flares through stochastic acceleration.     

Thermalization of solar wind

Three kinetic equations describing the linear and non-linear wave-particle interaction for an anisotropic solar wind plasma have been developed. These equations have been solved numerically to find the variation inT /T with respect to time, whereT andT are the perpendicular and parallel temperatures with respect to the ambient magnetic field of the solar wind. For wave energy greater than a critical value (strong turbulence), non-linear wave-particle interactions are important but do not lead to thermalization. On the other hand, weak nonlinear interactions tend to increaseT /T , but make only a negligible contribution in the quantitative sense. Thus, only the linear wave-particle interaction remains as the significant contributer to the increase ofT /T .     

Electric fields in the solar atmosphere: A review

Macroscopic electric fields in the solar atmosphere have received much less attention than magnetic fields, although they must play a role of comparable importance in plasma heating, and in charged particle acceleration and transport. We review various remote sensing techniques that have been developed, whose sensitivity is now 5–10 V cm ^–1 for measurement of the electric field component transverse to the line-of-sight. Our review of the processes most likely to produce observable fields in the solar atmosphere indicates that quasi-static, macroscopic values of E (the electric field component parallel to the magnetic vector) well above this detection threshold are predicted by the discharge model of flares, by models of return currents associated with flare particle beams, and by models of neutral sheets associated with two-ribbon flares and post-flare loops. In addition, both E and E components may be detectable in time dependent electric fields associated with MHD and plasma waves, and with plasma turbulence. The emission measures and time-scales associated with these electrified plasma volumes are as highly uncertain as our present understanding of the volumes, plasma conditions and processes involved in the liberation of flare energy. Observations of electric field vector intensities, orientations, time-behaviour and spatial distribution at the presently attained electric field sensitivity levels could provide new, direct information of great interest in the electrodynamics of solar magnetic structures.     

Thermal properties of the solar wind plasma

In the solar wind, electrostatic ion cyclotron waves can be excited by electrons when the flow velocity becomes supersonic. The waves reduce the proton temperature anisotropy and heat the protons effectively. Temperature equations for T _e ,T _p and T _p are solved numerically in the region from 1 AU to the Sun, with the non-thermal proton heating rate included as a parameter. Distributions of T _e ,T _p , T _p and the proton heating rate are determined and found to be in good agreement with the proton heating rate expected from the linear growth rate of electrostatic ion cyclotron waves. The electron thermal conductivity is reduced approximately 2–3 times smaller than the usual collisional one due to the plasma wave instabilities. Effective energy exchange rates from proton-proton and electron-proton interactions are 1–10 and 10–100 times larger than the Coulomb collision rates v ppand v ep,respectively.     

Lorentz Force: A Possible Driving Force for Sunspot Rotation

Zhao and Kosovichev (Astrophys. J. 591, 446, 2003) found two opposite sub-photospheric vortical flows in the depth range of 0 – 12 Mm around a fast rotating sunspot. So far there is no theoretical model explaining such flow motions. In this paper, we try to explain this phenomenon from the point of view of magnetic flux tubes interacting with large-scale vortical motions of plasma. In the deeper zone under the photosphere, the magnetic force may be less than the nonmagnetic force of plasma. The vortical flow located there twists the flux tube and magnetic free energy is built up in the tube. In the shallower zone under the photosphere, the magnetic force may be greater than the nonmagnetic force. Thus, part of the stored magnetic free energy is released to drive the plasma to rotate in two opposite directions, e.g., in the depth ranges of 0 – 3(5) and 9 – 12 Mm. In addition, we also define a vector of nonpotential magnetic stress τ, which can be related to flare occurrence. It is calculated for the active region NOAA 10930 on 11 December 2006. We find that: i) the integral of its line-of-sight (LOS) stress successively increases around the magnetic neutral line (MNL) prior to and during the flare and decreases to a minimum after the flare; ii) the integral of its transverse stress exceeds the integral of its LOS component by one order of magnitude over the whole field of view; iii) the transverse stress first points toward the MNL, then along it, and finally it points away from it. We need other data to verify whether or not the magnetic energy is transported in the horizontal direction to the neutral line, and then partly changes into the energy in LOS direction before and during the flare.     

Excitation of high-frequency Alfvén waves by plasma outflows from coronal reconnection events

We study a kinetic excitation mechanism for high-frequency dispersive Alfvén waves in the solar corona by magnetic reconnection events. The ion-cyclotron and Cerenkov kinetic effects are important for these waves which we call the ion-cyclotron kinetic Alfvén waves (IC KAWs). The plasma outflowing from the reconnection site sets up a neutralized proton beam in the surrounding plasma, providing free energy for the excitation of waves. The dependence of the phase velocity of the IC KAW on the parallel wavenumber is different from that on the perpendicular wavenumber. The phase velocity is an increasing function of the perpendicular wavenumber and overtakes the Alfvén velocity for sufficiently large values of k . However, the phase velocity is a decreasing function of k , and sufficiently large values of k result in a phase velocity below the Alfvén velocity. As a result, the IC KAWs can undergo the Cerenkov resonance with both super- and sub-Alfvénic particles, and for the waves to be excited the outflow velocity does not need to be super-Alfvénic, as for KAWs, but the beam/Alfvén velocity ratio can span a wide range of values. High growth rates of the order of 10⁴ s^–1 are found for the values of the plasma parameters typical for the low solar corona. The waves excited by (sub-)Alfvénic beams are damped mainly due to kinetic wave-particle interactions with ions at the cyclotron resonance (ion-cyclotron damping), and with ions and electrons at the Cerenkov resonance (Landau damping). Therefore, IC KAWs can heat the plasma species of the corona in both the parallel and perpendicular direction, giving rise to an anisotropic heating of the ions. The observational consequences of the processes under study are discussed.     

Nonlinear excitation of small-scale Alfvén waves by fast waves and plasma heating in the solar atmosphere

We study a nonlinear mechanism for the excitation of kinetic Alfvén waves (KAWs) by fast magneto-acoustic waves (FWs) in the solar atmosphere. Our focus is on the excitation of KAWs that have very small wavelengths in the direction perpendicular to the background magnetic field. Because of their small perpendicular length scales, these waves are very efficient in the energy exchange with plasmas and other waves. We show that the nonlinear coupling of the energy of the finite-amplitude FWs to the small-scale KAWs can be much faster than other dissipation mechanisms for fast wave, such as electron viscous damping, Landau damping, and modulational instability. The nonlinear damping of the FWs due to decay FW = KAW + KAW places a limit on the amplitude of the magnetic field in the fast waves in the solar corona and solar-wind at the level B/B ₀∼10⁻². In turn, the nonlinearly excited small-scale KAWs undergo strong dissipation due to resistive or Landau damping and can provide coronal and solar-wind heating. The transient coronal heating observed by Yohkoh and SOHO may be produced by the kinetic Alfvén waves that are excited by parametric decay of fast waves propagating from the reconnection sites.     

Kinetic theory of MHD instabilities in a nonuniform plasma

Most of the MHD instabilities originating from the nonuniformity of a plasma excite MHD surface wave. When the excited wave has a frequency _s which corresponds to the local shear Alfvén wave resonance (_s = k v _a (x), where v _a is the Alfvén speed and k is the wave number in the direction of the magnetic field), the surface wave resonantly mode converts to the kinetic Alfvén wave, the Alfvén wave having a perpendicular wavelength comparable to the ion gyroradius and being able to propagate across the magnetic field. We discuss various linear and nonlinear effects of this kinetic Alfvén wave on the plasma including particle acceleration and heating. A specific example for the case of a MHD Kelvin-Helmholtz instability is given.     

Stokes Diagnostics of 2D MHD-Simulated Solar Magnetogranulation

The properties of solar magnetic fields on scales less than the spatial resolution of solar telescopes are studied. A synthetic infrared spectropolarimetric diagnostic based on a 2D MHD simulation of magnetoconvection is used for this. Analyzed are two time sequences of snapshots that likely represent two regions of the network fields with their immediate surroundings on the solar surface with unsigned magnetic flux densities of 300 and 140 G. In the first region from the probability density functions of the magnetic field strength it is found that the most probable field strength at log τ ₅=0 is equal to 250 G. Weak fields (B<500 G) occupy about 70% of the surface, whereas stronger fields (B>1000 G) occupy only 9.7% of the surface. The magnetic flux is −28 G and its imbalance is −0.04. In the second region, these parameters are correspondingly equal to 150 G, 93.3%, 0.3%, −40 G, and −0.10. The distribution of line-of-sight velocities on the surface of log τ ₅=−1 is estimated. The mean velocity is equal to 0.4 km s⁻¹ in the first simulated region. The average velocity in the granules is −1.2 km s⁻¹ and in the intergranules it is 2.5 km s⁻¹. In the second region, the corresponding values of the mean velocities are equal to 0, −1.8, and 1.5 km s⁻¹. In addition the asymmetry of synthetic Stokes V profiles of the Fe i 1564.8 nm line is analyzed. The mean values of the amplitude and area asymmetry do not exceed 1%. The spatially smoothed amplitude asymmetry is increased to 10% whereas the area asymmetry is only slightly varied.     

The Electron Firehose and Ordinary-Mode Instabilities in Space Plasmas

Self-generated wave fluctuations are particularly interesting in the solar wind and magnetospheric plasmas, where Coulomb collisions are rare and cannot explain the observed states of quasi-equilibrium. Linear theory predicts that firehose and ordinary-mode instabilities can develop under the same conditions, which makes it challenging to separate the role of these instabilities in conditioning the space-plasma properties. The hierarchy of these two instabilities is reconsidered here for nonstreaming plasmas with an electron-temperature anisotropy T _∥>T _⊥, where ∥ and ⊥ denote directions with respect to the local mean magnetic field. In addition to the previously reported comparative analysis, here the entire 3D wave-vector spectrum of the competing instabilities is investigated, with a focus on the oblique firehose instability and the relatively poorly known ordinary-mode instability. Results show a dominance of the oblique firehose instability with a threshold lower than the parallel firehose instability and lower than the ordinary-mode instability. For stronger anisotropies, the ordinary mode can grow faster, with maximum growth rates exceeding those of the oblique firehose instability. In contrast to previous studies that claimed a possible activity of the ordinary-mode in the low β [<?1] regimes, here it is rigorously shown that only the high β [>?1] regimes are susceptible to these instabilities.     

Drift instabilities associated with the particles in ring current and inner edge of plasma sheet

Electromagnetic waves propagating transverse to the magnetic field, containing inhomogenous and loss cone plasma, may become unstable due to the excitation of resonant proton, resonant electron and drift cyclotron instabilities. Resonant proton instability gets excited in inhomogenous plasma, irrespective of the presence of temperature anisotropy, loss cone or temperature gradient. However, the growth rate of this instability is much smaller than the other two instabilities. The maximum growth rates of resonant electron instability are enhanced with the increase of loss cone index, gradients in transverse temperature and magnetic field, and with the decrease of temperature anisotropy and gradients in density and parallel temperature. The drift cyclotron instability exists in a bounded range of wave numbers and its growth rate increases with the increase of electron temperature, density and magnetic field gradient, and with the decrease of proton temperature and temperature anisotropy. In the region of ring current for beyond plasmapause the resonant proton and resonant electron instabilities have the characterstic frequencies around 0.1Ω_p and growth rates ~10^?6Ω_p and 10^?3Ω_p, respectively. In the ring current region the drift cyclotron instability is not excited whereas in the plasma sheet region the frequency and growth rate of this instability are around Ω_p and 10^?2Ω_p, respectively. These instabilities can accelerate the ring current particles along the magnetic field lines and dump them into the auroral region.     

The Effects of Altitude-Dependent Wave Particle Interactions on the Polar Wind Plasma

The energization of charged particles, due to interaction with the ambient electromagnetic turbulence, has a significant influence on the plasma transport in space. The effect of wave-particle interactions on the outflow characteristics of polar wind plasma was investigated. The theoretical model included gravitational acceleration (g), polarization electrostatic field (Ep), and divergence of the geomagnetic field. Within the simulation region (1.7 to 10 earth radii, Re) the ions were assume to be collisionless and the electrons to obey a Boltzmann relation. Profiles of altitude-dependent diffusion coefficients [D⊥ (O+) and D⊥ (H+)] were computed from the wave spectral density (S) observed by the Plasma Wave Instrument (PWI) on board DE-1. The effects of WPI were introduced via a Monte Carlo technique, and an iterative approach was used in order to converge to self-consistent results. The main conclusions of this study were the following. As a result of perpendicular heating, the temperature anisotropy (T| /T⊥) was reduced and even reversed (T| < T⊥) at high altitudes. The O+ velocity distribution function developed a conic shape at high altitudes. The altitudes above which the WPI influences the O+ ions were lower than those for the H+ ions. The escape flux of O+ could be enhanced by more than an order of magnitude while the H+ flux remains constant. The O+ ions are heated more efficiently than the H+ ions, especially at low altitudes due to the 'pressure cooker' effect. As the ions are heated and move to higher altitudes, the ion's Larmor radius a _L may become comparable to the perpendicular wavelength λ⊥. As the ratio aL /λ⊥ becomes > 1, the heating rate becomes self-limited and the ion distribution displays toroidal features. This result is consistent with the observation of O+ toroidal distribution in the high altitude ionosphere. Finally, the large variability in the wave spectral density S was studied. This variability was found to change our results only in a quantitative manner, while our conclusions remained qualitatively unchanged. This revised version was published online in July 2006 with corrections to the Cover Date.     

The basic properties of galaxy groups from the volume-limited main galaxy sample of SDSS DR6

Using the Berlind et al. algorithm and the Davis et al. algorithm, we find that the mean velocity dispersion, virial radius, and virial mass of the group catalogs identified in the volume-limited sample are much smaller than those identified in the flux-limited sample. Our study shows that these properties of groups are heavily influenced by the relative values of the linking parameters and that there may be no values for the linking lengths that will work perfectly for every sample. In addition, we note that the luminosity distribution of member galaxies of groups identified by different algorithms is nearly the same, but member galaxies of groups identified using the linking lengths b_⊥ = 0.14, b_∥ = 0.75 have a higher proportion of blue galaxies and a lower proportion of red galaxies than member galaxies of groups identified using the linking length b = 0.2, and the early-type fraction of member galaxies of groups identified using the linking length b = 0.2 is higher than that of member galaxies of groups identified using the linking lengths b_⊥ = 0.14, b_∥ = 0.75. Published in Astrofizika, Vol. 51, No. 4, pp. 555–566 (August 2008).     

On 17 – 22 January 2005 Events in Space Weather

This contribution is a follow-up to the recent paper of Kuznetsov et al. (Contrib. Astron. Obs. Skalnaté Pleso 36, 85, 2006) on the ground level enhancement (GLE) on 20 January 2005. We focused on a study of Forbush decrease (FD) of 17 – 18 and 21 – 22 January 2005, respectively. The data from the neutron monitor at Lomnicky Štít (1 min counts) and from the Geomagnetic Observatory in Hurbanovo, both in Slovakia, were used as the basis for our investigation. The data on magnetic field and solar wind from GOES 10 and 12, SOHO-CELIAS, ACE and WIND satellites were used for better understanding of the global evolution of the event. The magnetic field is transformed to the RTN (Radial – Tangential – Normal) system where only the disturbed part of the field is compared, i.e., daily variations and a constant part are subtracted. The field reduction method is described. Our results are temporal vector diagrams of variation of all parameters at all positions from where we used the data. The amplitudes of |B| exceed 100 nT and variations during the arrival of the wavefront of CME take place simultaneously at the ground-based station and at GOES satellites. The character of the variations is as if there would be regions with the dominant electric charge of opposite signs, or electric currents with different orientations in the CME. On the basis of the values v _p and n _p and using certain assumptions we determined the mass of CME on 17 January and 21 January, respectively, of 10¹² kg. A decrease of the cosmic ray level runs suddenly (during 10 minutes), starting, however, about two hours after a sudden change of the magnetic field.     

Large kinetic power in FRII radio jets

We investigate the total kinetic powers (L _j) and ages (t _age) of powerful jets of four FR II radio sources (Cygnus A, 3C 223, 3C 284, and 3C 219) by the detail comparison of the dynamical model of expanding cocoons with observed ones. It is found that these sources have quite large kinetic powers with the ratio of L _j to the Eddington luminosity (L _Edd) resides in 0.02<L _j/L _Edd<10. Reflecting the large kinetic powers, we also find that the total energy stored in the cocoon (E _c) exceed the energy derived from the minimum energy condition (E _min): 2<E _c/E _min<160. This implies that a large amount of kinetic power is carried by invisible components such as thermal leptons (electron and positron) and/or protons.     

Oscillations In Galaxies

Oscillations in galaxies have been investigated by numerical simulations. The various models used have density distributions corresponding to that of polytrope of index n in the range 0 ≤ n ≤ 4 and their evolution has been followed for more than 70 crossing times. The kinetic energy shows regular and smooth oscillations for models with n = 0, 1 and 2 whereas in other models it shows noisy oscillation. The oscillation in kinetic energy is observed to have a period of 3 crossing time irrespective of the density and size of the galaxy. The amplitude of oscillation is seen to decrease as the central density of the galaxy increases. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

Statistical Comparison of Magnetic Clouds with Interplanetary Coronal Mass Ejections for Solar Cycle 23

We compare the number and characteristics of interplanetary coronal mass ejections (ICMEs) to those of magnetic clouds (MCs) by using in-situ solar wind plasma and magnetic field observations made at 1 AU during solar cycle 23. We found that ≈ 28% of ICMEs appear to contain MCs, since 103 magnetic clouds (MCs) occurred during 1995  – 2006, and 307 ICMEs occurred during 1996 – 2006. For the period between 1996 and 2006, 85 MCs are identified as part of ICMEs, and six MCs are not associated with ICMEs, which conflicts with the idea that MCs are usually a subset of ICMEs. It was also found that solar wind conditions inside MCs and ICMEs are usually similar, but the linear correlation between geomagnetic storm intensity (Dst _min) and relevant solar wind parameters is better for MCs than for ICMEs. The differences between average event duration (Δt) and average proton plasma β (〈β〉) are two of the major differences between MCs and ICMEs: i) the average duration of ICMEs (29.6 h) is 44% longer than for MCs (20.6 hours), and ii) the average of 〈β〉 is 0.01 for MCs and 0.24 for ICMEs. The difference between the definition of a MC and that for an ICME is one of the major reasons for these average characteristics being different (i.e., listed above as items i) and ii)), and it is the reason for the frequency of their occurrences being different.