Search For Pulsations Research Articles

Parkes Radio and NuSTAR X-Ray Observations of the Composite Supernova Remnant B0453–685 in the Large Magellanic Cloud

Gamma-ray emission is observed coincident in position to the evolved, composite supernova remnant (SNR) B0453–685. Prior multiwavelength investigations of the region indicate that the pulsar wind nebula (PWN) within the SNR is the most likely origin for the observed gamma rays, with a possible pulsar contribution that becomes significant at energies below E ∼ 5 GeV. Constraints on the PWN hard X-ray spectrum are important for the most accurate broadband representation of PWN emission and determining the presence of a gamma-ray pulsar component. The results of Parkes radio and Nuclear Spectroscopic Telescope Array (NuSTAR) X-ray observations are presented on PWN B0453–685. We perform a search for the central pulsar in the new Parkes radio data, finding an upper limit of 12 μJy. A pulsation search in the new NuSTAR observations additionally provides a 3σ upper limit on the hard X-ray pulsed fraction of 56%. The PWN is best characterized with a photon index ΓX = 1.91 ± 0.20 in the 3–78 keV NuSTAR data, and the results are incorporated into existing broadband models. Last, we characterize a serendipitous source detected by Chandra and NuSTAR that is considered a new high-mass X-ray binary candidate.

An Image-based Search for Pulsar Candidates in the MeerKAT Bulge Survey

We report on the results of an image-based search for pulsar candidates toward the Galactic bulge. We used mosaic images from the MeerKAT radio telescope that were taken as part of a 173 deg2 survey of the bulge and Galactic center of our Galaxy at L band (856–1712 MHz) in all four Stokes I, Q, U, and V. The image rms noise levels of 12–17 μJy ba−1 represent a significant increase in sensitivity over past image-based pulsar searches. Our primary search criterion was circular polarization, but we used other criteria, including linear polarization, in-band spectral index, compactness, variability, and multiwavelength counterparts to select pulsar candidates. We first demonstrate the efficacy of this technique by searching for polarized emission from known pulsars and comparing our results with measurements from the literature. Our search resulted in a sample of 75 polarized sources. Bright stars or young stellar objects were associated with 28 of these sources, including a small sample of highly polarized dwarf stars with pulsar-like steep spectra. Comparing the properties of this sample with the known pulsars, we identified 30 compelling candidates for pulsation follow-up, including two sources with both strong circular and linear polarization. The remaining 17 sources are either pulsars or stars, but we cannot rule out an extragalactic origin or image artifacts among the brighter, flat-spectrum objects.

TRAPUM pulsar and transient search in the Sextans A and B galaxies and discovery of background FRB 20210924D

ABSTRACT The Small and Large Magellanic Clouds are the only galaxies outside our own in which radio pulsars have been discovered to date. The sensitivity of the MeerKAT radio interferometer offers an opportunity to search for a population of more distant extragalactic pulsars. The TRAPUM (TRansients And PUlsars with MeerKAT) collaboration has performed a radio-domain search for pulsars and transients in the dwarf star-forming galaxies Sextans A and B, situated at the edge of the Local Group 1.4 Mpc away. We conducted three 2-h multibeam observations at L band (856–1712 MHz) with the full array of MeerKAT. No pulsars were found down to a radio pseudo-luminosity upper limit of 7.9$\pm$0.4 Jy kpc$^{2}$ at 1400 MHz, which is 28 times more sensitive than the previous limit from the Murriyang telescope. This luminosity is 30 per cent greater than that of the brightest known radio pulsar and sets a cut-off on the luminosity distributions of the entire Sextans A and B galaxies for unobscured radio pulsars beamed in our direction. A fast radio burst was detected in one of the Sextans A observations at a dispersion measure (DM) of 737 pc cm$^{-3}$. We believe this is a background event not associated with the dwarf galaxy due to its large DM and its signal-to-noise ratio being strongest in the wide-field incoherent beam of MeerKAT.

Simulation-based Inference of Radio Millisecond Pulsars in Globular Clusters

