Source author record

Paul S. Ray

Paul S. Ray 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

36works
10topics
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

36 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

preprint2023arXiv

The Radius of PSR J0740+6620 from NICER with NICER Background Estimates

We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the best available NICER background estimates to constrain the number of photons detected from the source. This approach eliminates any potential issues in the cross-calibration between these two instruments, providing thus an independent check of the robustness of the analysis. The obtained neutron star parameter constraints are compatible with the already published results, with a slight dependence on how conservative the imposed background limits are. A tighter lower limit causes the inferred radius to increase, and a tighter upper limit causes it to decrease. We also extend the study of the inferred emission geometry to examine the degree of deviation from antipodality of the hot regions. We show that there is a significant offset to an antipodal spot configuration, mainly due to the non-half-cycle azimuthal separation of the two emitting spots. The offset angle from the antipode is inferred to be above 25 degrees with 84% probability. This seems to exclude a centered-dipolar magnetic field in PSR J0740+6620.

preprint2022arXiv

A New Flaring Black Widow Candidate and Demographics of Black Widow Millisecond Pulsars in the Galactic Field

We present the discovery of a new optical/X-ray source likely associated with the Fermi $γ$-ray source 4FGL J1408.6-2917. Its high-amplitude periodic optical variability, large spectroscopic radial velocity semi-amplitude, evidence for optical emission lines and flaring, and X-ray properties together imply the source is probably a new black widow millisecond pulsar binary. We compile the properties of the 41 confirmed and suspected field black widows, finding a median secondary mass of $0.027\pm0.003\,M_{\odot}$. Considered jointly with the more massive redback millisecond pulsar binaries, we find that the "spider" companion mass distribution remains strongly bimodal, with essentially zero systems having companion masses between $\sim0.07-0.1\,M_{\odot}$. X-ray emission from black widows is typically softer and less luminous than in redbacks, consistent with less efficient particle acceleration in the intrabinary shock in black widows, excepting a few systems that appear to have more efficient "redback-like" shocks. Together black widows and redbacks dominate the census of the fastest-spinning field millisecond pulsars in binaries with known companion types, making up $\gtrsim$80% of systems with $P_{\rm{spin}}<2\,\rm{ms}$. Similar to redbacks, the neutron star masses in black widows appear on average significantly larger than the canonical $1.4\,M_{\odot}$, and many of the highest-mass neutron stars claimed to date are black widows with $M_{\rm{NS}}\gtrsim2.1\,M_{\odot}$. Both of these observations are consistent with an evolutionary picture where spider millisecond pulsars emerge from short orbital period progenitors that had a lengthy period of mass transfer initiated while the companion was on the main sequence, leading to fast spins and high masses.

preprint2022arXiv

Discovery, Timing, and Multiwavelength Observations of the Black Widow Millisecond Pulsar PSR J1555-2908

We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the radio and gamma-ray pulsation discovery and timing, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for ~10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system.

preprint2022arXiv

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.

preprint2022arXiv

The GMRT High Resolution Southern Sky Survey for pulsars and transients -- III: searching for long period pulsars

Searching for periodic non-accelerated signals in presence of ideal white noise using the fully phase-coherent Fast Folding Algorithm (FFA) is theoretically established as a more sensitive search method than the Fast Fourier Transform (FFT) search with incoherent harmonic summing. In this paper, we present a comparison of the performance of an FFA search implementation using RIPTIDE and an FFT search implementation using PRESTO, over a range of signal parameters with white noise and with real telescope noise from the GHRSS survey with the uGMRT. We find that FFA search with appropriate de-reddening of time series, performs better than FFT search with spectral whitening for long period pulsars in real GHRSS noise conditions. We describe an FFA search pipeline implemented for the GHRSS survey looking for pulsars over a period range of 0.1 s to 100 s and up to dispersion measure of 500 pc cm$^{-3}$. We processed GHRSS survey data covering $\sim$ 1500 degree$^2$ of the sky with this pipeline. We re-detected 43 known pulsars with better signal-to-noise in the FFA search than in the FFT search. We also report discovery of two new pulsars including a long period pulsar having a short duty-cycle with this FFA search pipeline. The population of long period pulsars with periods of several seconds or higher can help to constrain the pulsar death-line.

