Observations of energetic water-group ions at comet giacobini-zinner: Implications for ion acceleration processes

Observations of energetic water-group pick-up ions made by the EPAS instrument during the ICE fly-by of comet P/Giacobini-Zinner are investigated for evidence concerning the processes which accelerate the ions from initial pick-up energies of around 10 keV up to energies of a few hundred kilo-electronvolts. The form of the ion spectrum in the ion rest frame is first investigated and compared with theoretical suggestions that exponential energy distributions might be produced by either first or second order Fermi acceleration in the cometary environment. It is shown that such distributions do not fit the data at all well, but that rather (over the EPAS ion bulk rest frame energy range of ∼30 to <300 keV) the observed distribution functions closely approximate an exponential in ion speed. The higher energy (> 50 keV) data also fit a power law distribution very well, but at lower energies the data tend to show a flattening below the power law form. It is also shown that the observed spectra are much softer than those calculated by Ip and Axford (1986, Planet. Space Sci. 34, 1061) for the G-Z ion pick-up region upstream of the bow shock, indicating that if these distributions are indeed due to second order Fermi acceleration as they propose, then the ion mean free path must be rather larger than used in their calculations (i.e. rather larger than 5 water-group ion gyroradii). Overall, it is concluded that at the present stage of theoretical development it is premature to draw firm conclusions about acceleration mechanisms from studies of the ion spectrum alone. However, from the general spectral analysis it is also possible to investigate the variations of ion intensity and spectral hardness which take place during the comet encounter, and to gain an indication of the degree of isotropy of the ion distribution in the rest frame of the flow. We find that a ∼.5 × 10⁴km-wide region (∼ 35 min along the spacecraft trajectory) of rapid ion intensification and spectral hardening occurs immediately upstream from the turbulent mass-loaded region, suggestive of a first order Fermi process in which ions are successively reflected between the outer layer of the slowed, turbulent region and waves in the faster upstream flow. It is shown that within the turbulent region the ion distribution becomes rapidly isotropized in its own bulk rest frame, indicative of a sudden reduction in ion mean free path. Additional processes (possibly second order Fermi acceleration) must also occur in order to account for the further modest intensifications of the energetic ion fluxes observed within the turbulent mass-loaded region, as well as the presence of ions in the outer pick-up region which have energies well in excess of the local pick-up energy.     

Interaction of interstellar pick-up ions with the solar wind

The interaction of interstellar pick-up ions with the solar wind is studied by comparing a model for the velocity distribution function of pick-up ions with actual measurements of He⁺ ions in the solar wind. The model includes the effects of pitch-ang'e diffusion due to interplanetary Alfvén waves, adiabatic deceleration in the expanding solar wind and the radial variation of the source function. It is demonstrated that the scattering mean free path is in the range 0.1 AU and that energy diffusion can be neglected as compared with adiabatic deceleration. The effects of adiabatic focusing, of the radial variation of the neutral density and of a variation of the solar wind velocity with distance from the Sun are investigated. With the correct choice of these parameters we can model the measured energy spectra of the pick-up ions reasonably well. It is shown that the measured differential energy density of the pick-up ions does not vary with the solar wind velocity and the direction of the interplanetary magnetic field for a given local neutral gas density and ionization rate. Therefore, the comparison of the model distributions with the measurements leads to a quantitative determination of the local interstellar gas density.Paper dedicated to Professor Hannes Alfvén on the occasion of his 80th birthday, 30 May 1988.     

Ejection of nitrogen from Titan's atmosphere by magnetospheric ions and pick-up ions

A 3-D Monte Carlo model is used to describe the ejection of N and N₂ from Titan due to the interaction of Saturn's magnetospheric N⁺ ions and molecular pick-up ions with its N₂ atmosphere. Based on estimates of the ion flux into Titan's corona, atmospheric sputtering is an important source of both atomic and molecular nitrogen for the neutral torus and plasma in Saturn's outer magnetosphere, a region now being studied by the Cassini spacecraft.     

Observations of energetic ions near the Venus ionopause

Ions (primarily O⁺) with spacecraft rest frame energies >40 eV have been observed by the Pioneer Venus Neutral Mass Spectrometer. The signature occurs in about 13% of the 700 orbits examined, primarily near the ionopause and at all solar zenith angles. The energetic ions coincide in location with superthermal ions observed by the Ion Mass Spectrometer and more rarely occur in some of the plasma clouds observed by the Electron Temperature Probe. These observations in conjunction with measurements by the Plasma Wave Instrument near the ionopause suggest that the ions are accelerated out of ionospheric plasma by the shocked solar wind through plasma waveparticle interactions.     

The spatial distribution of energetic ions and electrons in the magnetotail