Millisecond pulsars (MSPs) are abundant in globular clusters (GCs), which offer favorable environments for their creation. While the advent of recent, powerful facilities led to a rapid increase in MSP discoveries in GCs through pulsation searches, detection biases persist. In this work, we investigate the ability of current and future detections in GCs to constrain the parameters of the MSP population in GCs through a careful study of their luminosity function. Parameters of interest are the number of MSPs hosted by a GC, as well as the mean and the width of their luminosity function, which are typically affected by large uncertainties. While, as we show, likelihood-based studies can lead to ill-behaved posteriors on the size of the MSP population, we introduce a novel, likelihood-free analysis, based on marginal neural ratio estimation, which consistently produces well-behaved posteriors. We focus on the GC Terzan 5 (or Ter 5), which currently counts 48 detected MSPs. We find that 158 −104+294 MSPs should be hosted in this GC, but the uncertainty on this number remains large. We explore the performance of our new method on simulated Terzan 5-like data sets mimicking possible future observational outcomes. We find that significant improvement on the posteriors can be obtained by adding a reliable measurement of the diffuse radio emission of the GC to the analysis or by improving the detection threshold of current radio pulsation surveys by at least a factor of 2.

Multiwavelength identification of millisecond pulsar candidates in the Galactic bulge

Context. The existence of a population of millisecond pulsars in the Galactic bulge is supported, along with other evidence, by the Fermi GeV excess, an anomalous γ-ray emission detected almost 15 years ago in the direction of the Galactic center. However, radio surveys searching for pulsations have not yet revealed bulge millisecond pulsars. Aims. Identifying promising bulge millisecond pulsar candidates is key to motivating pointed radio pulsation searches. Candidates are often selected among steep-spectrum or polarized radio sources, but multiwavelength information can also be exploited: The aim of this work is to pinpoint strong candidates among the yet unidentified X-ray sources. Methods. We investigated the multiwavelength counterparts of sources detected by the Chandra X-ray observatory that have spectral properties expected for millisecond pulsars in the Galactic bulge. We considered that ultraviolet, optical, and strong infrared counterparts indicate that an X-ray source is not a bulge pulsar, while a radio or a faint infrared counterpart makes it a promising candidate. Results. We identify a large population of more than a thousand X-ray sources without optical, ultraviolet, or strong infrared counterparts. Among them, five are seen for the first time in unpublished radio imaging data from the Very Large Array. We provide the list of promising candidates, for most of which follow-up pulsation searches are ongoing.

Identifying Pulsar Candidates in Interferometric Radio Images Using Scintillation

Pulsars have been primarily detected by their narrow pulses or periodicity in time domain data. Interferometric surveys for pulsars are challenging due to the trade-off between beam sensitivity and beam size and the corresponding tradeoff between survey sensitivity (depth), sky coverage, and computational efforts. The detection sensitivity of time domain searches for pulsars is affected by dispersion smearing, scattering, and rapid orbital motion of pulsars in binaries. We have developed a new technique to select pulsar candidates in interferometric radio images by identifying scintillating sources and measuring their scintillation bandwidth and timescale. Identifying likely candidates allows sensitive, focused time-domain searches, saving computational effort. Pulsar scintillation is independent of its timing properties and hence offers a different selection of pulsars compared to time-domain searches. Candidates identified from this method could allow us to find hard-to-detect pulsars, such as submillisecond pulsars and pulsars in very compact, highly accelerated binary orbits. We use upgraded Giant Metrewave Radio Telescope (uGMRT) observations in the fields of PSR B1508+55, PSR J0437−4715, and PSR B0031−07 as test cases for our technique. We demonstrate that the technique correctly differentiates between the pulsar and other nonscintillating point sources and show that the extracted dynamic spectrum of the pulsar is equivalent to that extracted from the uGMRT phased array beam. We show the results from our analysis of known pulsar fields and discuss challenges in dealing with interference and instrumental effects.

Detection of Two Totally Eclipsing B-type Binaries with Extremely Low Mass Ratios

