Physical and chemical characteristics of the ISM inside and outside the heliosphere

This paper summarizes some of the discussions of working group 8–9 during the ISSI Conference on The Heliosphere in the Local Interstellar Medium. Because the subject of these working groups has become significantly broader during the last ten years, we have selected three topics for which recent observations have modified and improved our knowledge of the heliosphere and the surrounding interstellar medium. These topics are the number densities and ISM ionization states of hydrogen and helium, the newly discovered hot gas from the H wall seen in absorption, and the comparison between ISM and heliospheric minor element abundances. Papers from this volume in which more details on these topics can be found are quoted throughout the report.     

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     

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.     

The local interstellar medium

The local interstellar medium can be probed in different ways: by analyzing low energy X-ray data in the range 0.1–0.4 keV, where the radiation is absorbed by the interstellar gas at column densities in excess of about 10²⁰ cm^-2 — and can therefore be regarded as local, by determining the absorption of stellar emission spectra from nearby stars along their lines of sight by intervening gas and by directin situ measurements of those components which penetrate the heliosphere sufficiently far, provided they can be distinguished from interplanetary material. The current status of these different investigations gives the following picture: the solar system is surrounded by a bubble of hot gas (density 0.005cm^-3, temperature 10⁶ K) out to several tens of parsecs. More locally it is embedded in a small warm cloud of density 0.07cm^-3, temperature 7000 K, column density 5 × 10¹⁷ cm^-2 — which gives a mass of about 0.1M . The transition to the heliosphere is governed by solar UV ionization, snowploughing of the interstellar gas by the outwardly expanding solar wind and the bow shock. The heliosphere is the region inside the solar wind terminal shock. Classically it would be regarded as not yet affected by (or aware of) the obstacle ahead. Practically, the existence of the interstellar medium makes itself felt even far inside the heliosphere by the penetration of neutral gas, dust, plasma waves, shock accelerated particles and cosmic rays. These are the local probes of the interstellar medium.     

Modelling of the interstellar hydrogen distribution in the heliosphere

The detailed knowledge of the distribution of neutral interstellar hydrogen in the interplanetary space is necessary for a reliable interpretation of optical and H⁺ pickup ions observations. In the paper, we review the status of the modelling efforts with the emphasis on recent improvements in that field. We discuss in particular the role of the nonstationary, solar cycle-related effects and the consequences of hydrogen filtration through the heliospheric interface region for its distribution in the inner Solar System. We demonstrate also that the use of the simple cold model, neglecting the thermal character of the hydrogen gas (T 8000 K), is generally incorrect for the whole region of the inner heliosphere (R < 5 AU) since it leads to a substantial underestimation of the local hydrogen density and thus influences the derivation of the H properties in the outer heliosphere/LISM. Referring to recent Ulysses measurements, we point out also the need to consider in the modelling the effects of the latitudinal asymmetry of the ionization rate.     

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.     

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.     

The abundance of atomic 1H, 4He and 3He in the local interstellar cloud from pickup ion observations with SWICS on Ulysses

Pickup ions measured deep inside the heliosphere open a new way to determine the absolute atomic density of a number of elements and isotopes in the local interstellar cloud (LIC). We derive the atomic abundance of hydrogen and the two isotopes of helium from the velocity and spatial distributions of interstellar pickup protons and ionized helium measured with the Solar Wind Ion Composition Spectrometer (SWICS) on the Ulysses spacecraft between 2 and 5 AU. The atomic hydrogen density near the termination shock derived from interstellar pickup ion measurements is 0.115±0.025 cm^–3 and the atomic H/He ratio from these observations is found to be 7.7 ± 1.3 in the outer heliosphere. Comparing this value with the standard universal H/He ratio of 10 we conclude that filtration of hydrogen is small and that the ionization fraction of hydrogen in the LIC is low.     

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.     

Voyager Observations of the Magnetic Field in the Distant Heliosphere

