Search

Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the line of sight. We report the observation of annual and orbital variations in the curvature of scintillation arcs over a period of 16 yr for the bright millisecond pulsar, PSR J0437-4715. These variations are the signature of the relative transverse motions of Earth, the pulsar, and the scattering medium, which we model to obtain precise measurements of parameters of the pulsar's binary orbit and the scattering medium itself. We observe two clear scintillation arcs in most of our >5000 observations, and we show that they originate from scattering by thin screens located at distances D 1 = 89.8 ± 0.4 pc and D 2 = 124 ± 3 pc from Earth. The best-fit scattering model we derive for the brightest arc yields the pulsar's orbital inclination angle, i = 137.°1 ± 0.°3, and longitude of ascending node, Ω = 206.°3 ± 0.°4. Using scintillation arcs for precise astrometry and orbital dynamics can be superior to modeling variations in the diffractive scintillation timescale, because the arc curvature is independent of variations in the level of turbulence of interstellar plasma. This technique can be used in combination with pulsar timing to determine the full three-dimensional orbital geometries of binary pulsars and provides parameters essential for testing theories of gravity and constraining neutron star masses.

We report here the results obtained from a systematic optical photometric survey aimed at finding new compact binary millisecond pulsars (also known as “spiders”): the COmpact BInary PULsar SEarch. We acquired multiband optical images over 1 yr around 33 unidentified Fermi Large Area Telescope sources, selected as pulsar candidates based on their curved GeV spectra and steady γ -ray emission. We present the discovery of four optical variables coinciding with the Fermi sources 3FGL J0737.2−3233, 3FGL J2117.6+3725 (two systems in this field), and 3FGL J2221.6+6507, which we propose as new candidate spider systems. Indeed, they all show optical flux modulation consistent with orbital periods of 0.3548(5), 0.25328(6), 0.441961(2), and 0.165(4) days, respectively, with amplitudes ≳0.3 mag and colors compatible with companion star temperatures of 5000–6000 K. These properties are consistent with the “redback” subclass of spider pulsars. If confirmed as a millisecond pulsar, 3FGL J0737.2−3233 will be the closest known spider to Earth ( $D={659}_{-20}^{+16}\ \mathrm{pc}$ , from Gaia-DR3 parallax). We searched and did not find any X-ray sources matching our four candidates, placing 3 σ upper limits of ∼10 ^31 –10 ^32 erg s ^−1 (0.3–10 keV) on their soft X-ray luminosities. We also present and discuss other multiwavelength information on our spider candidates, from infrared to X-rays.

Pulsar timing array experiments have recently uncovered evidence for a nanohertz gravitational wave background by precisely timing an ensemble of millisecond pulsars. The next significant milestones for these experiments include characterizing the detected background with greater precision, identifying its source(s), and detecting continuous gravitational waves from individual supermassive black hole binaries. To achieve these objectives, generating accurate and precise times of arrival of pulses from pulsar observations is crucial. Incorrect polarization calibration of the observed pulsar profiles may introduce errors in the measured times of arrival. Further, previous studies have demonstrated that robust polarization calibration of pulsar profiles can reduce noise in the pulsar timing data and improve timing solutions. In this paper, we investigate and compare the impact of different polarization calibration methods on pulsar timing precision using three distinct calibration techniques: the Ideal Feed Assumption (IFA), Measurement Equation Modeling (MEM), and Measurement Equation Template Matching (METM). Three NANOGrav pulsars—PSRs J1643−1224, J1744−1134, and J1909−3744—observed with the 800 MHz and 1.5 GHz receivers at the Green Bank Telescope (GBT) are utilized for our analysis. Our findings reveal that all three calibration methods enhance timing precision compared to scenarios where no polarization calibration is performed. Additionally, among the three calibration methods, the IFA approach generally provides the best results for timing analysis of pulsars observed with the GBT receiver system. We attribute the comparatively poorer performance of the MEM and METM methods to potential instabilities in the reference noise diode coupled to the receiver and temporal variations in the profile of the reference pulsar, respectively.

Recent constraints on neutron star mass and radius have advanced our understanding of the equation of state (EOS) of cold dense matter. Some of them have been obtained by modeling the pulses of three millisecond X-ray pulsars observed by the Neutron Star Interior Composition Explorer (NICER). Here, we present a Bayesian parameter inference for a fourth pulsar, PSR J1231−1411, using the same technique with NICER and XMM-Newton data. When applying a broad mass-inclination prior from radio timing measurements and the emission region geometry model that can best explain the data, we find likely converged results only when using a limited radius prior. If limiting the radius to be consistent with the previous observational constraints and EOS analyses, we infer the radius to be 12.6 ± 0.3 km and the mass to be ${1.04}_{-0.03}^{+0.05}$ M _⊙ , each reported as the posterior credible interval bounded by the 16% and 84% quantiles. If using an uninformative prior but limited between 10 and 14 km, we find otherwise similar results, but ${R}_{\mathrm{eq}}={13.5}_{-0.5}^{+0.3}$ km for the radius. In both cases, we find a nonantipodal hot region geometry where one emitting spot is at the equator or slightly above, surrounded by a large colder region, and where a noncircular hot region lies close to southern rotational pole. If using a wider radius prior, we only find solutions that fit the data significantly worse. We discuss the challenges in finding the better fitting solutions, possibly related to the weak interpulse feature in the pulse profile.

