Source author record

Roger W. Romani

Roger W. Romani appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2026arXiv

X-ray Polarization of the Intrabinary Shock in Redback Pulsar J1723$-$2837

The intrabinary shocks (IBS) in spider pulsars emit non-thermal synchrotron X-rays from accelerated electrons and positrons in the shocked pulsar wind, likely energized by magnetic reconnection. The double-peaked X-ray light curves from these shocks have been well characterized in several spider systems. In this paper, we analyze Imaging X-ray Polarimetry Explorer (IXPE) observations of the redback pulsar J1723$-$2837 to examine the expected synchrotron polarization. Using advanced extraction methods that include spatial, temporal, and particle background weights, we constrain the polarization of the IBS. We compare different models for the magnetic field in the radiation zone and find that the best fit prefers a striped pulsar wind model over other polarized models, with maximum polarization degree of the IBS emission component $Π_{\rm IBS}=36^{+16}_{-15}\%$, in addition to an unpolarized non-IBS component. Since this is only 2.4$σ$, we cannot claim strong preference over an unpolarized model; we report a $99\%$ confidence level upper limit on the total polarization of both IBS and non-IBS components $Π_{99}<36\%$, which is improved over the $50\%$ limit obtained in previous work. The best-fit polarization of the IBS component is consistent with numerical simulations. Detailed tests of such models are accessible to future measurements.

preprint2023arXiv

Determination of X-ray pulsar geometry with IXPE polarimetry

Using observations of X-ray pulsar Her X-1 by the Imaging X-ray Polarimetry Explorer, we report on a highly significant ($>17σ$) detection of the polarization signal from an accreting neutron star. The observed degree of the linear polarization of $\sim$10\% is found to be far below theoretical expectations for this object, and stays low throughout the spin cycle of the pulsar. Both the polarization degree and the angle exhibit variability with pulse phase, which allowed us to measure the pulsar spin position angle 57(2) deg and the magnetic obliquity 12(4) deg, which is an essential step towards detailed modelling of the intrinsic emission of X-ray pulsars. Combining our results with the optical polarimetric data, we find that the spin axis of the neutron star and the angular momentum of the binary orbit are misaligned by at least $\sim$20 deg, which is a strong argument in support of the models explaining stability of the observed super-orbital variability with the precession of the neutron star.

preprint2022arXiv

An Optical Study of the Black Widow Population

