Source author record

Matthew Kerr

Matthew Kerr 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

43works
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

43 published item(s)

preprint2026arXiv

Detection of giant pulses from the transitional millisecond pulsar J1227$-$4853

We report the discovery of giant pulse (GP) emission from the transitional millisecond pulsar (tMSP) PSR J1227$-$4853, using 174 hours of single-pulse data from the upgraded Giant Metrewave Radio Telescope (uGMRT). This marks the first detection of GPs from a transitional MSP and adds to the small number of millisecond pulsars known to exhibit such extreme variability. A total of 235 GPs were detected across observations at 550-750 MHz, with widths as narrow as $1.28~μs$ and flux densities up to $\sim 10^4$ times the pulsar's mean flux density. The GPs are strongly localized in pulse phase, originating predominantly from the second and third main-pulse components, and are absent in the inter-pulse region. The arrival times of the GPs deviate significantly from Poisson statistics, with the waiting-time distribution well described by a Weibull model having a shape parameter of $k = 0.30$, indicative of strong temporal clustering. During an epoch of enhanced activity, the GP rate increased by nearly two orders of magnitude to $124~\mathrm{hr}^{-1}$, with a corresponding shape parameter of $k = 0.47$. This value is similar to that reported for a burst storm from the repeating fast radio burst FRB 20200120E, suggesting possible phenomenological parallels between GPs from compact binary systems and repeating FRBs.

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.

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

Consistency of the Parkes Pulsar Timing Array Signal with a Nanohertz Gravitational-wave Background

Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA) collaboration that the same models used in recent analyses resulted in strong evidence for a common-spectrum process in simulations where none is present. In this work, we introduce a methodology to distinguish pulsar power spectra with the same amplitude from noise power spectra of similar but distinct amplitudes. The former is the signature of a spatially uncorrelated pulsar term of a nanohertz gravitational-wave background, whereas the latter could represent ensemble pulsar noise properties. We test the methodology on simulated data sets. We find that the reported common process in PPTA pulsars is indeed consistent with the spectral feature of a pulsar term. We recommend this methodology as one of the validity tests that the real astrophysical and cosmological backgrounds should pass, as well as for inferences about the spatially uncorrelated component of the background.

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

Evaluating the prevalence of spurious correlations in pulsar timing array datasets

Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pulsar Timing Array data indicate that such a detection may arise spuriously. In this paper, we use simulated pulsar timing array datasets to further test the robustness of the inference methods for spectral and spatial correlations from a GWB. Expanding on our previous results, we find strong support (Bayes factors exceeding $10^5$) for the presence of a common-spectrum noise process in datasets where no common process is present, under a wide range of timing noise prescriptions per pulsar. We show that these results are highly sensitive to the choice of Bayesian priors on timing noise parameters, with priors that more closely match the injected distributions of timing noise parameters resulting in diminished support for a common-spectrum noise process. These results emphasize shortcomings in current methods for inferring the presence of a common-spectrum process, and imply that the detection of a common process is not a reliable precursor to detection of the GWB. Future searches for the nanohertz GWB should remain focussed on detecting spatial correlations, and make use of more tailored specifications for a common-spectrum noise process.

preprint2022arXiv

Evidence for PeV Proton Acceleration from Fermi-LAT Observations of SNR G106.3+2.7

The existence of a "knee" at energy ~1 PeV in the cosmic-ray spectrum suggests the presence of Galactic PeV proton accelerators called "PeVatrons". Supernova Remnant (SNR) G106.3+2.7 is a prime candidate for one of these. The recent detection of very high energy (0.1-100 TeV) gamma rays from G106.3+2.7 may be explained either by the decay of neutral pions or inverse Compton scattering by relativistic electrons. We report an analysis of 12 years of Fermi-LAT gamma-ray data which shows that the GeV-TeV gamma-ray spectrum is much harder and requires a different total electron energy than the radio and X-ray spectra, suggesting it has a distinct, hadronic origin. The non-detection of gamma rays below 10 GeV implies additional constraints on the relativistic electron spectrum. A hadronic interpretation of the observed gamma rays is strongly supported. This observation confirms the long-sought connection between Galactic PeVatrons and SNRs. Moreover, it suggests that G106.3+2.7 could be the brightest member of a new population of SNRs whose gamma-ray energy flux peaks at TeV energies. Such a population may contribute to the cosmic-ray knee and be revealed by future very high energy gamma-ray detectors.

