Topology of magnetodynamic flows: Kinetic and internal energy of a one-fluid plasma

A mathematical formalism describing the relation between the structure of a one-fluid, collisionless plasma and the topological features of the vector lines of the magnetic, velocity and electric fields is suggested and a rigorous approach to the problem of determining the local features of the plasma from the properties of the velocity field v, observed at a chosen point of space, is developed. The limits to drawing conclusions from the observational data (such as encountered in magnetospheric physics) are outlined in Section 2.Energization of a plasma during a stationary convection is an example of a process characterized by the properties of the tensor $\text{V}\text{= grad}\text{v}$. In particular, if the plasma density and the velocity and pressure components in a volume element of the plasma are known, then the time rates of change of the kinetic and internal energies are described by stretchings of the element and by its shearing in the surface defined by the vectors of magnetic field and perpendicular velocity. Criteria according to which plasma gains or loses kinetic and internal energies are derived.     

Magnetic Reconnection and Turbulent Mixing: From ISM to Clusters of Galaxies

Magnetic reconnection, or the ability of the magnetic field lines that are frozen in plasma to change their topology, is a fundamental problem of magnetohydrodynamics (MHD). Webriefly examine the problem starting with the well-known Sweet-Parker scheme, discuss effectsof tearing modes, anomalous resistivity and the concept of hyperresistivity. We show that the field stochasticity by itself provides a way toenable fast reconnection even if, at the scale of individual turbulent wiggles,the reconnection happens at the slow Sweet-Parker rate. We show that fast reconnectionallows efficient mixing of magnetic field in the direction perpendicular tothe local direction of magnetic field. While the idea of stochastic reconnection still requiresnumerical confirmation,our numerical simulations testify that mixing motions perpendicular to the local magnetic field are upto high degree hydrodynamical. This suggests that the turbulent heattransport should be similar to that in non-magnetized turbulent fluid, namely,should have a diffusion coefficient ～V _L L, whereV _L is the amplitude of the turbulent velocity and L is the scale of the turbulent motions. We present numericalsimulations which support this conclusion. The applicationof this idea to thermal conductivity in clusters of galaxies shows that thismechanism may dominate the diffusion of heat and may be efficient enoughto prevent cooling flow formation.     

On the rotation of a moonlet embedded in planetary rings

We examine the rotation of a small moonlet embedded in planetary rings caused by impacts of ring particles, using analytic calculation and numerical orbital integration for the three-body problem. Taking into account the Rayleigh distribution of particles' orbital eccentricities and inclinations, we evaluate both systematic and random components of rotation, where the former arises from an average of a large number of small impacts and the latter is contribution from large impacts. Calculations for parameter values corresponding to inner parts of Saturn's rings show that a moonlet would spin slowly in the prograde direction if most impactors are small particles whose velocity dispersion is comparable to or smaller than the moonlet's escape velocity. However, we also find that the effect of the random component can be significant, if the velocity dispersion of particles is larger and/or impacts of large particles comparable to the moonlet's size are common: in this case, both prograde and retrograde rotations can be expected. In the case of a small moonlet embedded in planetary rings of equal-sized particles, we find that the systematic component dominates the moonlet rotation when m/M?0.1 (m and M are the mass of a particle and a moonlet, respectively), while the random component is dominant when m/M?0.3. We derive the condition for the random component to dominate moonlet rotation on the basis of our results of three-body orbital integration, and confirm agreement with N-body simulation.     

Potential and inductive electric fields in the magnetosphere during auroras

During quiescent auroras the large-scale electric field is essentially irrotational. The volume formed by the plasma sheet and its extension into the auroral oval is connected to an external source by electric currents, which enter and leave the volume at different electric potentials and which supply sufficient energy to support the auroral activity. The location of the actual acceleration of particles depends on the internal distribution of electric fields and currents. One important feature is the energization of the carriers of the cross-tail current and another is the acceleration of electrons precipitated through relatively low-altitude magnetic-field-aligned potential drops.Substorm auroras depend on rapid and (especially initially) localized release of energy that can only be supplied by tapping stored magnetic energy. The energy is transmitted to the charged particle via electric inductive fields.The primary electric field due to changing electric currents is redistributed in a complicated way—but never extinguished—by polarization of charges. As a consequence, any tendency of the plasma to suppress magnetic-field-aligned components of the electric fields leads to a corresponding enhancement of the transverse component.     

Radiation Induced Dynamics of the Galactic Halo

