Neutral particle collisional relaxation of suprathermal distributions in aeronomy

There are numerous situations in laboratory experiments and in atmospheric science that are characterized by distributions of energetic neutral species with extended high-energy tails. This paper provides a detailed analysis of the relaxation of such isotropic nonequilibrium distributions of neutral species. We consider a minor constituent, referred to as the test-particle dilutely dispersed in a second component that acts as a heat bath at equilibrium. A hard sphere cross section is assumed for the collisions of the test particles and the heat bath particles, and collisions between test particles are not included in the analysis. We study the approach to equilibrium with a finite difference method of solution of the Boltzmann equation. The solution of the Boltzmann equation is also presented as an expansion in the eigenfunctions of the Boltzmann collision operator. The main objective of the paper is the calculation of the energy-dependent relaxation times for the distribution function. It is anticipated that these relaxation times will not exhibit a strong energy dependence owing to the energy-independent hard sphere cross section. This relaxation behaviour is important in the characterization of nonthermal populations of energetic atoms in aeronomy. The results are compared with a similar analysis for Coulomb collisions described in the preceeding paper.     

The Fokfer-Planck equapions for magnetic monopoles and charged particles

We investigate the Fokker-Planck equations (non-relativistic case) for the interaction between magnetic monopoles and charged particles. We find that, the collision integral in this case is also logarithmically divergent, showing that the main effect is still produced by far encounters. We give, for thermal equilibrium, expressions for the three relaxation times. Compared with the well-known expressions for charged particles, these are generally amplified by a factor (light velocity/thermal velocity)², showing that, under ordinary conditions, the interaction between monopoles and charged particles is negligible but that, under certain astrophysical conditions, it is not. In the latter case, the MHD equations must be modified accordingly and the modified equations are given.     

Suprathermal particle distributions in space physics: Kappa distributions and entropy

The energization of a charged test-particle of mass m in contact with a large ensemble of charged particles of mass M at equilibrium is studied with the Fokker-Planck equation for Coulomb collisions and a quasi-linear diffusion operator for wave-particle interactions. The features of the nonequilibrium steady state velocity distribution of the test-particle system is studied as a function of the mass ratio m/M, and the relative strengths of the wave-particle interactions and Coulomb collisions. It is shown that the steady distribution function is not necessarily a Kappa distribution. The temperature of heavy minor ions given by the model is shown to vary linearly with the mass ratio as observed in the solar wind. The time evolution of the distribution function with and without the energization by wave-particle interactions is calculated and it is demonstrated that the Kullback relative entropy rather than the Tsallis nonextensive entropy rationalizes the results obtained.     

Shock waves in a dusty plasma having q-nonextensive electron velocity distribution

The dust-acoustic shock waves have been theoretically investigated using reductive perturbation technique. An unmagnetized four-component dusty plasma system consisting of nonextensive q-distributed electrons, Boltzmann distributed ions, and negatively as well as positively charged dust particles has been considered. The solution of Burgers equation in planar geometry is numerically analyzed. It has been observed that the nonextensive q-distribution of electrons has a significant role in the formation of shock waves. The relevance of our results to astrophysics as well as laboratory plasmas are briefly discussed.     

Nonlinear dynamics of Landau damped kinetic Alfvén waves

In the present paper we have studied the nonlinear dynamical equation of Landau damped kinetic Alfvén wave (KAW) to investigate the nonlinear evolution of KAW and the resulting turbulent spectra in solar wind plasmas. We have introduced a parameter g which governs the coupling between the amplitude of the pump KAW and the density perturbation. The numerical solution has been carried out to see the dependence on the parameter g in the nonlinear part of our equation. Our results reveal the formation of damped localized structures of KAW as well as steepening of the turbulent spectra by increasing g when damping is taken into account. The power spectra of magnetic field fluctuations indicate the redistribution of energy among the higher wave numbers. Each power spectrum with and without damping splits up into two different scaling ranges, Kolmogorov scaling followed by a steeper scaling. The steepening in the power spectra with Landau damping is more than without Landau damping case (for the same value of g). This type of steeper spectra has also been observed in the solar wind and is attributed to the Landau damping effects.     

Collisional balance of the meteoritic complex