preprint2022arXiv

Exploring the MeV Sky with a Combined Coded Mask and Compton Telescope: The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds on and improves past and currently proposed missions at such energies. We outline the motivations for combining a coded mask and a Compton telescope and we define the scientific goals of such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which enable identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

preprint2022arXiv

Single Pulse Variability in Gamma-ray Pulsars

The Fermi Large Area Telescope receives $\ll$1 photon per rotation from any $γ$-ray pulsar. However, out of the billions of monitored rotations of the bright pulsars Vela (PSR~J0835$-$4510) and Geminga (PSR~J0633$+$1746), a few thousand have $\geq$2 pulsed photons. These rare pairs encode information about the variability of pulse amplitude and shape. We have cataloged such pairs and find the observed number to be in good agreement with simple Poisson statistics, limiting any amplitude variations to $<$19% (Vela) and $<$22% (Geminga) at 2$σ$ confidence. Using an array of basis functions to model pulse shape variability, the observed pulse phase distribution of the pairs limits the scale of pulse shape variations of Vela to $<$13% while for Geminga we find a hint of $\sim$20% single-pulse shape variability most associated with the pulse peaks. If variations last longer than a single rotation, more pairs can be collected, and we have calculated upper limits on amplitude and shape variations for assumed coherence times up to 100 rotations, finding limits of $\sim$1% (amplitude) and $\sim$3% (shape) for both pulsars. Because a large volume of the pulsar magnetosphere contributes to $γ$-ray pulse production, we conclude that the magnetospheres of these two energetic pulsars are stable over one rotation and very stable on longer time scales. All other $γ$-ray pulsars are too faint for similar analyses. These results provide useful constraints on rapidly-improving simulations of pulsar magnetospheres, which have revealed a variety of large-scale instabilities in the thin equatorial current sheets where the bulk of GeV $γ$-ray emission is thought to originate.

preprint2022arXiv

The Future of Gamma-Ray Experiments in the MeV-EeV Range

Gamma-rays, the most energetic photons, carry information from the far reaches of extragalactic space with minimal interaction or loss of information. They bring messages about particle acceleration in environments so extreme they cannot be reproduced on earth for a closer look. Gamma-ray astrophysics is so complementary with collider work that particle physicists and astroparticle physicists are often one in the same. Gamma-ray instruments, especially the Fermi Gamma-ray Space Telescope, have been pivotal in major multi-messenger discoveries over the past decade. There is presently a great deal of interest and scientific expertise available to push forward new technologies, to plan and build space- and ground-based gamma-ray facilities, and to build multi-messenger networks with gamma rays at their core. It is therefore concerning that before the community comes together for planning exercises again, much of that infrastructure could be lost to a lack of long-term planning for support of gamma-ray astrophysics. Gamma-rays with energies from the MeV to the EeV band are therefore central to multiwavelength and multi-messenger studies to everything from astroparticle physics with compact objects, to dark matter studies with diffuse large scale structure. These goals and new discoveries have generated a wave of new gamma-ray facility proposals and programs. This paper highlights new and proposed gamma-ray technologies and facilities that have each been designed to address specific needs in the measurement of extreme astrophysical sources that probe some of the most pressing questions in fundamental physics for the next decade. The proposed instrumentation would also address the priorities laid out in the recent Astro2020 Decadal Survey, a complementary study by the astrophysics community that provides opportunities also relevant to Snowmass.

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

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.

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

Gamma-ray heartbeat powered by the microquasar SS 433

