Search

The critical ionization velocity (CIV) mechanism in a finite size cloud is studied with a series of electrostatic particle-in-cell simulations. It is observed that an initial seed ionization, produced by non-CIV mechanisms, generates a cross-field ion beam which excites a modified beam-plasma instability (MBPI) with frequency in the range of the lower hybrid frequency. The excited waves accelerate electrons along the magnetic field up to the ion drift energy that exceeds the ionization energy of the neutral atoms. The heated electrons in turn enhance the ion beam by electron-neutral impact ionization, which establishes a positive feedback loop in maintaining the CIV process. It is also found that the efficiency of the CIV mechanism depends on the finite size of the gas cloud in the following ways: (1) Along the ambient magnetic field the finite size of the cloud, L (sub parallel), restricts the growth of the fastest growing mode, with a wavelength lambda (sub m parallel), of the MBPI. The parallel electron heating at wave saturation scales approximately as (L (sub parallel)/lambda (sub m parallel)) (exp 1/2); (2) Momentum coupling between the cloud and the ambient plasma via the Alfven waves occurs as a result of the finite size of the cloud in the direction perpendicular to both the ambient magnetic field and the neutral drift. This reduces exponentially with time the relative drift between the ambient plasma and the neutrals. The timescale is inversely proportional to the Alfven velocity. (3) The transvers e charge separation field across the cloud was found to result in the modulation of the beam velocity which reduces the parallel heating of electrons and increases the transverse acceleration of electrons. (4) Some energetic electrons are lost from the cloud along the magnetic field at a rate characterized by the acoustic velocity, instead of the electron thermal velocity. The loss of energetic electrons from the cloud seems to be larger in the direction of plasma drift relative to the neutrals, where the loss rate is characterized by the neutral drift velocity. It is also shown that a factor of 4 increase in the ambient plasma density, increases the CIV ionization yield by almost 2 orders of magnitude at the end of a typical run. It is concluded that a larger ambient plasma density can result in a larger CIV yield because of (1) larger seed ion production by non-CIV mechanisms, (2) smaller Alfven velocity and hence weak momentum coupling, and (3) smaller ratio of the ion beam density to the ambient ion density, and therefore a weaker modulation of the beam velocity. The simulation results are used to interpret various chemical release experiments in space.

Observations of thermal ions in Saturn's inner magnetosphere suggest distributed local sources rather than diffusive mass loading from a source located further out. We suggest that the plasma is produced and maintained mainly by 'self-sputtering' of E ring dust. Sputtered particles are 'picked up' by the planetary magnetospheric field and accelerated to corotation energies (of the order of 8 eV/amu). The sputter yield for oxygen on ice at, for example, 120 eV is about 5, which implies that an avalanche of self-sputtering occurs. The plasma density is built up until it is balanced by local losses, presumably pitch angle scattering into the loss cone and absorption in the planet's ionosphere. The plasma density determines the distribution of dust in the E ring through plasma drag. Thus a feedback mechanism between the plasma and the E ring dust is established. The model accounts for the principal plasma observations and simultaneously the radial optical depth profile of the E ring.

A simple analytic model for the directly driven, adiabatic solar wind-magnetosphere-ionosphere coupling is derived on the basis of frontside reconnection and force balance in the tail. Magnetosphere-ionosphere coupling occurs through Alfven waves which are driven by changes of electric fields in the magnetosphere. The model provides an algorithm which calculates the directly driven part of the auroral electrojets from an observed solar wind variation. For two days in May 1979 the cross correlation between the prediction and observation of AE is greater than 0.9. This suggests that most of the geomagnetic activity during these two days was directly driven by the solar wind.

The critical ionization velocity (CIV) for a neutral barium (Ba) gas cloud moving across the static magnetic field is studied numerically using quasi-linear equations and a parameter range which is typical for the shaped-charge Ba gas release experiments in space. For consistency the charge exchange between the background oxygen ions and neutral atoms and its reverse process, as well as the excitation of the neutral Ba atoms, are included. The numerical results indicate that when the ionization rate due to CIV becomes comparable to the charge exchange rate the energy lost to the ionization and excitation collisions by the superthermal electrons exceeds the energy gain from the waves that are excited by the ion beam. This results in a CIV yield less than the yield by the charge exchange process.

