Voyager energetic particle observations at interplanetary shocks and upstream of planetary bow shocks: 1977–1990

The Voyager 1 and 2 spacecraft include instrumentation that makes comprehensive ion (E 28 keV) and electron (E 22 keV) measurements in several energy channels with good temporal, energy, and compositional resolution. Data collected over the past decade (1977–1988), including observations upstream and downstream of four planetary bow shocks (Earth, Jupiter, Saturn, Uranus) and numerous interplanetary shocks to 30 AU, are reviewed and analyzed in the context of the Fermi and shock drift acceleration (SDA) models. Principal findings upstream of planetary bow shocks include the simultaneous presence of ions and electrons, detection of tracer ions characteristic of the parent magnetosphere (O, S, O⁺), power-law energy spectra extending to 5 MeV, and large (up to 100:1) anisotropies. Results from interplanetary shocks include observation of acceleration to the highest energies ever seen in a shock ( 22 MeV for protons, 220 MeV for oxygen), the saturation in energy gain to 300 keV at quasi-parallel shocks, the observation of shock-accelerated relativistic electrons, and separation of high-energy (upstream) from low-energy (downstream) populations to within 1 particle gyroradius in a near-perpendicular shock. The overall results suggest that ions and electrons observed upstream of planetary bow shocks have their source inside the parent magnetosphere, with first order Fermi acceleration playing a secondary role at best. Further, that quasi-perpendicular interplanetary shocks accelerate ions and electrons most efficiently to high energies through the shock-drift process. These findings suggest that great care must be exercised in the application of concepts developed for heliosphere shocks to cosmic ray acceleration through shocks at supernova remnants.     

High-latitude observations of energetic ions during the first Ulysses polar pass

We present observations of energetic ions from the Ulysses COSPIN Low Energy Telescope in the mid and high-latitude regions of the heliosphere prior to and during the first polar pass of the Ulysses spacecraft. After the encounter with Jupiter, Ulysses started on its journey out-of-the-ecliptic. Between 13°S and 29°S the spacecraft sampled the solar wind from both the streamer belt and the polar coronal hole. Here, co-rotating magnetic structures with forward and reverse shocks and containing accelerated energetic ions were observed.At latitudes greater than 29°S, Ulysses was completely immersed in the solar wind from the polar coronal hole. Here the co-rotating magnetic structures were weaker, and in general had only reverse shocks, but were still capable of accelerating the energetic ions, albeit with reduced intensity. The most recent results show that beyond 50°S, very few if any, reverse shocks are observed. However, accelerated ions from magnetic interaction regions are still observed. We report also on an intensity enhancement at 50°S due to the passage of a high-latitude CME.     

Solar wind charge states measured by ULYSSES/SWICS in the south polar hole

The Ulysses mission now has an extensive data base covering several passes of the south polar coronal hole as the spacecraft proceeds to higher latitudes. Using composition measurements from the SWICS experiment on the Ulysses spacecraft, we have obtained charge state distributions, and hence inferred coronal ionization temperatures, for several solar wind species. In particular, we present an overview of Oxygen ionization temperature measurements, based on the O⁷⁺/O⁶⁺ ratio, for the period January 1993 until April 1994 (23°S to 61°S heliographic latitude), and detailed Oxygen, Silicon and Iron charge state distributions of the south polar hole during a two month period of nearly continuous hole coverage, Dec 1993–Jan 1994 (45°S to 52°S heliographic latitude).     

The Freja F3C Cold Plasma Analyzer