Microquasars, the local siblings of extragalactic quasars, are binary systems comprising a compact object and a companion star. By accreting matter from their companions, microquasars launch powerful winds and jets, influencing the interstellar environment around them. Steady gamma-ray emission is expected to rise from their central objects, or from interactions between their outflows and the surrounding medium. The latter prediction was recently confirmed with the detection of SS 433 at high (TeV) energies. In this report, we analyze more than ten years of GeV gamma-ray data from the Fermi Gamma-ray Space Telescope on this source. Detailed scrutiny of the data reveal emission in the SS 433 vicinity, co-spatial with a gas enhancement, and hints for emission possibly associated with a terminal lobe of one of the jets. Both gamma-ray excesses are relatively far from the central binary, and the former shows evidence for a periodic variation at the precessional period of SS 433, linking it with the microquasar. This result challenges obvious interpretations and is unexpected from any previously published theoretical models. It provides us with a chance to unveil the particle transport from SS 433 and to probe the structure of the local magnetic field in its vicinity.

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.

preprint2019arXiv

Discovery of Millisecond Pulsars in the Globular Cluster Omega Centauri

The globular cluster Omega Centauri is the most massive and luminous cluster in the Galaxy. The $γ$-ray source FL8Y J1326.7$-$4729 is coincident with the core of the cluster, leading to speculation that hitherto unknown radio pulsars or annihilating dark matter may be present in the cluster core. Here we report on the discovery of five millisecond pulsars in Omega Centauri following observations with the Parkes radio telescope. Four of these pulsars are isolated with spin periods of 4.1, 4.2, 4.6 and 6.8 ms. The fifth has a spin period of 4.8 ms and is in an eclipsing binary system with an orbital period of 2.1 hours. Deep radio continuum images of the cluster centre with the Australian Telescope Compact Array reveal a small population of compact radio sources making it likely that other pulsars await discovery. We consider it highly likely that the millisecond pulsars are the source of the $γ$-ray emission. The long-term timing of these pulsars opens up opportunities to explore the dynamics and interstellar medium of the cluster.

preprint2019arXiv

Multi-scale time- and frequency-domain likelihood analysis with photon weights

We present an unbinned likelihood analysis formalism employing photon weights -- the probabilities that events are associated with a particular source. This approach is applicable to any photon-resolving instrument, and thus well suited to high-energy observations; we focus here on GeV $γ$-ray data from the Fermi Large Area Telescope. Weights connect individual photons to the outputs of a detailed, expensive likelihood analysis of a much larger data set. The weighted events can be aggregated into arbitrary time spans ranging from microseconds to years. Such retrospective grouping permits time- and frequency-domain analysis over a wide range of scales and enables characterization of disparate phenomena like blazar flares, $γ$-ray bursts, pulsar pulses, novae, $γ$-ray binaries, and other variable sources. To demonstrate the formalism, we incorporate photon weights into the Bayesian blocks algorithm and perform a hierarchical time scale analysis of 3C 279 activity. We analyze pulsar pulse profiles and estimate the unpulsed emission level and the optimal division of the data into on- and off-pulse intervals. We extend the formalism to Fourier analysis and derive estimators for power spectra, used to search for and characterize periodic sources. We show how the Fast Fourier transform can be used to probe orbital periods as short as a minute and we discuss mitigation of spurious signals. Our final example combines time- and frequency- domain analysis to jointly characterize the flares and orbital modulation of Cygnus X-3, yielding the strongest detection of the orbital signal ($>$13$σ$) to date. Finally, we discuss extensions of the work to other GeV sources and to X-ray and TeV observations.

preprint2019arXiv

Probing the emission states of PSR J1107-5907

The emission from PSR J1107-5907 is erratic. Sometimes the radio pulse is undetectable, at other times the pulsed emission is weak, and for short durations the emission can be very bright. In order to improve our understanding of these state changes, we have identified archival data sets from the Parkes radio telescope in which the bright emission is present, and find that the emission never switches from the bright state to the weak state, but instead always transitions to the off state. Previous work had suggested the identification of the off state as an extreme manifestation of the weak state. However, the connection between the off and bright emission reported here suggests that the emission can be interpreted as undergoing only two emission states: a bursting state consisting of both bright pulses and nulls as well as the weak-emission state.

preprint2018arXiv

Theoretically motivated search and detection of non-thermal pulsations from PSRs J1747-2958, J2021+3651, and J1826-1256

