Population Of Millisecond Pulsars Research Articles

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 days and a Mc, med ∼ 0.15 M⊙. Neither of these low-mass systems exhibit eclipses. M62J has only been detected in the two UHF 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.

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.

A Fermi-LAT Study of Globular Cluster Dynamical Evolution in the Milky Way: Millisecond Pulsars as the Probe

Using archival Fermi-LAT data with a time span of ∼12 yr, we study the population of Millisecond Pulsars (MSPs) in Globular Clusters (GlCs) and investigate their dependence on cluster dynamical evolution in the Milky Way. We show that the γ-ray luminosity (L γ ) and emissivity (i.e., ϵ γ = L γ /M, with M the cluster mass) are good indicators of the population and abundance of MSPs in GlCs, and they are highly dependent on the dynamical evolution history of the host clusters. Specifically speaking, the dynamically older GlCs with more compact structures are more likely to have larger L γ and ϵ γ , and these trends can be summarized as strong correlations with cluster stellar encounter rate Γ and the specific encounter rate (Λ = Γ/M), with L γ ∝ Γ0.70±0.11 and ϵ γ ∝ Λ0.73±0.13 for dynamically normal GlCs. However, as GlCs evolve into deep core collapse, these trends are found to be reversed, implying that strong encounters may have lead to the disruption of Low-Mass X-ray Binaries and ejection of MSPs from core-collapsed systems. Besides, the GlCs are found to exhibit larger ϵ γ with increasing stellar mass function slope (ϵ γ ∝ 10(0.57±0.1)α ), decreasing tidal radius ( ϵγ∝Rt−1.0±0.22 ) and distances from the Galactic Center (GC, ϵγ∝Rgc−1.13±0.21 ). These correlations indicate that, as GlCs losing kinetic energy and spiral in toward the GC, tidal stripping and mass segregation have a preference in leading to the loss of normal stars from GlCs, while MSPs are more likely to concentrate to cluster center and be deposited into the GC. Moreover, we gauge ϵ γ of GlCs is ∼10–1000 times larger than the Galactic bulge, the latter is thought to reside thousands of unresolved MSPs and may be responsible for the GC γ-ray excess, which supports that GlCs are generous contributors to the population of MSPs in the GC.

Probing the Pulsar Explanation of the Galactic-Center GeV Excess Using Continuous Gravitational-Wave Searches.

Over 10years ago, Fermi observed an excess of GeV gamma rays from the Galactic Center whose origin is still under debate. One explanation for this excess involves annihilating dark matter, another requires an unresolved population of millisecond pulsars concentrated at the Galactic Center. In this work, we use the results from LIGO and Virgo's most recent all-sky search for quasimonochromatic, persistent gravitational-wave signals from isolated neutron stars, which is estimated to be about 20%-50% of the population, to determine whether unresolved millisecond pulsars could actually explain this excess. First, we choose a luminosity function that determines the number of millisecond pulsars required to explain the observed excess. Then, we consider two models for deformations on millisecond pulsars to determine their ellipticity distributions, which are directly related to their gravitational-wave radiation. Lastly, based on null results from the O3 frequency-Hough all-sky search for continuous gravitational waves, we find that a large set of the parameter space in the pulsar luminosity function can be excluded. We also evaluate how these exclusion regions may change with respect to various model choices. Our results are the first of their kind and represent a bridge between gamma-ray astrophysics, gravitational-wave astronomy, and dark-matter physics.

Improved Galactic diffuse emission model strengthens the case for a Millisecond Pulsar explanation of the Fermi GeV excess

After more than a decade since its discovery, the Galactic center gamma-ray excess - discovered with the Fermi Large Area Telescope - remains puzzling. While the spectrum of the signal can be explained by either dark matter or an unresolved population of millisecond pulsars, the spatial morphology of this excess seems to hold the key to separate the two theories. In this contribution, we present the results of a recent study in which we use bleeding edge models for interstellar gas, inverse Compton emission, and stellar mass models to reanalyze the Galactic center excess. We find that the spatial morphology of the excess is highly correlated with stellar matter in the Galactic bulge, providing strong support for the millisecond pulsar hypothesis.