The optical study of the heated substellar companions of `Black Widow' (BW) millisecond pulsars (MSP) provides unique information on the MSP particle and radiation output and on the neutron star mass. Here we present analysis of optical photometry and spectroscopy of a set of relatively bright BWs, many newly discovered in association with Fermi $γ$-ray sources. Interpreting the optical data requires sophisticated models of the companion heating. We provide a uniform analysis, selecting the preferred heating model and reporting on the companion masses and radii, the pulsar heating power and neutron star mass. The substellar companions are substantially degenerate, with average densities $15-30\times$ Solar, but are inflated above their zero temperature radii. We find evidence that the most extreme recycled BW pulsars have both large $>0.8M_\odot$ accreted mass and low $<10^8$G magnetic fields. Examining a set of heavy BWs, we infer that neutron star masses larger than $2.19 M_\odot$ ($1σ$ confidence) or $2.08 M_\odot$ ($3σ$ confidence) are required; these bounds exclude all but the stiffest equations of state in standard tabulations.

preprint2022arXiv

Limits on X-ray Polarization at the Core of Centaurus A as Observed with the Imaging X-ray Polarimetry Explorer

We present measurements of the polarization of X-rays in the $2-8 \thinspace \mathrm{keV}$ band from the nucleus of the radio galaxy Centaurus A (Cen A), using a 100ks observation from the Imaging X-ray Polarimetry Explorer (IXPE). Nearly simultaneous observations of Cen A were also taken with the Swift, NuSTAR, and INTEGRAL observatories. No statistically significant degree of polarization is detected with IXPE. These observations have a minimum detectable polarization at $99 \%$ percent confidence (MDP$_{99}$) of $6.5 \%$ using a weighted, spectral model-independent calculation in the $2-8 \thinspace \mathrm{keV}$ band. The polarization angle $ψ$ is consequently unconstrained. Spectral fitting across three orders of magnitude in X-ray energy ($0.3-400 \thinspace \mathrm{keV}$) demonstrates that the SED of Cen A is well described by a simple power law with moderate intrinsic absorption ($N_H \sim 10^{23} \thinspace \mathrm{cm}^{-2}$) and a Fe K$α$ emission line, although a second unabsorbed power law is required to account for the observed spectrum at energies below $2 \thinspace \mathrm{keV}$. This spectrum suggests that the reprocessing material responsible for this emission line is optically thin and distant from the central black hole. Our upper limits on the X-ray polarization are consistent with the predictions of Compton scattering, although the specific seed photon population responsible for production of the X-rays cannot be identified. The low polarization degree, variability in the core emission, and the relative lack of variability in the Fe K$α$ emission line support a picture where electrons are accelerated in a region of highly disordered magnetic fields surrounding the innermost jet.

preprint2022arXiv

PSR J0952-0607: The Fastest and Heaviest Known Galactic Neutron Star

We describe Keck-telescope spectrophotometry and imaging of the companion of the ``black widow" pulsar PSR~J0952$-$0607, the fastest known spinning neutron star (NS) in the disk of the Milky Way. The companion is very faint at minimum brightness, presenting observational challenges, but we have measured multicolor light curves and obtained radial velocities over the illuminated ``day" half of the orbit. The model fits indicate system inclination $i=59.8\pm 1.9^\circ$ and a pulsar mass $M_{NS} = 2.35\pm 0.17 M_\odot$, the largest well-measured mass found to date. Modeling uncertainties are small, since the heating is not extreme; the companion lies well within its Roche lobe and a simple direct-heating model provides the best fit. If the NS started at a typical pulsar birth mass, nearly $1 M_\odot$ has been accreted; this may be connected with the especially low intrinsic dipole surface field, estimated at $6\times 10^7$G. Joined with reanalysis of other black widow and redback pulsars, we find that the minimum value for the maximum NS mass is $M_{\rm max} > 2.19 M_\odot$$(2.09 M_\odot)$ at $1σ$$(3σ)$ confidence. This is $\sim 0.15 M_\odot$ heavier than the lower limit on $M_{\rm max}$ implied by the white-dwarf--pulsar binaries measured via radio Shapiro-delay techniques.

preprint2022arXiv

Simultaneous space and phase resolved X-ray polarimetry of the Crab Pulsar and Nebula

The Crab pulsar and its nebula are among the most studied astrophysical systems, and constitute one of the most promising environments where high energy processes and particle acceleration can be investigated. They are the only objects for which previous X-ray polarisation has been reported. We present here the first Imaging X-ray Polarimetry Explorer (IXPE) observation of the Crab pulsar and nebula. The total pulsar pulsed emission in the [2-8] keV energy range is unpolarised. Significant polarisation up to 15% is detected only in the core of the main peak. The nebula has a total space integrated polarised degree of 20% and polarisation angle of 145deg. The polarised maps show a large variation in the local polarisation, and regions with polarised degree up to 45-50%. The polarisation pattern suggests a predominantly toroidal magnetic field.

preprint2022arXiv

Testing High-Energy Emission Models for Blazars with X-ray Polarimetry

Both leptonic and hadronic emission processes may contribute to blazar jet emission; which dominates in blazars's high energy emission component remains an open question. Some intermediate synchrotron peaked blazars transition from their low to high energy emission components in the X-ray band making them excellent laboratories to probe both components simultaneously, and good targets for the newly launched Imaging X-ray Polarimetry Explorer. We characterize the spectral energy distributions for three such blazars: CGRaBS~J0211+1051, TXS~0506+056, and S5~0716+714, predicting their X-ray polarization behavior by fitting a multizone polarized leptonic jet model. We find that a significant detection of electron synchrotron dominated polarization is possible with a 300~ks observation for S5~0716+714 and CGRaBS~J0211+1051 in their flaring states, while even 500~ks observations are unlikely to measure synchrotron self-Compton polarization. Importantly, non-leptonic emission processes like proton synchrotron are marginally detectable for our brightest ISP, S5~0716+714, during a flaring state. Improved {\it IXPE} data reduction methods or next generation telescopes like {\it eXTP} are needed to confidently measure SSC polarization.

preprint2022arXiv

The Long Filament of PSR J2030+4415

New X-ray and optical observations shed light on the remarkable X-ray filament of the Gamma-ray pulsar PSR~J2030+4415. Images of the associated H$α$ bow shock's evolution over the past decade compared with its velocity structure provide an improved kinematic distance of $\sim$0.5kpc. These velocities also imply that the pulsar spin axis lies $\sim 15^\circ$ from the proper motion axis which is close to the plane of the sky. The multi-bubble shock structure indicates that the bow shock stand-off was compressed to a small value $\sim 20-30$y ago when the pulsar broke through the bow shock to its present bubble. This compression allowed multi-TeV pulsar $e^\pm$ to escape to the external ISM, lighting up an external magnetic field structure as the `filament'. The narrow filament indicates excellent initial confinement and the full $15^\prime$ ($2.2$~pc=7~lt-y) projected length of the filament indicates rapid $e^\pm$ propagation to its end. Spectral variation along the filament suggests that the injected particle energy evolved during the break-through event.