Taking into account meteoroid measurements by in situ experiments, zodiacal light observations, and oblique angle hypervelocity impact studies, it is found that the observed size distributions of lunar microcraters usually do not represent the interplanetary meteoroid flux for particles with masses ?10^?10g. From the steepest observed lunar crater size distribution a “lunar flux” is derived which is up to 2 orders of magnitude higher than the interplanetary flux at the smallest particle masses. New models of the “lunar” and “interplanetary” meteoroid fluxes are presented. The spatial mass density of interplanetary meteoritic material at 1 AU is ～10^?16g/m³. A large fraction of this mass is in particles of 10^?6 to 10^?4 g. A detailed analysis of the effects of mutual collisions (i.e., destruction of meteoroids and production of fragment particles) and of radiation pressure has been performed which yielded a new picture of the balance of the meteoritic complex. It has been found that the collisional lifetime at 1 AU is shortest (～10⁴years) for meteoroids of 10^?4 to 1 g mass. For particles with masses m > 10^?5g, Poynting-Robertson lifetimes are considerably larger than collisional lifetimes. The collisional destruction rate of meteoroids with masses $\text{m ? 10}{}^{\mathrm{?3}}\text{g}$ is about 10 times larger than the rate of collisional production of fragment particles in the same mass range. About 9 tons/sec of these “meteor-sized” (m > 10^?5g) particles are lost inside 1 AU due to collisions and have to be replenished by other sources, e.g., comets. Under steady-state conditions, most of these large particles are “young”; i.e., they have not been fragmented by collisions and their initial orbits are not altered much by radiation pressure drag. Many more micrometeoroids of masses m ? 10^?5g are generated by collisions from more massive particles than are destroyed by collisions. The net collisional production rate of intermediate-sized particles 10^?10g ? m ? 10^?5g is found to be about 16 times larger at 1 AU than the Poynting-Robertson loss rate. The total Poynting-Robertson loss rate inside 1 AU is only about 0.26 tons/sec. The smallest fragment particles (m ? 10^?10g) will be largely injected into hyperbolic trajectories under the influence of radiation pressure (β meteoroids). These particles provide the most effecient loss mechanism from the meteoritic complex. When it is assumed that meteoroids fragment similarly to experimental impact studies with basalt, then it is found that interplanetary meteoroids in the mass range 10^?10g ? m ? 10^?5g cannot be in temporal balance under collisions and Poynting-Robertson drag but their spatial density is presently increasing with time.     

The vertical structure and thickness of Saturn's rings

We have considered the steady state vertical structure of Saturn's rings with regard to whether collapse to a monolayer due to collisions between particles, the end state predicted by Jeffreys (1947a), may be prevented by any of a variety of mechanisms. Given a broad distribution of particle sizes such as a typical power law n(R) = n₀R^?3, it is found that gravitational scattering of small particles by large particles maintains a true ring thickness of several times the radius of the largest particles, or many times the radius of the smallest particles. Thus the “many-particle-thick” condition which best satisfies optical observations, such as the opposition effect, may be reconciled with ongoing particle collisions. If we consider the obvious sources of energy available for such a process, we find that a ring thickness of only tens of meters may be sustained over the lifetime of the solar system. This implies a maximum particle size on the order of a few meters.     

The Fokker-Planck equation for charged particles moving in random electromagnetic fields

The Fokker-Planck equation which describes the motion of charged particles in a random electromagnetic field is derived from the Liouville equation by a new method. The size of the perturbing magnetic field, for the Fokker-Planck equation to be valid, is calculated in a regime appropriate for cosmic-ray diffusion.     

Roles of statistical mechanics in determining the density profile of mild relaxation

The r ^?4 law of cold collapse has been explained in other work. Here we try to explain the density profile of mild relaxation by statistical mechanics. In this paper we first generate many kinds of initial conditions with the same mass and energy to see whether there are other initial factors that can change the density profile of an isolated equilibrium self-gravitating system; then for a more general initial condition we discuss the role of mass and energy in determining the final density profile. Next we use our previous results obtained from statistical mechanics to fit these simulations, and find that when the masses of the particles in clumps are less than 5 % of the total mass, or the initial density is shallower than r ^?2, the whole virialized density profile (VDF) can be fitted well by our equation of state with three parameters, and some other cases can be explained by the theory with the r ^?4 law. We conclude that statistical physics may play an important role in determining the shape of VDF in the mild relaxation, mass and energy can control the values of the central density and the system’s radius, but there are still other initial configurations that can affect the VDF.     

High-energy scalarons in R gravity as a model for Dark Matter in galaxies

