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.     

Determination of the characteristics of the gamma-ray telescope Gamma-1

The Gamma-1 telescope has been developed through a collaboration of scientists in the USSR and France in order to conduct -ray astronomical observations within the energy range from 50 to 5000 MeV. The major characteristics of the telescope were established by Monte-Carlo simulations and calibrations made with the aid of electron and tagged -ray beams produced by an accelerator, and these have been found to be as follows: the effective area for photons coming along the instrument's axis varies from about 50 cm² at E = 50 MeV to approximately 230 cm² at E 300 MeV; the angular resolution (half opening of the cone embracing 68% events) is equal to 2.7° at E = 100 MeV, and 1.8° at E = 300 MeV; the energy resolution (FWHM) varies from 70% to 35% as the energy of the detected photons increases from 100 to 550 MeV; the telescope's field-of-view at the half-sensitivity level is 300–450 square degrees depending upon the spectrum of the detected radiation, and the event selection logic. Proceeding from the thus obtained characteristics it is demonstrated that a point source producing a photon flux J (E 100 MeV) = 3 × 10^-7 cm^-2 s^-1, can be detected with a 5 significance by observing it during 10⁶ s at the level of the Cygnus background, and a source having intensity J (E 100 MeV) = 10^-6 cm^-2 s^-1 can be detected to within a mean square positional accuracy of about 15.     

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.     

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.     

Plasma waves at the dayside magnetopause

Experimental investigations of plasma waves at the magnetopause, including recent results from the AMPTE/IRM satellite, show that both E and B fluctuations typically have a featureless spectrum which monotonically decreases with frequency; integrated rms amplitudes are typically a few mV m^-1 for E and 10 nT for B, though in particular E can be as much as an order of magnitude larger in exceptional cases. Surveys show a lack of correlation between wave parameters and the magnetopause parameters. Under the assumption that crossing the diffusion region would give a pronounced signature in the waves, the survey data allow an upper limit to be placed on the latitudinal extent of the diffusion region, which is about 1000 km — implying that it is not surprising that the wave data surveys have so far failed to detect it. The observed wave turbulence levels have been used to estimate diffusion coefficients under different assumptions for the wave mode, but the resulting diffusion coefficient is always too small to explain either reconnection or boundary layer formation. Recent work of Galeev et al. (1986) indicates that the dominant diffusion process may be magnetic field migration, which is a macroscopic process involving the interaction of tearing mode islands. Assuming this mode to be present at the observed level of B, a particle diffusion coefficient of nearly 10⁹ m² s^-1 is obtained. Another macroscopic diffusive process which could occur at the magnetopause is stochastic E × B scattering, which also implies a diffusion coefficient the order of 10⁹ m² s^-1 if the observed E spectrum is assumed to be a turbulent cascade consisting of convective cells.     

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     

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.     

The effect of rotation on the hydrodynamics of stellar collapse

Summary Using values of d, _min, and _max that Van Riper (1978) has found most promising for a hydrodynamic envelope ejection, we have shown that even a small amount of rotation in the initial core can stop its collapse before nuclear densities are reached. We expected _i > 0.02 to produce significant deviations from a spherically symmetric collapse, but have found that _i as much as ten times smaller than this will not allow the core to reach densities as high as in the spherical collapse. In no case, however, does the core flatten very much, nor does the value of become very large. Low final 's preclude the formation of an axisymmetric torus. They also indicate that deformation of an iron core into a triaxial configuration or fragmentation of the core during its collapse is an extremely unlikely event. (Note: Classically, must exceed 0.27 before a dynamic instability to non-axisymmetric perturbations is encountered.)The small degree of flattening of the core also suggests that the reduced moment of inertia I of the core will always be relatively small in magnitude and hence that the third time derivative of I, which is proportional to the energy emitted in gravity wave radiation, will not be very significant. Numerically calculated estimates of I- during some of these model evolutions supports this suspicion. If the _min and used here are found to be realistic values after the detailed physics of the core collapse is well understood, it is clear that gravitational radiation from a core collapse will be difficult to measure.Finally, we should point out that it is the relatively large values of Y_min (near 4/3) combined with values of d near unity that (a) prevented the core from flattening significantly in these models and (b) prevented the core from reaching high configurations. If realistic values of either one (or both) of these parameters are found to be much smaller in more complete models of the core collapse, then the core will have to become flatter (and denser) before pressure gradients will support it along the rotation axis. All of the conclusions drawn here would be modified accordingly under those circumstances. It should also be noted that in general relativistic models, the critical for spherical collapse is somewhat larger than 4/3 (Van Riper, 1979). Therefore, we predict that when fully general relativistic core collapses are performed including rotation, a given choice of _min and _i will produce a slightly flatter and slightly denser core than the corresponding model that has been presented here.     