preprint2021arXiv

PSR J1810+1744: Companion Darkening and a Precise High Neutron Star Mass

Keck-telescope spectrophotometry of the companion of PSR J1810+1744 shows a flat, but asymmetric light-curve maximum and a deep, narrow minimum. The maximum indicates strong gravity darkening near the L_1 point, along with a heated pole and surface winds. The minimum indicates a low underlying temperature and substantial limb darkening. The gravity darkening is a consequence of extreme pulsar heating and the near-filling of the Roche lobe. Light-curve modeling gives a binary inclination i=65.7+/-0.4deg. With the Keck-measured radial-velocity amplitude K_c=462.3+/-2.2km/s, this gives an accurate neutron star mass M_NS=2.13+/-0.04M_o, with important implications for the dense-matter equation of state. A classic direct-heating model, ignoring the L_1 gravitational darkening, would predict an unphysical M_NS>3M_o. A few other ``spider" pulsar binaries have similar large heating and fill factor; thus, they should be checked for such effects.

preprint2020arXiv

Atmospheric Circulation on Black Widow Companions

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.

preprint2020arXiv

Orbital Modulation of Gamma Rays from PSR~J2339$-$0533

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.

preprint2020arXiv

PSR J1709-4429's Proper Motion and its Relationship to SNR G343.1-2.3

We have obtained a deep 670 ks CXO ACIS image of the remarkable pulsar wind nebula (PWN) of PSR J1709-4429, in 4 epochs during 2018-2019. Comparison with an archival 2004 data set provides a pulsar proper motion mu = 13 \pm 3 mas/yr at a PA of 86 \pm 9 degree (1 sigma combined statistical and systematic uncertainties), precluding birth near the center of SNR G343.1-2.3. At the pulsar's characteristic age of 17 kyr, the association can be preserved through a combination of progenitor wind, birth kick and PWN outflow. Associated TeV emission may, however, indicate explosion in an earlier supernova. Inter-epoch comparison of the X-ray images shows that the PWN is dynamic, but we are unable to conclusively measure flow speeds from blob motion. The pulsar has generated a radio/X-ray wind bubble, and we argue that the PWN's long narrow jets are swept back by shocked pulsar wind venting from this cavity. These jets may trace the polar magnetic field lines of the PWN flow, an interesting challenge for numerical modeling.

preprint2020arXiv

PSR J2030+4415's Remarkable Bow Shock, PWN and Filament

We report on new X-ray and optical observations of PSR J2030+4415, a Gamma-ray Pulsar with an H$α$ 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$α$ velocity structure this allows us to estimate the distance as $0.75$ kpc.

preprint2016arXiv

Binary Black Holes, Gas Sloshing, and Cold Fronts in the X-ray Halo Hosting 4C+37.11

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.

preprint2016arXiv

Deep Chandra Observations of the Pulsar Wind Nebula Created by PSR B0355+54

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.

preprint2016arXiv

Intra-Binary Shock Heating of Black Widow Companions

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.

preprint2015arXiv

A Spectroscopic Study of the Extreme Black Widow PSR J1311-3430

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.

preprint2015arXiv

Keck Spectroscopy of Millisecond Pulsar J2215+5135: a Moderate-M_NS, High-Inclination Binary

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.

preprint2015arXiv

X-ray analysis of the proper motion and pulsar wind nebula for PSR J1741-2054

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 $μ=109 \pm 10 {\rm mas\ yr}^{-1}$ in a direction consistent with the symmetry axis of the observed H$α$ 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 $Γ$=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 $Γ$ = 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.

preprint2014arXiv

2FGL J1653.6-0159: A New Low in Evaporating Pulsar Binary Periods

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.