The F3C Cold Plasma Analyzer (CPA) instrument on theFreja spacecraft is designed to measure the energy per unit charge (E/Q) of ions oe electrons in the range 0<E/Q<200 V and complements the observations made by the F3H Hot Plasma Experiment. The CPA sensor, which is deployed on a boom, is an electrostatic analyzer which produces angle/energy images of particles incident on the sensor in a plane perpendicular to the boom axis. Charged particles incident normal to the CPA sensor housing axis of symmetry, which coincides with the boom axis, pass through collimators and enter a semi-spherical electrostatic analyzer which disperses particles in energy and azimuthal angle of arrival onto an imaging MCP detector thus producing images of the particle distributions in a plane perpendicular to the boom axis. Measurements are transmitted either as discrete 16×16 (angle/energy) images or as parameters related to the incident particle distribution function. Pixels in the discrete images are separated approximately equally in azimuthal angle while the 16 energy bins are separated approximately geometrically in energy. The ratio of the maximum to minimum energy imaged is programmable up to a maximum of more than a factor of ten, and the energy range itself is also under the control of the processor and can be varied by more than an order of magnitude. The density dynamic range of the sensor is increased by the introduction of an electrostatic gating system between the entrance aperture and the analyzer which can be used to duty-cycle low-energy electrons into the sensor thus keeping the count rate within appropriate levels. To reduce the effects of spacecraft induced perturbations on the lower-energy particle distributions, the sensor portion of the instrument is deployed on a 2 m long boom, perpendicular to the spacecraft spin axis. Spacecraft rotation is used to recover complete (4) angle/energy distributions every half spin period. In addition, the sensor skin may be biased with respect to the spacecraft ground to offset effects due to spacecraft charging. Current to the skin is monitored, making the exterior of the sensor equivalent to a large cylindrical Langmuir probe. Two separate processing paths for signals from the MCP anode may be chosen; slow and rast. The slow pulse processing path provides discrete angle/energy images at a nominal rate of 10 images per second and a peak burst mode rate of 100 images per second. The fast analog or current mode path provides crude parameterized estimates of densities, temperatures and drift velocities at nominal rates of up to 1000 parameters per second with a burst rate near 6000 parameters per second. Observations of cold ions and electrons in an unperturbed ionospheric plasma are presented which demonstrate the functionality of the instrument. Suprathermal ion observations in a transverse ion energization or acceleration region are also shown which demonstrate many of the small-scale features of these events.The Canadian Government's right to retain a non-exclusive, royalty free licence in and to any copyright is acknowledge.     

Gamma-ray line production in astronomical objects

This review focuses on the conditions for -ray line production in the most interesting astronomical objects, in light of the planned experiments: Gamma-1, GRO, Sigma, GRASP, and others. Among these objects are the Sun, the galactic center region, molecular and dust clouds, accreting and exploding stars.     

Simulated observations of the electron coronal density irregularity\overline {n^2 } /(\bar n)^2

A technique to derive the coronal density irregularity factor , wheren is the electron density, has been proposed by Fineschi and Romoli (1993). This technique will exploit the unique UVCS capability of cotemporal and cospatial measurements of both UV line radiation and K-coronal polarized brightness,pB.The ratio of the measured H I Lyman (Ly-) line intensity to the resonant-scattering dominated H I Lyman (Ly-) intensity can be used to extract the collisional component of the Ly-. This component yields an estimate of . The quantity is then obtained from the UVCS white-light K-coronal measurements.We present simulated observations of the UVCS for coronal atmosphere models with different filling factors and electron density profiles, and for different coronal structures (e.g., coronal holes, streamers). These simulations will show how the proposed technique may be used to probe inhomogeneities of the solar corona.     

The Medium Energy Instrument on Exosat

The Medium Energy Instrument on EXOSAT, although conceived as the main instrument for occultations, has been made sufficiently versatile to provide a significant advance over previous large area proportional counters when used for individual source studies of timing and spectra. The energy range is 1.2 to 50 keV, with E/E of 0.2 at 6 keV, sufficient to detect iron lines. The effective area of 1800 cm² and narrow field of view (3/4° × 3/4°) make it suitable for the detailed study of sources down to the 0.3 mCrab confusion limit. The unique facility provided by EXOSAT, allowing uninterrupted observations of X-ray sources for periods of up to 80 hours, backed up by a high capacity data link and on-board processing, enables timing studies to be performed over the range from milliseconds to days. Sophisticated background discrimination techniques giving a rejection efficiency of99% will control the background count rate to a suitably low value in the environment of the 200,000 km orbit.     

The solar grazing incidence (GRIST) and optical (SOT) telescopes joint accommodation

Some of the problems foreseen for the joint accommodation and operation of the Grazing Incidence Solar Telescope (GRIST) under study by ESA to operate in the extreme ultraviolet region (90 < < 1700 Å), and the Solar Optical Telescope (SOT), developed by NASA to operate in the ultraviolet, optical and infrared region (A > 1100 Å) on a Spacelab mission are described.Proceedings of the Conference Solar Physics from Space, held at the Swiss Federal Institute of Technology Zurich (ETHZ), 11–14 November 1980.     

Nonlinear theory of the two-point correlation function