We present a new method to combine multimass equilibrium dynamical models and pulsar timing data to constrain the mass distribution and remnant populations of Milky Way globular clusters (GCs). We first apply this method to 47 Tuc, a cluster for which there exists an abundance of stellar kinematic data and which is also host to a large population of millisecond pulsars. We demonstrate that the pulsar timing data allow us to place strong constraints on the overall mass distribution and remnant populations even without fitting on stellar kinematics. Our models favor a small population of stellar-mass black holes (BHs) in this cluster (with a total mass of ${446}_{-72}^{+75}{{M}}_{\odot }$ ), arguing against the need for a large (>2000 M _⊙ ) central intermediate-mass BH. We then apply the method to Terzan 5, a heavily obscured bulge cluster that hosts the largest population of millisecond pulsars of any Milky Way GC and for which the collection of conventional stellar kinematic data is very limited. We improve existing constraints on the mass distribution and structural parameters of this cluster and place stringent constraints on its black hole content, finding an upper limit on the mass in BHs of ∼4000 M _⊙ . This method allows us to probe the central dynamics of GCs even in the absence of stellar kinematic data and can be easily applied to other GCs with pulsar timing data, for which data sets will continue to grow with the next generation of radio telescopes.

Pulse profile modeling (PPM), the technique used to infer mass, radius, and geometric parameters for rotation-powered millisecond pulsars using data from the Neutron Star Interior Composition Explorer (NICER), relies on relativistic ray-tracing of thermal X-ray photons from hot spots on the neutron star surface to the observer. To verify our ray-tracing codes we have in the past conducted cross tests for simple hot spot geometries, focusing primarily on the implementation of the spacetime model. In this paper, we present verification for test problems that explore the more complex hot spot geometries that are now being employed in the NICER PPM analyses. We conclude that the accuracy of our computed waveforms is in general sufficiently high for analyses of current NICER data sets. We have however identified some extreme configurations where extra care may be needed.

We present phase-connected timing ephemerides, polarization pulse profiles, and Faraday rotation measurements of 12 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in the Commensal Radio Astronomy FAST Survey. The observational data for each pulsar span at least 1 yr. Among them, PSR J1840+2843 shows subpulse drifting, and five pulsars are detected to exhibit pulse nulling phenomena. PSR J0640−0139 and PSR J2031−1254 are isolated millisecond pulsars (MSPs) with stable spin-down rates ( $\dot{P}$ ) of 4.8981(6) × 10 ^−20 s s ^−1 and 6.01(2) × 10 ^−21 s s ^−1 , respectively. Additionally, one pulsar (PSR J1602−0611) is in a neutron star−white dwarf (WD) binary system with an 18.23-day orbit and a companion of ≤0.65 M _⊙ . PSR J1602−0611 has a spin period, companion mass, and orbital eccentricity that are consistent with the theoretical expectations for MSP−helium WD (He WD) systems. Therefore, we believe that it might be an MSP−He WD binary system. The locations of PSR J1751−0542 and PSR J1840+2843 on the $P-\dot{P}$ diagram are beyond the traditional death line. This indicates that FAST has discovered some low- $\dot{E}$ pulsars, contributing new samples for testing pulsar radiation theories. We estimated the distances of these 12 pulsars based on NE2001 and YMW16 electron density models, and our work enhances the data set for investigating the electron density model of the Galaxy.

PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be 2.08 ± 0.07 M _⊙ . Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. and Riley et al. reported measurements of the radius of PSR J0740+6620 based on Neutron Star Interior Composition Explorer (NICER) observations accumulated through 2020 April 17, and an exploratory analysis utilizing NICER background estimates and a data set accumulated through 2021 December 28 was presented in Salmi et al. Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 2022 April 21, totaling ∼1.1 Ms additional exposure compared to the data set analyzed in Miller et al. and Riley et al., and to data from XMM-Newton observations. We find that the equatorial circumferential radius of PSR J0740+6620 is ${12.92}_{-1.13}^{+2.09}$ km (68% credibility), a fractional uncertainty ∼83% the width of that reported in Miller et al., in line with statistical expectations given the additional data. If we were to require the radius to be less than 16 km, as was done in Salmi et al., then our 68% credible region would become $R={12.76}_{-1.02}^{+1.49}$ km, which is close to the headline result of Salmi et al. Our updated measurements, along with other laboratory and astrophysical constraints, imply a slightly softer equation of state than that inferred from our previous measurements.

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.

The companion to PSR J1622-0315, one of the most compact known redback millisecond pulsars, shows extremely low irradiation despite its short orbital period. We model this system to determine the binary parameters, combining optical observations from the New Technology Telescope in 2017 and the Nordic Optical Telescope in 2022 with the binary modeling code ICARUS . We find a best-fit neutron star mass of 2.3 ± 0.4 M _⊙ , and a companion mass of 0.15 ± 0.02 M _⊙ . We detect for the first time low-level irradiation from asymmetry in the minima as well as a change in the asymmetry of the maxima of its light curves over five years. Using starspot models, we find better fits than those from symmetric direct heating models, with consistent orbital parameters. We discuss an alternative scenario where the changing asymmetry is produced by a variable intrabinary shock. In summary, we find that PSR J1622-0315 combines low irradiation with variable light-curve asymmetry and a relatively high neutron star mass.

