Search for: All records

ABSTRACT We present the deepest and highest-resolution radio continuum imaging of the Galactic globular cluster Terzan 5, one of the most crowded locations in the radio sky. In these new 2–4 GHz Karl G. Jansky Very Large Array images, we detect 38 of the 49 confirmed pulsars, including extensive multi-frequency eclipse mapping of the luminous redback Ter5A. Nonetheless, there is still a large amount of diffuse residual flux from pulsars that are fainter than our 2.5 GHz continuum detection limit of $$\sim 11\, \mu$$Jy. Using a range of approaches including image-based simulations, we model the fluxes of the detected pulsars together with the residual flux. We find a minimum total population of $$N\sim 250$$ detectable pulsars in Terzan 5 and perhaps substantially more, though the luminosity function remains very uncertain. Consideration of the $$\gamma$$-ray properties of the cluster, though also not unambiguous to interpret, leads to consistent conclusions. These pulsar population estimates are larger than inferred from previous work and highlight Terzan 5 as a keystone target for next-generation radio facilities. 

Abstract We present 15 yr of Nançay and Green Bank radio telescope timing observations for PSR J1231−1411. This millisecond pulsar is a primary science target for the Neutron Star Interior Composition Explorer telescope (NICER, which discovered its X-ray pulsations), has accumulated near-continuousγ-ray data since the Fermi-Large Area Telescope’s launch, and has been studied extensively with the Green Bank and Nançay radio telescopes. We have undertaken a campaign with the Green Bank Telescope targeting specific orbital phases designed to improve our constraint on the pulsar’s mass through the detection of a relativistic Shapiro delay. Both frequentist and Bayesian techniques—the latter incorporating priors from white dwarf binary evolution models—are applied to 15 yr of radio observations, yielding relatively weak constraints on the companion and pulsar masses of $0.23_{-0.06}^{+0.09}$M_⊙and $1.87_{-0.67}^{+1.11}$M_⊙, respectively (68.3% CI from Bayesian fits); however, the orbital inclination is measured to better relative precision ( $79.80_{-4.70}^{+3.47}$°). Restricting the maximum allowed pulsar mass to 3M_⊙improves the constraint and lowers the measured mass to $1.71_{-0.56}^{+0.70}$M_⊙. A fully generalized Bayesian fit that simultaneously samples the noise and timing models yields a pulsar mass in close agreement with this value. While our radio-derived inclination result has informed recent NICER X-ray studies of J1231−1411, the lessons learned from this troublesome pulsar will also bolster future high-precision mass measurement campaigns and resulting constraints on the neutron star interior equation of state. 

Abstract We present timing solutions spanning nearly two decades for five redback (RB) systems found in globular clusters (GCs), created using a novel technique that effectively “isolates” the pulsar. By accurately measuring the time of passage through periastron (T₀) at points over the timing baseline, we use a piecewise-continuous binary model to get local solutions of the orbital variations that we pair with long-term orbital information to remove the orbital timing delays. The isolated pulse times of arrival can then be fit to describe the spin behavior of the millisecond pulsar (MSP). The results of our timing analyses via this method are consistent with those of conventional timing methods for binaries in GCs as demonstrated by analyses of NGC 6440D. We also investigate the observed orbital phase variations for these systems. Quasiperiodic oscillations in Terzan 5P’s orbit may be the result of changes to the gravitational quadruple moment of the companion as prescribed by the Applegate model. We find a striking correlation between the standard deviation of the phase variations as a fraction of a system’s orbit ( ${\sigma }_{\mathrm{\Delta }{T}_0}$) and the MSP’s spin frequency, as well as a potential correlation between ${\sigma }_{\mathrm{\Delta }{T}_0}$and the binary’s projected semimajor axis. While long-term RB timing is fraught with large systematics, our work provides a needed alternative for studying systems with significant orbital variations, especially when high-cadence monitoring observations are unavailable. 

Abstract The Fourier domain acceleration search (FDAS) and Fourier domain jerk search (FDJS) are proven matched-filtering techniques for detecting binary pulsar signatures in time-domain radio astronomy data sets. Next-generation radio telescopes such as the SPOTLIGHT project at the Giant Metrewave Radio Telescope (GMRT) produce data at rates that mandate real-time processing, as storage of the entire captured data set for subsequent offline processing is infeasible. The computational demands of FDAS and FDJS make them challenging to implement in real-time detection pipelines, requiring costly high-performance computing facilities. To address this, we propose Pulscan, an unmatched-filtering approach that achieves order-of-magnitude improvements in runtime performance compared to FDAS while being able to detect both accelerated and some jerked binary pulsars. We profile the sensitivity of Pulscan using a distribution (N= 10,955) of synthetic binary pulsars (simulated post–radio-frequency interference mitigation) and compare its performance with FDAS and FDJS. Our implementation of Pulscan includes an OpenMP version for multicore CPU acceleration, a version for heterogeneous CPU/GPU environments such as NVIDIA Grace Hopper, and a fully optimized NVIDIA GPU implementation for integration into an AstroAccelerate pipeline, which will be deployed in the SPOTLIGHT project at the GMRT. Our results demonstrate that unmatched filtering in Pulscan can serve as an efficient data reduction step prior to FDAS or FDJS, selecting data sets for further analysis and focusing subsequent computational resources on likely binary pulsar signatures. 

Abstract The existence of compact stellar remnants in the mass range 2–5M_⊙has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5M_⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte CarloN-body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc⁻³yr⁻¹. 