A recent effort to digitize the plasma density by using the low frequency cutoff of trapped continuum radiation in the vicinity of the Jovian plasma sheet has revealed the existence of sharply defined density structures in the plasma sheet. These structures typically have a plasma density which is relatively constant but of order 50 percent greater or less than in the surrounding plasma. At the boundaries of these structures, the transitions from low to high density occur on time scales of about ten seconds, which correspond to spatial dimensions on the order of a few ion Larmor radii. The structures themselves last for intervals from less than a minute to more than five minutes, corresponding to size scales from a fraction of a Jovian radius to more than a Jovian radius, depending on the velocity of the structure relative to the spacecraft. In view of the importance of near corotation plasma flows, these structures are likely to be limited in both the longitudinal and radial dimensions and, therefore, could represent flux tubes with greatly varying plasma content. These observations are presented as among the first to directly address the theoretically proposed interchange instability.

The field line resonance expansion is studied with a magnetohydrodynamic (MHD) model applicable to earth's magnetotail. Using a Frobenius series solution, this expansion is found to be inconsistent. The inconsistency is seen by evaluating the terms in the model equations that the expansion neglects as small. Their effect is to provide a solution more singular than that determined by the expansion procedure. A tentative alternative expansion procedure, retaining these neglected terms approximately, yields a nonsingular solution, suggesting that resonant mode conversion rather than field line resonance is the appropriate phenomenon.

The present study reinvestigates the evidence for nonadiabatic transport in the quiet central plasma sheet, using AMPTE IRM data from the plasma sheet boundary layer and active times selected on the basis of large AE values. It is found that as the plasma is transported from the plasma sheet boundary layer into the central plasma sheet, both its temperature and its density (n) increase. The plasma obeys the relation p varies as n exp 4/3 for quiet times (AE is less than 100 nT) and p varies as n exp 5/3 for AE greater than 300 nT. The temperature in the quiet plasma sheet is usually less than 6 keV, and high-temperature values are more likely to be observed in what is defined as the active plasma sheet. It is suggested that the plasma sheet contains a mixture of high-entropy 'bubbles' and low-entropy 'blobs.' It is argued that these either merge or are lost from the tail before they are convected into the near-earth tail.

A recent effort to digitize the plasma density by using the low frequency cutoff of trapped continuum radiation in the vicinity of the Jovian plasma sheet has revealed the existence of sharply defined density structures in the plasma sheet. These structures typically have a plasma density which is relatively constant but of order 50 percent greater or less than in the surrounding plasma. At the boundaries of these structures, the transitions from low to high density occur on time scales of about ten seconds, which correspond to spatial dimensions on the order of a few ion Larmor radii. The structures themselves last for intervals from less than a minute to more than five minutes, corresponding to size scales from a fraction of a Jovian radius to more than a Jovian radius, depending of the velocity of the structure relative to the spacecraft. In view of the importance of near corotation plasma flows, these structures are likely to be limited in both the longitudinal and radial dimensions and, therefore, could represent flux tubes with greatly varying plasma content. These observations are presented as among the first to directly address the theoretically proposed interchange instability.

A new code is described with open boundary conditions for simulating electrostatic instabilities driven by field-aligned currents in narrow channels whose widths are smaller than the plasma skin depth. The reflection of Alfven waves at regions containing plasma waves coupling the electron and ion motion is found to cause a reduction of the field-aligned current and limits it to a critical value which is marginal for the excitation of electrostatic ion-cyclotron waves at a height of 1-2 R(e). Electron heating due to current-driven plasma turbulence causes the overall plasma density in the system to decrease. The plasma density in the region of upward current is enhanced relative to that in the region of downward current, and the current density in upward-directed current channels can be larger than in downward-directed ones.

Results are reported from a stochastic analysis of a 5-day time series of the geomagnetic AE index during an active period. The original data, the power spectrum, and the autocorrelation function are shown, and the steps in the analysis are described in detail. It is found that the autocorrelation time scale is about 50 min, giving a correlation dimension (for the construction of a time series of m-dimensional vectors) of 2.4. This result is consistent with either colored-noise or deterministic-chaos magnetosphere models, indicating the need for further investigation.