The detection of O- and B-type stars with extremely low-mass companions is very important for understanding the formation and evolution of binary stars. However, finding them remains a challenge because the low-mass components in such systems contribute such small flux to the total. During our search for pulsations among O- and B-type stars using the Transiting Exoplanet Survey Satellite (TESS) data, we found two short-period and B-type (B9) eclipsing binaries with orbital periods of 1.61613 and 2.37857 days. Photometric solutions of the two close binaries were derived by analyzing the TESS light curves with the Wilson–Devinney method. It is discovered that both of them are detached binaries with extremely low mass ratios of 0.067(2) for TIC 260342097 and 0.140(3) for TIC 209148631. The determined mass ratio indicates that TIC 260342097 is one of the lowest mass ratios among known B-type binary systems. We showed that the two systems have total eclipses with a broad and flat secondary minimum, suggesting that the photometric parameters could be derived reliably. The absolute parameters of the two binaries are estimated and it is found that the secondary components in the two systems are overluminous and oversize when compared with the normal low-mass and cool main-sequence (MS) stars. These findings may imply that the two systems are composed of a B-type MS primary and a cool pre-MS secondary with orbital periods shorter than 2.5 days. They are valuable targets to test theories of binary star formation and evolution.

Searching for pulsars, magnetars, and fast radio bursts in the sculptor galaxy using MeerKAT

Abstract The Sculptor Galaxy (NGC 253), located in the Southern Hemisphere, far off the Galactic Plane, has a relatively high star-formation rate of about 7 M⊙ yr−1 and hosts a young and bright stellar population, including several super star clusters and supernova remnants. It is also the first galaxy, apart from the Milky Way Galaxy to be associated with two giant magnetar flares. As such, it is a potential host of pulsars and/or fast radio bursts in the nearby Universe. The instantaneous sensitivity and multibeam sky coverage offered by MeerKAT therefore make it a favourable target. We searched for pulsars, radio-emitting magnetars, and fast radio bursts in NGC 253 as part of the TRAPUM large survey project with MeerKAT. We did not find any pulsars during a four-hour observation, and derive a flux density limit of 4.4 μJy at 1400 MHz, limiting the pseudo-luminosity of the brightest putative pulsar in this galaxy to 54 Jy kpc2. Assuming universality of pulsar populations between galaxies, we estimate that detecting a pulsar as bright as this limit requires NGC 253 to contain a pulsar population of ⪞20 000. We also did not detect any single pulses and our single pulse search flux density limit is 62 mJy at 1284 MHz. Our search is sensitive enough to have detected any fast radio bursts and radio emission similar to the brighter pulses seen from the magnetar SGR J1935+2154 if they had occurred during our observation.

Investigating the nature of the 2.4 hr-period eclipsing cataclysmic variable W2 in 47 Tuc

) is a cataclysmic variable (CV) in the Galactic globular cluster 47 Tucanae. Its modulation was discovered within the CATS@BAR project. The source shows all the properties of magnetic CVs, but whether it is a polar or an intermediate polar is still a matter of debate. This paper investigates the spectral and temporal properties of the source, using all archival X-ray data from and Early Data Release, to establish whether the source falls within the category of polars or intermediate polars. We fitted archival spectra with three different models: a power law, a bremsstrahlung and an optically thin thermal plasma. We also explored the temporal properties of the source with searches for pulsations with a power spectral density analysis and a Rayleigh test ($Z_n^2$). displays a mean luminosity of $ over a 20-year span, despite lower values in a few epochs. The source is not detected in the latest observation, taken with in 2022, and we infer an X-ray luminosity $ The source spectral shape does not change over time and can be equally well fitted with each of the three models, with a best-fit photon index of 1.6 for the power law and best-fit temperatures of 10 keV for both the bremsstrahlung and the thermal plasma models. We confirm the previously detected period of 8649 s, ascribed to the binary orbital period, and found a cycle-to-cycle variability associated with this periodicity. No other significant pulsation is detected. Considering the source orbital period, luminosity, spectral characteristics, long-term evolution and strong cycle-to-cycle variability, we suggest that is a magnetic CV of the polar type.

A Diamond Ring in X-Ray: Investigating the Diffuse Emission of Supernova Remnant G28.8+1.5

Supernova remnants are products of supernova explosions, which are usually observed in the X-ray band. However, the diffuse emission from the galactic supernova remnant G28.8+1.5 has been overlooked following initial searches for pulsars and pulsar wind nebulae. This paper presents background-subtracted imaging and spatial-resolved spectral analyses using archival XMM-Newton data to reveal unprecedented details of the diffuse emission, providing updated estimates of distance and age. For the first time, an iron overabundance is revealed and the progenitor is suggested to be 9–15 M ⊙. These findings refine our understanding of G28.8+1.5 and contribute to a broader comprehension of the Galactic supernova remnant family.

A VLITE Search for Millisecond Pulsars in Globular Clusters: Discovery of a Pulsar in GLIMPSE-C01