Based on a theoretical selection of pulsars as candidates for detection at X-ray energies, we present an analysis of archival X-ray observations done with Chandra and XMM-Newton of PSR J1747-2958 (the pulsar in the "Mouse" nebula), PSR J2021+3651 (the pulsar in the "Dragonfly" nebula), and PSR J1826-1256. X-ray pulsations from PSR J1747-2958 and PSR J1826-1256 are detected for the first time, and a previously reported hint of an X-ray pulsation from PSR J2021+3651 is confirmed with a higher significance. We analyze these pulsars' spectra in regards to the theoretically predicted energy distribution, finding a remarkable agreement, and provide here a refined calculation of the model parameters taking into account the newly derived X-ray spectral data.

preprint2016arXiv

Gravitational-wave cosmology across 29 decades in frequency

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $Ω_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

preprint2016arXiv

On the difference between gamma-ray-detected and non-gamma-ray-detected pulsars

We compare radio profile widths of young, energetic gamma-ray-detected and non-gamma-ray-detected pulsars. We find that the latter typically have wider radio profiles, with the boundary between the two samples exhibiting a dependence on the rate of rotational energy loss. We also find that within the sample of gamma-ray-detected pulsars, radio profile width is correlated with both the separation of the main gamma-ray peaks and the presence of narrow gamma-ray components. These findings lead us to propose that these pulsars form a single population where the main factors determining gamma ray detectability are the rate of rotational energy loss and the proximity of the line of sight to the rotation axis. The expected magnetic inclination angle distribution will be different for radio pulsars with and without detectable gamma rays, naturally leading to the observed differences. Our results also suggest that the geometry of existing radio and outer-magnetosphere gamma-ray emission models are at least qualitatively realistic, implying that information about the viewing geometry can be extracted from profile properties of pulsars.

preprint2016arXiv

Young and Millisecond Pulsar GeV Gamma-ray Fluxes from the Galactic Center and Beyond

Gamma-ray observations have shown pulsars to be efficient converters of rotational energy into GeV photons and it is of wide-ranging interest to determine their contribution to the gamma-ray background. We arrive at flux predictions from both the young (<~ Myr) and millisecond (~Gyr) Galactic pulsar populations. We find that unresolved pulsars can yield both a significant fraction of the excess GeV gamma rays near the Galactic Center and an inverse Compton flux in the inner kpc similar to that inferred by Fermi. We compare models of the young pulsar population and millisecond pulsar population to constraints from gamma-ray and radio observations. Overall, we find that the young pulsars should outnumber millisecond pulsars as unassociated gamma-ray point sources in this region. The number of young radio pulsars discovered near the Galactic Center is in agreement with our model of the young pulsar population. Deeper radio observations at higher latitudes can constrain the total gamma-ray emission from both young and millisecond pulsars from the inner galaxy. While this is a step towards better understanding of pulsars, cosmic rays in the Milky Way, and searches for dark matter, we also discuss a few interesting puzzles that arise from the underlying physics of pulsar emission and evolution.

preprint2015arXiv

A study of spatial correlations in pulsar timing array data

Pulsar timing array experiments search for phenomena that produce angular correlations in the arrival times of signals from millisecond pulsars. The primary goal is to detect an isotropic and stochastic gravitational wave background. We use simulated data to show that this search can be affected by the presence of other spatially correlated noise, such as errors in the reference time standard, errors in the planetary ephemeris, the solar wind and instrumentation issues. All these effects can induce significant false detections of gravitational waves. We test mitigation routines to account for clock errors, ephemeris errors and the solar wind. We demonstrate that it is non-trivial to find an effective mitigation routine for the planetary ephemeris and emphasise that other spatially correlated signals may be present in the data.

preprint2015arXiv

Detection and localization of single-source gravitational waves with pulsar timing arrays

