We present a method for enhancing the sensitivity of X-ray telescopic observations with imaging polarimeters, with a focus on the gas pixel detectors (GPDs) to be flown on the Imaging X-ray Polarimetry Explorer (IXPE). Our analysis measures photoelectron directions, X-ray absorption points and X-ray energies for 2-8keV event tracks, with estimates for both the statistical and systematic (reconstruction) uncertainties. We use a weighted maximum likelihood combination of predictions from a deep ensemble of ResNet convolutional neural networks, trained on Monte Carlo event simulations. We define a figure of merit to compare the polarization bias-variance trade-off in track reconstruction algorithms. For power-law source spectra, our method improves on current state-of-the-art (and previous deep learning approaches), providing ~45% increase in effective exposure times. For individual energies, our method produces 20-30% absolute improvements in modulation factor for simulated 100% polarized events, while keeping residual systematic modulation within 1 sigma of the finite sample minimum. Absorption point location and photon energy estimates are also significantly improved. We have validated our method with sample data from real GPD detectors.

We report on new X-ray and optical observations of PSR J2030+4415, a Gamma-ray Pulsar with an H$\alpha$ bow shock. These data reveal the velocity structure of the bow shock apex and resolve unusual X-ray structure in its interior. In addition the system displays a very long, thin filament, extending at least $5^\prime$ at $\sim 130^\circ$ to the pulsar motion vector. Careful astrometry, compared with a short archival exposure, detects the pulsar proper motion at 85 mas yr$^{-1}$. With the H$\alpha$ velocity structure this allows us to estimate the distance as $0.75$ kpc.

We report on orbital modulation of the 100-600 MeV gamma-ray emission of the $P_{\rm B}=4.6$ hr millisecond pulsar binary PSR J2339$-$0533 using 11 yr of Fermi Large Area Telescope data. The modulation has high significance (chance probability $p\approx 10^{-7}$), is approximately sinusoidal, peaks near pulsar superior conjunction, and is detected only in the low-energy 100-600 MeV band. The modulation is confined to the on-pulse interval, suggesting that the variation is in the 2.9-ms pulsed signal itself. This contrasts with the few other known systems exhibiting GeV orbital modulations, as these are unpulsed and generally associated with beamed emission from an intrabinary shock. The origin of the modulated pulsed signal is not yet clear, although we describe several scenarios, including Compton upscattering of photons from the heated companion. This would require high coherence in the striped pulsar wind.

We present a model for atmospheric wind circulation in binary millisecond pulsar (MSP) companions, showing how the optical light curve (LC) and radial velocities (RV) are sensitive to the wind flow, causing LC orbital phase shifts and asymmetries, as observed for several `spider' MSP. Velocity widths of spectral lines offer additional opportunities for measuring surface wind speed. As examples, we fit optical data for a black widow pulsar J1959+2048 and a redback pulsar J2215+5135; the wind heating models (WH) are statistically strongly preferred over direct heating (DH) for both objects, although the latter is even better fit with a heated spot. In general, WH effects tend to increase the inferred orbital inclination $i$ and decrease the inferred companion center-of-mass radial velocity amplitude $K_\mathrm{c}$; both effects decrease the inferred neutron star mass. Even with such a decrease, we find large masses for the two neutron stars: $2.18\pm0.09M_\odot$ and $2.28^{+0.10}_{-0.09}M_\odot$, respectively (for the modest surface speeds fit from the bulk heat flow; supersonic photospheric winds can slightly change these values). These are among the highest masses known, and our improved modeling increases confidence that the results are important for understanding the dense matter equation of state.

Since the birth of X-ray Astronomy, spectacular advances have been seen in the imaging, spectroscopic and timing studies of the hot and violent X-ray Universe, and further leaps forward are expected in the future. On the other hand, polarimetry is very much lagging behind: after the measurements of the Crab Nebula and Scorpius X-1, obtained by OSO-8 in the 70s, no more observations have been performed in the classical X-ray band, even if some interesting results have been obtained in hard X-rays and in soft gamma-rays. The NASA/ASI mission IXPE, scheduled for the launch in 2021, is going to provide for the first time imaging X-ray polarimetry in the 2-8 keV band thanks to its photoelectric polarimeter, coupled with ~25'' angular resolution X-ray mirrors. Its orders of magnitude improvement in sensitivity with respect to the OSO-8 Bragg polarimeter implies scientifically meaningful polarimetric measurements for at least the brightest specimens of most classes of X-ray sources. In 2027, the Chinese-led mission eXTP should also be launched. In addition to timing and spectroscopic instruments, eXTP will have on board photoelectric polarimeters very similar to those of IXPE, but with a total effective area 2-3 times larger. Building on IXPE results, eXTP will increase the number of sources for which significant polarimetric measurements could be obtained. However, further progresses, such as exploring a broader energy range, considering a larger effective area, improving the angular resolution, and performing wide-field polarization measurements, are needed to reach a mature phase for X-ray polarimetry. In the first part of this White Paper we will discuss a few scientific cases in which a next generation X-ray Polarimetry mission can provide significant advances. In the second part, a possible concept for a medium-class Next Generation X-ray Polarimetry (NGXP) mission will be sketched.

We describe a geometric model for synchrotron and synchrotron self-Compton (SSC) radiation from blazar jets, involving multiple emission zones with turbulent magnetic fields and fully self-consistent seed photon mixing for SSC. Including the effects of jet divergence, particle cooling and the Relativistic PA rotation (RPAR) to the observer frame, we find that the multi-zone model recovers simple predictions for SSC polarization, but describes new dependencies on jet viewing geometry and zone multiplicity. Increasing the zone number decreases both synchrotron and SSC polarization, but with different scaling. A rise in synchrotron polarization fraction $\Pi_{\rm Sync}$ at high energies is guaranteed by basic relativity considerations, and strengthened by jet non-uniformity. Finite light travel time effects can suppress the synchrotron polarization at energies well below the $\nu_{\rm Sync}$ peak. In general $\Pi_{\rm Sync}$ and $\Pi_{\rm SSC}$ are correlated with $\Pi_{\rm SSC}/\Pi_{\rm Sync}\approx 0.3$, but individual realizations can lie far from this trend. This study lets us estimate $\Pi$ across the SED, leading to predictions in the X-ray band helpful for planning observations with {\it IXPE} and other upcoming X-ray polarization missions.

X-ray polarization should provide new probes of magnetic field geometry and acceleration physics near the base of blazar jets, but near-future missions will have limited sensitivity. We thus use existing lower energy data and X-ray variability measurements in the context of a basic synchro-Compton model to predict the X-ray polarization level and the probability of detection success for individual sources, listing the most attractive candidates for an IXPE campaign. We find that, as expected, several high-peak blazars such as Mrk 421 can be easily measured in 100 ksec exposures. Most low peak sources should only be accessible to triggered campaigns during bright flares. Surprisingly, a few intermediate peak sources can have anomalously high X-ray polarization and thus are attractive targets.

We collect new and archival optical observations of nine "black-widow" millisecond pulsar binaries. New measurements include direct imaging with the Keck, Gemini-S, MDM, and LCO 2~m telescopes. This is supplemented by synthesized colors from Keck long-slit spectra. Four black-widow optical companions are presented here for the first time. Together these data provide multicolor photometry covering a large fraction of the orbital phase. We fit these light curves with a direct (photon) heating model using a version of the ICARUS light-curve modeling code. The fits provide distance and fill-factor estimates, inclinations, and heating powers. We compare the heating powers with the observed GeV luminosities, noting that the ratio is sensitive to pulsar distance and to the gamma-ray beaming. We make a specific correction for "outer-gap" model beams, but even then some sources are substantially discrepant, suggesting imperfect beaming corrections and/or errors in the fit distance. The fits prefer large metal abundance for half of the targets, a reasonable result for these wind-stripped secondaries. The companion radii indicate substantial Roche-lobe filling, $f_c \approx 0.7-1$ except for PSR J0952$-$0607, which with $f_c< 0.5$ has a companion density $\rho \approx 10\,{\rm g\,cm^{-3}}$, suggesting unusual evolution. We note that the direct-heating fits imply large heating powers and rather small inclinations, and we speculate that unmodeled effects can introduce such bias.

We model millisecond pulsars winds colliding with radiatively-driven companion winds in black widow and redback systems. For the redbacks, the geometry of this intrabinary shock (IBS) is quite sensitive to the expected equatorial concentration in the pulsar outflow. We thus analytically extend IBS thin-shock models to $\sim \sin^{2n}\theta$ pulsar winds. We compute the synchrotron emission from such shocks, including the build-up and cooling of the particle population as it accelerates along the IBS. For reasonable parameters, this IBS flux dominates the binary emission in the X-ray band. The modeling shows subtle variation in spectrum across the IBS peak, accessible to sensitive X-ray studies. As example applications, we fit archival CXO/XMM data from the black widow pulsar J1959+2048 and the redback PSR J2339-0533, finding that the model reproduces well the orbital light curve profiles and energy spectra. The results show a very hard injected electron spectrum, indicating likely dominance by reconnection. The light curve fitting is sensitive to the geometric parameters, including the very important orbital inclination $i$. Coupled with optical fits of the companion star, such IBS X-ray light curve modeling can strongly constrain the binary geometry and the energetics of the MSP wind.

Even with several thousand {\it Fermi}-LAT blazar detections, the $\gamma$-ray emission mechanism is poorly understood. We explore correlated optical/$\gamma$-ray flux variations for 178 {\it Fermi}-LAT blazars regularly monitored by KAIT, SMARTS, and the Steward Observatory. Out of the 178 sources, 121 show a measurable ($>1\sigma$) discrete correlation function peak. Using the derived time-lags and Bayesian block light-curve decompositions, we measure the fraction of common and orphan flares between the two bands. After accounting for sampling and sensitivity limitations we quantify for the first time the true orphan flare rates of optical and $\gamma$-ray flares: 54.5\% of optical and 20\% of $\gamma$-ray flares are orphan events. Both the intraband temporal relation and the small orphan $\gamma$-ray flare fraction point toward leptonic processes as the likely mechanism for the high-energy emission. Motivated to discriminate between synchrotron self-Compton and external-Compton dominance in individual sources, we use the flux-flux variations to determine the slope $m$ of the ${\rm log}\,f_{\rm opt} - {\rm log}\,f_{\gamma}$ dependence. The slope distribution suggests a bimodal population with high and intermediate synchrotron-peak objects showing larger $m$ than low synchrotron-peak objects. We find that $m$ is naturally decreased through pollution from the orphan (typically optical) flares and develop a method to statistically recover, given the sources' measured orphan flare rate, the intrinsic $m$. While source classes show composite behavior, the majority of BL Lac objects favor $m=2$, indicating a synchrotron self-Compton origin for the $\gamma$-rays. No preference for either $m$ is found in flat spectrum radio quasars.

We present the first limits on inelastic electron-scattering dark matter and dark photon absorption using a prototype SuperCDMS detector having a charge resolution of 0.1 electron-hole pairs (CDMS HVeV, a 0.93 gram CDMS HV device). These electron-recoil limits significantly improve experimental constraints on dark matter particles with masses as low as 1 MeV/$\mathrm{c^2}$. We demonstrate a sensitivity to dark photons competitive with other leading approaches but using substantially less exposure (0.49 gram days). These results demonstrate the scientific potential of phonon-mediated semiconductor detectors that are sensitive to single electronic excitations.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

We report on high-energy properties of the black widow pulsar PSR J2241$-$5236 in the X-ray and the Fermi-LAT (GeV gamma-ray) bands. In the LAT band, the phase-averaged gamma-ray light curve shows orbital modulation below $\sim$1 GeV with a chance probability ($p$) monotonically decreasing with time to $p\sim 10^{-5}$. The peak of the light curve is near the superior conjunction of the pulsar (binary phase $\phi_{\rm B}\approx 0.25$). We attribute the modulation to the intra-binary shock (IBS) emission and search for IBS signatures in the archival X-ray data. We find that the X-ray spectral fit requires a non-thermal component, which implies a possible IBS origin of the X-rays. We discuss our observations in the context of IBS scenarios.

Blazars are known for their energetic multiwavelength flares from radio wavelengths to high-energy $\gamma$-rays. In this work, we study radio, optical, and $\gamma$-ray light curves of 145 bright blazars spanning up to 8~yr, to probe the flaring activity and interband correlations. Of these, 105 show $>1\sigma$ correlations between one or more wavebands, 26 of which have a $>3\sigma$ correlation in at least one wavelength pair, as measured by the discrete correlation function. The most common and strongest correlations are found between the optical and $\gamma$-ray bands, with fluctuations simultaneous within our $\sim 30$~d resolution. The radio response is usually substantially delayed with respect to the other wavelengths with median time lags of $\sim 100$--160~d. A systematic flare identification via Bayesian block analysis provides us with a first uniform sample of flares in the three bands, allowing us to characterise the relative rates of multiband and "orphan" flares. Multiband flares tend to have higher amplitudes than "orphan" flares.

We describe a geometric model for synchrotron radiation from blazar jets, involving multiple emission zones with turbulent magnetic fields and a transient core with a helical B field. Including the effects of jet divergence, particle cooling and the Relativistic PA rotation (RPAR) to the observer frame, we find polarization behavior consistent with recent data from monitoring campaigns. We predict that under some circumstances multi-$\pi$ rotation phases should exhibit relativistically-induced steps in rate $d{\rm PA}/dt$ and modulation in polarization $\Pi$ that can be helpful in pinning down the jet $\Gamma$ and $\theta_{\rm obs}$. Also, RPAR enhances waveband differences that will be particularly interesting for comparing radio, optical and, soon, X-ray PA and $\Pi$ variations.

We present an analysis of archival Gemini g', r', K and Keck H, K_s imaging of this nearby short period binary (P_B=95.8min) 2.87ms pulsar. The heated companion is clearly detected. Direct pulsar heating provides an acceptable model at the revised >700pc parallax distance. The relatively shallow light curve modulation prefers an inclination i<40deg; this high latitude view provides a likely explanation for the lack of radio signatures of wind dispersion or eclipse. It also explains the low minimum companion mass and, possibly, the faintness of the source in X- and gamma-rays.

We report on Gemini, NuSTAR and 8-year Fermi observations of the most distant blazar QSO~J0906$+$6930 ($z=5.48$). We construct a broadband spectral energy distribution (SED) and model the SED using a synchro-Compton model. The measurements find a $\sim 4 \times 10^9 M_\odot$ mass for the black hole and a spectral break at $\sim$4 keV in the combined fit of the new NuSTAR and archival Chandra data. The SED fitting constrains the bulk Doppler factor $\delta$ of the jet to $9^{+2.5}_{-3}$ for QSO~J0906$+$6930. Similar, but weaker $\delta$ constraints are derived from SED modeling of the three other claimed $z>5$ blazars. Together, these extrapolate to $\sim620$ similar sources, fully 20\% of the optically bright, high mass AGN expected at $5<z<5.5$. This has interesting implications for the early growth of massive black holes.

We assess the science reach and technical feasibility of a satellite mission based on precision atomic sensors configured to detect gravitational radiation. Conceptual advances in the past three years indicate that a two-satellite constellation with science payloads consisting of atomic sensors based on laser cooled atomic Sr can achieve scientifically interesting gravitational wave strain sensitivities in a frequency band between the LISA and LIGO detectors, roughly 30 mHz to 10 Hz. The discovery potential of the proposed instrument ranges from from observation of new astrophysical sources (e.g. black hole and neutron star binaries) to searches for cosmological sources of stochastic gravitational radiation and searches for dark matter.

We describe new measurements of the remarkable H$\alpha$/UV/X-ray bowshock and pulsar wind nebula of the isolated millisecond pulsar PSR J2124$-$3358. {\it CXO} imaging shows a one-sided jet structure with a softer equatorial outflow. KOALA IFU spectroscopy shows that non-radiative emission dominates the bow shock and that the H$\alpha$ nebula is asymmetric about the pulsar velocity with an elongation into the plane of the sky. We extend analytic models of the contact discontinuity to accommodate such shapes and compare these to the data. Using {\it HST} UV detections of the pulsar and bow shock, radio timing distance and proper motion measurements and the {\it CXO}-detected projected spin axis we model the 3-D PWN momentum flux distribution. The integrated momentum flux depends on the ionization of the ambient ISM, but for an expected ambient WNM we infer $I = 2.4 \times 10^{45} {\rm g\,cm^2}$. This implies $M_{NS} = 1.6-2.1 M_\odot$, depending on the equation of state, which in turn suggests that the MSP gained significant mass during recycling and then lost its companion. However, this conclusion is at present tentative, since lower ionization allows $\sim 30\%$ lower masses and uncertainty in the parallax allows up to 50% error.

We have studied the variability of the black-widow type binary millisecond pulsar PSR J1311$-$3430 from optical to gamma-ray energies. We confirm evidence for orbital modulation in the weak off-pulse $\ge$200-MeV emission, with a peak at $\phi_B\approx 0.8$, following pulsar inferior conjunction. The peak has a relatively hard spectrum, extending above $\sim 1$\,GeV. XMM-Newton and Swift UV observations also show that this source's strong X-ray flaring activity is associated with optical/UV flares. With a duty cycle $\sim 7-19$%, this flaring is quite prominent with an apparent power-law intensity distribution. Flares are present at all orbital phases, with a slight preference for $\phi_B=0.5-0.7$. We explore possible connections of these variabilities with the intrabinary shock and magnetic activity on the low mass secondary.

The companions of evaporating binary pulsars (black widows and related systems) show optical emission suggesting strong heating. In a number of cases large observed temperatures and asymmetries are inconsistent with direct radiative heating for the observed pulsar spindown power and expected distance. Here we describe a heating model in which the pulsar wind sets up an intrabinary shock (IBS) against the companion wind and magnetic field, and a portion of the shock particles duct along this field to the companion magnetic poles. We show that a variety of heating patterns, and improved fits to the observed light curves, can be obtained at expected pulsar distances and luminosities, at the expense of a handful of model parameters. We test this `IBS-B' model against three well observed binaries and comment on the implications for system masses.

