Source author record

A. K. Harding

A. K. Harding appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph.GA astro-ph.SR astro-ph.IM

Catalog footprint

What is connected

73works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background

After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

preprint2022arXiv

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.

preprint2022arXiv

Pulse Peak Migration during the Outburst Decay of the Magnetar SGR 1830-0645: Crustal Motion and Magnetospheric Untwisting

Magnetars, isolated neutron stars with magnetic field strengths typically $\gtrsim10^{14}$~G, exhibit distinctive months-long outburst epochs during which strong evolution of soft X-ray pulse profiles, along with nonthermal magnetospheric emission components, is often observed. Using near-daily NICER observations of the magnetar SGR 1830-0645 during the first 37 days of a recent outburst decay, a pulse peak migration in phase is clearly observed, transforming the pulse shape from an initially triple-peaked to a single-peaked profile. Such peak merging has not been seen before for a magnetar. Our high-resolution phase-resolved spectroscopic analysis reveals no significant evolution of temperature despite the complex initial pulse shape. Yet the inferred surface hot spots shrink during the peak migration and outburst decay. We suggest two possible origins for this evolution. For internal heating of the surface, tectonic motion of the crust may be its underlying cause. The inferred speed of this crustal motion is $\lesssim100$~m~day$^{-1}$, constraining the density of the driving region to $ρ\sim10^{10}$~g~cm$^{-3}$, at a depth of $\sim200$~m. Alternatively, the hot spots could be heated by particle bombardment from a twisted magnetosphere possessing flux tubes or ropes, somewhat resembling solar coronal loops, that untwist and dissipate on the 30-40~day timescale. The peak migration may then be due to a combination of field-line footpoint motion (necessarily driven by crustal motion) and evolving surface radiation beaming. These novel dataset paints a vivid picture of the dynamics associated with magnetar outbursts, yet it also highlights the need for a more generic theoretical picture where magnetosphere and crust are considered in tandem.

preprint2022arXiv

Search for new cosmic-ray acceleration sites within the 4FGL catalog Galactic plane sources

Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Supernova Remnants (SNRs), IC 443, W44, W49B and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5 degrees from the Galactic plane. Using 8 years of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton-proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS~I +61 303; the colliding wind binary eta Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites.

preprint2022arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

preprint2022arXiv

X-ray burst and persistent emission properties of the magnetar SGR 1830-0645 in outburst

We report on NICER X-ray monitoring of the magnetar SGR 1830-0645 covering 223 days following its October 2020 outburst, as well as Chandra and radio observations. We present the most accurate spin ephemerides of the source so far: $ν=0.096008680(2)$~Hz, $\dotν=-6.2(1)\times10^{-14}$~Hz~s$^{-1}$, and a significant second and third frequency derivative terms indicative of non-negligible timing noise. The phase-averaged 0.8--7~keV spectrum is well fit with a double-blackbody (BB) model throughout the campaign. The BB temperatures remain constant at 0.46 and 1.2 keV. The areas and flux of each component decreased by a factor of 6, initially through a steep decay trend lasting about 46 days followed by a shallow long-term one. The pulse shape in the same energy range is initially complex, exhibiting three distinct peaks, yet with clear continuous evolution throughout the outburst towards a simpler, single-pulse shape. The rms pulsed fraction is high and increases from about 40% to 50%. We find no dependence of pulse shape or fraction on energy. These results suggest that multiple hotspots, possibly possessing temperature gradients, emerged at outburst-onset, and shrank as the outburst decayed. We detect 84 faint bursts with \nicer, having a strong preference for occurring close to the surface emission pulse maximum the first time this phenomenon is detected in such a large burst sample. This likely implies a very low altitude for the burst emission region, and a triggering mechanism connected to the surface active zone. Finally, our radio observations at several epochs and multiple frequencies reveal no evidence of pulsed or burst-like radio emission.

preprint2021arXiv

Broadband X-ray Burst Spectroscopy of the FRB-Emitting Galactic Magnetar

Magnetars are young, magnetically-powered neutron stars possessing the strongest magnetic fields in the Universe. Fast Radio Bursts (FRBs) are extremely intense millisecond-long radio pulses of primarily extragalactic origin, and a leading attribution for their genesis focuses on magnetars. A hallmark signature of magnetars is their emission of bright, hard X-ray bursts of sub-second duration. On April 27th 2020, the Galactic magnetar SGR J1935+2154 emitted hundreds of X-ray bursts in a few hours. One of these temporally coincided with an FRB, the first detection of an FRB from the Milky Way. Here we present spectral and temporal analyses of 24 X-ray bursts emitted 13 hours prior to the FRB and seen simultaneously with the NASA NICER and Fermi/GBM missions in their combined energy range, 0.2 keV-30 MeV. These broadband spectra permit direct comparison with the spectrum of the FRB-associated X-ray burst (FRB-X). We demonstrate that all 24 NICER/GBM bursts are very similar temporally, albeit strikingly different spectrally, from FRB-X. The singularity of the FRB-X burst is perhaps indicative of an uncommon locale for its origin. We suggest that this event originated in quasi-polar open or closed magnetic field lines that extend to high altitudes.

preprint2021arXiv

Einstein@Home Discovery of the Gamma-ray Millisecond Pulsar PSR J2039-5617 Confirms Its Predicted Redback Nature

The Fermi Large Area Telescope gamma-ray source 3FGL J2039.6$-$5618 contains a periodic optical and X-ray source that was predicted to be a "redback" millisecond pulsar (MSP) binary system. However, the conclusive identification required the detection of pulsations from the putative MSP. To better constrain the orbital parameters for a directed search for gamma-ray pulsations, we obtained new optical light curves in 2017 and 2018, which revealed long-term variability from the companion star. The resulting orbital parameter constraints were used to perform a targeted gamma-ray pulsation search using the Einstein@Home distributed volunteer computing system. This search discovered pulsations with a period of 2.65 ms, confirming the source as a binary MSP now known as PSR J2039$-$5617. Optical light curve modelling is complicated, and likely biased, by asymmetric heating on the companion star and long-term variability, but we find an inclination $i > 60°$, for a low pulsar mass between $1.1 M_{\odot} < M_{\rm psr} < 1.6 M_{\odot}$ and a companion mass of 0.15--0.22 $M_{\odot}$, confirming the redback classification. Timing the gamma-ray pulsations also revealed significant variability in the orbital period, which we find to be consistent with quadrupole moment variations in the companion star, suggestive of convective activity. We also find that the pulsed flux is modulated at the orbital period, potentially due to inverse Compton scattering between high-energy leptons in the pulsar wind and the companion star's optical photon field.

preprint2020arXiv

A radiatively-quiet glitch and anti-glitch in the magnetar 1E 2259+586

