Global Ambipolar Potentials and Electric Fields at Mars Inferred From MAVEN Observations
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The parallel electric field in the solar wind is much smaller than the V x B motional electric field, yet in the proper dimensions units it is very 'large'. At the orbit of earth it is within a few percent of being at the Dreicer limit. This 'large' electric field is required for quasi-neutrality; it will be shown to have interesting consequences for the electron velocity distribution function and the description of transport of heat. Interestingly, a similar dimensionless situation also occurs at the base of the transition region, while below the transition region the dimensionless electric field is very weak. These facts suggest a new way to look at the thermal-suprathermal dichotomy in velocity distributions as the response of a plasma where charge neutrality requires such large dimensionless electric fields.
This paper investigates the properties of a one-dimensional fluid model of plasma convection in the equatorial F region ionosphere. The model equations are similar in form to Burgers equation except for additional higher-order spatial derivatives. Like Burgers equation, solution to the model have the form of propagating, shocklike structures. Numerical simulations of the model closely resemble the steepened structures observed by sounding rocket plasma density probes within equatorial spread F. Simulated denstiy power spectra, like the spectra computed from in situ data, seem to possess power law forms with a break at wavelengths of about 100 m. The precise wavenumber of the spectral break is determined by the ambipolar diffusion coefficient. The model predicts that electric field fluctuations perpendicular to the direction of plasma steepening should be proportional to the plasma density fluctuations. Electric field fluctuations parallel to the steepening will be due primarily to the ambipolar field and have a Boltzmann relationship with density (square of the absolute value of delta E) approximately equal to (K(exp 2))(square of the absolute value of (delta n/n)). At wavelengths less than about 300 m, the ambipolar field should be the dominant component of the total field intensity.
The dependence of ionospheric O+ escape flux on electromagnetic energy flux and electron precipitation into the ionosphere is derived for a hypothetical ambipolar pick-up process, powered the relative motion of plasmas and neutral upper atmosphere, and by electron precipitation, at heights where the ions are magnetized but influenced by photo-ionization, collisions with gas atoms, ambipolar and centrifugal acceleration. Ion pick-up by the convection electric field produces "ring-beam" or toroidal velocity distributions, as inferred from direct plasma measurements, from observations of the associated waves, and from the spectra of incoherent radar echoes. Ring-beams are unstable to plasma wave growth, resulting in rapid relaxation via transverse velocity diffusion, into transversely accelerated ion populations. Ion escape is substantially facilitated by the ambipolar potential, but is only weakly affected by centrifugal acceleration. If, as cited simulations suggest, ion ring beams relax into non-thermal velocity distributions with characteristic speed equal to the local ion-neutral flow speed, a generalized "Jeans escape" calculation shows that the escape flux of ionospheric O+ increases with Poynting flux and with precipitating electron density in rough agreement with observations.
NASA’s Endurance rocket mission resolved a subtle energy shift in photoelectrons escaping from Earth’s polar caps, from which the first successful measurement of Earth’s ambipolar electrostatic field was made. Preflight models predicted that coulomb collisions with thermal electrons could add an extra energy shift of up to 1 eV. However, no such additional shift was reported in the initial analysis. Here, we identify a previously unreported effect related to the spacecraft wake, which introduced a systematic +0.25eV energy shift that countered and masked the shift from coulomb collisions. After applying an in-flight calibration for this wake effect, we find a total photopeak energy shift consistent with the combined influence of the ambipolar field and the subtler -0.37eV coulomb shift, predicted by the fully kinetic STET model. Further analysis of these recalibrated measurements yields an ambipolar field estimate in close agreement with the initial analysis.
New calculations of boundary-layer electron concentration profiles for entry of a blunt-nosed slender body into the earth's atmosphere are compared with previous calculations in which ambipolar diffusion was neglected. The old and new results agree in those flight regimes where ambipolar diffusion is unimportant, but large differences are noted in both peak electron concentration and profile shape at the higher altitudes, where diffusion effects are greatest. The new results are also compared with flight-measured profiles and with calculated profiles for a viscous-shock-layer theory which was recently reported in the literature. The boundary-layer results and the data agree in most respects. Differences which occur between predicted results and the data in the outer parts of the profile are discussed in terms of the effects of aerodynamic heating of the probes.
A theory is developed for quasi-neutral and nonneutral ambipolar diffusion in a weakly ionized multiconstituent plasma. Magnetic field effects are disregarded. The effects of negative ions and various positive and negative ion diffusion coefficients are identified. It is found that the presence of negative ions leads to enhanced electron diffusion and a retarded diffusion of positive ions. In some cases electrons may be forcibly propelled from the area of perturbations in ionic concentrations. The ambipolar electric field is observed to cause negative ions to flow in a countergradient fashion where the gradient is small. Electrons tend toward a more homogeneous spatial distribution in the nonneutral case than in the quasi-neutral case.
A theory is developed for the ambipolar diffusion, drift, distortion, and production of ionization inhomogeneities within weakly ionized plasma in the presence of the geomagnetic field and motions of the neutral gas. The theory developed includes spatially varying velocity fields of the neutral gas, with a special emphasis on investigating the length scale restrictions necessary for the plasma to be treated as spatially homogeneous. Several effects of the magnetic field on ionization inhomogeneities are discussed. These include the anisotropy of ambipolar diffusion, the distortion and drift due to a gravitational field, the production of irregularities due to the turbulent acceleration field, the drift relative to the neutral gas, and the distortion of ionization inhomogeneities due to the effects of a uniform wind. Pertinent equations are derived for a multispecies plasma containing positive or negative ions, and for the case of electrons and one species of positive ions.
The combined magnetic braking-ambipolar diffusion problem in weakly ionized, rigidly rotating disks is studied. An analytical solution is presented for a disk whose angular velocity and magnetic yield vectors are aligned with the symmetry axis, illustrating the effects of the relative azimuthal drift of neutrals and ions. The effects of radial drift are added, commenting on the ratio of the characteristic ambipolar diffusion and magnetic braking time scales in high-mass and low-mass disks. A numerical calculation is used to show the combined action of these two processes.
A time-independent kinetic plasma model was developed on the basis of the Li et al. (1988) semikinetic plasma model and was used to study the interaction of the Saturnian ionosphere and ring plasma. The model includes the gravitational magnetic mirror and centripetal and ambipolar electric forces, and the effect of the mixing of two plasma populations. The results obtained indicate that the density, temperature, and composition of plasma near the rings changing in the direction from the inner C ring to the outer A ring, due to the fact that the predominant source of plasma changes from the ionosphere to the rings. The model results also suggest that the outflow of hydrogen from the ionosphere to the rings may be shut off for field lines passing through the outer B and A ring, due to the ambipolar electric field set up by the warm ring plasma trapped near the ring plane by the centipetal force. In these regions, there will be a net flux of O(+) ions from the rings to the ionosphere.
The formation and evolution of polarization charge layers on cylindrical plasma streams moving in vacuum are investigated using analytic theory and 2D electrostatic particle-in-cell computer simulations. It is shown that the behavior of the electron charge layer goes through three stages. An early time expansion is driven by electrostatic repulsion of electrons in the charge layer. At the intermediate stage, the simulations show that the electron-charge-layer expansion is halted by the positively charged plasma stream. Electrons close to the stream are pulled back to the stream and a second electron expansion follows in time. At the late stage, the expansion of the ion charge layer along the magnetic field lines accompanies the electron expansion to form an ambipolar expansion. It is found that the velocities of these electron-ion expansions greatly exceed the velocities of ambipolar expansions which are driven by plasma temperatures.
A Monte Carlo code is used to explore the full effects of the convection electric field on distributions of upflowing O(+) ions from the cusp/cleft ionosphere. Trajectories of individual ions/neutrals are computed as they undergo multiple charge-exchange collisions. In the ion state, the trajectories are computed in realistic models of the magnetic field and the convection, corotation, and ambipolar electric fields. The effects of ion-ion collisions are included, and the trajectories are computed with and without simultaneous stochastic heating perpendicular to the magnetic field by a realistic model of broadband, low frequency waves. In the neutral state, ballistic trajectories in the gravitational field are computed. The initial conditions of the ions, in addition to ambipolar electric field and the number densities and temperatures of O(+), H(+), and electrons as a function of height in the cusp/cleft region were obtained from the results of Gombosi and Killeen (1987), who used a hydrodynamic code to simulate the time-dependent frictional-heating effects in a magnetic tube during its motion though the convection throat. The distribution of the ion fluxes as a function of height are constructed from the case histories.
The effects of self-consistently including superthermal electrons in the definition of the ambipolar electric field are investigated for the case of plasmaspheric refilling after a geomagnetic storm. By using the total electron population in the hydrodynamic equations, a method for incorporating superthermal electron parameters in the electric field and electron temperature calculation is developed. Also, the ambipolar electric field is included in the kinetic equation for the superthermal electrons through a change of variables using the total energy and the first adiabatic invariant. Calculations based on these changes are performed by coupling time-dependent models of the thermal plasma and superthermal electrons. Results from this treatment of the electric field and the self-consistent development of the solution are discussed in detail. Specifically, there is a decreased thermal electron density in the plasmasphere during the first few minutes of refilling, a slightly accelerated proton shock front, and a decreased superthermal electron flux due to the deceleration by the electric field. The timescales of plasmaspheric refilling are discussed and determined to be somewhat shorter than previously calculated for the thermal plasma and superthermal electron population due to the effects of the field-aligned potential.
Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field.
A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.
Kinetic effects due to wave-particle interactions and suprathermal electrons have been suggested in the literature as possible solar wind acceleration mechanisms. Ion cyclotron resonant heating, in particular, has been associated with some qualitative features observed in the solar wind. In terms of solar wind acceleration, however, it is interesting to compare the kinetic effects of suprathermal electrons with those due to the wave-particle interactions. The combined effects of the two acceleration mechanisms on the fast solar wind have been studied by Tam and Chang (1999a,b). In this study. we investigate the role of the suprathermal electron population in the acceleration of the solar wind. Our model follows the global kinetic evolution of the fast solar wind under the influence of ion cyclotron resonant heating, while taking into account Coulomb collisions, and the ambipolar electric field that is consistent with the particle distributions themselves. The kinetic effects due to the suprathermal electrons, which we define to be the tail of the electron distributions, can be included in the model as an option. By comparing the results with and without the inclusion of the suprathermal electron effects, we determine the relative importance of suprathermal electrons and wave-particle interactions in driving the solar wind. We find that although suprathermal electrons enhance the ambipolar electric potential in the solar wind considerably, their overall influence as an acceleration mechanism is relatively insignificant in a wave-driven solar wind.