Search

We describe and model high-velocity (approximately 15,000 km S(exp -1)) Ly Alpha and H Alpha emission from the supernova remnant SNR 1987A seen in 1997 September and October with the Space Telescope Imaging Spectrograph. Part of this emission comes from a reverse shock located at approximately 75% of the radius of the inner boundary ofthe innercircumstellar ring and confined within + or - 30 degrees of the equatorial plane. Departure from axisymmetry in the Ly Alpha and H Alpha emission correlates with that seen in nonthermal radio emission and reveals an asymmetry in the circumstellar gas distribution. We also see diffuse high-velocity Ly-Alpha emission from supernova debris inside the reverse shock that may be due to excitation by nonthermal particles accelerated by the shock.

The high-velocity Ly(Alpha) emission from SN 1987A observed with the Space Telescope Imaging Spectrograph (STIS) evidently comes from a reverse shock formed where the outer envelope of SN 1987A strikes ionized gas inside the inner circumstellar ring. The observations can be explained by a simple kinematic model, in which the Ly(Alpha) emission comes from hydrogen atoms with radial velocity approximately 15,000 km s(exp -1) crossing a reverse shock in the shape of a slightly prolate ellipsoid with equatorial radius 4.8 x 10(exp 17) cm or approximately 80% of the distance to the inner surface of the inner ring. N v double Lambda 1239, 1243 emission, if present, has a net luminosity approximately less than 30% times that of the Ly(Alpha) emission. Future STIS observations should enable us to predict the time of impact with the inner ring and to determine unambiguously whether or not N v emission is present. These observations will offer a unique opportunity to probe the structure of SN 1987A's circumstellar environment and the hydrodynamics and kinetics of very fast shocks.

By the time the expanding envelope of a Type 2 supernova becomes transparent in the optical continuum, most of the gamma-ray luminosity produced by radioactive Fe/Co/Ni clumps propagates into the hydrogen/helium envelope and is deposited there, if at all. The resulting fast electrons excite He 1 and H 1, the two- photon continua of which are the dominant internal sources of ultraviolet radiation. The UV radiation is blocked by scattering in thousands of resonance lines of metals and converted by fluorescence into optical and infrared emission lines that escape freely. We describe results of Monte Carlo calculations that simulate non-LTE scattering and fluorescence in more than five million allowed lines of Ca, Sc, Ti, V, Cr, Mn, Fe, Co, and Ni. For a model approximating conditions in the envelope of SN 1987A, the calculated emergent spectrum resembles the observed one. For the first 2 yr after explosion, the ultraviolet radiation (lambda less than or approximately equals 3000) is largely blocked and converted into a quasi continuum of many thousands of weak optical and infrared emission lines and some prominent emission features, such as the Ca 2 lambdalambda8600 triplet. Later, as the envelope cools and expands, it becomes more transparent, and an increasing fraction of the luminosity emerges in the UV band.

The observations of CO and SiO in the infrared spectrum of SN 1987A clearly indicate that molecules can form in the debris of a supernova explosion. Since H2 is not easily observable we compute its abundance theoretically. For conditions typical of the inner (v less than 2500 km/s) envelope of SN 1987A, the fraction of H that is in molecular form rises to approx. 1% by t approx. 800 days. For t less than 500 days the formation is dominated by the gas-phase reactions H + H(+) yields H2(+) + hv; H2(+) + H yields H2 + H(+). Thereafter, the formation is dominated by the reactions H + e yields H(-) + hv; H(-) + H yields H2 + e. At early times the H(-) may absorb approx. 10%-30% of visible photons, contributing to the apparent paucity of H alpha emission. For t greater than 1000 days the abundance of H2 'freezes out' due to the slowing of all reactions. The opacity of the supernova envelope in the range 912 less than lambda less than or approx. equal to 1400 A (the upper limit depending on temperature) is dominated by resonance scattering in the Lyman and Werner bands of H2. The resulting fluorescence emission bands of H2 in the range 1150 less than lambda less than 1650 A may be observable in the UV spectra of supernovae at late times.

We model the evolution of the He I lines, especially He I 1.083 micrometers and He I 2.058 micrometers, from the envelope of SN 1987A during the first 3 yr after the explosion. The He I metastable 2(sup 3)S and 2(sup 1)S states are populated by fast electrons resulting from energy deposition by gamma rays and depopulated by electron collisions and two-photon decay. The 2.058 micrometer line is excited mainly by nonthermal electrons at all times. The 1.083 micrometer line is excited mainly by thermal electron impact for t less than or approximately = 450 days and by recombination afterward. We can fit the observed evolution of He I 2.058 micrometers with a model in which approximately 3 solar mass of nearly pure helium and approximately 11 solar mass of hydrogen mixed with primordial helium occupy the volume emitting most of the light (v less than or approximately = 2500 km/s). For t less than or approximately = 450 days, trace carbon in the inner part of the helium suppresses the He I 2.058 micrometer line by photoabsorbing He I lambda 584 and depopulating the 2(sup 1) P(sup 0) state. The observed 1.08 micrometer line is a blend of He I 1.083 micrometers and (S I) 1.082 micrometers. We estimate that He I 1.083 micrometers contributes 40% of the observed flux for t less than or approximately = 550 days and approximately 60% thereafter.