This paper demonstrates the impact of state-of-the-art instrumental calibration techniques on the precision of arrival times obtained from 9.6 yr of observations of millisecond pulsars using the Murriyang 64 m CSIRO Parkes Radio Telescope. Our study focuses on 21 cm observations of 25 high-priority pulsars that are regularly observed as part of the Parkes Pulsar Timing Array project, including those predicted to be the most susceptible to calibration errors. We employ measurement equation template matching (METM) for instrumental calibration and matrix template matching (MTM) for arrival time estimation, resulting in significantly improved timing residuals with up to a sixfold reduction in white noise compared to arrival times estimated using scalar template matching and conventional calibration based on the ideal feed assumption. The median relative reduction in white noise is 33%, and the maximum absolute reduction is 4.5 μ s. For PSR J0437−4715, METM and MTM reduce the best-fit power-law amplitude (2.7 σ ) and spectral index (1.7 σ ) of the red noise in the arrival time residuals, which can be tentatively interpreted as mitigation of 1/ f noise due to otherwise unmodeled steps in polarimetric response. These findings demonstrate the potential to directly enhance the sensitivity of pulsar timing array experiments through more accurate methods of instrumental calibration and arrival time estimation.

This study presents an orbital-phase-dependent analysis of three black widow spider millisecond pulsars (BW MSPs), aiming to investigate the magnetic field within the eclipse environment. The ultrawide-bandwidth low-frequency receiver of the Parkes Murriyang radio telescope is utilized for full polarization observations covering frequencies from 704−4032 MHz. Depolarization of pulsed emission is observed during the eclipse phase of three BW MSPs, namely PSR J0024−7204J, PSR J1431−4715, and PSR J1959+2048, consistent with previous studies of other BW MSPs. We estimated orbital-phase-dependent rotation measure values for these MSPs. The wide bandwidth observations also provided constraints on the eclipse cutoff frequency for these BW MSPs. For PSR J0024−7204J, we report temporal variation of the eclipse cutoff frequency coupled with changes in the electron column density within the eclipse medium across six observed eclipses. Moreover, the eclipse cutoff frequency for PSR J1431−4715 is determined to be 1251 ± 80 MHz, leading to the conclusion that synchrotron absorption is the primary mechanism responsible for the eclipsing. Additionally, for PSR J1959+2048, the estimated cutoff frequency exceeded 1400 MHz, consistent with previous studies. With this investigation, we have doubled the sample size of BW MSPs with orbital-phase-resolved studies, allowing a better probe of the eclipse environment.

Millisecond pulsars are subject to accelerations in globular clusters (GCs) that manifest themselves in both the first and second spin period time derivatives, and can be used to explore the mass distribution of the potentials they inhabit. Here we report on over 20 yr of pulsar timing observations of five millisecond radio pulsars in the core of the core-collapse GC NGC 6752 with the Parkes (Murriyang) and MeerKAT radio telescopes, which have allowed us to measure the proper motions, positions, and first and second time derivatives of the pulsars. The pulsar timing parameters indicate that all the pulsars in the core experience accelerations and jerks that can be explained only if an amount of nonluminous mass of at least 2.56 × 10 ^3 M _⊙ is present in the core of NGC 6752. On the other hand, our studies highly disfavor the presence of an intermediate-mass black hole at the center of the cluster, with a mass equal to or greater than ∼3000 M _⊙ .

We present the phase-connected timing solutions of all five pulsars in globular cluster M3 (NGC 5272), namely PSRs M3A to F (PSRs J1342+2822A to F), with the exception of PSR M3C, from FAST archival data. In these timing solutions, those of PSRs M3E and F are obtained for the first time. We find that PSRs M3E and F have low-mass companions and are in circular orbits with periods of 7.1 and 3.0 days, respectively. For PSR M3C, we have not detected its signal in all 41 observations. We found no X-ray counterparts for these pulsars in archival Chandra images in the band of 0.2–20 keV. From the autocorrelation function analysis of M3A and M3B’s dynamic spectra, the scintillation timescale ranges from 7.0 ± 0.3 to 60.0 ± 0.6 minutes, and the scintillation bandwidth ranges from 4.6 ± 0.2 to 57.1 ± 1.1 MHz. The measured scintillation bandwidths from the dynamic spectra indicate strong scintillation, and the scattering medium is anisotropic. From the secondary spectra, we captured a scintillation arc only for PSR M3B with a curvature of 649 ± 23 m ^−1 mHz ^−2 .

Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio frequency. However, pulse arrival times are also affected by radio-frequency-dependent “chromatic” noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect the achromatic noise properties in each pulsar. The models we compare include existing DM models used by the North American Nanohertz Observatory for Gravitational waves (NANOGrav) and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar data set. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from ${\mathrm{log}}_{10}{A}_{\mathrm{RN}}=-{14.1}_{-0.1}^{+0.1}$ to $-{14.7}_{-0.5}^{+0.3}$ , and the spectral index broadens from ${\gamma }_{\mathrm{RN}}={2.6}_{-0.4}^{+0.5}$ to ${\gamma }_{\mathrm{RN}}={3.5}_{-0.9}^{+1.2}$ . We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs.