We report on the timing and spectral properties of the soft X-ray emission from the magnetar 1E 2259+586 from January 2013, $\sim 8$ months after the detection of an anti-glitch, until September 2019, using the Neil Gehrels Swift and NICER observatories. During this time span, we detect two timing discontinuities. The first, occurring around 5 years after the April 2012 anti-glitch, is a relatively large spin-up glitch with a fractional amplitude $Δν/ν=1.24(2)\times10^{-6}$. We find no evidence for flux enhancement or change in the spectral or pulse profile shape around the time of this glitch. This is consistent with the picture that a significant number of magnetar spin-up glitches are radiatively-quiet. Approximately 1.5 years later in April 2019, 1E 2259+586 exhibited an anti-glitch with spin-down of a fractional amplitude $Δν/ν=-5.8(1)\times10^{-7}$; similar to the fractional change detected in 2012. We do not, however, detect any change to the pulse-profile shape or increase in the rms pulsed flux of the source, nor do we see any possible bursts from its direction around the time of the anti-glitch; all of which occurred during the 2012 event. Hence, similar to spin-up glitches, anti-glitches can occur silently. This may suggest that these phenomena originate in the neutron star interior, and that their locale and triggering mechanism do not necessarily have to be connected to the magnetosphere. Lastly, our observations suggest that the occurrence rate of spin-up and spin-down glitches is about the same in 1E 2259+586, with the former having a larger net fractional change.

preprint2020arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's rotation, and from the $l = 2$, $m = 1$ mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the $l = m = 2$ mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars' spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711-6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of $10^{38}$ kg m$^2$ for the star's moment of inertia, and imply a gravitational-wave-derived upper limit on the star's ellipticity of $1.2\!\times\!10^{-8}$. We also place new limits on the emission amplitude at the rotation frequency of the pulsars.

preprint2016arXiv

3FGL Demographics Outside the Galactic Plane using Supervised Machine Learning: Pulsar and Dark Matter Subhalo Interpretations

Nearly 1/3 of the sources listed in the Third Fermi Large Area Telescope (LAT) catalog (3FGL) remain unassociated. It is possible that predicted and even unanticipated gamma-ray source classes are present in these data waiting to be discovered. Taking advantage of the excellent spectral capabilities achieved by the Fermi LAT, we use machine learning classifiers (Random Forest and XGBoost) to pinpoint potentially novel source classes in the unassociated 3FGL sample outside the Galactic plane. Here we report a total of 34 high-confidence Galactic candidates at |b| > 5 degrees. The currently favored standard astrophysical interpretations for these objects are pulsars or low-luminosity globular clusters hosting millisecond pulsars (MSPs). Yet, these objects could also be interpreted as dark matter annihilation taking place in ultra-faint dwarf galaxies or dark matter subhalos. Unfortunately, Fermi LAT spectra are not sufficient to break degeneracies between the different scenarios. Careful visual inspection of archival optical images reveals no obvious evidence for low-luminosity globular clusters or ultra-faint dwarf galaxies inside the 95% error ellipses. If these are pulsars, this would bring the total number of MSPs at |b| > 5 degrees to 106. We find this number to be in excellent agreement with predictions from a new population synthesis of MSPs that predicts 100-126 high-latitude 3FGL MSPs depending on the choice of high-energy emission model. If, however, these are dark matter substructures, we can place upper limits on the number of Galactic subhalos surviving today and on dark matter annihilation cross sections. These limits are beginning to approach the canonical thermal relic cross section for dark matter particle masses below ~100 GeV in the bottom quark annihilation channel.

preprint2016arXiv

A New, Low Braking Index For the LMC Pulsar B0540-69

We report the results of a 16-month monitoring campaign using the Swift satellite of PSR B0540-69, a young pulsar in the Large Magellanic Cloud. Phase connection was maintained throughout the campaign so that a reliable ephemeris could be determined, and the length of the campaign is adequate to accurately determine the spin frequency and its first and second derivatives. The braking index is 0.031 +/- 0.013 (90% confidence), a value much lower than previously reported for B0540-69 and almost all other young pulsars. We use data from the extensive monitoring campaign with RXTE to show that timing noise is unlikely to significantly affect the measurement. This is the first measurement of the braking index in the pulsar's recently discovered high spin-down state. We discuss possible mechanisms for producing the low braking index.

preprint2016arXiv

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center.

preprint2016arXiv

Implementation of a goodness-of-fit test for finding optimal concurrent radio and γ-ray pulsar light curves

Since the launch of the Fermi Large Area Telescope in 2008 the number of known $γ$-ray pulsars has increased immensely to over 200, many of which are also visible in the radio and X-ray bands. Seyffert et al. (2011) demonstrated how constraints on the viewing geometries of some of these pulsars could be obtained by comparing their observed radio and $γ$-ray light curves by eye to light curves from geometric models. While these constraints compare reasonably well with those yielded by more rigorous single-wavelength approaches, they are still a somewhat subjective representation of how well the models reproduce the observed radio and $γ$-ray light curves. Constructing a more rigorous approach is, however, made difficult by the large uncertainties associated with the $γ$-ray light curves as compared to those associated with the radio light curves. Naively applying a $χ^{2}$-like goodness-of-fit test to both bands invariably results in constraints dictated by the radio light curves. A number of approaches have been proposed to address this issue. In this paper we investigate these approaches and evaluate the results they yield. Based on what we learn, we implement our own version of a goodness-of-fit test, which we then use to investigate the behaviour of the geometric models in multi-dimensional phase space.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Measurement of the high-energy gamma-ray emission from the Moon with the Fermi Large Area Telescope

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

preprint2016arXiv

Pulsar emission in the very-high-energy regime

The vast majority of the pulsars detected by the Fermi Large Area Telescope (LAT) display spectra with exponential cutoffs falling in a narrow range around a few GeV. Early spectral modelling predicted spectral cutoff energies of up to 100 GeV. More modern studies estimated spectral cutoff energies in the 1-20 GeV range. It was therefore not expected that pulsars would be visible in the very-high-energy (VHE; >100 GeV) regime. The VERITAS detection (confirmed by MAGIC) of pulsed emission from the Crab pulsar up to 400 GeV (and now possibly up to 1 TeV) therefore raised important questions about our understanding of the electrodynamics and local environment of pulsars. H.E.S.S. has now detected pulsed emission from the Vela pulsar in the 20-120 GeV range, making this the second pulsar detected by a ground-based Cherenkov telescope. We will review the latest developments in VHE pulsar science, including an overview of recent observations and refinements to radiation models and magnetic field structures. This will assist us in interpreting the VHE emission detected from the Crab and Vela pulsars, and predicting the level of VHE emission expected from other pulsars, which will be very important for the upcoming CTA.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2016arXiv

Very-high energy emission from pulsars

The vast majority of pulsars detected by the Fermi Large Area Telescope (LAT) display exponentially cutoff spectra with cutoffs falling in a narrow band around a few GeV. Early spectral modelling predicted spectral cutoffs at energies of up to 100 GeV, assuming curvature radiation. It was therefore not expected that pulsars would be visible in the very-high energy (VHE) regime (>100 GeV). The VERITAS announcement of the detection of pulsed emission from the Crab pulsar at energies up to 400 GeV (and now up to 1.5 TeV as detected by MAGIC) therefore raised important questions about our understanding of the electrodynamics and local environment of pulsars. H.E.S.S. has now detected pulsed emission from the Vela pulsar down to tens of GeV, making this the second pulsar detected by a ground-based Cherenkov telescope. Deep upper limits have also been obtained by VERITAS and MAGIC for the Geminga pulsar. We will review the latest developments in VHE pulsar science, including an overview of the latest observations, refinements, and extensions to radiation models and magnetic field structures, and the implementation of new radiation mechanisms. This will assist us in understanding the VHE emission detected from the Crab pulsar, and predicting the level of VHE emission expected from other pulsars, which is very important for the upcoming CTA.

