Discovery and timing of three millisecond pulsars in radio and gamma-rays with the GMRT and Fermi-LAT

We performed deep observations to search for radio pulsations in the directions of 375 unassociated Fermi Large Area Telescope (LAT) gamma-ray sources using the Giant Metrewave Radio Telescope (GMRT) at 322 and 607 MHz. In this paper we report the discovery of three millisecond pulsars (MSPs), PSR J0248+4230, PSR J1207$-$5050 and PSR J1536$-$4948. We conducted follow up timing observations for around 5 years with the GMRT and derived phase coherent timing models for these MSPs. PSR J0248$+$4230 and J1207$-$5050 are isolated MSPs having periodicities of 2.60 ms and 4.84 ms. PSR J1536-4948 is a 3.07 ms pulsar in a binary system with orbital period of around 62 days about a companion of minimum mass 0.32 solar mass. We also present multi-frequency pulse profiles of these MSPs from the GMRT observations. PSR J1536-4948 is an MSP with an extremely wide pulse profile having multiple components. Using the radio timing ephemeris we subsequently detected gamma-ray pulsations from these three MSPs, confirming them as the sources powering the gamma-ray emission. For PSR J1536-4948 we performed combined radio-gamma-ray timing using around 11.6 years of gamma-ray pulse times of arrivals (TOAs) along with the radio TOAs. PSR J1536-4948 also shows evidence for pulsed gamma-ray emission out to above 25 GeV, confirming earlier associations of this MSP with a >10 GeV point source. The multi-wavelength pulse profiles of all three MSPs offer challenges to models of radio and gamma-ray emission in pulsar magnetospheres.

A Luminous and Highly-variable Gamma-ray Flare Following the 2017 Periastron of PSR B1259-63/LS 2883

Three periastron passages of the PSR B1259$-$63/LS 2883 binary system, consisting of a 48 ms rotation-powered pulsar and a $sim30$ M$_{odot}$ Be star, have been observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, in 2010, 2014, and 2017. During the most-recent periastron passage, sustained low-level gamma-ray emission was observed over a $sim3$-week long interval immediately after periastron, which was followed by an interval of no emission. Sporadic flares were detected starting 40 days post-periastron and lasted approximately 50 days, during which the emission displayed significant spectral curvature, variability on timescales as short as 1.5 minutes, and peak flux levels well in excess of the pulsar spin-down power. By contrast, during the 2010 and 2014 periastron passages, significant gamma-ray emission was not observed with the LAT until 30 and 32 days post-periastron, respectively. The previous flares did not exhibit spectral curvature, showed no short term variability, and did not exceed the pulsar spin-down power. The high flux and short timescales observed in 2017 suggest significant beaming of the emission is required and constrain the size of the emission region. The flares occur long enough after periastron that the neutron star should already have passed through the extended disk-like outflow, thus constraining options for target material and seed photon sources for inverse Compton models.

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

The anatomy of $gamma$-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 $gamma$-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 $gamma$-ray model LCs proved difficult though due to the increased complexity of the geometric models for the $gamma$-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.

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.

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

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

Fermi 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 $alpha$ and the viewing angle $zeta$ to be close to 90$^circ$, which is consistent with independent studies of the radio emission from PSR J0737-3039A. A small misalignment angle between the pulsars spin axis and the systems 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.

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

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 Nanc{c}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 4sigma$) 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.

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