We use cyclic spectroscopy to perform high-frequency resolution analyses of multihour baseband Arecibo observations of the millisecond pulsar PSR B1937+21. This technique allows for the examination of scintillation features in far greater detail than is otherwise possible under most pulsar timing array observing setups. We measure scintillation bandwidths and timescales in each of eight subbands across a 200 MHz observing band in each observation. Through these measurements we obtain intra-epoch estimates of the frequency scalings for scintillation bandwidth and timescale. Thanks to our high-frequency resolution and the narrow scintles of this pulsar, we resolve scintillation arcs in the secondary spectra due to the increased Nyquist limit, which would not have been resolved at the same observing frequency with a traditional filterbank spectrum using NANOGrav’s current time and frequency resolutions, and the frequency-dependent evolution of scintillation arc features within individual observations. We observe the dimming of prominent arc features at higher frequencies, possibly due to a combination of decreasing flux density and the frequency dependence of the plasma refractive index of the interstellar medium. We also find agreement with arc curvature frequency dependence predicted by Stinebring et al. in some epochs. Thanks to the frequency-resolution improvement provided by cyclic spectroscopy, these results show strong promise for future such analyses with millisecond pulsars, particularly for pulsar timing arrays, where such techniques can allow for detailed studies of the interstellar medium in highly scattered pulsars without sacrificing the timing resolution that is crucial to their gravitational-wave detection efforts.

This paper presents an investigation of the spectral properties of 10 millisecond pulsars (MSPs) discovered by the uGMRT, observed from 2017 to 2023 using band 3 (300–500 MHz) and 4 (550–750 MHz) of the uGMRT. For these MSPs, we have reported a range of spectral indices from ∼0 to −4.8, while averaging the full observing band and all the observing epochs. For every MSP, we calculated the mean flux density across 7–8 subbands, each with approximately 25 MHz bandwidth spanning band 3 and band 4. We computed their modulation indices as well as average and maximum-to-median flux density within each subband. Using the temporal variations of the flux density we calculated the refractive scintillation timescales and estimated a structure function with a time lag for eight MSPs in the sample. We note a significant temporal evolution of the in-band spectra, classified into three categories based on the nature of the best-fit power-law spectra, having single positive spectral indices, multiple broken power laws, and single negative spectral indices. Additionally, indications of a low-frequency turnover and temporal variations of the turnover frequency (to the extent that turnover was observed for some of the epochs while not seen for the rest) were noted for all the MSPs. To the best of our knowledge, this is the first systematic investigation probing temporal changes in MSP spectra as well as in the turnover frequency. Future exploration with dense monitoring combined with the modeling of spectra can provide vital insight into the intrinsic emission properties of MSPs and interstellar medium properties.

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.

Neutron stars represent unique laboratories, offering insights into the physics of supranuclear-density matter and serving as potential hosts for dark matter. This study explores the impact of dark matter cores on rapidly rotating neutron stars through the two-fluid approximation, assuming minimal interaction between baryonic matter and dark matter. The investigation employs phenomenological models for fermionic and bosonic dark matter, revealing that universal relations governing mass and radius changes due to rotation remain largely unaffected in the presence of a dark matter core. Specifically, for a 5% dark matter mass fraction, the percent deviations in total mass ( M _tot ), the baryonic equatorial radius ( R _Be ), and polar-to-equatorial baryonic radius ratio ( R _ratioB ) are within 3.9%, 1.8%, and 1.4%, respectively. These findings suggest that the universal relations governing neutron star shape can be utilized to infer constraints on the properties of dark matter cores even in cases where the dark matter significantly softens the neutron star’s equation of state.

The Fermi All-sky Variability Analysis (FAVA) provides a photometric alternative for identifying week-long gamma-ray flares across the entire sky while being independent of any diffuse Galactic or isotropic emission model. We reviewed 779 weeks of Fermi Large Area Telescope data analyzed by FAVA to estimate the rate and origin of Galactic gamma-ray flares, and to search for new variable Galactic gamma-ray transients. We report an estimated yearly rate of ∼8.5 Galactic gamma-ray flares yr ^–1 , with ∼1 flare yr ^–1 coming from unknown sources. Out of the known gamma-ray sources that are spatially coincident with these detected flares, we report gamma-ray flares for six of them for the first time. All six are classified as pulsars, or a source of unknown nature but which positionally overlaps with known supernova remnants or pulsar wind nebulae (PWNe). This potentially means these sites are tentative candidates to be the second known site of a variable gamma-ray PWN, after the famous Crab Nebula’s PWN. Additionally, we identify nine unassociated flares that are unlikely to have originated from known gamma-ray sources.

At 327 MHz, the observed emission of PSR B1937+21 is greatly affected by scattering in the interstellar medium, on a timescale of order the pulse period. We use the bright impulsive giant pulses emitted by the pulsar to measure the impulse response of the interstellar medium and then recover the intrinsic emission of the pulsar by deconvolution—revealing fine structure on timescales not normally observable. We find that the intrinsic widths of the main pulse and interpulse in the pulse profile are similar to those measured at higher frequencies. We detect 60,270 giant pulses, which typically appear as narrow, ∼100 ns bursts consisting of one to a few nanoshots with widths ≲ 10 ns. However, about 10% of the giant pulses exhibit multiple bursts that seem to be causally related to each other. We also report the first detection of giant micropulses in PSR B1937+21, primarily associated with the regular main pulse emission. These are distinct from giant pulses not only in the phases at which they occur, but also in their larger widths, of order a microsecond, and steeper energy distribution. These measurements place useful observational constraints on emission mechanisms for giant pulses as well as the regular radio emission of millisecond pulsars.

