Millisecond Pulsars Research Articles

The Relationships between Emission Geometry and Jitter Noise in Millisecond Pulsars

The relationships between several parameters of an emission geometry and jitter noise in 17 millisecond pulsars (MSPs) are investigated. By assuming the jitter noise is due only to a pulse variation in phase, the former can be modeled as changes in the plasma flow rate leading to variation in the measured pulse arrival time relative to the predicted time. In the model for pulsar magnetospheres with multiple emission states, the plasma flow is associated with the emission states, and a change in the emission state corresponds to a change in the plasma flow causing variation in the pulse arrival time. These can be specified in an emission geometry defined by the obliquity and viewing angles, measured from the rotation axis to the magnetic axis and to the line of sight, respectively. We calculate the maximum change in the emission state based on the reported jitter noise for each of the MSPs. Using the results, we show that the MSPs possess relatively large obliquity angles, which is consistent with observation, and the jitter noise exhibits dependency on frequency. We find that the jitter noise in our sample displays an exponential decay as a function that combines the obliquity angle and the rotation period, revealing the correlation among the three parameters. This suggests that the magnitude of the jitter noise is likely specific to an MSP. We discuss how jitter noise may be related to the evolution of an MSP.

Non-radial oscillations in anisotropic dark energy stars

We study the non-radial f-mode oscillations of both isotropic and anisotropic dark energy stars by using the modified Chaplygin prescription of dark energy to model the stellar matter. The anisotropic pressure in the system is modeled with Bowers–Liang prescription. By solving the stellar structure equations in presence of anisotropy, we study the global properties of the dark energy star and compare the mass-radius profiles with data from GW events and milli-second pulsars. The stellar radial and anisotropic pressure profiles have also been investigated. We proceed to determine the prominent non-radial l=2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$l=2$$\\end{document}f-mode frequencies of the anisotropic dark energy star by employing the Cowling approximation and analyse and quantify the spectra by varying the anisotropic parameter. We report that f-mode spectra of dark energy star have distinctly different behaviour compared to neutron star and quark star, and this may possibly help in its future identification. Further, the tidal deformability factors of the anisotropic dark energy stars have also been analyzed.

A large population of neutron star low-mass X-ray binaries with long outburst recurrence time?

ABSTRACT Low-mass X-ray binaries (LMXBs) with neutron stars show quite different features that depend on the rate of mass transfer from the donor star. With a high transfer rate, the Z sources are in a persistent soft spectral state, and with a moderate transfer rate the transient Atoll sources have outburst cycles like the black hole X-ray binaries. The observations document very long outburst recurrence times for quite a number of sources. We follow with our computations the evolution of the accretion disc until the onset of the ionization instability. For sources with a low mass transfer rate, the accumulation of matter in the disc is essentially reduced due to the continuous evaporation of matter from the disc to the coronal flow. Different mass transfer rates result in nearly the same amount of matter accumulated for the outburst, which means that the outburst properties are similar for sources with short outburst cycles and sources with long outburst cycles, contrary to some expectations. Then, for systems with long recurrence time, less sources will be detected and the total population of LMXBs could be larger than it appears. This would relieve the apparent problem that the observed number of LMXBs as progenitors of millisecond pulsars (MSPs) is too small compared to the number of MSPs. Concerning the few quasi-persistent sources with year-long soft states, we argue that these states are not outbursts, but quasi-stationary hot states as in Z sources.

Heating of millisecond pulsars by magnetic field decay

AbstractMillisecond pulsars (MSPs) are believed to be very old neutron stars (NSs) whose age may exceed significantly years. Although cooling scenarios of isolated NSs predict for that age a surface temperature K, observations of the nearest MSP J0437‐4715 indicate well above that value. Besides the heating of the polar cap surface by backflowing charged particles, Joule heating in the crust can contribute to the overall heat budget of MSPs. Because the dipolar field component, derived from and measurements, is much too weak for remarkable heating, smaller‐scale structures should be analyzed to understand whether they can supply the demanded heat. For this purpose, we study the small‐scale field structure of radio pulsars. Magnetic field components, significantly stronger than the dipolar one, may exist especially at the surface of MSPs. We assign upper limits to the strength of single‐field components up to a multipolarity of and the corresponding deviations from axial symmetry . Arguments are provided that the decay of the small‐scale components with or of the crustal magnetic field may cause the relatively high surface temperature of isolated MSPs.

A 1.9 solar-mass neutron star candidate in a 2-year orbit