We show that in the framework of R² gravity and in the linearized approach it is possible to obtain spherically symmetric stationary states that can be used as a model for galaxies. Such approach could represent a solution to the Dark Matter Problem. In fact, in the model, the Ricci curvature generates a high energy term that can in principle be identified as the dark matter field making up the galaxy. The model can also help to have a better understanding on the theoretical basis of Einstein-Vlasov systems. Specifically, we discuss, in the linearized R² gravity, the solutions of a Klein-Gordon equation for the spacetime curvature. Such solutions describe high energy scalarons, a field that in the context of galactic dynamics can be interpreted like the no-light-emitting galactic component. That is, these particles can be figured out like wave-packets showing stationary solutions in the Einstein-Vlasov system. In such approximation, the energy of the particles can be thought of as the galactic dark matter component that guarantees the galaxy equilibrium. Thus, because of the high energy of such particles the coupling constant of the R²-term in the gravitational action comes to be very small with respect to the linear term R. In this way, the deviation from standard General Relativity is very weak, and in principle the theory could pass the Solar System tests. As pertinent to the issue under analysis in this paper, we present an analysis on the gravitational lensing phenomena within this framework.Although the main goal of this paper is to give a potential solution to the Dark Matter Problem within galaxies, we add a section where we show that an important property of the Bullet Cluster can in principle be explained in the scenario introduced in this work.To the end, we discuss the generic prospective to give rise to the Dark Matter component of most galaxies within extended gravity.     

The time-scale for core collapse in spherical star clusters

The collapse time for a cluster of equal-mass stars is usually stated to be either 330 central relaxation times (t_rc) or 12-19 half-mass relaxation times (t_rh). But the first of these times applies only to the late stage of core collapse, and the second only to low-concentration clusters. To clarify how the time depends on the density profile, the Fokker-Planck equation is solved for the evolution of a variety of isotropic cluster models, including King models, models with power-law density cusps of ρ ∼ r^−γ, and models with nuclei. The collapse times for King models vary considerably with the cluster concentration when expressed in units of t_rc or t_rh, but vary much less when expressed in units of t_rc divided by a dimensionless measure of the temperature gradient in the core. Models with cusps have larger temperature gradients and evolve faster than King models, but not all of them collapse: those with 0 < γ < 2 expand because they start with a temperature inversion. Models with nuclei collapse or expand as the nuclei would in isolation if their central relaxation times are short; otherwise their evolution is more complicated. Suggestions are made for how the results can be applied to globular clusters, galaxies, and clusters of dark objects in the centers of galaxies.Scott D. Tremaine     

Ion-acoustic solitons in pair-ion plasma with non-thermal electrons

Properties of fully nonlinear ion-acoustic solitary waves in an unmagnetized and collisionless pair-ion (PI) plasma containing superthermal electrons obeying Cairns distribution have been analyzed. A linear biquadratic dispersion relation has been derived, which yields the fast (supersonic) and slow (subsonic) modes in a pair-ion-electron plasma with nonthermal electrons. For nonlinear analysis, Korteweg-de Vries equation is obtained using the reductive perturbation technique. It is found that in case of slow mode, both electrostatic hump and dip type structures are formed depending on the temperature difference between positively and negatively charged ions, whereas, only dip type solitary structures have been observed for fast mode. The present work may be employed to explore and to understand the formation of solitary structures in the space (especially, the Earth’s ionosphere where two distinct pair ion species (H ^±) are present) and laboratory produced pair-ion plasmas with nonthermal electrons.     

Charged shell supported by a phantom energy

This paper discusses the evolution of a thin spherically symmetric self gravitating phantom shell around the charged shell. The general equations describing the motion of shell with a general form of equation of state are derived. The different types of space-time R _± and T _± regions and shell motion are classified depending on the parameters of the problem. The mechanical stability analysis of this spherically symmetric thin shell with charge in Reissner-Nordstrom (RN) to linearized spherically symmetric perturbation about static equilibrium solution is carried out.     

Head-on collision between positron acoustic waves in homogeneous and inhomogeneous plasmas

The head-on collision between positron acoustic solitary waves (PASWs) as well as the production of rogue waves (RWs) in homogeneous and PASWs in inhomogeneous unmagnetized plasma systems are investigated deriving the nonlinear evolution equations. The plasmas are composed of immobile positive ions, mobile cold and hot positrons, and hot electrons, where the hot positrons and hot electrons are assumed to follow the Kappa distributions. The evolution equations are derived using the appropriate coordinate transformation and the reductive perturbation technique. The effects of concentrations, kappa parameters of hot electrons and positrons, and temperature ratios on the characteristics of PASWs and RWs are examined. It is found that the kappa parameters and temperature ratios significantly modify phase shifts after head-on collisions and RWs in homogeneous as well as PASWs in inhomogeneous plasmas. The amplitudes of the PASWs in inhomogeneous plasmas are diminished with increasing kappa parameters, concentration and temperature ratios. Further, the amplitudes of RWs are reduced with increasing charged particles concentration, while it enhances with increasing kappa- and temperature parameters. Besides, the compressive and rarefactive solitons are produced at critical densities from KdV equation for hot and cold positrons, while the compressive solitons are only produced from mKdV equation for both in homogeneous and inhomogeneous plasmas.     

