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.     

The synthesis of helium and CNO isotopes in massive stars

We present helium and CNO isotopic yields for massive mass-losing stars in the initial mass range 15M M _i 50M . We investigate their dependence on assumptions about mass loss rates, internal mixing processes, and metallicity, and specify the contributions from stellar winds and from supernova ejecta.     

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.     

Anomalous Cosmic Rays and Solar Modulation

We review aspects of anomalous cosmic rays (ACRs) that bear on the solar modulation of energetic particles in the heliosphere. We show that the latitudinal and radial gradients of these particles exhibit a 22-year periodicity in concert with the reversal of the Sun's magnetic field. The power-law index of the low energy portion of the energy spectrum of ACRs at the shock in 1996 appears to be -1.3, suggesting that the strength of the solar wind termination shock at the helioequatorial plane is relatively weak, with s 2.8. The rigidity dependence of the perpendicular interplanetary mean free path in the outer heliosphere for particles with rigidities between 0.2 and 0.7 GV varies approximately as R2, where R is particle rigidity. There is evidence that ACR oxygen is primarily multiply charged above 20 MeV/nuc and primarily singly-charged below 16 MeV/nuc. The location of the termination shock was at 65 AU in 1987 and 85 AU in 1994.     

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.     

Optical identification of X-ray sources in the Einstein observatory medium and deep surveys

Methods are discussed for establishing the optical identification of X ray sources in the medium and deep X-ray surveys of the Einstein Observatory. Of the 63 X-ray sources with a statistical significance of 5 in the medium survey (Maccacaro et al. 1981), optical identification work is summarized for 51, of which identifications have been made with 30 active galactic nuclei. The optical properties of some of these X-ray selected objects are briefly discussed.The Einstein deep survey of Pavo (Griffiths et al. 1981) is used to illustrate the problems and methods used for securing optical identifications for X-ray sources in the deep survey fields. Identifications have been made with 4 QSOs at the bright end of the optical candidate distribution (together with 3 G stars) and it is shown that a further 7 fainter objects are also likely to be QSOs.     

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.     

Characteristics and interpretation of the photometric variability of eta carinae and its nebula

A short review is given on the history of the peculiar variable object Car and on a number of relevant references describing and discussing its physical characteristics and behaviour, based on different types of observational techniques. The star is known to be variable since the 17th century. The excessive mass loss to which it was subject during the 19th century is now visible as an ellipsoidal reflection nebula of 15 diameter: the so-called homunculus. The remainder of the paper is spent on different kinds of problems partly based on the results of a decade of photometric monitoring in the VBLUW photometric system of Walraven. Foreground reddening and reddening by dust in the homunculus are determined and amount to E(B - V) _J = 0 50 and < 6, respectively. Scanning of the homunculus revealed an estimate for the photometric characteristics of the central object, which presumably consists of a massive hot star surrounded by a cooler gas envelope. The total luminosity is derived using fluxes of various sources in the wavelength region 0.15 < < 175 n resulting in M _bol = - 12 3 ± 0 2. The total observed flux corrected for foreground extinction corresponds to a star with R 96 R if T _eff 30 000 K. The mass may be near 150 M . The excess luminosity in 1843, when the star was presumably bolometrically at least 2 5 brighter than at present, may have been caused by envelope-energized pulsations when the star's luminosity was close to its Eddington limit. The temperature should then have been 50 000 K. The mass loss rate, during the excess luminosity phase lasting 30 yr, is estimated to amount to M 4 × 10^-3 M yr^-1. At present the mass loss may be M 10^-4M yr^-1. Since the homunculus is mainly built up from material expelled in the 30 yr interval (from 1830 to 1860), its total mass amounts to M _hom 0.15 M . The historical observations of the colours of Car and a comparison with the characteristics of S Dor type stars, suggest that Car was subject to a number of S Dor type phases similar to those of P Cyg (in the 17th century), S Dor and others. A satisfactory explanation is found for the complete historical light curve. The decrease in light after the 1843 maximum by 9 ^m , presumably consists of a fading of the luminosity excess and the S Dor effect by 2 5 and 3^m, respectively, and a 3 5 extinction by circumstellar dust. The small amplitude light variations which Car showed during the last decade, were studied with the aid of the variations of the Balmer jump. They are presumably caused by temperature variations of the central star.Percy and Welch (1983) (Publ. Astron. Soc. Pacific 95, 491) have observed P Cyg on a number of nights in 1982 and found for the photometric variations a time scale of 30 to 50 days and an amplitude of 0 _. ^m 15.Based partly on observations collected at the ESO, La Silla, Chile.     

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     

