The future of space plasma simulation

Space plasma simulation is a subject which is in its infancy, but which is already having an important impact on space science. Its growth is being spurred onward by the continuing increase in capacity (speed and memory) of computers and by advances in the sophistication of numerical models. These advances are making it possible to simulate more realistic situations using more complex models. Already significant three-dimensional MHD calculations of the magnetosphere and its interaction with the solar wind have been carried out. In addition multi-dimensional particle simulations are illuminating many of the microscopic physics processes which go on (instabilities, saturation levels and wave nonlinearity, shock structure, etc.). Notwithstanding these advances, the surface has only been scratched; many challenges and opportunities are provided by simulation both for the space physicist and the model builders (also for computer designers). In MHD models more physics need to be included (Hall effect, gyroviscosity, accurate models of boundaries, how do we put microscopic physics effects into macroscopic codes, etc.). For model builders correct treatments of systems containing a large range of important space and time scales, magnetic field strengths, Alfven wave velocity, etc. present real challenges. What are the best ways to diagnose these complex models and obtain meaningful information? What quantities should be looked at? How should they be displayed? A discussion of the promises, the prospects, and the challenges of the above topics will be given with examples taken from recent work.     

EXOSAT medium energy observations of Cyg X-2

On July 5.–6. 1983, during the EXOSAT performance verification (PV) and calibration phase, a raster scan of Cygnus X-2 was performed. In contrast to the previously observed smooth intensity variations on timescales of hours, the source revealed a behaviour unknown until now: active periods with high energy flares recurring on time scales of 300–500 s were interrupted by quiet periods of several hours. At all intensity levels the source spectra clearly require a two component continuum (blackbody + thermal bremsstrahlung). In addition, a weak iron emission line with equivalent widths between 39 an 70 eV was detected. The source has a much harder spectrum during the flares than during quiet periods, indicating drastic temperature changes within the emission region, while the absolute iron line flux does not vary. From the spectral characteristics it becomes clear that self-comptonization of the thermal bremsstrahlung spectrum plays an important role. The time variability and spectral behaviour in this peculiar state allow Cyg X-2 to be classified as a Low Mass X-ray Binary System (LMXB) very similar to the prototype of this class, Sco X-1.     

The long term variability of 4U 1258-61 (GX304-1)

We present optical spectroscopy and photometry of 4U 1258-61 obtained over the course of six years. Evidence for the evolution of a circumstellar envelope throughout the course of the observations is presented.     

An X-ray study of M51 (NGC 5194) and its companion (NGC 5195)

Observations of NGC 5194/95 with the Einstein HRI show a very strong nuclear X-ray source, surrounded by a diffuse flux, three point sources and the companion. The diffuse flux, which correlates well with the radio continuum, is likely to originate from the disk population with age 2·10⁹ yrs. The large luminosity from the nuclear source, together with optical and radio observations, shows that it belongs to the low luminosity active nuclei, thus extending this class to luminosities less than 10⁴⁰ erg/s.     

Spectral computation of Tandem mirror equilibria with Finite Larmor Radius effects

A direct iterative method of solving for Tandem equilibria by moving magnetic field lines in a manner to satisfy the linearized equilibrium equations converges much more rapidly than standard relaxation techniques, typically in under a fifty iterations. At the highest 's the number of iterations increase, but is still far less than other methods. In quadrupole tandem mirror equilibrium, octupole and higher distortions of the flux surfaces are important which forces us to abandon finite differences in the angle-like flux variable and resort to a spectral decomposition to solve the equilibrium equations. We display equilibria at the high expected for MFTF-B and show how Finite Larmor Radius (FLR) effects strongly suppress these azimuthal distortions.     

The Broad Band X-Ray Telescope

The Broad-Band X-Ray Telescope (BBXRT) has been designed to perform high sensitivity, moderate resolution spectrophototnetry of X-ray sources in the 0.3–12 keV band from the Shuttle. It consists of a coaligned pair of high throughput, conical X-ray imaging mirrors, with a cryogenically-cooled, multiple element, Si(Li) spectrometer at the focus of each. On axis, BBXRT will have an effective area of 580 cm² at 2 keV and 250 cm² at 7 keV, and a spectral resolution of 110 eV at 2 keV and 150 eV at 7 keV. A 10⁴ s observation with BBXRT will allow a determination of the continuum spectral shape for sources near the Einstein deep survey limit.     

Some crucial X-ray observations — One or two SNR(s) in Crux?

A new EXOSAT (LE/CMA) observation of the region in Crux (R.A. 11h 45m, Dec. -62°) where Markert et al. (1981) reported the existence of two x-ray SNR's is presented. After cleaning the CMA field from the point source component, due to the UV emission of the numerous stars in the field, the smoothed x-ray contours are compared to the 408 MHz radio map of Caswell et al. (1983). The existence of two, well-separate x-ray emission regions is confirmed by EXOSAT, and the current x-ray/radio picture is not sufficent to distinguish clearly between the assumption of one or two (possibly interacting) SNR's in the region.     

Long term observations of Cyg X-3 with Vela 5B

We present long-term (1969–1979) observations of Cygnus X-3, obtained by the Vela 5B satellite. The 3–12 keV light curve has 10 day time resolution. Cyg X-3 is a peculiar high-luminosity X-ray source, radiating from the radio region to hard gamma rays of more than 10¹⁶ eV. It has a 4.8 hour period, probably orbital, which is not resolved by our present analysis. Long term periodicities of 17, 20, and 33–34 days have been reported by several authors, and explained as the effects of apsidal motion, precession, or an eccentric orbit. We do not observe the 17 and 33–34 day variations, and set upper limits significantly lower than the reported amplitude of the 33–34 day variation. There is weak evidence for a 20 day flux variation. The light curve shows high and low states which alternate with a characteristic timescale of 1 year. There is no counterpart, at this time resolution, of the giant radio outburst of 1972 September.     

X-ray scattering by intergalactic dust

The imaging capabilities of the Exosat and Einstein satellites at soft X-ray wavelengths have begun to show that suitable Galactic X-ray sources have extended ( 10 arcmin) haloes due to scattering of soft X-rays by interstellar dust. A simple argument suggests that similar haloes, due to scattering by intergalactic dust, should exist around distant (z 1) quasars and detailed analysis confirms this conclusion. A search for such haloes around suitable X-ray quasars could provide valuable, model-independent, constraints on the amount and origin of intergalactic dust.     

Soft X-ray emission from radio Loops I and III compared

We examine X-ray maps of the whole sky covering energies from 0.1 to 6 keV. We model Loop I X-ray emission as being caused by an old supernova that exploded into an already warm interstellar medium. After comparison with Loop III we deduce that there may be a general temperature gradient in the gas as we move away from the plane in the northern galactic hemisphere.