We present the results of a high angular resolution study of the BL Lac object Markarian 501 in the radio band. We consider data taken at 14 different epochs, ranging between 1.6 GHz and 22 GHz in frequency, and including new Space VLBI observations obtained on 2001 March 5 and 6 at 1.6 and 5 GHz. We study the kinematics of the parsec-scale jet and estimate its bulk velocity and orientation with respect to the line of sight. Limb brightened structure in the jet is clearly visible in our data and we discuss its possible origin in terms of velocity gradients in the jet. Quasi-simultaneous multi-wavelength observations allow us to map the spectral index distribution and to compare it to the jet morphology. Finally, we estimate the physical parameters of the parsec-scale jet.

We report EGRET upper limits on the high-energy gamma-ray emission from clusters of galaxies. EGRET observations between 1991 and 2000 were analyzed at positions of 58 individual clusters from a flux-limited sample of nearby X-ray bright galaxy clusters. Subsequently, a coadded image from individual galaxy clusters has been analyzed using an adequately adapted diffuse gamma-ray foreground model. The resulting 2 sigma upper limit for the average cluster is \~ 6 x 10^{-9} cm^{-2} s^{-1} for E > 100 MeV. Implications of the non--detection of prominent individual clusters and of the general inability to detect the X-ray brightest galaxy clusters as a class of gamma-ray emitters are discussed. We compare our results with model predictions on the high-energy gamma-ray emission from galaxy clusters as well as with recent claims of an association between unidentified or unresolved gamma-ray sources and Abell clusters of galaxies and find these contradictory.

We present the results of deep X-ray and gamma-ray observations of the Geminga pulsar obtained in the final years of the ASCA and CGRO missions, and an upper limit from RXTE. A phase-connected ephemeris from the gamma-rays is derived that spans the years 1973-2000, after allowing for a minor glitch in frequency of Delta f/f = 6.2 x 10^-10 in late 1996. ASCA observations of the hard X-ray pulse profile in 1994 and 1999 confirm this glitch. An improved characterization of the hard X-ray pulse profile and spectrum from the long ASCA observation of 1999 confirms that there is a non-thermal X-ray component that is distinct from the gamma-ray spectrum as measured by EGRET. It can be parameterized as a power-law of photon index Gamma = 1.72 +/- 0.10 with a flux of 2.62 x 10^-13 ergs/cm^2/s in the 0.7-5 keV band and pulsed fraction 0.54 +/- 0.05, similar to, but more precise than values measured previously. An extrapolation of this spectrum into the energy band observed by the RXTE PCA is consistent with the non-detection of pulsed emission from Geminga with that instrument. These results are interpreted in the context of outer-gap models, and motivations for future X-ray observations of Geminga are given.

We perform X-ray studies of the shell-type supernova remnant (SNR) gamma-Cygni associated with the brightest EGRET unidentified source 3EG J2020+4017. In addition to the thermal emissions with characteristic temperature of kT = 0.5-0.9 keV, we found an extremely hard X-ray component from several clumps localized in the northern part of the remnant. This component is described by a power-law with a photon index of 0.8-1.5. Both the absolute flux and the spectral shape of the nonthermal X-rays cannot be explained by the synchrotron or inverse-Compton mechanisms. We argue that the unusually hard X-ray spectrum can be naturally interpreted in terms of nonthermal bremsstrahlung from Coulomb-loss-flattened electron distribution in dense environs with the gas density about 10 to 100 cm^-3 . For given spectrum of the electron population, the ratio of the bremsstrahlung X- and gamma-ray fluxes depends on the position of the ``Coulomb break'' in the electron spectrum. The bulk of gamma-rays detected by EGRET would come from the radio-bright and X-ray dim cloud at southeast, where very dense gas and strong magnetic field would illuminate the cloud in the radio and gamma-ray bands, but suppress the bremsstrahlung X-ray emission due to the shift of the ``Coulomb break'' in the electron spectrum towards higher energies.

New optical and radio data on the BL Lacertae object AO 0235+16 have been collected in the last four years by a wide international collaboration, which confirm the intense activity of this source. The optical data also include the results of the Whole Earth Blazar Telescope (WEBT) first-light campaign organized in November 1997. The optical spectrum is observed to basically steepen when the source gets fainter. We have investigated the existence of typical variability time scales and of possible correlations between the optical and radio emissions by means of visual inspection, Discrete Correlation Function analysis, and Discrete Fourier Transform technique. The major radio outbursts are found to repeat quasi-regularly with a periodicity of about 5.7 years; this period is also in agreement with the occurrence of some of the major optical outbursts, but not all of them.