Collapse,equilibrium, and fragmentation of rotating,adiabatic clouds

Numerical calculations of the collapse of adiabatic clouds from uniform density and rotation initial conditions show that when restricted to axisymmetry, the clouds form either near-equilibrium spheroids or rings. Rings form in the collapse of low thermal energy clouds and have = T/ ¦W¦ 0.43. When the axisymmetric constraint is removed and an initial m=2 density variation is introduced, clouds either collapse to form near-equilibrium ellipsoids or else fragment into binary systems through a bar phase. Ellipsoids form in the collapse of high thermal energy clouds and have 3 0.27. The results are consistent with the critical values of for instabilities in Maclaurin spheroids, and suggest that protostellar clouds may undergo a dynamic fragmentation in the nonisothermal collapse regime.National Academy of Sciences — National Research Council Resident Research Associate.     

The population of massive stars in R136 from HST/FOC UV observations

New ultraviolet (1300 A, 3400 A),HST FOC observations have been used to derive the UV color-magnitude diagram (CMD) of R136, with the main scientific goal of studying the upper end of the stellar mass function at ultraviolet wavelengths where the color degeneracy encountered in visual CMDs is less severe. The CMD has been compared to a set of theoretical isochrones, which have been computed using the latest generation of evolutionary models and model atmospheres for early type stars. Wolf-Rayet stars are included. Comparison of theTheoretical andobserved CMD suggests that there are no stars brighter than M₁₃₀–11. We use the observed main sequence turn-off and the known spectroscopic properties of the stellar population to derive constraints on the most probable age of R136. The presence of WNL stars and the lack of red supergiants suggests a most likely age of 3±1 Myr. A theoretical isochrone of 3±1 Myr is consistent with the observed stellar content of R136 if the most massive stars have initial masses around 50 M.Bases on Observations with the NASA/ESA Hubble Space Telescope, obtained at the STScI, which is operated by AURA, Inc., under NASA contract NAS5-26555.Astrophysics Division, Space Science Department, ESA     

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 evolution of massive stars to explosion

We review the possible evolutionary paths from massive stars to explosive endpoints as various types of supernovae associated with Population I and hence with massive stars: Type II-P, Type II-L, Type Ib, Type Ic, and the hybrid events SN 1987K and SN 1993J. We identify SN 1954A as another hybrid event from the evidence for both H and He in its spectrum with velocities nearly the same as SN 1983J. Evidence for ejected⁵⁶Ni mass of 0.07 M suggests that SN II-P underwent standard iron core collapse, not collapse of an O–Ne–Mg core nor thermonuclear explosion of a C–O core. Most SN II-P presumably arise in single stars or wide binaries of 10–20 M. There may be indirect evidence for duplicity in some cases in the form of strong Ba II lines, such as characterized SN 1987A. SN II-L are recognizably distinct from typical SN II-P and must undergo a significantly different evolution. Despite indications that SN II-L have small envelopes that may be helium enriched, they are also distinct from events like SN 1993J that must have yet again a different evolution. The SN II-L that share a common Luminosity seem to have ejected a small nickel mass and hence may come from stars with O–Ne–Mg cores. The amount of nickel ejected by the exceptionally bright events, SN 1980K and SN 1979C, remains controversial. SN Ib require the complete loss of the H envelope, either to a binary companion or to a wind. The few identified have relatively large ejecta masses. It is not clear what evolutionary processes distinguish SN Ib's evolving in binary systems from hybrid events that retain some H in the envelope. SN Ic events are both H and He deficient. Binary models that can account for transfer of an extended helium envelope from low mass helium cores, 2 to 4 M, imply C–O core masses that are roughly consistent with that deduced from the ejecta mass plus a neutron star, 2 to 3 M. It is possible that the hybrid events are the result of Roche lobe overflow and that the pure events, SN Ib or SN Ic, result from common envelope evolution.     