preprint2016arXiv

Young and middle age pulsar light-curve morphology: Comparison of Fermi observations with gamma-ray and radio emission geometries

Thanks to the huge amount of gamma-ray pulsar photons collected by the Fermi Large Area Telescope since June 2008, it is now possible to constrain gamma-ray geometrical models by comparing simulated and observed light-curve morphological characteristics. We assumed vacuum-retarded dipole pulsar magnetic field and tested simulated and observed morphological light-curve characteristics in the framework of two pole emission geometries, Polar Cap (PC), radio, and Slot Gap (SG), and Outer Gap (OG)/One Pole Caustic (OPC) emission geometries. We compared simulated and observed/estimated light-curve morphological parameters as a function of observable and non-observable pulsar parameters. The PC model gives the poorest description of the LAT pulsar light-curve morphology. The OPC best explains both the observed gamma-ray peak multiplicity and shape classes. The OPC and SG models describe the observed gamma-ray peak-separation distribution for low- and high-peak separations, respectively. This suggests that the OPC geometry best explains the single-peak structure but does not manage to describe the widely separated peaks predicted in the framework of the SG model as the emission from the two magnetic hemispheres. The OPC radio-lag distribution shows higher agreement with observations suggesting that assuming polar radio emission, the gamma-ray emission regions are likely to be located in the outer magnetosphere. The larger agreement between simulated and LAT estimations in the framework of the OPC suggests that the OPC model best predicts the observed variety of profile shapes. The larger agreement between observations and the OPC model jointly with the need to explain the abundant 0.5 separated peaks with two-pole emission geometries, calls for thin OPC gaps to explain the single-peak geometry but highlights the need of two-pole caustic emission geometry to explain widely separated peaks.

preprint2015arXiv

Cosmic-ray positrons from millisecond pulsars

Observations by the Fermi Large Area Telescope of gamma-ray millisecond pulsar light curves imply copious pair production in their magnetospheres, and not exclusively in those of younger pulsars. Such pair cascades may be a primary source of Galactic electrons and positrons, contributing to the observed enhancement in positron flux above ~10 GeV. Fermi has also uncovered many new millisecond pulsars, impacting Galactic stellar population models. We investigate the contribution of Galactic millisecond pulsars to the flux of terrestrial cosmic-ray electrons and positrons. Our population synthesis code predicts the source properties of present-day millisecond pulsars. We simulate their pair spectra invoking an offset-dipole magnetic field. We also consider positrons and electrons that have been further accelerated to energies of several TeV by strong intrabinary shocks in black widow and redback systems. Since millisecond pulsars are not surrounded by pulsar wind nebulae or supernova shells, we assume that the pairs freely escape and undergo losses only in the intergalactic medium. We compute the transported pair spectra at Earth, following their diffusion and energy loss through the Galaxy. The predicted particle flux increases for non-zero offsets of the magnetic polar caps. Pair cascades from the magnetospheres of millisecond pulsars are only modest contributors around a few tens of GeV to the lepton fluxes measured by AMS-02, PAMELA, and Fermi, after which this component cuts off. The contribution by black widows and redbacks may, however, reach levels of a few tens of percent at tens of TeV, depending on model parameters.

preprint2015arXiv

Discovery of a Spin-Down State Change in the LMC Pulsar B0540-69

We report the discovery of a large, sudden, and persistent increase in the spin-down rate of B0540-69, a young pulsar in the Large Magellanic Cloud, using observations from the Swift and RXTE satellites. The relative increase in the spin-down rate of 36% is unprecedented for B0540-69. No accompanying change in the spin rate is seen, and no change is seen in the pulsed X-ray emission from B0540-69 following the change in the spin-down rate. Such large relative changes in the spin-down rate are seen in the recently discovered class of 'intermittent pulsars', and we compare the properties of B0540-69 to such pulsars. We consider possible changes in the magnetosphere of the pulsar that could cause such a large change in the spin-down rate.

preprint2015arXiv

Discovery of Gamma-ray Pulsations from the Transitional Redback PSR J1227-4853

The 1.69 ms spin period of PSR J1227-4853 was recently discovered in radio observations of the low-mass X-ray binary XSS J12270-4859 following the announcement of a possible transition to a rotation-powered millisecond pulsar state, inferred from decreases in optical, X-ray, and gamma-ray flux from the source. We report the detection of significant (5$σ$) gamma-ray pulsations after the transition, at the known spin period, using ~1 year of data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The gamma-ray light curve of PSR J1227-4853 can be fit by one broad peak, which occurs at nearly the same phase as the main peak in the 1.4 GHz radio profile. The partial alignment of light-curve peaks in different wavebands suggests that at least some of the radio emission may originate at high altitude in the pulsar magnetosphere, in extended regions co-located with the gamma-ray emission site. We folded the LAT data at the orbital period, both pre- and post-transition, but find no evidence for significant modulation of the gamma-ray flux. Analysis of the gamma-ray flux over the mission suggests an approximate transition time of 2012 November 30. Continued study of the pulsed emission and monitoring of PSR J1227-4853, and other known redback systems, for subsequent flux changes will increase our knowledge of the pulsar emission mechanism and transitioning systems.

preprint2015arXiv

Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113

We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

preprint2015arXiv

PSR J1906+0722: An Elusive Gamma-ray Pulsar