preprint2021arXiv

A Detection of Red Noise in PSR J1824$-$2452A and Projections for PSR B1937+21 using NICER X-ray Timing Data

We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long time-scale, correlated variations ("red noise") in the pulse times of arrival from the millisecond pulsars PSR J1824$-$2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency dependent propagation effects of the interstellar medium. Our Bayesian search methodology yields strong evidence (natural log Bayes factor of $9.634 \pm 0.016$) for red noise in PSR J1824$-$2452A, but is inconclusive for PSR B1937+21. In the interest of future X-ray missions, we devise and implement a method to simulate longer and higher precision X-ray datasets to determine the timing baseline necessary to detect red noise. We find that the red noise in PSR B1937+21 can be reliably detected in a 5-year mission with a time-of-arrival (TOA) error of 2 microseconds and an observing cadence of 20 observations per month compared to the 5 microsecond TOA error and 11 observations per month that NICER currently achieves in PSR B1937+21. We investigate detecting red noise in PSR B1937+21 with other combinations of observing cadences and TOA errors. We also find that an injected stochastic gravitational wave background (GWB) with an amplitude of $A_{\rm GWB}=2\times10^{-15}$ and spectral index of $γ_{\rm GWB}=13/3$ can be detected in a pulsar with similar TOA precision to PSR B1937+21, but with no additional red noise, in a 10-year mission that observes the pulsar 15 times per month and has an average TOA error of 1 microsecond.

preprint2021arXiv

An Empirical Background Model for the NICER X-ray Timing Instrument

NICER has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4-12 keV; 50 detectors) of 0.87 c/s, but in 5 percent (1 percent) of cases, the rate exceeds 10 (300) c/s. Model residuals persist at 20-30 percent of the initial rate for the brightest GTIs, implying one or more missing model parameters. Filtering criteria are given to flag GTIs likely to have unsatisfactory background predictions. With such filtering, we estimate a detection limit, 1.20 c/s (3 sigma, single GTI) at 0.4-12 keV, equivalent to 3.6e-12 erg/cm^2/s for a Crab-like spectrum. The corresponding limit for soft X-ray sources is 0.51 c/s at 0.3-2.0 keV, or 4.3e-13 erg/cm^2/s for a 100 eV blackbody. Faint-source filtering selects 85 percent of the background GTIs, and higher rates are expected for targets scheduled more favorably. An application of the model to 1 s timescale makes it possible to distinguish source flares from possible surges in the background.

preprint2021arXiv

Bayesian Solar Wind Modeling with Pulsar Timing Arrays

Using Bayesian analyses we study the solar electron density with the NANOGrav 11-year pulsar timing array (PTA) dataset. Our model of the solar wind is incorporated into a global fit starting from pulse times-of-arrival. We introduce new tools developed for this global fit, including analytic expressions for solar electron column densities and open source models for the solar wind that port into existing PTA software. We perform an ab initio recovery of various solar wind model parameters. We then demonstrate the richness of information about the solar electron density, $n_E$, that can be gleaned from PTA data, including higher order corrections to the simple $1/r^2$ model associated with a free-streaming wind (which are informative probes of coronal acceleration physics), quarterly binned measurements of $n_E$ and a continuous time-varying model for $n_E$ spanning approximately one solar cycle period. Finally, we discuss the importance of our model for chromatic noise mitigation in gravitational-wave analyses of pulsar timing data and the potential of developing synergies between sophisticated PTA solar electron density models and those developed by the solar physics community.

preprint2021arXiv

NICER Discovery of Millisecond X-ray Pulsations and an Ultracompact Orbit in IGR J17494-3030