We consider the influence of the nonlinear stage of gravitational instability on the two-point correlation functions of gravitationally bound objects. Based on the theory of nonlinear gravitational contraction of a single density peak of dissipationless matter (Gurevich and Zybin, 1988a,b; 1990) we develop a method for calculating the two-point correlation functions of different objects of any mass. The method works good in the region of strong correlations and can be easily extended to calculate higher correlation functions. We show that the main contribution to the correlation function _i in the region of strong correlations _i 1 is made by pair systems located outside large clusters of objects. In this region the shape of _i is determined only by the nonlinear dynamics of gravitational contraction of dissipationless matter and has the form _i C ^–, where 1.8 is a universal parameter.     

Some results of ionospheric investigations by means of coherent radiowaves emitted by satellites

In this paper we discuss theoretical expressions, determining the difference of Doppler shifts of various coherent radiowave frequencies emitted by a radiator moving in the ionosphere or interplanetary medium. The rotating Doppler effect (Faraday effect) caused by the Doppler shifts ±H of the ordinary and extraordinary waves is also considered. In a three-dimensional inhomogeneous ionosphere, stationary in time (N/t = 0), is determined in the general case, by an equation with three variables. The equation for proper depends only on the local value of the electron concentration N _c around the radiator and on integral values, determining, by means of additional calculations, the angle of refraction or its components, the horizontal gradients of electron concentration N/x and N/y, and in some cases, the integral electron concentration ₀ ^zcN dz. We describe the analysis of the measurements, made with the satellites Cosmos I, II and partially XI, assuming that N/t = N/y = 0, with a two variables equation. The expected errors are considered. The results coincide well for different points (Moscow, The Crimea, Sverdlovsk) and thus agree with the measurements of _H and with height-frequency ionospheric characteristics. The curve giving electron concentration versus height N (z) in the outer ionosphere (above the maximum of F2), shows a new maximum higher than the main maximum of the ionosphere N _MF2 at 120–140 km. At this maximum the value of N (z) is (0.9–0.95) N _MF2. The new data on the large-scale horizontal inhomogeneities of the ionosphere, exceed the previous ones by about a factor 10. By means of the irregular variations of the spectrum W() of the inhomogenous formation is determined. Three unknown constant maxima with values 16 to 18 km, 28 to 32 km and 100 to 120 km are found. The spectrum W () mainly characterizes the local properties of the ionosphere along the orbit of the satellite.     

Relative ionizations in the nearest interstellar gas

We compare CLOUDY predictions for the equilibrium ionization in the interstellar cloud surrounding the solar system with pick-up ion data. The incident radiation field includes contributions from hot stars, the emission from the conductive cloud boundary and the diffuse FUV back-ground. To within the observational uncertainties, CLOUDY predictions for the ratios n(He)/n(O), n(N)/n(O), n(Ne)/n(O), and n(He)/n(Ne) are consistent with pick-up ion data, provided that O and N are filtered by 50% in the heliopause region and the outer heliosphere as predicted by others. Thus, the steady-state ionization model and assumed radiation field appear approximately valid. However, the youth and low intervening column density towards the Vela pulsar leave open the possibility that the parent supernova explosion 10,500 years ago, and 200 pc distant, may also have affected LISM ionization, although the mechanism is uncertain. Support for this last possibility is provided by the apparent signature of the Vela explosion in the terrestrial geological record.Abbreviations ISM Interstellar Medium - FUV Far Ultraviolet - EUV Extreme Ultraviolet - SNR SN remnant - SXRB SXR Background - LISM Local Interstellar Matter     

A New Mechanism of Coronal Heating

The interaction between network magnetic fields and emerging intranetwork fields may lead to magnetic reconnection and microflares, which generate fast shocks with an Alfvén Mach number M _A<2. Protons and less abundant ions in the solar corona are then heated and accelerated by fast shocks. Our study of shock heating shows that (a) the nearly nondeflection of ion motion across the shock ramp leads to a large perpendicular thermal velocity (v _th), which is an increasing function of the mass/charge ratio; (b) the heating by subcritical shocks with 1.1 M_A 1.5 leads to a large temperature anisotropy with T/T 50 for O⁵⁺ ions and a mild anisotropy with T_/T 1.2 for protons; (c) the large perpendicular thermal velocity of He⁺⁺ and O⁵⁺ ions can be converted to the radial outflow velocity (u) in the divergent coronal field lines; and (d) the heating and acceleration by shocks with 1.1 M_A 1.5 can lead to u(O⁵⁺) v _th(O⁵⁺) 460 km s^–1 for O⁵⁺ ions, u(He⁺⁺) v _th(He⁺⁺) 360 km s^–1 for He⁺⁺ ions, and u(H⁺) v _th(H⁺) 240 km s^–1 for protons at r=3–4 R . Our results can explain recent SOHO observations of the heating and acceleration of protons and heavier ions in the solar corona.     