Hydrogen collisions in planetary atmospheres, ionospheres, and magnetospheres

Hydrogen is the most abundant element in the universe. Molecular hydrogen is the dominant chemical species in the atmospheres of the giant planets. Because of their low masses, neutral and ionized hydrogen atoms are the dominant species in the high atmospheres of many planets. Finally, protons are the principal heavy component of the solar wind.Here we present a critical evaluation of the current state of understanding of the chemical reaction rates and collision cross sections for several important hydrogen collision processes in planetary atmospheres, ionospheres, and magnetospheres. Accurate ab initio quantum theory will play an important role. The collision processes are grouped as follows:

(a)

H⁺+H charge transfer,

(b)

H⁺+H₂(v) charge transfer and vibrational relaxation, and

(c)

H₂(v,J)+H₂ vibrational, rotational, and ortho-para relaxation.

In each case we provide explicit representations as tabulations or compact formulas. Particularly important conclusions are that H⁺+H₂(v) collisions are more likely to result in vibrational relaxation than charge transfer and H₂ ortho-para conversion is at least an order-of-magnitude faster than previously assumed.     

Transport equations for aeronomy

In this paper we present results for a general system of transport equations appropriate to a multi-constituent gas mixture. This system includes a continuity, momentum, internal energy, pressure tensor and heat flow equation for each species. The results can be applied to both collision dominated and collisionless plasmas with there being explicit limits derived for the validity of the various expressions. In the limit of very frequent collisions the pressure tensor and heat flow equations give the usual Navier-Stokes results for the viscous stress tensor and heat flow vector. Furthermore, the momentum equation includes thermal diffusion and thermoelectric transport coefficients equivalent to the second approximation of Chapman and Cowling. The basic system of equations has been applied to different regions of the ionosphere and neutral atmosphere. It is found that: (1) The viscous stress tensor and heat flow expressions used in previous studies of the neutral thermosphere may not be appropriate; (2) The transport coefficients normally used for mid-latitude F2-region and topside studies seem to be adequate; (3) The high speed flow of plasma in the polar topside ionosphere is likely to be strongly affected by stresses and heat flow; and (4) E- and F-region ionization at high latitudes is substantially affected by stresses and heat flow.     

Physical processes in Jupiter's ring: Clues to its origin by Jove!

The particles making up the Jovian ring may be debris which has been excavated by micrometeoroids from the surfaces of many unseen (R ? 1 km) parent bodies (or “mooms” as we will occasionally call them) residing in the ring. A distribution of particle sizes exists: large objects are sources for the small visible ring particles and also account for the absorption of charged particles noted by Pioneer; the small grains are generated by micrometeoroid impacts, by jostling collisions among different-sized particles, and by self-fracturing due to electrostatic stresses. The latter are most effective in removing surface asperities to thereby produce smooth and crudely equidimensional grains. The presence of intermediate-sized (radius of several to several hundred microns) objects is also expected; these particles will have a total area comparable to the area of the visible ring particles. The nominal size (?2 μm) of the visible particles derived from their forward-scattering characteristics is caused, at least in part, by a selection effect but may also reflect a fundamental grain size or the preferential generation of certain sizes along with the destruction of others. The tiny ring particles have short lifetimes (?10²?10³ years) limited by erosion due to sputtering and meteoroid impacts. Plasma drag significantly modifies orbits in ～10² years but Poynting-Robertson drag is not effective (T_PR ～ 10⁵ years) in removing debris. The ring width is influenced by the distribution of source satellites, by the initial ejection velocity off them, by electromagnetic scattering, and by solar radiation forces. In the absence of electromagnetic forces, debris will reimpact a mother satellite or collide with another particle in about 10 years. A relative drift between different-sized particles, caused by a lessened effective gravity due to the Lorentz force, will substantially shorten these times to less than a month. The ring thickness is determined by a balance between initial conditions (abetted perhaps by electromagnetic scattering) and collisional damping; existence of the “halo” over the diffuse disk compared to its relative absence over the bright ring indicates the presence of mooms in the bright ring but not in the faint disk. Small satellites (R ? 1 km) will not reaccumulate colliding dust grains whereas satellites having the size of J14 or J16 may be able to do so, depending upon their precise shape, size, density, and location. Visible ring structure could indicate separate source satellites. The particles in the faint inner disk are delivered from the bright ring by orbital evolution principally under plasma drag. The halo is comprised of small particles (～0.1 μm) partially drawn out of the faint disk by interactions with the tilted Jovian magnetic field.     