We report the detection of 376.05 Hz (2.66 ms) coherent X-ray pulsations in NICER observations of a transient outburst of the low-mass X-ray binary IGR J17494-3030 in 2020 October/November. The system is an accreting millisecond X-ray pulsar in a 75 minute ultracompact binary. The mass donor is most likely a $\simeq 0.02 M_\odot$ finite-entropy white dwarf composed of He or C/O. The fractional rms pulsed amplitude is 7.4%, and the soft (1-3 keV) X-ray pulse profile contains a significant second harmonic. The pulsed amplitude and pulse phase lag (relative to our mean timing model) are energy-dependent, each having a local maximum at 4 keV and 1.5 keV, respectively. We also recovered the X-ray pulsations in archival 2012 XMM-Newton observations, allowing us to measure a long-term pulsar spin-down rate of $\dotν= -2.1(7)\times10^{-14}$ Hz/s and to infer a pulsar surface dipole magnetic field strength of $\simeq 10^9$ G. We show that the mass transfer in the binary is likely non-conservative, and we discuss various scenarios for mass loss from the system.

preprint2021arXiv

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature.

preprint2020arXiv

A NICER View of Spectral and Profile Evolution for Three X-ray Emitting Millisecond Pulsars

We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821-24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separation between pulse components of PSR J0218+4232 decreases with increasing energy at $\gt 3σ$ confidence. The 95% upper limit on pulse separation evolution for PSRs B1937+21 and B1821-24 is less than 2 milliperiods per keV. Our phase-resolved spectral results provide updated constraints on the non-thermal X-ray emission of these three pulsars. The photon indices of the modeled X-ray emission spectra for each pulse component of PSR B1937+21 are inconsistent with each other at the 90% confidence level, suggesting different emission origins for each pulse. We find that the PSR B1821-24 and PSR J0218+4232 emission spectra are invariant with phase at the 90% confidence level. We describe the implications of our profile and spectral results in the context of equatorial current sheet emission models for these three MSPs with non-thermal, magnetospheric X-ray emission.

preprint2020arXiv

Bright X-ray and Radio Pulses from a Recently Reactivated Magnetar

Magnetars are young, rotating neutron stars that possess larger magnetic fields ($B$ $\approx$ $10^{13}$-$10^{15}$ G) and longer rotational periods ($P$ $\approx$ 1-12 s) than ordinary pulsars. In contrast to rotation-powered pulsars, magnetar emission is thought to be fueled by the evolution and decay of their powerful magnetic fields. They display highly variable radio and X-ray emission, but the processes responsible for this behavior remain a mystery. We report the discovery of bright, persistent individual X-ray pulses from XTE J1810-197, a transient radio magnetar, using the Neutron star Interior Composition Explorer (NICER) following its recent radio reactivation. Similar behavior has only been previously observed from a magnetar during short time periods following a giant flare. However, the X-ray pulses presented here were detected outside of a flaring state. They are less energetic and display temporal structure that differs from the impulsive X-ray events previously observed from the magnetar class, such as giant flares and short X-ray bursts. Our high frequency radio observations of the magnetar, carried out simultaneously with the X-ray observations, demonstrate that the relative alignment between the X-ray and radio pulses varies on rotational timescales. No correlation was found between the amplitudes or temporal structure of the X-ray and radio pulses. The magnetar's 8.3 GHz radio pulses displayed frequency structure, which was not observed in the pulses detected simultaneously at 31.9 GHz. Many of the radio pulses were also not detected simultaneously at both frequencies, which indicates that the underlying emission mechanism producing these pulses is not broadband. We find that the radio pulses from XTE J1810-197 share similar characteristics to radio bursts detected from fast radio burst (FRB) sources, some of which are now thought to be produced by active magnetars.

preprint2020arXiv

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

preprint2019arXiv

Timing the pulsations of the accreting millisecond pulsar SAX J1808.4-3658 during its 2019 outburst

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

preprint2015arXiv

A Highly Eccentric 3.9-Millisecond Binary Pulsar in the Globular Cluster NGC 6652