We show results of numerical simulations of a three component plasma consisting of electrons, ions and dust with external gravitation and radiation fields. We perform simulation runs, starting from an analytic halo equilibrium, balancing pressure, gravitational, and radiative forces. Within these the equilibrium is perturbed by the radiation of a typical OB-star association. The perturbation has a total energy input of 10⁷ L _⊙ and a duration of 30 Myrs. After switching off the perturbation, the simulations are continued to further investigate the dynamics induced. We start with a self consistent one-fluid MHD model without background magnetic field and show for an asymmetric case that the system approaches a new equilibrium after switching on the perturbation. Later it relaxes into the starting configuration again, when the additional radiation is turned off. We then show, first by including a disk-parallel magnetic field and then by redoing the simulations with a full three-fluid code, the influence of magnetic fields and species separation on the plasma dynamics. With our computations we demonstrate that these features can be important for the explanation of the structures of galactic halos and large scale mass flows. This revised version was published online in July 2006 with corrections to the Cover Date.     

Production of pulses and interpulses in pulsars

Pulsars accelerate the charged particles moving along their magnetic field lines due to their rapidly spinning motion. Particles gain maximum energy from pulsars within the light cylinder when they are moving along the field lines perpendicular to the rotation velocity. In pulsars with non-aligned rotation and magnetic axes, the production of two intense and sharp pulses (main pulse and interpulse) separated by 180° longitude occur at the two regions near the light cylinder where the rotation velocity is perpendicular to the magnetic field. Since the radiating particles move radially along the relativistically compressed magnetic field lines, the observer in the stationary frame receives beamed and transversely compressed radiation pulse. Near the light cylinder position angle varies smoothly during pulsar rotation in a way as Radhakrishnan and Cook (1969) expect its variation near the magnetic pole, as the field lines experience relativistic compression in the direction of rotation. The motion of two charge species along the field lines produce orthogonal modes at each pulse longitude.     

Hydromagnetic waves driven by velocity shear instability in the magnetospheric boundary layer

The bulk flow of the solar wind plasma in the flank-side of the magnetospheric boundary layer, where the magnetic field lines are closed, has a component transverse to the ambient field. There is quite a strong velocity shear. The theoretical model ignores inhomogeneities in the ambient field and the mass density which occur at the magnetopause on about the same length scale as that of the velocity shear.Consideration is restricted to hydromagnetic waves which have a k-vector nearly normal to the B_o-V_o plane, i.e., approximately the magnetopause surface (k_x >k_z ～ k_y′k_xL_B > 1 and L_B = 0.1 ～ 1.0 R_E where L_B is a characteristic length of the boundary layer). It is found that a long-period (T ? 40 sec) hydromagnetic wave [the Alfvén-like wave ($\text{Ω}{}_{\mathrm{<mtext>A</mtext>}}$)] driven by velocity shear instability can be excited in the shear plasma. It is also found that the group velocity of the HM-wave is directed almost along the magnetic field line and that the magnetic variance in the shear plasma tends to be parallel to the B_o-V_o plane. The velocity shear instability in the magnetospheric boundary layer is judged to be a likely source of long-period magnetic pulsations.     

Analysis of HST ultraviolet spectra for T Tauri stars: RU Lup

The ultraviolet spectra of the star RU Lup obtained with the Hubble Space Telescope are analyzed. Emission lines are identified. The presence of absorption components with a nearly zero residual intensity in the Mg II resonance doublet lines is indicative of mass outflow with a velocity V _∞?300 km s^?1. These lines also exhibit a broad (?1400 km s^?1 at the base) component originating in the star itself. The profiles of the (optically thin) Si II] and Si III]1892 Å lines for the first time unequivocally prove that these lines originate in an accretion shock wave rather than in the chromosphere, with the gas infall velocity being V ₀?400 km s^?1. The intensity ratio of the C IV 1550 Å and Si IV 1400 Å resonance doublet components was found to be close to unity, suggesting a high accreted-gas density, logN ₀>12.5. Molecular H₂ Lyman lines formed in the stellar wind were detected. The H I Lα luminosity of RU Lup was found from their intensities to exceed 10% of L _bol. Radiation pressure in the Lα line on atomic hydrogen may play a significant role in the initial acceleration of stellar-wind matter, but the effect of Lα emission on the dynamics of molecular gas is negligible.     

Compound diffusion of energetic particles: a Kappa model for the parallel distribution function

One can assume that energetic particles follow magnetic field lines while they propagate through a magnetized plasma. The latter scenario is usually described by the so-called field line random walk limit. This limit, however, is only valid if parallel diffusion is suppressed. As soon as the latter effect is taken into account, perpendicular transport becomes subdiffusive. This physical scenario is usually called compound diffusion or compound subdiffusion and can be described by a Chapman-Kolmogorov equation. In the latter equation the parallel distribution function is an essential ingredient. In the present paper we replace the standard Gaussian model by a Kappa distribution to compute distribution functions and mean square displacements across the field.     