We have searched for radio pulsations toward 49 Fermi Large Area Telescope (LAT) 1FGL Catalog γ -ray sources using the Green Bank Telescope at 350 MHz. We detected 18 millisecond pulsars (MSPs) in blind searches of the data; 10 of these were discoveries unique to our survey. 16 are binaries, with eight having short orbital periods P _B < 1 day. No radio pulsations from young pulsars were detected, although three targets are coincident with apparently radio-quiet γ -ray pulsars discovered in LAT data. Here, we give an overview of the survey and present radio and γ -ray timing results for the 10 MSPs discovered. These include the only isolated MSP discovered in our survey and six short- P _B binary MSPs. Of these, three have very-low-mass companions ( M _c ≪ 0.1 M _⊙ ) and hence belong to the class of black widow pulsars. Two have more massive, nondegenerate companions with extensive radio eclipses and orbitally modulated X-ray emission consistent with the redback class. Significant γ -ray pulsations have been detected from nine of the discoveries. This survey and similar efforts suggest that the majority of Galactic γ -ray sources at high Galactic latitudes are either MSPs or relatively nearby nonrecycled pulsars, with the latter having on average a much smaller radio/ γ -ray beaming ratio as compared to MSPs. It also confirms that past surveys suffered from an observational bias against finding short- P _B MSP systems.

Noise characterization for pulsar-timing applications accounts for interstellar dispersion by assuming a known frequency dependence of the delay it introduces in the times of arrival (TOAs). However, calculations of this delay suffer from misestimations due to other chromatic effects in the observations. The precision in modeling dispersion is dependent on the observed bandwidth. In this work, we calculate the offsets in infinite-frequency TOAs due to misestimations in the modeling of dispersion when using varying bandwidths at the Green Bank Telescope. We use a set of broadband observations of PSR J1643−1224, a pulsar with unusual chromatic timing behavior. We artificially restricted these observations to a narrowband frequency range, then used both the broad- and narrowband data sets to calculate residuals with a timing model that does not account for time variations in the dispersion. By fitting the resulting residuals to a dispersion model and comparing the fits, we quantify the error introduced in the timing parameters due to using a reduced frequency range. Moreover, by calculating the autocovariance function of the parameters, we obtained a characteristic timescale over which the dispersion misestimates are correlated. For PSR J1643−1224, which has one of the highest dispersion measures (DM) in the NANOGrav pulsar timing array, we find that the infinite-frequency TOAs suffer from a systematic offset of ∼22 μ s due to incomplete frequency sampling, with correlations over about one month. For lower-DM pulsars, the offset is ∼7 μ s. This error quantification can be used to provide more robust noise modeling in the NANOGrav data, thereby increasing the sensitivity and improving the parameter estimation in gravitational wave searches.

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.

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.

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.

The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between 2021 April 16 and 17 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multifrequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment and the 100 m Green Bank Telescope in a 3 yr period encompassing the shape change event, between 2020 February and 2023 February. As of 2023 February, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying time-of-arrival residuals display a strong nonmonotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency, ν ) nor a change in dispersion measure alone (which would produce a delay proportional to ν ^−2 ). However, it does bear some resemblance to the two previous “chromatic timing events” observed in J1713+0747, as well as to a similar event observed in PSR J1643−1224 in 2015.

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.

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 statistical F -test and the quadratic dependence of the reduced χ ^2 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.

We investigated the multiband emission from the pulsar binaries XSS J12270−4859, PSR J2039−5617, and PSR J2339−0533, which exhibit orbital modulation in the X-ray and gamma-ray bands. We constructed the sources’ broadband spectral energy distributions and multiband orbital light curves by supplementing our X-ray measurements with published gamma-ray results, and we modeled the data using intrabinary shock (IBS) scenarios. While the X-ray data were well explained by synchrotron emission from electrons/positrons in the IBS, the gamma-ray data were difficult to explain with the IBS components alone. Therefore, we explored other scenarios that had been suggested for gamma-ray emission from pulsar binaries: (1) inverse-Compton emission in the upstream unshocked wind zone and (2) synchrotron radiation from electrons/positrons interacting with the kilogauss magnetic field of the companion. Scenario (1) requires that the bulk motion of the wind substantially decelerates to ∼1000 km s ^−1 before reaching the IBS for increased residence time, in which case the formation of a strong shock is untenable, inconsistent with the X-ray phenomenology. Scenario (2) can explain the data if we assume the presence of electrons/positrons with a Lorentz factor ∼ 10 ^8 (∼0.1 PeV) that pass through the IBS and tap a substantial portion of the pulsar voltage drop. These findings raise the possibility that the orbitally modulating gamma-ray signals from pulsar binaries can provide insights into the flow structure and energy conversion within pulsar winds and particle acceleration nearing PeV energies in pulsars. These signals may also yield greater understanding of kilogauss magnetic fields potentially hosted by the low-mass stars in these systems.

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.

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) × 10 ^9 and (1.81 ± 0.01) × 10 ^9 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.

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.

We present results from long-term monitoring of frequency-dependent eclipses of the radio emission from PSR J1544+4937, which is a black widow spider millisecond pulsar (MSP) in a compact binary system. The majority of such systems often exhibit relatively long-duration radio eclipses caused by ablated material from their companion stars. With the wide spectral bandwidth of the upgraded Giant Metrewave Radio Telescope, we present the first systematic study of temporal variation of eclipse cutoff frequency. With decade-long monitoring of 39 eclipses for PSR J1544+4937, we notice significant changes in the observed cutoff frequency ranging from 343 ± 7 to ≥740 MHz. We also monitored changes in eclipse cutoff frequency on timescales of tens of days and observed a maximum change of ≥315 MHz between observations that were separated by 22 days. In addition, we observed a change of ∼47 MHz in eclipse cutoff frequency between adjacent orbits, i.e., on timescales of ∼2.9 hr. We infer that such changes in the eclipse cutoff frequency depict an eclipse environment for the PSR J1544+4937 system that is dynamically evolving, where, along with the change in electron density, the magnetic field could also be varying. We also report a significant correlation between the eclipse cutoff frequency and the mass-loss rate of the companion. This study provides the first direct evidence of mass-loss rate affecting the frequency-dependent eclipsing in a spider MSP.

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 create N -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.