Fourier parameters of heliospheric current sheet and their significance

If the path of the neutral line on the coronal source surface is expressible as a singlevalued function (colatitude vs longitude ), then Fourier analysis of ctn with respect to leads to a simple algorithm for realistically mapping the neutral line outward to model the heliospheric current sheet (HCS) at distancesr1 AU. To be compatible with MHD, the source surface used for this mapping should be prolate (aligned with dipole axis) rather than spherical. Orientation of the Sun's magnetic-dipole moment is indicated by them=1 Fourier amplitude (a ₁ sin +b ₁ cos ) of ctn on the source surface. Physical features (including the neutral line) on a prolate source surface intrinsically map to lower dipole latitudes atr1 AU in the heliosphere, and Ulysses observations of a unipolar field at latitudes beyond 30°S (when the neutral line on the source surface still reached 39°S) confirm the expected geometry.     

Non-BBN Constraints on the Key Cosmological Parameters

Since the baryon-to-photon ratio 10 is in some doubt at present, we ignore the constraints on 10 from big bang nucleosynthesis (BBN) and fit the three key cosmological parameters (h, M, 10) to four other observational constraints: Hubble parameter (ho), age of the universe (to), cluster gas (baryon) fraction (fo fGh3/2), and effective shape parameter (o). We consider open and flat CDM models and flat CDM models, testing goodness of fit and drawing confidence regions by the 2 method. CDM models with M = 1 (SCDM models) are accepted only because we allow a large error on ho, permitting h < 0.5. Open CDM models are accepted only for _M 0.4. CDM models give similar results. In all of these models, large 10 ( 6) is favored strongly over small 10 ( 2), supporting reports of low deuterium abundances on some QSO lines of sight, and suggesting that observational determinations of primordial 4He may be contaminated by systematic errors. Only if we drop the crucial o constraint are much lower values of M and 10 permitted.     

Extragalactic Abundances of Hydrogen,Deuterium and Helium New Steps,Missteps and Next Steps

Estimates of the deuterium abundance in quasar absorbers are reviewed, including a brief account of incorrect claims published by the author and a brief review of the problem of hydrogen contamination. It is concluded that the primordial abundance may be universal with a value (D/H)P 10-4, within about a factor of two, corresponding to _Bh _0.7 ² 0.0fs2 or 10 2.7 in the Standard Big Bang. This agrees with current limits on primordial helium, YP 0.243, which are shown to be surprisingly insensitive to models of stellar enrichment. It also agrees with a tabulated sum of the total density of baryons in observed components. Much lower primordial deuterium ( 2 × 10^-5) is also possible but disagrees with currently estimated helium abundances; the larger baryon density in this case fits better with current models of the Lyman- forest but requires the bulk of the baryons to be in some currently uncounted form.     

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 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.     

LISM structure — Fragmented superbubble shell?

Small scale structure in local interstellar matter (LISM) is considered. Overall morphology of the local cloud complex is inferred from Ca II absorption lines and observations of H I in white dwarf stars. Clouds with column densities ranging from 2–100 × 10¹⁷ cm^–2 are found within 20 pc of the Sun. Cold (50 K) dense (10⁵ cm^–3) small (5–10 au) clouds could be embedded and currently undetected in the upwind gas. The Sun appears to be embedded in a filament of gas with thickness 0.7 pc, and cross-wise column density 2 × 10¹⁷ cm^–2. The local magnetic field direction is parallel to the filament, suggesting that the physical process causing the filamentation is MHD related. Enhanced abundances of refractory elements and LISM kinematics indicate outflowing gas from the Scorpius-Centaurus Association. The local flow vector and Sco data are consistent with a 4,000,000 year old superbubble shell at –22 km s^–1, which is a shock front passing through preshock gas at –12 km s^–1, and yielding cooled postshock gas at –26 km s^–1in the upwind direction. A preshock magnetic field strength of 1.6 G, and postshock field strength of 5.2 G embedded in the superbubble shell, are consistent with the data.Abbreviations LISM Local ISM - SIC Surrounding Interstellar Cloud - LIC Local Interstellar Cloud     

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.