High ionization star forming (SF) galaxies are easily identified with strong emission line techniques such as the BPT diagram, and form an obvious ionization sequence on such diagrams. We use a locally optimally emitting cloud model to fit emission line ratios that constrain the excitation mechanism, spectral energy distribution, abundances and physical conditions along the star-formation ionization sequence. Our analysis takes advantage of the identification of a sample of pure star-forming galaxies, to define the ionization sequence, via mean field independent component analysis. Previous work has suggested that the major parameter controlling the ionization level in SF galaxies is the metallicity. Here we show that the observed SF- sequence could alternatively be interpreted primarily as a sequence in the distribution of the ionizing flux incident on gas spread throughout a galaxy. Metallicity variations remain necessary to model the SF-sequence, however, our best models indicate that galaxies with the highest and lowest observed ionization levels (outside the range -0.37 < log [O III]/H\b{eta} < -0.09) require the variation of an additional physical parameter other than metallicity, which we determine to be the distribution of ionizing flux in the galaxy.

To identify extragalactic magnetar flares, we have searched for their periodic tails by generating Lomb periodograms of the emission following short bursts detected by the Burst and Transient Source Experiment (BATSE). Out of 358 short bursts examined, one has a significant tail periodicity (T = 13.8 s, P = 4 x 10^-5). The most probable host galaxy for this burst is "The Fireworks Galaxy" NGC 6946 (d = 5.9 Mpc). At this distance, the energy of the spike, (2.7 +/- 0.3) x 10^44 ergs, is akin to those of the galactic magnetar giant flares, as are its duration (~ 0.4 s) and temperature (250 +/- 60 keV). For the tail emission, however, our estimated temperature of 60 +/- 5 keV is harder and the energy release of (4.3 +/- 0.8) x 10^45 ergs is larger than those of the galactic magnetar flares. Regardless of the host, such a large ratio of tail-to-spike energy would imply that magnetar flare tails might be detectable out to further distances than previously thought.

We suggest a simple fitting formula to represent Comptonized X- and gamma-ray spectra from a hot ($kT_e \gtrsim 10$ keV), Thomson thick ($\tau_T \gtrsim 5$) hybrid thermal/nonthermal plasma in spherical geometry with homogeneous soft-photon injection throughout the Comptonizing region. We have used this formula to fit a large data base of Monte-Carlo generated photon spectra, and provide correlations between the physical parameters of the plasma and the fit parameters of our analytic fit function. These allow us to construct Thomson thick Comptonization spectra without performing computer intensive Monte Carlo simulations of the high-Thomson-depth hybrid-plasma Comptonization problem. Our formulae can easily be used in data analysis packages such as XSPEC, thus rendering rapid $\chi^2$ fitting of such spectra to real data feasible.

The recent simultaneous detection of optical, X-ray and gamma-ray photons from GRB990123 during the burst provides the first broadband multi-wavelength characterization of the burst spectrum and evolution. Here we show that a direct correlation exists between the time-varying gamma-ray spectral shape and the prompt optical emission. This combined with the unique signatures of the time-resolved spectra of GRB990123 convincingly supports earlier predictions of the saturated Comptonization model. Contrary to other suggestions, we find that the entire continuum from optical to gamma-rays can be generated from a single source of leptons (electrons and pairs). The optical flux only appears to lag the gamma-ray flux due to the high initial Thomson depth of the plasma. Once the plasma has completely thinned out, the late time afterglow behavior of our model is the same as in standard models based on the Blandford-McKee (1976) solution.

GX 339-4 is an unusual black hole candidate in that it is a persistent source, being detected by X-ray telescopes most of the time, but it also has nova-like flaring states. In 1996 we performed radio, optical, X-ray, and gamma-ray observations of GX 339-4 when it was in a hard state (= soft X-ray low state). Here we present a brief summary of some of the results.

As part of our multiwavelength campaign of GX 339-4 observations in 1996 we present our radio, X-ray, and gamma-ray observations made in July, when the source was in a hard state (= soft X-ray low state). The radio observations were made at the time when there was a possible radio jet. We show that the radio spectrum was flat and significantly variable, and that the radio spectral shape and amplitude at this time were not anomalous for this source. Daily light curves from our pointed observation July 9-23 using OSSE, from BATSE, and from the ASM on RXTE also show that there was no significant change in the X- and gamma-ray flux or hardness during the time the possible radio jet-like feature was seen. The higher energy portion of our pointed RXTE observation made July 26 can be equally well fit using simple power law times exponential (PLE) and Sunyaev-Titarchuk (ST) functions. An additional soft component is required, as well as a broad emission feature centered on 6.4 keV. This may be an iron line that is broadened by orbital Doppler motions and/or scattering off a hot medium. Its equivalent width is 600 eV. Our simplistic continuum fitting does not require an extra reflection component. Both a PLE and a ST model also fit our OSSE spectrum on its own. Although the observations are not quite simultaneous, combining the RXTE and CGRO spectra we find that the PLE model easily fits the joint spectrum. However, the ST model drops off too rapidly with increasing energies to give an acceptable joint fit.