The spatial distributions of energetic ion and electron bursts observed on the IMP 7 and 8 satellites in the Earth's magnetotail were studied. It was found that the ion bursts were more frequently detected in the dusk than in the dawn quarter of the neutral sheet whereas the electron bursts, more frequently in the dawn than the dusk quarter. The degree of dawn-dusk asymmetry is however energy dependent; the distribution for higher energy particle bursts exhibits higher degree of asymmetry. The morphologies of the distributions manifest themselves as seasonal variations of the most probable solar ecliptic latitudes at which the ion and electron bursts were observed. The amplitudes of the variations are about 25° with the seasonal variation for ions leading that for electrons by about 2 months.     

Electron capture of nuclei in the region surrounding magnetar

Effects of ultra-strong magnetic field on electron capture rates for ⁵⁷Fe, ⁵⁸Co and ⁵⁹Ni have been analyzed in the nuclear shell model and under the Landau energy levels quantized approximation in the ultra-strong magnetic field, the result increase about 3 orders magnitude. The rate of change of electron abundance, $\dot{Y}_{e}$ , for every nuclide and total $\dot{Y}_{e}$ in the condition without magnetic field and B=4.414×10¹⁵ G have been calculated, and exceed about 6 orders of magnitude generally. These conclusions play an important role in future studying the evolution of magnetar.     

Observations of energetic ions in inverted V events

Simultaneous measurements of keV ions and electrons with the ESRO 1A satellite have shown the following ion characteristics among others. Ions of about 6 keV energy are strongly field-aligned on the flanks of the inverted V events (mainly through the disappearance of the ion flux near 90° pitch angle). Field-aligned electron fluxes are often found in the same regions of the inverted V events where the ions are field-aligned. At the centre of inverted V events isotropization occurs (except in some small events). The 1 keV ion flux at large pitch angles (80°) is generally not reduced very much when the 6 keV, 80° ion flux shows strongly decreased values. The ratio of the 1 to 6 keV ion flux has a maximum near the centre of an inverted V event where the electron spectrum is hardest and the 6 keV ions are isotropic (or nearly isotropic).The observations are interpreted in terms of a model with two oppositely directed field-aligned electrostatic potential drops: one upper accelerating electrons downward and one lower, produced by the electron influx, which accelerates ions downward. Ion scattering in turbulent wave fields is proposed to be responsible for the observation that the 1 keV ion flux at large pitch angles does not decrease strongly where the 6 keV ion flux does and as an explanation of the isotropization at the centre of the event. The source problem for the ions is eliminated by the precipitating electrons ionizing continuously the thin neutral atmosphere even at altitudes of a few thousand kilometers.     

Representative abundances of a few positive ions in the innermost coma of comet Halley

The production rate of H₂O molecules at a heliocentric distance of 1 AU for comet Halley and the abundance ratio with respect to water (H₂O) of parent molecules at the cometary nucleus from the paper of Yamamoto (1987) have been used to compute the number densities of positive ions viz. H₃O⁺, H₃S⁺, H₂CN⁺, H₃CO⁺, CH₃OH ₂ ⁺ and NH ₄ ⁺ at various cometocentric distances within 600 kms from the nucleus.The role of proton transfer reactions in producing major ionic species is discussed. A major finding of the present investigation is that NH ₄ ⁺ ion which may be produced through proton transfer reactions is the most abundant ion near the nucleus of a comet unless the abundance of NH₃ as a parent is abnormally low. Using the quoted value of Q(NH₃)/Q(H₂O) for comet Halley and the life times of NH₃ and H₂O molecules, the abundance ratio N(NH₃)/N(H₂O) is found to be one-third of that used in the present paper. The consequent proportionate decrease in the NH ₄ ⁺ ions does not, however, affect its superiority in number density over other ions near the nucleus.The number density of the next most abundant ion viz. H₃O⁺ is found to be 4 × 10⁴ cm^-3 at the nucleus of comet Halley and decreases by a factor of two only upto a distance of 600 K ms from the nucleus. The ionic mass peak recorded by VEGA and GIOTTO spacecrafts atm/q = 18 is most probably composite of the minor ionic species H₂O⁺, as its number density = 10² cm^-3 remains virtually constant in the inner coma and of NH ₄ ⁺ , the number density of which at large cometocentric distances may add to the recorded peak atmlq = 18. The number densities of other major ions produced through proton transfer from H₃O⁺ are also discussed in the region within 600 K ms from the nucleus of comet Halley.     

3C 120 and the surrounding region of sky

Image processing performed on a series of photographs of the superluminal Seyfert galaxy, 3C 120, shows the outer optical disc to consist of fragmented segments generally pointing toward the centre. One long arm of peculiar, separated knots comes off to the W and SW. A peculiar companion is seen along the line of the NW radio jet. In the interior, optical jets are detected which are aligned along the direction of the outer radio jets. A region of the sky 45 ×; 25 degrees around 3C120 is investigated. It is found that:

