Ion-driven electrostatic ion-cyclotron instability in the two-electron component plasma

Recently various waves and particle observations of auroral field lines have been carried out by the Swedish magnetospheric research satellite Viking. In order to study wave-particle interaction extensively, waves and subsequent ion beams are measured simultaneously. Both narrowband EIC waves and wideband EIC waves of up to several times the proton gyro-frequency have been detected. Along with the observed wave characteristics, a cool low-density electron population has been detected in ion-beam regions sometimes. Therefore, we try to find an effect on EIC wave instability by the presence of cold electrons using computer simulation.     

Analysis methods for multi-component wave measurements on board the DEMETER spacecraft

We describe analysis methods to estimate parameters of electromagnetic waves based on the multi-component measurements of the DEMETER spacecraft. Using the fact that the wave magnetic field is perpendicular to the wave vector, the wave normal direction can be estimated by different methods. We use these plane-wave estimates to interpret measurements of the observed wave emissions. For instance, we use the recently developed singular value decomposition (SVD) technique. The results of the plane-wave analysis have an advantage that they often allow a straightforward interpretation. These different methods have been successfully tested with the data of previous spacecraft. All these methods are also implemented in the analysis tools designed for the analysis of the DEMETER wave measurements.We show the first results of these analysis techniques for different types of wave emissions observed on board DEMETER. Obliquely propagating right-hand polarized electromagnetic waves at a few hundreds of Hz are usually connected with a multi-ion mode structure below the local proton cyclotron frequency and with a sharp lower cutoff of left-hand polarized waves, as well as with right-hand polarized waves tunelling below the multi-ion cross-over frequency. Electron and proton whistlers are also very frequently observed on DEMETER. An unusual narrow-band emission at 140 Hz (well below the local proton cyclotron frequency) serves us as another case for a detailed analysis. We find that these waves are right-hand polarized and obliquely propagating.Using this example case, we also present analysis methods to estimate continuous distribution of wave energy density as a function of wave vector directions. These techniques of wave distribution function (WDF) analysis need both wave and particle measurements. In the analyzed case, two different methods of WDF analysis give similar results consistent with the results of the plane-wave techniques. To identify the source region we use the backward ray-tracing method. The wave normal direction obtained by the analysis of multi-component data is used for a simulation of wave propagation from the point of measurement. By this procedure, we obtain an inverse trajectory of the wave ray. We can thus follow the ray path back to the anticipated source region which is in our case located a few degrees of latitude to the South from the spacecraft position.     

Electrostatic electron cyclotron harmonic instability near Ganymede

Jupiter’s moon—Ganymede—is the largest satellite in our solar system. Galileo spacecraft made six close flybys to explore Ganymede. More information was acquired about particle population, magnetic field and plasma waves during these encounters. In this paper, our aim is to study the generation of electrostatic electron cyclotron harmonic (ECH) emissions in the vicinity of Ganymede using the observed particle data. The calculated ECH wave’s growth rates are analyzed in the light of observations of plasma waves along the path of Galileo near Ganymede. Dispersion relation for electrostatic mode is solved to obtain the temporal growth rates. A new electron distribution function, fitted to distribution observed near Ganymede, is used in the calculations. A parametric study is performed to evaluate the effect of loss-cone angle and the ratio of plasma to gyro-frequency on growth rates. It is found that ECH waves growth rates generally decrease as the loss-cone angle is increased. However, the ratio plasma to gyro-frequency has almost no effect on the growth rates. These parameters vary considerably along the Galileo trajectory near Ganymede. This is the first study which relates the occurrence of ECH waves with the particle and magnetic field data in the vicinity of Ganymede. The study of ECH wave growth rate near Ganymede is important for the calculation of pitch angle scattering rates of low-energy electrons and their subsequent precipitation into the thin atmosphere of Ganymede producing ultraviolet emissions. Results of the present study may also be relevant for the upcoming JUNO and JUICE missions to Jupiter.     