The envleope of SN 1987A will strike its circumstellar ring in 12 +/- 3 yr after explosion (A.D. 1999+/-3), the exact time depending weakly on the uncertain density of diffuse gas between the supernova and the ring. The impact will drive a radiative shock into the ring with velocity approximatley 200-400 km s(exp -1). The shocked ring will become a bright optical and ultraviolet emsiison-line source. A bright arc will suddenly appear at the near side of the ring and grow into an entire ring about 11 months later. The luminosities of the brightest lines, H-alpha H-alpha, N v lambda lambda 1238, 1242, and O VI lambda lambda 1032, 1038 will rise rapidly to approxminately 10(exp 36)-10(exp 37) ergs s(exp -1) and remain bright for several years after impact. The emission lines from the shocked ring will have FWHM approximatley 300-600 km s(exp -1) and complex profiles that will depend on position and will be sensitive to the details of the density distribution of gas in the ring. Strong EUV radiation from the shock will photoionize the unshocked ring, causing emission of narrow FWHM equivalent to 15 km s(exp -1) H-alpha, H-beta and (O III) lambda lambda 4959, 5007 lines with luminosities approximatley 10(exp 35) ergs s(exp -1).The EUV radiation will probably cause the nebulosity beyond the ring to become visable again. The EUV radiation may also illuminate the unshocked outer supernova envelope, causing visible emission of broad FWHM equivalent to 10(exp 4) km s(exp -1) H-alpha and H-beta lines.

During the first several weeks after the explosion of SN 1987A, Fe/Co/Ni clumps, containing approx. 1% of the mass of the supernova envelope, absorbed most of the energy release by Ni-56 and Co-56 decay. As a result, the clumps expanded relative to the substrate, forming a 'nickel bubble' of low-density Fe/Co/Ni. Later the clumps captured approx. 10% of the radioactive luminosity of gamma rays and positrons. Assuming that these elements are not mixed microscopically with other elements, we find that the clumps must occupy approx. greater than 30% of the volume of the emitting region (radial velocity approx. less than 2500 km/s). The result indicates that the emission at late times is dominated by an extra source of heating and ionization, most likely photoionization by two-photon continuum from metastable helium in the gas surrounding the clumps. The resulting 'frothy' structure, consisting of bubbles of low-density Fe surrounded by higher-density filaments of H, He, and other elements, will persist and may be seen in the spectra and structure of supernova remnants.

The angular autocorrelation function and spectrum of the soft X-ray background is studied below a discrete source detection limit, using two deep images from the Rosat X-ray satellite. The average spectral shape of pointlike sources, which account for 40 to 60 percent of the background intensity, is determined by using the autocorrelation function. The background spectrum, in the 0.5-0.9 keV band (M band), is decomposed into a pointlike source component characterized by a power law and a diffuse component represented by a two-temperature plasma. These pointlike sources cannot contribute more than 60 percent of the X-ray background intensity in the M band without exceeding the total observed flux in the R7 band. Spectral analysis has shown that the local soft diffuse component, although dominating the background intensity at energies not greater than 0.3 keV, contributes only a small fraction of the M band background intensity. The diffuse component may represent an important constituent of the interstellar or intergalactic medium.

We study the evolution of the forbidden Ca II 7300 A and Ca II 8600 A emission lines in the spectrum SN 1987A. We can explain the observations with a model in which (1) the lines come from about 1.7 x 10 exp -4 solar mass of calcium distributed with filling factor about 0.1 in a sphere expanding with velocity 2500 km/s; (2) the electron number density is given by 4 x 10 exp 9/cu cm/(t/100d)-cubed; (3) the gas temperature decreases from T roughly 6400 K at 200 days to T roughly 3900 K at 500 days; and (4) ultraviolet pumping in the Ca II 3934 A, 3969 A (K, H) lines is required to explain the line ratios for t greater than about 350 days. The mass of calcium in our model is consistent with the LMC abundance of calcium in about 5 solar mass of hydrogen in the inner envelope of SN 1987A. If more calcium is produced by supernova nucleosynthesis, it must be in dense clumps that are too cool to emit the forbidden Ca II 7300 A and Ca II 8600 A lines.

