Emission In Pulsar Magnetosphere Research Articles

Discovery and Timing of Three Millisecond Pulsars in Radio and Gamma-Rays with the Giant Metrewave Radio Telescope and Fermi Large Area Telescope

Abstract We performed deep observations to search for radio pulsations in the directions of 375 unassociated Fermi Large Area Telescope γ-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 ∼5 yr 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 an orbital period of ∼62 days about a companion of a minimum mass of 0.32 M ⊙. We also present multifrequency 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 γ-ray pulsations from these three MSPs, confirming them as the sources powering the γ-ray emission. For PSR J1536–4948 we performed combined radio–γ-ray timing using ∼11.6 yr of γ-ray pulse times of arrival (TOAs) along with the radio TOAs. PSR J1536–4948 also shows evidence for pulsed γ-ray emission out to above 25 GeV, confirming earlier associations of this MSP with a ≥10 GeV point source. The multiwavelength pulse profiles of all three MSPs offer challenges to models of radio and γ-ray emission in pulsar magnetospheres.

General-relativistic pulsar magnetospheric emission

Most current pulsar emission models assume photon production and emission within the magnetosphere. Low frequency radiation is preferentially produced in the vicinity of the polar caps whereas the high-energy tail is shifted to regions closer but still inside the light-cylinder. We conducted a systematic study of the merit of several popular radiation sites like the polar cap, the outer gap and the slot gap. We computed sky maps emanating from each emission site according to a prescribed distribution function for the emitting particles made of an electron/positron mixture. Calculations are performed using a three dimensional integration of the plasma emissivity in the vacuum electromagnetic field of a rotating centred general-relativistic dipole. We compare Newtonian electromagnetic fields to their general-relativistic counterpart. In the latter case, light bending is also taken into account. As a typical example, light-curves and sky maps are plotted for several power-law indices of the particle distribution function. The detailed pulse profiles strongly depend on the underlying assumption about the fluid motion subject to strong electromagnetic fields. This electromagnetic topology enforces the photon propagation direction directly, or indirectly, from aberration effects. We also discuss the implication of a net stellar electric charge on to the sky maps. Taking into account the electric field strongly affects the light-curves originating close to the light-cylinder where the electric field strength becomes comparable to the magnetic field strength.

The Einstein@Home Gamma-ray Pulsar Survey. II. Source Selection, Spectral Analysis, and Multiwavelength Follow-up

Abstract We report on the analysis of 13 gamma-ray pulsars discovered in the Einstein@Home blind search survey using Fermi Large Area Telescope (LAT) Pass 8 data. The 13 new gamma-ray pulsars were discovered by searching 118 unassociated LAT sources from the third LAT source catalog (3FGL), selected using the Gaussian Mixture Model machine-learning algorithm on the basis of their gamma-ray emission properties being suggestive of pulsar magnetospheric emission. The new gamma-ray pulsars have pulse profiles and spectral properties similar to those of previously detected young gamma-ray pulsars. Follow-up radio observations have revealed faint radio pulsations from two of the newly discovered pulsars and enabled us to derive upper limits on the radio emission from the others, demonstrating that they are likely radio-quiet gamma-ray pulsars. We also present results from modeling the gamma-ray pulse profiles and radio profiles, if available, using different geometric emission models of pulsars. The high discovery rate of this survey, despite the increasing difficulty of blind pulsar searches in gamma rays, suggests that new systematic surveys such as presented in this article should be continued when new LAT source catalogs become available.

Fermi Gamma-Ray Pulsars: Understanding the High-energy Emission from Dissipative Magnetospheres

Abstract Based on the Fermi observational data, we reveal meaningful constraints for the dependence of the macroscopic conductivity (σ) of dissipative pulsar magnetosphere models on the corresponding spin-down rate, . Our models are refinements of the FIDO (Force-free Inside, Dissipative Outside) models, which have dissipative regions that are restricted on the equatorial current sheet outside of the the light-cylinder. Taking into account the observed cutoff energies of all of the Fermi pulsars and assuming that (a) the corresponding γ-ray pulsed emission is due to curvature radiation at the radiation-reaction-limit regime, and (b) this emission is produced at the equatorial current sheet near the light cylinder, we show that the Fermi data provide clear indications about the corresponding accelerating electric-field components. A direct comparison between the Fermi cutoff energies and the model ones reveals that σ increases with for high -values, while it saturates for low ones. This comparison indicates also that the corresponding gap width increases toward low -values. Assuming the Goldreich–Julian flux for the emitting particles, we calculate the total γ-ray luminosity (L γ ). A comparison between the dependence of the Fermi L γ -values and the model ones on indicates an increase of the emitting particle multiplicity with . Our modeling, guided by the Fermi data alone, enhances our understanding of the physical mechanisms behind the high-energy emission in pulsar magnetospheres.