Wave packet propagation in an amplifying medium and its application to the dispersion characteristics and to the generation mechanisms of Pc 1 events

The formulae which give the propagation characteristics of a wave packet in a dispersive and amplifying medium, are established. Application is made to the propagation of Pc 1 elements through a magnetosphere constituted of a cold plasma and a high energy proton population. It is shown that the spectral shape, in a frequency-time coordinate system, of the Pc 1 elements is related to two terms : v = d²ω/dk², which represents the variation of the group velocity with frequency and which depends only on the cold plasma characteristics, and μ = -d²γ/dk², in which γ is the amplification coefficient depending on the frequency and which is related to the high energy particle distribution function. When v ? μ, only the usual dispersion effects occur, but a new method is found for determining the line of force on which the micropulsations are generated, without making any assumption about the cold plasma density distribution inside the magnetosphere. It is also possible to deduce some characteristics about the high energy proton distribution. Theoretical computations are presented, which give the frequency variation of the amplification coefficient as a function of the e-folding energy and the anisotropy factor of these high energy protons. Applications are made to ~30 pearl events which are analysed in detail according to this theory. When μ ? v, other effects do appear. After a preliminary phase, the pearl elements can become parallel for a while, or even re-erect before lying again; the duration of each element gives an indication about the number of interacting particles. The conditions for the validity of the quasi-linear theory, and some other non-linear effects related with the interpretation of Pc 1 micropulsations are also discussed.     

A Kinetic Model of Coronal Heating and Acceleration by Ion-Cyclotron Waves: Preliminary Results

We present a kinetic model of the heating and acceleration of coronal protons by outward-propagating ion-cyclotron waves on open, radial magnetic flux tubes. In contrast to fluid models which typically insist on bi-Maxwellian distributions and which spread the wave energy and momentum over the entire proton population, this model follows the kinetic evolution of the collisionless proton distribution function in response to the combination of the resonant wave-particle interaction and external forces. The approximation is made that pitch-angle scattering by the waves is faster than all other processes, resulting in proton distributions which are uniform over the resonant surfaces in velocity space. We further assume, in this preliminary version, that the waves are dispersionless so these resonant surfaces are portions of spheres centered on the radial sum of the Alfvén speed and the proton bulk speed. We incorporate the fact that only those protons with radial speeds less than the bulk speed will be resonant with outward-propagating waves, so this rapid interaction acts only on the sunward half of the distribution. Despite this limitation, we find that the strong perpendicular heating of the resonant particles, coupled with the mirror force, results in substantial outward acceleration of the entire distribution. The proton distribution evolves towards an incomplete shell in velocity space, and appears vastly different from the distributions assumed in fluid models. Evidence of these distinctive distributions should be observable by instruments on Solar Probe.     

TWO-FLUID MODEL FOR HEATING OF THE SOLAR CORONA AND ACCELERATION OF THE SOLAR WIND BY HIGH-FREQUENCY ALFVÉN WAVES

A model of the solar corona and wind is developed which includes for the first time the heating and acceleration effects of high-frequency Alfvén waves in the frequency range between 1 Hz and 1 kHz. The waves are assumed to be generated by the small-scale magnetic activity in the chromospheric network. The wave dissipation near the gyro-frequency, which decreases with increasing solar distance, leads to strong coronal heating. The resulting heating function is different from other artificial heating functions used in previous model calculations. The associated thermal pressure-gradient force and wave pressure-gradient force together can accelerate the wind to high velocities, such as those observed by Helios and Ulysses. Classical Coulomb heat conduction is also considered and turns out to play a role in shaping the temperature profiles of the heated protons. The time-dependent two-fluid (electrons and protons) model equations and the time-dependent wave-spectrum equation are numerically integrated versus solar distance out to about 0.3 AU. The solutions finally converge and settle on time-stationary profiles which are discussed in detail. The model computations can be made to fit the observed density profiles of a polar coronal hole and polar plume with the sonic point occurring at 2.4 R and 3.2 R , respectively. The solar wind speeds obtained at 63 R are 740 km s^-1 and 540 km s^-1; the mass flux is 2.1 and 2.2 × 10⁸ cm^-2 s^-1 (normalized to 1 AU), respectively. The proton temperature increases from a value of 4 × 10⁵ K at the lower boundary to 2 × 10⁶ K in the corona near 2 R .     