Effect of general loss-cone distribution function on the shear-driven electrostatic ion-cyclotron instability

The shear-driven electrostatic ion-cyclotron instability (EICI) is studied using the loss-cone distribution function by particle aspect analysis. The effect of the loss-cone distribution on the dispersion relation and growth rate of weak shear-driven EICI is studied. The whole plasma is considered to consist of resonant and non-resonant particles. The wave is assumed to propagate obliquely to the static magnetic field. It is found that the frequency of the EICI is Doppler shifted due to the transverse inhomogeneous flow in the direction of the magnetic field. It is also found that for anisotropic plasma the critical velocity shear needed to excite EICI depends upon the loss-cone distribution index (J). With the increasing values the loss-cone distribution indices (J), the critical value of normalized velocity shear needed to generate EICI in anisotropic plasma, decreases and is of the order of the weak shear. The loss-cone distribution acts as a source of free energy and generates the weak shear-driven EICI at longer perpendicular perturbations. It also lowers the transverse and parallel energy of the resonant ions. The study may explain the frequently observed EICI in the auroral acceleration region.     

The ‘Roche-Limit’ of ionospheric plasma and the formation of the plasmapause

The plasmapause position is determined by the innermost equipotential surface which is tangent to the ‘Roche-Limit’ surface of the ionospheric plasma filling the magnetosphere. When the thermal particles corotate with the Earth's angular velocity, the ‘Roche-Limit’ equatorial distance is L_c=5.78 [R_E]. When the angular convection velocity is evaluated from the quiet time electric field distribution E3 of McIIwain (1972), L_c depends on the local time. Its minimum value is then L_C=4.5Near 2400 LT, and the plasmapause shape and position satisfactorily fit the observations. The diffusive equilibrium dnesity distribution appropriated inside the plasmasphere, becomes convectively unstable beyond L = L_c, where the collisions type of model satisfactorily represents the observations. In the intermediate region between the plasmapause and the last closed magnetic field line, contimues ionization fluxes are expected to flow out of the midlatitude ionosphere     

Electron heating resulting from interhemispherical transport of photoelectrons

Ambient electron heating rates along several magnetic field lines have been determined for subsequent studies of electron and ion temperatures. Use is made of the modified diffusion method for computing the heating of the ambient plasma, and the escape fluxes from both hemispheres are coupled by self-consistent upper boundary conditions supplied by interhemispheric fluxes degraded in energy along the magnetic field tubes. Heating rates and fluxes are presented for several low L-shells appropriate for noon solstice conditions when both hemispheres are illuminated. The opacity of the field tubes as a function of L is expected to go through a minimum due to the transition from large collective effects of coulomb small angle scattering and energy loss for high L-shells, to a domination by neutral scattering all along the field lines of low L-shells.     

Effect of parallel electric field on electrostatic ion-cyclotron instability in anisotropic plasma in the presence of ion beam and general distribution function—Particle aspect analysis

Dispersion relation, resonant energy transferred, growth rate and marginal instability criteria for the electrostatic ion-cyclotron wave with general loss-cone distribution in low-β anisotropic, homogeneous plasma in the auroral acceleration region are discussed by investigating the trajectories of the charged particles. Effects of the parallel electric field, ion beam velocity, steepness of the loss-cone distribution and temperature anisotropy on resonant energy transferred and growth rate of the instability are discussed. It is found that the effect of the parallel electric field is to stabilize the wave and enhance the transverse acceleration of ions whereas the effect of steepness of loss-cone, ion beam velocity and the temperature anisotropy is to enhance the growth rate and decrease the transverse acceleration of ions. The steepness of the loss-cone also introduces a peak in the growth rate which shifts towards the lower side of the perpendicular wave number with the increasing steepness of the loss-cone.     

Criticism of reconnection models of the magnetosphere