We report discovery and characterization of a main-sequence G star orbiting a dark object with mass 1.90±0.04M⊙. The system was discovered via Gaia astrometry and has an orbital period of 731 days. We obtained multi-epoch RV follow-up over a period of 639 days, allowing us to refine the Gaia orbital solution and precisely constrain the masses of both components. The luminous star is a ≳12,Gyr-old, low-metallicity halo star near the main-sequence turnoff (,K; ; ; M≈0.79M⊙) with a highly enhanced lithium abundance. The RV mass function sets a minimum companion mass for an edge-on orbit of M2>1.67M⊙, well above the Chandrasekhar limit. The Gaia inclination constraint, i=68.7±1.4,deg, then implies a companion mass of M2=1.90±0.04M⊙. The companion is most likely a massive neutron star: the only viable alternative is two massive white dwarfs in a close binary, but this scenario is disfavored on evolutionary grounds. The system’s low eccentricity ( e=0.122±0.002) disfavors dynamical formation channels and implies that the neutron star likely formed with little mass loss ( ≲1M⊙) and with a weak natal kick (). Stronger kicks with more mass loss are not fully ruled out but would imply that a larger population of similar systems with higher eccentricities should exist. The current orbit is too small to have accommodated the neutron star progenitor as a red supergiant or super-AGB star. The simplest formation scenario – isolated binary evolution – requires the system to have survived unstable mass transfer and common envelope evolution with a donor-to-accretor mass ratio >10. The system, which we call Gaia NS1, is likely a progenitor of symbiotic X-ray binaries and long-period millisecond pulsars. Its discovery challenges binary evolution models and bodes well for Gaia’s census of compact objects in wide binaries.

New Timing Results of MSPs from NICER Observations

Millisecond pulsars (MSPs) are known for their long-term stability. Using six years of observations from the Neutron Star Interior Composition Explorer (NICER), we have conducted an in-depth analysis of the X-ray timing results for six MSPs: PSRs B1937+21, B1821-24, J0437-4715, J0030+0451, J0218+4232, and J2124-3358. The timing stability parameter σz has been calculated, revealing remarkable timing precision on the order of 10−14 for PSRs B1937+21 and J0437-4715, and 10−13 for PSRs B1821-24, J0218+4232, and J0030+0451 over a timescale of 1000 days. These findings underscore the feasibility of autonomous in-orbit timekeeping using X-ray observations of MSPs. In addition, the consistency of long-term spin-down noise in the X-ray and radio bands has been investigated by comparison with IPTA radio data.

The thousand-pulsar-array programme on MeerKAT XIII: Timing, flux density, rotation measure and dispersion measure timeseries of 597 pulsars

Abstract We report here on the timing of 597 pulsars over the last four years with the MeerKAT telescope. We provide Times-of-Arrival, pulsar ephemeris files and per-epoch measurements of the flux density, dispersion measure (DM) and rotation measure (RM) for each pulsar. In addition we use a Gaussian process to model the timing residuals to measure the spin frequency derivative at each epoch. We also report the detection of 11 glitches in 9 individual pulsars. We find significant DM and RM variations in 87 and 76 pulsars respectively. We find that the DM variations scale approximately linearly with DM, which is broadly in agreement with models of the ionised interstellar medium. The observed RM variations seem largely independent of DM, which may suggest that the RM variations are dominated by variations in the interstellar magnetic field on the line of sight, rather than varying electron density. We also find that normal pulsars have around 5 times greater amplitude of DM variability compared to millisecond pulsars, and surmise that this is due to the known difference in their velocity distributions.

Generic gravito-magnetic clock effects

ABSTRACT General relativity predicts that two counter-orbiting clocks around a spinning mass differ in the time required to complete the same orbit. The difference in these two values for the orbital period is generally referred to as the gravito-magnetic (GM) clock effect. It has been proposed to measure the GM clock effect using atomic clocks carried by satellites in prograde and retrograde orbits around the Earth. The precision and stability required for satellites to accurately perform this measurement remains a challenge for current instrumentation. One of the most accurate clocks in the Universe is a millisecond pulsar, which emits periodic radio pulses with high stability. Timing of the pulsed signals from millisecond pulsars has proven to be very successful in testing predictions of general relativity and the GM clock effect is potentially measurable in binary systems. In this work, we derive the generic GM clock effect by considering a slowly spinning binary system on an elliptical orbit, with both arbitrary mass ratio and arbitrary spin orientations. The spin–orbit interaction introduces a perturbation to the orbit, causing the orbital plane to precess and nutate. We identify several different contributions to the clock effects: the choice of spin supplementary condition and the observer-dependent definition of a full revolution and ‘nearly identical’ orbits. We discuss the impact of these subtle definitions on the formula for GM clock effects and show that most of the existing formulae in the literature can be recovered under appropriate assumptions.