Head-on collision of ion-acoustic solitary waves in magnetized plasmas with nonextensive electrons and positrons

This article presents the first study of the head-on collision of two ion-acoustic solitary waves (IASWs) in magnetized plasmas with nonextensive electrons and positrons using the extended Poincaré-Lighthill-Kuo (PLK) method. The effects of the ion gyro-frequency to ion plasma frequency ratio, the positron to ion number density ratio, the electrons temperature to positrons temperature ratio, and the nonextensive parameter q on the phase shifts are investigated. It is shown that these factors significantly modify the phase shifts.     

Modification of Proton Velocity Distributions by Alfvénic Turbulence in the Solar Wind

In the present paper, the proton velocity distribution function (VDF) in the solar wind is determined by numerically solving the kinetic evolution equation. We compare the results obtained when considering the effects of external forces and Coulomb collisions with those obtained by adding effects of Alfvén wave turbulence. We use Fokker–Planck diffusion terms to calculate the Alfvénic turbulence, which take into account observed turbulence spectra and kinetic effects of the finite proton gyroradius. Assuming a displaced Maxwellian for the proton VDF at the simulation boundary at 14 solar radii, we show that the turbulence leads to a fast (within several solar radii) development of the anti-sunward tail in the proton VDF. Our results provide a natural explanation for the nonthermal tails in the proton VDFs, which are often observed in-situ in the solar wind beyond 0.3 AU.     

Location and parameters of a microwave millisecond spike event

A typical microwave millisecond spike event on November 2, 1997 was observed by the radio spectrograph of National Astronomical Observatories (NAOs) at 2.6–3.8 GHz with high time and frequency resolution. This event was also recorded by Nobeyama Radio Polarimeters (NoRP) at 1–35 GHz and Radio Heliograph (NoRH) at 17 GHz. The source at 17 GHz is located in one foot-point of a small bright coronal loop of YOHKOH SXT and SOHO EIT images with strong photospheric magnetic field in SOHO MDI magnetograph. It is assumed that the electron cyclotron maser instability and gyro-resonance absorption dominate, respectively, the rising and decay phase of the spike event. For different harmonic number of gyro-frequency or magnetic field strength, a fitting program with free plasma parameters is used to minimize the difference between the observational and theoretical values of the exponential growth and decay rates for a given spike. The plasma parameters at third harmonic number are more comparable to their typical values in solar corona. Hence, it is able to provide a diagnosis for the source parameters (magnetic field, density, and temperature), the properties of radiations (wave vector and propagation angle), and the properties of non-thermal electrons (density, pitch angle, and energy). The results are also comparable with the diagnosis of the gyro-synchrotron radiation model, the frequency drift rates and a dipole magnetic field model, as well as the YOHKOH SXT and SOHO MDI data. This study is supported by the NFSC project nos. 10333030 and 10273025, and “973” program with no. G2000078403.     

Comparative study of ion cyclotron waves at Mars, Venus and Earth