We present the Robert C. Byrd Green Bank Telescope discovery of the highly eccentric binary millisecond pulsar PSR J1835$-$3259A in the Fermi Large Area Telescope-detected globular cluster NGC 6652. Timing over one orbit yields the pulse period 3.89 ms, orbital period 9.25 d, eccentricity $\sim 0.95$, and an unusually high companion mass of $0.74\,M_{\odot}$ assuming a $1.4\,M_{\odot}$ pulsar. We caution that the lack of data near periastron prevents a precise measurement of the eccentricity, and that further timing is necessary to constrain this and the other orbital parameters. From tidal considerations, we find that the companion must be a compact object. This system likely formed through an exchange encounter in the dense cluster environment. Our initial timing results predict the measurements of at least two post-Keplerian parameters with long-term phase-connected timing: the rate of periastron advance $\dotω \sim 0.1^{\circ}\,$yr$^{-1}$, requiring 1 yr of phase connection; and the Einstein delay $γ_{\mathrm{GR}} \sim 10\,$ms, requiring 2-3 yr of timing. For an orbital inclination $i > 50^{\circ}$, a measurement of $\sin{i}$ is also likely. PSR J1835$-$3259A thus provides an opportunity to measure the neutron star mass with high precision; to probe the cluster environment; and, depending on the nature of the companion, to investigate the limits of general relativity.

preprint2015arXiv

The GMRT High Resolution Southern Sky Survey for pulsars and transients -I. Survey description and initial discoveries

We are conducting a survey for pulsars and transients using the Giant Metrewave Radio Telescope (GMRT). The GMRT High Resolution Southern Sky (GHRSS) survey is an off-Galactic-plane (|b|>5) survey in the declination range -40 deg to -54 deg at 322 MHz. With the high time (up to 30.72 micro-sec) and frequency (up to 0.016275 MHz) resolution observing modes, the 5-sigma detection limit is 0.5 mJy for a 2 ms pulsar with 10% duty cycle at 322 MHz. Total GHRSS sky coverage of 2866 square-deg, will result from 1953 pointing, each covering 1.8 square-deg. The 10-sigma detection limit for a 5 milli-sec transient burst is 1.6 Jy for the GHRSS survey. In addition, the GHRSS survey can reveal transient events like rotating radio transients or fast radio bursts. With 35% of the survey completed (i.e. 1000 square-deg), we report the discovery of 10 pulsars, one of which is a millisecond pulsar (MSP), one of the highest pulsar per square degree discovery rates for any off-Galactic plane survey. We re-detected 23 known in-beam pulsars. Utilising the imaging capability of the GMRT we also localised 4 of the GHRSS pulsars (including the MSP) in the gated image plane within +/- 10 arc-second. We demonstrated rapid convergence in pulsar timing with a more precise position than is possible with single dish discoveries. We also exhibited that we can localise the brightest transient sources with simultaneously obtained lower time resolution imaging data, demonstrating a technique that may have application in the SKA.

preprint2015arXiv

X-Ray Studies of Redbacks

We consider the X-ray properties of the redback class of eclipsing millisecond pulsars. These are transitional systems between accreting low-mass X-ray binaries and binary millisecond pulsars orbiting white dwarfs, and hence their companions are non-degenerate and nearly Roche-lobe filling. The X-ray luminosity seems to scale with the fraction of the pulsar sky subtended by the companion, suggesting the shock region is not much larger than the companion, which is supported by modeling of the orbital light curves. The typical X-ray photon spectral index is $\sim 1$ and the typical 0.3-8 keV X-ray efficiency, assuming a shock size on the order of the companion's Roche lobe cross-section, is on the order of 10%. We present an overview of previous investigations, and present new observations of two redbacks, a Chandra observation of PSR J1628$-$3205 and a XMM-Newton observation of PSR J2129$-$0429. The latter shows a clearly double peaked orbital light curve with variation of the non-thermal flux by a factor of $\sim 11$, with peaks around orbital phases 0.6 and 0.9. We suggest the magnetic field of the companion plays a significant role in the X-ray emission from intrabinary shocks in redbacks.

preprint2014arXiv

Discovery of PSR J1227-4853: A transition from a low-mass X-ray binary to a redback millisecond pulsar

