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.Comment: 51 pages, 23 figures, 3 tables; low-resolution figures; accepted for publication by Ap

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.Comment: 10 pages, 5 figures, "To appear in the Proceedings of the 3rd Annual Conference on High Energy Astrophysics in Southern Africa (HEASA2015), 18-20 June 2015, Johannesburg, South Africa, Eds. M. Boettcher, D. Buckley, S. Colafrancesco, P. Meintjes and S. Razzaque.

Predictions of Gamma-ray Emission from Globular Cluster Millisecond Pulsars Above 100 MeV

The recent Fermi detection of the globular cluster (GC) 47 Tucanae highlighted the importance of modeling collective gamma-ray emission of millisecond pulsars (MSPs) in GCs. Steady flux from such populations is also expected in the very high energy (VHE) domain covered by ground-based Cherenkov telescopes. We present pulsed curvature radiation (CR) as well as unpulsed inverse Compton (IC) calculations for an ensemble of MSPs in the GCs 47 Tucanae and Terzan 5. We demonstrate that the CR from these GCs should be easily detectable for Fermi, while constraints on the total number of MSPs and the nebular B-field may be derived using the IC flux components.Comment: 12 pages, 2 figures, accepted for publication in ApJ