Ion cyclotron waves are generated in the solar wind when it picks up freshly ionized planetary exospheric ions. These waves grow from the free energy of the highly anisotropic distribution of fresh pickup ions, and are observed in the spacecraft frame with left-handed polarization and a wave frequency near the ion’s gyrofrequency. At Mars and Venus and in the Earth’s polar cusp, the solar wind directly interacts with the planetary exospheres. Ion cyclotron waves with many similar properties are observed in these diverse plasma environments. The ion cyclotron waves at Mars indicate its hydrogen exosphere to be extensive and asymmetric in the direction of the interplanetary electric field. The production of fast neutrals plays an important role in forming an extended exosphere in the shape and size observed. At Venus, the region of exospheric proton cyclotron wave production may be restricted to the magnetosheath. The waves observed in the solar wind at Venus appear to be largely produced by the solar-wind-Venus interaction, with some waves at higher frequencies formed near the Sun and carried outward by the solar wind to Venus. These waves have some similarity to the expected properties of exospherically produced proton pickup waves but are characterized by magnetic connection to the bow shock or by a lack of correlation with local solar wind properties respectively. Any confusion of solar derived waves with exospherically derived ion pickup waves is not an issue at Mars because the solar-produced waves are generally at much higher frequencies than the local pickup waves and the solar waves should be mostly absorbed when convected to Mars distance as the proton cyclotron frequency in the plasma frame approaches the frequency of the solar-produced waves. In the Earth’s polar cusp, the wave properties of ion cyclotron waves are quite variable. Spatial gradients in the magnetic field may cause this variation as the background field changes between the regions in which the fast neutrals are produced and where they are re-ionized and picked up. While these waves were discovered early in the magnetospheric exploration, their generation was not understood until after we had observed similar waves in the exospheres of Mars and Venus.     

Ion cyclotron instability in the solar wind

An ion cyclotron instability, arising because of the relative drift between the beam and the main components of the proton distribution function in the solar wind at 1 AU, is studied. The instability is excited in a bounded range of wave numbers provided the relative drift exceeds a certain minimum value called instability threshold. For 1, the instability threshold is smaller than or equal to the threshold of magnetosonic and Alfvén instabilities. The growth rates are enhanced by increasing relative drift and ratio of beam to main proton number density and by decreasing the wave numbers.     

Generation of Type U Bursts in Solar Radio Emission

If plasma waves propagate in the direction of the plasma density decrease, their spectrum shifts to large wave numbers (to small phase velocities). This means that the spectrum of plasma waves excited by an electron beam concentrates near the distribution function (“plateau”) border, which shifts in the region of low velocities in the process of quasilinear relaxation. As the spectrum of excited plasma waves shifts in the region of large wave numbers, their frequency grows in accordance with the dispersion equation, which describes these waves. When the growth of the plasma wave frequency exceeds the decrease of the frequency owing to the regular inhomogeneity in the corona, the branch with positive frequency drift appears on the dynamic spectrum of the radio emission. Our computations allow us to estimate the density and energy of electron beams generating type U bursts.     

A theory for the solar type-I radio continuum

The type-I radio continuum may arise from the combination of two electrostatic waves, both directed nearly normal to the magnetic field. One wave, near the upper-hybrid frequency, is generated by gyroresonance with superthermal electrons and comes into equilibrium with these electrons. The other wave, at the lower-hybrid frequency, is generated by the loss-cone instability of trapped superthermal protons in those wave directions for which the lower-hybrid frequency is an exact multiple of the proton gyrofrequency. The brightness temperature of the continuum indicates both the energy of the superthermal electrons and the existance of at least a small number of superthermal protons.     

The effects of bounce resonances on wave growth rates in the magnetosphere

Expressions for the ULF wave growth rates are obtained by consistent study of the bounce and drift motions of geomagnetically trapped particles. Different cases depending on the relative values of bounce, drift and wave frequencies are examined. It is shown that the growth (damping) rate as a function of wave frequency has an oscillatory character with maxima associated with an averaged bounce frequency. The instability conditions are investigated. For wave frequencies exceeding the Larmor drift frequency, the instability appears if the derivative function with respect to the energy is positive over a sufficient period. If, however, the wave frequency is less than an effective Larmor drift frequency, the instability is possible even if this derivative is negative everywhere.     

Galileo plasma wave observations of iogenic hydrogen