The latitudinal structure of the heliospheric magnetic field during much of the solar cycle is determined by a "sector zone", in which both positive and negative magnetic polarities are observed, and by the unipolar regions above and below the sector zone. Distinct corotating streams and interactions regions are found primarily in the sector zone during the declining phase of the solar cycle. Within a few AU, the streams and interaction regions are distinct and are related to solar features. A restructuring of the solar wind occurs between 1 AU and 15 AU, in which the isolated streams, interaction regions and shocks merge to form compound streams and merged interaction regions ("MIRs"). Memory of the source conditions is lost in this process. In the region between 30 AU and the termination shock (the "distant heliosphere"), the pressure of interstellar pickup protons dominates that of the magnetic field and solar wind particles and largely controls the dynamical processes. During 1983 and 1994, corotating streams and corotating interaction regions were observed at 1 AU. Merged interaction regions were observed at 15 AU in 1983, but not at 45 AU during 1994. This result suggests a further restructuring of the solar wind in the distant heliosphere, but variations from one solar cycle to the next might also contribute to the result. Approaching solar minimum in 1996, the latitudinal extent of the sector zone decreased, and Voyager 2 gradually entered the unipolar region below it. The speed was lower in the sector zone than below it. At Voyagers 1 and 2, the change in cosmic ray intensity is related to the magnetic field strength during each year from 1983 through 1996. The magnetic field strength has a multifractal distribution throughout the heliosphere. This fundamental symmetry of the heliosphere has not been incorporated explicitly in cosmic ray propagation models.     

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.     

Cosmic ray modulation in the 3-D heliosphere

The basic physical processes that lead to the long-term modulation of cosmic rays by the solar wind have been known for many years. However our knowledge of the structure of the heliosphere, which determines which processes are most important for the modulation, and of the variation of this structure with time and solar activity level is still incomplete. Study of the modulation provides a tool for probing the scale and structure of the heliosphere. While the Pioneer and Voyager spacecraft are surveying the radial structure and extent of the heliosphere at modest heliographic latitudes, theUlysses mission is the first to undertake a nearly complete scan of the latitudinal structure of the modulated cosmic ray intensity in the inner heliosphere (R<5.4 AU).Ulysses will reach latitudes of 80°S in September 1994 and 80°N in July 1995 during the approach to minimum activity in the 11 year solar cycle. We present a first report of measurements extending to latitudes of 52°S, which show surprisingly little latitudinal effect in the modulated intensities and suggest that at this time modulation in the inner heliosphere may be much more spherically symmetric than had generally been believed based upon models and previous observations.     

Detection limits in the far infrared

The advent of far infrared arrays will change fundamentally the means of analyzing observations in this spectral region. Sources much fainter than traditional confusion limits will be extracted from images by using computer algorithms similar to CLEAN or DAOPHOT. We have conducted numerical experiments to evaluate these techniques and show that they will permit long integrations (10,000 sec at 60 m, 200 sec at 100 m) to achieve nearly photon-background-limited performance and hence very deep detection limits. The dominant noise sources—photon noise, confusion by distant galaxies, and confusion by IR cirrus — scale with nearly the same power of the telescope aperture. As a result, the integration times required to reach confusion limits are nearly aperture-independent.     

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.     

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.     

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     

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     

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.     

Latitude and solar-cycle dependence of radial IMF intensity

I describe a simple procedure for extrapolating the observed solar magnetic field into the heliosphere, which averages the asymptotic fields computed using the standard source surface and current sheet models. The resultant field is characterized by strong latitudinal gradients (maintained by volume currents outside the source surface) and by abrupt reversals in direction at the current sheets. The model yields good agreement with the observed long-term variation of the radial IMF component in the ecliptic, and is used to predict the variation of |B _r | along the latitudinal trajectory of Ulysses during 1990–1994. As found in earlier studies, the magnitude ofB _r at any latitude is determined largely by the strength and relative orientation of the Sun's dipole moment.     

Observed Densities in the Universe

Baryons observed in Ly absorbers contribute to the density parameter 0 by bar 0.06 in close agreement with the value of 0.06 from primordial nucleosynthesis (H0=55 km s-1 Mpc-1, = 0 assumed throughout). A number of methods are known to measure 0 from density fluctuations; bound structures tend to yield lower values (m 0.2-0.4), field galaxies over large scales higher, but still undercritical values (m 0.6 ± 0.2). The best compromise value is 0 0.5, but the present methods are blind to diffusely distributed, exotic matter which still could make 0 = 1. A satisfactory solution of 0 (and ) will only come from a fundamental cosmological test (e.g. the Hubble diagram of [evolution-corrected] supernovae type Ia) in combination with the CMB fluctuation spectrum.