Updates on research on SN 1987A are provided, with emphasis on the aftermath of the explosion: the interpretation of the light curve and spectrum, the physics of the envelope, and the interaction of the supernova with its circumstellar envelope. UVOIR light curves, gamma rays and hard X-rays, and dust formation are examined, and explosion dynamics is explored. Attention is given to optical and IR spectrum development, with emphasis on the photospheric spectrum, nebular spectrum, ionization, emission-line profiles, line strengths, and iron, cobalt, and nickel. The circumstellar environment is discussed, with special reference to progenitor mass loss, circumstellar ring and nebulosity, radio emission, and soft X-rays.

The radio and X-ray images of the supernova remnant CTB 109 have the morphology of a semicircular shell. It is shown that such a structure is a natural result of a supernova explosion that occurred 3,000 yr ago near an interface between the diffuse interstellar medium and a dense molecular cloud. The calculated X-ray, infrared, and radio fluxes agree fairly well with the observed values.

A general theory is presented for hydrogen recombination line formation in an expanding medium in which some of the lines are optically thick. This theory is used to calculate the time evolution of the hydrogen lines of SN 1987A at t equal to or greater than 150 days, assuming that the supernova envelope is a homologously expanding uniform sphere. The theoretical luminosities and ratios of the recombination lines agree remarkably well with the observations. For the first 2 yr, the supernova envelope is optically thick to Balmer continuum. For t equal to or less than 400 days, hydrogen is ionized primarily from the n = 2 level by Balmer continuum photons, which are provided partly by the two-photon decay of the 2s state and partly by emission lines of heavy elements.

An hourglass-shaped circumstellar shell, with radius about 0.65 lt-yr at the equator and about 2.4 lt-yr at the poles, surrounds SN 1987A. The shell is a bubble blown by the wind of the blue giant progenitor of SN 1987A into an equatorially concentrated wind of a prior red giant stage. The bubble was ionized by the initial UV flash of the supernova, and the emission-line region is concentrated in a ring at the waist of the hourglass. In about 2002 AD the supernova envelope will strike the shell. The resulting young supernova remnant will become a luminous source of soft X-rays and broad optical and ultraviolet emission lines.

An equation is derived for calculating the energy distribution of fast electrons in a partially ionized gas and a method is provided to solve for the electron degradation spectrum and the energy deposition in different forms (ionization, excitation, or heating). As an example, the energy degradation of fast electrons in a gas of pure hydrogen is calculated, considering excitations to the lowest 10 atomic levels. The Bethe approximation and the continuous slowing-down approximation are discussed and it is concluded that these approximations are accurate to the order of 20 percent for electrons with initial energy of greater than about keV. The method and results can be used to determine heating, excitations, and ionizations by high-energy photoelectrons or cosmic-ray particles in various astrophysical circumstances, such as the interstellar medium, supernova envelopes, and QSO emission-line clouds.

A dense circumstellar shell, with radius of about 0.5 lt-yr, surrounds SN 1987 A. In 16 yr or less after outburst, the expanding debris of SN 1987 A will strike this shell. When it does, the hot gas and relativistic electrons resulting from the forward and reverse shocks will radiate X-rays, infrared radiation, and nonthermal radio waves. The remnant of SN 1987 A will then become one of the brightest X-ray and radio sources in the LMC.

An analytic method is developed for calculating the emergent spectrum of gamma-rays and X-rays scattered in a homogeneous medium with low-temperature electrons. The Klein-Nishina corrections of the scattering cross section and absorption processes are taken in account. The wavelength relaxation and the spatial diffusion problems are solved separately, and the emergent spectrum is calculated by convolving the evolution function of the spectrum in an infinite medium with the photon luminosity resulting from the spatial diffusion in a finite sphere. The analytic results are compared with that of Monte Carlo calculations and it is concluded that the analytic result is quite accurate.

Recommendations are given regarding National Science Foundation (NSF) astronomy programs and the NASA Space Astrophysics program. The role of ground based astronomy is reviewed. The role of National Optical Astronomy Observatories (NOAO) in ground-based night-time astronomical research is discussed. An enhanced Explored Program, costs and management of small and moderate space programs, the role of astrophysics within NASA's space exploration initiative, suborbital and airborne astronomical research, the problems of the Hubble Space Telescope, and astronomy education are discussed. Also covered are policy issues related to the role of science advisory committees, international cooperation and competition, archiving and distribution of astronomical data, and multi-wavelength observations of variable sources.

A stellar wind from a massive OB or Wolf-Rayet star in a binary system will strike the surface or stellar wind of its companion, forming shocked gas that can radiate X-rays. The X-ray spectrum from the shocked winds will vary in a predictable way with orbital phase, owing to photoelectric absorption by the stellar winds. Detailed models are calculated for the hydrodynamics and X-ray emission from two such systems. In one of these systems (HD 165052), the winds are nearly identical in strength. In the other (V444 Cygni), the wind of the Wolf-Rayet star overwhelms and crushes that of its companion. The calculated X-ray luminosities agree fairly well with the observed values for HD 165052 and for V444 Cygni. These results can be scaled to other such systems.