We have analyzed the evolution of the spectral hardness parameter Epk as a function of fluence in gamma-ray bursts. We fit 41 pulses within 26 bursts with the trend reported by Liang & Kargatis (1996) which found that Epk decays exponentially with respect to photon fluence. We also fit these pulses with a slight modification of this trend, where Epk decays linearly with energy fluence. In both cases, we found the set of 41 pulses to be consistent with the trend. For the latter trend, which we believe to be more physical, the distribution of the decay constant is roughly log-normal, with a mean of 1.75 +/- 0.07 and a FWHM of 1.0 +/- 0.1. Regarding an earlier reported invariance in the decay constant among different pulses in a single burst, we found probabilities of 0.49 to 0.84 (depending on the test used) that such invariance would occur by coincidence, most likely due to the narrow distribution of decay constant values among pulses.

If gamma-ray bursts are associated with dense star-forming regions in young galaxies, photoelectric absorption by the dense circumburster material (CBM) will occur. As the burst evolves, the surrounding material is photoionized, leading to fluorescence line emission and reduced photoelectric absorption opacity. We have analyzed this process in detail, accounting for the time-dependent photoelectric absorption, photoionization and fluorescence line emission from the CBM. We find that even if GRBs are hosted in dense star-forming regions, photoionization of the GRB environment leads to a constant, but very weak level of delayed fluorescence line emission on timescales of weeks to years after the burst. A temporally evolving iron~K edge absorption feature can serve as diagnostic tool to reveal the density structure of the CBM and may provide an opportunity for redshift measurements. We also investigated whether photoelectric absorption could be responsible for the spectral evolution of the low-energy slopes of some bright BATSE $\gamma$-ray bursts displaying extremely hard spectra below the peak energy, inconsistent with the optically-thin synchrotron shock model. We find that a very strong metal enrichment (~ 100 times solar-system abundances) in the gamma-ray burst environment and a rather peculiar spatial distribution of the CBM would be necessary in order to account for the observed hard spectra below a few 100 keV and their temporal evolution.

Recent studies have found that the spectral evolution of pulses within gamma-ray bursts (GRBs) is consistent with simple radiative cooling. Perhaps more interesting was a report that some bursts may have a single cooling rate for the multiple pulses that occur within it. We determine the probability that the observed "cooling rate invariance" is purely coincidental by sampling values from the observed distribution of cooling rates. We find a 0.1-26% probability that we would randomly observe a similar degree of invariance based on a variety of pulse selection methods and pulse comparison statistics. This probability is sufficiently high to warrant skepticism of any intrinsic invariance in the cooling rate.

The time-resolved spectra of gamma-ray bursts (GRBs) remain in conflict with many proposed models for these events. After proving that most of the bursts in our sample show evidence for spectral "shape-shifting", we discuss what restrictions that BATSE time-resolved burst spectra place on current models. We find that the synchrotron shock model does not allow for the steep low-energy spectral slope observed in many bursts, including GRB 970111. We also determine that saturated Comptonization with only Thomson thinning fails to explain the observed rise and fall of the low-energy spectral slope seen in GRB 970111 and other bursts. This implies that saturated Comptonization models must include some mechanism which can cause the Thomson depth to increase intially in pulses.

Using Monte Carlo methods to simulate the inverse Compton scattering of soft photons, we model the spectrum of the Galactic black hole candidate Cyg X-1, which shows evidence of a nonthermal tail extending beyond a few hundred keV. We assume an ad hoc sphere of leptons, whose energy distribution consists of a Maxwellian plus a high energy power-law tail, and inject 0.5 keV blackbody photons. The spectral data is used to constrain the nonthermal plasma fraction and the power-law index assuming a reasonable Maxwellian temperature and Thomson depth. A small but non-negligible fraction of nonthermal leptons is needed to explain the power-law tail.

We report evidence that the asymptotic low-energy power law slope alpha (below the spectral break) of BATSE gamma-ray burst photon spectra evolves with time rather than remaining constant. We find a high degree of positive correlation exists between the time-resolved spectral break energy E_pk and alpha. In samples of 18 "hard-to-soft" and 12 "tracking" pulses, evolution of alpha was found to correlate with that of the spectral break energy E_pk at the 99.7% and 98% confidence levels respectively. We also find that in the flux rise phase of "hard-to-soft" pulses, the mean value of alpha is often positive and in some bursts the maximum value of alpha is consistent with a value > +1. BATSE burst 3B 910927, for example, has a alpha_max equal to 1.6 +/- 0.3. These findings challenge GRB spectral models in which alpha must be negative of remain constant.