We present results from a search for pulsars in globular clusters, including the discovery of a new millisecond pulsar in the stellar cluster GLIMPSE-C01. We searched for low-frequency radio sources within 97 globular clusters using images from the VLA Low-band Ionosphere and Transient Experiment (VLITE) and epochs 1 and 2 of the VLITE Commensal Sky Survey. We discovered 10 sources in our search area, four more than expected from extragalactic source counts at our sensitivity limits. The strongest pulsar candidate was a point source found in GLIMPSE-C01 with a spectral index ∼ − 2.6, and we present additional measurements at 0.675 and 1.25 GHz from the GMRT and 1.52 GHz from the VLA that confirm the spectral index. Using archival Green Bank Telescope S-band data from 2005, we detect a binary pulsar with a spin period of 19.78 ms within the cluster. Although we cannot confirm that this pulsar is at the same position as the steep-spectrum source using the existing data, the pulse flux is consistent with the predicted flux density from other frequencies, making it a probable match. The source also shows strong X-ray emission, indicative of a higher magnetic field than most millisecond pulsars, suggesting that its recycling was interrupted. We demonstrate that low-frequency searches for steep-spectrum sources are an effective way to identify pulsar candidates, in particular on sightlines with high dispersion.

Enhancing Pulsar Candidate Identification with Self-tuning Pseudolabeling Semisupervised Learning

In the field of astronomy, machine-learning technologies are becoming increasingly crucial for identifying radio pulsars. However, the process of acquiring labeled data, which is both time-consuming and potentially biased, poses a significant limitation to current methodologies. In response to these challenges, this study proposes and validates a self-tuning pseudolabeling semisupervised learning approach. This approach synthesizes a vast unlabeled data set with a considerably smaller set of labeled data, markedly enhancing classifier performance and effectuating a transition from traditional fully supervised learning methods to more efficient radio pulsar detection strategies. Our experimental outcomes demonstrate that even with a training set comprised of only 100 labeled pulsar candidates, this method can attain a recall rate of 92.35% and an F1 score of 93.89%. When the number of labeled examples is increased to 800, we observe a further improvement in performance, with the recall rate rising to 97.50% and the F1 score reaching 97.16%. The utility of the semisupervised learning approach is evident even with minimal labeled data, which is a common scenario in the search for pulsars, including in environments like globular clusters. What stands out is the method’s capacity to detect pulsar candidates effectively with only a limited number of labeled examples. This emphasizes the robust potential of our approach to facilitate early-stage pulsar surveys and highlights its capability to yield substantial results even when labeled data are in short supply.

Short-period pulsating hot subdwarf stars observed by TESS

We present the results of an extension of our Transiting Exoplanet Survey Satellite (TESS) search for short-period pulsations in compact stellar objects observed during the second and fourth years of the TESS mission, which targeted the northern ecliptic hemisphere. For many of the targets, we exploited unpublished spectroscopic data to confirm or re-evaluate the object’s spectral classification. From the TESS photometry, we identified 50 short-period hot-subdwarf pulsators, including 35 sdB and 15 sdOB stars. The sample contains 26 pulsators that were unknown prior to the TESS mission. Nine stars show signals at both low and high frequencies and have been categorized as “hybrid” pulsators. For each pulsator, we report the list of prewhitened frequencies, along with and their amplitude spectra calculated from the TESS data. We attempt to identify possible multiplets caused by stellar rotation and we report five candidates with rotation periods between 11 and 46 d. With the search for p-mode pulsating hot subdwarfs in TESS Sectors 1–60 complete, we discuss the completeness of the study, as well as the instability strip and the evolutionary status of the stars we found. We also compare the distribution of pulsation periods as a function of effective temperature and surface gravity with theoretical predictions. We find that the percentage of undetected pulsators in the TESS mission increases with decreasing brightness measurements of stars, reaching 25% near the 15th magnitude. When comparing the distribution of hot subdwarfs in the log g − Teff plane with stellar models, we underline the importance of a proper treatment of the hydrogen-rich envelope composition (strongly affected by microscopic diffusion processes). We also emphasize that the stellar mass is a significant factor in understanding the instability strip. The p-mode instability strip is confirmed to be narrower than predicted by prior non-adiabatic calculations based on models incorporating equilibrium between gravitational settling and radiative levitation for iron. This implies that competing mixing processes ignored in these models must play a role in reducing the amount of levitating iron in the stellar envelope. Interestingly, we find that the coolest p-mode pulsators with Teff ≲ 30 000 K (including the hybrid ones) tend to cluster around the terminal age of the extreme horizontal branch of canonical mass (TAEHB at ∼0.47 M⊙). This trend is expected from the non-adiabatic pulsation calculations. Otherwise, the overall pulsation period distributions tend to reproduce the predicted trends in Teff and log g.