We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($Δf / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this.

preprint2015arXiv

Search for Early Gamma-ray Production in Supernovae Located in a Dense Circumstellar Medium with the Fermi LAT

Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We search for a gamma-ray excess at each SNe location in a one year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months and 3 months). For the most promising source of the sample, SN 2010jl (PTF10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in gamma rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the gamma-ray luminosity and the ratio of gamma-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic gamma-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at 95% confidence level (CL) for the source SN 2010jl (PTF10aaxf).

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 anatomy of $γ$-ray pulsar light curves

We previously obtained constraints on the viewing geometries of 6 Fermi LAT pulsars using a multiwavelength approach (Seyffert et al., 2011). To obtain these constraints we compared the observed radio and $γ$-ray light curves (LCs) for those 6 pulsars by eye to LCs predicted by geometric models detailing the location and extent of emission regions in a pulsar magnetosphere. As a precursor to obtaining these constraints, a parameter study was conducted to reinforce our qualitative understanding of how the underlying model parameters effect the LCs produced by the geometric models. Extracting useful trends from the $γ$-ray model LCs proved difficult though due to the increased complexity of the geometric models for the $γ$-ray emission relative to those for the radio emission. In this paper we explore a second approach to investigating the interplay between the model parameters and the LC atlas. This approach does not attempt to understand how the set of model parameters influences the LC shapes directly, but rather, more fundamentally, investigates how the set of model parameters effects the sky maps from which the latter are extracted. This allows us to also recognise structure within the atlas itself, as we are now able to attribute certain features of the LCs to specific features on the sky map, meaning that we not only understand how the structure of single LCs come about, but also how their structure changes as we move through the geometric solution space.

preprint2015arXiv

The Contribution of Millisecond Pulsars to the Galactic Cosmic-Ray Lepton Spectrum

Pulsars are believed to be sources of relativistic electrons and positrons. The abundance of detections of gamma-ray millisecond pulsars by Fermi Large Area Telescope coupled with their light curve characteristics that imply copious pair production in their magnetospheres, motivated us to investigate this old pulsar population as a source of Galactic electrons and positrons and their contribution to the enhancement in cosmic-ray positron flux at GeV energies. We use a population synthesis code to predict the source properties (number, position, and power) of the present-day Galactic millisecond pulsars, taking into account the latest Fermi and radio observations to calibrate the model output. Next, we simulate pair cascade spectra from these pulsars using a model that invokes an offset-dipole magnetic field. We assume free escape of the pairs from the pulsar environment. We then compute the cumulative spectrum of transported electrons and positrons at Earth, following their diffusion and energy losses as they propagate through the Galaxy. Our results indicate that the predicted particle flux increases for non-zero offsets of the magnetic polar caps. Comparing our predicted spectrum and positron fraction to measurements by AMS-02, PAMELA, and Fermi, we find that millisecond pulsars are only modest contributors at a few tens of GeV, after which this leptonic spectral component cuts off. The positron fraction is therefore only slightly enhanced above 10 GeV relative to a background flux model. This implies that alternative sources such as young, nearby pulsars and supernova remnants should contribute additional primary positrons within the astrophysical scenario.

preprint2015arXiv

The effect of an offset-dipole magnetic field on the Vela pulsar's gamma-ray light curves

Over the past six years, the Fermi Large Area Telescope has detected more than 150 gamma-ray pulsars, discovering a variety of light curve trends and classes. Such diversity hints at distinct underlying magnetospheric and/or emission geometries. We implemented an offset-dipole magnetic field, with an offset characterised by parameters epsilon and magnetic azimuthal angle, in an existing geometric pulsar modelling code which already includes static and retarded vacuum dipole fields. We use these different magnetic field solutions in conjunction with standard emission geometries, namely the two-pole caustic and outer gap models (the latter only for non-offset dipoles), and construct intensity maps and light curves for several pulsar parameters. We compare our model light curves to the Vela data from the second pulsar catalogue of Fermi. We use a refined chi-square grid search method for finding best-fit light curves for each of the different models. Our best fit is for the retarded vacuum dipole field and the outer gap model.

preprint2015arXiv

The Effect of Different Magnetospheric Structures on Predictions of Gamma-ray Pulsar Light Curves

The second pulsar catalogue of the Fermi Large Area Telescope (LAT) will contain in excess of 100 gamma-ray pulsars. The light curves (LCs) of these pulsars exhibit a variety of shapes, and also different relative phase lags with respect to their radio pulses, hinting at distinct underlying emission properties (e.g., inclination and observer angles) for the individual pulsars. Detailed geometric modelling of the radio and gamma-ray LCs may provide constraints on the B-field structure and emission geometry. We used different B-field solutions, including the static vacuum dipole and the retarded vacuum dipole, in conjunction with an existing geometric modelling code, and constructed radiation sky maps and LCs for several different pulsar parameters. Standard emission geometries were assumed, namely the two-pole caustic (TPC) and outer gap (OG) models. The sky maps and LCs of the various B-field and radiation model combinations were compared to study their effect on the resulting LCs. As an application, we compared our model LCs with Fermi LAT data for the Vela pulsar, and inferred the most probable configuration in this case, thereby constraining Vela's high-altitude magnetic structure and system geometry.

preprint2014arXiv

Constraints on the Emission Geometries and Spin Evolution of Gamma-ray Millisecond Pulsars

Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using a maximum likelihood technique. We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II) or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.

preprint2014arXiv

Implementation of an offset-dipole magnetic field in a pulsar modelling code

The light curves of gamma-ray pulsars detected by the Fermi Large Area Telescope show great variety in profile shape and position relative to their radio profiles. Such diversity hints at distinct underlying magnetospheric and/or emission geometries for the individual pulsars. We implemented an offset-dipole magnetic field in an existing geometric pulsar modelling code which already includes static and retarded vacuum dipole fields. In our model, this offset is characterised by a parameter epsilon (with epsilon = 0 corresponding to the static dipole case). We constructed sky maps and light curves for several pulsar parameters and magnetic fields, studying the effect of an offset dipole on the resulting light curves. A standard two-pole caustic emission geometry was used. As an application, we compared our model light curves with Fermi data for the bright Vela pulsar.

preprint2014arXiv

Inferred cosmic-ray spectrum from ${\it Fermi}$-LAT $γ$-ray observations of the Earth's limb

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $γ$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

preprint2014arXiv

Light-curve modelling constraints on the obliquities and aspect angles of the young Fermi pulsars

In more than four years of observation the Large Area Telescope on board the Fermi satellite has identified pulsed $γ$-ray emission from more than 80 young pulsars, providing light curves with high statistics. Fitting the observations with geometrical models can provide estimates of the magnetic obliquity $α$ and aspect angle $ζ$, yielding estimates of the radiation beaming factor and luminosity. Using $γ$-ray emission geometries (Polar Cap, Slot Gap, Outer Gap, One Pole Caustic) and radio emission geometry, we fit $γ$-ray light curves for 76 young pulsars and we jointly fit their $γ$-ray plus radio light curves when possible. We find that a joint radio plus $γ$-ray fit strategy is important to obtain ($α$, $ζ$) estimates that can explain simultaneous radio and $γ$-ray emission. The intermediate-to-high altitude magnetosphere models, Slot Gap, Outer Gap, and One pole Caustic, are favoured in explaining the observations. We find no evolution of $α$ on a time scale of a million years. For all emission geometries our derived $γ$-ray beaming factors are generally less than one and do not significantly evolve with the spin-down power. A more pronounced beaming factor vs. spin-down power correlation is observed for Slot Gap model and radio-quiet pulsars and for the Outer Gap model and radio-loud pulsars. For all models, the correlation between $γ$-ray luminosity and spin-down power is consistent with a square root dependence. The $γ$-ray luminosities obtained by using our beaming factors not exceed the spin-down power. This suggests that assuming a beaming factor of one for all objects, as done in other studies, likely overestimates the real values. The data show a relation between the pulsar spectral characteristics and the width of the accelerator gap that is consistent with the theoretical prediction for the Slot Gap model.

preprint2014arXiv

Modelling the $γ$-ray and radio light curves of the double pulsar system