Presupernova models and supernovae

Present status of the theories for presupernova evolution and triggering mechanisms of supernova explosions are summarized and discussed from the standpoint of the theory of stellar structure and evolution. It is not intended to collect every detail of numerical results thus far obtained, but to extract physically clear-cut understanding from complexities of the numerical stellar models. For this purpose the evolution of stellar cores is discussed in a generalized fashion. The following types of the supernova explosions are discussed. The carbon deflagration supernova of intermediate mass star which results in the total disruption of the star. Massive star evolves into a supernova triggered by photo-dissociation of iron nuclei which results in a formation of a neutron star or a black hole depending on its mass. These two are typical types of the sueprnovae. Between them there remains a range of mass for which collapse of the stellar core is triggered by electron captures, which has been recently shown to leave a neutron star despite oxygen deflagration competing with the electron captures. Also discussed are combustion and detonation of helium or carbon which take place in accreting white dwarfs, and the collapse which is triggered by electron-pair creation in very massive stars.Appendix: Notations A mass number of atomic nucleus - B _v(a, b) incomplete beta function - c _p specific heat at constant pressure - c _p sound velocity - c(sub) center of the star - E _e mean energy of an electron captured by nucleus - E _n nuclear energy release from unit mass of the nuclear fuel specified by n - E _thr threshold energy (9.3) - E _thr,0 energy difference between the ground states of daughter nucleus and parent nucleus (9.1) - E energy of gamma ray emitted from daughter nucleus (9.1) - E _v mean energy of a neutrino emitted by electron capture (9.1) - f flatness parameter (2.17) - g local gravitational acceleration (2.16) - H atomic mass unit - H _p scale height of pressure (2.22) - H (sub) hydrogen-burning shell - k Boltzmann constant - l mixing length of convection - L _cr(M _r ) local Eddington's critical luminosity (4.3) - L _n integrated nuclear energy generation rate by nuclear fuel specified by n - L _v neutrino luminosity - L _{v, cr}(M _r ) local Eddington's critical neutrino luminosity (11.2) - M (current) mass of a star - m _M core mass contained interior to the carbon-burning shell - M _Ch Chandrasekhar's limiting mass (9.6) - M _H core mass contained interior to the hydrogen-burning shell - M _He core mass contained interior to the helium-burning shell - M _ms mass of a star at its zero-age min-sequence - M _O core mass contained interior to the oxygen-burning shell - M _r mass contained interior to a shell at r - M _Si core mass contained interior to the silicon-burning shell - M _WD mass of white dwarf (7.1) - M ₀ normalization factor to the non-dimensional mass (3.3) - M ₁ core mass (3.6) - N polytropic index between pressure and density (2.3) - n polytropic index between pressure and temperature (10.1) - N _A Avogadro number - N _ad adiabatic polytropic index - N _e number of electrons in unit mass of matter - NSE nuclear statistical equilibrium - P pressure - ph (sub) photosphere - Q _e mass fraction of the envelope exterior of the shell e (2.14) - R stellar radius - r radial distance of a shell - r ₀ normalization factor to the non-dimensional radius (3.2) - s specific entropy - S _i specific entropy of ions - T temperature - U homology invariant defined by (2.1) - u _gas specific internal energy of gas - u _rad energy of the radiation field per volume in which unit mass of gas is contained (6.4) - V homology invariant defined by (2.2) - _def velocity of deflagration front (6.10) - X concentration by weight of hydrogen - Y concentration by weight of helium - Y _e mole number of electrons in one gram of matter (9.7) - Y _v mole number of neutrinos in one gram of matter - Z concentration by weight of the elements other than hydrogen and helium - z shock strength (6.6) - 1 (sub) usually denotes the core edge (2.13) - ratio of the mixing length to the scale height of pressure (l/H _p ) - ratio of gas pressure to the total pressure - ratio of the specific heats - gD locus of singularity in U-V plane (2.5) - M(H _p ) mass contained within unit scale height of pressure (4.4) - _ec energy rate by electron captures (9.5) - _n nuclear energy generation rate by the nuclear fuel specified by n - _v neutrino loss rate - L _v ^(D) neutrino loss rate excluding the neutrinos from the electron captures (9.4) - non-dimensional density (3.1) - P/, not the non-dimensional temperature (2.7) - _W Weinberg's angle (5.8) - opacity - _v neutrino opacity (11.2) - describes the effect of electron degeneracy in equation of state (2.19) - _ec rate of electron capture - mean molecular weight - _e mean molecular weight of electrons - _e chemical potential of an electron excluding the rest mass (8.1) - _i mean molecular weight of ions - non-dimensional radius (3.1) - non-dimensional pressure (3.1) - matter density - _cr ^GR critical density above which the general relativistic instability sets in - _cr critical density for reimplosion of the core by beta processes (Section 5) - _ign density at the ignition - _nse density above which the deflagrated matter results in NSE composition - _e non-dimensional entropy of electron-per one electron in units of k(9.2) - _ff timescale of free fall (6.2) - _h (H _p ) timescale of heat transport over unit scale height of pressure (4.4) - _n nuclear timescale for a change in temperature (6.1) - non-dimensional mass (3.1) - _e chemical potential of an electron in units of kT (8.1)     