The radio galaxy 0402+379 is believed to host a supermassive black hole binary (SMBHB). The two compact core sources are separated by a projected distance of 7.3 pc, making it the most (spatially) compact resolved SMBHB known. We present new multi-frequency VLBI observations of 0402+379 at 5, 8, 15 and 22 GHz, and combine with previous observations spanning 12 years. A strong frequency dependent core shift is evident, which we use to infer magnetic fields near the jet base. After correcting for these shifts we detect significant relative motion of the two cores at $\beta=v/c=0.0054 \pm 0.0003$ at $PA= -34.4^\circ$. With some assumptions about the orbit, we use this measurement to constrain the orbital period $P\approx 3 \times 10^4$ y and SMBHB mass $M \approx 15 \times 10^9\ M_\odot$. While additional observations are needed to confirm this motion and obtain a precise orbit, this is apparently the first black hole system resolved as a visual binary.

We present broadband spectral energy distributions (SEDs) and light curves of the gamma-ray binary 1FGL J1018.6$-$5856 measured in the X-ray and the gamma-ray bands. We find that the orbital modulation in the low-energy gamma-ray band is similar to that in the X-ray band, suggesting a common spectral component. However, above a GeV the orbital light curve changes significantly. We suggest that the GeV band contains significant flux from a pulsar magnetosphere, while the X-ray to TeV light curves are dominated by synchrotron and Compton emission from an intrabinary shock (IBS). We find that a simple one-zone model is inadequate to explain the IBS emission, but that beamed Synchrotron-self Compton radiation from adiabatically accelerated plasma in the shocked pulsar wind can reproduce the complex multiband light curves, including the variable X-ray spike coincident with the gamma-ray maximum. The model requires inclination $\sim$50$^\circ$ and orbital eccentricity $\sim$0.35, consistent with the limited constraints from existing optical observations. This picture motivates searches for pulsations from the energetic young pulsar powering the wind shock.

We report on six new Chandra observations of the Geminga pulsar wind nebula (PWN). The PWN consists of three distinct elongated structures - two $\approx 0.2 d_{250}$ pc long lateral tails and a segmented axial tail of $\approx 0.05 d_{250}$ pc length, where $d_{250}=d/(250 {\rm pc})$. The photon indices of the power law spectra of the lateral tails, $\Gamma \approx 1$, are significantly harder than those of the pulsar ($\Gamma \approx 1.5$) and the axial tail ($\Gamma \approx 1.6$). There is no significant diffuse X-ray emission between the lateral tails -- the ratio of the X-ray surface brightness between the south tail and this sky area is at least 12. The lateral tails apparently connect directly to the pulsar and show indication of moving footpoints. The axial tail comprises time-variable emission blobs. However, there is no evidence for constant or decelerated outward motion of these blobs. Different physical models are consistent with the observed morphology and spectra of the Geminga PWN. In one scenario, the lateral tails could represent an azimuthally asymmetric shell whose hard emission is caused by the Fermi acceleration mechanism of colliding winds. In another scenario, the lateral tails could be luminous, bent polar outflows, while the blobs in the axial tail could represent a crushed torus. In a resemblance to planetary magnetotails, the blobs of the axial tail might also represent short-lived plasmoids which are formed by magnetic field reconnection in the relativistic plasma of the pulsar wind tail.

We report on Chandra X-ray Observatory (CXO) observations of the pulsar wind nebula (PWN) associated with PSR B0355+54 (eight observations with a 395 ks total exposure, performed over an 8 month period). We investigated the spatial and spectral properties of the emission coincident with the pulsar, compact nebula (CN), and extended tail. We find that the CN morphology can be interpreted in a way that suggests a small angle between the pulsar spin axis and our line-of-sight, as inferred from the radio data. On larger scales, emission from the 7' (2 pc) tail is clearly seen. We also found hints of two faint extensions nearly orthogonal to the direction of the pulsar's proper motion. The spectrum extracted at the pulsar position can be described with an absorbed power-law + blackbody model. The nonthermal component can be attributed to magnetospheric emission, while the thermal component can be attributed to emission from either a hot spot (e.g., a polar cap) or the entire neutron star surface. Surprisingly, the spectrum of the tail shows only a slight hint of cooling with increasing distance from the pulsar. This implies either a low magnetic field with fast flow speed, or particle re-acceleration within the tail. We estimate physical properties of the PWN and compare the morphologies of the CN and the extended tail with those of other bow shock PWNe observed with long CXO exposures.

The low mass companions of evaporating binary pulsars (black widows and their ilk) are strongly heated on the side facing the pulsar. However in high-quality photometric and spectroscopic data the heating pattern does not match that expected for direct pulsar illumination. Here we explore heating mediated by an intra-binary shock (IBS). We develop a simple analytic model and implement it in the popular `ICARUS' light curve code. The model is parameterized by the wind momentum ratio beta and velocity v_Rel v_orb and assumes that the reprocessed pulsar wind emits prompt particles or radiation to heat the companion surface. We illustrate an interesting range of light curve asymmetries controlled by these parameters. The code also computes the IBS synchrotron emission pattern, and thus can model black widow X-ray light curves. As a test we apply the results to the high quality asymmetric optical light curves of PSR J2215+5135; the resulting fit gives a substantial improvement upon direct heating models and produces an X-ray light curve consistent with that seen. The IBS model parameters imply that, at the present loss rate, the companion evaporation has a characteristic timescale tau_evap ~150My. Still, the model is not fully satisfactory, indicating additional unmodeled physical effects.

We analyzed deep $Chandra$ ACIS-I exposures of the cluster-scale X-ray halo surrounding the radio source 4C+37.11. This remarkable system hosts the closest resolved pair of super-massive black hole and an exceptionally luminous elliptical galaxy, the likely product of a series of past mergers. We characterize the halo with $r_{500} = 0.95$ Mpc, $M_{500} = (2.5 \pm 0.2) \times 10^{14} \ M_{\rm{\odot}}$, $ kT = 4.6\pm 0.2$ keV, and a gas mass of $M_{\rm g,500} = (2.2 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.09 \pm 0.01$. The entropy profile shows large non-gravitational heating in the central regions. We see several surface brightness jumps, associated with substantial temperature and density changes, but approximate pressure equilibrium, implying that these are sloshing structures driven by a recent merger. A residual intensity image shows core spiral structure closely matching that seen for the Perseus cluster, although at $z=0.055$ the spiral pattern is less distinct. We infer the most recent merger occurred $1-2$ Gyr ago and that the event that brought the two observed super-massive black holes to the system core is even older. Under that interpretation, this black hole binary pair has, unusually, remained at pc-scale separation for more than 2 Gyr.

We present multi-frequency Very Long Baseline Array (VLBA) follow-up observations of VLBA Imaging and Polarimetry Survey sources identified as likely compact symmetric objects (CSOs) or super-massive binary black holes (SBBHs). We also present new spectroscopic redshifts for 11 sources observed with the Hobby-Eberly Telescope. While no new SBBHs can be confirmed from these observations, we have identified 24 CSOs in the sample, 15 of which are newly designated, and refuted 52 candidates leaving 33 unconfirmed candidates. This is the first large uniform sample of CSOs which can be used to elicit some of the general properties of these sources, including morphological evolution and environmental interaction. We have detected polarised emission from two of these CSOs the properties of which are consistent with Active Galactic Nuclei unification schemes.

We have identified an optical/X-ray binary with orbital period P_b=5.47h as the likely counterpart of the Fermi source 2FGL J2039.6-5620. GROND, SOAR and DES observations provide an accurate orbital period and allow us to compare with the light curve of an archival XMM exposure. Like many short-period optical X-ray binaries associated with LAT sources this may be a interacting (black widow/redback) millisecond pulsar binary. The X-ray light curve is consistent with the emission associated with an intrabinary shock. The optical light curve shows evidence of companion heating, but has a peculiar asymmetric double peak. The nature of this optical structure is not yet clear; additional optical studies and, especially, detection of an orbital modulation in a gamma-ray pulsar are needed to elucidate the nature of this peculiar source.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10{\deg}), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

We present Keck spectroscopic measurements of the millisecond pulsar binary J2215+5135. These data indicate a neutron-star (NS) mass M_NS=1.6Mo, much less than previously estimated. The pulsar heats the companion face to T_D~9000K; the large heating efficiency may be mediated by the intrabinary shock dominating the X-ray light curve. At the best-fit inclination i=88.8deg, the pulsar should be eclipsed. We find weak evidence for such eclipses in the pulsed gamma-rays; an improved radio ephemeris allows use of up to 5 times more Fermi-LAT gamma-ray photons for a definitive test of this picture. If confirmed, the gamma-ray eclipse provides a novel probe of the dense companion wind and the pulsar magnetosphere.

We observed the Fermi-discovered gamma-ray binary 1FGL J1018.6-5856 at 20 epochs over 50 days using the CHIRON spectrograph, obtaining spectra at R~25,000 covering 4090-8908A. The average spectrum confirms an O6 V((f)) spectral type and extinction E(B-V) = 1.35+/-0.04. Variable absorption line equivalent widths suggest substantial contamination by wind line features. The limited S/N ratio hindered accurate continuum definition and prevented measurement of a high quality radial velocity curve. Nevertheless, the best data indicate a radial velocity amplitude <40 km/s for the He II lines and substantially lower for H I. We argue that this indicates a most likely compact object mass <2.2Msun. While black hole solutions are not excluded, a neutron star source of the gamma-ray emission seems preferred.

One of the most intriguing results from the gamma-ray instruments in orbit has been the detection of powerful flares from the Crab Nebula. These flares challenge our understanding of pulsar wind nebulae and models for particle acceleration. We report on the portion of a multiwavelength campaign using Keck, HST, and Chandra concentrating on a small emitting region, the Crab's inner knot, located a fraction of an arcsecond from the pulsar. We find that the knot's radial size, tangential size, peak flux, and the ratio of the flux to that of the pulsar are correlated with the projected distance of the knot from the pulsar. A new approach, using singular value decomposition for analyzing time series of images, was introduced yielding results consistent with the more traditional methods while some uncertainties were substantially reduced. We exploit the characterization of the knot to discuss constraints on standard shock-model parameters that may be inferred from our observations assuming the inner knot lies near to the shocked surface. These include inferences as to wind magnetization, shock shape parameters such as incident angle and poloidal radius of curvature, as well as the IR/optical emitting particle enthalpy fraction. We find that while the standard shock model gives good agreement with observation in many respects, there remain two puzzles: (a) The observed angular size of the knot relative to the pulsar--knot separation is much smaller than expected; (b) The variable, yet high degree of polarization reported is difficult to reconcile with a highly relativistic downstream flow.

We report on a series of spectroscopic observations of PSR J1311-3430, an extreme black-widow gamma-ray pulsar with a helium-star companion. In a previous study we estimated the neutron star mass as M_NS= 2.68+/-0.14M_Sun (statistical error), based on limited spectroscopy and a basic (direct heating) light curve model; however, much larger model-dependent systematics dominate the mass uncertainty. Our new spectroscopy reveals a range of complex source behavior. The variable He I companion wind emission lines can dominate broad-band photometry, especially in red filters or near minimum brightness, and the wind flux should complete companion evaporation in a spin-down time. The heated companion face also undergoes dramatic flares, reaching 40,000K over 20% of the star; this is likely powered by a magnetic field generated in the companion. The companion center-of-light radial velocity is now well measured with K_CoL = 615.4+/-5.km/s. We detect non-sinusoidal velocity components due to the heated face flux distribution. Using our spectra to excise flares and wind lines, we generate substantially improved light curves for companion continuum fitting. We show that the inferred inclination and neutron star mass, however, remain sensitive to the poorly constrained heating pattern. The neutron star's mass, M_NS, is likely less than the direct heating value and could range as low as 1.8M_Sun for extreme equatorial heating concentration. While we cannot yet pin down M_NS, our data imply that an intrabinary shock reprocesses the pulsar emission and heats the companion. Improved spectra and, especially, models that include such shock heating are needed for precise parameter measurement.

We obtained six observations of PSR J1741-2054 using the $Chandra$ ACIS-S detector totaling $\sim$300 ks. By registering this new epoch of observations to an archival observation taken 3.2 years earlier using X-ray point sources in the field of view, we have measured the pulsar proper motion at $\mu =109 \pm 10 {\rm mas\ yr}^{-1}$ in a direction consistent with the symmetry axis of the observed H$\alpha$ nebula. We investigated the inferred past trajectory of the pulsar but find no compelling association with OB associations in which the progenitor may have originated. We confirm previous measurements of the pulsar spectrum as an absorbed power law with photon index $\Gamma$=2.68$\pm$0.04, plus a blackbody with an emission radius of (4.5$^{+3.2}_{-2.5})d_{0.38}$ km, for a DM-estimated distance of $0.38d_{0.38}$ kpc and a temperature of $61.7\pm3.0$ eV. Emission from the compact nebula is well described by an absorbed power law model with a photon index of $\Gamma$ = 1.67$\pm$0.06, while the diffuse emission seen as a trail extending northeast of the pulsar shows no evidence of synchrotron cooling. We also applied image deconvolution techniques to search for small-scale structures in the immediate vicinity of the pulsar, but found no conclusive evidence for such structures.

The origin of the extragalactic $\gamma$-ray background (EGB) has been debated for some time. { The EGB comprises the $\gamma$-ray emission from resolved and unresolved extragalactic sources, such as blazars, star-forming galaxies and radio galaxies, as well as radiation from truly diffuse processes.} This letter focuses on the blazar source class, the most numerous detected population, and presents an updated luminosity function and spectral energy distribution model consistent with the blazar observations performed by the {\it Fermi} Large Area Telescope (LAT). We show that blazars account for 50$^{+12}_{-11}$\,\% of the EGB photons ($>$0.1\,GeV), and that {\it Fermi}-LAT has already resolved $\sim$70\,\% of this contribution. Blazars, and in particular low-luminosity hard-spectrum nearby sources like BL Lacs, are responsible for most of the EGB emission above 100\,GeV. We find that the extragalactic background light, which attenuates blazars' high-energy emission, is responsible for the high-energy cut-off observed in the EGB spectrum. Finally, we show that blazars, star-forming galaxies and radio galaxies can naturally account for the amplitude and spectral shape of the background in the 0.1--820\,GeV range, leaving only modest room for other contributions. This allows us to set competitive constraints on the dark-matter annihilation cross section.

We have been using the 0.76-m Katzman Automatic Imaging Telescope (KAIT) at Lick Observatory to optically monitor a sample of 157 blazars that are bright in gamma rays, being detected with high significance ($\ge 10\sigma$) in one year by the Large Area Telescope (LAT) on the {\it Fermi Gamma-ray Space Telescope}. We attempt to observe each source on a 3-day cadence with KAIT, subject to weather and seasonal visibility. The gamma-ray coverage is essentially continuous. KAIT observations extend over much of the 5-year {\it Fermi} mission for several objects, and most have $>100$ optical measurements spanning the last three years. These blazars (flat-spectrum radio quasars and BL~Lac objects) exhibit a wide range of flaring behavior. Using the discrete correlation function (DCF), here we search for temporal relationships between optical and gamma-ray light curves in the 40 brightest sources in hopes of placing constraints on blazar acceleration and emission zones. We find strong optical--gamma-ray correlation in many of these sources at time delays of $\sim 1$ to $\sim 10$ days, ranging between $-40$ and +30 days. A stacked average DCF of the 40 sources verifies this correlation trend, with a peak above 99% significance indicating a characteristic time delay consistent with 0 days. These findings strongly support the widely accepted leptonic models of blazar emission. However, we also find examples of apparently uncorrelated flares (optical flares with no gamma-ray counterpart and gamma-ray flares with no optical counterpart) that challenge simple, one-zone models of blazar emission. Moreover, we find that flat-spectrum radio quasars tend to have gamma rays leading the optical, while intermediate and high synchrotron peak blazars with the most significant peaks have smaller lags/leads.

The {\gamma}-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse {\gamma}-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

We have identified an optical binary with orbital period P_b=4488s as the probable counterpart of the Fermi source 2FGL J1653.6-0159. Although pulsations have not yet been detected, the source properties are consistent with an evaporating millisecond pulsar binary; this P_b=75min is the record low for a spin-powered system. The heated side of the companion shows coherent radial velocity variations, with amplitude K=666.9+/-7.5 km/s for a large mass function of f(M)=1.60+/-0.05 M_sun. This heating suggests a pulsar luminosity ~3x10^34 erg/s. The colors and spectra show additional hard emission dominating at binary minimum. Its origin is, at present, unclear. This system is similar to PSR J1311-3430, with a low mass H-depleted companion, a dense shrouding wind and, likely, a large pulsar mass.

In order to determine the location of the gamma-ray emission site in blazars, we investigate the time-domain relationship between their radio and gamma-ray emission. Light-curves for the brightest detected blazars from the first 3 years of the mission of the Fermi Gamma-ray Space Telescope are cross-correlated with 4 years of 15GHz observations from the OVRO 40-m monitoring program. The large sample and long light-curve duration enable us to carry out a statistically robust analysis of the significance of the cross-correlations, which is investigated using Monte Carlo simulations including the uneven sampling and noise properties of the light-curves. Modeling the light-curves as red noise processes with power-law power spectral densities, we find that only one of 41 sources with high quality data in both bands shows correlations with significance larger than 3-sigma (AO 0235+164), with only two more larger than even 2.25-sigma (PKS 1502+106 and B2 2308+34). Additionally, we find correlated variability in Mrk 421 when including a strong flare that occurred in July-September 2012. These results demonstrate very clearly the difficulty of measuring statistically robust multiwavelength correlations and the care needed when comparing light-curves even when many years of data are used. This should be a caution. In all four sources the radio variations lag the gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. Continuous simultaneous monitoring over a longer time period is required to obtain high significance levels in cross-correlations between gamma-ray and radio variability in most blazars.

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

We have obtained WIYN and SOAR i' images of BL Lacertae objects (BL Lacs) and used these to detect or constrain the flux of the host galaxy. Under common standard candle assumptions these data provide estimates of, or lower bounds on, the redshift. Our targets are a set of flat-spectrum radio counterparts of high flux Fermi Large Area Telescope (LAT) sources, with sensitive spectral observations showing them to be continuum-dominated BL Lacs. In this sample 5 of 11 BL Lacs yielded significant host detections, with standard candle redshifts z=0.13-0.58. Our estimates and lower bounds are generally in agreement with other redshifts estimates, although our z=0.374 estimate for J0543-5532 implies a significantly sub-luminous host.

We report on a survey for H$\alpha$ bow shock emission around nearby $\gamma$-detected energetic pulsars. This survey adds three Balmer-dominated neutron star bow shocks to the six previously confirmed examples. In addition to the shock around {\it Fermi} pulsar PSR J1741$-$2054, we now report H$\alpha$ structures around two additional $\gamma$-ray pulsars, PSR J2030+4415 and PSR J1509$-$5850. These are the first known examples of H$\alpha$ nebulae with pre-ionization halos. With new measurements, we show that a simple analytic model can account for the angular size and flux of the bow shocks' apices. The latter, in particular, provides a new pulsar probe and indicates large moments of inertia and smaller distances than previously assumed, in several cases. In particular we show that the re-measured PSR J0437$-$4715 shock flux implies $I = (1.7\pm 0.2) \times 10^{45}/(f_{HI} {\rm sin}i) {\rm g\,cm^2}$. We also derive a distance $d\approx 0.72$kpc for the $\gamma$-ray only pulsar PSR J2030+4415 and revised distances for PSRs J1959+2048 (1.4kpc) and J2555+6535 ($\sim 1$kpc), smaller than the conventional DM-estimated values. Finally we report upper limits for 94 additional LAT pulsars. An estimate of the survey sensitivity indicates that for a warm neutral medium filling factor $\phi_{WNM}\sim 0.3$ there should be a total of $\sim$ nine H$\alpha$ bow shocks in our LAT-targeted survey; given that seven such objects are now known, a much larger $\phi_{WNM}$ seems problematic.

This paper describes the data release of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey conducted between 2005 and 2007. Light curves, spectra, classifications, and ancillary data are presented for 10,258 variable and transient sources discovered through repeat ugriz imaging of SDSS Stripe 82, a 300 deg2 area along the celestial equator. This data release is comprised of all transient sources brighter than r~22.5 mag with no history of variability prior to 2004. Dedicated spectroscopic observations were performed on a subset of 889 transients, as well as spectra for thousands of transient host galaxies using the SDSS-III BOSS spectrographs. Photometric classifications are provided for the candidates with good multi-color light curves that were not observed spectroscopically. From these observations, 4607 transients are either spectroscopically confirmed, or likely to be, supernovae, making this the largest sample of supernova candidates ever compiled. We present a new method for SN host-galaxy identification and derive host-galaxy properties including stellar masses, star-formation rates, and the average stellar population ages from our SDSS multi-band photometry. We derive SALT2 distance moduli for a total of 1443 SN Ia with spectroscopic redshifts as well as photometric redshifts for a further 677 purely-photometric SN Ia candidates. Using the spectroscopically confirmed subset of the three-year SDSS-II SN Ia sample and assuming a flat Lambda-CDM cosmology, we determine Omega_M = 0.315 +/- 0.093 (statistical error only) and detect a non-zero cosmological constant at 5.7 sigmas.

4C+37.11, at z=0.055 shows two compact radio nuclei, imaged by VLBI at 7mas separation, making it the closest known resolved super-massive black hole binary (SMBHB). An important question is whether this unique object is young, caught on the way to a gravitational in-spiral and merger, or has `stalled' at 7pc. We describe new radio/optical/X-ray observations of the massive host and its surrounding X-ray halo. These data reveal X-ray/optical channels following the radio outflow and large scale edges in the X-ray halo. These structures are promising targets for further study which should elucidate their relationship to the unique SMBHB core.