Pulsar timing arrays (PTAs) can be used to search for very low frequency ($10^{-9}$--$10^{-7}$ Hz) gravitational waves (GWs). In this paper we present a general method for the detection and localization of single-source GWs using PTAs. We demonstrate the effectiveness of this new method for three types of signals: monochromatic waves as expected from individual supermassive binary black holes in circular orbits, GWs from eccentric binaries and GW bursts. We also test its implementation in realistic data sets that include effects such as uneven sampling and heterogeneous data spans and measurement precision. It is shown that our method, which works in the frequency domain, performs as well as published time-domain methods. In particular, we find it equivalent to the $\mathcal{F}_{e}$-statistic for monochromatic waves. We also discuss the construction of null streams -- data streams that have null response to GWs, and the prospect of using null streams as a consistency check in the case of detected GW signals. Finally, we present sensitivities to individual supermassive binary black holes in eccentric orbits. We find that a monochromatic search that is designed for circular binaries can efficiently detect eccentric binaries with both high and low eccentricities, while a harmonic summing technique provides greater sensitivities only for binaries with moderate eccentricities.

preprint2015arXiv

Improving Pulsar Timing Precision with Single Pulses

The measurement error of pulse times of arrival (TOAs) in the high S/N limit is dominated by the quasi-random variation of a pulsar's emission profile from rotation to rotation. Like measurement noise, this noise is only reduced as the square root of observing time, posing a major challenge to future pulsar timing campaigns with large aperture telescopes, e.g. the Five-hundred-metre Aperture Spherical Telescope and the Square Kilometre Array. We propose a new method of pulsar timing that attempts to approximate the pulse-to-pulse variability with a small family of 'basis' pulses. If pulsar data are integrated over many rotations, this basis can be used to measure sub-pulse structure. Or, if high-time resolution data are available, the basis can be used to 'tag' single pulses and produce an optimal timing template. With realistic simulations, we show that these applications can dramatically reduce the effect of pulse-to-pulse variability on TOAs. Using high-time resolution data taken from the bright PSR J0835-4510 (Vela), we demonstrate a 25-40% improvement in TOA precision. Crucially for pulsar timing applications, we further establish that these techniques produce TOAs with gaussian residuals. Improvements of this level halve the telescope time required to reach a desired TOA precision. Although some gains can be achieved with existing data, the greatest improvements result from the 'tagging' approach, which in turn requires online or posthoc analysis of single pulses, an important consideration for the design of future instrumentation.

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

Limits on Planet Formation Around Young Pulsars and Implications for Supernova Fallback Disks

We have searched a sample of 151 young, energetic pulsars for periodic variation in pulse time-of-arrival arising from the influence of planetary companions. We are sensitive to objects with masses two orders of magnitude lower than those detectable with optical transit timing, but we find no compelling evidence for pulsar planets. For the older pulsars most likely to host planets, we can rule out Mercury analogues in one third of our sample and planets with masses $>0.4 M_{\oplus}$ and periods $P_b<1$ yr in all but 5% of such systems. If pulsar planets form primarily from supernova fallback disks, these limits imply that such disks do not form, are confined to $<0.1$ AU radii, are disrupted, or form planets more slowly ($>2$ Myr) than their protoplanetary counterparts.

preprint2015arXiv

Periodic modulation in pulse arrival times from young pulsars: a renewed case for neutron star precession

In a search for periodic variation in the arrival times of pulses from 151 young, energetic pulsars, we have identified seven cases of modulation consistent with one or two harmonics of a single fundamental with time-scale 0.5-1.5 yr. We use simulations to show that these modulations are statistically significant and of high quality (sinusoidal) even when contaminated by the strong stochastic timing noise common to young pulsars. Although planetary companions could induce such modulation, the large implied masses and 2:1 mean motion resonances challenge such an explanation. Instead, the modulation is likely to be intrinsic to the pulsar, arising from quasi-periodic switching between stable magnetospheric states, and we propose that precession of the neutron star may regulate this switching.

preprint2015arXiv

The 1st Fermi Lat Supernova Remnant Catalog

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

preprint2015arXiv

The Binary Nature of PSR J2032+4127

PSR J2032+4127 is a gamma-ray and radio-emitting pulsar which has been regarded as a young luminous isolated neutron star. However, its recent spin-down rate has extraordinarily increased by a factor of two. We present evidence that this is due to its motion as a member of a highly-eccentric binary system with a 15-solar-mass Be star, MT91~213. Timing observations show that, not only are the positions of the two stars coincident within 0.4 arcsec, but timing models of binary motion of the pulsar fit the data much better than a model of a young isolated pulsar. MT91~213, and hence the pulsar, lie in the Cyg~OB2 stellar association, which is at a distance of only 1.4-1.7 kpc. The pulsar is currently on the near side of, and accelerating towards, the Be star, with an orbital period of 20-30 years. The next periastron is well-constrained to occur in early 2018, providing an opportunity to observe enhanced high-energy emission as seen in other Be-star binary systems.