Electric gaps in pulsar magnetospheres

Two contrasting models of the magnetosphere are studied: a classical, quasi-mhd model, in which it is the componentE that accelerates electrons to relativistic energies; and a model with e⁺-e^– production, in which the componentE is locally crucial, both in an acceleration domain near the star and in a dissipation domain beyond the light-cylinder.     

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.     

Recent results on the parameters of the interstellar helium from the ULYSSES/GAS experiment

Velocity and direction of the flow of the interstellar helium and its temperature and density have been determined from the measurements of the ULYSSES/GAS experiment for two different epochs: during the in-ecliptic path of ULYSSES, representing solar maximum conditions, and during the south to the north pole transition (11/94-6/95), close to the solar minimum conditions. Within the improved error bars the values are consistent with results published earlier.The determination of the density n of the interstellar helium at the heliospheric boundary from observations in the inner solar system requires knowledge about the loss processes experienced by the particles on their way to the observer. The simultaneous observation of the helium particles arriving on direct and indirect orbits at the observer provides a tool to directly determine the effects of the loss processes assumed to be predominantly photoionization and — for particles travelling close to the Sun — electron impact ionization by high-energy solar wind electrons.Such observations were obtained with the ULYSSES/GAS instrument in February 1995, before the spaceprobe passed its perihelion. From these measurements values for the loss rates and the interstellar density could be derived. Assuming photoionization to be the only loss process reasonable fits to the observations were obtained for an ionization rate = 1.1 · 10^–7 s^–1 and a density n 1.7 · 10^–2 cm^–3. Including, in addition, electron impact ionization, a photoionization = 0.6 · 10^–7 s^–1 was sufficient to fit both observations, resulting in a density n 1.4 · 10^–2 cm^–3.On leave from Space Research Centre, Warsaw, Poland.     

Some problems of gamma-astronomy

The modern state of gamma-ray astronomy is reviewed, the paper being mainly devoted to the theoretical models that describe generation of gamma-ray emission under astrophysical conditions. Basic information on the processes of generation and absorption of gamma-rays, as well as the results of observations for various gamma-ray photon energies are reported.In the region of soft gamma-ray emission (i.e., for energies less than tens of MeV), where emission in gamma-ray lines dominates, we also discuss the nature of gamma-ray bursts, the origin of gamma-ray emission from the galactic centre, etc.Discrete sources and, in particular, the mysterious source Cyg X-3 are discussed in the region of very high (E > 10¹² eV) and ultra-high (E > 10¹⁵ eV) energy gamma-ray emission.A larger portion of the review is devoted to the analysis of cosmic-ray origin on the basis of the available gamma-ray data in the region from several tens of MeV to several GeV. The peculiarity of this energy range is, in particular, in the fact that the diffuse galactic emission was observed mainly there. We also discuss the problem of determination of the cosmic-ray density gradient from the gamma-ray data.The origin of high-latitude gamma-ray emission, the problem of galactic gamma-ray halo, etc., are discussed.The theoretical models explaining the nature of unidentified gamma-ray sources, as well as the results of measurements and theoretical estimations of a gamma-ray flux from SN1987A are analysed.List of Notations m electron mass, m = 9.108 × 10^–28 g, - M proton mass, M = 1.672 × 10^–24 g, - e electron charge, e = 4.803 × 10^–10 CGS - c velocity of light, c = 2.9979 × 10¹⁰ cm s^–1, - k Boltzmann constant, k = 1.380 × 10^–16 erg grad^–1, - e electron - p proton - gamma-ray photon - p antiproton - ^{0 0}-meson - -lepton - e ⁺ positron - r, , x radio-frequency, gamma-ray, and X-ray emission bands - E total energy of a particle - E _k kinetic energy - p particle momentum - spectral index for particles - spectral index for emission - n particle density (concentration) - H magnetic field strength - T temperature - _ph energy of low-energy photons - emission frequency - r _H Larmor radius of relativistic particles - k wave number - , z cylindric coordinates, in this case the coordinate (radius) along the galactic disk, z perpendicular to the galactic disk - M solar mass, M = 1.99 × 10³³ g.     