We report a 5.4\sigma\ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 \mu s, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density

Fermi has provided the largest sample of gamma-ray selected blazars to date. In this work we use a uniformly selected set of 211 BL Lacertae (BL Lac) objects detected by it Fermi during its first year of operation. We have obtained redshift constraints for 206 out of the 211 BL Lacs in our sample making it the largest and most complete sample of BL Lacs available in the literature. We use this sample to determine the luminosity function of BL Lacs and its evolution with cosmic time. We find that for most BL Lac classes, the evolution is positive with a space density peaking at modest redshift (z~1.2). The low-luminosity, high-synchrotron peaked (HSP) BL Lacs are an exception, showing strong negative evolution, with number density increasing for z$\lesssim$0.5. Since this rise corresponds to a drop-off in the density of flat-spectrum radio quasars (FSRQs), a possible interpretation is that these HSPs represent an accretion-starved end-state of an earlier merger-driven gas-rich phase. We additionally find that the known BL Lac correlation between luminosity and photon spectral index persists after correction for the substantial observational selection effects with implications for the so called `blazar sequence'. Finally, estimating the beaming corrections to the luminosity function, we find that BL Lacs have an average Lorentz factor of $\gamma=6.1^{+1.1}_{-0.8}$, and that most are seen within 10$^{\circ}$ of the jet axis.

We report on Keck spectra of ten Fermi blazars. J0622+3326, previously unobserved, is shown to be a flat-spectrum radio quasar at redshift z=1.062. The others are known BL Lac-type objects that have resisted previous attempts to secure redshifts. Using a photometric monitoring campaign with the 0.76 m Katzman Automatic Imaging Telescope at Lick Observatory, we identified epochs when the relativistic jet emission was fainter than usual, thus triggering the Keck spectroscopy. This strategy gives improved sensitivity to stars and ionized gas in the host galaxy, thereby providing improved redshift constraints for seven of these sources.

The Vela supernova remnant is the closest supernova remnant to Earth containing an active pulsar, the Vela pulsar (PSR B0833-45). This pulsar is the archetype of the middle-aged pulsar class and powers a bright pulsar wind nebula (PWN), Vela-X, spanning a region of 2{\deg} x 3{\deg} south of the pulsar and observed in the radio, X-ray and very high energy gamma-ray domains. The detection of the Vela-X PWN by the Fermi Large Area Telescope (LAT) was reported in the first year of the mission. Subsequently, we have re-investigated this complex region and performed a detailed morphological and spectral analysis of this source using 4 years of Fermi-LAT observations. This study lowers the threshold for morphological analysis of the nebula from 0.8 GeV to 0.3 GeV, allowing inspection of distinct energy bands by the LAT for the first time. We describe the recent results obtained on this PWN and discuss the origin of the newly detected spatial features.

Using the Giant Metrewave Radio Telescope (GMRT) we performed deep observations to search for radio pulsations in the directions of unidentified Fermi Large Area Telescope (LAT) gamma-ray sources. We report the discovery of an eclipsing black-widow millisecond pulsar, PSR J1544+4937, identified with the un-cataloged gamma-ray source Fermi J1544.2+4941. This 2.16 ms pulsar is in a 2.9 hours compact circular orbit with a very low-mass companion (Mc > 0.017 Msun). At 322 MHz this pulsar is found to be eclipsing for 13% of its orbit, whereas at 607 MHz the pulsar is detected throughout the low-frequency eclipse phase. Variations in the eclipse ingress phase are observed, indicating a clumpy and variable eclipsing medium. Moreover, additional short-duration absorption events are observed around the eclipse boundaries. Using the radio timing ephemeris we were able to detect gamma-ray pulsations from this pulsar, confirming it as the source powering the gamma-ray emission.

We report the Fermi Large Area Telescope discovery of gamma-ray pulsations from the 22.7 ms pulsar A in the double pulsar system J0737-3039A/B. This is the first mildly recycled millisecond pulsar (MSP) detected in the GeV domain. The 2.7 s companion object PSR J0737-3039B is not detected in gamma rays. PSR J0737-3039A is a faint gamma-ray emitter, so that its spectral properties are only weakly constrained; however, its measured efficiency is typical of other MSPs. The two peaks of the gamma-ray light curve are separated by roughly half a rotation and are well offset from the radio and X-ray emission, suggesting that the GeV radiation originates in a distinct part of the magnetosphere from the other types of emission. From the modeling of the radio and the gamma-ray emission profiles and the analysis of radio polarization data, we constrain the magnetic inclination $\alpha$ and the viewing angle $\zeta$ to be close to 90$^\circ$, which is consistent with independent studies of the radio emission from PSR J0737-3039A. A small misalignment angle between the pulsar's spin axis and the system's orbital axis is therefore favored, supporting the hypothesis that pulsar B was formed in a nearly symmetric supernova explosion as has been discussed in the literature already.

We present results from our analysis of Chandra X-ray Observatory, W. M. Keck Observatory, and Karl G. Jansky Very Large Array (VLA) images of the Crab Nebula that were contemporaneous with the gamma-ray flare of 2011 April. Despite hints in the X-ray data, we find no evidence for statistically significant variations that pinpoint the specific location of the flares within the Nebula. The Keck observations extend this conclusion to the "inner knot", i.e., the feature within an arcsecond of the pulsar. The VLA observations support this conclusion. We also discuss theoretical implications of the gamma-ray flares and suggest that the most dramatic gamma-ray flares are due to radiation-reaction-limited synchrotron emission associated with sudden, dissipative changes in the current system sustained by the central pulsar.

We report on spectroscopic observations covering most of the 475 BL Lacs in the 2nd Fermi LAT catalog of AGN. Including archival measurements (correcting several erroneous literature values) we now have spectroscopic redshifts for 44% of the BL Lacs. We establish firm lower redshift limits via intervening absorption systems and statistical lower limits via searches for host galaxies for an additional 51% of the sample leaving only 5% of the BL Lacs unconstrained. The new redshifts raise the median spectroscopic z from 0.23 to 0.33 and include redshifts as large as z=2.471. Spectroscopic redshift minima from intervening absorbers have ~ z= 0.70, showing a substantial fraction at large z and arguing against strong negative evolution. We find that detected BL Lac hosts are bright ellipticals with black hole masses M_\bullet ~ 10^{8.5-9}, substantially larger than the mean of optical AGN and LAT Flat Spectrum Radio Quasar samples. A slow increase in M_\bullet with z may be due to selection bias. We find that the power-law dominance of the optical spectrum extends to extreme values, but this does not strongly correlate with the gamma-ray properties, suggesting that strong beaming is the primary cause of the range in continuum dominance.

We report the detection of radio emission from PSR J1311-3430, the first millisecond pulsar discovered in a blind search of Fermi Large Area Telescope (LAT) gamma-ray data. We detected radio pulsations at 2 GHz, visible for <10% of ~4.5-hrs of observations using the Green Bank Telescope (GBT). Observations at 5 GHz with the GBT and at several lower frequencies with Parkes, Nancay, and the Giant Metrewave Radio Telescope resulted in non-detections. We also report the faint detection of a steep spectrum continuum radio source (0.1 mJy at 5 GHz) in interferometric imaging observations with the Jansky Very Large Array. These detections demonstrate that PSR J1311-3430, is not radio quiet and provides additional evidence that the radio beaming fraction of millisecond pulsars is very large. The radio detection yields a distance estimate of 1.4 kpc for the system, yielding a gamma-ray efficiency of 30%, typical of LAT-detected MSPs. We see apparent excess delay in the radio pulsar as the pulsar appears from eclipse and we speculate on possible mechanisms for the non-detections of the pulse at other orbital phases and observing frequencies.

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

We have obtained initial spectroscopic observations and additional photometry of the newly discovered Pb=94min gamma-ray black-widow pulsar PSR J1311-3430. The Keck spectra show a He-dominated, nearly H-free photosphere and a large radial-velocity amplitude of 609.5+/-7.5km/s. Simultaneous seven-color GROND photometry further probes the heating of this companion, and shows the presence of a flaring infrared excess. We have modeled the quiescent light curve, constraining the orbital inclination and masses. Simple heated light-curve fits give M_NS=2.7Msun, but show systematic light-curve differences. Adding extra components allows a larger mass range to be fit, but all viable solutions have M_NS>2.1Msun. If confirmed, such a large M_NS substantially constrains the equation of state of matter at supernuclear densities.

The pulsar wind nebula associated with PSR J1826-1334, HESS J1825-137, is a bright very high energy source with an angular extent of ~1 degree and spatially-resolved spectroscopic TeV measurements. The gamma-ray spectral index is observed to soften with increasing distance from the pulsar, likely the result of cooling losses as electrons traverse the nebula. We describe analysis of X-ray data of the extended nebula, as well as 3-D time-dependent spectral energy distribution modeling, with emphasis on the spatial variations within HESS J1825-137. The multi-wavelength data places significant constraints on electron injection, transport, and cooling within the nebula. The large size and high nebular energy budget imply a relatively rapid initial pulsar spin period of 13 \pm 7 ms and an age of 40 \pm 9 kyr. The relative fluxes of each VHE zone can be explained by advective particle transport with a radially decreasing velocity profile with $v(r) \propto r^{-0.5}$. The evolution of the cooling break requires an evolving magnetic field which also decreases radially from the pulsar, $B(r,t) \propto r^{-0.7} \dot{E}(t)^{1/2}$. Detection of 10 TeV flux ~80 pc from the pulsar requires rapid diffusion of high energy particles with $\tau_{esc} \approx 90 (R / 10 pc)^2 (E_e/100 TeV)^{-1}$ year, contrary to the common assumption of toroidal magnetic fields with strong magnetic confinement. The model predicts a rather uniform Fermi LAT surface brightness out to ~1 degree from the pulsar, in good agreement with the recently discovered LAT source centered 0.5 degree southwest of PSR J1826-1334 with extension 0.6 \pm 0.1 degree.

We report on a new Chandra exposure of PSR J1809-2332, the recently discovered pulsar powering the bright EGRET source 3EG J1809-2328. By registration of field X-ray sources in an archival exposure, we measure a significant proper motion for the pulsar point source over an ~11 year baseline. The shift of 0.30+/-0.06" (at PA= 153.3+/-18.4) supports an association with proposed SNR parent G7.5-1.7. Spectral analysis of diffuse emission in the region also supports the interpretation as a hard wind nebula trail pointing back toward the SNR.

We have found an optical/X-ray counterpart candidate for the bright, but presently unidentified, Fermi source 2FGL J1311.7-3429. This counterpart undergoes large amplitude quasi-sinusoidal optical modulation with a 1.56h (5626s) period. The modulated flux is blue at peak, with T_eff ~14,000K, and redder at minimum. Superimposed on this variation are dramatic optical flares. Archival X-ray data suggest modest binary modulation, but no eclipse. With the gamma-ray properties, this appears to be another black-widow-type millisecond pulsar. If confirmation pulses can be found in the GeV data, this binary will have the shortest orbital period of any known spin-powered pulsar. The flares may be magnetic events on the rapidly rotating companion or shocks in the companion-stripping wind. While this may be a radio-quiet millisecond pulsar, we show that such objects are a small subset of the gamma-ray pulsar population.

Traditional pulsar polarization sweep analysis starts from the point dipole rotating vector model (RVM) approximation. If augmented by a measurement of the sweep phase shift, one obtains an estimate of the emission altitude (Blaskiewicz, Cordes, & Wasserman). However, a more realistic treatment of field line sweepback and finite altitude effects shows that this estimate breaks down at modest altitude ~ 0.1R_{LC}. Such radio emission altitudes turn out to be relevant to the young energetic and millisecond pulsars that dominate the \gamma-ray population. We quantify the breakdown height as a function of viewing geometry and provide simple fitting formulae that allow observers to correct RVM-based height estimates, preserving reasonable accuracy to R ~ 0.3R_{LC}. We discuss briefly other observables that can check and improve height estimates.

We present a summary of the Fermi Pulsar Search Consortium (PSC), an international collaboration of radio astronomers and members of the Large Area Telescope (LAT) collaboration, whose goal is to organize radio follow-up observations of Fermi pulsars and pulsar candidates among the LAT gamma-ray source population. The PSC includes pulsar observers with expertise using the world's largest radio telescopes that together cover the full sky. We have performed very deep observations of all 35 pulsars discovered in blind frequency searches of the LAT data, resulting in the discovery of radio pulsations from four of them. We have also searched over 300 LAT gamma-ray sources that do not have strong associations with known gamma-ray emitting source classes and have pulsar-like spectra and variability characteristics. These searches have led to the discovery of a total of 43 new radio millisecond pulsars (MSPs) and four normal pulsars. These discoveries greatly increase the known population of MSPs in the Galactic disk, more than double the known population of so-called `black widow' pulsars, and contain many promising candidates for inclusion in pulsar timing arrays.