We examine the coincidence of times of high $\gamma$-ray flux and ejections of superluminal components from the core in EGRET blazars based on a VLBA monitoring program at 22 and 43 GHz from November 1993 to July 1997. In 23 cases of $\gamma$-ray flares for which sufficient VLBA data exist, 10 of the flares (in 8 objects) fall within 1$\sigma$ uncertainties of the extrapolated epoch of zero separation from the core of a superluminal radio component. In each of two sources (0528+134 and 1730-130) two successive $\gamma$-ray flares were followed by the appearance of new superluminal components. We carried out statistical simulations which show that if the number of coincidences $\ge$ 10 the radio and $\gamma$-ray events are associated with each other at >99.999% confidence. Our analysis of the observed behavior, including variability of the polarized radio flux, of the sources before, during, and after the $\gamma$-ray flares suggests that the $\gamma$-ray events occur in the superluminal radio knots. This implies that the $\gamma$-ray flares are caused by inverse Compton scattering by relativistic electrons in the parsec-scale regions of the jet rather than closer to the central engine.

We present the results of a program to monitor the structure of the radio emission in 42 $\gamma$-ray bright blazars (31 quasars and 11 BL Lac objects) with the VLBA at 43, 22, and occasionally 15 and 8.4 GHz, over the period from November 1993 to July 1997. We determine proper motions in 33 sources and find that the apparent superluminal motions in $\gamma$-ray sources are much faster than for the general population of bright compact radio sources. This follows the strong dependence of the $\gamma$-ray flux on the level of relativistic beaming for both external-radiation Compton and synchrotron self-Compton emission. There is a positive correlation (correlation coefficient $r$=0.45) between the flux density of the VLBI core and the $\gamma$-ray flux and a moderate correlation (partial correlation coefficient $r$=0.31) between $\gamma$-ray apparent luminosity and superluminal velocities of jet components, as expected if the $\gamma$-ray emission originates in a very compact region of the relativistic jet and is highly beamed. In 43% of the sources the jet bends by more than 20$^\circ$ on parsec scales, which is consistent with amplification by projection effects of modest actual changes in position angle. In 27 of the sources in the sample there is at least one non-core component that appears to be stationary during our observations. Different characteristics of stationary features close to and farther from the core lead us to suggest two different classes of stationary components: those within about 2 milliarcseconds (mas) of the core, probably associated with standing hydrodynamical compressions, and those farther down the jet, which tend to be associated with bends in the jet.

We conducted a multi-frequency campaign for the TeV blazar Mrk~421 in 1998 April. The campaign started from a pronounced high amplitude flare recorded by SAX and Whipple; ASCA observation started three days later. In the X-ray data, we detected multiple flares, occuring on time scales of about one day. ASCA data clearly reveal spectral variability. The comparison of the data from ASCA, EUVE and RXTE indicates that the variability amplitudes in the low energy synchrotron component are larger at higher photon energies. In TeV Gamma-rays, large intra-day variations -- which were correlated with the X-ray flux -- were observed when results from three Cherenkov telescopes are combined. The RMS variability of TeV Gamma--rays was similar to that observed in hard X-rays, above 10 keV. The X-ray light curve reveals flares which are almost symmetric for most of cases, implying the dominant time scale is the light crossing time through the emitting region. The structure function analysis based on the continuous X-ray light curve of seven days indicates that the characteristic time scale is ~0.5 day. The analysis of ASCA light curves in various energy bands appears to show both soft (positive) and hard (negative) lags. These may not be real, as systematic effects could also produce these lags, which are all much smaller than an orbit. If the lags of both signs are real, these imply that the particle acceleration and X-ray cooling time scales are similar.

We report the variable X-ray emission from BL Lacertae detected in the ASCA ToO observation conducted during the EGRET and RXTE pointings, coincident with the 1997 July outburst. The source showed a historically high state of X-ray, optical, and gamma-ray emission, with its 2-10 keV flux peaking at ~3.3 E-11 erg/cm2/s. This is more than 3 times higher than the value measured by ASCA in 1995. We detected two rapid flares which occured only in the soft X-ray band, while the hard X-ray flux also increased, but decayed with a much longer time scale. Together with the requirement of a very steep and varying power law dominating the soft X-ray band in addition to the hard power law, we suggest that both the high energy end of the synchrotron spectrum and the hard inverse Compton spectrum were visible in this source during the outburst. We discuss the possible origins of the observed variability time scales, and interpret the short time scales of the soft X-ray variability as reflecting the size of the emission region.