Theoretical modelling of Algol disks

A brief review of various theoretical approaches to model accretion disks is presented. Emphasis is given to models that determine self-consistently the structure of a disk together with the radiation field. It is argued that a proper treatment of the vertical structure is essential for calculating theoretical spectra to be compared with observations. In particular, it is shown that hot layers above an accretion disk (sometimes called disk chromospheres or coronae), whose presence is indicated by recent UV observations of strong emission lines of highly ionized species, may be explained using simple energy balance arguments.1987–88 JILA Visiting Fellow.This work was in part supported by a NASA grant ADP U-003-88 (Plavec and Hubeny). I also wish to thank the organizers of the IAU Colloquium 107 for the travel grant which enabled me to attend the meeting.     

Near Earth Current Meander (Necm) Model of Substorms

We propose that the appropriate instability to trigger a substorm is a tailward meander (in the equatorial plane) of the strong current filament that develops during the growth phase. From this single assumption follows the entire sequence of events for a substorm. The main particle acceleration mechanism in the plasma sheet is curvature drift with a dawn-dusk electric field, leading to the production of auroral arcs. Eventually the curvature becomes so high that the ions cannot negotiate the sharp turn at the field-reversal region, locally, at a certain time. The particle motion becomes chaotic, causing a local outward meander of the cross-tail current. An induction electric field is produced by Lenz's law, E ^ind=–A/t. An outward meander with B _z>0 will cause E×B flow everywhere out from the disturbance; this reaction is a macroscopic instability which we designate the electromotive instability. The response of the plasma is through charge separation and a scalar potential, E ^es=–. Both types of electric fields have components parallel to B in a realistic magnetic field. For MHD theory to hold the net E must be small; this usually seems to happen (because MHD often does hold), but not always. Part of the response is the formation of field-aligned currents producing the well-known substorm current diversion. This is a direct result of a strong E ^ind (the cause) needed to overcome the mirror force of the current carriers; this enables charge separation to produce an opposing electrostatic field E ^es (the effect). Satellite data confirm the reality of a strong E in the plasma sheet by counter-streaming of electrons and ions, and by the inverse ion time dispersion, up to several 100 keV. The electron precipitation is associated with the westward traveling surge (WTS) and the ion with omega () bands, respectively. However, with zero curl, E ^es cannot modify the emf =Edl=–d^M/dt of the inductive electric field E ^ind (a property of vector fields); the charge separation that produces a reduction of E must enhance the transverse component E . The new plasma flow becomes a switch for access to the free energy of the stressed magnetotail. On the tailward side the dusk-dawn electric field with EJ<0 will cause tailward motion of the plasma and a plasmoid may be created; it will move in the direction of least magnetic pressure, tailward. On the earthward side the enhanced dawn-dusk induction electric field with EJ>0 will cause injection into the inner plasma sheet, repeatedly observed at moderate energies of 1–50 keV. This same electric field near the emerging X-line will accelerate particles non-adiabatically to moderate energies. With high magnetic moments in a weak magnetic field, electrons (ions) can benefit from gradient and curvature drift to attain high energies (by the ratio of the magnetic field magnitude) in seconds (minutes).     

Gamma Rays From Molecular Clouds

High energy -rays from individual giant molecular clouds contain unique information about the hidden sites of acceleration of galactic cosmic rays, and provide a feasible method for study of propagation of cosmic rays in the galactic disk on scales 100 pc. I discuss the spectral features of ⁰-decay -radiation from clouds/targets located in proximity of relatively young proton accelerators, and speculate that such `accelerator+target systems in our Galaxy can be responsible for a subset of unidentified EGRET sources. Also, I argue that the recent observations of high energy -rays from the Orion complex contain evidence that the level of the `sea of galactic cosmic rays may differ significantly from the flux and the spectrum of local (directly detected) particles.     

Temperature Fluctuations,H II Region Abundances,and Primordial Helium

There is evidence for temperature fluctuations in Planetary Nebulae and in some Galactic H II regions. If such fluctuations occur in the low metallicity, extragalactic H II regions used to probe the primordial helium abundance, the derived 4He mass fraction, YP, could be systematically different from the true primordial value. Although this effect could be large, there are no data which allow us to estimate the size of the temperature fluctuations for the extragalactic H II regions. Therefore, we have explored this effect via Monte Carlo simulations of the data in which the abundances derived from a fiducial data set are modified by T chosen from a distribution with 0 T Tmax where Tmax is varied from 500 K to 4000 K.     