The Galileo plasma wave instrument has detected intense electromagnetic wave emissions approximately centered on the second and fourth harmonics of the local proton gyrofrequency during the close equatorial flyby of Io on 7 December 1995. Their frequencies suggest these emissions are likely generated locally by an instability driven by non thermal protons. Given that this process occurs close to Io, we suggest that hydrogen-bearing compounds, escaping from Io, are broken up/ionized near this moon, thereby releasing protons. Newly-created protons are thus injected in the Jovian corotating plasma with the corotation velocity, leading to the formation of a ring in velocity space. Several electromagnetic wave–particle instabilities can be driven by a ring of newborn protons. Given that the corotating plasma is sub-Alfvénic relative to Io, the magnetosonic mode cannot be destabilized by this proton ring. The full dispersion relation is studied using the WHAMP program (Rönmark, 1982. Rep. 179. Kiruna Geophys. Inst., Kiruna, Sweden) as well as a new algorithm that allows us to fit the distribution function of newborn protons in a more realistic way. This improvement in the ring model is necessary to explain the relative narrowness of the observed spectral peaks. The measured E/B ratio is also used to identify the relevant instability and wave mode: this mode results from the coupling between the ion Bernstein and the ion cyclotron mode (IBCW). To our knowledge this mode has not yet been studied. From the instability threshold an estimate of the density of newborn protons around Io is thus given; at about 2 Io radii from the surface and 40°W longitude from the sub-Jupiter meridian, this density is found to be ≥0.5% of the local plasma density (∼4000 cm⁻³), namely ≥20 cm⁻³. Assuming a stationary pickup process and a r⁻ⁿ distribution of pickup protons within several Io radii of Io’s wake, this implies that more than 10²⁶ protons/s are created around Io. The ultimate origin of these protons is an open issue.     

Thin layered field aligned current in the solar wind

The distribution of thin layered field aligned currents (current layers) in the leading edge of the fast stream at 1 a.u. is studied by comparing solar wind plasma parameters (bulk velocity, density proton temperature and magnetic field intensity) and occurrence frequency of the current layer. Each leading edge studied is either a stream interface or an interplanetary shock. It is found that the occurrence frequency is related best to the proton temperature variations; the occurrence frequency increases with the rise of proton temperature. Possible mechanisms which cause high occurrence frequency of current layers in the leading edge are discussed.     

Effects of the magnetic field model and wave polarisation on the estimation of proton number densities in the magnetosphere using field line resonances

The cold, core plasma mass density in the Earth's magnetosphere may be deduced from the resonant behaviour of ultra-low frequency (ULF; 1–100 mHz), magnetohydrodynamic (MHD) waves. Ground-based magnetometers are the most widely used instruments for recording the signature of ULF wave activity in the magnetosphere. For a suitable model of the background magnetic field and a functional form for the variation of the proton number density with radial distance, the resonant frequencies of ULF waves provide estimates of the equatorial plasma mass density. At high latitudes, the magnetic field model becomes critical when estimating the plasma mass density from FLR data. We show that a dipole field model is generally inadequate for latitudes greater than ∼65° geomagnetic compared with models that are keyed to magnetic activity, interplanetary magnetic field and solar wind properties. Furthermore, the method often relies on the detection of the fundamental ULF resonance, which changes frequency depending on the polarisation of the oscillation. Using idealised toroidal and poloidal oscillation modes, the range of the derived densities as the ULF wave polarisation changes is of the same order as changing the density function from a constant value throughout the magnetosphere to assuming constant Alfven speed in a dipole geometry.     

The origin of pulsar glitches

It is usually assumed that pulsar glitches are caused by the large-scale unpinning of superfluid neutron vortices in the solid crust of a neutron star and that vortex motion relative to the crust is highly dissipative at low velocities, owing to the excitation of long-wavelength Kelvin waves. The force per unit length acting on a vortex as a result of Kelvin wave excitation has been calculated for a polycrystalline structure using the free-vortex Green function. An approximate upper limit for the maximum pinning force has been obtained which, for the form of structure anticipated, is many orders of magnitude too small for consistency with the observed size and frequency of glitches. The corollary is that glitches do not originate in the crust: the necessary pinning may be given by the interaction between neutron and proton vortices in the liquid core of the star.     