Guillemot et al. recently reported the discovery of $γ$-ray pulsations from the 22.7ms pulsar (pulsar A) in the famous double pulsar system J0737-3039A/B. The $γ$-ray light curve (LC) of pulsar A has two peaks separated by approximately half a rotation, and these are non-coincident with the observed radio and X-ray peaks. This suggests that the $γ$-ray emission originates in a part of the magnetosphere distinct from where the radio and X-ray radiation is generated. Thus far, three different methods have been applied to constrain the viewing geometry of pulsar A (its inclination and observer angles $α$ and $ζ$): geometric modelling of the radio and $γ$-ray light curves, modelling of the position angle sweep in phase seen in the radio polarisation data, and independent studies of the time evolution of the radio pulse profile of pulsar A. These three independent, complementary methods have yielded consistent results: pulsar A's rotation axis is likely perpendicular to the orbital plane of the binary system, and its magnetic axis close to lying in the orbital plane (making this pulsar an orthogonal rotator). The observer is furthermore observing emission close to the magnetic axis. Thus far, however, current models could not reproduce all the characteristics of the radio and $γ$-ray light curves, specifically the large radio-to-$γ$ phase lag. In this paper we discuss some preliminary modelling attempts to address this problem, and offer ideas on how the LC fits may be improved by adapting the standard geometric models in order to reproduce the profile positions more accurately.

preprint2013arXiv

Broadband Pulsations from PSR B1821-24: Implications for Emission Models and the Pulsar Population of M28

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

preprint2013arXiv

Constraints on the Galactic Population of TEV Pulsar Wind Nebulae Using Fermi Large Area Telescope Observations

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV gamma-ray emitters. Since launch, the Fermi Large Area Telescope (LAT)identified five high-energy (100MeV <E< 100 GeV) gamma-ray sources as PWNe, and detected a large number of PWNe candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV gamma-ray unidentifiedsources (UNIDs) are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58TeV PWNe and UNIDs within 5deg of the Galactic Plane to establish new constraints on PWNe properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their gamma-rayfluxes for energies above 10 GeV. The spectral energy distributions (SED) andupper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e. between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWNe candidates are described in detail and compared with existing models. A population study of GeV PWNe candidates as a function of the pulsar/PWN system characteristics is presented.

preprint2013arXiv

Detection of the Characteristic Pion-Decay Signature in Supernova Remnants

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

preprint2013arXiv

Fermi LAT pulsed detection of PSR J0737-3039A in the double pulsar system

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

preprint2013arXiv

GMRT discovery of PSR J1544+4937, an eclipsing black-widow pulsar identified with a Fermi LAT source

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

preprint2013arXiv

Modeling the non-recycled Fermi gamma-ray pulsar population

We use Fermi Gamma-ray Space Telescope detections and upper limits on non-recycled pulsars obtained from the Large Area Telescope (LAT) to constrain how the gamma-ray luminosity L depends on the period P and the period derivative \dot{P}. We use a Bayesian analysis to calculate a best-fit luminosity law, or dependence of L on P and \dot{P}, including different methods for modeling the beaming factor. An outer gap (OG) magnetosphere geometry provides the best-fit model, which is L \propto P^{-a} \dot{P}^{b} where a=1.36\pm0.03 and b=0.44\pm0.02, similar to but not identical to the commonly assumed L \propto \sqrt{\dot{E}} \propto P^{-1.5} \dot{P}^{0.5}. Given upper limits on gamma-ray fluxes of currently known radio pulsars and using the OG model, we find that about 92% of the radio-detected pulsars have gamma-ray beams that intersect our line of sight. By modeling the misalignment of radio and gamma-ray beams of these pulsars, we find an average gamma-ray beaming solid angle of about 3.7π for the OG model, assuming a uniform beam. Using LAT-measured diffuse fluxes, we place a 2σ upper limit on the average braking index and a 2σ lower limit on the average surface magnetic field strength of the pulsar population of 3.8 and 3.2 X 10^{10} G, respectively. We then predict the number of non-recycled pulsars detectable by the LAT based on our population model. Using the two-year sensitivity, we find that the LAT is capable of detecting emission from about 380 non-recycled pulsars, including 150 currently identified radio pulsars. Using the expected five-year sensitivity, about 620 non-recycled pulsars are detectable, including about 220 currently identified radio pulsars. We note that these predictions significantly depend on our model assumptions.

preprint2013arXiv

The First Pulse of the Extremely Bright GRB 130427A: A Test Lab for Synchrotron Shocks

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

preprint2013arXiv

The Vela-X Pulsar Wind Nebula revisited with 4 years of Fermi-Large Area Telescope observations

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

preprint2012arXiv

Binary Millisecond Pulsar Discovery via Gamma-Ray Pulsations

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

preprint2012arXiv

Constraining gamma-ray pulsar gap models with a simulated pulsar population

With the large sample of young gamma-ray pulsars discovered by the Fermi Large Area Telescope (LAT), population synthesis has become a powerful tool for comparing their collective properties with model predictions. We synthesised a pulsar population based on a radio emission model and four gamma-ray gap models (Polar Cap, Slot Gap, Outer Gap, and One Pole Caustic) normalizing to the number of detected radio pulsars in select group of surveys. The luminosity and the wide beams from the outer gaps can easily account for the number of Fermi detections in 2 years of observations. The wide slot-gap beams requires an increase by a factor of ~10 of the predicted luminosity to produce a reasonable number of gamma-ray pulsars. Such large increases in the luminosity may be accommodated by implementing offset polar caps. The narrow polar-cap beams contribute at most only a handful of LAT pulsars. Standard distributions in birth location and pulsar spin-down power (Edot) fail to reproduce the LAT findings: all models under-predict the number of LAT pulsars with high Edot, and they cannot explain the high probability of detecting both the radio and gamma-ray beams at high Edot. The beaming factor remains close to 1 over 4 decades in Edot evolution for the slot gap whereas it significantly decreases with increasing age for the outer gaps. The evolution of the slot-gap luminosity with Edot is compatible with the large dispersion of gamma-ray luminosity seen in the LAT data. The stronger evolution predicted for the outer gap, which is linked to the polar cap heating by the return current, is apparently not supported by the LAT data. The LAT sample of gamma-ray pulsars therefore provides a fresh perspective on the early evolution of the luminosity and beam width of the gamma-ray emission from young pulsars, calling for thin and more luminous gaps.

preprint2012arXiv

Constraining viewing geometries of pulsars with single-peaked gamma-ray profiles using a multiwavelength approach

Since the launch of the Large Area Telescope (LAT) on board the Fermi spacecraft in June 2008, the number of observed gamma-ray pulsars has increased dramatically. A large number of these are also observed at radio frequencies. Constraints on the viewing geometries of 5 of 6 gamma-ray pulsars exhibiting single-peaked gamma-ray profiles were derived using high-quality radio polarization data (Weltevrede et al., 2010). We obtain independent constraints on the viewing geometries of 6 by using a geometric emission code to model the Fermi LAT and radio light curves (LCs). We find fits for the magnetic inclination and observer angles by searching the solution space by eye. Our results are generally consistent with those previously obtained (Weltevrede et al., 2010), although we do find small differences in some cases. We will indicate how the gamma-ray and radio pulse shapes as well as their relative phase lags lead to constraints in the solution space. Values for the flux correction factor corresponding to the fits are also derived (with errors).

preprint2012arXiv