The parsec-scale radio properties of blazars detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been investigated using observations with the Very Long Baseline Array (VLBA). Comparisons between LAT and non-LAT detected samples were made using contemporaneous data. In total, 232 sources were used in the LAT-detected sample. This very large, radio flux-limited sample of active galactic nuclei (AGN) provides insights into the mechanism that produces strong gamma-ray emission. It has been found that LAT-detected BL Lac objects are very similar to the non-LAT BL Lac objects in most properties, although LAT BL Lac objects may have longer jets. The LAT flat spectrum radio quasars (FSRQs) are significantly different from non-LAT FSRQs and are likely extreme members of the FSRQ population. Contemporaneous observations showed a strong correlation, whereas no correlation is found using archival radio data. Most of the differences between the LAT and non-LAT populations are related to the cores of the sources, indicating that the gamma-ray emission may originate near the base of the jets (i.e., within a few pc of the central engine). There is some indication that LAT-detected sources may have larger jet opening angles than the non-LAT sources. Strong core polarization is significantly more common among the LAT sources, suggesting that gamma-ray emission is related to strong, uniform magnetic fields at the base of the jets of the blazars. Observations of sources in two epochs indicate that core fractional polarization was higher when the objects were detected by the LAT. Included in our sample are several non-blazar AGN such as 3C84, M82, and NGC 6251.

Whether a correlation exists between the radio and gamma-ray flux densities of blazars is a long-standing question, and one that is difficult to answer confidently because of various observational biases which may either dilute or apparently enhance any intrinsic correlation between radio and gamma-ray luminosities. We introduce a novel method of data randomization to evaluate quantitatively the effect of these biases and to assess the intrinsic significance of an apparent correlation between radio and gamma-ray flux densities of blazars. The novelty of the method lies in a combination of data randomization in luminosity space (to ensure that the randomized data are intrinsically, and not just apparently, uncorrelated) and significance assessment in flux space (to explicitly avoid Malmquist bias and automatically account for the limited dynamical range in both frequencies). The method is applicable even to small samples that are not selected with strict statistical criteria. For larger samples we describe a variation of the method in which the sample is split in redshift bins, and the randomization is applied in each bin individually; this variation is designed to yield the equivalent to luminosity-function sampling of the underlying population in the limit of very large, statistically complete samples. We show that for a smaller number of redshift bins, the method yields a worse significance, and in this way it is conservative in that it does not assign a stronger, artificially enhanced significance. We demonstrate how our test performs as a function of number of sources, strength of correlation, and number of redshift bins used, and we show that while our test is robust against common-distance biases and associated false positives for uncorrelated data, it retains the power of other methods in rejecting the null hypothesis of no correlation for correlated data.

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the Crab nebula lasting for approximately nine days in April of 2011. The source, which at optical wavelengths has a size of ~11 ly across, doubled its gamma-ray flux within eight hours. The peak photon flux was (186 +- 6) 10-7 cm-2 s-1 above 100 MeV, which corresponds to a 30-fold increase compared to the average value. During the flare, a new component emerged in the spectral energy distribution, which peaked at an energy of (375 +- 26) MeV at flare maximum. The observations imply that the emission region was likely relativistically beamed toward us and that variations in its motion are responsible for the observed spectral variability.

We report on optical spectroscopy of 165 Flat Spectrum Radio Quasars (FSRQs) in the Fermi 1LAC sample, which have helped allow a nearly complete study of this population. Fermi FSRQ show significant evidence for non-thermal emission even in the optical; the degree depends on the gamma-ray hardness. They also have smaller virial estimates of hole mass than the optical quasar sample. This appears to be largely due to a preferred (axial) view of the gamma-ray FSRQ and non-isotropic (H/R ~ 0.4) distribution of broad-line velocities. Even after correction for this bias, the Fermi FSRQ show higher mean Eddington ratios than the optical population. A comparison of optical spectral properties with Owens Valley Radio Observatory radio flare activity shows no strong correlation.

The radio properties of blazars detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been observed contemporaneously by the Very Long Baseline Array (VLBA). In total, 232 sources were observed with the VLBA. Ninety sources that were previously observed as part of the VLBA Imaging and Polarimetry Survey (VIPS) have been included in the sample, as well as 142 sources not found in VIPS. This very large, flux-limited sample of active galactic nuclei (AGN) provides insights into the mechanism that produces strong gamma-ray emission. In particular, we see that gamma-ray emission is related to strong, uniform magnetic fields in the cores of the host AGN. Included in this sample are non-blazar AGN such as 3C84, M82, and NGC 6251. For the blazars, the total VLBA radio flux density at 5 GHz correlates strongly with gamma-ray flux. The LAT BL Lac objects tend to be similar to the non-LAT BL Lac objects, but the LAT flat-spectrum radio quasars (FSRQs) are significantly different from the non-LAT FSRQs. Strong core polarization is significantly more common among the LAT sources, and core fractional polarization appears to increase during LAT detection.

In a radio search with the Green Bank Telescope of three unidentified low Galactic latitude Fermi-LAT sources, we have discovered the middle-aged pulsar J2030+3641, associated with 1FGL J2030.0+3641 (2FGL J2030.0+3640). Following the detection of gamma-ray pulsations using a radio ephemeris, we have obtained a phase-coherent timing solution based on gamma-ray and radio pulse arrival times that spans the entire Fermi mission. With a rotation period of 0.2 s, spin-down luminosity of 3e34 erg/s, and characteristic age of 0.5 Myr, PSR J2030+3641 is a middle-aged neutron star with spin parameters similar to those of the exceedingly gamma-ray-bright and radio-undetected Geminga. Its gamma-ray flux is 1% that of Geminga, primarily because of its much larger distance, as suggested by the large integrated column density of free electrons, DM=246 pc/cc. We fit the gamma-ray light curve, along with limited radio polarimetric constraints, to four geometrical models of magnetospheric emission, and while none of the fits have high significance some are encouraging and suggest that further refinements of these models may be worthwhile. We argue that not many more non-millisecond radio pulsars may be detected along the Galactic plane that are responsible for LAT sources, but that modified methods to search for gamma-ray pulsations should be productive -- PSR J2030+3641 would have been found blindly in gamma rays if only >0.8 GeV photons had been considered, owing to its relatively flat spectrum and location in a region of high soft background.

We have measured dramatic flux and spectral variations through the 0.193d orbit of the optical counterpart of the unidentified gamma-ray source 0FGL J2339.8-0530. This compact object companion is strongly heated, with T_eff varying from ~6900K (superior conjunction) to <3000K at minimum. A combined fit to the light curve and radial velocity amplitudes implies M_1 ~ 0.075M_sun, M_2 ~ 1.4 M_sun and inclination i ~ 57deg. Thus this is a likely `black widow' system with Edot~10^{34-34.5} erg/s pulsar driving companion mass loss. This wind, also suggested by the X-ray light curve, may prevent radio pulse detection. Our measurements constrain the pulsar's reflex motion, increasing the possibility of a pulse detection in the gamma-ray signal.

Fermi has provided the largest sample of {\gamma}-ray selected blazars to date. In this work we use a complete sample of FSRQs detected during the first year of operation to determine the luminosity function (LF) and its evolution with cosmic time. The number density of FSRQs grows dramatically up to redshift \sim0.5-2.0 and declines thereafter. The redshift of the peak in the density is luminosity dependent, with more luminous sources peaking at earlier times; thus the LF of {\gamma}-ray FSRQs follows a luminosity-dependent density evolution similarly to that of radio-quiet AGN. Also using data from the Swift Burst Alert Telescope we derive the average spectral energy distribution of FSRQs in the 10 keV-100 GeV band and show that there is no correlation of the peak {\gamma}-ray luminosity with {\gamma}-ray peak frequency. The coupling of the SED and LF allows us to predict that the contribution of FSRQs to the Fermi isotropic {\gamma}-ray background is 9.3(+1.6/-1.0) (\pm3% systematic uncertainty) in the 0.1-100GeV band. Finally we determine the LF of unbeamed FSRQs, finding that FSRQs have an average Lorentz factor of {\gamma} = 11.7(+3.3/-2.2), that most are seen within 5\circ of the jet axis, and that they represent only ~0.1 % of the parent population.

We report the detection of gamma-ray pulsations from the high-magnetic-field rotation-powered pulsar PSR J1119-6127 using data from the Fermi Large Area Telescope. The gamma-ray light curve of PSR J1119-6127 shows a single, wide peak offset from the radio peak by 0.43 pm 0.02 in phase. Spectral analysis suggests a power law of index 1.0 pm 0.3 with an energy cut-off at 0.8 pm 0.2 GeV. The first uncertainty is statistical and the second is systematic. We discuss the emission models of PSR J1119-6127 and demonstrate that despite the object's high surface magnetic field---near that of magnetars---the field strength and structure in the gamma-ray emitting zone are apparently similar to those of typical young pulsars. Additionally, we present upper limits on the \gam-ray pulsed emission for the magnetically active PSR J1846-0258 in the supernova remnant Kesteven 75 and two other energetic high-B pulsars, PSRs J1718-3718 and J1734-3333. We explore possible explanations for the non-detection of these three objects, including peculiarities in their emission geometry.

We report the probable identification of the X-ray counterpart to the gamma-ray pulsar PSR J2021+4026 using imaging with the Chandra X-ray Observatory ACIS and timing analysis with the Fermi satellite. Given the statistical and systematic errors, the positions determined by both satellites are coincident. The X-ray source position is R.A. 20h21m30.733s, Decl. +40 deg 26 min 46.04sec (J2000) with an estimated uncertainty of 1.3 arsec combined statistical and systematic error. Moreover, both the X-ray to gamma-ray and the X-ray to optical flux ratios are sensible assuming a neutron star origin for the X-ray flux. The X-ray source has no cataloged infrared-to-visible counterpart and, through new observations, we set upper limits to its optical emission of i' >23.0 mag and r' > 25.2mag. The source exhibits an X-ray spectrum with most likely both a powerlaw and a thermal component. We also report on the X-ray and visible light properties of the 43 other sources detected in our Chandra observation.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

Most pulsars observed by the Fermi LAT have gamma-ray luminosities scaling with spindown power Edot as L_gamma (Edot x 10^33 erg/s)^{1/2}. However, there exist one detection and several upper limits an order of magnitude or more fainter than this trend. We describe these `sub-luminous' gamma-ray pulsars, and discuss the case for this being an orientation effect. Of the 12 known young radio pulsars with Edot>10^34 erg/s and d<2kpc several are substantially sub-luminous. The limited available geometrical constraints favor aligned geometries for these pulsars, although no one case for alignment is compelling. In this scenario GeV emission detected from such sub-luminous pulsars can be due to a lower altitude, lower-power accelerator gap.

We announce the discovery of 1 - 100 GeV gamma-ray emission from the archetypal TeV pulsar wind nebula HESS J1825-137 using 20 months of survey data from the Fermi Large Area Telescope (LAT). The gamma-ray emission detected by the LAT is significantly spatially extended, with a best-fit rms extension of sigma = 0.56{\deg} $\pm$ 0.07{\deg} for an assumed Gaussian model. The 1 - 100 GeV LAT spectrum of this source is well described by a power-law with a spectral index of 1.38 $\pm$ 0.12 $\pm$ 0.16 and an integral flux above 1 GeV of (6.50 $\pm$ 0.21 $\pm$ 3.90) x 10^{-9} cm^{-2} s^{-1}. The first errors represent the statistical errors on the fit parameters, while the second ones are the systematic uncertainties. Detailed morphological and spectral analyses bring new constraints on the energetics and magnetic field of the pulsar wind nebula system. The spatial extent and hard spectrum of the GeV emission are consistent with the picture of an inverse Compton origin of the GeV-TeV emission in a cooling-limited nebula powered by the pulsar PSR J1826-1334.

We report the discovery of two millisecond pulsars in a search for radio pulsations at the positions of \emph{Fermi Large Area Telescope} sources with no previously known counterparts, using the Nan\c{c}ay radio telescope. The two millisecond pulsars, PSRs J2017+0603 and J2302+4442, have rotational periods of 2.896 and 5.192 ms and are both in binary systems with low-eccentricity orbits and orbital periods of 2.2 and 125.9 days respectively, suggesting long recycling processes. Gamma-ray pulsations were subsequently detected for both objects, indicating that they power the associated \emph{Fermi} sources in which they were found. The gamma-ray light curves and spectral properties are similar to those of previously-detected gamma-ray millisecond pulsars. Detailed modeling of the observed radio and gamma-ray light curves shows that the gamma-ray emission seems to originate at high altitudes in their magnetospheres. Additionally, X-ray observations revealed the presence of an X-ray source at the position of PSR J2302+4442, consistent with thermal emission from a neutron star. These discoveries along with the numerous detections of radio-loud millisecond pulsars in gamma rays suggest that many \emph{Fermi} sources with no known counterpart could be unknown millisecond pulsars.