Comments on the proton whistler generation process

The problem of the propagation of an electromagnetic wave originating for instance in a lightning flash through the ionospheric medium is analysed in order to understand the formation at high ionospheric altitudes of the so-called proton whistler. It is shown that the accessibility of the hydrodynamic (or kinetic) proton resonance at the satellite altitude requires that a mode conversion process must take place slightly above the transition region separating the one ion (O⁺) from the two ion (O⁺ + H⁺) component plasmas. Moreover, the transformation conditions in the wave conversion region imply that the magnetic field should be (almost) perpendicular to the density gradient. Otherwise, the incident electromagnetic wave will never reach the satellite altitude in the frequency range of the proton whistler. However, some former proton whistler theories have postulated that the signal is the result of simple ionospheric propagation effects, in contradiction with the above results. These former proton whistler theories are reviewed and it is shown that the basic flaw in these theories lies in that the incident electromagnetic wave has been supposed from the beginning to have reached the high ionospheric altitudes where is located the satellite without being influenced by the lower ionospheric layers. Some various aspects, like the high variability of the wave electric to magnetic field ratio and the harmonics bands as observed by Injun are analysed in the light of the obtained results. Finally, numerical solutions of the wave dispersion relation for both the fast hydrodynamic mode (the extraordinary mode) and the slow ion kinetic mode are presented which shows that a coupling process between the two modes may take place at various frequencies between the O⁺ and the H⁺ gyrofrequencies.     

Proton cyclotron frequency phenomena in the topside ionosphere

Proton cyclotron echoes and spurs are phenomena related to the proton cyclotron frequency discovered on topside sounder ionograms from Canadian Alouette satellites. The echoes and spurs appears on the ionograms at apparent ranges which lead to a frequency close to the proton cyclotron frequency; the frequency is obtained by taking the reciprocal of the time elapsed between the transmission of the sounder pulse and the reception of the signal at the satellite. Aloutte II and ISIS and II ionograms for about sixty satellite passes were scaled to study the charateristics of these phenomena. Generally, proton cyclotron echoes and spurs occured on the ionograms at frequencies below the electron plasma frequency f_N, the echoes predominantly slightly above the electron cyclotron frequency f_H and the spurs just below f_N. They appeared most often when a harmonic of the electron cyclotron frequency nf_H(n = 1, 2, 3, 4) was approximately equal to one of the other characteristic frequencies, that is: (1) nf_H ≈ f_N, (2) nf_H ≈ f_zS, the frequency of the Z wave at the heght of the satellite, and (3) nf_H ≈ f_T, the upper hybrid resonance frequency.Proton cyclotron echoes, spurs and protein cyclotron wave patterns have many features in common in addition to their fundamental relationship with the proton cyclotron frequency. The echoes and spurs are observed most often when when nf_H overlaps one of the other characteristic frequencies, that is: nf_H ≈ f_N, nf_H ≈ f_zS, and nf_H ≈ f_T. The proton cyclotron wave pattern is observed under the first of the three conditions. It appears that the occurence of the phenomena is related to the plama conditions, the geographic location not being important in itself except that reflects different plasma conditions. Although proton cyclotron echoes and spurs were observed more often near the geomagnetic equator, consistent with the results of Matuura and Nishizaki,(8) they still observed at high latitudes even near the north geomagnetic pole.The echoes and spurs occur at frequencies below f_N, the echoes predominantly slightly above f_H and the spurs just below f_N. Generally it is easy to distinguish between the two since usually they appear separately or, if together, often an echo would terminate and a spur begin at a slightly different apparent range. But it is not always easy since sometimes it appeared that a proton cyclotron echo and a spur formed a continuous trace, suggesting that perhaps they may be different manifestations of the same phenomenon. Work is continuing in an attemp to understand the origin of proton cyclotron echoes, spurs, and proton cyclotron wave patterns.