preprint2014arXiv

A Survey for H$α$ Pulsar Bow Shocks

We report on a survey for H$α$ bow shock emission around nearby $γ$-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$α$ structures around two additional $γ$-ray pulsars, PSR J2030+4415 and PSR J1509$-$5850. These are the first known examples of H$α$ 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 $γ$-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 $ϕ_{WNM}\sim 0.3$ there should be a total of $\sim$ nine H$α$ bow shocks in our LAT-targeted survey; given that seven such objects are now known, a much larger $ϕ_{WNM}$ seems problematic.

preprint2014arXiv

Imaging Redshift Estimates for Fermi BL Lacs

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.

preprint2014arXiv

Temporal Correlations Between Optical and Gamma-ray Activity in Blazars

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σ$) 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.

preprint2013arXiv

A Multi-wavelength Study of the Host Environment of SMBHB 4C+37.11

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.

preprint2013arXiv

Chandra, Keck and VLA Observations of the Crab Nebula during the 2011-April Gamma-ray Flare

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.

preprint2013arXiv

Photometrically Triggered Keck Spectroscopy of Fermi BL Lac Objects

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.

preprint2013arXiv

Spectroscopy of The Largest Ever Gamma-ray Selected BL Lac Sample

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.

preprint2012arXiv

2FGL J1311.7-3429 Joins the Black Widow Club

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.

preprint2012arXiv

A Chandra Proper Motion for PSR J1809-2332

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.

preprint2012arXiv

Altitude Limits for Rotating Vector Model Fitting of Pulsar Polarization

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 γ-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.

preprint2012arXiv

Multi-Zone Modeling of The Pulsar Wind Nebula HESS J1825-137

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 $τ_{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.

preprint2012arXiv

PSR J1311-3430: A Heavyweight Neutron Star with a Flyweight Helium Companion

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.

preprint2012arXiv

Spectroscopy of Broad Line Blazars from 1LAC

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.

preprint2011arXiv

Contemporaneous VLBA 5 GHz Observations of LAT-Detected Blazars

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.

preprint2011arXiv

The Identification of the X-ray Counterpart to PSR J2021+4026

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.

preprint2011arXiv

The Orbit and Companion of Probable gamma-ray Pulsar J2339-0533

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.

preprint2010arXiv

Blazars in the Fermi Era: The OVRO 40-m Telescope Monitoring Program

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.

preprint2010arXiv

Constraining Pulsar Magnetosphere Geometry with Gamma-Ray Light Curves

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.

preprint2010arXiv

Imaging Redshift Estimates for BL Lacertae Objects

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.

preprint2010arXiv

The Balmer-dominated Bow Shock and Wind Nebula Structure of Gamma-ray Pulsar PSR J1741-2054

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.

preprint2010arXiv

The Extended X-ray Nebula of PSR J1420-6048

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.

preprint2010arXiv

The Galactic Population of Young Gamma-ray Pulsars

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.

preprint2009arXiv

Filling in the Gaps in the 4.85 GHz Sky

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.

preprint2009arXiv

Optical Spectroscopy of Bright Fermi LAT Blazars

We report on HET and Palomar 5 m spectroscopy of recently identified $γ$-ray blazars in the {\it Fermi} LAT Bright Source List. These data provide identifications for 10 newly discovered $γ$-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.

preprint2008arXiv

An Atlas For Interpreting Gamma-Ray Pulsar Light Curves

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.

preprint2007arXiv

Birth Kick Distributions and the Spin-Kick Correlation of Young Pulsars

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 Ω$ ~120rad/s and iv) the models suggest that the anisotropy emerges on a timescale $τ$ ~1-3s.

preprint2005arXiv

A Northern Survey of Gamma-Ray Blazar Candidates

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.

preprint2002arXiv

Chandra Imaging of the Gamma-Ray Source GeV J1809-2327

We report on Chandra imaging observations of the Galactic Unidentified $γ$-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 $γ$-ray emitter. We also detect numerous X-ray sources from the young stellar association in the adjacent HII region S32.

preprint2001arXiv

Phase-resolved Crab Studies with a Cryogenic TES Spectrophotometer

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.

preprint2001arXiv

X-Ray and Radio Observations of Bright GeV Sources

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.

preprint1995arXiv

Ultra-Low Frequency Gravitational Radiation from Massive Black Hole Binaries

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.