Radio Polarization of Millisecond Pulsars with Multipolar Magnetic Fields

NICER has observed a few millisecond pulsars where the geometry of the X-ray-emitting hotspots on the neutron star have been analyzed in order to constrain the mass and radius from X-ray light-curve modeling. One example, PSR J0030 + 0451, has been shown to possibly have significant multipolar magnetic fields at the stellar surface. Using force-free simulations of the magnetosphere structure, it has been shown that the radio, X-ray, and γ-ray light curves can be modeled simultaneously with an appropriate field configuration. An even more stringent test is to compare predictions of the force-free magnetosphere model with observations of radio polarization. This paper attempts to reproduce the radio polarization of PSR J0030 + 0451 using a force-free magnetospheric solution. As a result of our modeling, we can reproduce certain features of the polarization well.

Steering the Frequency of Hydrogen Master Based on Pulsar Observations

Hydrogen masters have high short-term stability, using them as master clocks can generate high-precision local standard time in a short period of time. However, the hydrogen maser has a frequency drift phenomenon, which leads to its poor long-term stability, thus affecting the accuracy of the local time. The rotation of the millisecond pulsar is highly stable. With the help of its high long-term stability, it can regularly control the frequency of the hydrogen atomic clock, thus control the real-time signal. The stability of four millisecond pulsars in the second data set released by the International Pulsar Timing Array (IPTA) is analyzed. At the same time, the frequency stability of a hydrogen maser of the National Time Service Center of the Chinese Academy of Sciences (NTSC) is analyzed by Hadamard variance. Finally, a strategy for steering the frequency of hydrogen master using pulsars is given.

The GMRT High-resolution Southern Sky Survey for Pulsars and Transients. VII. Timing of the Spider Millisecond Pulsar PSR J1242–4712, a Bridge between Redback and Black Widow Pulsars

We present a timing solution for the 5.31 ms spider millisecond pulsar (MSP) J1242−4712, discovered with the GMRT. PSR J1242−4712 orbits a companion of minimum mass 0.08 M ⊙ with an orbital period of 7.7 hr and occupies a relatively unexplored region in the orbital period versus companion mass space. We did not detect gamma-ray pulsations for this MSP, and also could not identify the optical counterpart for PSR J1242–4712 in available optical/near-infrared data. The profile of J1242−4712 evolves with frequency, showing a clear single component at lower frequencies and a three-component profile at 650 MHz. PSR J1242−4712 eclipses for a very short duration near superior conjunction (orbital phase ∼ 0.23−0.25) below 360 MHz. Moreover, significant dispersion measure delays and errors in the pulse times of arrivals are observed near inferior conjunction (orbital phase ∼ 0.7), along with an observed eclipse in one epoch at 650 MHz. Observed eclipses and significant orbital period variability suggest that PSR J1242−4712 is possibly not a helium star−white dwarf binary, but has a semi- or nondegenerate companion, indicating that this is a “spider” MSP lying in a region between typical black widows and redbacks. This system may represent a distinct category of spider MSPs, displaying characteristics that bridge the gap between known black widow and redback MSPs.

A targeted radio pulsar survey of redback candidates with MeerKAT

ABSTRACT Redbacks are millisecond pulsar binaries with low-mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localization regions of unidentified but pulsar-like Fermi-LAT gamma-ray sources. However, these candidates remain unconfirmed until pulsations are detected. As part of the TRAPUM project, we searched for radio pulsations from six of these redback candidates with MeerKAT. We discovered three new radio millisecond pulsars, PSRs J0838−2827, J0955−3947, and J2333−5526, confirming their redback nature. PSR J0838−2827 remained undetected for 2 yr after our discovery despite repeated observations, likely due to evaporated material absorbing the radio emission for long periods of time. While, to our knowledge, this system has not undergone a transition to an accreting state, the disappearance, likely caused by extreme eclipses, illustrates the transient nature of spider pulsars and the heavy selection bias in uncovering their radio population. Radio timing enabled the detection of gamma-ray pulsations from all three pulsars, from which we obtained 15-yr timing solutions. All of these sources exhibit complex orbital period variations consistent with gravitational quadrupole moment variations in the companion stars. These timing solutions also constrain the binary mass ratios, allowing us to narrow down the pulsar masses. We find that PSR J2333−5526 may have a neutron star mass in excess of 2 M⊙.

Characterizing the Gamma-Ray Emission Properties of the Globular Cluster M5 with the Fermi-LAT