The discovery and timing follow up of millisecond pulsars (MSPs) are necessary not just for their usefulness in pulsar timing arrays (PTAs) but also for investigating their own intriguing properties. In this work, we provide the findings of the decade-long timing of four MSPs discovered by the Giant Meterwave Radio Telescope (GMRT), including their timing precision, model parameters, and newly detected proper motions. We compare the timing results for these MSPs before and after the GMRT upgrade in 2017 and characterize the improvement in timing precision due to the bandwidth upgrade. We discuss the suitability of these four GMRT MSPs as well as the usefulness of the decade-long timing data for PTA experiments. These data may aid in the global effort to improve the signal-to-noise ratios of recently detected signature of gravitational waves in cross-correlation statistics of residuals of MSPs.

We present the results of an ultradeep radio continuum survey, containing ∼480 hr of observations, of the Galactic globular cluster 47 Tucanae with the Australia Telescope Compact Array. This comprehensive coverage of the cluster allows us to reach rms noise levels of 1.19 μ Jy beam ^−1 at 5.5 GHz, 940 nJy beam ^−1 at 9 GHz, and 790 nJy beam ^−1 in a stacked 7.25 GHz image. This is the deepest radio image of a globular cluster and the deepest image ever made with the Australia Telescope Compact Array. We identify ATCA J002405.702-720452.361, a faint (6.3 ± 1.2 μ Jy at 5.5 GHz, 5.4 ± 0.9 μ Jy at 9 GHz), flat-spectrum ( α = −0.31 ± 0.54) radio source that is positionally coincident with the cluster center and potentially associated with a faint X-ray source. No convincing optical counterpart was identified. We use radio, X-ray, optical, and UV data to show that explanations involving a background active galactic nucleus, a chromospherically active binary, or a binary involving a white dwarf are unlikely. The most plausible explanations are that the source is an undiscovered millisecond pulsar or a weakly accreting black hole. If the X-ray source is associated with the radio source, the fundamental plane of black-hole activity suggests a black hole mass of ∼54–6000 M _⊙ , indicating an intermediate-mass black hole or a heavy stellar-mass black hole.

We report detection results of nine millisecond pulsars (MSPs) at 8600 MHz using simultaneous 2250 and 8600 MHz observations conducted with the Shanghai Tian Ma Radio Telescope. Mainly benefiting from updated ephemerids with 2250 MHz observations, integrated profiles of all nine MSPs at 8600 MHz are successfully obtained by coherently adding multi-epoch (3–83 epochs) observation data spanning from 19–1210 days, which significantly increases the number of MSPs with published profiles (from 4 to 11) above 8000 MHz, as seven of our target MPSs had no related results previously. Combining our new flux density and pulse width measurements with previous low-frequency results, we study their integrated profile evolution and spectral behaviors in a wider frequency range. We find their component separations and pulse widths remain almost constant, which is consistent with previous findings. While dramatic evolution of integrated profiles exists at low frequencies, we observe a potential end of the related evolution around 5 GHz in eight MSPs. The spectra of four MSPs are found to deviate from a single power law, and we fit them with a broken power law. The change in the profile of PSR J1713+0747, which started around MJD 59320−59321, seems to be more prominent as the observation frequency increases. Compared with the effects of the interstellar medium, we prefer to explain this event as some changes in the magnetosphere. We also find its integrated profile possibly had not recovered to the pre-event state until MJD 59842–59857.

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems colocated with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to ≤11 known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with six undetected in radio. Overall, ≥236 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot{E}$ decreases to its observed minimum near 10 ^33 erg s ^−1 , approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties, and fit results to guide and be compared with modeling results.

We studied the broadband X-ray timing and spectral behaviors of the newly confirmed accreting millisecond X-ray pulsar MAXI J1816−195 during its 2022 outburst. We used data from the Insight-HXMT Medium Energy (ME) and High Energy (HE) telescopes, NICER, and NuSTAR that cover the energy range between 0.8 and 210 keV. A coherent timing analysis of solely the Insight-HXMT HE data across the full outburst revealed a complex behavior of the timing residuals, also prominently visible in the independent Insight-HXMT ME and NICER data, particularly at the rising part of the outburst and at the very end in the NICER data. Therefore, we broke down the full outburst into a (noisy) rising part, covering only about five days, from MJD 59737.0 to 59741.9, and a decaying part, lasting for 19 days across MJD 59741.9–59760.6. Fitting for the decaying part, a timing model including a frequency ν and frequency time derivative $\dot{\nu }$ component yielded a value of (+9.0 ± 2.1) × 10 ^−14 Hz s ^−1 for $\dot{\nu }$ , which could be interpreted as a spinup under our model assumptions. We detected X-ray pulsations up to ∼95 keV in a combination of Insight-HXMT HE observations. The pulse profiles were quite stable over the whole outburst and could be well described by a truncated Fourier series using two harmonics, the fundamental and the first overtone. Both components kept alignment in the range 0.8–64 keV. The joint and time-averaged NICER and Insight-HXMT spectra in the energy range 1–150 keV were well fitted by the absorbed Comptonization model compps plus disk blackbody with two additional Gaussian components. Using the bolometric flux and spinup values both evaluated during the decay phase, we determined a magnetic field strength of (0.2–2) × 10 ^8 G for MAXI J1816−195.