Hubble Space Telescope Detection of the Millisecond Pulsar J2124−3358 and its Far-ultraviolet Bow Shock Nebula

Abstract We observed the nearby millisecond pulsar J2124–3358 with the Hubble Space Telescope in broad far-UV (FUV) and optical filters. The pulsar is detected in both bands with fluxes F(1250–2000 Å) = (2.5 ± 0.3) × 10−16 erg s−1 cm−2 and F(3800–6000 Å) = (6.4 ± 0.4) × 10−17 erg s−1 cm−2, which corresponds to luminosities of ≈5.8 × 1027 and 1.4 × 1027 erg s−1, for d = 410 pc and E(B − V) = 0.03. The optical-FUV spectrum can be described by a power-law model, , with slope α = 0.18–0.48 for a conservative range of color excess, E(B − V) = 0.01–0.08. Since a spectral flux rising with frequency is unusual for pulsar magnetospheric emission in this frequency range, it is possible that the spectrum is predominantly magnetospheric (power law with α < 0) in the optical, while it is dominated by thermal emission from the neutron star surface in the FUV. For a neutron star radius of 12 km, the surface temperature would be between 0.5 × 105 and 2.1 × 105 K for α ranging from −1 to 0, E(B − V) = 0.01–0.08, and d = 340–500 pc. In addition to the pulsar, the FUV images reveal extended emission that is spatially coincident with the known Hα bow shock, making PSR J2124–3358 the second pulsar (after PSR J0437−4715) with a bow shock detected in the FUV.

Polarized curvature radiation in pulsar magnetosphere

The propagation of polarized emission in pulsar magnetosphere is investigated in this paper. The polarized waves are generated through curvature radiation from the relativistic particles streaming along curved magnetic field lines and co-rotating with the pulsar magnetosphere. Within the 1/{\deg} emission cone, the waves can be divided into two natural wave mode components, the ordinary (O) mode and the extraord nary (X) mode, with comparable intensities. Both components propagate separately in magnetosphere, and are aligned within the cone by adiabatic walking. The refraction of O-mode makes the two components separated and incoherent. The detectable emission at a given height and a given rotation phase consists of incoherent X-mode and O-mode components coming from discrete emission regions. For four particle-density models in the form of uniformity, cone, core and patches, we calculate the intensities for each mode numerically within the entire pulsar beam. If the co-rotation of relativistic particles with magnetosphere is not considered, the intensity distributions for the X-mode and O-mode components are quite similar within the pulsar beam, which causes serious depolarization. However, if the co-rotation of relativistic particles is considered, the intensity distributions of the two modes are very different, and the net polarization of out-coming emission should be significant. Our numerical results are compared with observations, and can naturally explain the orthogonal polarization modes of some pulsars. Strong linear polarizations of some parts of pulsar profile can be reproduced by curvature radiation and subsequent propagation effect.

FERMILARGE AREA TELESCOPE OBSERVATIONS OF THE CRAB PULSAR AND NEBULA