Theoretical studies of the flux and energy spectrum of gamma radiation from the Sun

The purpose of this work is to study the various -ray-production mechanisms in solar flares and to calculate the flux, the spectrum, and the decay curves of radiation. Using the continuity equation and taking into account the energy losses for solar-flare-accelerated particles, we obtain the time-dependent particle distribution and thus the time behavior of the resulting rays. The important processes for producing rays in solar flares are found to be nonthermal electron bremsstrahlung, decay of neutral mesons, positron annihilation, neutron capture, and decay of excited nuclei. The results are applied to several known solar flares. For a large flare such as the class 3⁺ on February 23, 1956, continuous rays with energies up to 100 MeV from electron bremsstrahlung and neutral meson decays are observable at the orbit of the Earth by existing -ray detectors. Line rays from positron annihilation (0.51 MeV), neutron capture (2.23 MeV), and deexcitation of excited nuclei O¹⁶ (6.14 and 7.12 MeV) and C¹² (4.43 MeV) are particularly strong and well above the continuous -ray background due to electron bremsstrahlung. These lines can be detected at the Earth.NASA-NRC Resident Research Associate.     

Pulsars as gamma ray sources: Nebular shocks and magnetospheric gaps

I summarize the results of recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is primarily electrons and positrons with an admixture of heavy ions. Shocks which contain heavy ions that are a minority constituent by number but which carry most of the energy density in the upstream medium put 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) E ^-2, where N(E)dE is the number of particles with energy between E and E+dE. Synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front, provides the mechanism of thermalization and non-thermal particle acceleration. The maximum energy achievable by the pairs is _± m _± c ² = m _i c ² ₁/Z _i, where ₁ is the Lorentz factor of the upstream flow and Z _i is the atomic number of the ions. The shock's spatial structure contains a series of overshoots in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. These overshoots provide a new interpretation of the structure of the inner regions of the Crab Nebula, in particular of the wisps, surface brightness enhancements near the pulsar. The wisps appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure ₁ 4 × 10⁶, the upstream magnetic field B ₁ to be 3 × 10^-5 Gauss, as well as to show that the total ion flow is 3 × 10³⁴ elementary charges/sec, in good agreement with the total current flow predicted by the early Goldreich and Julian (1969) model. The total pair outflow is shown to be about 5 × 10³⁷ pairs per second, in good agreement with the particle flux required to explain the nebular X—ray source.The energetics of particle acceleration within the magnetospheres of rotation powered pulsars and the consequences for pulsed gamma ray emission are also briefly discussed. The gamma ray luminosity above 100 MeV is shown to scale in proportion to _R ^1/2 , as is in accord with some of the simplest ideas about polar cap models. Models based on acceleration in the outer magnetosphere are also briefly discussed.     

The importance of anisotropic interstellar turbulence and molecular-cloud magnetic mirrors for galactic cosmic-ray propagation

Recent studies suggest that when magnetohydrodynamic (MHD) turbulence is excited by stirring a plasma at large scales, the cascade of energy from large to small scales is anisotropic, in the sense that small-scale fluctuations satisfy the inequality k k , where k and k are, respectively, the components of a fluctuations wave vector and to the background magnetic field. Such anisotropic fluctuations are very inefficient at scattering cosmic rays. Results based on the quasilinear approximation for scattering of cosmic rays by anisotropic MHD turbulence are presented and explained. The important role played by molecular-cloud magnetic mirrors in confining and isotropizing cosmic rays when scattering is weak is also discussed.     

Ariel-6 medium energy spectral observations of active galaxies

X-ray spectra of the BL Lac type object Mkn 421 and several Seyfert type 1 galaxies; IIIZw2, MCG8-11-11 and NGC 4151, have been obtained using the Leicester University instrument on board the Ariel-6 satellite. The Mkn 421 spectrum is best represented by two powerlaw components, the soft component having 3.4 whilst the hard flux has 1.0. In MCG8-11-11 there is clear evidence for spectral variability between our observation in late 1979 and that of HEAO-1/A2 in 1977. The Ariel-6 spectrum of MCG8-11-11 can be fitted by a powerlaw of index 2.1 together with an iron line at 6.2 keV with an equivalent width of 1.6 keV. The first X-ray spectrum of IIIZw2 is also presented, fitting with a powerlaw we find an index of 1.7. With the exception of NGC 4151 there is no evidence for a significant column of cool material along the line of sight.