1. A nebulous filament about 3/4 degree in length points to 3C 120.
2. Hydrogen clouds of redshiftz = ?130 and ?210 km s^?1 are situated at 3 and 1 degrees on either side of 3C 120.
3. Eleven low-surface-brightness galaxies with 4500 <z < 5300 km s^?1 fall within a radius of 8 degrees.
4. Seven quasars withz ? 1.35 and radio fluxesS _b ? 0.3 fall within a radius of 10 degrees.

It is concluded that the concentration of these objects in the vicinity of this unique, active galaxy has a negligible chance of being accidental and that all those objects of diverse redshift are at the same nearby distance. This smaller distance reduces the supposed superluminal motions in 3C 120 to quite precedented ejection velocities.     

On the origin and excitation of the extended nebula surrounding NGC 1275

We use line-of-sight velocity information on the filamentary emission-line nebula of NGC 1275 to infer a dynamical model of the nebula's flow through the surrounding intracluster gas. We detect outflowing gas and flow patterns that match simulations of buoyantly rising bubbles from which we deduce that some of the nebula filaments have been drawn out of NGC 1275. We find a radial gradient of the ratio [N  ii ]λ6584/Hα which may be due to a variation in metallicity, interactions with the surrounding intracluster medium or a hardening of the excitation mechanism. We find no preferred spatial correlation of stellar clusters within the filaments and there is a notable lack of [O  iii ]λ5007 emission, therefore it is unlikely that the filaments are ionized by stellar ultraviolet.     

Long-period comet flux in the planetary region: Dynamical evolution from the Oort cloud

This study analyzes the evolution of 2 × 10⁵ orbits with initial parameters corresponding to the orbits of comets of the Oort cloud under the action of planetary, galactic, and stellar perturbations over 2 × 10⁹ years. The dynamical evolution of comets of the outer (orbital semimajor axes a > 10⁴ AU) and inner (5 × 10³ < a (AU) < 10⁴) parts of the comet cloud is analyzed separately. The estimates of the flux of “new” and long-period comets for all perihelion distances q in the planetary region are reported. The flux of comets with a > 10⁴ AU in the interval 15 AU < q < 31 AU is several times higher than the flux of comets in the region q < 15 AU. We point out the increased concentration of the perihelia of orbits of comets from the outer cloud, which have passed several times through the planetary system, in the Saturn-Uranus region. The maxima in the distribution of the perihelia of the orbits of comets of the inner Oort cloud are located in the Uranus-Neptune region. “New” comets moving in orbits with a < 2 × 10⁴ AU and arriving at the outside of the planetary system (q > 25 AU) subsequently have a greater number of returns to the region q < 35 AU. The perihelia of the orbits of these comets gradually drift toward the interior of the Solar System and accumulate beyond the orbit of Saturn. The distribution of the perihelia of long-period comets beyond the orbit of Saturn exhibits a peak. We discuss the problem of replenishing the outer Oort cloud by comets from the inner part and their subsequent dynamical evolution. The annual rate of passages of comets of the inner cloud, which replenish the outer cloud, in the region q < 1 AU in orbits with a > 10⁴ AU (～ 5.0 × 10^?14 yr^?1) is one order of magnitude lower than the rate of passage of comets from the outer Oort cloud (～ 9.1 × 10^?13 yr^?1).     

Energization of positive ions in the cometary foreshock region

The energization of positive ions in front of a cometary bow shock is investigated. Ions produced by ionization of the cometary neutrals interact with the solar wind protons to produce, among other waves, large amplitude oscillations of the ambient magnetic field. Such oscillations are convected towards the comet at the unperturbed solar wind speed far from the shock and at a lower speed closer to the shock (due to the solar wind mass loading) ; hence, they can energize the suprathermal ions by Fermi acceleration. The spatial extension of the acceleration region is of the order of 10⁶ km and the resulting ion energy spectrum is harder than in the Earth's bow shock case. The energization of cometary ions produces an additional deceleration of the solar wind, such that the cometary bow shock of Halley-type comet may be regarded as a “cosmic ray shock”.     

Gasdynamic interpretation of ice measurements at the Giacobini-Zinner encounter

Plasma measurements of the ICE spacecraft at its encounter with comet Giacobini-Zinner are found to be in qualitative agreement with corresponding predictions of a gasdynamic model on the solar wind interaction with the cometary neutral gas.     

The dust coma of Comet P/Giacobini-Zinner in the infrared

