Millisecond Pulsar PSR J1023 Research Articles

Constraining the magnetic field geometry of the millisecond pulsar PSR J0030+0451 from joint radio, thermal X-ray, andγ-ray emission

Context.With the advent of multi-wavelength electromagnetic observations of neutron stars – spanning many decades in photon energies – from radio wavelengths up to X-rays andγ-rays, it has become possible to significantly constrain the geometry and the location of the associated emission regions.Aims.In this work, we use results from the modelling of thermal X-ray observations of PSR J0030+0451 from the Neutron Star Interior Composition Explorer (NICER) mission and phase-aligned radio andγ-ray pulse profiles to constrain the geometry of an off-centred dipole that is able to reproduce the light curves in these respective bands simultaneously.Methods.To this aim, we deduced a configuration with a simple dipole off-centred from the location of the centre of the thermal X-ray hot spots. We show that the geometry is compatible with independent constraints from radio andγ-ray pulsations only, leading to a fixed magnetic obliquity ofα ≈ 75° and a line-of-sight inclination angle ofζ ≈ 54°.Results.We demonstrate that an off-centred dipole cannot be rejected by accounting for the thermal X-ray pulse profiles. Moreover, the crescent shape of one spot is interpreted as the consequence of a small-scale surface dipole on top of the large-scale off-centred dipole.

Search for continuous gravitational waves from PSR J0437−4715 with a hidden Markov model in O3 LIGO data

Results are presented for a semi-coherent search for continuous gravitational waves from the millisecond pulsar PSR J0437$-$4715, using a hidden Markov model to track spin wandering, in LIGO data from the third LIGO-Virgo observing run. This is the first search for PSR J0437$-$4715 to cover a wide frequency range from $60$ Hz to $500$ Hz and simultanously accommodate random spin deviations from the secular radio ephemeris. Two searches are performed with plausible coherence times of $10$ days and $30$ days, as the frequency wandering time-scale of the gravitational-wave-emitting quadrupole is unknown. The former analysis yields no surviving candidates, while the latter yields five candidates after the veto procedure. The detection statistic of each of the five survivors is mapped as a function of sky position, in preparation for follow-up analyses in the future, e.g. during LIGO-Virgo-KAGRA fourth observing run.

Measurements of pulse jitter and single-pulse variability in millisecond pulsars using MeerKAT

ABSTRACT Using the state-of-the-art SKA precursor, the MeerKAT radio telescope, we explore the limits to precision pulsar timing of millisecond pulsars achievable due to pulse stochasticity (jitter). We report new jitter measurements in 15 of the 29 pulsars in our sample and find that the levels of jitter can vary dramatically between them. For some, like the 2.2 ms pulsar PSR J2241−5236, we measure an implied jitter of just ∼4 ns h−1, while others, like the 3.9 ms PSR J0636−3044, are limited to ∼100 ns h−1. While it is well known that jitter plays a central role to limiting the precision measurements of arrival times for high signal-to-noise ratio observations, its role in the measurement of dispersion measure (DM) has not been reported, particularly in broad-band observations. Using the exceptional sensitivity of MeerKAT, we explored this on the bright millisecond pulsar PSR J0437−4715 by exploring the DM of literally every pulse. We found that the derived single-pulse DMs vary by typically 0.0085 cm−3 pc from the mean, and that the best DM estimate is limited by the differential pulse jitter across the band. We postulate that all millisecond pulsars will have their own limit on DM precision which can only be overcome with longer integrations. Using high-time resolution filterbank data of 9 μs, we also present a statistical analysis of single-pulse phenomenology. Finally, we discuss optimization strategies for the MeerKAT pulsar timing program and its role in the context of the International Pulsar Timing Array.

Three-dimensional Modeling of the Magnetothermal Evolution of Neutron Stars: Method and Test Cases