XSS J12270-4859 is an X-ray binary associated with the Fermi LAT gamma-ray source 1FGL J1227.9-4852. In 2012 December, this source underwent a transition where the X-ray and optical luminosity dropped and the spectral signatures of an accretion disc disappeared. We report the discovery of a 1.69 millisecond pulsar (MSP), PSR J1227-4853, at a dispersion measure of 43.4 pc cm$^{-3}$ associated with this source, using the GMRT at 607 MHz. This demonstrates that, post-transition, the system hosts an active radio MSP. This is the third system after PSR J1023+0038 and PSR J1824-2452I showing evidence of state switching between radio MSP and low-mass X-ray binary (LMXB) states. We report timing observations of PSR J1227-4853 with the GMRT and Parkes, which give a precise determination of the rotational and orbital parameters of the system. The companion mass measurement of 0.17 to 0.46 M$_{sun}$ suggests that this is a redback system. PSR J1227-4853 is eclipsed for about 40% of its orbit at 607 MHz; with additional short-duration eclipses at all orbital phases. We also find that the pulsar is very energetic, with a spin-down luminosity of ~ 10$^{35}$ erg s$^{-1}$. We report simultaneous imaging and timing observations with the GMRT, which suggests that eclipses are caused by absorption, rather than dispersion smearing or scattering.

preprint2014arXiv

Intrabinary Shock Emission from Black Widows and Redbacks

Eclipsing millisecond pulsars in close ($P_b < 1$~day) binary systems provide a different view of pulsar winds and shocks than do isolated pulsars. Since 2009, the numbers of these systems known in the Galactic field has increased enormously. We have been systematically studying many of these newly discovered systems at multiple wavelengths. Typically, the companion is nearly Roche-lobe filling and heated by the pulsar which drives mass loss from the companion. The pulsar wind shocks with this material just above the surface of the companion. We discuss various observational properties of this shock, including radio eclipses, orbitally modulated X-ray emission, and the potential for $γ$-ray emission. Redbacks, whose companions are likely non-degenerate and significantly more massive, generally have more luminous shocks than black widows which have very low mass companions. This is expected since the more massive redback companions intercept a greater fraction of the pulsar wind. We also compare these systems to accreting millisecond pulsars, which may be progenitors of black widows and in some cases can pass back and forth between redback and accretion phases.

preprint2014arXiv

Understanding the Neutron Star Population with the SKA

Since their discovery in the late 1960's the population of known neutron stars (NSs) has grown to ~2500. The last five decades of observations have yielded many surprises and demonstrated that the observational properties of NSs are remarkably diverse. The surveys that will be performed with SKA (the Square Kilometre Array) will produce a further tenfold increase in the number of Galactic NSs known. Moreover, the SKA's broad spectral coverage, sub-arraying and multi-beaming capabilities will allow us to characterise these sources with unprecedented efficiency, in turn enabling a giant leap in the understanding of their properties. Here we review the NS population and outline our strategies for studying each of the growing number of diverse classes that are populating the "NS zoo". Some of the main scientific questions that will be addressed by the much larger statistical samples and vastly improved timing efficiency provided by SKA include: (i) the spin period and spin-down rate distributions (and thus magnetic fields) at birth, and the associated information about the SNe wherein they are formed; (ii) the radio pulsar-magnetar connection; (iii) the link between normal radio pulsars, intermittent pulsars and rotating radio transients; (iv) the slowest possible spin period for a radio pulsar (revealing the conditions at the pulsar death-line); (v) proper motions of pulsars (revealing SN kick physics); (vi) the mass distribution of NSs (vii) the fastest possible spin period for a recycled pulsar (constraining magnetosphere-accretion disc interactions, gravitational wave radiation and the equation-of-state); (viii) the origin of high eccentricity millisecond pulsars (MSPs); (ix) the formation channels for recently identified triple systems; and finally (x) how isolated MSPs are formed. We expect that the SKA will break new ground unveiling exotic systems that will challenge... [abridged]

preprint2013arXiv

Looking for Stars and Finding the Moon: Effects of Lunar Gamma-ray Emission on Fermi LAT Light Curves