Pulsations of the R Coronae Borealis stars

The radial pulsations of very luminous, low-mass models (L/M 10⁴, solar units), which are possible representatives of the R CrB stars, have been examined. These pulsations are extremely nonadiabatic. We find that there are in some cases at least one extra (strange) mode which makes interpretation difficult. The blue instability edges are also peculiar, in that there is an abrupt excursion of the blue edge to the blue for L/M sufficiently large. The range of periods of the model encompasses observed periods of the Cepheid-like pulsations of actual R CrB stars.     

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.     

Relations between ISM inside and outside the heliosphere

Thanks to remarkable new tools, such as the Goddard High Resolution Spectrograph (GHRS) on board the HST and the EUVE spectrometer on the interstellar side, and Ulysses particle detectors on the heliospheric side, it is possible now to begin to compare abundances and physical properties of the interstellar matter outside the heliosphere (from absorption features in the stellar spectra), and inside the heliosphere (from in situ or remote detection of the interstellar neutrals or their derivatives, the pick-up ions or the Anomalous Cosmic Rays detected by the two Voyager spacecraft).Ground-based and UV spectra of nearby stars show that the Sun is located between two volumes of gas of different heliocentric velocities V and temperatures T (see also Linsky et al, this issue). One of these clouds has the same velocity (V= 25.6 km s^–1 from = 255 and =8) and temperature (6700 K) as the heliospheric helium of interstellar origin probed by Ulysses, and is certainly surrounding our star (and then the Local Interstellar Cloud or LIC). This Identification allows comparisons between interstellar constituents on both sides of the heliospheric interface.Ly-alpha background data (absorption cell and recent HST-GHRS spectra) suggest that the heliospheric neutral H velocity is smaller by 5–6 km s^–1 than the local cloud velocity, and therefore that H is decelerated at its entrance into the heliosphere, in agreement with interaction models between the heliosphere and the ISM which include the coupling with the plasma. This is in favor of a non negligible electron density (at least 0.05 cm³). There are other indications of a rather large ionization of the ambient ISM, such as the ionization equilibrium of interstellar magnesium and of sodium. However the resulting range for the plasma density is still broad.The heliospheric neutral hydrogen number density (0.08–0.16 cm^–3) is now less precisely determined than the helium density (0.013–0.017 cm^–3, see Gloeckler, Witte et al, Mobius, this issue). The comparison between the neutral hydrogen to neutral helium ratios in the ISM (recent EUVE findings) and in the heliosphere, suggests that 15 to 70% of H does not enter the heliosphere. The comparison between the interstellar oxygen relative abundance (with respect to H and He) in the ISM and the heliospheric abundance deduced from pick-up ions is also in favor of some filtration, and thus of a non-negligible ionization.For a significant ISM plasma density, one expects a Hydrogen wall to be present as an intermediate state of the interstellar H around the interface between inside and outside. Since 1993, the two UVS instruments on board Voyager 1 and 2 indeed reveal clearly the existence of an additional Ly-alpha emission, probably due to a combination of light from the compressed H wall, and from a galactic source. On the other hand, the decelerated and heated neutral hydrogen of this H wall has recently been detected in absorption in the spectra of nearby stars (see Linsky, this issue).     

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     

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.     

Large scale structure of the universe

The current situation with the cosmological model and fundamental constants is briefly reviewed. Here, we concentrate on evolutionary effects of large-scale structure formation, in particular, the relationship with the quasar distribution and dynamics is discussed. We argue that groups of bright quasars with few or more than dozen of members within regions l _LS(100–150)h ^–1 Mpc found atz<2 may belong to concentrations of young rich clusters of galaxies, and thus be distant Great Attractors like the local GA or the Shapley concentration. These early large-scale galactic structures (i) provide a natural way to bias the distribution of Abell clusters, and (ii) suggest that the spectrum of primordial density perturbations is nearly flat on scales encompassing both the cluster and GAs,l=k ^–1(10,100)h ^–1 Mpc: _k ² k ³ P(k) k , =1 _–0.4 ^+0.6 , whereP(k) is the power spectrum of density perturbations.     

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.