We analyzed the globular cluster M5 (NGC 5904) using 15 yr of gamma-ray data from the Fermi Large Area Telescope (LAT). Using rotation ephemerides generated from Arecibo and FAST radio telescope observations, we searched for gamma-ray pulsations from the seven millisecond pulsars (MSPs) identified in M5. We detected no significant pulsations from any of the individual pulsars. In addition, we searched for possible variations of the gamma-ray emission as a function of orbital phase for all six MSPs in binary systems, but we did not detect any significant modulations. The gamma-ray emission from the direction of M5 is well described by an exponentially cutoff power-law spectral model, although other models cannot be excluded. The phase-averaged emission is consistent with being steady on a timescale of a few months. We estimate the number of MSPs in M5 to be between 1 and 10, using the gamma-ray conversion efficiencies for well-characterized gamma-ray MSPs in the Third Fermi-LAT Catalog of Gamma-ray Pulsars, suggesting that the sample of known MSPs in M5 is (nearly) complete, even if it is not currently possible to rule out a diffuse component of the observed gamma rays from the cluster.

Discoveries and timing of pulsars in M62

ABSTRACT Using MeerKAT, we have discovered three new millisecond pulsars (MSPs) in the bulge globular cluster M62: M62H, M62I, and M62J. All three are in binary systems, which means all ten known pulsars in the cluster are in binaries. M62H has a planetary-mass companion with a median mass Mc, med ∼ 3 MJ and a mean density of ρ ∼ 11 g cm−3. M62I has an orbital period of 0.51 d and a Mc, med ∼ 0.15 M⊙. Neither of these low-mass systems exhibit eclipses. M62J has only been detected in the two Ultra High Frequency band (816 MHz) observations with a flux density S816 = 0.08 mJy. The non-detection in the L-band (1284 MHz) indicates it has a relatively steep spectrum (β < −3.1). We also present 23-yr-long timing solutions obtained using data from the Parkes ‘Murriyang’, Effelsberg, and MeerKAT telescopes for the six previously known pulsars. For all these pulsars, we measured the second spin-period derivatives and the rate of change of orbital period caused by the gravitational field of the cluster, and their proper motions. From these measurements, we conclude that the pulsars’ maximum accelerations are consistent with the maximum cluster acceleration assuming a core-collapsed mass distribution. Studies of the eclipses of the redback M62B and the black widow M62E at four and two different frequency bands, respectively, reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies. The presence of only binary MSPs in this cluster challenges models which suggest that the MSP population of core-collapsed clusters should be dominated by isolated MSPs.

Pulse Jitter and Single-pulse Variability in Millisecond Pulsars

Understanding the jitter noise resulting from single-pulse phase and shape variations is important for the detection of gravitational waves using pulsar timing arrays. We present measurements of the jitter noise and single-pulse variability of 12 millisecond pulsars that are part of the International Pulsar Timing Array sample using the Five-hundred-meter Aperture Spherical radio Telescope. We find that the levels of jitter noise can vary dramatically among pulsars. A moderate correlation with a correlation coefficient of 0.57 between jitter noise and pulse width is detected. To mitigate jitter noise, we perform matrix template matching using all four Stokes parameters. Our results reveal a reduction in jitter noise ranging from 6.7% to 39.6%. By performing longitude-resolved fluctuation spectrum analysis, we identify periodic intensity modulations in 10 pulsars. In PSR J0030+0451, we detect single pulses with energies more than 10 times the average pulse energy, suggesting the presence of giant pulses. We also observe a periodic mode-changing phenomenon in PSR J0030+0451. We examine the achievable timing precision by selecting a subset of pulses with a specific range of peak intensity, but no significant improvement in timing precision is achievable.

Research on Pulsar Time Steered Atomic Time Algorithm Based on DPLL

In today’s society, there is a wide demand for high-precision and high-stability time service in the fields of electric power, communication, transportation and finance. At present, the time standard in various countries is mainly based on atomic clocks, but the frequency drift of atomic clocks will affect the long-term stability performance. Compared with atomic clocks, millisecond pulsars have better long-term stability and can complement with the excellent short-term stability of atomic clocks. In order to improve the long-term stability of the atomic timescale, and then improve the timing accuracy, this paper proposes an algorithm for steering the atomic clock ensemble (ACE) by ensemble pulsar time (EPT) based on digital phase locked loop (DPLL). First, the ACE and EPT are generated by the ALGOS algorithm, then the ACE is steered by EPT based on DPLL to calibrate the long-term frequency drift of the atomic clock, so that the generated steered atomic time follows both the short-term stability characteristics of ACE and the long-term stability characteristics of EPT, and finally, the steered atomic time is used to calibrate the local cesium clock. The experimental results show that the long-term stability of atomic time after steering is improved by 2 orders of magnitude compared with that before steering, and the daily drift of a local cesium clock after calibration is less than 9.47 ns in 3 yr, 3 orders of magnitude higher than that before calibration on accuracy.