Discovery of the millisecond pulsar PSR J2043+1711 in a Fermi source with the Nancay Radio Telescope

We report the discovery of the millisecond pulsar PSR J2043+1711 in a search of a Fermi Large Area Telescope (LAT) source with no known associations, with the Nancay Radio Telescope. The new pulsar, confirmed with the Green Bank Telescope, has a spin period of 2.38 ms, is relatively nearby (d <~ 2 kpc), and is in a 1.48 day orbit around a low mass companion, probably a He-type white dwarf. Pulsed gamma-ray emission was detected in the data recorded by the Fermi LAT. The gamma-ray light curve and spectral properties are typical of other gamma-ray millisecond pulsars seen with Fermi. X-ray observations of the pulsar with Suzaku and the Swift/XRT yielded no detection. At 1.4 GHz we observe strong flux density variations because of interstellar diffractive scintillation, however a sharp peak can be observed at this frequency during bright scintillation states. At 327 MHz the pulsar is detected with a much higher signal-to-noise ratio and its flux density is far more steady. However, at that frequency the Arecibo instrumentation cannot yet fully resolve the pulse profile. Despite that, our pulse time-of-arrival measurements have a post-fit residual rms of 2 \mus. This and the expected stability of this system has made PSR J2043+1711 one of the first new Fermi-selected millisecond pulsars to be added to pulsar gravitational wave timing arrays. It has also allowed a significant measurement of relativistic delays in the times of arrival of the pulses due to the curvature of space-time near the companion, but not yet with enough precision to derive useful masses for the pulsar and the companion. A mass for the pulsar between 1.7 and 2.0 solar masses can be derived if a standard millisecond pulsar formation model is assumed. In this article we also present a comprehensive summary of pulsar searches in Fermi LAT sources with the Nancay Radio Telescope to date.

preprint2012arXiv

Five New Millisecond Pulsars From a Radio Survey of 14 Unidentified Fermi-LAT Gamma-ray Sources

We have discovered five millisecond pulsars (MSPs) in a survey of 14 unidentified Fermi-LAT sources in the southern sky using the Parkes radio telescope. PSRs J0101-6422, J1514-4946, and J1902-5105 reside in binaries, while PSRs J1658-5324 and J1747-4036 are isolated. Using an ephemeris derived from timing observations of PSR J0101-6422 (P =2.57 ms, DM=12 pc cm-3), we have detected γ-ray pulsations and measured its proper motion. Its γ-ray spectrum (a power law of Γ = 0.9 with a cutoff at 1.6 GeV) and efficiency are typical of other MSPs, but its radio and γ-ray light curves challenge simple geometric models of emission. The high success rate of this survey-enabled by selecting γ-ray sources based on their detailed spectral characteristics-and other similarly successful searches indicate that a substantial fraction of the local population of MSPs may soon be known.

preprint2012arXiv

Gamma-ray Activity in the Crab Nebula: The Exceptional Flare of April 2011

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

preprint2012arXiv

Gamma-ray observations of the Orion Molecular Clouds with the Fermi Large Area Telescope

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1° \times1° pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212°), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

preprint2012arXiv

Measurement of separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope

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

preprint2012arXiv

Radio Searches of Fermi LAT Sources and Blind Search Pulsars: The Fermi Pulsar Search Consortium

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

preprint2012arXiv

SGR J1550-5418 bursts detected with the Fermi Gamma-ray Burst Monitor during its most prolific activity

We have performed detailed temporal and time-integrated spectral analysis of 286 bursts from SGR J1550-5418 detected with the Fermi Gamma-ray Burst Monitor (GBM) in January 2009, resulting in the largest uniform sample of temporal and spectral properties of SGR J1550-5418 bursts. We have used the combination of broadband and high time-resolution data provided with GBM to perform statistical studies for the source properties. We determine the durations, emission times, duty cycles and rise times for all bursts, and find that they are typical of SGR bursts. We explore various models in our spectral analysis, and conclude that the spectra of SGR J1550-5418 bursts in the 8-200 keV band are equally well described by optically thin thermal bremsstrahlung (OTTB), a power law with an exponential cutoff (Comptonized model), and two black-body functions (BB+BB). In the spectral fits with the Comptonized model we find a mean power-law index of -0.92, close to the OTTB index of -1. We show that there is an anti-correlation between the Comptonized Epeak and the burst fluence and average flux. For the BB+BB fits we find that the fluences and emission areas of the two blackbody functions are correlated. The low-temperature BB has an emission area comparable to the neutron star surface area, independent of the temperature, while the high-temperature blackbody has a much smaller area and shows an anti-correlation between emission area and temperature. We compare the properties of these bursts with bursts observed from other SGR sources during extreme activations, and discuss the implications of our results in the context of magnetar burst models.

preprint2011arXiv

Discovery of millisecond pulsars in radio searches of southern Fermi LAT sources

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

preprint2011arXiv

Discovery of two millisecond pulsars in Fermi sources with the Nancay Radio Telescope

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

preprint2011arXiv

Geometric modelling of radio and gamma-ray light curves of 6 Fermi LAT pulsars

The Fermi Large Area Telescope (LAT) has recently reported the detection of pulsed gamma-rays from 6 young pulsars (J0631+1036, J0659+1414, J0742-2822, J1420-6048, J1509-5850, and J1718-3825), all exhibiting single-peaked pulse profiles (Weltevrede et al., 2010). High-quality radio polarization data are also available for 5 of these pulsars, allowing derivation of constraints on their viewing geometries. We obtain independent constraints on the viewing geometries of these pulsars by using a geometric pulsar emission code to model the Fermi LAT and radio light curves. We find fits for the magnetic inclination and observer angles alpha and zeta with typical errors of ~ 5deg. Our results are generally consistent with those obtained by Weltevrede et al. (2010), although we do find differences in some cases. Our model may lastly provide a framework to constrain the radio emission altitude.

preprint2011arXiv

Modeling Light Curves of the Phase-Aligned Gamma-ray Millisecond Pulsar Subclass

The gamma-ray population of millisecond pulsars (MSPs) detected by the Fermi Large Area Telescope (LAT) has been steadily increasing. A number of the more recent detections, including PSR J0034-0534, PSR J1939+2134 (B1937+21; the first MSP ever discovered), PSR J1959+2048 (B1957+20; the first black widow system), and PSR J2214+3000 exhibit an unusual phenomenon: nearly phase-aligned radio and gamma-ray light curves (LCs). To account for the phase alignment, we explore geometric models where both the radio and gamma-ray emission originate either in the outer magnetosphere near the light cylinder (R_LC) or near the polar caps (PCs). We obtain reasonable fits for the first three of these MSPs in the context of "altitude-limited" outer gap (alOG) and two-pole caustic (alTPC) geometries. The outer magnetosphere phase-aligned models differ from the standard outer gap (OG) / two-pole caustic (TPC) models in two respects: first, the radio emission originates in caustics at relatively high altitudes compared to the usual low-altitude conal radio beams; second, we allow the maximum altitude of the gamma-ray emission region as well as both the minimum and maximum altitudes of the radio emission region to vary within a limited range. Alternatively, there also exist phase-aligned LC solutions for emission originating near the stellar surface in a slot gap (SG) scenario ("low-altitude slot gap" (laSG) models). We find best-fit LCs using a Markov chain Monte Carlo (MCMC) maximum likelihood approach. Our fits imply that the phase-aligned LCs are likely of caustic origin, produced in the outer magnetosphere, and that the radio emission may come from close to R_LC. We lastly constrain the emission altitudes with typical uncertainties of ~0.3R_LC. Our results describe a third gamma-ray MSP subclass, in addition to the two (with non-aligned LCs) previously found.