preprint2015arXiv

Timing Gamma-ray Pulsars with the Fermi Large Area Telescope: Timing Noise and Astrometry

We have constructed timing solutions for 81 gamma-ray pulsars covering more than five years of Fermi data. The sample includes 37 radio-quiet or radio-faint pulsars which cannot be timed with other telescopes. These timing solutions and the corresponding pulse times of arrival are prerequisites for further study, e.g. phase-resolved spectroscopy or searches for mode switches. Many gamma-ray pulsars are strongly affected by timing noise, and we present a new method for characterizing the noise process and mitigating its effects on other facets of the timing model. We present an analysis of timing noise over the population using a new metric for characterizing its strength and spectral shape, namely its time-domain correlation. The dependence of the strength on spin frequency and spin-down rate is in good agreement with previous studies. We find that noise process power spectra $S(f)$ for unrecycled pulsars are steep, with strong correlations over our entire data set and spectral indices $S(f)\propto f^{-α}$ of 5 to 9. One possible explanation for these results is the occurrence of unmodelled, episodic 'microglitches'. Finally, we show that our treatment of timing noise results in robust parameter estimation, and in particular we measure a precise timing position for each pulsar. We extensively validate our results with multi-wavelength astrometry, and using our updated position, we firmly identify the X-ray counterpart of PSR J1418-6058.

preprint2015arXiv

Young Pulsars and the Galactic Center GeV Gamma-ray Excess

Studies of Fermi data indicate an excess of GeV gamma rays around the Galactic center (GC), possibly due to dark matter. We show that young gamma-ray pulsars can yield a similar signal. First, a high concentration of GC supernovae naturally leads to a population of kicked pulsars symmetric about the GC. Second, while very-young pulsars with soft spectra reside near the Galactic plane, pulsars with spectra that have hardened with age accumulate at larger angles. This combination, including unresolved foreground pulsars, traces the morphology and spectrum of the Excess.

preprint2014arXiv

The Three Discrete Nulling Timescales of PSR J1717-4054

PSR J1717-4054 is one of the small class of pulsars which null on intermediate (~hour) timescales. Such pulsars represent an important link between classical nullers, whose emission vanishes for a few rotations, and the intermittent pulsars which null for months and years. Using the Parkes radio telescope and the Australia Telescope Compact Array, we have studied the emission from PSR J1717-4054 over intervals from single pulses to years. We have identified and characterised nulling at three discrete timescales: the pulsar emits during "active states" separated by nulls lasting many thousands of rotations, while active states themselves are interrupted by nulls with a bimodal distribution of durations; one to two rotations, or tens of rotations. We detect no subpulse structure or flux density variations during active states, and we place stringent limits (<0.1% of the mean active flux density) on pulsed emission during inactive states. Finally, our high-quality data have also allowed us to measure for the first time many important properties of PSR J1717-4054, including its position, spindown rate, spectrum, polarization characteristics, and pulse broadening in the interstellar medium.

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

Chandra Observations of Black-Widow Pulsars

We describe the first X-ray observations of binary millisecond pulsars PSRs J0023+0923, J1810+1744, J2215+5135, and J2256-1024. All four are Fermi gamma-ray sources and three are 'black-widow' pulsars, with companions of mass < 0.1 solar masses. Data were taken using the Chandra X-Ray Observatory and covered a full binary orbit for each pulsar. Two pulsars, PSRs J2215+5135 and J2256-1024, show significant orbital variability and X-ray flux minima at the times of eclipses observed at radio wavelengths. This phenomenon is consistent with intrabinary shock emission characteristic of black-widow pulsars. The other two pulsars, PSRs J0023+0923 and J1810+1744, do not demonstrate significant variability, but are fainter than the other two sources. Spectral fits yield power-law indices that range from 1.4 to 2.3 and blackbody temperatures in the hundreds of eV. The spectrum for PSR J2215+5135 shows a significant hard X-ray component (41% of counts are above 2 keV), which is additional evidence for the presence of intrabinary shock emission.