We are conducting a search for new gamma-ray binaries by making high signal-to-noise light curves of all cataloged Fermi LAT sources and searching for periodic variability using appropriately weighted power spectra. The light curves are created using a variant of aperture photometry where photons are weighted by the probability that they came from the source of interest. From this analysis we find that the light curves of a number of sources near the ecliptic plane are contaminated by gamma-ray emission from the Moon. This shows itself as modulation on the Moon's sidereal period in the power spectra. We demonstrate that this contamination can be removed by excluding times when the Moon was too close to a source. We advocate that this data screening should generally be used when analyzing LAT data from a source located close to the path of the Moon.

preprint2012arXiv

A wide-band, active antenna system for long wavelength radio astronomy

We describe an "active" antenna system for HF/VHF (long wavelength) radio astronomy that has been successfully deployed 256-fold as the first station (LWA1) of the planned Long Wavelength Array. The antenna system, consisting of crossed dipoles, an active balun/preamp, a support structure, and a ground screen has been shown to successfully operate over at least the band from 20 MHz (15 m wavelength) to 80 MHz (3.75 m wavelength) with a noise figure that is at least 6 dB better than the Galactic background emission noise temperature over that band. Thus, the goal to design and construct a compact, inexpensive, rugged, and easily assembled antenna system that can be deployed many-fold to form numerous large individual "stations" for the purpose of building a large, long wavelength synthesis array telescope for radio astronomical and ionospheric observations was met.

preprint2012arXiv

The Long Wavelength Array Software Library

The Long Wavelength Array Software Library (LSL) is a Python module that provides a collection of utilities to analyze and export data collected at the first station of the Long Wavelength Array, LWA1. Due to the nature of the data format and large-N ($\gtrsim$100 inputs) challenges faced by the LWA, currently available software packages are not suited to process the data. Using tools provided by LSL, observers can read in the raw LWA1 data, synthesize a filter bank, and apply incoherent de-dispersion to the data. The extensible nature of LSL also makes it an ideal tool for building data analysis pipelines and applying the methods to other low frequency arrays.

preprint2011arXiv

A 350-MHz GBT Survey of 50 Faint Fermi Gamma-ray Sources for Radio Millisecond Pulsars

We have used the Green Bank Telescope at 350MHz to search 50 faint, unidentified Fermi Gamma-ray sources for radio pulsations. So far, these searches have resulted in the discovery of 10 millisecond pulsars, which are plausible counterparts to these unidentified Fermi sources. Here we briefly describe this survey and the characteristics of the newly discovered MSPs.

preprint2011arXiv

A Radio Pulsar Search of the Gamma-ray Binaries LS I +61 303 and LS 5039

LS I +61 303 and LS 5039 are exceptionally rare examples of HMXBs with MeV-TeV emission, making them two of only five known or proposed "gamma-ray binaries". There has been disagreement within the literature over whether these systems are microquasars, with stellar winds accreting onto a compact object to produce high energy emission and relativistic jets, or whether their emission properties might be better explained by a relativistic pulsar wind colliding with the stellar wind. Here we present an attempt to detect radio pulsars in both systems with the Green Bank Telescope. The upper limits of flux density are between 4.1-14.5 uJy, and we discuss the null results of the search. Our spherically symmetric model of the wind of LS 5039 demonstrates that any pulsar emission will be strongly absorbed by the dense wind unless there is an evacuated region formed by a relativistic colliding wind shock. LS I +61 303 contains a rapidly rotating Be star whose wind is concentrated near the stellar equator. As long as the pulsar is not eclipsed by the circumstellar disk or viewed through the densest wind regions, detecting pulsed emission may be possible during part of the orbit.

preprint2011arXiv

A Search for Millisecond Pulsars in Fermi LAT-Detected Globular Clusters