Stochastic Wave Dark Matter with Fermi-LAT γ-Ray Pulsar Timing Array

Pulsar timing arrays (PTAs) can detect disturbances in the fabric of spacetime on a galactic scale by monitoring the arrival time of pulses from millisecond pulsars (MSPs). Recent advancements have enabled the use of γ-ray radiation emitted by MSPs, in addition to radio waves, for PTA experiments. Wave dark matter (DM), a prominent class of DM candidates, can be detected with PTAs due to its periodic perturbations of the spacetime metric. In response to this development, we perform in this Letter a first analysis of applying the γ-ray PTA to detect the ultralight axion-like wave DM, with the data of Fermi Large Area Telescope (Fermi-LAT). Despite its much smaller collecting area, the Fermi-LAT γ-ray PTA demonstrates a promising sensitivity potential. We show that the upper limits not far from those of the dedicated radio-PTA projects can be achieved. Moreover, we initiate a cross-correlation analysis using the data of two Fermi-LAT pulsars. The cross-correlation of phases, while carrying key information on the source of the spacetime perturbations, has been ignored in the existing data analyses for the wave DM detection with PTAs. Our analysis indicates that taking this information into account can improve the sensitivity to wave DM by ≳50% at masses below 10−23 eV.

On the Spin Period Distribution of Millisecond Pulsars

Spin period distribution provides important clues to understand the formation of millisecond pulsars (MSPs). To uncover the intrinsic period distribution, we analyze three samples of radio MSPs in the Galactic field and globular clusters. The selection bias due to pulse broadening has been corrected but turns out to be negligible. We find that all the samples can be well described by a Weibull distribution of spin frequencies. Considering MSPs in the Galactic field or globular clusters and in isolation or binary systems, we find no significant difference in the spin distribution among these subpopulations. Based on the current known population of MSPs, we find that submillisecond pulsars are unlikely to be discovered by the Square Kilometre Array, although up to ∼10 discoveries of pulsars that spin faster than the current record holder of P = 1.4 ms are expected.

Formation of PSR J1012+5307 with an Extremely Low-mass White Dwarf: Testing Magnetic Braking Models

PSR J1012+5307 is a millisecond pulsar with an extremely low-mass white dwarf (WD) companion in an orbit of 14.5 hr. Magnetic braking (MB) plays an important role in influencing the orbital evolution of binary systems with a low-mass (≲1–2 M ⊙) donor star. At present, there exist several different MB descriptions. In this paper, we investigate the formation of PSR J1012+5307 as a probe to test plausible MB models. Employing a detailed stellar evolution model with the MESA code, we find that the convection- and rotation-boosted MB and “Intermediate” MB models can reproduce the WD mass, WD radius, WD surface gravity, neutron star mass, and orbital period observed in PSR J1012+5307. However, our simulated WD has a higher effective temperature than the observations. Three other MB mechanisms, including the standard MB model, are too weak to account for the observed orbital period in the Hubble time. A long cooling timescale caused by H-shell flashes of the WD may alleviate the discrepancy between the simulated effective temperature and the observed value.

Effects of scattering in the accretion funnel on the pulse profiles of accreting millisecond pulsars

The hotspot emission of accreting millisecond pulsars (AMPs) undergoes scattering in the accretion flow between the disk inner radius and the neutron star surface. The scattering optical depth of the flow depends on the photon emission angle, which is a function of the pulse phase, and reaches its maximum when the hotspot is closest to the observer. At sufficiently large optical depths the observed pulse profile should develop a secondary minimum, the depth of which depends on the accretion rate and the emission geometry. Such a dip evolving with the accretion rate might explain the phase shift and pulse profile evolution observed in AMPs during outbursts. Accounting for scattering is important for accurate modeling of the AMP pulse profiles in order to improve the accuracy of determination of the neutron star parameters, such as their masses and radii. In this paper we present a simplified analytical model for the Thomson optical depth of the accretion funnel, and apply it to simulating the pulse profiles. We show that scattering in the accretion funnel has a significant effect on the pulse profiles at accretion rates of M ≳ 10−10 M⊙yr−1. Our model predicts a gradual evolution of the pulse profile with the accretion rate that appears to be consistent with the observations.