TRAPUM search for pulsars in supernova remnants and pulsar wind nebulae − I. Survey description and initial discoveries

ABSTRACT We present the description and initial results of the TRAPUM (TRAnsients And PUlsars with MeerKAT) search for pulsars associated with supernova remnants (SNRs), pulsar wind nebulae, and unidentified TeV emission. The list of sources to be targeted includes a large number of well-known candidate pulsar locations but also new candidate SNRs identified using a range of criteria. Using the 64-dish MeerKAT radio telescope, we use an interferometric beamforming technique to tile the potential pulsar locations with coherent beams which we search for radio pulsations, above a signal-to-noise of 9, down to an average flux density upper limit of 30 μJy. This limit is target-dependent due to the contribution of the sky and nebula to the system temperature. Coherent beams are arranged to overlap at their 50 per cent power radius, so the sensitivity to pulsars is not degraded by more than this amount, though realistically averages around 65 per cent if every location in the beam is considered. We report the discovery of two new pulsars; PSR J1831−0941 is an adolescent pulsar likely to be the plerionic engine of the candidate PWN G20.0+0.0, and PSR J1818−1502 appears to be an old and faint pulsar that we serendipitously discovered near the centre of a SNR already hosting a compact central object. The survey holds importance for better understanding of neutron star birth rates and the energetics of young pulsars.

A Millisecond Pulsar Binary Embedded in a Galactic Center Radio Filament

The Galactic center is host to a population of extraordinary radio filaments, thin linear structures that trace out magnetic field lines running perpendicular to the Galactic plane. Using Murriyang, the 64 m Parkes radio telescope, we conducted a search for pulsars centered on the position of a compact source in the filament G359.0−0.2. We discovered a millisecond pulsar (MSP), PSR J1744−2946, with a period P = 8.4 ms, that is bound in a 4.8 hr circular orbit around an M c > 0.05 M ⊙ companion. The pulsar dispersion measure of 673.7 ± 0.1 pc cm−3 and Faraday rotation measure of 3011 ± 3 rad m−2 are the largest of any known MSP. Its radio pulses are moderately scattered due to multipath propagation through the interstellar medium, with a scattering timescale of 0.87 ± 0.08 ms at 2.6 GHz. Using MeerKAT, we localized the pulsar to a point source embedded in a low-luminosity radio filament, the “Sunfish”, that is unrelated to G359.0−0.2. Our discovery of the first MSP within 1° of the Galactic center hints at a large population of these objects detectable via high-frequency surveys. The association with a filament points to pulsars as the energy source responsible for illuminating the Galactic center radio filaments.

A Needle in a Cosmic Haystack: A Review of FRB Search Techniques

Ephemeral Fast Radio Bursts (FRBs) must be powered by some of the most energetic processes in the Universe. That makes them highly interesting in their own right, and as precise probes for estimating cosmological parameters. This field thus poses a unique challenge: FRBs must be detected promptly and immediately localised and studied based only on that single millisecond-duration flash. The problem is that the burst occurrence is highly unpredictable and that their distance strongly suppresses their brightness. Since the discovery of FRBs in single-dish archival data in 2007, detection software has evolved tremendously. Pipelines now detect bursts in real time within a matter of seconds, operate on interferometers, buffer high-time and frequency resolution data, and issue real-time alerts to other observatories for rapid multi-wavelength follow-up. In this paper, we review the components that comprise a FRB search software pipeline, we discuss the proven techniques that were adopted from pulsar searches, we highlight newer, more efficient techniques for detecting FRBs, and we conclude by discussing the proposed novel future methodologies that may power the search for FRBs in the era of big data astronomy.

Study on tiered storage algorithm based on heat correlation of astronomical data