Abstract Neutron stars harbor extremely strong magnetic fields within their solid outer crust. The topology of this field strongly influences the surface temperature distribution and, hence, the star’s observational properties. In this work, we present the first realistic simulations of the coupled crustal magnetothermal evolution of isolated neutron stars in three dimensions accounting for neutrino emission, obtained with the pseudo-spectral code parody. We investigate both the secular evolution, especially in connection with the onset of instabilities during the Hall phase, and the short-term evolution following episodes of localized energy injection. Simulations show that a resistive tearing instability develops in about a Hall time if the initial toroidal field exceeds G. This leads to crustal failures because of the huge magnetic stresses coupled with the local temperature enhancement produced by dissipation. Localized heat deposition in the crust results in the appearance of hot spots on the star surface, which can exhibit a variety of patterns. Because the transport properties are strongly influenced by the magnetic field, the hot regions tend to drift away and get deformed following the magnetic field lines while cooling. The shapes obtained with our simulations are reminiscent of those recently derived from NICER X-ray observations of the millisecond pulsar PSR J0030+0451.

A Numerical Model for the Multiwavelength Lightcurves of PSR J0030+0451

Abstract Recent modeling of Neutron Star Interior Composition Explorer (NICER) observations of the millisecond pulsar PSR J0030+0451 suggests that the magnetic field of the pulsar is non-dipolar. We construct a magnetic field configuration where foot points of the open field lines closely resemble the hotspot configuration from NICER observations. Using this magnetic field as input, we perform force-free simulations of the magnetosphere of PSR J0030+0451, showing the three-dimensional structure of its plasma-filled magnetosphere. Making simple and physically motivated assumptions about the emitting regions, we are able to construct the multiwavelength lightcurves that qualitatively agree with the corresponding observations. The agreement suggests that multipole magnetic structures are the key to modeling this type of pulsar, and can be used to constrain the magnetic inclination angle and the location of radio emission.

A NICER View of PSR J0030+0451: Evidence for a Global-scale Multipolar Magnetic Field

Abstract Recent modeling of Neutron Star Interior Composition Explorer observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar’s surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451’s magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission.

A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation

Abstract We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030+0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar’s surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius R eq are, respectively, 1.34 − 0.16 + 0.15 M ⊙ and 12.71 − 1.19 + 1.14 km , while the compactness GM / R eq c 2 = 0.156 − 0.010 + 0.008 is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

A NICER View of PSR J0030+0451: Implications for the Dense Matter Equation of State