Abstract We report the discovery of CHIME J1634+44, a long-period radio transient (LPT) unique for two aspects: it is the first known LPT to emit fully circularly polarized radio bursts, and it is the first LPT with a significant spin-up. Given that high circular polarization (>90%) has been observed in FRB 20201124A and in some giant pulses of PSR B1937+21, we discuss the implications of the high circular polarization of CHIME J1634+44 and conclude its emission mechanism is likely to be “pulsar-like.” While CHIME J1634+44 has a pulse period of 841 s, its burst arrival patterns are indicative of a secondary 4206 s period, probably associated with binary activity. The timing properties suggest it has a significantly negative period derivative of $\dot{P}=-9.03\left(0.11\right)\times {10}^{-12}$s s⁻¹. Few systems have been known to spin up, most notably transitional millisecond pulsars and cataclysmic binaries, both of which seem unlikely progenitors for CHIME J1634+44. If the period was only associated with the spin of the object, then the spin-up is likely generated by accretion of material from a companion. If, however, the radio pulse period and the orbital period are locked, as appears to be the case for two other LPTs, the spin-up of CHIME J1634+44 could be driven by gravitational-wave radiation. 

We present a 34 yr timing solution of the redback pulsar system Terzan 5A (Ter5A). Ter5A, also known as B1744−24A or J1748−2446A, has a 11.56 ms pulse period, a ~0.1 Msun dwarf companion star, and an orbital period of 1.82 hr. Ter5A displays highly variable eclipses and orbital perturbations. Using new timing techniques, we have determined a phase-connected timing solution for this system over 34 yr. This is the longest ever published for a redback pulsar. We find that the pulsar’s spin variability is much larger than most globular cluster pulsars. In fact, of the nine redback pulsars with published or in-preparation long-term timing solutions, Ter5A is by far the noisiest. We see no evidence of strong correlations between orbital and spin variability of the pulsar. We also find that long-term astrometric timing measurements are likely too contaminated by this variability to be usable, and therefore they require careful short-term timing to determine reasonable positions. Finally, we measure an orbital period contraction of  −2.5(3) x 10^-13, which is likely dominated by the general relativistic orbital decay of the system. The effects of the orbital variability due to the redback nature of the pulsar are not needed to explain the observed orbital period derivative, but they are constrained to less than ~30% of the observed value. 

Abstract Alongside the population of several hundred radio millisecond pulsars currently known in Milky Way globular clusters, a subset of six slowly spinning pulsars (spin periods 0.3–4 s) are also observed. With inferred magnetic fields ​​​​​​≳10¹¹G and characteristic ages ≲​​​​​​10⁸yr, explaining the formation of these apparently young pulsars in old stellar populations poses a major challenge. One popular explanation is that these objects are not actually young but instead have been partially spun up via accretion from a binary companion. In this scenario, accretion in a typical low-mass X-ray binary (LMXB) is interrupted by a dynamical encounter with a neighboring object in the cluster. Instead of complete spin-up to millisecond spin periods, the accretion is halted prematurely, leaving behind a “partially recycled” neutron star. In this Letter, we use a combination of analytic arguments motivated by LMXB evolution andN-body simulations to show that this partial recycling mechanism is not viable. Realistic globular clusters are not sufficiently dense to interrupt mass transfer on the short timescales required to achieve such slow spin periods. We argue that collapse of massive white dwarfs and/or neutron star collisions are more promising ways to form slow pulsars in old globular clusters. 

Abstract Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we createN-body models that reproduce the ratio of single-to-binary pulsars in Milky Way–like GCs. We focus especially on NGC 6752, a typical core-collapsed cluster with many observed millisecond pulsars. Previous studies suggested that an increased rate of neutron star binary disruption in the densest clusters could explain the overabundance of single pulsars in these systems. Here, we demonstrate that binary disruption is ineffective and instead we propose that two additional dynamical processes play dominant roles: (1) tidal disruption of main-sequence stars by neutron stars and (2) gravitational collapse of heavy white dwarf binary merger remnants. Neutron stars formed through these processes may also be associated with fast radio bursts similar to those observed recently in an extragalactic GC. 

Abstract The Fourier domain acceleration search (FDAS) is an effective technique for detecting faint binary pulsars in large radio astronomy data sets. This paper quantifies the sensitivity impact of reducing numerical precision in the graphics processing unit (GPU)-accelerated FDAS pipeline of the AstroAccelerate (AA) software package. The prior implementation used IEEE-754 single-precision in the entire binary pulsar detection pipeline, spending a large fraction of the runtime computing GPU-accelerated fast Fourier transforms. AA has been modified to use bfloat16 (and IEEE-754 double-precision to provide a “gold standard” comparison) within the Fourier domain convolution section of the FDAS routine. Approximately 20,000 synthetic pulsar filterbank files representing binary pulsars were generated using SIGPROC with a range of physical parameters. They have been processed using bfloat16, single-precision, and double-precision convolutions. All bfloat16 peaks are within 3% of the predicted signal-to-noise ratio of their corresponding single-precision peaks. Of 14,971 “bright” single-precision fundamental peaks above a power of 44.982 (our experimentally measured highest noise value), 14,602 (97.53%) have a peak in the same acceleration and frequency bin in the bfloat16 output plane, while in the remaining 369 the nearest peak is located in the adjacent acceleration bin. There is no bin drift measured between the single- and double-precision results. The bfloat16 version of FDAS achieves a speedup of approximately 1.6× compared to single-precision. A comparison between AA and the PRESTO software package is presented using observations collected with the GMRT of PSR J1544+4937, a 2.16 ms black widow pulsar in a 2.8 hr compact orbit.