We present precise phase-connected pulse timing solutions for 16 gamma-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one very faint radio pulsar (PSR J1124-5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multi-wavelength follow up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124-5916, and J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw component. Another Chandra ACIS image of PSR J1732-3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the gamma-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124-5916 to establish the gamma-ray to radio phase offset.

Using the Parkes radio telescope we have carried out deep observations of eleven unassociated gamma-ray sources. Periodicity searches of these data have discovered two millisecond pulsars, PSR J1103-5403 (1FGL J1103.9-5355) and PSR J2241-5236 (1FGL J2241.9-5236), and a long period pulsar, PSR J1604-44 (1FGL J1604.7-4443). In addition we searched for but did not detect any radio pulsations from six gammaray pulsars discovered by the Fermi satellite to a level of - 0.04 mJy (for pulsars with a 10% duty cycle). Timing of the millisecond pulsar PSR J1103-5403 has shown that its position is 9' from the centroid of the gamma-ray source. Since these observations were carried out, independent evidence has shown that 1FGL J1103.9-5355 is associated with the flat spectrum radio source PKS 1101-536. It appears certain that the pulsar is not associated with the gamma-ray source, despite the seemingly low probability of a chance detection of a radio millisecond pulsar. We consider that PSR J1604-44 is a chance discovery of a weak, long period pulsar and is unlikely to be associated with 1FGL J1604.7-4443. PSR J2241-5236 has a spin period of 2.2 ms and orbits a very low mass companion with a 3.5 hour orbital period. The relatively high flux density and low dispersion measure of PSR J2241-5236 makes it an excellent candidate for high precision timing experiments. The gamma-rays of 1FGL J2241.9-5236 have a spectrum that is well modelled by a power law with exponential cutoff, and phasebinning with the radio ephemeris results in a multi-peaked gamma-ray pulse profile. Observations with Chandra have identified a coincident X-ray source within 0.1" of the position of the pulsar obtained by radio timing

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

We searched for radio pulsars in 25 of the non-variable, unassociated sources in the Fermi LAT Bright Source List with the Green Bank Telescope at 820 MHz. We report the discovery of three radio and gamma-ray millisecond pulsars (MSPs) from a high Galactic latitude subset of these sources. All of the pulsars are in binary systems, which would have made them virtually impossible to detect in blind gamma-ray pulsation searches. They seem to be relatively normal, nearby (<=2 kpc) millisecond pulsars. These observations, in combination with the Fermi detection of gamma-rays from other known radio MSPs, imply that most, if not all, radio MSPs are efficient gamma-ray producers. The gamma-ray spectra of the pulsars are power-law in nature with exponential cutoffs at a few GeV, as has been found with most other pulsars. The MSPs have all been detected as X-ray point sources. Their soft X-ray luminosities of ~10^{30-31} erg/s are typical of the rare radio MSPs seen in X-rays.

We have simulated a Galactic population of young pulsars and compared with the Fermi LAT sample, constraining the birth properties, beaming and evolution of these spin-powered objects. Using quantitative tests of agreement with the distributions of observed spin and pulse properties, we find that short birth periods P_0 ~ 50ms and gamma-ray beams arising in the outer magnetosphere, dominated by a single pole, are strongly preferred. The modeled relative numbers of radio-detected and radio-quiet objects agree well with the data. Although the sample is local, extrapolation to the full Galaxy implies a gamma-ray pulsar birthrate 1/(59 yr). This is shown to be in good agreement with the estimated Galactic core collapse rate and with the local density of OB star progenitors. We give predictions for the numbers of expected young pulsar detections if Fermi LAT observations continue 10 years. In contrast to the potentially significant contribution of unresolved millisecond pulsars, we find that young pulsars should contribute little to the Galactic gamma-ray background.

The Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope provides an unprecedented opportunity to study gamma-ray blazars. To capitalize on this opportunity, beginning in late 2007, about a year before the start of LAT science operations, we began a large-scale, fast-cadence 15 GHz radio monitoring program with the 40-m telescope at the Owens Valley Radio Observatory (OVRO). This program began with the 1158 northern (declination>-20 deg) sources from the Candidate Gamma-ray Blazar Survey (CGRaBS) and now encompasses over 1500 sources, each observed twice per week with a ~4 mJy (minimum) and 3% (typical) uncertainty. Here, we describe this monitoring program and our methods, and present radio light curves from the first two years (2008 and 2009). As a first application, we combine these data with a novel measure of light curve variability amplitude, the intrinsic modulation index, through a likelihood analysis to examine the variability properties of subpopulations of our sample. We demonstrate that, with high significance (7-sigma), gamma-ray-loud blazars detected by the LAT during its first 11 months of operation vary with about a factor of two greater amplitude than do the gamma-ray-quiet blazars in our sample. We also find a significant (3-sigma) difference between variability amplitude in BL Lacertae objects and flat-spectrum radio quasars (FSRQs), with the former exhibiting larger variability amplitudes. Finally, low-redshift (z<1) FSRQs are found to vary more strongly than high-redshift FSRQs, with 3-sigma significance. These findings represent an important step toward understanding why some blazars emit gamma-rays while others, with apparently similar properties, remain silent.

The radio properties of blazars detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been observed as part of the VLBA Imaging and Polarimetry Survey (VIPS). This large, flux-limited sample of active galactic nuclei (AGN) provides insights into the mechanism that produces strong gamma-ray emission. At lower flux levels, radio flux density does not directly correlate with gamma-ray flux. We find that the LAT-detected BL Lacs tend to be similar to the non-LAT BL Lacs, but that the LAT-detected FSRQs are often significantly different from the non-LAT FSRQs. The differences between the gamma-ray loud and quiet FSRQs can be explained by Doppler boosting; these objects appear to require larger Doppler factors than those of the BL Lacs. It is possible that the gamma-ray loud FSRQs are fundamentally different from the gamma-ray quiet FSRQs. Strong polarization at the base of the jet appears to be a signature for gamma-ray loud AGN.

Pulsars PSR J0248+6021 (rotation period P=217 ms and spin-down power Edot = 2.13E35 erg/s) and PSR J2240+5832 (P=140 ms, Edot = 2.12E35 erg/s) were discovered in 1997 with the Nancay radio telescope during a northern Galactic plane survey, using the Navy-Berkeley Pulsar Processor (NBPP) filter bank. GeV gamma-ray pulsations from both were discovered using the Fermi Large Area Telescope. Twelve years of radio and polarization data allow detailed investigations. The two pulsars resemble each other both in radio and in gamma-ray data. Both are rare in having a single gamma-ray pulse offset far from the radio peak. The high dispersion measure for PSR J0248+6021 (DM = 370 pc cm^-3) is most likely due to its being within the dense, giant HII region W5 in the Perseus arm at a distance of 2 kpc, not beyond the edge of the Galaxy as obtained from models of average electron distributions. Its high transverse velocity and the low magnetic field along the line-of-sight favor this small distance. Neither gamma-ray, X-ray, nor optical data yield evidence for a pulsar wind nebula surrounding PSR J0248+6021. The gamma-ray luminosity for PSR J0248+6021 is L_ gamma = (1.4 \pm 0.3)\times 10^34 erg/s. For PSR J2240+5832, we find either L_gamma = (7.9 \pm 5.2) \times 10^34 erg/s if the pulsar is in the Outer arm, or L_gamma = (2.2 \pm 1.7) \times 10^34 erg/s for the Perseus arm. These luminosities are consistent with an L_gamma ~ sqrt(Edot) rule. Comparison of the gamma-ray pulse profiles with model predictions, including the constraints obtained from radio polarization data, favor emission in the far magnetosphere. These two pulsars differ mainly in their inclination angles and acceleration gap widths, which in turn explains the observed differences in the gamma-ray peak widths.

We have detected an Halpha bow shock nebula around PSR J1741-2054, a pulsar discovered through its GeV gamma-ray pulsations. The pulsar is only ~1.5" behind the leading edge of the shock. Optical spectroscopy shows that the nebula is non-radiative, dominated by Balmer emission. The Halpha images and spectra suggest that the pulsar wind momentum is equatorially concentrated and implies a pulsar space velocity ~150km/s, directed 15+/-10deg out of the plane of the sky. The complex Halpha profile indicates that different portions of the post-shock flow dominate line emission as gas moves along the nebula and provide an opportunity to study the structure of this unusual slow non-radiative shock under a variety of conditions. CXO ACIS observations reveal an X-ray PWN within this nebula, with a compact ~2.5" equatorial structure and a trail extending several arcmin behind. Together these data support a close (<0.5kpc) distance, a spin geometry viewed edge-on and highly efficient gamma-ray production for this unusual, energetic pulsar.

We report the discovery of eight gamma-ray pulsars in blind frequency searches using the LAT, onboard the Fermi Gamma-ray Space Telescope. Five of the eight pulsars are young (tau_c<100 kyr), energetic (Edot>10^36 erg/s), and located within the Galactic plane (|b|<3 deg). The remaining three are older, less energetic, and located off the plane. Five pulsars are associated with sources included in the LAT bright gamma-ray source list, but only one, PSR J1413-6205, is clearly associated with an EGRET source. PSR J1023-5746 has the smallest characteristic age (tau_c=4.6 kyr) and is the most energetic (Edot=1.1E37 erg/s) of all gamma-ray pulsars discovered so far in blind searches. PSRs J1957+5033 and J2055+25 have the largest characteristic ages (tau_c~1 Myr) and are the least energetic (Edot~5E33 erg/s) of the newly-discovered pulsars. We present the timing models, light curves, and detailed spectral parameters of the new pulsars. We used recent XMM observations to identify the counterpart of PSR J2055+25 as XMMU J205549.4+253959. In addition, publicly available archival Chandra X-ray data allowed us to identify the likely counterpart of PSR J1023-5746 as a faint, highly absorbed source, CXOU J102302.8-574606. The large X-ray absorption indicates that this could be among the most distant gamma-ray pulsars detected so far. PSR J1023-5746 is positionally coincident with the TeV source HESS J1023-575, located near the young stellar cluster Westerlund 2, while PSR J1954+2836 is coincident with a 4.3 sigma excess reported by Milagro at a median energy of 35 TeV. Deep radio follow-up observations of the eight pulsars resulted in no detections of pulsations and upper limits comparable to the faintest known radio pulsars, indicating that these can be included among the growing population of radio-quiet pulsars in our Galaxy being uncovered by the LAT, and currently numbering more than 20.

The \emph{Fermi} Large Area Telescope (LAT) discovered a new gamma-ray source near the Galactic plane, \object{Fermi J0109+6134}, when it flared brightly in 2010 February. The low Galactic latitude (b =-1.2\degr) indicated that the source could be located within the Galaxy, which motivated rapid multi-wavelength follow-up including radio, optical, and X-ray observations. We report the results of analyzing all 19 months of LAT data for the source, and of X-ray observations with both \emph{Swift} and the \emph{Chandra X-ray Observatory}. We determined the source redshift, z =0.783, using a Keck LRIS observation. Finally, we compiled a broadband spectral energy distribution (SED) from both historical and new observations contemporaneous with the 2010 February flare. The redshift, SED, optical line width, X-ray absorption, and multi-band variability indicate that this new GeV source is a blazar seen through the Galactic plane. Because several of the optical emission lines have equivalent width >5\AA, this blazar belongs in the flat-spectrum radio quasar category.

We demonstrate a method for quantitatively comparing gamma-ray pulsar light curves with magnetosphere beaming models. With the Fermi LAT providing many pulsar discoveries and high quality pulsar light curves for the brighter objects, such comparison allows greatly improved constraints on the emission zone geometry and the magnetospheric physics. Here we apply the method to Fermi LAT light curves of a set of bright pulsars known since EGRET or before. We test three approximate models for the magnetosphere structure and two popular schemes for the location of the emission zone, the Two Pole Caustic (TPC) model and the Outer Gap (OG) model. We find that OG models and relatively physical B fields approximating force-free dipole magnetospheres are preferred at high statistical significance. An application to the full LAT pulsar sample will allow us to follow the emission zone's evolution with pulsar spindown.

We have obtained high dynamic range, good natural seeing i' images of BL Lacertae objects (BL Lacs) to search for the AGN host and thus constrain the source redshift. These objects are drawn from a sample of bright flat-spectrum radio sources that are either known (via recent Fermi LAT observations) gamma-ray emitters or similar sources that might be detected in continuing gamma-ray observations. All had spectroscopic confirmation as BL Lac sources, but no redshift solution. We detected hosts for 25/49 objects. As these galaxies have been argued to be standard candles, our measured host magnitudes provide redshift estimates (ranging from 0.2--1.0). Lower bounds are established on the redshifts of non-detections. The mean of the fit redshifts (and lower limits) is higher than those of spectroscopic solutions in the radio- and gamma-ray- loud parent samples, suggesting corrections may be needed for the luminosity function and evolution of these sources.

We present multifrequency Very Long Baseline Array (VLBA) observations of J16021+3326. These observations, along with variability data obtained from the Owens Valley Radio Observatory (OVRO) candidate gamma-ray blazar monitoring program, clearly indicate this source is a blazar. The peculiar characteristic of this blazar, which daunted previous classification attempts, is that we appear to be observing down a precessing jet, the mean orientation of which is aligned with us almost exactly.

The vicinity of the unidentified EGRET source 3EG J1420-6038 has undergone extensive study in the search for counterparts, revealing the energetic young pulsar PSR J1420-6048 and its surrounding wind nebula as a likely candidate for at least part of the emission from this bright and extended gamma-ray source. We report on new Suzaku observations of PSR J1420-6048, along with analysis of archival XMM Newton data. The low background of Suzaku permits mapping of the extended X-ray nebula, indicating a tail stretching ~8' north of the pulsar. The X-ray data, along with archival radio and VHE data, hint at a pulsar birthsite to the North, and yield insights into its evolution and the properties of the ambient medium. We further explore such properties by modeling the spectral energy distribution (SED) of the extended nebula.

We present results from the first two years of our fast-cadence 15 GHz gamma-ray blazar monitoring program, part of the F-GAMMA radio monitoring project. Our sample includes the 1158 blazars north of -20 degrees declination from the Candidate Gamma-Ray Blazar Survey (CGRaBS), which encompasses a significant fraction of the extragalactic sources detected by the Fermi Gamma-ray Space Telescope. We introduce a novel likelihood analysis for computing a time series variability amplitude statistic that separates intrinsic variability from measurement noise and produces a quantitative error estimate. We use this method to characterize our radio light curves. We also present results indicating a statistically significant correlation between simultaneous average 15 GHz radio flux density and gamma-ray photon flux.

The Large Area Telescope (LAT), Fermi's main instrument, is providing a new view of the local energetic pulsar population. In addition to identifying a pulsar origin of a large fraction of the bright unidentified Galactic EGRET sources, the LAT results provide a great opportunity to study a sizable population of high-energy pulsars. Correlations of their physical properties, such as the trend of the luminosity versus the rotational energy loss rate, help identify global features of the gamma-ray pulsar population. Several lines of evidence, including the light curve and spectral features, suggest that gamma-ray emission from the brightest pulsars arises largely in the outer magnetosphere.

We present measurements of the Hubble diagram for 103 Type Ia supernovae (SNe) with redshifts 0.04 < z < 0.42, discovered during the first season (Fall 2005) of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey. These data fill in the redshift "desert" between low- and high-redshift SN Ia surveys. We combine the SDSS-II measurements with new distance estimates for published SN data from the ESSENCE survey, the Supernova Legacy Survey, the Hubble Space Telescope, and a compilation of nearby SN Ia measurements. Combining the SN Hubble diagram with measurements of Baryon Acoustic Oscillations from the SDSS Luminous Red Galaxy sample and with CMB temperature anisotropy measurements from WMAP, we estimate the cosmological parameters w and Omega_M, assuming a spatially flat cosmological model (FwCDM) with constant dark energy equation of state parameter, w. For the FwCDM model and the combined sample of 288 SNe Ia, we find w = -0.76 +- 0.07(stat) +- 0.11(syst), Omega_M = 0.306 +- 0.019(stat) +- 0.023(syst) using MLCS2k2 and w = -0.96 +- 0.06(stat) +- 0.12(syst), Omega_M = 0.265 +- 0.016(stat) +- 0.025(syst) using the SALT-II fitter. We trace the discrepancy between these results to a difference in the rest-frame UV model combined with a different luminosity correction from color variations; these differences mostly affect the distance estimates for the SNLS and HST supernovae. We present detailed discussions of systematic errors for both light-curve methods and find that they both show data-model discrepancies in rest-frame $U$-band. For the SALT-II approach, we also see strong evidence for redshift-dependence of the color-luminosity parameter (beta). Restricting the analysis to the 136 SNe Ia in the Nearby+SDSS-II samples, we find much better agreement between the two analysis methods but with larger uncertainties.

We present ugriz light curves for 146 spectroscopically confirmed or spectroscopically probable Type Ia supernovae from the 2005 season of the SDSS-II Supernova survey. The light curves have been constructed using a photometric technique that we call scene modelling, which is described in detail here; the major feature is that supernova brightnesses are extracted from a stack of images without spatial resampling or convolution of the image data. This procedure produces accurate photometry along with accurate estimates of the statistical uncertainty, and can be used to derive photometry taken with multiple telescopes. We discuss various tests of this technique that demonstrate its capabilities. We also describe the methodology used for the calibration of the photometry, and present calibrated magnitudes and fluxes for all of the spectroscopic SNe Ia from the 2005 season.

We report on HET and Palomar 5 m spectroscopy of recently identified $\gamma$-ray blazars in the {\it Fermi} LAT Bright Source List. These data provide identifications for 10 newly discovered $\gamma$-ray flat spectrum radio quasars (FSRQ) and six new BL Lacs plus improved spectroscopy for six additional BL Lacs. We substantially improve the identification completeness of the bright LAT blazars and give new redshifts and $z$ constraints, new estimates of the black hole masses and new measurements of the optical SED.

Sixteen pulsars have been discovered so far in blind searches of photons collected with the Large Area Telescope on the Fermi Gamma-ray Space Telescope. We here report the discovery of radio pulsations from two of them. PSR J1741-2054, with period P=413ms, was detected in archival Parkes telescope data and subsequently has been detected at the Green Bank Telescope (GBT). Its received flux varies greatly due to interstellar scintillation and it has a very small dispersion measure of DM=4.7pc/cc, implying a distance of ~0.4kpc and possibly the smallest luminosity of any known radio pulsar. At this distance, for isotropic emission, its gamma-ray luminosity above 0.1GeV corresponds to 25% of the spin-down luminosity of dE/dt=9.4e33erg/s. The gamma-ray profile occupies 1/3 of pulse phase and has three closely-spaced peaks with the first peak lagging the radio pulse by delta=0.29P. We have also identified a soft Swift source that is the likely X-ray counterpart. In many respects PSR J1741-2054 resembles the Geminga pulsar. The second source, PSR J2032+4127, was detected at the GBT. It has P=143ms, and its DM=115pc/cc suggests a distance of ~3.6kpc, but we consider it likely that it is located within the Cyg OB2 stellar association at half that distance. The radio emission is nearly 100% linearly polarized, and the main radio peak precedes by delta=0.15P the first of two narrow gamma-ray peaks that are separated by Delta=0.50P. Faint, diffuse X-ray emission in a Chandra image is possibly its pulsar wind nebula. PSR J2032+4127 likely accounts for the EGRET source 3EG J2033+4118, while its pulsar wind is responsible for the formerly unidentified HEGRA source TeV J2032+4130.

We describe a 4.85 GHz survey of bright, flat-spectrum radio sources conducted with the Effelsberg 100 m telescope in an attempt to improve the completeness of existing surveys, such as CRATES. We report the results of these observations and of follow-up 8.4 GHz observations with the VLA of a subset of the sample. We comment on the connection to the WMAP point source catalog and on the survey's effectiveness at supplementing the CRATES sky coverage.

Compact Symmetric Objects (CSOs) are small (less than 1 kpc) radio sources which have symmetric double lobes or jets. The dominant theory for the small size of these objects is that they are young radio sources which could grow into larger radio galaxies, but the currently small population of known CSOs makes it difficult to definitively determine whether or not this is the case. While a greater number of Gigahertz peaked sources can be identified by sifting through spectral surveys, this yields none of the dynamics of the sources, and also brings Quasars into the sample, which although interesting are peaked around 1 Gigahertz for very different reasons. We have used the 5 GHz VLBA Imaging and Polarization Survey (VIPS) to identify 103 CSO candidates morphologically, and are following up on these sources with multifrequency VLBA observations to confirm CSO identifications and to study their dynamics. The identification of candidates from within the survey will be discussed, as well as preliminary results from the follow-up observations.

We highlight recent theoretical and observational progress in several areas of neutron star astrophysics, and discuss the prospect for advances in the next decade.

We have simulated a population of young spin-powered pulsars and computed the beaming pattern and lightcurves for the three main geometrical models: polar cap emission, two-pole caustic ("slot gap") emission and outer magnetosphere emission. The light curve shapes depend sensitively on the magnetic inclination alpha and viewing angle zeta. We present the results as maps of observables such as peak multiplicity and gamma-ray peak separation in the (alpha, zeta) plane. These diagrams can be used to locate allowed regions for radio-loud and radio-quiet pulsars and to convert observed fluxes to true all-sky emission.

We describe a comprehensive pulsar monitoring campaign for the Large Area Telescope (LAT) on the {\em Fermi Gamma-ray Space Telescope} (formerly GLAST). The detection and study of pulsars in gamma rays give insights into the populations of neutron stars and supernova rates in the Galaxy, into particle acceleration mechanisms in neutron star magnetospheres, and into the "engines" driving pulsar wind nebulae. LAT's unprecedented sensitivity between 20 MeV and 300 GeV together with its 2.4 sr field-of-view makes detection of many gamma-ray pulsars likely, justifying the monitoring of over two hundred pulsars with large spin-down powers. To search for gamma-ray pulsations from most of these pulsars requires a set of phase-connected timing solutions spanning a year or more to properly align the sparse photon arrival times. We describe the choice of pulsars and the instruments involved in the campaign. Attention is paid to verifications of the LAT pulsar software, using for example giant radio pulses from the Crab and from PSR B1937+21 recorded at Nan\c{c}ay, and using X-ray data on PSR J0218+4232 from XMM-Newton. We demonstrate accuracy of the pulsar phase calculations at the microsecond level. Data Table 1 is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ .

We report the Suzaku observations of the high luminosity blazar SWIFT J0746.3+2548 (J0746) conducted in November 2005. This object, with z = 2.979, is the highest redshift source observed in the Suzaku Guaranteed Time Observer (GTO) period, is likely to show high gamma-ray flux peaking in the MeV range. As a result of the good photon statistics and high signal-to-noise ratio spectrum, the $Suzaku$ observation clearly confirms that J0746 has an extremely hard spectrum in the energy range of 0.3-24 keV, which is well represented by a single power-law with a photon index of 1.17 and Galactic absorption. The multiwavelength spectral energy distribution of J0746 shows two continuum components, and is well modeled assuming that the high-energy spectral component results from Comptonization of the broad-line region photons. In this paper we search for the bulk Compton spectral features predicted to be produced in the soft X-ray band by scattering external optical/UV photons by cold electrons in a relativistic jet. We discuss and provide constraints on the pair content resulting from the apparent absence of such features.

We present a measurement of the rate of type Ia supernovae (SNe Ia) from the first of three seasons of data from the SDSS-II Supernova Survey. For this measurement, we include 17 SNe Ia at redshift $z\le0.12$. Assuming a flat cosmology with $\Omega_m = 0.3=1-\Omega_\Lambda$, we find a volumetric SN Ia rate of $[2.93^{+0.17}_{-0.04}({\rm systematic})^{+0.90}_{-0.71}({\rm statistical})] \times 10^{-5} {\rm SNe} {\rm Mpc}^{-3} h_{70}^3 {\rm year}^{-1}$, at a volume-weighted mean redshift of 0.09. This result is consistent with previous measurements of the SN Ia rate in a similar redshift range. The systematic errors are well controlled, resulting in the most precise measurement of the SN Ia rate in this redshift range. We use a maximum likelihood method to fit SN rate models to the SDSS-II Supernova Survey data in combination with other rate measurements, thereby constraining models for the redshift-evolution of the SN Ia rate. Fitting the combined data to a simple power-law evolution of the volumetric SN Ia rate, $r_V \propto (1+z)^{\beta}$, we obtain a value of $\beta = 1.5 \pm 0.6$, i.e. the SN Ia rate is determined to be an increasing function of redshift at the $\sim 2.5 \sigma$ level. Fitting the results to a model in which the volumetric SN rate, $r_V=A\rho(t)+B\dot \rho(t)$, where $\rho(t)$ is the stellar mass density and $\dot \rho(t)$ is the star formation rate, we find $A = (2.8 \pm 1.2) \times 10^{-14} \mathrm{SNe} \mathrm{M}_{\sun}^{-1} \mathrm{year}^{-1}$, $B = (9.3^{+3.4}_{-3.1})\times 10^{-4} \mathrm{SNe} \mathrm{M}_{\sun}^{-1}$.

We describe recent Chandra ACIS observations of the Vela-like pulsar PSR J2021+3651 and its pulsar wind nebula (PWN). This `Dragonfly Nebula' displays an axisymmetric morphology, with bright inner jets, a double-ridged inner nebula, and a ~30" polar jet. The PWN is embedded in faint diffuse emission: a bow shock-like structure with standoff ~1' brackets the pulsar to the east and emission trails off westward for 3-4'. Thermal (kT=0.16 +/-0.02 keV) and power law emission are detected from the pulsar. The nebular X-rays show spectral steepening from Gamma=1.5 in the equatorial torus to Gamma=1.9 in the outer nebula, suggesting synchrotron burn-off. A fit to the `Dragonfly' structure suggests a large (86 +/-1 degree) inclination with a double equatorial torus. Vela is currently the only other PWN showing such double structure. The >12 kpc distance implied by the pulsar dispersion measure is not supported by the X-ray data; spectral, scale and efficiency arguments suggest a more modest 3-4 kpc.

This paper presents spectroscopy of supernovae discovered in the first season of the Sloan Digital Sky Survey-II Supernova Survey. This program searches for and measures multi-band light curves of supernovae in the redshift range z = 0.05 - 0.4, complementing existing surveys at lower and higher redshifts. Our goal is to better characterize the supernova population, with a particular focus on SNe Ia, improving their utility as cosmological distance indicators and as probes of dark energy. Our supernova spectroscopy program features rapid-response observations using telescopes of a range of apertures, and provides confirmation of the supernova and host-galaxy types as well as precise redshifts. We describe here the target identification and prioritization, data reduction, redshift measurement, and classification of 129 SNe Ia, 16 spectroscopically probable SNe Ia, 7 SNe Ib/c, and 11 SNe II from the first season. We also describe our efforts to measure and remove the substantial host galaxy contamination existing in the majority of our SN spectra.

We have applied the torus fitting procedure described in Ng & Romani (2004) to PWNe observations in the Chandra data archive. This study provides quantitative measurement of the PWN geometry and we characterize the uncertainties in the fits, with statistical errors coming from the fit uncertainties and systematic errors estimated by varying the assumed fitting model. The symmetry axis $\Psi$ of the PWN are generally well determined, and highly model-independent. We often derive a robust value for the spin inclination $\zeta$. We briefly discuss the utility of these results in comparison with new radio and high energy pulse measurements

The Sloan Digital Sky Survey-II Supernova Survey has identified a large number of new transient sources in a 300 sq. deg. region along the celestial equator during its first two seasons of a three-season campaign. Multi-band (ugriz) light curves were measured for most of the sources, which include solar system objects, Galactic variable stars, active galactic nuclei, supernovae (SNe), and other astronomical transients. The imaging survey is augmented by an extensive spectroscopic follow-up program to identify SNe, measure their redshifts, and study the physical conditions of the explosions and their environment through spectroscopic diagnostics. During the survey, light curves are rapidly evaluated to provide an initial photometric type of the SNe, and a selected sample of sources are targeted for spectroscopic observations. In the first two seasons, 476 sources were selected for spectroscopic observations, of which 403 were identified as SNe. For the Type Ia SNe, the main driver for the Survey, our photometric typing and targeting efficiency is 90%. Only 6% of the photometric SN Ia candidates were spectroscopically classified as non-SN Ia instead, and the remaining 4% resulted in low signal-to-noise, unclassified spectra. This paper describes the search algorithm and the software, and the real-time processing of the SDSS imaging data. We also present the details of the supernova candidate selection procedures and strategies for follow-up spectroscopic and imaging observations of the discovered sources.

We describe a uniform all-sky survey of bright blazars, selected primarily by their flat radio spectra, that is designed to provide a large catalog of likely gamma-ray AGN. The defined sample has 1625 targets with radio and X-ray properties similar to those of the EGRET blazars, spread uniformly across the |b| > 10 deg sky. We also report progress toward optical characterization of the sample; of objects with known R < 23, 85% have been classified and 81% have measured redshifts. One goal of this program is to focus attention on the most interesting (e.g., high redshift, high luminosity, ...) sources for intensive multiwavelength study during the observations by the Large Area Telescope (LAT) on GLAST.

We examine the radio properties of EGRET-detected blazars observed as part of the VLBA Imaging and Polarimetry Survey (VIPS). VIPS has a flux limit roughly an order of magnitude below the MOJAVE survey and most other samples that have been used to study the properties of EGRET blazars. At lower flux levels, radio flux density does not directly correlate with gamma-ray flux density. We do find that the EGRET-detected blazars tend to have higher brightness temperatures, greater core fractions, and possibly larger than average jet opening angles. A weak correlation is also found with jet length and with polarization. All of the well-established trends can be explained by systematically larger Doppler factors in the gamma-ray loud blazars, consistent with the measurements of higher apparent velocities found in monitoring programs carried out at radio frequencies above 10 GHz.

The Sloan Digital Sky Survey-II (SDSS-II) has embarked on a multi-year project to identify and measure light curves for intermediate-redshift (0.05 < z < 0.35) Type Ia supernovae (SNe Ia) using repeated five-band (ugriz) imaging over an area of 300 sq. deg. The survey region is a stripe 2.5 degrees wide centered on the celestial equator in the Southern Galactic Cap that has been imaged numerous times in earlier years, enabling construction of a deep reference image for discovery of new objects. Supernova imaging observations are being acquired between 1 September and 30 November of 2005-7. During the first two seasons, each region was imaged on average every five nights. Spectroscopic follow-up observations to determine supernova type and redshift are carried out on a large number of telescopes. In its first two three-month seasons, the survey has discovered and measured light curves for 327 spectroscopically confirmed SNe Ia, 30 probable SNe Ia, 14 confirmed SNe Ib/c, 32 confirmed SNe II, plus a large number of photometrically identified SNe Ia, 94 of which have host-galaxy spectra taken so far. This paper provides an overview of the project and briefly describes the observations completed during the first two seasons of operation.

We present extensive ugrizYHJK photometry and optical spectroscopy of SN 2005gj obtained by the SDSS-II and CSP Supernova Projects, which give excellent coverage during the first 150 days after the time of explosion. These data show that SN 2005gj is the second clear case, after SN 2002ic, of a thermonuclear explosion in a dense circumstellar environment. Both the presence of singly and doubly ionized iron-peak elements (FeIII and weak SII, SiII) near maximum light as well as the spectral evolution show that SN 2002ic-like events are Type Ia explosions. Independent evidence comes from the exponential decay in luminosity of SN 2005gj, pointing to an exponential density distribution of the ejecta. The interaction of the supernova ejecta with the dense circumstellar medium is stronger than in SN 2002ic: (1) the supernova lines are weaker; (2) the Balmer emission lines are more luminous; and (3) the bolometric luminosity is higher close to maximum light. The velocity evolution of the Halpha components suggest that the CSM around SN 2005gj is clumpy and it has a flatter density distribution compared with the steady wind solution, in agreement with SN 2002ic. An early X-ray observation with Chandra gives an upper-limit on the mass loss rate from the companion of < 2x10^{-4} Msun/yr.

We have assembled an 8.4 GHz survey of bright, flat-spectrum (alpha > -0.5) radio sources with nearly uniform extragalactic (|b|>10 deg) coverage for sources brighter than S_{4.8 GHz} = 65 mJy. The catalog is assembled from existing observations (especially CLASS and the Wright et al. PMN-CA survey), augmented by reprocessing of archival VLA and ATCA data and by new observations to fill in coverage gaps. We refer to this program as CRATES, the Combined Radio All-sky Targeted Eight GHz Survey. The resulting catalog provides precise positions, sub-arcsecond structures, and spectral indices for some 11,000 sources. We describe the morphology and spectral index distribution of the sample and comment on the survey's power to select several classes of interesting sources, especially high energy blazars. Comparison of CRATES with other high-frequency surveys also provides unique opportunities for identification of high-power radio sources.

Evidence from pulsar wind nebula symmetry axes and radio polarization observations suggests that pulsar motions correlate with the spin directions. We assemble this evidence for young isolated pulsars and show how it can be used to quantitatively constrain birth kick scenarios. We illustrate by computing several plausible, but idealized, models where the momentum thrust is proportional to the neutrino cooling luminosity of the proto-neutron star. Our kick simulations include the effects of pulsar acceleration and spin-up and our maximum likelihood comparison with the data constrains the model parameters. The fit to the pulsar spin and velocity measurements suggests that: i) the anisotropic momentum required amounts to ~10% of the neutrino flux, ii) while a pre-kick spin of the star is required, the preferred magnitude is small 10-20rad/s, so that for the best-fit models iii) the bulk of the spin is kick-induced with $\bar \Omega$ ~120rad/s and iv) the models suggest that the anisotropy emerges on a timescale $\tau$ ~1-3s.