Reconnection involves singular lines called X-lines on the day and night sides of the magnetosphere, and the reconnection rate is proportional to the component of the electric field along the X-line. Although there is some indirect support for this model, nevertheless direct support is totally lacking. However, there are two distinct pieces of clearly contradictory observational evidence on the dayside. First is the failure to account for the implied energy dissipation by the magnetopause current, over 10¹¹ W, which should be easily observable as heating or enhanced flow of the plasma near the magnetopause. In marked contrast to this prediction, HEOS-2 satellite data reveal a plasma with decreased energy density and reduced flow. Second, the boundary of closed magnetic field lines is in the wrong location. In the reconnection process the plasma outflow would cut across open field lines toward higher latitudes; there should be a band of open field lines equatorward of the cleft. Observations of trapped energetic particles indicate closed field lines within the entry layer and cleft. Either one of these pieces of evidence is sufficient by itself to require drastic revision, even rejection, of the reconnection model. There is also contradictory evidence on the night side. The last closed field line capable of trapping energetic particles is poleward of auroral arcs. The implication is that the X-line is at the distant magnetopause, and not in the plasma sheet. Consequently, even if the reconnection process were operative at the nightside X-line, it would be isolated from steady state plasma sheet and auroral processes. On the other hand, substorm phenomena, in which stored magnetic energy is converted into particle kinetic energy, necessarily involve an induced electric field; that is excluded in theories of the reconnection process in which it is assumed that curl E = 0. Nevertheless, the observed easy access of energetic solar flare particles to the polar caps, and especially the preservation of interplanetary anisotropies as differences between the two polar caps, argues strongly for an open magnetosphere, with interconnection between geomagnetic and inter-planetary magnetic field lines. It is suggested that the resolution of this apparent paradox involves electric fields parallel to the magnetic field lines somewhere on the dawn and dusk sides of the magnetosphere, with an equipotential dayside magnetopause.     

Stability of Electrostatic Electron Cyclotron Waves in a Multi-Ion Plasma

We have studied the stability of the electrostatic electron cyclotron wave in a plasma composed of hydrogen, oxygen and electrons. To conform to satellite observations in the low latitude boundary layer we model both the ionic components as drifting perpendicular to the magnetic field. Expressions for the frequency and the growth rate of the wave have been derived. We find that the plasma can support electron cyclotron waves with a frequency slightly greater than the electron cyclotron frequency ω _ce ; these waves can be driven unstable when the drift velocities of both the ions are greater than the phase velocity of the wave. We thus introduce another source of instability for these waves namely multiple ion beams drifting perpendicular to the magnetic field.     

Dissipation of MHD waves in an inhomogeneous streaming plasma

The investigation of instabilities adopting the point of view of inhomogeneous mass flow, physically corresponds to consideration of stability of the perturbations whose wavelengths in the direction of plasma inhomogeneities are much larger than the characteristic plasma scale length. The dissipation of hydromagnetic-waves and instabilities takes place due to the inhomogeneous plasma flow. Both the velocity and plasma density vary in the direction perpendicular to the magnetic field. It is found that the Alfvén wave branch and magnetosonic branch may be driven unstable by the velocity shear. Instability, oscillatory modes, marginal instability and overstability are worked out.     

Gravitational instability of a thermally-conducting plasma flowing through a porous medium in the presence of suspended particles

Magnetogravitational instability of a thermally-conducting, rotating plasma flowing through a porous medium with finite conductivity and finite Larmor radius in the presence of suspended particles has been investigated. The wave propagation has been considered for both parallel and perpendicular axes of rotation. Magnetic field is being taken in the vertical direction. A general dispersion relation has been derived through relevant linearized perturbation equations. It has been observed that the condition of instability is determined by the Jeans's criterion in its modifed form. Thermal conductivity replaces the adiabatic velocity of sound by the isothermal one. Rotation decreases the Larmor radius. Porosity decreases the Alfvén velocity. In case of a viscous medium the effects of FLR, rotation, and suspended particles are not observed in the Jeans's condition, for transverse propagation for rotational axis parallel to the magnetic field. The effects of rotation and FLR are decreased by the porosity and the suspended particles. Finite conductivity removes the Alfvén velocity from Jeans's condition.     

Slow convection of a magnetized plasma and the earth plasma sheet

Stationary convection of an isotropic, infinitely conducting plasma in a magnetic field with non-trivial geometry is discussed under the assumption that the inertial term in the equation of motion may be ignored. The energy gained or lost by a volume element of plasma per unit time does not vary along the field-lines. Simple relations between the components of the current density, depending on the field-line geometry, exist. Similar relations hold for the components of the plasma velocity.The theoretical analysis is applied to the geomagnetically-quiet plasma sheet and a qualitative physical picture of the sheet is suggested. The observed structure of the sheet is compatible with Axford-Hines type of convection perhaps combined with a low-speed flow from a distant neutral point. The magnetic-field-aligned currents are driven by the deformations of the closed field-lines which are enforced by the solar wind.     

Viscous solutions of the two-fluid solar wind equations

The two-fluid equations describing the ideal, steady, viscous solar wind are examined, and supersonic solutions are sought in which the electron and ion temperatures vary as inverse powers of the radial distance from the Sun. Just two solutions are found, and these are analogous to those found by Whang et al. (1966) and Dahlberg (1970) in one-fluid theory.