We report on gamma-ray observations of the Crab Pulsar and Nebula using 8 months of survey data with the Fermi Large Area Telescope (LAT). The high quality light curve obtained using the ephemeris provided by the Nancay and Jodrell Bank radio telescopes shows two main peaks stable in phase with energy. The first gamma-ray peak leads the radio main pulse by (281 \pm 12 \pm 21) mus, giving new constraints on the production site of non-thermal emission in pulsar magnetospheres. The improved sensitivity and the unprecedented statistics afforded by the LAT enable precise measurement of the Crab Pulsar spectral parameters: cut-off energy at E_c = (5.8 \pm 0.5 \pm 1.2) GeV, spectral index of Gamma = (1.97 \pm 0.02 \pm 0.06) and integral photon flux above 100 MeV of (2.09 \pm 0.03 \pm 0.18) x 10^{-6} cm^{-2} s^{-1}. The first errors represent the statistical error on the fit parameters, while the second ones are the systematic uncertainties. Pulsed gamma-ray photons are observed up to ~ 20 GeV which precludes emission near the stellar surface, below altitudes of around 4 to 5 stellar radii in phase intervals encompassing the two main peaks. The spectrum of the nebula in the energy range 100 MeV - 300 GeV is well described by the sum of two power-laws of indices Gamma_{sync} = (3.99 \pm 0.12 \pm 0.08) and Gamma_{IC} = (1.64 \pm 0.05 \pm 0.07), corresponding to the falling edge of the synchrotron and the rising edge of the inverse Compton components, respectively. This latter, which links up naturally with the spectral data points of Cherenkov experiments, is well reproduced via inverse Compton scattering from standard Magnetohydrodynamics (MHD) nebula models, and does not require any additional radiation mechanism.

An empirical model for the beams of radio pulsars

Motivated by recent results on the location of the radio emission in pulsar magnetospheres, we have developed a model which can account for the large diversity found in the average profile shapes of pulsars. At the centre of our model lies the idea that radio emission at a particular frequency arises from a wide range of altitudes above the surface of the star, and that it is confined to a region close to the last open field lines. We assert that the radial height range over which emission occurs is responsible for the complex average pulse shapes rather than the transverse (longitudinal) range proposed in most current models. By implementing an abrupt change in the height range to discriminate between young, short-period, highly energetic pulsars and their older counterparts, we obtain the observed transition between the simple and complex average pulse profiles observed in each group respectively. Monte Carlo simulations are used to demonstrate the match of our model to real observations.

The Shadow of a Pulsar and the Inward Radio Emission in Pulsar Magnetosphere

We discuss observational facts that can be interpreted in terms of inward emission in pulsar magnetosphere. These include the main-pulse/interpulse anticorrelation in B1822-09 and the locations of emission and absorption features in the pulse profile of B0950+08. Weaknesses of geometrical models employing the inward emission are carefully discussed.

The induced turbulence effect on propagation of radio emission in pulsar magnetospheres

The effect of photon-beam-induced turbulence on propagation of radio emission in a pulsar magnetosphere is discussed. Beamed radio emission with a high brightness temperature can generate low-frequency plasma waves in the pulsar magnetosphere and these waves scatter the radio beam. We consider this effect on propagation of radio emission both in the open field line region and in the closed field line region. The former is applicable to most cases of pulsar radio emission where the propagation is confined to the polar region; it is shown that the induced process is not effective for radio emission of moderately high brightness temperature but can have a severe effect on giant pulses. For giant pulses not to be affected by this process, they must be emitted very close to the light cylinder. We show that the induced process is efficient in the closed field line region, inhibiting propagation of the radio emission in this region.

Modulated electrostatic modes in pair plasmas: Modulational stability profile and envelope excitations

A pair plasma consisting of two types of ions, possessing equal masses and opposite charges, is considered. The nonlinear propagation of modulated electrostatic wave packets is studied by employing a two-fluid plasma model. Considering propagation parallel to the external magnetic field, two distinct electrostatic modes are obtained, namely a quasiacoustic lower moddfe and a Langmuir-like, as optic-type upper one, in agreement with experimental observations and theoretical predictions. Considering small yet weakly nonlinear deviations from equilibrium, and adopting a multiple-scale technique, the basic set of model equations is reduced to a nonlinear Schrödinger equation for the slowly varying electric field perturbation amplitude. The analysis reveals that the lower (acoustic) mode is stable and may propagate in the form of a dark-type envelope soliton (a void) modulating a carrier wave packet, while the upper linear mode is intrinsically unstable, and may favor the formation of bright-type envelope soliton (pulse) modulated wave packets. These results are relevant to recent observations of electrostatic waves in pair-ion (fullerene) plasmas, and also with respect to electron-positron plasma emission in pulsar magnetospheres.

Centrifugal acceleration of plasma in pulsar magnetosphere