We present the first results of the VLBA Imaging and Polarimetry Survey (VIPS), a 5 GHz VLBI survey of 1,127 sources with flat radio spectra. Through automated data reduction and imaging routines, we have produced publicly available I, Q, and U images and have detected polarized flux density from 37% of the sources. We have also developed an algorithm to use each source's I image to automatically classify it as a point-like source, a core-jet, a compact symmetric object (CSO) candidate, or a complex source. The mean ratio of the polarized to total 5 GHz flux density for VIPS sources with detected polarized flux density ranges from 1% to 20% with a median value of about 5%. We have also found significant evidence that the directions of the jets in core-jet systems tend to be perpendicular to the electric vector position angles (EVPAs). The data is consistent with a scenario in which ~24% of the polarized core-jets have EVPAs that are anti-aligned with the directions of their jet components and which have a substantial amount of Faraday rotation. In addition to these initial results, plans for future follow-up observations are discussed.

We report on VLBA astrometry and CXO imaging of PSR J0538+2817 in the supernova remnant S147. We measure a parallax distance of $1.47^{+0.42}_{-0.27}$kpc along with a high-precision proper motion, giving a transverse velocity $V_\perp=400^{+114}_{-73}$km/s. A small extended wind nebula is detected around the pulsar; the symmetry axis of this structure suggests that the spin axis lies $12\arcdeg\pm4\arcdeg$ from the velocity vector (2-D), but the emission is too faint for robust model independent statements. The neutron star is hot, consistent with the young ~40kyr kinematic age. The pulsar progenitor is likely a runaway from a nearby cluster, with NGC 1960 (M36) a leading candidate.

We describe further observations of QSO J0906+6930, a z=5.48 blazar likely to be detected in gamma-rays. New radio and X-ray data place significant constraints on any kpc-scale extension of the VLBA-detected jet. Improved optical spectroscopy detects absorption from an intervening galaxy at z=1.849 and raise the possibility that this distant, bright source is lensed. We combine the new data into an improved SED for the blazar core and comment on the Compton keV-GeV flux component.

We report on the discovery of a supermassive binary black hole system in the radio galaxy 0402+379, with a projected separation between the two black holes of just 7.3 pc. This is the closest black hole pair yet found by more than two orders of magnitude. These results are based upon recent multi-frequency observations using the Very Long Baseline Array (VLBA) which reveal two compact, variable, flat-spectrum, active nuclei within the elliptical host galaxy of 0402+379. Multi-epoch observations from the VLBA also provide constraints on the total mass and dynamics of the system. Low spectral resolution spectroscopy using the Hobby-Eberly Telescope indicates two velocity systems with a combined mass of the two black holes of ~1.5 x 10^8 solar masses. The two nuclei appear stationary while the jets emanating from the weaker of the two nuclei appear to move out and terminate in bright hot spots. The discovery of this system has implications for the number of close binary black holes that might be sources of gravitational radiation. Green Bank Telescope observations at 22 GHz to search for water masers in this interesting system are also presented.

It has been suggested that the Crab pulsar's proper motion is well aligned with the symmetry axis of the pulsar wind nebula. We have re-visited this question, examining over 6 years of F547M WFPC2 chip 3 images to obtain a best-fit value of $\mu_\ast = 14.9 \pm 0.8$mas/yr at PA $278\arcdeg \pm 3\arcdeg$. At $26\arcdeg\pm 3\arcdeg$ to the nebula axis, this substantially relaxes constraints on the birth kick of this pulsar. Such misalignment allows the momentum to be imparted over $\sim$1s timescales.