A new pulsar candidate in 47 Tucanae discovered with MeerKAT imaging

ABSTRACT MeerKAT imaging of the globular cluster 47 Tucanae (47 Tuc) reveals 1.28 GHz continuum emission at the locations of 20 known millisecond pulsars (MSPs). We use time series and spectral imaging to investigate the image-domain characteristics of the MSPs, and search for previously unknown sources of interest. The MSPs exhibit a range of differences in their temporal and spectral properties compared the general background radio source population. Temporal variability differs strongly from pulsar to pulsar, some appearing to vary randomly on 15 min time-scales, others varying coherently by factors of >10 on time-scales of hours. The error in the typical power law fit to the spectrum emerges as a powerful parameter for indentifying the MSPs. This behaviour is likely due to differing diffractive scintillation conditions along the sightlines to the MSPs. One MSP exhibits tentative periodic variations that are consistent with modulation due the orbit of an eclipsing binary system. One radio source has spectro-temporal properites closely resembling those of the MSP population in the cluster, and we report its position as a candidate new MSP, or alternatively an interferometric localization of one of six MSPs which do not yet have an accurate position from the timing solutions.

The status of the galactic center gamma-ray excess

The Galactic Center Gamma-Ray Excess has a spectrum, angular distribution, and overall intensity that agree remarkably well with that expected from annihilating dark matter particles in the form of a m_X \simmX∼ 50 GeV thermal relic. Previous claims that these photons are clustered on small angular scales or trace the distribution of known stellar populations once appeared to favor interpretations in which this signal originates from a large population of unresolved millisecond pulsars. More recent work, however, has overturned these conclusions, finding that the observed gamma-ray excess does not contain discernible small scale power, and is distributed with approximate spherical symmetry, not tracing any known stellar populations. In light of these results, it now appears significantly more likely that the Galactic Center Gamma-Ray Excess is produced by annihilating dark matter.

Discovery and Timing Analysis of New Pulsars in Globular Cluster NGC 5024: New Observations from FAST

The most distant globular cluster (GC) with known pulsars is NGC 5024 (M53). In this paper, we report the discovery of a new binary millisecond pulsar, PSR J1312+1810E (M53E), and present the new timing solutions for M53B–E based on 22 observations from the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These discoveries and timing work benefit from FAST’s high sensitivity. We find that M53C is the only isolated millisecond pulsar known in this distant GC, with a spin period of 12.53 ms and spin period derivative of 5.26 × 10−20 s s−1. Our results reveal orbital periods of 47.7, 5.8, and 2.4 days for M53B, D, and E, respectively. The companions, with masses of 0.25, 0.27, and 0.18 M ⊙, respectively, are likely to be white dwarf stars; if they are extended objects, they do not eclipse the pulsars. We find no X-ray counterparts for these millisecond pulsars in archival Chandra images in the band of 0.3–8 keV. The characteristics of this pulsar population are similar to the population of millisecond pulsars in the Galactic disk, as expected from the low stellar density of M53.

A binary pulsar in a 53-minute orbit

Spider pulsars are neutron stars that have a companion star in a close orbit. The companion star sheds material to the neutron star, spinning it up to millisecond rotation periods, while the orbit shortens to hours. The companion is eventually ablated and destroyed by the pulsar wind and radiation1,2. Spider pulsars are key for studying the evolutionary link between accreting X-ray pulsars and isolated millisecond pulsars, pulsar irradiation effects and the birth of massive neutron stars3–6. Black widow pulsars in extremely compact orbits (as short as 62 minutes7) have companions with masses much smaller than 0.1 M⊙. They may have evolved from redback pulsars with companion masses of about 0.1–0.4 M⊙ and orbital periods of less than 1 day8. If this is true, then there should be a population of millisecond pulsars with moderate-mass companions and very short orbital periods9, but, hitherto, no such system was known. Here we report radio observations of the binary millisecond pulsar PSR J1953+1844 (M71E) that show it to have an orbital period of 53.3 minutes and a companion with a mass of around 0.07 M⊙. It is a faint X-ray source and located 2.5 arcminutes from the centre of the globular cluster M71.

Mind the gap: the discrepancy between simulation and reality drives interpretations of the Galactic Center Excess

