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.     

Energetic particles at high latitudes

We review work on diffusion coefficients of energetic particles with an attempt to extract implications on their behaviour at high latitudes. In the ecliptic plane results from solar energetic particle propagation between the Sun and about 5 AU can be described by an effective radial mean free path _r which is approximately constant as a function of distancer. When particle propagation in three dimensions in the heliosphere is considered it is not sufficient to consider _r only. Jovian electrons can be used as probes to determine the parameters of three-dimensional diffusion. In the polar regions diffusion is dominated by its parallel component. Some predictions how should vary with latitude are discussed. For different choices of this variation we present expectations for intensity-time profiles of solar particle events during the Ulysses polar passages.     

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.     

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.     

SOHO: The Solar and Heliospheric Observatory

The Solar and Heliospheric Observatory (SOHO), together with the Cluster mission, constitutes ESA's Solar Terrestrial Science Programme (STSP), the first Cornerstone of the Agency's long-term programme Space Science — Horizon 2000. STSP, which is being developed in a strong collaborative effort with NASA, will allow comprehensive studies to be made of the both the Sun's interior and its outer atmosphere, the acceleration and propagation of the solar wind and its interaction with the Earth. This paper gives a brief overview of one part of STSP, the SOHO mission.     

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.     

Electron coronal density irregularity\overline {n^2 } /(\bar n)^2 from measurements of K-coronal and lyman lines brightnesses

We propose a technique to derive the coronal density irregularity factor , wheren is the electron density. The absolute photometric comparison between the intensity of UV lines and the white-light K-coronal polarized brightness (pB) provides an unique constraint on the inhomogeneity of the corona. 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 collisonal component of the Ly-. This component yields an estimate of . The quantity is then obtained from white-light K-coronal measurements. The use of lines of the same atomic species minimizes the effects due to outflow velocities (i.e., Doppler dimming), and reduces the errors introduced by the uncertainties in the ionization balance, the atomic parameters, and the solar abundances. The UVCS/SOHO unique capability of performing cotemporal and cospatial measurements of the Ly- and Ly- lines, and ofpB makes this instrument ideal for implementing this technique.     

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.     

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.     

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.     

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.     

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.     

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     

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

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.     

Primordial Helium and ΔY/ΔZ from H II Regions and from Fine Structure in the Main Sequence Based on Hipparcos Parallaxes

The primordial helium abundance YP is important for cosmology and the ratio Y/Z of the changes relative to primordial abundances constrains models of stellar evolution. While the most accurate estimates of YP come from emission lines in extragalactic H II regions, they involve an extrapolation to zero metallicity which itself is closely tied up with the slope Y/Z. Recently certain systematic effects have come to light in this exercise which make it useful to have an independent estimate of Y/Z from fine structure in the main sequence of nearby stars. We derive such an estimate from Hipparcos data for stars with Z Z and find values between 2 and 3, which are consistent with stellar models, but still have a large uncertainty.     

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.     

Reverse Algols

An Algol is a binary system having a semidetached configuration where the less massive component is in contact with the critical equipotential surface. A reverse Algol is a binary system having a semidetached configuration where the more massive component is in contact with the critical equipotential surface. In 1985, Leung suggested 5 reverse Algol systems at the Beijing Colloquium. Two more such systems have been discovered recently. The spectral types of these systems range from early B to mid G. There is also a wide spread in mass ratio among these systems. There appear to be two types of reverse Algols, hot and cool systems. The hot systems have their more massive components as the hotter stars and the cool systems their more massive components as the cooler stars. The mass-radius relation of the reverse Algols is very similar to that of the contact and near-contact systems. It is believed that reverse Algols represent the pre-mass-reversal semidetached phase of close binary evolution. Since selection effects apply to both the regular Algols and reversed Algols in a similar manner, the evolutionary time scale between them would be simply the ratio of the number of confirmed systems of these two types of Algols.     

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.