Influence of inelastic collisions with hydrogen atoms on the formation of AlI and SiI lines in stellar spectra

We have performed calculations by abandoning the assumption of local thermodynamic equilibrium (within the so-called non-LTE approach) for Al I and Si I with model atmospheres corresponding to stars of spectral types F–G–Kwith differentmetal abundances. To take into account inelastic collisions with hydrogen atoms, for the first time we have applied the cross sections calculated by Belyaev et al. using model approaches within the formalism of the Born–Oppenheimer quantum theory. We show that for Al I non-LTE leads to higher ionization (overionization) than in LTE in the spectral line formation region and to a weakening of spectral lines, which is consistent with earlier non-LTE studies. However, our results, especially for the subordinate lines, differ quantitatively from the results of predecessors. Owing to their large cross sections, the ion-pair production and mutual neutralization processes Al I(nl) + HI(1s) ? Al II(3s ²) + H^? provide a close coupling of highly excited Al I levels with the Al II ground state, which causes the deviations from the equilibrium level population to decrease compared to the calculations where the collisions only with electrons are taken into account. For three moderately metal-deficient dwarf stars, the aluminum abundance has been determined from seven Al I lines in different models of their formation. Under the assumption of LTE and in non-LTE calculations including the collisions only with electrons, the Al I 3961 ?A resonance line gives a systematically lower abundance than the mean abundance from the subordinate lines, by 0.25–0.45 dex. The difference for each star is removed by taking into account the collisions with hydrogen atoms, and the rms error of the abundance derived from all seven Al I lines decreases by a factor of 1.5–3 compared to the LTE analysis. We have calculated the non- LTE corrections to the abundance for six subordinate Al I lines as a function of the effective temperature (4500 K ≤ T _eff ≤ 6500 K), surface gravity (3.0 ≤ log g ≤ 4.5), and metal abundance ([M/H] = 0, ?1, ?2, and ?3). For Si I including the collisions with HI leads to the establishment of equilibrium populations in the spectral line formation region even in hot metal-deficient models and to vanishingly small departures from LTE in spectral lines.     

Soliton and shocks in pair ion plasma in presence of superthermal electron

Solitons and shocks are addressed in a pair ion plasma in the presence of a kappa distribution. The dissipation is taken care of through the kinematic viscosity of both positive and negative ions in the plasma. The Kadomtsev–Petviashvili–Burger (KPB) equation is derived using the small amplitude expansion method. The Abel equation is obtained from the KPB equation and a solution is obtained by using the factorization method. The effect of the parameters κ and β (temperature ratio of ion species) is observed. Analytically we can find both solitons and shocks. The change of profile from soliton to shocks is shown in the figures. This study may be of wide relevance for the study of the formation of shocks and solitons in laboratory-produced pair ion plasmas.     

Power distributions for self-gravitating astrophysical systems based on nonextensive Tsallis kinetics

The long-time development of self-gravitating gaseous astrophysical systems (in particular, the evolution of the protoplanet accretion disk) is mainly determined by relatively fast processes of the collision relaxation of particles. However, slower dynamical processes related to force (Newton or Coulomb) interactions between particles should be included (as q-collisions) in the nonextensive kinetic theory as well. In the present paper, we propose a procedure to include the Newton self-gravity potential and the centrifugal potential in the near-equilibrium power-like q-distribution in the phase space, obtained (in the framework of nonextensive statistics) by means of the modified Boltzmann equation averaged with respect to an unnormalized distribution. We show that if the power distribution satisfies the stationary q-kinetic equation, then the said equation imposes clear restrictions on the character of the long-term force field and on the possible dependence of hydrodynamic parameters of the coordinates: it determines those parameters uniquely. We provide a thermodynamic stability criterion for the equilibrium of the nonextensive system. The results allow us to simulate the evolution of gaseous astrophysical systems (in particular, the gravitational stability of rotating protoplanet accretion disks) more adequately.