We present 1-20 micrometers photometry of P/Giacobini-Zinner obtained at the NASA Infrared Telescope Facility, during 1985 June-September (r = 1.57-1.03 AU). A broad, weak 10 micrometers silicate emission feature was detected on August 26.6; a similar weak emission feature could have been hidden in the broadband photometry on other dates. The total scattering and emitting cross section of dust in the inner coma was similar to that in other short-period comets, but a factor of 10 (r = 1.56 AU) to 100 (r = 1.03 AU) lower than the amount of dust in Comet Halley. The thermal emission continuum can be fit with models weighted toward either small or large absorbing grains. The dust production rate near perihelion was approximately 10(5) g/s (small-grain model) to approximately 10(6) g/s (large-grain model). The corresponding dust/gas mass ratio on August 26 was approximately 0.1-1. A silicate-rich heterogeneous grain model with an excess of large particles is compatible with the observed spectrum of Giacobini-Zinner on August 26. Thus, weak or absent silicate emission does not necessarily imply an absence of silicates in the dust, although the abundance of silicate particles < or = 1 micrometer radius must have been lower than in Comet Halley.     

Reflection of pre-accelerated pick-up ions at the solar wind termination shock: The seed for anomalous cosmic rays

It has been hypothesized for quite some time that interplanetary pick-up ions due to energization taking place in the region close to the solar wind termination shock, at some fraction and as an outcome of a complicated chain of processes, eventually are converted into species of the anomalous cosmic-ray particles. For the actual conversion efficiency it is of great importance to know the energy distribution of these pick-up ions upon their arrival at the shock. It turns out that pre-acceleration of these ions during their passage through the heliosphere shall substantially increase their chances to become reflected at the shock into the upstream direction which is a prerequisite for a further climb-up in energy by virtue of Fermi-1 acceleration processes. In this paper we start out from stochastically pre-accelerated pick-up ions and investigate their behaviour at the shock. With the use of adiabatic approaches in the de Hoffman-Teller frame of the shock, we calculate the energy distribution function of the reflected part of pick-up ions. From the calculated distribution functions it turns out that the reflected ions in the average suffer an energy increase by about a factor of 10, still not enough to let them move off the shock by spatial diffusion in the upstream direction. Thus, since converted back into the shock, they can undergo repeated reflections and energy gains till the diffusion-convection limit is reached. As we show in addition, the reflection probability for pick-up ions is about a factor of 10 higher than expected from the present literature and strongly varies with the off-axis angle, pointing to the fact that the termination shock represents a surface with a three-dimensionally varying source strength for the production of anomalous cosmic rays. The ACR source pattern is also expected to vary during the solar cycle and the relevant injection energies are expected to be larger by factors of 10 to 100 than the canonically adopted 1 keV nucl^–1.Institute for Problems of Mechanics of the Russian Academy of Sciences, Prospect Vernadskogo 101, 117526, Moscow, Russia.     

Steady-state observations of geomagnetically trapped energetic heavy ions and their implications for theory

Measurements of energetic heavy ions using the Explorer 45 and ATS-6 satellites are reviewed and the resulting implications for theory are evaluated. The measured ions are basically protons and helium ions in the energy range from 0.1 to 1 MeV/nucleon. The equatorial energetic ion distributions inside L = 4.5 are found to be very stable for extended periods of time. These ions are very closely confined to the equatorial plane and are sharply peaked as a function of L around a value designated as L_max. Beyond L = 5.0 the fluxes of these ions are more variable with order of magnitude variations being observed at L = 6.6 on the time scales of minutes, hours, or days. The region inside L = 4.5 appears to be well described by radial diffusive transport driven by fluctuations in the geomagnetic field coupled with losses due to charge exchange and Coulomb interactions with ambient hydrogen geocorona and terrestrial plasma environment. From an analysis relating the position in L-value of the maximum intensity, L_max, observed for a given ion species and energy, it is argued that the influence of fluctuations in the convection electric field as discussed by Cornwall (1972) are not effective in radially diffusing in L ions with energies greater than a few hundred kiloelectron volts per nucleon. The source of these ions remains basically undetermined and its determination must await further measurements.     

Pickup interstellar helium ions in the region of the solar gravitational cone

Our numerical analyses of the velocity and spatial distributions of pickup interstellar helium ions in the region of the solar gravitational cone in the ecliptic plane at a distance of 1 AU show that the ion density maximum must be displaced relative to the neutral helium cone axis in the direction of the Earth’s revolution around the Sun. The solar wind parameters in the numerical model correspond to their observed values during the crossing of the helium cone by the ACE spacecraft in 1998. At these parameters, the calculated angular displacement is 5°. The absence of a similar displacement in the ACE measurements is shown to stem from the fact that the spectrometer onboard ACE records and identifies only a fraction of the pickup helium ions with fairly high magnitudes and certain directions of the velocities.