The optical light curve of the energetic gamma-ray burst GRB 991216 is consistent with jet-like behavior in which a power-law decay steepens from t**(-1.22 +/- 0.04) at early times to t**(-1.53 +/- 0.05) in a gradual transition at around 2 d. The derivation of the late-time decay slope takes into account the constant contribution of a host or intervening galaxy which was measured 110 d after the event at R = 24.56 +/- 0.14, although the light curve deviates from a single power law whether or not a constant term is included. The early-time spectral energy distribution of the afterglow can be described as F_nu ~ nu**(-0.74 +/- 0.05) or flatter between optical and X-ray, which, together with the slow initial decay, is characteristic of standard adiabatic evolution in a uniformly dense medium. Assuming that a reported absorption-line redshift of 1.02 is correct, the apparent isotropic energy of 6.7 x 10**53 erg is reduced by a factor of ~ 200 in the jet model, and the initial half-opening angle is ~ 6 deg. GRB 991216 is the third good example of a jet-like afterglow (following GRB 990123 and GRB 990510), supporting a trend in which the apparently most energetic gamma-ray events have the narrowest collimation and a uniform ISM environment. This, plus the absence of evidence for supernovae associated with jet-like afterglows, suggests that these events may originate from a progenitor in which angular momentum plays an important role but a massive stellar envelope or wind does not, e.g., the coalescence of a compact binary.

Because of the scheduled GLAST mission by NASA, there is strong scientific justification for preparation for very extensive blazar monitoring in the optical bands to exploit the opportunity to learn about blazars through the correlation of variability of the gamma-ray flux with flux at lower frequencies. Current optical facilities do not provide the required capability.Developments in technology have enabled astronomers to readily deploy automatic telescopes. The effort to create an Automatic Telescope Network (ATN) for blazar monitoring in the GLAST era is described. Other scientific applications of the networks of automatic telescopes are discussed. The potential of the ATN for science education is also discussed.

We present results of Compton Gamma-Ray Observatory EGRET observations of the unidentified high-energy gamma-ray sources 2EG J1049-5847 (GEV J1047-5840, 3EG J1048-5840) and 2EG J1103-6106 (3EG J1102-6103). These sources are spatially coincident with the young, energetic radio pulsars PSRs B1046-58 and J1105-6107, respectively. We find evidence for an association between PSR B1046-58 and 2EG J1049-5847. The gamma-ray pulse profile, obtained by folding time-tagged photons having energies above 400 MeV using contemporaneous radio ephemerides, has probability of arising by chance of 1.2E-4 according to the binning-independent H-test. A spatial analysis of the on-pulse photons reveals a point source of equivalent significance 10.2 sigma. Off-pulse, the significance drops to 5.8 sigma. Archival ASCA data show that the only hard X-ray point source in the 95% confidence error box of the gamma-ray source is spatially coincident with the pulsar within the 1' uncertainty (Pivovaroff, Kaspi & Gotthelf 1999). The double peaked gamma-ray pulse morphology and leading radio pulse are similar to those seen for other gamma-ray pulsars and are well-explained in models in which the gamma-ray emission is produced in charge-depleted gaps in the outer magnetosphere. The inferred pulsed gamma-ray flux above 400 MeV, (2.5 +/- 0.6) x 10E-10 erg/cm^2/s, represents 0.011 +/- 0.003 of the pulsar's spin-down luminosity, for a distance of 3 kpc and 1 sr beaming. For PSR J1105-6107, light curves obtained by folding EGRET photons using contemporaneous radio ephemerides show no significant features. We conclude that this pulsar converts less than 0.014 of its spin-down luminosity into E > 100 MeV gamma-rays beaming in our direction (99% confidence), assuming a distance of 7 kpc, 1 sr beaming and a duty cycle of 0.5.

