Energy dissipation and ion heating at the heliospheric termination shock

The Los Alamos hybrid simulation code is used to examine heating and the partition of dissipation energy at the perpendicular heliospheric termination shock in the presence of pickup ions. The simulations are one-dimensional in space but three-dimensional in field and velocity components, and are carried out for a range of values of pickup ion relative density. Results from the simulations show that because the solar wind ions are relatively cold upstream, the temperature of these ions is raised by a relatively larger factor than the temperature of the pickup ions. An analytic model for energy partition is developed on the basis of the Rankine-Hugoniot relations and a polytropic energy equation. The polytropic index gamma used in the Rankine-Hugoniot relations is varied to improve agreement between the model and the simulations concerning the fraction of downstream heating in the pickup ions as well as the compression ratio at the shock. When the pickup ion density is less than 20%, the polytropic index is about 5/3, whereas for pickup ion densities greater than 20%, the polytropic index tends toward 2.2, suggesting a fundamental change in the character of the shock, as seen in the simulations, when the pickup ion density is large. The model and the simulations both indicate for the upstream parameters chosen for Voyager 2 conditions that the pickup ion density is about 25% and the pickup ions gain the larger share ( approximately 90%) of the downstream thermal pressure, consistent with Voyager 2 observations near the shock

Particle Energization in an Expanding Magnetized Relativistic Plasma

Using a 2-1/2-dimensional particle-in-cell (PIC) code to simulate the relativistic expansion of a magnetized collisionless plasma into a vacuum, we report a new mechanism in which the magnetic energy is efficiently converted into the directed kinetic energy of a small fraction of surface particles. We study this mechanism for both electron-positron and electron-ion (mi/me=100, me is the electron rest mass) plasmas. For the electron-positron case the pairs can be accelerated to ultra-relativistic energies. For electron-ion plasmas most of the energy gain goes to the ions.Comment: 7 pages text plus 5 figures, accepted for publication by Physical Review Letter

The moving boundary problem in the presence of a dipole magnetic field

An exact analytic solution is obtained for a uniformly expanding, neutral, infinitely conducting plasma sphere in an external dipole magnetic field. The electrodynamical aspects related to the radiation and transformation of energy were considered as well. The results obtained can be used in analyzing the recent experimental and simulation data.Comment: 17 pages, 1 figure, Submitted to J. Phys. A, Math. and Genera

An exact solution of the moving boundary problem for the relativistic plasma expansion in a dipole magnetic field

An exact analytic solution is obtained for a uniformly expanding, neutral, highly conducting plasma sphere in an ambient dipole magnetic field with an arbitrary orientation of the dipole moment in the space. Based on this solution the electrodynamical aspects related to the emission and transformation of energy have been considered. In order to highlight the effect of the orientation of the dipole moment in the space we compare our results obtained for parallel orientation with those for transversal orientation. The results obtained can be used to treat qualitatively experimental and simulation data, and several phenomena of astrophysical and laboratory significance.Comment: 7 pages, 2 figures. arXiv admin note: substantial text overlap with arXiv:physics/060323

Simulations of Electron Acceleration at Collisionless Shocks: The Effects of Surface Fluctuations

Energetic electrons are a common feature of interplanetary shocks and planetary bow shocks, and they are invoked as a key component of models of nonthermal radio emission, such as solar radio bursts. A simulation study is carried out of electron acceleration for high Mach number, quasi-perpendicular shocks, typical of the shocks in the solar wind. Two dimensional self-consistent hybrid shock simulations provide the electric and magnetic fields in which test particle electrons are followed. A range of different shock types, shock normal angles, and injection energies are studied. When the Mach number is low, or the simulation configuration suppresses fluctuations along the magnetic field direction, the results agree with theory assuming magnetic moment conserving reflection (or Fast Fermi acceleration), with electron energy gains of a factor only 2 - 3. For high Mach number, with a realistic simulation configuration, the shock front has a dynamic rippled character. The corresponding electron energization is radically different: Energy spectra display: (1) considerably higher maximum energies than Fast Fermi acceleration; (2) a plateau, or shallow sloped region, at intermediate energies 2 - 5 times the injection energy; (3) power law fall off with increasing energy, for both upstream and downstream particles, with a slope decreasing as the shock normal angle approaches perpendicular; (4) sustained flux levels over a broader region of shock normal angle than for adiabatic reflection. All these features are in good qualitative agreement with observations, and show that dynamic structure in the shock surface at ion scales produces effective scattering and can be responsible for making high Mach number shocks effective sites for electron acceleration.Comment: 26 pages, 12 figure

Long-range attraction between particles in dusty plasma and partial surface tension of dusty phase boundary

Effective potential of a charged dusty particle moving in homogeneous plasma has a negative part that provides attraction between similarly charged dusty particles. A depth of this potential well is great enough to ensure both stability of crystal structure of dusty plasma and sizable value of surface tension of a boundary surface of dusty region. The latter depends on the orientation of the surface relative to the counter-ion flow, namely, it is maximal and positive for the surface normal to the flow and minimal and negative for the surface along the flow. For the most cases of dusty plasma in a gas discharge, a value of the first of them is more than sufficient to ensure stability of lenticular dusty phase void oriented across the counter-ion flow.Comment: LATEX, REVTEX4, 7 pages, 6 figure

Non-Maxwellian Proton Velocity Distributions in Nonradiative Shocks

The Balmer line profiles of nonradiative supernova remnant shocks provide the means to measure the post-shock proton velocity distribution. While most analyses assume a Maxwellian velocity distribution, this is unlikely to be correct. In particular, neutral atoms that pass through the shock and become ionized downstream form a nonthermal distribution similar to that of pickup ions in the solar wind. We predict the H alpha line profiles from the combination of pickup protons and the ordinary shocked protons, and we consider the extent to which this distribution could affect the shock parameters derived from H alpha profiles. The Maxwellian assumption could lead to an underestimate of shock speed by up to about 15%. The isotropization of the pickup ion population generates wave energy, and we find that for the most favorable parameters this energy could significantly heat the thermal particles. Sufficiently accurate profiles could constrain the strength and direction of the magnetic field in the shocked plasma, and we discuss the distortions from a Gaussian profile to be expected in Tycho's supernova remnant.Comment: 13 pages, 6 figure