We present a relativistic model for the centrifugal acceleration of plasma bunches and the coherent radio emission in pulsar magnetosphere. We find that rotation broadens the width of leading component compared to the width of trailing component. We explain this difference in the component widths using the nested cone emission geometry. We estimate the effect of pulsar spin on the Stokes parameters, and find that the inclination between the rotation and magnetic axes can introduce an asymmetry in the circular polarization of the conal components. We analyse the single pulse polarization data of PSR B032954 at 606 MHz, and find that in its conal components, one sense of circular polarization dominates in the leading component while the other sense dominates in the trailing component. Our simulation shows that changing the sign of the impact parameter changes the sense of circular polarization as well as the swing of polarization angle

A search for the optical counterpart to PSR B1821-24 in M 28

We have analysed archival HST/WFPC2 images in both the F555W & F814W bands of the core field of the globular cluster M 28 in an attempt to identify the optical counterpart of the magnetospherically active millisecond pulsar PSR B1821-24. Examination of the radio derived error circle yielded several potential candidates, down to a magnitude of V ~ 24.5 (V 0 ~ 23.0). Each were further investigated, both in the context of the CMD of M 28, and also with regard to phenomenological models of pulsar magnetospheric emission. The latter was based on both luminosity-spindown correlations and known spectral flux density behaviour in this regime from the small population of optical pulsars observed to date. None of the potential candidates exhibited emission expected from a magnetospherically active pulsar. The fact that the magnetic field & spin coupling for PSR B1821-24 is of a similar magnitude to that of the Crab pulsar in the vicinity of the light cylinder has suggested that the millisecond pulsar may well be an efficient nonthermal emitter. ASCA's detection of a strong synchrotron-dominated X-ray pulse fraction encourages such a viewpoint. We argue that only future dedicated 2-d high speed photometry observations of the radio error-circle can finally resolve this matter.

A new mechanism for pulsar gamma-ray emission

Gamma-ray emission in pulsar magnetospheres is attributed to synchrotron radiation, which tends to decrease the pitch angle of the particle, being balanced by plasma processes tending to increase the pitch angle. The plasma processes are non-resonant instabilities that drive non-resonant quasilinear diffusion (NQD), thereby pumping energy from waves and the parallel motion of the particle into the perpendicular motion of the particle. It is shown that NQD can maintain the pitch angles for particles near the light-cylinder such that they radiate synchrotron radiation at MeV energies. Compared to conventional emission mechanisms (such as polar cap or outer gap models), the resulting spectrum has a relatively low upper cut-off from about a few to 100 MeV. Possible observational consequences of this mechanism are discussed.

Pulsar Magnetospheric Emission Mapping: Images and Implications of Polar Cap Weather

The beautiful sequences of ``drifting'' subpulses observed in some radio pulsars have been regarded as among the most salient and potentially instructive characteristics of their emission, not least because they have appeared to represent a system of subbeams in motion within the emission zone of the star. Numerous studies of these ``drift'' sequences have been published, and a model of their generation and motion articulated long ago by Ruderman & Sutherland (1975); but efforts thus far have failed to establish an illuminating connection between the drift phemomenon and the actual sites of radio emission. Through a detailed analysis of a nearly coherent sequence of ``drifting'' pulses from pulsar B0943+10, we have in fact identified a system of subbeams circulating around the magnetic axis of the star. A mapping technique, involving a ``cartographic'' transform and its inverse, permits us to study the character of the polar-cap emission ``map'' and then to confirm that it, in turn, represents the observed pulse sequence. On this basis, we have been able to trace the physical origin of the ``drifting-subpulse'' emission to a stably rotating and remarkably organized configuration of emission columns, in turn traceable possibly to the magnetic polar-cap ``gap'' region envisioned by some theories.

Detection of Magnetospheric X-Ray Pulsation from the Millisecond Pulsar B1821−24

We report the first clear detection of X-ray pulsation of magnetospheric origin from a millisecond pulsar, PSR B1821-24, with the Advanced Satellite for Cosmology and Astrophysics. The photon arrival time shows a periodicity of 3.05 ms period, as observed in the radio band. The observed X-ray pulse is double peaked. The pulses are characterized by a sharp temporal profile and hard power-law spectrum with a phase-averaged photon index of ~1.9. These two features are quite similar to the X-ray/γ-ray pulses from the Crab pulsar and characteristic of the nonthermal emission from the magnetosphere of the neutron star, but contradictory to the thermal emission from its surface. Since this pulsar has physical parameters completely different from those of young pulsars, we expect these observations to provide significant constraints on models of pulsar magnetospheric emission, which have been developed mostly based on the observations of young pulsars.