Abstract Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA’s Neutron Star Interior Composition Explorer (NICER) mission. In this Letter we study the implications of the mass–radius inference reported for this source by Riley et al. for the dense matter equation of state (EoS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EoS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EoS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EoS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass–radius likelihood function derived by Riley et al. using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.

PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter

Abstract Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state (EoS) of this dense matter is to measure both a star’s equatorial circumferential radius R e and its gravitational mass M. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are km and (68%). The independent analysis reported in the companion paper by Riley et al. explores different emitting spot models, but finds spot shapes and locations and estimates of R e and M that are consistent with those found in this work. We show that our measurements of R e and M for PSR J0030+0451 improve the astrophysical constraints on the EoS of cold, catalyzed matter above nuclear saturation density.

Neutron star radius measurement from the ultraviolet and soft X-ray thermal emission of PSR J0437−4715

ABSTRACT We analysed the thermal emission from the entire surface of the millisecond pulsar PSR J0437−4715 observed in the ultraviolet and soft X-ray bands. For this, we calculated non-magnetized, partially ionized atmosphere models of hydrogen, helium, and iron compositions and included plasma frequency effects that may affect the emergent spectrum. This is particularly true for the coldest atmospheres composed of iron (up to a few per cent changes in the soft X-ray flux). Employing a Markov chain Monte Carlo method, we found that the spectral fits favour a hydrogen atmosphere, disfavour a helium composition, and rule out iron atmosphere and blackbody models. By using a Gaussian prior on the dust extinction, based on the latest 3D map of Galactic dust, and accounting for the presence of hot polar caps found in the previous work, we found that the hydrogen atmosphere model results in a well-constrained neutron star radius ${R_{\rm NS}}= 13.6^{+0.9}_{-0.8}{\, {\rm km}}$ and bulk surface temperature ${T_{\rm eff}^{\infty }}=\left(2.3\pm 0.1\right){\times 10^{5}}{\, {\rm K}}$. This relatively large radius favours a stiff equation of state and disfavours a strange quark composition inside neutron stars.

The broad X-ray emission of the millisecond pulsar PSR J0437–4715

AbstractThe nearest millisecond pulsar PSR J0437–4715 is the ideal target to constrain the dense matter equation of state using the lightcurve modelling method. Our analysis combining XMM-Newton, NuSTAR, and ROSAT observations removed ambiguities in the spectral modelling of the surface emission from PSR J0437–4715. Furthermore, the NuSTAR observation demonstrated that the non-thermal hard tail emission was pulsed at the pulsar spin period. These features are crucial to model the lightcurve and to measure the radius of the neutron star. This conference proceeding is based on the publication Guillot et al. (2016).

SCINTILLATION ARCS IN LOW-FREQUENCY OBSERVATIONS OF THE TIMING-ARRAY MILLISECOND PULSAR PSR J0437–4715

ABSTRACT Low-frequency observations of pulsars provide a powerful means for probing the microstructure in the turbulent interstellar medium (ISM). Here we report on high-resolution dynamic spectral analysis of our observations of the timing-array millisecond pulsar PSR J0437–4715 with the Murchison Widefield Array (MWA), enabled by our recently commissioned tied-array beam processing pipeline for voltage data recorded from the high time resolution mode of the MWA. A secondary spectral analysis reveals faint parabolic arcs akin to those seen in high-frequency observations of pulsars with the Green Bank and Arecibo telescopes. Data from Parkes observations at a higher frequency of 732 MHz reveal a similar parabolic feature with a curvature that scales approximately as the square of the observing wavelength (λ 2) to the MWA's frequency of 192 MHz. Our analysis suggests that scattering toward PSR J0437–4715 predominantly arises from a compact region about 115 pc from the Earth, which matches well with the expected location of the edge of the Local Bubble that envelopes the local Solar neighborhood. As well as demonstrating new and improved pulsar science capabilities of the MWA, our analysis underscores the potential of low-frequency pulsar observations for gaining valuable insights into the local ISM and for characterizing the ISM toward timing-array pulsars.

Rotochemical heating of millisecond and classical pulsars with anisotropic and density-dependent superfluid gap models

When a rotating neutron star loses angular momentum, the progressive reduction of the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium, inducing reactions that release heat (rotochemical heating). This effect has previously been studied by Fern\'andez & Reisenegger (2005) for non-superfluid neutron stars and by Petrovich & Reisenegger (2010) for superfluid millisecond pulsars. Both studies found that pulsars reach a quasi-steady state in which the compression driving the matter out of beta equilibrium is balanced by the reactions trying to restore the equilibrium. We extend previous studies by considering the effect of density-dependence and anisotropy of the superfluid energy gaps, for the case in which the dominant reactions are the modified Urca processes, the protons are non-superconducting, and the neutron superfluidity is parametrized by models proposed in the literature. By comparing our predictions with the surface temperature of the millisecond pulsar PSR J0437-4715 and upper limits for twenty-one classical pulsars, we find the millisecond pulsar can be only explained by the models with the effectively largest energy gaps (type B models), the classical pulsars require with the gap models that vanish for some angle (type C) and two different envelope compositions. Thus, no single model for neutron superfluidity can simultaneously account for the thermal emission of all available observations of non-accreting neutron stars, possibly due to our neglect of proton superconductivity.

THE LOW-FREQUENCY CHARACTERISTICS OF PSR J0437–4715 OBSERVED WITH THE MURCHISON WIDE-FIELD ARRAY

We report on the detection of the millisecond pulsar PSR J0437-4715 with the Murchison Widefield Array (MWA) at a frequency of 192 MHz. Our observations show rapid modulations of pulse intensity in time and frequency that arise from diffractive scintillation effects in the interstellar medium (ISM), as well as prominent drifts of intensity maxima in the time-frequency plane that arise from refractive effects. Our analysis suggests that the scattering screen is located at a distance of $\sim$80-120 pc from the Sun, in disagreement with a recent claim that the screen is closer ($\sim$10 pc). Comparisons with higher frequency data from Parkes reveals a dramatic evolution of the pulse profile with frequency, with the outer conal emission becoming comparable in strength to that from the core and inner conal regions. As well as demonstrating high time resolution science capabilities currently possible with the MWA, our observations underscore the potential to conduct low-frequency investigations of timing-array millisecond pulsars, which may lead to increased sensitivity for the detection of nanoHertz gravitational waves via the accurate characterisation of ISM effects.

Timing, polarimetry and physics of the bright, nearby millisecond pulsar PSR J0437−4715 – a single-pulse perspective

Single pulses from radio pulsars contain a wealth of information about emission and propagation in the magnetosphere and insight into their timing properties. It was recently demonstrated that single-pulse emission is responsible for limiting the timing stability of the brightest of millisecond pulsars. We report on an analysis of more than a million single-pulses from PSR J0437-4715 and present various statistical properties such as the signal-to-noise ratio (S/N) distribution, timing and polarimetry of average profiles integrated from subpulses with chosen S/N cut-offs, modulation properties of the emission, phase-resolved statistics of the S/N, and two dimensional spherical histograms of the polarization vector orientation. The last of these indicates the presence of orthogonally polarised modes (OPMs). Combined with the dependence of the polarisation fraction on the S/N and polarimetry of the brightest pulses, the existence of OPMs constrains pulsar emission mechanisms and models for the plasma physics in the magnetosphere.

EVIDENCE FOR GAMMA-RAY EMISSION FROM THE LOW-MASS X-RAY BINARY SYSTEM FIRST J102347.6+003841

The low-mass X-ray binary (LMXB) system FIRST J102347.6+003841 hosts a newly born millisecond pulsar (MSP) PSR J1023+0038 that was revealed as the first and only known rotation-powered MSP in a quiescent LMXB. While the system is shown to have an accretion disk before 2002, it remains unclear how the accretion disk has been removed in order to reveal the radio pulsation in 2007. In this Letter, we report the discovery of gamma-rays spatially consistent with FIRST J102347.6+003841, at a significance of 7 standard deviations, using data obtained by the Fermi Gamma-ray Space Telescope. The gamma-ray spectrum can be described by a power law (PL) with a photon index of 2.9+-0.2, resulting in an energy flux above 200 MeV of (5.5+-0.9)x10^{-12} erg cm^{-2}s^{-1}. The gamma-rays likely originate from the MSP PSR J1023+0038, but also possibly from an intrabinary shock between the pulsar and its companion star. To complement the gamma-ray study, we also re-investigate the XMM-Newton data taken in 2004 and 2008. Our X-ray spectral analysis suggests that a broken PL with two distinct photon indices describes the X-ray data significantly better than a single PL. This indicates that there exists two components and that both components appear to vary with the orbital phase. The evidence for gamma-ray emission conforms with a recent suggestion that gamma-rays from PSR J1023+0038 may be responsible for ejecting the disk material out of the system.

SEARCH FOR VERY HIGH ENERGY GAMMA-RAY EMISSION FROM PULSAR-PULSAR WIND NEBULA SYSTEMS WITH THE MAGIC TELESCOPE

The MAGIC collaboration has searched for high-energy gamma-ray emission of some of the most promising pulsar candidates above an energy threshold of 50 GeV, an energy not reachable up to now by other ground-based instruments. Neither pulsed nor steady gamma-ray emission has been observed at energies of 100 GeV from the classical radio pulsars PSR J0205+6449 and PSR J2229+6114 (and their nebulae 3C58 and Boomerang, respectively) and the millisecond pulsar PSR J0218+4232. Here, we present the flux upper limits for these sources and discuss their implications in the context of current model predictions.

PULSED GAMMA RAYS FROM THE MILLISECOND PULSAR J0030+0451 WITH THEFERMILARGE AREA TELESCOPE

We report the discovery of gamma-ray pulsations from the nearby isolated millisecond pulsar PSR J0030+0451 with the Large Area Telescope (LAT) on the \emph{Fermi} Gamma-ray Space Telescope (formerly GLAST). This discovery makes PSR J0030+0451 the second millisecond pulsar to be detected in gamma-rays after PSR J0218+4232, observed by the EGRET instrument on the Compton Gamma Ray Observatory. The spin-down power $\dot E = $ 3.5 $\times$ 10$^{33}$ ergs s$^{-1}$ is an order of magnitude lower than the empirical lower bound of previously known gamma-ray pulsars. The emission profile is characterized by two narrow peaks, respectively 0.07 $\pm$ 0.01 and 0.08 $\pm$ 0.02 wide, separated by 0.44 $\pm$ 0.02 in phase. The first gamma-ray peak falls 0.15 $\pm$ 0.01 after the main radio peak. The pulse shape is similar to that of the normal gamma-ray pulsars. An exponentially cut-off power-law fit of the emission spectrum leads to an integral photon flux above 100 MeV of (6.76 $\pm$ 1.05 $\pm$ 1.35) $\times 10^{-8}$ cm$^{-2}$ s$^{-1}$ with cut-off energy (1.7 $\pm$ 0.4 $\pm$ 0.5) GeV. Based on its parallax distance of $(300 \pm 90)$ pc, we obtain a gamma-ray efficiency $L_\gamma / \dot{E} \simeq 15%$ for the conversion of spin-down energy rate into gamma-ray radiation, assuming isotropic emission.

Precision Timing of PSR J0437−4715: An Accurate Pulsar Distance, a High Pulsar Mass, and a Limit on the Variation of Newton’s Gravitational Constant

Analysis of ten years of high-precision timing data on the millisecond pulsar PSR J0437-4715 has resulted in a model-independent kinematic distance based on an apparent orbital period derivative, Pbdot, determined at the 1.5% level of precision (Dk = 157.0 +/- 2.4 pc), making it one of the most accurate stellar distance estimates published to date. The discrepancy between this measurement and a previously published parallax distance estimate is attributed to errors in the DE200 Solar System ephemerides. The precise measurement of Pbdot allows a limit on the variation of Newton's gravitational constant, |Gdot/G| < 23 x 10^{-12} 1/yr. We also constrain any anomalous acceleration along the line of sight to the pulsar to |a(Sun)/c| < 1.5 x 10^{-18} 1/s at 95% confidence, and derive a pulsar mass, m(psr) = 1.76 +/- 0.20 M, one of the highest estimates so far obtained.

Empirical Constraints on the General Relativistic Electric Field Associated with PSR J0437-4715

We simulate the magnetosphere of the nearby millisecond pulsar PSR J0437-4715, which is expected to have an unscreened electric potential due to the lack of magnetic pair production. We incorporate general relativistic (GR) effects and study curvature radiation (CR) by primary electrons but neglect inverse Compton scattering of thermal X-ray photons by these electrons. We find that the CR spectrum cuts off at energies below ~17 GeV, well below the threshold of the High Energy Stereoscopic System (H.E.S.S.) telescope (100 GeV), while other models predict a much higher cutoff of 100 GeV. GR theory also predicts a relatively narrow pulse (βo ~ 0.2 phase width) centered on the magnetic axis. EGRET observations above 100 MeV significantly constrain the application of the Muslimov & Harding model for γ-ray production as a result of GR frame dragging and ultimately its polar cap (PC) current and accelerating potential. Whereas the standard prediction of this pulsar's γ-ray luminosity due to GR frame dragging is ~10% of the spin-down power, a nondetection by forthcoming H.E.S.S. observations will constrain it to 0.3%, enforcing an even more severe revision of the accelerating electric field and PC current.