Because of the dawn-dusk electric field E(dd), plasma in the magnetotail convects from the lobe toward the central plasma sheet (CPS). In the absence of space or velocity diffusion due to plasma turbulence, convection would yield a steady state distribution function f = V exp (-2/3) g(v exp 2 V exp 2/3), where V is the flux tube volume. Starting with such a distribution function and a plasma beta which varies from beta greater than 1 in the CPS to beta much smaller than 1 in the lobe, the evolution of the ion distribution function was studied considering the combined effects of ion diffusion by kinetic Alfven waves (KAW) in the ULF frequency range (1-10 mHz) and convection due to E(dd) x B drift in the plasma sheet boundary layer (PSBL) and outer central plasma sheet (OCPS). The results show that, during the early stages after launching the KAWs, a beamlike ion distribution forms in the PSBL and at the same time the plasma density and temperature decrease in the OCPS. Following this stage, ions in the beams convect toward the CPS resulting in an increase of the plasma temperature in the OCPS.

Voyager observations have revealed that the Saturnian B ring is much more structured than previously anticipated. Some of these structures occur on radial scale lengths of less than 15 km. In the visible the ring an optical depth minimum at a small distance outside of synchronous orbit. It is suggested in this paper that the electromagnetically induced radial transport of angular momentum associated with radial transport of charged submicron-size dust particles may explain these features. This mechanism induces an instability which produces, over geological times, significant radial structuring of the ring. Mass contained in particles less than about 1 cm in size is transported away from the synchronous orbit. This produces the optical depth minimum there. A simulation, starting with a radially uniform mass and number density, yields a radial optical depth profile which is similar to what is observed. It is found that it takes about 400-800 million yr to produce the observed structure. This estimate for the age of B ring agrees with some other independent estimates.

Current collection by outgassing probes in motion relative to a magnetized plasma may be significantly affected by plasma processes that cause electron heating and cross field transport. Simulations of a neutral gas cloud moving across a static magnetic field are discussed. The authors treat a low-Beta plasma and use a 2-1/2 D electrostatic code linked with the authors' Plasma and Neutral Interaction Code (PANIC). This study emphasizes the understanding of the interface between the neutral gas cloud and the surrounding plasma where electrons are heated and can diffuse across field lines. When ionization or charge exchange collisions occur a sheath-like structure is formed at the surface of the neutral gas. In that region the crossfield component of the electric field causes the electron to E times B drift with a velocity of the order of the neutral gas velocity times the square root of the ion to electron mass ratio. In addition a diamagnetic drift of the electron occurs due to the number density and temperature inhomogeneity in the front. These drift currents excite the lower-hybrid waves with the wave k-vectors almost perpendicular to the neutral flow and magnetic field again resulting in electron heating. The thermal electron current is significantly enhanced due to this heating.

Time-series observations of the magnetotail-lobe magnetic field have been Fourier analyzed to compute the frequency-weighted energy density Pfz in the range 1-30 mHz. Pfz is generally observed in the range 0.0001-0.01 gamma-squared Hz with a mean value of 0.0012 during substorm growth phases and 0.001 in the comparison intervals. No strong correlation of Pfz is found with the auroral electrojet index in either set of intervals, but during substorm growth phases Pfz may vary by an order of magnitude over time scales of 30 min, with a tendency for higher power levels to occur later in the growth phase. Increases in Pfz precede by about 10 min localized expansive phase activity observed in individual magnetograms.

The electrostatic charge of small dust grains in a plasma in which the temperature varies in time is discussed, pointing out that secondary electron emission might introduce charge separation. If the sign of the charge on small grains is opposite to that on big ones, enhanced coagulation can occur which will affect the size distribution of grains in a plasma. Two scenarios where this process might be relevant are considered: a hot plasma environment with temperature fluctuations and a cold plasma environment with transient heating events. The importance of the enhanced coagulation is uncertain, because the plasma parameters in grain-producing environments such as a molecular cloud or a protoplanetary disk are not known. It is possible, however, that this process is the most efficient mechanism for the growth of grains in the size range of 0.1-500 microns.