The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model (γEM) used in Galactic center analyses affects the interpretation of the GCE's nature.To address this issue, we construct an ultra-fast and powerful inference pipeline based on convolutional Deep Ensemble Networks. We explore the two main competing hypotheses for the GCE using a set of γEMs with increasing parametric freedom. We calculate the fractional contribution (f src) of a dim population of MSPs to the total luminosity of the GCE and analyze its dependence on the complexity of the γEM. For the simplest γEM, we obtain f src = 0.10 ± 0.07, while the most complex model yields f src = 0.79 ± 0.24.In conclusion, we find that the statement about the nature of the GCE (dark matter or not) strongly depends on the assumed γEM.The quoted results for f src do not account for the additional uncertainty arising from the fact that the observed gamma-ray sky is out-of-distribution concerning the investigated γEM iterations. We quantify the reality gap between our γEMs using deep-learning-based One-Class Deep Support Vector Data Description networks, revealing that all employed γEMs have gaps to reality.Our study casts doubt on the validity of previous conclusions regarding the GCE and dark matter, and underscores the urgent need to account for the reality gap and consider previously overlooked “out of domain” uncertainties in future interpretations.

Multi-messenger Approaches to Supermassive Black Hole Binary Detection and Parameter Estimation. II. Optimal Strategies for a Pulsar Timing Array

Pulsar timing arrays (PTAs) are Galactic-scale gravitational wave (GW) detectors consisting of precisely timed pulsars distributed across the sky. Within the decade, PTAs are expected to detect nanohertz GWs emitted by close-separation supermassive black hole binaries (SMBHBs), thereby opening up the low-frequency end of the GW spectrum for science. Individual SMBHBs which power active galactic nuclei are also promising multi-messenger sources; they may be identified via theoretically predicted electromagnetic (EM) signatures and be followed up by PTAs for GW observations. In this work, we study the detection and parameter estimation prospects of a PTA which targets EM-selected SMBHBs. Adopting a simulated Galactic millisecond pulsar population, we envisage three different pulsar timing campaigns which observe three mock sources at different sky locations. We find that an all-sky PTA which times the best pulsars is an optimal and feasible approach to observe EM-selected SMBHBs and measure their source parameters to high precision (i.e., comparable to or better than conventional EM measurements). We discuss the implications of our findings in the context of future PTA experiments with the planned Deep Synoptic Array-2000 and the multi-messenger studies of SMBHBs such as the well-known binary candidate OJ 287.

A Comparison of Millisecond Pulsar Populations between Globular Clusters and the Galactic Field

We have performed a systematic study of the rotational, orbital, and X-ray properties of millisecond pulsars (MSPs) in globular clusters (GCs) and compared their nature with those of the MSPs in the Galactic field (GF). We found that GC MSPs generally rotate slower than their counterparts in the GF. Different from the expectation of a simple recycling scenario, no evidence for the correlation between the orbital period and the rotation period can be found in the MSP binaries in GCs. There is also an indication that the surface magnetic field of GC MSPs is stronger than those in the GF. All these suggest dynamical interactions in GCs can alter the evolution of MSPs/their progenitors, which can leave an imprint on their X-ray emission properties. While the MSPs in both GF and GCs have similar distributions of X-ray luminosity and hardness, our sample supports the notion that these two populations follow different relations between the X-ray luminosity and spin-down power. We discuss this in terms of both the pulsar emission model and the observational bias.

Bayesian parameter estimation for targeted anisotropic gravitational-wave background

Extended sources of the stochastic gravitational backgrounds have been conventionally searched on the spherical harmonics bases. The analysis during the previous observing runs by the ground-based gravitational-wave detectors, such as LIGO and Virgo, have yielded the constraints on the angular power spectrum ${C}_{\ensuremath{\ell}}$, yet it lacks the capability of estimating other parameters such as a spectral index. In this paper, we introduce an alternative Bayesian formalism to search for such stochastic signals with a particular distribution of anisotropies on the sky. This approach provides a Bayesian posterior of model parameters and also enables selection tests among different signal models. While the conventional analysis fixes the highest angular scale a priori, here we show a more systematic and quantitative way to determine the cutoff scale based on a Bayes factor, which depends on the amplitude and the angular scale of observed signals. Also, we analyze the third observing runs of LIGO and Virgo for the population of millisecond pulsars and obtain the 95% constraints of the signal amplitude, $\ensuremath{\epsilon}<2.7\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}8}$.