preprint2011arXiv

Observations of Energetic High Magnetic Field Pulsars with the Fermi Large Area Telescope

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

preprint2011arXiv

Precise Gamma-Ray Timing and Radio Observations of 17 Fermi Gamma-Ray Pulsars

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

preprint2011arXiv

PSR J0007+7303 in the CTA1 SNR: New Gamma-ray Results from Two Years of Fermi-LAT Observations

One of the main results of the Fermi Gamma-Ray Space Telescope is the discovery of γ-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its γ-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of γ-ray emission in the off-pulse phase interval at the ~ 6σ level. The flux from this emission in the energy range E \geq 100 MeV is F_100 = (1.73\pm0.40)\times10^(-8) photons/cm^2/s and is best fitted by a power law with a photon index of Γ = 2.54\pm0.14. The pulsed γ-ray flux in the same energy range is F_100 = (3.95\pm0.07)\times10^(-7) photons/cm^2/s and is best fitted by an exponentially-cutoff power-law spectrum with a photon index of Γ = 1.41 \pm 0.23 and a cutoff energy E_c = 4.04 \pm 0.20 GeV. We find no flux variability neither at the 2009 May glitch nor in the long term behavior. We model the γ-ray light curve with two high-altitude emission models, the outer gap and slot gap, and find that the model that best fits the data depends strongly on the assumed origin of the off-pulse emission. Both models favor a large angle between the magnetic axis and observer line of sight, consistent with the nondetection of radio emission being a geometrical effect. Finally we discuss how the LAT results bear on the understanding of the cooling of this neutron star.

preprint2011arXiv

PSR J2030+3641: radio discovery and gamma-ray study of a middle-aged pulsar in the now identified Fermi-LAT source 1FGL J2030.0+3641

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

preprint2011arXiv

Pulsed Gamma Rays from the Original Millisecond and Black Widow Pulsars: a case for Caustic Radio Emission?

We report the detection of pulsed gamma-ray emission from the fast millisecond pulsars (MSPs) B1937+21 (also known as J1939+2134) and B1957+20 (J1959+2048) using 18 months of survey data recorded by the \emph{Fermi} Large Area Telescope (LAT) and timing solutions based on radio observations conducted at the Westerbork and Nançay radio telescopes. In addition, we analyzed archival \emph{RXTE} and \emph{XMM-Newton} X-ray data for the two MSPs, confirming the X-ray emission properties of PSR B1937+21 and finding evidence ($\sim 4σ$) for pulsed emission from PSR B1957+20 for the first time. In both cases the gamma-ray emission profile is characterized by two peaks separated by half a rotation and are in close alignment with components observed in radio and X-rays. These two pulsars join PSRs J0034-0534 and J2214+3000 to form an emerging class of gamma-ray MSPs with phase-aligned peaks in different energy bands. The modeling of the radio and gamma-ray emission profiles suggests co-located emission regions in the outer magnetosphere.

preprint2011arXiv

What can we learn from phase alignment of gamma-ray and radio pulsar light curves?

The Fermi Large Area Telescope (LAT) has revolutionized high-energy (HE) astronomy, and is making enormous contributions particularly to gamma-ray pulsar science. As a result of the many new pulsar discoveries, the gamma-ray pulsar population is now approaching 100. Some very famous millisecond pulsars (MSPs) have also been detected: J1939+2134 (B1937+21), the first MSP ever discovered, as well as J1959+2048 (B1957+20), the first black widow pulsar system. These, along with other MSPs such as PSR J0034-0534 and J2214+3000, are rare among the pulsar population in that they exhibit nearly phase-aligned radio and gamma-ray light curves (LCs). Traditionally, pulsar LCs have been modelled using standard HE models in conjunction with low-altitude conal beam radio models. However, a different approach is needed to account for phase-aligned LCs. We explored two scenarios: one where both the radio and gamma-ray emission originate in the outer magnetosphere, and one where the emission comes from near the polar caps (PCs) on the stellar surface. We find best-fit LCs using a Markov chain Monte Carlo (MCMC) technique for the first class of models. The first scenario seems to be somewhat preferred, as is also hinted at by the radio polarization data. This implies that the phase-aligned LCs are possibly of caustic origin produced in the outer magnetosphere, in contrast to the usual lower-altitude conal beam radio models. Lastly, we constrain the emission altitudes with typical uncertainties of ~10% of the light cylinder radius. The modelled pulsars are members of a third gamma-ray MSP subclass, in addition to two others with non-aligned radio and gamma-ray LCs.

preprint2010arXiv

Constraints on Cosmological Dark Matter Annihilation from the Fermi-LAT Isotropic Diffuse Gamma-Ray Measurement

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

preprint2010arXiv

Eight gamma-ray pulsars discovered in blind frequency searches of Fermi LAT data

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

preprint2010arXiv

Fermi LAT Search for Photon Lines from 30 to 200 GeV and Dark Matter Implications

Dark matter (DM) particle annihilation or decay can produce monochromatic $γ$-rays readily distinguishable from astrophysical sources. $γ$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $γ$-ray line analysis, and integrated over most of the sky. We obtain $γ$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

preprint2010arXiv

PSRs J0248+6021 and J2240+5832: Young Pulsars in the Northern Galactic Plane. Discovery, Timing, and Gamma-ray observations

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

preprint2010arXiv

Radio and Gamma-Ray Constraints on the Emission Geometry and Birthplace of PSR J2043+2740

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

preprint2010arXiv

Three Millisecond Pulsars in FERMI LAT Unassociated Bright Sources

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

preprint2009arXiv

Probing millisecond pulsar emission geometry using light curves from the Fermi Large Area Telescope

An interesting new high-energy pulsar sub-population is emerging following early discoveries of gamma-ray millisecond pulsars (MSPs) by the Fermi Large Area Telescope (LAT). We present results from 3D emission modeling, including the Special Relativistic effects of aberration and time-of-flight delays and also rotational sweepback of B-field lines, in the geometric context of polar cap (PC), outer gap (OG), and two-pole caustic (TPC) pulsar models. In contrast to the general belief that these very old, rapidly-rotating neutron stars (NSs) should have largely pair-starved magnetospheres due to the absence of significant pair production, we find that most of the light curves are best fit by TPC and OG models, which indicates the presence of narrow accelerating gaps limited by robust pair production -- even in these pulsars with very low spin-down luminosities. The gamma-ray pulse shapes and relative phase lags with respect to the radio pulses point to high-altitude emission being dominant for all geometries. We also find exclusive differentiation of the current gamma-ray MSP population into two MSP sub-classes: light curve shapes and lags across wavebands impose either pair-starved PC (PSPC) or TPC / OG-type geometries. In the first case, the radio pulse has a small lag with respect to the single gamma-ray pulse, while the (first) gamma-ray peak usually trails the radio by a large phase offset in the latter case. Finally, we find that the flux correction factor as a function of magnetic inclination and observer angles is typically of order unity for all models. Our calculation of light curves and flux correction factor for the case of MSPs is therefore complementary to the "ATLAS paper" of Watters et al. for younger pulsars.