We have searched for millisecond pulsars (MSPs) in two globular clusters detected by the Fermi Large Area Telescope. These clusters contained no known MSPs prior to their detections in gamma rays. The discovery of gamma ray emission from many MSPs and the prevalence of MSPs in globular clusters points to a population of MSPs as the likely source of the detected GeV emission, directing our search for new cluster MSPs. We observed NGC 6652 and NGC 6388 at 2 GHz with the Green Bank Ultimate Pulsar Processing Instrument pulsar backend at the Green Bank Telescope using the coherent dedispersion mode. We have discovered one MSP in the gamma ray error circle of NGC 6652. This pulsar is interesting because, while positionally coincident with the GC, it has a much lower dispersion measure than expected from the NE2001 galactic free electron density model. It is unclear whether the MSP is a foreground pulsar or a cluster member, and whether the pulsar, cluster, or both, is responsible for the gamma ray emission. Timing the MSP will give the pulsar position and a solid identification of the pulsar as a cluster member if it is within a few core radii of the cluster center, as well as the opportunity to search for gamma ray pulsations and determine the origin of the GeV emission.

preprint2011arXiv

Long-term X-ray monitoring of LS I +61 303: analysis of spectral variability and flares

We report on the full analysis of a Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) monitoring of the γ-ray binary system LS I +61 303. The data set covers 42 contiguous cycles of the system orbital motion. Analyzing this X-ray monitoring dataset, the largest to date for this source, we report on the variability of the orbital profile and the spectral distribution, and provide strong evidence for an anti-correlation between flux and spectral index (the higher the flux, the harder the spectral index). Furthermore, we present the analysis of two newly discovered ks-timescale flares, which present significant variability also on shorter timescales, and tend to occur at orbital phases between 0.6-0.9. However, a detailed timing analysis of the flares does not show any coherent or quasi-coherent (QPO) structure in their power spectra. We also investigated the possible appearance of the radio super-orbital modulation at X-rays energies, but we could not unambiguously detect such modulation in the system flux history, nor in the evolution of its orbital modulation fraction.

preprint2011arXiv

The Advanced X-ray Timing Array (AXTAR): A US MIDEX Mission Concept

AXTAR is a NASA MIDEX mission concept for X-ray timing of compact objects that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument is a collimated, thick Si pixel detector with 2-50 keV coverage and over 3 square meters effective area. For timing observations of accreting neutron stars and black holes, AXTAR provides at least a factor of five improvement in sensitivity over the RXTE PCA. AXTAR also carries a sensitive sky monitor that acts as a trigger for pointed observations of X-ray transients in addition to providing high duty cycle monitoring of the X-ray sky. We review the science goals and design choices that face a next generation timing mission. We then describe the technical concept for AXTAR and summarize a preliminary mission design study at the NASA/MSFC Advanced Concepts Office.

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.

preprint2010arXiv

Circularly polarized emission from the transient bursting radio source GCRT J1745-3009

We report detection of strong circularly polarized emission from the transient bursting source GCRT J1745-3009 based on new analysis of 325 MHz GMRT observations conducted on 28 September 2003. We place 8 Solar radius as the upper limit on the size of the emission region. The implied high brightness temperature required for an object beyond 1 pc and the high fraction of circular polarization firmly establish the emission as coherent. Electron cyclotron or plasma emission from a highly subsolar magnetically dominated dwarf located less than 4 kpc away could have given rise to the GCRT radio emission.

preprint2010arXiv

H-alpha Emission Variability in the gamma-ray Binary LS I +61 303

LS I +61 303 is an exceptionally rare example of a high mass X-ray binary (HMXB) that also exhibits MeV-TeV emission, making it one of only a handful of "gamma-ray binaries". Here we present H-alpha spectra that show strong variability during the 26.5 day orbital period and over decadal time scales. We detect evidence of a spiral density wave in the Be circumstellar disk over part of the orbit. The H-alpha line profile also exhibits a dramatic emission burst shortly before apastron, observed as a redshifted shoulder in the line profile, as the compact source moves almost directly away from the observer. We investigate several possible origins for this red shoulder, including an accretion disk, mass transfer stream, and a compact pulsar wind nebula that forms via a shock between the Be star's wind and the relativistic pulsar wind.

preprint2010arXiv

Pulsar Timing with the Fermi LAT