preprint2012arXiv

Pulsars in Gamma Rays: What Fermi Is Teaching Us

The 2nd Fermi-LAT pulsar catalog includes 117 gamma-ray pulsars, of which roughly one third are millisecond pulsars (MSPs) while the remaining two thirds split evenly into young radio-loud and radio-quiet pulsars. Although this large population will enable future, detailed studies of emission mechanisms and the evolution of the underlying neutron star population, some nearly-universal properties are already clear and unequivocal. We discuss some of these aspects below, including the altitude of the gamma-ray emission site and the shape of the gamma-ray spectrum and its implications for the radiation mechanism.

preprint2012arXiv

Search for Spatially Extended Fermi-LAT Sources Using Two Years of Data

Spatial extension is an important characteristic for correctly associating gamma-ray-emitting sources with their counterparts at other wavelengths and for obtaining an unbiased model of their spectra. We present a new method for quantifying the spatial extension of sources detected by the Large Area Telescope (LAT), the primary science instrument on the Fermi Gamma-ray Space Telescope (Fermi). We perform a series of Monte Carlo simulations to validate this tool and calculate the LAT threshold for detecting the spatial extension of sources. We then test all sources in the second Fermi-LAT catalog (2FGL) for extension. We report the detection of seven new spatially extended sources.

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

Improving Sensitivity to Weak Pulsations with Photon Probability Weighting

All gamma-ray telescopes suffer from source confusion due to their inability to focus incident high-energy radiation, and the resulting background contamination can obscure the periodic emission from faint pulsars. In the context of the Fermi Large Area Telescope, we outline enhanced statistical tests for pulsation in which each photon is weighted by its probability to have originated from the candidate pulsar. The probabilities are calculated using the instrument response function and a full spectral model, enabling powerful background rejection. With Monte Carlo methods, we demonstrate that the new tests increase the sensitivity to pulsars by more than 50% under a wide range of conditions. This improvement may appreciably increase the completeness of the sample of radio-loud gamma-ray pulsars. Finally, we derive the asymptotic null distribution for the H-test, expanding its domain of validity to arbitrarily complex light curves.

preprint2011arXiv

Likelihood Methods for the Detection and Characterization of Gamma-ray Pulsars with the Fermi Large Area Telescope

The sensitivity of the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope allows detection of thousands of new gamma-ray sources and detailed characterization of the spectra and variability of bright sources. Unsurprisingly, this increased capability leads to increased complexity in data analysis. Likelihood methods are ideal for connecting models with data, but the computational cost of folding the model input through the multi-scale instrument response function is appreciable. Both interactive analysis and large projects---such as analysis of the full gamma-ray sky---can be prohibitive or impossible, reducing the scope of the science possible with the LAT. To improve on this situation, we have developed pointlike, a software package for fast maximum likelihood analysis of LAT data. It is interactive by design and its rapid evaluation of the likelihood facilitates exploratory and large-scale, all-sky analysis. We detail its implementation and validate its performance on simulated data. We demonstrate its capability for interactive analysis and present several all-sky analyses. These include a search for new gamma-ray sources and the selection of LAT sources with pulsar-like characteristics for targeted radio pulsation searches. We conclude by developing sensitive periodicity tests incorporating spectral information obtained from pointlike.

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

Hunting for New Gamma-ray Binaries - Technique Development

There are only a few sources that are definitely known to be gamma-ray binaries. Two of these are listed as associations in the Fermi LAT Bright Source List. We are developing novel techniques to extract high signal-to-noise light curves of all cataloged Fermi sources and to search for periodic variability using appropriately weighted power spectra. The detection of periodic variability would be strong evidence for the detection of a new gamma-ray binary. The LAT's sensitivity provides the opportunity to open up completely new discovery space for additional binary systems, potentially involving novel astrophysics. We present here demonstrations of the sensitivity gains obtained through the use of these techniques.

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.