Tracking millisecond pulsars responsible for the Fermi GeV excess

More than 10 years ago, an excess of γ-ray photons coming from the Galactic center was discovered in the Fermi-LAT data. First attributed to dark matter, it has since been shown that it should have at least a partial stellar origin. One hypothesis is the presence of a population of millisecond pulsars (MSPs) confined in the Galactic bulge. We here present our recent progress in the selection of MSP candidates.

Millisecond Pulsars in Dense Star Clusters: Evolution, Scaling Relations, and the Galactic-center Gamma-Ray Excess

The number of millisecond pulsars (MSPs) observed in Milky Way globular clusters has increased explosively in recent years, but the underlying population is still uncertain due to observational biases. We use state-of-the-art N-body simulations to study the evolution of MSP populations in dense star clusters. These cluster models span a wide range in initial conditions, including different initial masses, metallicities, and virial radii, which nearly cover the full range of properties exhibited by the population of globular clusters in the Milky Way. We demonstrate how different initial cluster properties affect the number of MSPs, for which we provide scaling relations as a function of cluster age and mass. As an application, we use our formulae to estimate the number of MSPs delivered to the Galactic center from inspiralling globular clusters to probe the origin of the Galactic-center gamma-ray excess detected by Fermi. We predict about 400 MSPs in the Galactic center from disrupted globular clusters, which can potentially explain most of the observed gamma-ray excess.

The Andromeda gamma-ray excess: background systematics of the millisecond pulsars and dark matter interpretations

ABSTRACT Since the discovery of an excess in gamma rays in the direction of M31, its cause has been unclear. Published interpretations focus on dark matter or stellar related origins. Studies of a similar excess in the Milky Way centre motivate a correlation of the spatial morphology of the signal with the distribution of stellar mass in M31. However, a robust determination of the best theory for the observed excess emission is challenging due to uncertainties in the astrophysical gamma-ray foreground model. We perform a spectro-morphological analysis of the M31 gamma-ray excess using state-of-the-art templates for the distribution of stellar mass in M31 and novel astrophysical foreground models for its sky region. We construct maps for the old stellar populations of M31 based on data from the PAndAS survey and carefully remove the foreground stars. We also produce improved astrophysical foreground models via novel image inpainting techniques based on machine learning methods. Our stellar maps, mimicking the location of a population of millisecond pulsars in the bulge of M31, reach a 5.4σ significance, making them as strongly favoured as the simple phenomenological models usually considered in the literature, e.g. disc-like templates. This detection is robust to generous variations of the astrophysical foreground model. Once the stellar templates are included in the astrophysical model, we show that the dark matter annihilation interpretation of the signal is unwarranted. We demonstrate that about one million unresolved millisecond pulsars naturally explain the observed gamma-ray luminosity per stellar mass, energy spectrum, and stellar bulge-to-disc flux ratio.

Targeted search for the stochastic gravitational-wave background from the galactic millisecond pulsar population

