Search for: All records

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 Pulsar timing is a powerful tool that, by accounting for every rotation of a pulsar, precisely measures the spin frequency, spin frequency derivatives, astrometric position, binary parameters when applicable, properties of the interstellar medium, and potentially general relativistic effects. Typically, this process demands fairly stringent scheduling requirements for monitoring observations as well as deep domain knowledge to “phase connect” the timing data. We present an algorithm that automates the pulsar-timing process for binary pulsars, whose timing solutions have an additional level of complexity, although the algorithm works for isolated pulsars as well. Using the statisticalF-test and the quadratic dependence of the reducedχ²near a minimum, the global rotation count of a pulsar can be determined efficiently and systematically. We have used our algorithm to establish timing solutions for two newly discovered binary pulsars, PSRs J1748−2446aq and J1748−2446at, in the globular cluster Terzan 5, using ∼70 Green Bank Telescope observations from the last 13 yr. 

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 Precise localizations of a small number of repeating fast radio bursts (FRBs) using very long baseline interferometry (VLBI) have enabled multiwavelength follow-up observations revealing diverse local environments. However, the 2%–3% of FRB sources that are observed to repeat may not be representative of the full population. Here we use the VLBI capabilities of the full CHIME Outrigger array for the first time to localize a nearby (40 Mpc), bright (kJy), and apparently one-off FRB source, FRB 20250316A, to its environment on 13 pc scales. We use optical and radio observations to place deep constraints on associated transient emission and the properties of its local environment. We place a 5σupper limit ofL_{9.9 GHz} < 2.1 × 10²⁵erg s⁻¹Hz⁻¹on spatially coincident radio emission, a factor of 100 lower than any known compact persistent radio source associated with an FRB. Our Keck Cosmic Webb Imager observations allow us to characterize the gas density, metallicity, nature of gas ionization, dust extinction, and star formation rate through emission line fluxes. We leverage the exceptional brightness and proximity of this source to place deep constraints on the repetition of FRB 20250316A and find that it is inconsistent with all well-studied repeaters given the nondetection of bursts at lower spectral energies. We explore the implications of a measured offset of 190 ± 20 pc from the center of the nearest star formation region in the context of progenitor channels. FRB 20250316A marks the beginning of an era of routine localizations for one-off FRBs on tens of milliarcseconds scales, enabling large-scale studies of their local environments. 

Abstract Pulsar timing is a process of iteratively fitting pulse arrival times to constrain the spindown, astrometric, and possibly binary parameters of a pulsar, by enforcing integer numbers of pulsar rotations between the arrival times. Phase connection is the process of unambiguously determining those rotation numbers between the times of arrival while determining a pulsar timing solution. Pulsar timing currently requires a manual process of step-by-step phase connection performed by individuals. In an effort to quantify and streamline this process, we created the Algorithmic Pulsar Timer (APT), an algorithm that can accurately phase connect and time isolate pulsars. Using the statistical F-test and knowledge of parameter uncertainties and covariances, the algorithm decides what new data to include in a fit, when to add additional timing parameters, and which model to attempt in subsequent iterations. Using these tools, the algorithm can phase-connect timing data that previously required substantial manual effort. We tested the algorithm on 100 simulated systems, with a 99% success rate. APT combines statistical tests and techniques with a logical decision-making process, very similar to the manual one used by pulsar astronomers for decades, and some computational brute force, to automate the often tricky process of isolated pulsar phase connection, setting the foundation for automated fitting of binary pulsar systems. 

Abstract Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 ± 0.8 mm s^–1yr^–1, which deviates from the acceleration predicted by equilibrium Milky Way (MW) models at a 4σlevel. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the MW interacting with orbiting dwarf galaxies (∼1 mm s^–1yr^–1), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in Gaia DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 60 kyr, a companion mass of 0.3M_⊙(indicative of a white dwarf or low-mass star), and a semimajor axis of 1900 au. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios. 

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. 

Abstract We present the results of a deep study of the neutron star (NS) population in the globular cluster M28 (NGC 6626), using the full 330 ks 2002–2015 ACIS data set from the Chandra X-ray Observatory and coordinated radio observations taken with the Green Bank Telescope (GBT) in 2015. We investigate the X-ray luminosity ( L X ), spectrum, and orbital modulation of the seven known compact binary millisecond pulsars in the cluster. We report two simultaneous detections of the redback PSR J1824−2452I (M28I) and its X-ray counterpart at L X = [8.3 ± 0.9] × 10 31 erg s −1 . We discover a double-peaked X-ray orbital flux modulation in M28I during its pulsar state, centered around pulsar inferior conjunction. We analyze the spectrum of the quiescent NS low-mass X-ray binary to constrain its mass and radius. Using both hydrogen and helium NS atmosphere models, we find an NS radius of R = 9.2–11.5 km and R = 13.0–17.5 km, respectively, for an NS mass of 1.4 M ⊙ (68% confidence ranges). We also search for long-term variability in the 46 brightest X-ray sources and report the discovery of six new variable low-luminosity X-ray sources in M28.