preprint2009arXiv

Radio detection of LAT PSRs J1741-2054 and J2032+4127: no longer just gamma-ray pulsars

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

preprint2009arXiv

The Spectral Energy Distribution of Fermi bright blazars

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the γ-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $ν$ - Log $ν$ F$_ν$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $ν_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

A. K. Harding

What is connected

Connect this record

See the researcher in context

Building this map preview

73 published item(s)

A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Pulse Peak Migration during the Outburst Decay of the Magnetar SGR 1830-0645: Crustal Motion and Magnetospheric Untwisting

Search for new cosmic-ray acceleration sites within the 4FGL catalog Galactic plane sources

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

X-ray burst and persistent emission properties of the magnetar SGR 1830-0645 in outburst

Broadband X-ray Burst Spectroscopy of the FRB-Emitting Galactic Magnetar

Einstein@Home Discovery of the Gamma-ray Millisecond Pulsar PSR J2039-5617 Confirms Its Predicted Redback Nature

A radiatively-quiet glitch and anti-glitch in the magnetar 1E 2259+586

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

3FGL Demographics Outside the Galactic Plane using Supervised Machine Learning: Pulsar and Dark Matter Subhalo Interpretations

A New, Low Braking Index For the LMC Pulsar B0540-69

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Implementation of a goodness-of-fit test for finding optimal concurrent radio and γ-ray pulsar light curves

Localization and broadband follow-up of the gravitational-wave transient GW150914

Measurement of the high-energy gamma-ray emission from the Moon with the Fermi Large Area Telescope

Pulsar emission in the very-high-energy regime

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

Very-high energy emission from pulsars

Young and middle age pulsar light-curve morphology: Comparison of Fermi observations with gamma-ray and radio emission geometries

Cosmic-ray positrons from millisecond pulsars

Discovery of a Spin-Down State Change in the LMC Pulsar B0540-69

Discovery of Gamma-ray Pulsations from the Transitional Redback PSR J1227-4853

Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113

PSR J1906+0722: An Elusive Gamma-ray Pulsar

Search for Early Gamma-ray Production in Supernovae Located in a Dense Circumstellar Medium with the Fermi LAT

The 1st Fermi Lat Supernova Remnant Catalog

The anatomy of $γ$-ray pulsar light curves

The Contribution of Millisecond Pulsars to the Galactic Cosmic-Ray Lepton Spectrum

The effect of an offset-dipole magnetic field on the Vela pulsar's gamma-ray light curves

The Effect of Different Magnetospheric Structures on Predictions of Gamma-ray Pulsar Light Curves

Constraints on the Emission Geometries and Spin Evolution of Gamma-ray Millisecond Pulsars

Implementation of an offset-dipole magnetic field in a pulsar modelling code

Inferred cosmic-ray spectrum from ${\it Fermi}$-LAT $γ$-ray observations of the Earth's limb

Light-curve modelling constraints on the obliquities and aspect angles of the young Fermi pulsars

Modelling the $γ$-ray and radio light curves of the double pulsar system

Broadband Pulsations from PSR B1821-24: Implications for Emission Models and the Pulsar Population of M28

Constraints on the Galactic Population of TEV Pulsar Wind Nebulae Using Fermi Large Area Telescope Observations

Detection of the Characteristic Pion-Decay Signature in Supernova Remnants

Fermi LAT pulsed detection of PSR J0737-3039A in the double pulsar system

GMRT discovery of PSR J1544+4937, an eclipsing black-widow pulsar identified with a Fermi LAT source

Modeling the non-recycled Fermi gamma-ray pulsar population

The First Pulse of the Extremely Bright GRB 130427A: A Test Lab for Synchrotron Shocks

The Vela-X Pulsar Wind Nebula revisited with 4 years of Fermi-Large Area Telescope observations

Binary Millisecond Pulsar Discovery via Gamma-Ray Pulsations

Constraining gamma-ray pulsar gap models with a simulated pulsar population

Constraining viewing geometries of pulsars with single-peaked gamma-ray profiles using a multiwavelength approach

Discovery of the millisecond pulsar PSR J2043+1711 in a Fermi source with the Nancay Radio Telescope

Five New Millisecond Pulsars From a Radio Survey of 14 Unidentified Fermi-LAT Gamma-ray Sources

Gamma-ray Activity in the Crab Nebula: The Exceptional Flare of April 2011

Gamma-ray observations of the Orion Molecular Clouds with the Fermi Large Area Telescope

Measurement of separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope

Radio Searches of Fermi LAT Sources and Blind Search Pulsars: The Fermi Pulsar Search Consortium

SGR J1550-5418 bursts detected with the Fermi Gamma-ray Burst Monitor during its most prolific activity

Discovery of millisecond pulsars in radio searches of southern Fermi LAT sources

Discovery of two millisecond pulsars in Fermi sources with the Nancay Radio Telescope

Geometric modelling of radio and gamma-ray light curves of 6 Fermi LAT pulsars

Modeling Light Curves of the Phase-Aligned Gamma-ray Millisecond Pulsar Subclass

Observations of Energetic High Magnetic Field Pulsars with the Fermi Large Area Telescope

Precise Gamma-Ray Timing and Radio Observations of 17 Fermi Gamma-Ray Pulsars

PSR J0007+7303 in the CTA1 SNR: New Gamma-ray Results from Two Years of Fermi-LAT Observations

PSR J2030+3641: radio discovery and gamma-ray study of a middle-aged pulsar in the now identified Fermi-LAT source 1FGL J2030.0+3641

Pulsed Gamma Rays from the Original Millisecond and Black Widow Pulsars: a case for Caustic Radio Emission?

What can we learn from phase alignment of gamma-ray and radio pulsar light curves?

Constraints on Cosmological Dark Matter Annihilation from the Fermi-LAT Isotropic Diffuse Gamma-Ray Measurement

Eight gamma-ray pulsars discovered in blind frequency searches of Fermi LAT data

Fermi LAT Search for Photon Lines from 30 to 200 GeV and Dark Matter Implications

PSRs J0248+6021 and J2240+5832: Young Pulsars in the Northern Galactic Plane. Discovery, Timing, and Gamma-ray observations

Radio and Gamma-Ray Constraints on the Emission Geometry and Birthplace of PSR J2043+2740

Three Millisecond Pulsars in FERMI LAT Unassociated Bright Sources

Probing millisecond pulsar emission geometry using light curves from the Fermi Large Area Telescope

Radio detection of LAT PSRs J1741-2054 and J2032+4127: no longer just gamma-ray pulsars

The Spectral Energy Distribution of Fermi bright blazars