The millisecond pulsars, old-recycled objects spinning with high frequency $\mathcal{O}$ (kHz) sustaining the deformation from their spherical shape, may emit gravitational waves (GW). These are one of the potential candidates contributing to the anisotropic stochastic gravitational-wave background (SGWB) observable in the ground-based GW detectors. Here, we present the results from a likelihood-based targeted search for the SGWB due to millisecond pulsars in the Milky Way, by analyzing the data from the first three observing runs of Advanced LIGO and Advanced Virgo detector. We assume that the shape of SGWB power spectra and the sky distribution is known a priori from the population synthesis model. The information of the ensemble source properties, i.e., the in-band number of pulsars, ${N}_{\mathrm{obs}}$ and the averaged ellipticity, ${\ensuremath{\mu}}_{\ensuremath{\epsilon}}$ is encoded in the maximum likelihood statistic. We do not find significant evidence for the SGWB signal from the considered source population. The best Bayesian upper limit with 95% confidence for the parameters are ${N}_{\mathrm{obs}}\ensuremath{\le}8.8\ifmmode\times\else\texttimes\fi{}{10}^{4}$ and ${\ensuremath{\mu}}_{ϵ}\ensuremath{\le}1.4\ifmmode\times\else\texttimes\fi{}1{0}^{\ensuremath{-}6}$, which is comparable to the bounds on mean ellipticity with the GW observations of the individual pulsars. Finally, we show that for the plausible case of ${N}_{\mathrm{obs}}=40000$, with the one year of observations, the one-sigma sensitivity on ${\ensuremath{\mu}}_{\ensuremath{\epsilon}}$ might reach $1.5\ifmmode\times\else\texttimes\fi{}1{0}^{\ensuremath{-}7}$ and $4.1\ifmmode\times\else\texttimes\fi{}1{0}^{\ensuremath{-}8}$ for the second-generation detector network having $\mathrm{A}+$ sensitivity and third-generation detector network, respectively.

Observational Evidence for a Spin-up Line in the P– Diagram of Millisecond Pulsars

It is believed that millisecond pulsars attain their fast spins by accreting matter and angular momentum from companion stars. Theoretical modeling of the accretion process suggests a spin-up line in the period–period derivative (P-) diagram of millisecond pulsars, which plays an important role in population studies of radio millisecond pulsars and accreting neutron stars in X-ray binaries. Here we present observational evidence for such a spin-up line using a sample of 143 radio pulsars with P < 30 ms. We also find that PSRs J1823−3021A and J1824−2452A, located near the classic spin-up line, are consistent with the broad population of millisecond pulsars. Finally, we show that our approach of Bayesian inference can probe accretion physics, allowing constraints to be placed on the accretion rate and the disk–magnetosphere interaction.

Luminosity functions consistent with a pulsar-dominated Galactic Center excess

A new population of millisecond pulsars is a long-standing proposed explanation for the excess of GeV-scale gamma rays emanating from the region surrounding the center of the Milky Way (the “Galactic Center excess”). We examine several simple parameterizations of possible luminosity functions for this population, as well as several benchmark luminosity functions proposed in the literature, and compare the predicted populations of resolved point sources to the Fermi 4FGL-DR2 point source catalog and a sub-population recently identified using wavelet-based methods. We provide general results that can be used to translate upper limits on the number of resolved point sources associated with the excess, and the fraction of the flux in the excess that can be attributed to resolved sources, into limits on the luminosity function parameter space. We discuss a number of important systematic uncertainties, including in the detection threshold model and the total flux attributed to the excess. We delineate regions of parameter space (containing existing benchmark models) where there is no apparent tension with current data, and the number of total pulsars needed to explain the excess is in the range of \U0001d4aa(104-5). In the future, lowered point source detection thresholds could be achieved either by new analysis methods or new data. An order-of-magnitude reduction in the sensitivity threshold (which may already be achieved by novel analyses probing sub-threshold source populations) could hope to resolve more than 30% of the flux of the excess even in pessimistic scenarios.

Stellar core-merger-induced collapse: new formation pathways for black holes, Thorne–Żytkow objects, magnetars, and superluminous supernovae

ABSTRACT Most neutron stars (NSs) and black holes (BHs) are believed to be the final remnants in the evolution of massive stars. In this study, we propose a new formation channel for the formation of BHs and peculiar NSs [specifically, magnetars and Thorne–Żytkow objects (T$\dot{\rm Z}$Os)], which we refer to as the core-merger-induced collapse (CMIC) model. This model involves the merger during a common-envelope phase of an oxygen/neon/magnesium composition white dwarf and the core of a hydrogen-rich or helium-rich non-degenerate star, leading to the creation of peculiar new types of objects. The results of binary population synthesis simulations show that the CMIC channel could make important contributions to the populations of (millisecond) pulsars, T$\dot{\rm Z}$Os, magnetars, and BHs. The possibility of superluminous supernovae powered by T$\dot{\rm Z}$Os, magnetars, and BHs formed through the CMIC model is also being investigated. Magnetars with immediate matter surroundings formed after the CMIC might be good sources for fast radio bursts.