We present the discovery and timing solutions of four millisecond pulsars (MSPs) discovered in the Arecibo 327 MHz Drift-Scan Pulsar Survey. Three of these pulsars are in binary systems, consisting of a redback (PSR J2055+1545), a black widow (PSR J1630+3550), and a neutron star–white dwarf binary (PSR J2116+1345). The fourth MSP, PSR J2212+2450, is isolated. We present the multiyear timing solutions as well as polarization properties across a range of radio frequencies for each pulsar. We perform a multiwavelength search for emission from these systems and find an optical counterpart for PSR J2055+1545 in Gaia DR3, as well as a gamma-ray counterpart for PSR J2116+1345 with the Fermi-LAT telescope. Despite the close colocation of PSR J2055+1545 with a Fermi source, we are unable to detect gamma-ray pulsations, likely due to the large orbital variability of the system. This work presents the first two binaries found by this survey with orbital periods shorter than a day; we expect to find more in the 40% of the survey data that have yet to be searched.

M71E is a spider pulsar (i.e., a millisecond pulsar with a tight binary companion) with the shortest known orbital period of P = 53.3 minutes, which was discovered by Pan et al. Their favored evolutionary model suggests that it bridges between two types of spider pulsars, namely, it descended from a “redback” and will become a “black widow.” Using Hubble Space Telescope (HST) archival imaging data, we report the first optical identification of its companion COM-M71E. The HST and pulsar timing coordinates are in excellent agreement (within ∼10 mas). If M71E is associated with the globular cluster M71, our measured brightness of COM-M71E ( m _F606W ≈ 25.3) is broadly consistent with the expectation from Pan et al.'s preferred binary evolutionary model of a stripped dwarf companion, while it is also compatible with an ultralow-mass degenerate companion. Future multiwavelength photometric and spectroscopic observations can characterize the companion and test the evolutionary scenarios.

We present timing solutions for 21 pulsars discovered in 350 MHz surveys using the Green Bank Telescope (GBT). All were discovered in the Green Bank North Celestial Cap pulsar survey, with the exception of PSR J0957−0619, which was found in the GBT 350 MHz Drift-scan pulsar survey. The majority of our timing observations were made with the GBT at 820 MHz. With a spin period of 37 ms and a 528 days orbit, PSR J0032+6946 joins a small group of five other mildly recycled wide binary pulsars, for which the duration of recycling through accretion is limited by the length of the companion’s giant phase. PSRs J0141+6303 and J1327+3423 are new disrupted recycled pulsars. We incorporate Arecibo observations from the NANOGrav pulsar timing array into our analysis of the latter. We also observed PSR J1327+3423 with the Long Wavelength Array, and our data suggest a frequency-dependent dispersion measure. PSR J0957−0619 was discovered as a rotating radio transient, but is a nulling pulsar at 820 MHz. PSR J1239+3239 is a new millisecond pulsar (MSP) in a 4 days orbit with a low-mass companion. Four of our pulsars already have published timing solutions, which we update in this work: the recycled wide binary PSR J0214+5222, the noneclipsing black widow PSR J0636+5128, the disrupted recycled pulsar J1434+7257, and the eclipsing binary MSP J1816+4510, which is in an 8.7 hr orbit with a redback-mass companion.

The recent observation of a common red-noise process in pulsar timing arrays (PTAs) suggests that the detection of nanohertz gravitational waves might be around the corner. However, in order to confidently attribute this red process to gravitational waves, one must observe the Hellings–Downs curve—the telltale angular correlation function associated with a gravitational-wave background. This effort is complicated by the complex modeling of pulsar noise. Without proper care, misspecified noise models can lead to false-positive detections. Background estimation using “quasi-resampling” methods such as sky scrambles and phase shifts, which use the data to characterize the noise, are therefore important tools for assessing significance. We investigate the ability of current PTA experiments to estimate their background with “quasi-independent” scrambles—characterized by a statistical “match” below the fiducial value: ∣ M ∣ < 0.1. We show that sky scrambling is affected by “saturation” after ${ \mathcal O }(10)$ quasi-independent realizations; subsequent scrambles are no longer quasi-independent. We show that phase scrambling saturates after ${ \mathcal O }(100)$ quasi-independent realizations. With so few independent scrambles, it is difficult to make reliable statements about the ≳5 σ tail of the null distribution of the detection statistic. We sketch out various methods by which one may increase the number of independent scrambles. We also consider an alternative approach, wherein one reframes the background estimation problem so that the significance is calculated using statistically dependent scrambles. The resulting p -value is in principle well defined, but may be susceptible to failure if assumptions about the data are incorrect.

Black widows are millisecond pulsars ablating their companions. The material blown from the companion blocks the radio emission, resulting in radio eclipses. The properties of the eclipse medium are poorly understood. Here, we present direct evidence of the existence of magnetic fields in the eclipse medium of the black widow PSR J2051−0827 using observations made with the Five-hundred-meter Aperture Spherical radio Telescope. We detect a regular decrease in rotation measure (RM) in the egress of the eclipse, changing from 60 to −28.7 rad m ^−2 . The RM gradually changes back to normal when the line of sight moves away from the eclipse. The estimated line-of-sight magnetic field strength in the eclipse medium is ∼0.1 G. The RM reversal could be caused by a change in the magnetic field strength along the line of sight due to binary orbital motion. The RM reversal phenomenon has also been observed in some repeating fast radio bursts (FRBs), and the study of spider pulsars may provide additional information about the origin of FRBs.