We present a time-lapse sequence of 20 near-infrared (J- and K'-band) snapshots of the central 20x20 arcsec of the Crab pulsar wind nebula, taken at subarcsecond resolution with the Hokupa'a/QUIRC adaptive optics camera on the Gemini North Telescope, and sampled at intervals of 10 min and 24 hr. It is observed that the equatorial wisps and polar knots in the termination shock of the pulsar wind appear to fluctuate in brightness on kilosecond time-scales. Maximum flux variations of +/-24 +/-4 and +/-14 +/-4 per cent relative to the mean (in 1.2 ks) are measured for the wisps and knots respectively, with greatest statistical significance in J band where the nebula background is less prominent. The J and K' flux densities imply different near-infrared spectra for the nonthermal continuum emission from the wisps and outermost polar knot (sprite), giving F_nu ~ nu^{-0.56+/-0.12} and F_nu ~ nu^{-0.21+/-0.13) respectively. The data are compared with existing optical and UV photometry and applied to constrain theories of the variability of the wisps (relativistic ion-cyclotron instability) and knots (relativistic fire hose instability).

We report on Chandra ACIS imaging of the pulsar wind nebula (PWN) of the young Vela-like PSR B1706-44, which shows the now common pattern of an equatorial wind and polar jets. The structure is particularly rich, showing a relativistically boosted termination shock, jets with strong confinement, a surrounding radio/X-ray PWN and evidence for a quasi-static `bubble nebula'. The structures trace the pulsar spin geometry and illuminate its possible relation to SNR G343.1-2.3. We also obtain improved estimates of the pulsar flux and nebular spectrum, constraining the system age and energetics.

In preparation for the Supernova Survey of the Sloan Digital Sky Survey (SDSS) II, a proposed 3-year extension to the SDSS, we have conducted an early engineering and science run during the fall of 2004, which consisted of approximately 20 scheduled nights of repeated imaging of half of the southern equatorial stripe. Transient supernova-like events were detected in near real-time and photometric measurements were made in the five SDSS filter bandpasses with a cadence of ~2 days. Candidate type Ia supernovae (SNe) were pre-selected based on their colors, light curve shape, and the properties of the host galaxy. Follow-up spectroscopic observations were performed with the Astrophysical Research Consortium 3.5m telescope and the 9.2m Hobby-Eberly Telescope to confirm their types and measure the redshifts. The 2004 campaign resulted in 22 spectroscopically confirmed SNe, which includes 16 type Ia, 5 type II, and 1 type Ib/c. These SN Ia will help fill in the sparsely sampled redshift interval of z = 0.05 - 0.35, the so-called 'redshift desert', in the Hubble diagram. Detailed investigation of the spectral properties of these moderate-redshift SNe Ia will also provide a bridge between local SNe and high-redshift objects, and will help us understand the systematics for future cosmological applications that require high photometric precision. Finally, the large survey volume also provides the opportunity to select unusual supernovae for spectroscopic study that are poorly sampled in other surveys. We report on some of the early results from this program and discuss potential future applications.

We report on Chandra ACIS and XMM-Newton MOS/PN imaging observations of two pulsar wind nebulae (K3/PSR J1420-6048 and G313.3+0.1=`the rabbit') associated with the Galactic unidentified gamma-ray source GeV J1417-6100. With the excellent ACIS imaging, the very energetic pulsar PSR J1420-6048 is separated from its surrounding nebula. This nebula has surprisingly little compact structure, although a faint arc is seen near the pulsar. Similarly, two point sources are resolved in the rabbit nebula. The large XMM-Newton collecting area provides useful spectral constraints on the rabbit and the associated point sources. Based on spectra and X-ray morphology, we identify one point source as a plausible pulsar counterpart. Large backgrounds and low source counts limited pulse search sensitivities, but we report pulse upper limits and a candidate 108ms period for the rabbit pulsar based on the XMM-Newton data and an ACIS CC observation. Comparison of the X-ray images with high resolution ATCA radio maps shows that the non-thermal X-ray emission corresponds well with the radio structure.

PSR J0537-6910 is a young, energetic, rotation-powered X-ray pulsar with a spin period of 16 ms located in the Large Magellanic Cloud. We have searched for previously undetected radio pulsations (both giant and standard) from this pulsar in a 12-hour observation taken at 1400 MHz with the Parkes 64-m radio telescope. The very large value of the magnetic field at the light cylinder radius suggests that this pulsar might be emitting giant radio pulses like those seen in other pulsars with similar field strengths. No radio emission of either kind was detected from the pulsar, and we have established an upper limit of ~ 25 mJy kpc^2 for the average 1400-MHz radio luminosity of PSR J0537-6910. The 5-sigma single-pulse detection threshold was ~ 750 mJy for a single 80-us sample. These limits are likely to be the best obtainable until searches with greatly improved sensitivity can be made with next-generation radio instruments.

We describe time-tagged observations of the Vela pulsar (PSR B0833-45) with the HST STIS MAMA-UV detectors. Using an optimal extraction technique, we obtain a NUV light curve and crude FUV phase-resolved spectra. The pulse is dominated by complex non-thermal emission at all phases. We obtain approximate spectral indices for the brightest pulse components, relate these to pulsations observed in other wavebands, and constrain the thermal surface emission from the neutron star. An upper limit on the Rayleigh-Jeans component at UV pulse minimum suggests a rather low mean surface temperature T_eff^\infty <4.6 x 10^5 (d_{300}/R_{14.2})^2 K. If confirmed by improved phase-resolved spectroscopy, this observation has significant impact on our understanding of neutron star interiors.

In preparation for GLAST, we have compiled a sample of blazar candidates to increase the pool of well studied AGN from which GLAST counterparts will be drawn. Sources were selected with our Figure of Merit (FoM) ranking; thus, they have radio and X-ray properties very similar to the EGRET blazars. Spectroscopic confirmation of these candidates is in progress, and more than 70% of these objects have been identified as flat spectrum radio quasars and BL Lac objects. We present ~250 new optical blazar identifications based on McDonald Observatory spectroscopy, 224 with redshifts. Of these, 167 are in our FoM-selected set. To motivate the Gamma-ray nature of these objects, we analyzed the current release of the EGRET data for possible point sources at their radio positions. We develop two distinct methods to combine multiple EGRET observations of a sky position into a single detection significance. We report a detection of the signal of the set of blazar candidates in the EGRET data at the > 3 sigma level by both techniques. We predict that the majority of these blazar candidates will be found by GLAST due to its increased sensitivity, duty cycle and resolving power.

We observed the gamma-ray pulsar Geminga with the HST STIS/MAMA detectors to measure the Geminga's UV spectrum and pulsations. The slope of the far-UV (FUV) spectrum is close to that of a Rayleigh-Jeans spectrum, suggesting that the FUV radiation is dominated by thermal emission from the neutron star (NS) surface. The measured FUV flux, F_FUV=(3.7+/-0.2)x10^-15 ergs cm^-2 s^-1 in 1155-1702 A band, corresponds to a brightness temperature T_RJ=(0.3-0.4)(d_200/R_13)^2 MK. The soft thermal component of the Geminga's X-ray spectrum measured with the XMM-Newton observatory corresponds to a temperature T_s=0.49+/-0.01 MK and radius R_s = (12.9+/-1.0)d_200 km. Unlike other NSs detected in the UV-optical, the FUV spectrum of Geminga lies below the extrapolation of the soft thermal component, which might be associated with the Geminga's very low temperature. Surprisingly, the thermal FUV radiation is strongly pulsed, showing a narrow dip at a phase close to that of a broader minimum of the soft X-ray light curve. The strong pulsations might be attributed to partial occultations of the thermal UV radiation by regions of the magnetosphere filled with electron/positron plasma. In contrast with the FUV spectrum, the near-IR (NIR) through near-UV (NUV) spectrum of Geminga is clearly nonthermal. It can be described by a power-law model, F_nu \propto nu^{-alpha} with a slope alpha = 0.43+/-0.15, close to alpha = 0.56+/-0.24 of the hard X-ray (E>2.5 keV) magnetospheric component. The observed NIR-optical-NUV spectrum is marginally consistent with the extrapolation of the X-ray magnetospheric spectrum. The NUV pulsations, however, do not show a clear correlation with the hard X-ray pulsations.

A remarkable number of pulsar wind nebulae (PWN) are coincident with EGRET gamma-ray sources. X-ray and radio imaging studies of unidentified EGRET sources have resulted in the discovery of at least 6 new pulsar wind nebulae (PWN). Stationary PWN (SPWN) appear to be associated with steady EGRET sources with hard spectra, typical for gamma-ray pulsars. Their toroidal morphologies can help determine the geometry of the pulsar which is useful for constraining models of pulsed gamma-ray emission. Rapidly moving PWN (RPWN) with more cometary morphologies seem to be associated with variable EGRET sources in regions where the ambient medium is dense compared to what is typical for the ISM.

PSR B0540-69 is one of a small handful of pulsars known to emit giant pulses and the only extra-galactic pulsar known to do so. We observed the pulsar for a total of 72 hr over a 6 month interval and detected 141 giant pulses. We have obtained correct phasing between the radio arrival times of the giant pulses and the X-ray pulse profile. The giant pulses occur in two phase windows, located 6.7 ms before and 5.0 ms after the midpoint of the X-ray profile. We have detected the integrated profile of the pulsar at 1.4 GHz with a flux density of only 24 micro-Jy. The statistics of the giant pulses are clearly power-law and it is likely that the giant pulses contribute only a few per cent to the integrated pulse flux density, similar to other pulsars with giant pulse emission. Simultaneous X-ray and radio observations show no significant increase in the X-ray flux at the time of the radio giants. The relative enhancements in the X-ray emission must be at least 5.5 times smaller than in the radio.

The traditional view of radio pulsar wind nebulae (PWN), encouraged by the Crab nebula's X-ray and radio morphologies, is that they are a result of the integrated history of their pulsars' wind. The radio emission should therefore be largely unaffected by recent pulsar activity, and hence minimally correlated with structures in the X-ray nebulae. Observations of several PWN, both stationary (sPWN) and rapidly moving (rPWN), now show clear morphological relationships between structures in the radio and X-ray with radio intensity variations on the order of unity. We present high-resolution X-ray and radio images of several PWN of both types and discuss the morphological relationships between the two wavebands.

We report the discovery of a radio-loud flat-spectrum QSO at z=5.47 with properties similar to those of the EGRET gamma-ray blazars. This source is the brightest radio QSO at z>5, with a pc-scale radio jet and a black hole mass estimate >10^{10}M_\odot. It appears to be the most distant blazar discovered to date. High energy observations of this source can provide powerful probes of the background radiation in the early universe.

It has been noted that the Crab and Vela pulsar proper motions lie along the symmetry axes of their wind nebulae. In an effort to promote this observation to a serious test of kick physics, we are using CXO images and other data to estimate the angle between the proper motion and PWN (i.e. spin) axis for a number of pulsars. Here we give a progress report on this work and the constraints that these data provide on kick models. Present data suggest that a kick duration of tau_K~3s is sufficient to explain the alignment of most pulsars. This rules out E-M and hydrodynamic kick models, but is fairly consistent with proposed anisotropic neutrino emission. However, some objects, especially PSR J0538+2817 show such good alignment that even neutrino models are challenged.

Supplementing existing survey data with VLA observations, we have extended $\gamma-$ray counterpart identifications down to decl. = -40$^\circ$ using our Figure of Merit approach. We find blazar counterparts for $\sim$ 70% of EGRET sources above decl. = -40$^\circ$ away from the Galaxy. Spectroscopic confirmation is in progress, and spectra for $\sim$ two dozen sources are presented here. We find evidence that increased exposure in the bulge region allowed EGRET to detect relatively faint blazars; a clear excess of non-blazar objects in this region however argues for an additional (new) source class.

We present results from X-ray and radio observations of the recently discovered young Vela-like pulsar PSR J2021+3651, which is coincident with the EGRET gamma-ray source GeV 2020+3658. A 19.0-ks Chandra ACIS-S observation has revealed a ~20'' x 10'' pulsar wind nebula that is reminiscent of the equatorial tori seen around some young pulsars, along with thermal emission from an embedded point source (kT = 0.15 +/- 0.02 keV). We name the nebula G75.2+0.1. Its spectrum is well fit by an absorbed power-law model with photon index 1.7 +/- 0.3, hydrogen column density nH = 7.8 +/- 1.7 x 10^21 cm^-2, and an unabsorbed 0.3-10.0 keV flux of 1.9 +/- 0.3 x 10^-12 erg cm^-2 s^-1. We have spatially fit G75.2+0.1 with a model that assumes a toroidal morphology, and from this we infer that the torus is highly inclined 83 deg +/- 1 deg to the line of sight. A 20.8-ks Chandra observation in continuous-clocking mode reveals a possible pulse detection, with a pulsed fraction of ~37% and an H-test probability of occuring by chance of 1.2 x 10^-4. Timing observations with the Arecibo radio telescope spanning two years show that PSR J2021+3651 glitched sometime between MJDs 52616 and 52645 with parameters delta(v)/v = (2.587 +/- 0.002) x 10^-6 and delta(dot(v))/v = (6.2 +/- 0.3) x 10^-3, similar to those of the largest glitches observed in the Vela pulsar. PSR J2021+3651 is heavily scattered (T_sc = 17.7 ms +/- 0.9 ms at 1 GHz) and exhibits a significant amount of timing noise.

We observed PSR J0437-4715 with the FUV-MAMA detector of the Hubble Space Telescope Imaging Spectrometer (STIS) to measure the pulsar's spectrum and pulsations. For the first time, UV emission from a millisecond pulsar is detected. The measured flux, $(2.0\pm 0.2)\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the 1150-1700 \AA range, corresponds to the luminosity $L_{\rm FUV}=(4.7\pm 0.5)\times 10^{27}$ erg s$^{-1}$, for the distance of 140 pc and negligible interstellar extinction. The shape of the observed spectrum suggests thermal emission from the neutron star surface with a surprisingly high temperature of about $1\times 10^5$ K, above the upper limit on the surface temperature of the younger ``ordinary'' pulsar J0108-1431. For the few-Gyr-old J0437-4715, such a temperature requires a heating mechanism to operate. The spectrum of J0437-4715 shows marginal evidence of an emission line at 1372 \AA, which might be a gravitationally redshifted Zeeman component of the Hydrogen Ly$\alpha$ line in a magnetic field $\sim 7\times 10^8$ G. No pulsations are detected, with a $3\sigma$ upper limit of 50% on pulsed fraction.

CXO imaging has shown that equatorial tori, often with polar jets, are very common in young pulsar wind nebulae (PWNe). These structures are interesting both for what they reveal about the relativistic wind itself and for the (nearly) model-independent measurement of the neutron star spin orientation they provide. The later is a particularly valuable probe of pulsar emission models and of neutron star physics.We describe here a procedure for fitting simple 3-D torus models to the X-ray data which provides robust estimates of the geometric parameters. An application to 6 PWN tori gives orientations, PWN shock scales and post-shock wind speeds along with statistical errors. We illustrate the use of these data by commenting on the implications for kick physics and for high energy beaming models.

We report the discovery of the first giant pulses from an extragalactic radio pulsar. Observations of PSR B0540-69 in the Large Magellanic Cloud made with the Parkes radio telescope at 1.38 GHz show single pulses with energy more than 5000 times that of the average pulse energy. This is only the second young pulsar, after the Crab, to show giant pulse emission. Similar to the Crab pulsar, the giant pulses occur in two distinct phase ranges and have significant arrival time jitter within these ranges. The location of the giant pulses appears to lag the peak of the (sinusoidal) X-ray profile by 0.37 and 0.64 phase, although absolute timing between the radio and X-ray data is not yet secure. The dispersion measure of the giant pulses is 146.5 pc/cc, in agreement with the detection of the pulsar at 0.64 GHz by Manchester et al. (1993). The giant pulses are scatter broadened at 1.4 GHz with an exponential scattering time of 0.4 ms and have an emission bandwidth of at least 256 MHz. In 8 hr of integration we have failed to detect any integrated flux density from the pulsar to a level of 13 microJy, assuming a duty cycle of 10%. This implies the spectral index between 0.64 and 1.38 GHz is steeper than -4.4.

Using survey data, we have re-evaluated the correlation of flat spectrum radio sources with EGRET sources in the Northern sky. A likelihood analysis incorporating the radio and X-ray properties and the Gamma-ray source localization is used to gauge the reliability of associations and to search for counterparts of previously unidentified EGRET sources. Above |b|=10deg, where the classification is complete, we find that 70% of the Northern EGRET sources have counterparts similar to the bright EGRET blazars. For several of these we identify known blazar counterparts more likely than the earlier proposed 3EG association; for ~20 we have new identifications. Spectroscopic confirmation of these candidates is in progress and we have found flat spectrum radio quasars and BL Lac counterparts with redshifts as high as 4. We also find strong evidence for a set of 28 objects with no plausible counterpart like the known EGRET Blazars. These thus represent either a new extragalactic population or a population of Galactic objects with a large scale height. The survey has been extended into the plane, where we find several new blazar candidates; the bulk of the sources are, however, Galactic. Looking ahead to the GLAST era we predict that several of the present 3EG sources are composite and that higher resolution data will break these into multiple Blazar IDs.

We find evidence for a faint wind nebula surrounding PSR J0538+2817 in CXO-ACIS imaging. This object is particularly interesting, as the pulsar spindown age is largest for any such X-ray PWN. If interpreted as an equatorial torus, the PWN supports the claimed association with the S147 supernova remnant and implies good alignment between the pulsar spin and space velocity. Comparison of the SNR age, X-ray cooling age and characteristic age suggests a birth spin period of >130ms. In turn, if we accept as causal the alignment of the linear and angular momenta, this places strong constraints on the origin of the `kick' at the neutron star birth.

We have made radio continuum, HI and X-ray observations in the direction of the unidentified EGRET source 3EG J1410-6147, using the Australia Telescope Compact Array and the Chandra X-ray Observatory. The observations encompass the supernova remnant (SNR) G312.4-0.4 and the two young pulsars PSRs J1412-6145 and J1413-6141. We derive a lower distance limit of 6 kpc to the SNR, although interpretation of positive velocity features in the HI spectrum may imply the SNR is more distant than 14 kpc. PSR J1412-6145, with an age of 50 kyr, is the pulsar most likely associated with SNR G312.4-0.4. X-rays are not detected from either pulsar and diffuse X-ray emission near the bright western edge of the SNR is weak. Although there is circumstantial evidence that this western region is a pulsar wind nebula (PWN), the embedded pulsar PSR J1412-6145 is apparently not sufficiently powerful to explain the radio enhancement. The origin of the electron acceleration in this region and of the gamma-rays remain unidentified, unless the distance to PSR J1413-6141 is at least a factor of 3 lower than its dispersion measure distance.

Recent observations have detected a number of young pulsars from the power peak in the gamma-ray band to the incoherent photon peak in the optical/IR. We have made progress on the multiwavelength phenomenology of pulsar emission and beaming, but a wide variation of light curves between different objects and different energy bands makes the full story complex. I sketch here a `Unified Model' of pulsar beaming and summarize the radiation mechanisms and their interplay in outer magnetosphere models of pulsar high energy emission.

We are developing photon-counting cameras employing cryogenic arrays of energy-resolving TES (Transition Edge Sensor) pixels. These are being tested in ground-based instruments, but will have their greatest impact when employed on space platforms, where they can cover the 10micrometer-100nm range with high time- and moderate energy- resolution. Here we summarize briefly existing device performance, current directions in array camera development and anticipated capabilities.

We have examined the soft X-ray plus optical/UV spectrum of the nearby isolated neutron star RX J1856-3754, comparing with detailed models of a thermally emitting surface. Like previous investigators, we find the spectrum is best fit by a two-temperature blackbody model. In addition, our simulations constrain the allowed viewing geometry from the observed pulse fraction upper limits. These simulations show that RX J1856-3754 is very likely to be a normal young pulsar, with the non-thermal radio beam missing Earth's line of sight. The SED limits on the model parameter space put a strong constraint on the star's M/R. At the measured parallax distance, the allowed range for Mstar=1.5Msun is Rstar=13.7+/-0.6km. Under this interpretation, the EOS is relatively stiff near nuclear density and the `Quark Star' EOS posited in some previous studies is strongly excluded. The data also constrain the surface T distribution over the polar cap.

We present a new grid of LTE model atmospheres for weakly magnetic (B<=10e10G) neutron stars, using opacity and equation of state data from the OPAL project and employing a fully frequency- and angle-dependent radiation transfer. We discuss the differences from earlier models, including a comparison with a detailed NLTE calculation. We suggest heating of the outer layers of the neutron star atmosphere as an explanation for the featureless X-ray spectra of RX J1856.5-3754 and RX J0720.4-3125 recently observed with Chandra and XMM.

We have carried out a survey for `giant pulses' in 6 young, Vela-like pulsars. In no cases did we find single pulses with flux densities more than 10 times the mean flux density. However, in PSR B1706-44 we have detected giant micro-pulses very similar to those seen in the Vela pulsar. In PSR B1706-44 these giant micro-pulses appear on the trailing edge of the profile and have an intrinsic width of ~1 ms. The cumulative probability distribution of their intensities is best described by a power-law. If the power-law continues to higher intensities, then 3.7x10^6 rotations are required to obtain a pulse with 20 times the mean pulse flux. This number is similar to the giant pulse rate in PSR B1937+21 and PSR B1821-24 but significantly higher than that for the Crab.

We report on Chandra imaging observations of the Galactic Unidentified $\gamma$-ray source GEV J1809-2327, comparing the X-ray images with new VLA 1.46 GHz and 4.86 GHz maps. The X-ray images reveal a point source connected to a non-thermal X-ray/radio nebula, supporting a pulsar/wind model for the $\gamma$-ray emitter. We also detect numerous X-ray sources from the young stellar association in the adjacent HII region S32.

There is growing evidence that some variable high-energy gamma-ray sources are associated with pulsar wind nebulae. We review the current status of X-ray and radio studies of the sources most likely to be significant emitters above 100 MeV. The assumption of the gamma-radiation arising from synchrotron processes puts a lower limit on the magnetic field and suggests there may be significant doppler boosting causing the asymmetries seen in the morphology of these nebulae. However, these sources are also in close proximity to molecular clouds, suggesting that variable gamma-ray emission may arise from a ram pressure confined pulsar wind nebula in a dense medium.

We present X-ray, radio, and infrared observations of the 68 ms pulsar PSR J1420-6048 and its surrounding nebula, a possible counterpart of the gamma-ray source GeV J1417-6100/3EG J1420-6038. Pulsed X-ray emission at the radio period is marginally detected by ASCA from a source embedded in the hard spectrum X-ray nebula AX J1420.1-6049. At radio wavelengths, the pulsar is found to be strongly linearly and circularly polarized, and the polarization sweep is measured. A comparison of high resolution ATCA radio imaging of the Kookaburra's upper wing (G313.6+0.3), which contains the pulsar and the X-ray nebula, with infrared images suggests the radio emission is partly non-thermal.

We are developing time- and energy-resolved near-IR/optical/UV photon detectors based on sharp superconducting-normal transition edges in thin films. We report observations of the Crab pulsar made during prototype testing at the McDonald 2.7m telescope with a fiber-coupled transition-edge sensor (TES) system. These data show substantial (d[alpha]~0.3), rapid variations in the spectral index through the pulse profile, with a strong phase-varying IR break across our energy band. These variations correlate with X-ray spectral variations, but no single synchrotron population can account for the full Spectral Energy Distribution (SED). We also describe test spectrophotopolarimetry observations probing the energy dependence of the polarization sweep; this may provide a new key to understanding the radiating particle population.

We have carried out a survey for `giant pulses' in 5 millisecond pulsars. We detect individual pulses from the high Edot pulsar PSR B1821-24 with energies exceeding 50x the mean pulse energy. These giant pulses are concentrated in a narrow phase window coincident with the power-law non-thermal pulse seen in hard X-rays. This is the third example of the giant pulse phenomenon. It supports the idea that large B fields in the outer magnetosphere are critical to the formation of such pulses and further suggests a direct connection between giant pulses and high energy emission.

We present X-ray and radio studies of sources which are brightabove 1 GeV (F_{>1GeV} > 4e-8 ph/cm^2/s. Only 11 out of ~30 of these gamma-ray sources have been identified with lower energy counterparts: 5 blazars and 6 pulsars. Three of these pulsars are surrounded by radio pulsar wind nebulae (PWN), two of which are also seen as bright, extended X-ray synchrotron nebulae. The ASCA X-ray telescope has observed 28 of the bright GeV sources, revealing an excess of F_{2-10keV} > 10e-12 ergs/cm^2/s sources within the {\it EGRET} error contours of the unidentified sources. Although several supernova remnants are positionally coincident with these sources, we find no X-ray evidence of high energy particle production in SNR shell shocks consistent with the GeV positions. We also present initial results from follow on radio imaging studies of several fields containing unidentified sources. We have discovered new X-ray/radio nebulae in three of these fields which are strong candidates for PWN. These sources, along with a similar nebula in CTA 1 and the PWN around PSR B1853+01 in W44, are all positionally coincident with variable EGRET sources. This suggests a class of variable gamma-ray sources associated with synchrotron emitting regions powered by the winds of young pulsars.

We present a catalog of 2-10 keV ASCA Gas Imaging Spectrometer images of fields containing bright sources of GeV emission. The images cover $\sim 85%$ of the 95% confidence position contour for 28 of the 30 sources with $> 1$ GeV fluxes above $4\times 10^{-8}$ ph ${\rm cm}^{-2} {\rm s}^{-1}$. We find an excess of hard X-ray sources with $F_{2-10 keV} \ga 10^{-12}$ ergs ${\rm cm}^{-2} {\rm s}^{-1}$ positionally coincident with unidentified sources of GeV emission. We comment on radio-loud and radio-quiet pulsar candidates, as well as several SNR and massive binaries as possible sources of the GeV emission. We also find evidence for a class of variable $\gamma-$ray sources associated with extended regions of hard X-ray emission, and propose that for these sources, a significant percentage of the $\gamma-$ray emission is generated in synchrotron nebulae surrounding fast, young pulsars.

We model the formation of magnetospheric components of millisecond pulsar light curves, deriving an approximate model for the curved space, `swept-back' dipole field and following photon emission, propagation, and scattering. Magnetospheric pulse components are strongly affected by rapid rotation and Schwarzschild effects.

We consider the effect of rapid rotation on the light curves of neutron stars with hot polar caps. For $P \approx 3$ms spin periods, the pulse fractions can be as much as an order of magnitude larger than with simple slowly-rotating (Schwarzschild) estimates. Doppler boosting, in particular, leads to characteristic distortion and ``soft lags'' in the pulse profiles, which are easily measurable in light curves with moderate energy resolution. With $\sim 10^5$ photons it should also be possible to isolate the more subtle distortions of light travel time variations and frame dragging. Detailed analysis of high quality millisecond pulsar data from upcoming X-ray missions must include these effects.

CO J=1-0 mapping observations of the dark nebula Lynds 227, which is possibly associated with the unidentified EGRET source 2EG J1811-2339, are presented. We detected a large amount of molecular gas along Lynds 227 with a total mass of (1-2)E4 Msun surrounding an X-ray synchrotron nebula, which was detected by the ASCA satellite within the error circle of 2EG J1811-2339. Molecular gas along Lynds 227 shows spatial anti-correlation with the X-ray synchrotron nebula, suggesting that the synchrotron nebula is interacting with the Lynds 227 cloud. We propose a $\gamma$-ray emission mechanism for 2EG J1811-2339: High energy electrons are injected from a rotation-powered pulsar and accelerated further in shock waves generated by the interaction with ambient matter. The high energy electrons move into the molecular cloud at the Lynds 227 and collide with dense interstellar matter to yield high energy $\gamma$-ray photons, mainly through relativistic bremsstrahlung.

We report on the first astronomical observations with a photon counting pixel detector that provides arrival time- (delta t = 100ns) and energy- (delta E_gamma < 0.15eV) resolved measurements from the near IR through the near UV. Our test observations were performed by coupling this Transition Edge Sensor (TES) device to a 0.6m telescope; we have obtained the first simultaneous optical near-IR phase-resolved spectra of the Crab pulsar. A varying infrared turnover gives evidence of self-absorption in the pulsar plasma. The potential of such detectors in imaging arrays from a space platform are briefly described.

We have mapped the radio emission in the error ellipse of GeV J1417-6100 (2EGS J1418-6049) at 13cm and 20cm using the Australia Telescope Compact Array. We find a large shell with extended wings, at the edge of which is a non-thermal, polarized structure with a center filled morphology (the `Rabbit'), coincident with an extended, hard X-ray source. We discuss the various sources seen within the ellipse as potential counterparts of the gamma-ray source. We conclude that the most likely scenario is that the Rabbit is a wind nebula surrounding a radio-quiet gamma-ray pulsar.

We have detected weak soft X-ray emission from the Pulsar Wind Nebula trailing the high velocity star PSR 2224+65 (the `Guitar Nebula'). This X-ray flux gives evidence of \gamma~10^7 eV particles in the pulsar wind and constrains the properties of the post-shock flow. The X-ray emission is most easily understood if the shocked pulsar wind is partly confined in the nebula and if magnetic fields in this zone can grow to near equipartition values.

We revisit the association of unidentified Galactic plane EGRET sources with tracers of recent massive star formation and death. Up-to-date catalogs of OB associations, SNRs, young pulsars, HII regions and young open clusters were used in finding counterparts for a recent list of EGRET sources. It has been argued for some time that EGRET source positions are correlated with SNRs and OB associations as a class; we extend such analyses by finding additional counterparts and assessing the probability of individual source identifications. Among the several scenarios relating EGRET sources to massive stars, we focus on young neutron stars as the origin of the gamma-ray emission. The characteristics of the candidate identifications are compared to the known gamma-ray pulsar sample and to detailed Galactic population syntheses using our outer gap pulsar model of gamma-ray emission. Both the spatial distribution and luminosity function of the candidates are in good agreement with the model predictions; we infer that young pulsars can account for the bulk of the excess low latitude EGRET sources. We show that with this identification, the gamma-ray point sources provide an important new window into the history of recent massive star death in the solar neighborhood.

Using a Hartree-Fock formalism, we estimate energy levels and photon cross sections for atomic iron in magnetic fields B ~ 10^13 G. Computing ionization equilibrium and normal mode opacities with these data, we construct LTE neutron star model atmospheres at 5.5 < Log(T_eff) < 6.5 and compute emergent spectra. We examine the dependence of the emergent spectra on T_eff and B. We also show the spectral variation with the angle between the magnetic field and the atmosphere normal and describe the significant limb darkening in the X-ray band. These results are compared with recent detailed computations of neutron star H model atmospheres in high fields and with low field Fe and H model atmospheres constructed from detailed opacities. The large spectral differences for different surface compositions may be discernible with present X-ray data; we also note improvements needed to allow comparison of Fe models with high quality spectra.

We compute model atmospheres and emergent spectra for low field (B<10^10 G) neutron stars, using new opacity and equation of state data from the OPAL project. These computations, incorporating improved treatments of flux transport and convective stability, provide spectra for hydrogen, solar abundance and iron atmospheres. We compare our results to high field magnetic atmospheres, available only for hydrogen. An application to apparently thermal flux from the low field millisecond pulsar PSR J0437--4715 shows that H atmospheres fit substantially better than Fe models. We comment on extension to high fields and the implication of these results for neutron star luminosities and radii.

For massive black hole binaries produced in galactic mergers, we examine the possibility of inspiral induced by interaction with field stars. We model the evolution of such binaries for a range of galaxy core and binary parameters, using numerical results from the literature to compute the binary's energy and angular momentum loss rates due to stellar encounters and including the effect of back-action on the field stars. We find that only a small fraction of binary systems can merge within a Hubble time via unassisted stellar dynamics. External perturbations may, however, cause efficient inspiral. Averaging over a population of central black holes and galaxy mergers, we compute the expected background of gravitational radiation with periods Pw ~ 1-10 yr. Comparison with sensitivities from millisecond pulsar timing suggests that the strongest sources may be detectable with modest improvements to present experiments.

We compute the variation of the beaming fraction with the efficiency of high energy gamma-ray production in the outer gap pulsar model of Romani and Yadigaroglu. This allows us to correct the fluxes observed for pulsars in the EGRET band and to derive a simple estimate of the variation of efficiency with age. Integration of this model over the population of young neutron stars gives the expected number of gamma-ray pulsars along with their distributions in age and distance. This model also shows that many of the unidentified EGRET plane sources should be pulsars, and predicts the gamma-ray fluxes of known radio pulsars. The contribution of unresolved pulsars to the background flux in the EGRET band is found to be about 5 %. For an animation of our pulsar model see http://geminga.stanford.edu/users/ion/home.html .

The computer animation illustrates the geometries described in a paper by the same authors. The preprint is available as number 9401045. The opening scene shows dipole field lines emanating from the polar caps of a rotating neutron star. The dipole axis is inclined along the green rods. The field lines shown are defined from the condition that they be tangent to the light cylinder (the cylindrical radius at which the tangential velocity of rotation reaches the speed of light). The static dipole field lines are smoothly morphed into the correct retarted-potential vacuum solutions. A red surface spanning these field lines is painted. In the next scene the blue surfaces represent the outer gaps above the surface of last closed field lines. High energy emission (blue) is produced in these outer gaps, and is beamed tangentially along the field lines. The radio emission (green) originates close to the surface of the star and is beamed along the dipole axes. The inclination angle $\alpha$ of the dipole and the viewing angle $\zeta$ are chosen to match the Crab parameters; $\alpha$ = 70, $\zeta$ = 65. The corresponding light curve is computed and shown for these angles, and the red dot traces rotation phase. The next scene shows the situation for angles appropriate to PSR1706-44; $\alpha$ = 45, $\zeta$ = 65. The final scene is a possibility for Geminga; $\alpha$ = 20, $\zeta$ = 75. These angles are poorly constrained as there is no radio emission.

There are now a half dozen young pulsars detected in high energy photons by the Compton GRO, showing a variety of emission efficiencies and pulse profiles. We present here a calculation of the pattern of high energy emission on the sky in a model which posits $\gamma$-ray production by charge depleted gaps in the outer magnetosphere. This model accounts for the radio to $\gamma$-ray pulse offsets of the known pulsars, as well as the shape of the high energy pulse profiles. We also show that $\sim 1/3$ of emitting young radio pulsars will not be detected due to beaming effects, while $\sim 2.5 \times$ the number of radio-selected $\gamma$-ray pulsars will be viewed only high energies. Finally we compute the polarization angle variation and find that the previously misunderstood optical polarization sweep of the Crab pulsar arises naturally in this picture. These results strongly support an outer-magnetosphere location for the $\gamma-$ray emission.

We describe a refined calculation of high energy emission from rotation-powered pulsars based on the Outer Gap model of Cheng, Ho \&~Ruderman (1986a,b). We have improved upon previous efforts to model the spectra from these pulsars (e. g. Cheng, et al. 1986b; Ho 1989) by following the variation in particle production and radiation properties with position in the outer gap. Curvature, synchrotron and inverse-Compton scattering fluxes vary significantly over the gap and their interactions {\it via} photon-photon pair production build up the radiating charge populations at varying rates. We have also incorporated an approximate treatment of the transport of particle and photon fluxes between gap emission zones. These effects, along with improved computations of the particle and photon distributions, provide very important modifications of the model gamma-ray flux. In particular, we attempt to make specific predictions of pulse profile shapes and spectral variations as a function of pulse phase and suggest further extensions to the model which may provide accurate computations of the observed high energy emissions.