The exchange of energy between the plasma mantle and the plasma sheet boundary layer (PSBL) is examined with a one-dimensional magnetotail model. The energy exchange occurs via Poynting flux generated by the localized mode conversion of a surface wave to an Alfven wave. This Poynting flux propagates through the lobe and into the PSBL where it is absorbed by two processes. The first arises from a gradient in the plasma beta causing a smooth absorption of Poynting flux. The second process results from the localized mode conversion of the decaying surface wave to an Alfven wave, causing a localized absorption of energy. A numerical solution of the linearized ideal MHD equations is obtained by assuming an adiabatic equation of state.

The dc electric field acceleration of electrons out of a thermal plasma and the evolution of the runaway tail are studied numerically, using a relativistic quasi-linear code based on the Ritz-Galerkin method and finite elements. A small field-aligned electric field is turned on at a certain time. The resulting distribution function from the runaway process is used to calculate the synchrotron emission during the evolution of the runaway tail. It is found that, during the runaway tail formation, which lasts a few tens of seconds for typical solar flare conditions, the synchrotron emission level is low, almost ot the same order as the emission from the thermal plasma, at the high-frequency end of the spectrum. However, the emission is enhanced explosively in a few microseconds by several orders of magnitude at the time the runaway tail stops growing along the magnetic field and tends toward isotropy due to the pitch-angle scattering of the fast particles. Results indicate that, in order to account for the observed synchrotron emission spectrum of a typical solar flare, the electric field acceleration phase must be accompanied or preceded by a heating phase which yields an enhanced electron temperature of about 2-15 keV in the flare region if the electric field is 0.1-0.2 times the Dreicer field and cyclotron-to-plasma frequency ratios are of order 1-2.

An expression is derived that describes the coherent Cerenkov radiated power from a group of test particles in a plasma medium moving parallel to a magnetic field. In this analysis, each particle has an arbitrary position and velocity along a field line and, as a consequence, both the spatial and temporal coherence of the radiation are considered. As an example, it is demonstrated that a monoenergetic electron beam consisting of small pulses can generate wave powers well above incoherent levels if the pulse spacing is comparable to an integer number of emission wavelengths. It is also shown that, if the beam particles have a velocity spread, Delta-V, the wave powers will decrease in time due to the reduced temporal coherence of the particle radiators, where this coherence scales as 1/Delta-V. This latter effect applies to any charged particle beam propagating in a magnetoplasma, because even an initially monoenergetic beam becomes thermalized by electrostatic wave-particle interactions reducing the radiated power.

Magnetic-field models based on pressure equilibrium in the quiet magnetotail require nonadiabatic cooling of the plasma as it convects inward or a decrease of the flux tube content. Recent in situ observations of plasma density and temperature indicate that, during quiet convection, the flux tube content may actually increase. Thus the plasma must be cooled during quiet times. The earth plasma sheet is generally significantly hotter after the expansion phase of a substorm than before the plasma sheet thinning begins and cools during the recovery phase. Heating mechanisms such as reconnection, current sheet acceleration, plasma expansion, and resonant absorption of surface waves are discussed. It seems that all mechanisms are active, albeit in different regions of the plasma sheet. Near-earth tail signatures of substorms require local heating as well as a decrease of the flux tube content. It is shown that the resonant absorption of surface waves can provide both.

The critical ionization velocity (CIV) phenomenon involves the transfer of energy from neutrals moving relative to a plasma to electrons. Three mechanisms are discussed and illustrated with numerical simulations. If the elastic electron neutral collision frequency is large (comparable to the electron cyclotron frequency) the electrons can gain enough energy from the neutrals to ionize them. Newly created ions form beams in the plasma which are unstable to lower hybrid waves propagating obliquely to the magnetic field. These waves have an electric field component parallel to the background magnetic field which can accelerate electrons. However, the efficiency of this heating is significantly reduced by electromagnetic effects which occur when the relative velocity exceeds the Alfven speed or when the plasma beta is large. A third mechanism involves the E x B drift of electrons in the secondary electric field due to charge separation at the front of a neutral gas cloud. This E field is opposite to the neutral gas velocity and can be much larger than the polarization field. The electron drift excites the modified two-stream instability which causes rapid electron heating and thus ionization in the front.