We present an overview of precise pulsar timing using data from the Large Area Telescope (LAT) on Fermi. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. In addition to determining the spindown behavior of the pulsars and detecting glitches and timing noise, such timing analyses allow the precise determination of the pulsar position, thus enabling detailed multiwavelength follow up.

preprint2009arXiv

The Orbits of the Gamma-ray Binaries LS I +61 303 and LS 5039

LS I +61 303 and LS 5039 are two of only a handful of known high mass X-ray binaries (HMXBs) that exhibit very high energy emission in the MeV-TeV range, and these "gamma-ray binaries" are of renewed interest due to the recent launch of the Fermi Gamma-ray Space Telescope. Here we present new radial velocities of both systems based on recent red and blue optical spectra. Both systems have somewhat discrepant orbital solutions available in the literature, and our new measurements result in improved orbital elements and resolve the disagreements. The improved geometry of each orbit will aid in studies of the high energy emission region near each source.

preprint2008arXiv

GCRT J1742-3001: A New Radio Transient Towards the Galactic Center

We report the detection of a new transient radio source, GCRT J1742-3001, located ~1 degree from the Galactic center. The source was detected ten times from late 2006 to 2007 May in our 235 MHz transient monitoring program with the Giant Metrewave Radio Telescope (GMRT). The radio emission brightened in about one month, reaching a peak observed flux density of ~100 mJy on 2007 January 28, and decaying to ~50 mJy by 2007 May when our last monitoring observation was made. Two additional faint, isolated 235 MHz detections were made in mid-2006, also with the GMRT. GCRT J1742-3001 is unresolved at each epoch, with typical resolutions of ~20 arcsec x 10 arcsec. No polarization information is available from the observations. Based on nondetections in observations obtained simultaneously at 610 MHz, we deduce that the spectrum of GCRT J1742-3001 is very steep, with a spectral index less than about -2. Follow-up radio observations in 2007 September at 330 MHz and 1.4 GHz, and in 2008 February at 235 MHz yielded no detections. No X-ray counterpart is detected in a serendipitous observation obtained with the X-ray telescope aboard the Swift satellite during the peak of the radio emission in early 2007. We consider the possibilities that GCRT J1742-3001 is either a new member of an existing class of radio transients, or is representative of a new class having no associated X-ray emission.

preprint2007arXiv

The long-term evolution of the spin, pulse shape, and orbit of the accretion-powered millisecond pulsar SAX J1808.4-3658

We present a 7 yr timing study of the 2.5 ms X-ray pulsar SAX J1808.4-3658, an X-ray transient with a recurrence time of ~2 yr, using data from the Rossi X-ray Timing Explorer covering 4 transient outbursts (1998-2005). We verify that the 401 Hz pulsation traces the spin frequency fundamental and not a harmonic. Substantial pulse shape variability, both stochastic and systematic, was observed during each outburst. Analysis of the systematic pulse shape changes suggests that, as an outburst dims, the X-ray "hot spot" on the pulsar surface drifts longitudinally and a second hot spot may appear. The overall pulse shape variability limits the ability to measure spin frequency evolution within a given X-ray outburst (and calls previous nudot measurements of this source into question), with typical upper limits of |nudot| < 2.5x10^{-14} Hz/s (2 sigma). However, combining data from all the outbursts shows with high (6 sigma) significance that the pulsar is undergoing long-term spin down at a rate nudot = (-5.6+/-2.0)x10^{-16} Hz/s, with most of the spin evolution occurring during X-ray quiescence. We discuss the possible contributions of magnetic propeller torques, magnetic dipole radiation, and gravitational radiation to the measured spin down, setting an upper limit of B < 1.5x10^8 G for the pulsar's surface dipole magnetic field and and Q/I < 5x10^{-9} for the fractional mass quadrupole moment. We also measured an orbital period derivative of Pdot = (3.5+/-0.2)x10^{-12} s/s. This surprising large Pdot is reminiscent of the large and quasi-cyclic orbital period variation observed in the so-called "black widow" millisecond radio pulsars, supporting speculation that SAX J1808.4-3658 may turn on as a radio pulsar during quiescence. In an appendix we derive an improved (0.15 arcsec) source position from optical data.