With the surge in astronomical data volume, modern astronomical research faces significant challenges in data storage, processing, and access. The I/O bottleneck issue in astronomical data processing is particularly prominent, limiting the efficiency of data processing. To address this issue, this paper proposes a tiered storage algorithm based on the access characteristics of astronomical data. The C4.5 decision tree algorithm is employed as the foundation to implement an astronomical data access correlation algorithm. Additionally, a data copy migration strategy is designed based on tiered storage technology to achieve efficient data access. Preprocessing tests were conducted on 418GB NSRT (Nanshan Radio Telescope) formaldehyde spectral line data, showcasing that tiered storage can potentially reduce data processing time by up to 38.15%. Similarly, utilizing 802.2 GB data from FAST (Five-hundred-meter Aperture Spherical radio Telescope) observations for pulsar search data processing tests, the tiered storage approach demonstrated a maximum reduction of 29.00% in data processing time. In concurrent testing of data processing workflows, the proposed astronomical data heat correlation algorithm in this paper achieved an average reduction of 17.78% in data processing time compared to centralized storage. Furthermore, in comparison to traditional heat algorithms, it reduced data processing time by 5.15%. The effectiveness of the proposed algorithm is positively correlated with the associativity between the algorithm and the processed data. The tiered storage algorithm based on the characteristics of astronomical data proposed in this paper is poised to provide algorithmic references for large-scale data processing in the field of astronomy in the future.

Applying Hybrid Clustering in Pulsar Candidate Sifting with Multi-modality for FAST Survey

Pulsar search is always the basis of pulsar navigation, gravitational wave detection and other research topics. Currently, the volume of pulsar candidates collected by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) shows an explosive growth rate that has brought challenges for its pulsar candidate filtering system. Particularly, the multi-view heterogeneous data and class imbalance between true pulsars and non-pulsar candidates have negative effects on traditional single-modal supervised classification methods. In this study, a multi-modal and semi-supervised learning based on a pulsar candidate sifting algorithm is presented, which adopts a hybrid ensemble clustering scheme of density-based and partition-based methods combined with a feature-level fusion strategy for input data and a data partition strategy for parallelization. Experiments on both High Time Resolution Universe Survey II (HTRU2) and actual FAST observation data demonstrate that the proposed algorithm could excellently identify pulsars: On HTRU2, the precision and recall rates of its parallel mode reach 0.981 and 0.988 respectively. On FAST data, those of its parallel mode reach 0.891 and 0.961, meanwhile, the running time also significantly decreases with the increment of parallel nodes within limits. Thus, we can conclude that our algorithm could be a feasible idea for large scale pulsar candidate sifting for FAST drift scan observation.

The NANOGrav 12.5 yr Data Set: A Computationally Efficient Eccentric Binary Search Pipeline and Constraints on an Eccentric Supermassive Binary Candidate in 3C 66B

The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05 yr period and low redshift (∼0.02) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using pulsar timing array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it to simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5 yr data set employing PTA signal models containing Earth term-only as well as Earth+pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity e 0 < 0.5 and the symmetric mass ratio η > 0.1. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of 88.1 ± 3.7 ns for the Earth term-only and 81.74 ± 0.86 ns for the Earth+pulsar term searches for e 0 < 0.5 and η > 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 109 and (1.81 ± 0.01) × 109 M ☉. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5 yr data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations.

The Green Bank North Celestial Cap Survey. IX. Timing Follow-up for 128 Pulsars

The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large data sets of archival observations and connect them to recent, high-cadence observations taken using the Canadian Hydrogen Intensity Mapping Experiment telescope. This pipeline handles data for 128 pulsars and has produced measurements of spin, positional, and orbital parameters that connect data over observation gaps as large as 2000 days. We have also measured glitches in the timing residuals for five of the pulsars included and proper motion for 19 sources (13 new). We include updates to orbital parameters for 19 pulsars, including nine previously unpublished binaries. For two of these binaries, we provide updated measurements of post-Keplerian binary parameters, which result in much more precise estimates of the total masses of both systems. For PSR J0509+3801, the much improved measurement of the Einstein delay yields much improved mass measurements for the pulsar and its companion, 1.399(6) M ⊙ and 1.412(6) M ⊙, respectively. For this system, we have also obtained a measurement of the orbital decay due to the emission of gravitational waves, ṖB=−1.37(7)×10−12 , which is in agreement with the rate predicted by general relativity for these masses.