We report optical observations of the millisecond pulsar binary system PSR J1622–0315 with the Lulin 1 m telescope in Taiwan and the Lijiang 2.4 m telescope in China between 2019 and 2021. The companion of the pulsar, which is of V ∼ 19 mag, showed ellipsoidal-distorted orbital variations in its light curves. The best-fit model to the light curves, with the binary code PHOEBE , gives a companion mass of 0.122 ± 0.006 M _⊙ . This places PSR J1622–0315 in the spider-system subclass. We compared the properties of PSR J1622–0315 with other spider pulsar binaries for the scalings between the spin-down luminosity derived for the pulsar, irradiation luminosity of the companion, and X-ray luminosity of the binary. We find that pulsar irradiation in PSR J1622–0315 is insignificant and the irradiation luminosity of the transitional millisecond pulsars PSR J1023+0038 and PSR J1227–4853 are the highest among the redback systems.

Pulsar distances are notoriously difficult to measure, and play an important role in many fundamental physics experiments, such as pulsar timing arrays. Here, we perform a cross-match between International PTA pulsars (IPTA) and Gaia's Data Release 2 (DR2) and Data Release 3 (DR3). We then combine the IPTA pulsar’s parallax with its binary companion’s parallax, found in Gaia, to improve the distance measurement to the binary. We find seven cross-matched IPTA pulsars in Gaia DR2, and when using Gaia DR3 we find six IPTA pulsar cross-matches but with seven Gaia objects. Moving from Gaia DR2 to Gaia DR3, we find that the Gaia parallaxes for the successfully cross-matched pulsars improved by 53%, and pulsar distances improved by 29%. Finally, we find that binary companions with a

We report the discovery of a 2.11 ms binary millisecond pulsar during a targeted search of the redback optical candidate coincident with the γ -ray source 3FGL J0212.5+5320 using the Robert C. Byrd Green Bank Telescope (GBT) with the Breakthrough Listen backend at L band. Over a seven month period, five pointings were made near inferior conjunction of the pulsar in its 20.9 hr orbit, resulting in two detections, lasting 12 and 42 minutes. The pulsar dispersion measure (DM) of 25.7 pc cm ^−3 corresponds to a distance of 1.15 kpc in the NE2001 Galactic electron density model, consistent with the Gaia parallax distance of 1.16 ± 0.03 kpc for the companion star. We suspect the pulsar experiences wide orbital eclipses, similar to other redbacks, as well as scintillation and DM delays caused by its interaction with its companion and surroundings. Although the pulsar was only detected over ≈3.7% of the orbit, its measured acceleration is consistent with published binary parameters from optical radial velocity spectroscopy and light-curve modeling of the companion star, and it provides a more precise mass ratio and a projected semimajor axis for the pulsar orbit. We also obtained a refined optical photometric orbit ephemeris, and observed variability of the tidally distorted companion over 7 yr. A hard X-ray light curve from NuSTAR shows expected orbit-modulated emission from the intrabinary shock. The pulsar parameters and photometric ephemeris greatly restrict the parameter space required to search for a coherent timing solution including pulsar spin-down rate, either using Fermi γ -rays or further radio pulse detections.

It has been shown that the nonminimal coupling between geometry and matter can provide models for massive compact stars that are consistent with the conformal bound on the sound speed, $0\leqslant {c}_{s}^{2}\leqslant {c}^{2}/3$ , where the core density approaches a few times the nuclear saturation density. We impose the conformal upper bound on the sound speed on Rastall’s field equations of gravity, with Krori–Barua potentials in the presence of an anisotropic fluid as a matter source, to estimate the radius of the most massive pulsar ever observed, PSR J0952–0607. For its measured mass M = 2.35 ± 0.17 M _⊙ , we obtain a radius R = 14.087 ± 1.0186 km as inferred by the model. We investigate a possible connection between Rastall gravity and the MIT bag model with an equation of state, ${p}_{r}(\rho )\approx {c}_{s}^{2}\left(\rho -{\rho }_{{\rm{s}}}\right)$ , in the radial direction, with ${c}_{s}=c/\sqrt{3}$ and a surface density ρ _s slightly above the nuclear saturation density ρ _nuc = 2.7 × 10 ^14 g cm ^–3 . The corresponding mass–radius diagram is in agreement with our estimated value of the radius and with astrophysical observations of other pulsars at 68% confidence level.

We report on the discovery of the companion star to the millisecond pulsar PSR J1835−3259B in the Galactic globular cluster NGC 6652. Taking advantage of deep photometric archival observations acquired through the Hubble Space Telescope in near-UV and optical bands, we identified a bright and blue object at a position compatible with that of the radio pulsar. The companion is located along the helium-core white dwarf cooling sequence, and the comparison with binary evolution models provides a mass of 0.17 ± 0.02 M _⊙ , a surface temperature of 11,500 ± 1900 K, and a very young cooling age of only 200 ± 100 Myr. The mass and the age of the companion are compatible with a progenitor star of about 0.87 M _⊙ , which started transferring mass to the primary during its evolution along the subgiant branch and stopped during the early red giant branch phase. Combining together the pulsar mass function and the companion mass, we found that this system is observed at an almost edge-on orbit and hosts a neutron star with a mass of 1.44 ± 0.06 M _⊙ , thus suggesting a highly nonconservative mass accretion phase. The young age of the WD companion is consistent with the scenario of a powerful, relatively young MSP indicated by the earlier detection of gamma-rays from this system.