We present the results of a multiwavelength campaign for Mrk 501 performed in March 1996 with ASCA, EGRET, Whipple, and optical telescopes. In the X-ray band, a spectral break was observed around 2 keV. We report here for the first time the detection of high-energy gamma-ray flux from Mrk 501 with EGRET with 3.5 sigma significance (E>100 MeV). Higher flux was also observed in April/May 1996, with 4.0 sigma significance for E>100 MeV, and 5.2 sigma significance for E>500 MeV. The gamma-ray spectrum was measured to be flatter than most of the gamma-ray blazars. We find that the multiband spectrum in 1996 is consistent with that calculated from a one-zone SSC model in a homogeneous region. In the context of this model, we investigate the values of the magnetic field strength and the beaming factor allowed by the observational results. We compare the March 1996 multiwavelength spectrum with that in the flare state in April 1997. Between these two epochs, the TeV flux increase is well correlated with that observed in keV range. The keV and TeV amplitudes during the April 1997 flare are accurately reproduced, assuming that the population of synchrotron photons in 1996 are scattered by the newly injected relativistic electrons, having maximum energies of G_max = 6.0e6. However, the TeV spectrum observed during March 1996 campaign is flatter than predicted by our models. We find that this cannot be explained by either higher order Comptonization or the contribution of the `seed' IR photons from the host galaxy for the first-order external radiation Comptonization.

We report the discovery of a remarkable gamma-ray transient source near the Galactic plane, GRO J1838-04. This source was serendipitously discovered by EGRET in June 1995 with a peak intensity of 4 E(-6) ph/cm/cm/s (for photon energies larger than 100 MeV) and a 5.9 sigma significance. At that time, GRO J1838-04 was the second brightest gamma-ray source in the sky. A subsequent EGRET pointing in late September 1995 detected the source at a flux smaller than its peak value by a factor of 7. We determined that no radio-loud spectrally-flat blazar is within the error box of GRO J1838-04. We discuss the origin of the \ggg-ray transient source and show that interpretations in terms of AGNs or isolated pulsars are highly problematic. GRO J1838-04 provides strong evidence for the existence of a new class of variable gamma-ray sources.

We present a method to assess the reliability of the identification of EGRET sources with extragalactic radio sources. We verify that EGRET is detecting the blazar class of AGN. However, many published identifications are found to be questionable. We provide a table of 42 blazars which we expect to be robust identifications of EGRET sources. This includes one previously unidentified EGRET source, the lensed AGN PKS 1830-210 near the direction of the Galactic center. We provide the best available positions for 16 more radio sources which are also potential identifications for previously unidentified EGRET sources. All high Galactic latitude EGRET sources (b>3 degrees) which demonstrate significant variability can be identified with flat spectrum radio sources. This suggests that EGRET is not detecting any type of AGN other than blazars. This identification method has been used to establish with 99.998% confidence that the peak gamma-ray flux of a blazar is correlated with its average 5 GHz radio flux. An even better correlation is seen between gamma-ray flux and the 2.29 GHz flux density measured with VLBI at the base of the radio jet. Also, using high confidence identifications, we find that the radio sources identified with EGRET sources have larger correlated VLBI flux densities than the parent population of flat radio spectrum sources.

We report the results of a CGRO multi-instrument series of observations of the unidentified gamma-ray source 2CG 135+1 and of its possible counterpart, the peculiar radio source GT 0236+610 coincident with the Be star LSI+61 303. We monitored the time variable radio source GT 0236+610 during the CGRO observations. OSSE and COMPTEL data were obtained during the period May-June 1994, and BATSE data for the period April 1994-January 1995. We discuss the time variability of the gamma-ray emission and spectral properties of 2CG~135+1. Understanding the nature of 2CG 135+1 may be of crucial importance for the interpretation of a class of unidentified time variable gamma-ray sources in the Galaxy.

We report the results of a CGRO 3-week observation of the binary system containing the 47 ms pulsar PSR B1259-63 orbiting around a Be star companion in a very eccentric orbit. The PSR B1259-63 binary is a unique system for the study of the interaction of a rapidly rotating pulsar with time-variable nebular surroundings. CGRO observed the PSR B1259-63 system in coincidence with its most recent periastron passage (January 3-23, 1994). Unpulsed and non-thermal hard X-ray emission was detected up to 200 keV, with a photon index $1.8 \pm 0.2$ and a flux of ~4 mCrab, corresponding to a luminosity of a few 10^{34} erg/s at the distance of 2 kpc. The hard X-ray flux and spectrum detected by CGRO agrees with the X-ray emission measured by simultaneous ASCA observations. EGRET upper limits are significant, and exclude strong inverse Compton cooling in the PSR B1259-63 system. We interpret the observed non-thermal emission as synchrotron radiation of shocked electron/positron pairs of the relativistic pulsar wind interacting with the mass outflow from the Be star. Our results clearly indicate, for the first time in a binary pulsar, that high energy emission can be shock-powered rather than caused by accretion. The lack of X-ray/gamma-ray pulsations constrains models of high-energy emission from rapidly rotating pulsars.

We report the observation by the Compton Gamma Ray Observatory of a spectacular flare of radio source PKS 1622-297. A peak flux of 17E-6 cm^-2 s^-1 (E > 100 MeV) was observed. The corresponding isotropic luminosity is 2.9E49 erg/s. We find that PKS 1622-297 exhibits gamma-ray intra-day variability. A flux increase by a factor of at least 3.6 was observed to occur in less than 7.1 hours (with 99% confidence). Assuming an exponential rise, the corresponding doubling time is less than 3.8 hours. A significant flux decrease by a factor of ~2 in 9.7 hours was also observed. Without beaming, the rapid flux change and large isotropic luminosity are inconsistent with the Elliot-Shapiro condition (assuming that gas accretion is the immediate source of power for the gamma-rays). This inconsistency suggests that the gamma-ray emission is beamed. A minimum Doppler factor of 8.1 is implied by the observed lack of pair-production opacity (assuming x-rays are emitted co-spatially with the gamma-rays). Simultaneous observation by EGRET and OSSE finds a spectrum adequately fit by a power law with photon index of -1.9. Although the significance is not sufficient to establish this beyond doubt, the high-energy gamma-ray spectrum appears to evolve from hard to soft as a flare progresses.