We conduct the first-ever pulsar polarization array (PPA) analysis to detect the ultralight axionlike dark matter (ALDM) using the polarization data of 22 millisecond pulsars from the third data release of Parkes pulsar timing array. As one of the major dark matter candidates, the ultralight ALDM exhibits a pronounced wave nature on astronomical scales and offers a promising solution to small-scale structure issues within local galaxies. While the linearly polarized pulsar light travels through the ALDM galactic halo, its position angle (PA) can be subject to an oscillation induced by the ALDM Chern-Simons coupling with electromagnetic field. The PPA is thus especially suited for detecting the ultralight ALDM by correlating polarization data across the arrayed pulsars. To accomplish this task, we develop an advanced Bayesian analysis framework that allows us to construct pulsar PA residual time series, model noise contributions properly and search for pulsar cross-correlations. We find that for an ALDM density of ρ_{0}=0.4  GeV/cm^{3}, the Parkes PPA offers the best global limits on the ALDM Chern-Simons coupling, namely ≲10^{-13.5}-10^{-12.2}  GeV^{-1}, for the mass range of 10^{-22}-10^{-21}  eV. The crucial role of pulsar cross-correlation in recognizing the nature of the derived limits is also highlighted.

The birth-mass function of neutron stars encodes rich information about supernova explosions, double-star evolution and the properties of matter under extreme conditions. To date, it has remained poorly constrained by observations, however. Applying probabilistic corrections to account for mass accreted by recycled pulsars in binary systems to mass measurements of 90 neutron stars, we find that the birth masses of neutron stars can be described by a unimodal distribution that smoothly turns on at 1.1 M-circle dot and peaks at similar to 1.27 M-circle dot, before declining as a steep power law. Such a 'turn-on' power-law distribution is strongly favoured against the widely adopted empirical double-Gaussian model at the 3 sigma level. The power-law shape may be inherited from the initial mass function of massive stars, but the relative dearth of massive neutron stars implies that single stars with initial masses greater than similar to 18 M-circle dot do not form neutron stars, in agreement with the absence of massive red supergiant progenitors of supernovae.

Rotating radio transients (RRATs) are neutron stars that emit sporadic radio bursts. We detected 1955 single pulses from RRAT J1913+1330 using the 19 beam receiver of the Five-hundred-meter Aperture Spherical radio Telescope. These pulses were detected in 19 distinct clusters, with 49.4% of them occurring with a waiting time of one rotation period. The energy distribution of these individual pulses exhibited a wide range, spanning 3 orders of magnitude, reminiscent of repeating fast radio bursts (FRBs). Furthermore, we observed abrupt variations in pulse profile, width, peak flux, and fluence between adjacent sequential pulses. These findings suggest that this RRAT could be interpreted as a pulsar with extreme pulse-to-pulse modulation. The presence of sequential pulse trains during active phases, along with significant pulse variations in profile, fluence, flux, and width, should be intrinsic to a subset of RRATs. Our results indicate that J1913+1330 represents a peculiar source that shares certain properties with populations of nulling pulsars, giant pulses, and FRBs from different perspectives. The dramatic pulse-to-pulse variation observed in J1913+1330 could be attributed to unstable pair creation above the polar cap region and the variation of the site where streaming pairs emit coherently. Exploring a larger sample of RRATs exhibiting similar properties to J1913+1330 has the potential to significantly advance our understanding of pulsars, RRATs, and FRBs.

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 ⊙ .

>With their huge number of stars, many of which are members of binary systems, globular clusters provide a rich source of pulsar progenitors. Those in binary systems will evolve,often spinning up or “recycling” the pulsar through accretion of mass and angular momentum from the evolving companion star. If this accretion process goes through completion,

The birth mass function of neutron stars encodes rich information about supernova explosions, double star evolution, and properties of matter under extreme conditions. To date, it has remained poorly constrained by observations, however. Applying probabilistic corrections to account for mass accreted by recycled pulsars in binary systems to mass measurements of 90 neutron stars, we find that the birth masses of neutron stars can be described by a unimodal distribution that smoothly turns on at 1.1 M_⊙, peaks at 1.27 M_⊙, before declining as a steep power law. Such a “turn-on" power-law distribution is strongly favoured against the widely-adopted empirical double-Gaussian model at the 3 σ level. The power-law shape may be inherited from the initial mass function of massive stars, but the relative dearth of massive neutron stars implies that single stars with initial masses greater than ≈ 18 M_⊙ do not form neutron stars, in agreement with the absence of massive red supergiant progenitors to supernovae.

The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. The timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR J0437–4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Composition Explorer (NICER). The main aim of the NICER mission is to constrain the neutron star equation of state by inferring the compactness ( M p / R ) of the star. Direct measurements of the mass M p from pulsar timing therefore substantially improve constraints on the radius R and the equation of state. Here we use observations spanning 26 yr from Murriyang, the 64 m Parkes radio telescope, to improve the timing model for this pulsar. Among the new precise measurements are the pulsar mass M p = 1.418 ± 0.044 M ⊙ , distance D = 156.96 ± 0.11 pc, and orbital inclination angle i = 137.°506 ± 0.°016, which can be used to inform the X-ray pulse profile models inferred from NICER observations. We demonstrate that these results are consistent between multiple data sets from the Parkes Pulsar Timing Array (PPTA), each modeled with different noise assumptions. Using the longest available PPTA data set, we measure an apparent second derivative of the pulsar spin frequency and discuss how this can be explained either by kinematic effects due to the proper motion and radial velocity of the pulsar or excess low-frequency noise such as a gravitational-wave background.

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1\sigma$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA.

The Parkes 20 cm multibeam pulsar surveys have discovered nearly half of the known pulsars and revealed many distant pulsars with high dispersion measures. Using a sample of 1301 pulsars from these surveys, we have explored the spatial distribution and birth rate of normal pulsars. The pulsar distances used to calculate the pulsar surface density are estimated from the YMW16 electron-density model. When estimating the impact of the Galactic background radiation on our survey, we projected pulsars in the galaxy onto the Galactic plane, assuming that the flux density distribution of pulsars is uniform in all directions, and utilized the most up-to-date background temperature map. We also used an up-to-date version of the ATNF Pulsar Catalogue to model the distribution of pulsar flux densities at 1400 MHz. We derive an improved radial distribution for the pulsar surface density projected onto the Galactic plane, which has a maximum value at ∼4 kpc from the Galactic center. We also derive the local surface density and birth rate of pulsars, obtaining 47 ± 5 kpc ^−2 and ∼4.7 ± 0.5 kpc ^−2 Myr ^−1 , respectively. For the total number of potentially detectable pulsars in the galaxy, we obtain (1.1 ± 0.2) × 10 ^4 and (1.1 ± 0.2)×10 ^5 before and after applying the Tauris & Manchester beaming correction model. The radial distribution function is used to estimate the proportion of pulsars in each spiral arm and the Galactic center.

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 (FAST). We detect a regular decrease in rotation measure (RM) in the egress of eclipse, changing from $60\,\rm rad\,m^{-2}$ to $-28.7\,\rm 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 $\sim 0.1$ G. The RM reversal could be caused by a change of 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 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 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.

The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit.

We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes 'Murriyang' radio telescope. The data span is up to 18 yr with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim$3 yr of more recent data primarily obtained using an ultra-wide-bandwidth receiver system that operates between 704 and 4032 MHz. We provide calibrated pulse profiles, flux density dynamic spectra, pulse times of arrival, and initial pulsar timing models. We describe methods for processing such wide-bandwidth observations and compare this data release with our previous release.

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 co-located 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 $\leq 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 6 undetected in radio. Overall, >235 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.

The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR J0614−3329 and PSR J1902−5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of γ = 13/3, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission.

PSRJ1910-5959A (J1910A) is a binary millisecond pulsar in a 0.837 day circular orbit around a helium white dwarf (HeWD) companion. This pulsar is located 6.3 arcmin away from the centre of the globular cluster NGC6752. Given the large offset, the association of the pulsar to NGC6752 has been debated. We have made use of two decades of archival Parkes 64-m "Murriyang" telescope data and recently carried out observations with the MeerKAT telescope. We obtained Pulse times of arrival using standard data reduction techniques and analysed using Bayesian pulsar timing techniques. We analysed the pulsar's total intensity and polarisation profile, to study the interstellar scattering along the line of sight, and the pulsar's geometry by applying the rotating vector model. We obtain precise measurements of several post-Keplerian parameters: the range $r=0.202(6)T_\odot$ and shape s=0.999823(4) of the Shapiro delay, from which we infer the orbital inclination to be $88.9^{+0.15}_{-0.14}\deg$ and the masses of both the pulsar and the companion to be $1.55(7)M_{\odot}$ and $0.202(6)M_{\odot}$ respectively; a secular change in the orbital period $\dot{P}_{\rm b}=-53^{+7.4}_{-6.0}\times 10^{-15}$\,s\,s$^{-1}$ that proves the association to NGGC6752 and a secular change in the projected semi-major axis of the pulsar $\dot{x}= -40.7^{+7.3}_{-8.2}\times10^{-16}$\,s\,s$^{-1}$ that is likely caused by the spin-orbit interaction from a misaligned HeWD spin, at odds with the likely isolated binary evolution of the system. We also discuss some theoretical models for the structure and evolution of WDs in NS-WD binaries by using J1910A's companion as a test bed. J1910A is a rare system for which several parameters of both the pulsar and the HeWD companion can be accurately measured. As such, it is a test bed to discriminate between alternative models for HeWD structure and cooling.

Observing and timing a group of millisecond pulsars (MSPs) with high rotational stability enables the direct detection of gravitational waves (GWs). The GW signals can be identified from the spatial correlations encoded in the times-of-arrival of widely spaced pulsar-pairs. The Chinese Pulsar Timing Array (CPTA) is a collaboration aiming at the direct GW detection with observations carried out using Chinese radio telescopes. This short article serves as a `table of contents' for a forthcoming series of papers related to the CPTA Data Release 1 (CPTA DR1) which uses observations from the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Here, after summarizing the time span and accuracy of CPTA DR1, we report the key results of our statistical inference finding a correlated signal with amplitude $\log A_{\rm c}= -14.4 \,^{+1.0}_{-2.8}$ for spectral index in the range of $\alpha\in [-1.8, 1.5]$ assuming a GW background (GWB) induced quadrupolar correlation. The search for the Hellings-Downs (HD) correlation curve is also presented, where some evidence for the HD correlation has been found that a 4.6-$\sigma$ statistical significance is achieved using the discrete frequency method around the frequency of 14 nHz. We expect that the future International Pulsar Timing Array data analysis and the next CPTA data release will be more sensitive to the nHz GWB, which could verify the current results.

ABSTRACTThe radio polarization properties of the pulsar population are only superficially captured by the conventional picture of pulsar radio emission. We study the broadband polarization of 271 young radio pulsars, focusing particularly on circular polarization, using high-quality observations made with the Ultra-Wideband Low receiver on Murriyang, the Parkes radio telescope. We seek to encapsulate polarization behaviour on a population scale by defining broad categories for frequency- and phase-dependent polarization evolution, studying the co-occurrences of these categorizations and comparing them with average polarization measurements and spin-down energy ($\dot{E}$). This work shows that deviations of the linear polarization position angle from the rotating vector model are linked to the presence of circular polarization features and to frequency evolution of the polarization. Polarization fraction, circular polarization contribution, and profile complexity all evolve with $\dot{E}$ across the population, with the profiles of high-$\dot{E}$ pulsars being simple and highly linearly polarized. The relationship between polarization fraction and circular contribution is also seen to evolve such that highly polarized profiles show less variation in circular contribution with frequency than less strongly polarized profiles. This evolution is seen both across the population and across frequency for individual sources. Understanding pulsar radio polarization requires detailed study of individual sources and collective understanding of population-level trends. For the former, we provide visualizations of their phase- and frequency-resolved polarization parameters. For the latter, we have highlighted the importance of including the impact of circular polarization and of $\dot{E}$.

We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^{\alpha}$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $\alpha = -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $\alpha = -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability.

We report the discovery of a new 5.78 ms-period millisecond pulsar (MSP), PSR J1740-5340B (NGC 6397B), in an eclipsing binary system discovered with the Parkes radio telescope (now also known as Murriyang), Australia, and confirmed with the MeerKAT radio telescope in South Africa. The measured orbital period, 1.97 days, is the longest among all eclipsing binaries in globular clusters (GCs) and consistent with that of the coincident X-ray source U18, previously suggested to be a 'hidden MSP'. Our XMM-Newton observations during NGC 6397B's radio quiescent epochs detected no X-ray flares. NGC 6397B is either a transitional MSP or an eclipsing binary in its initial stage of mass transfer after the companion star left the main sequence. The discovery of NGC 6397B potentially reveals a subgroup of extremely faint and heavily obscured binary pulsars, thus providing a plausible explanation to the apparent dearth of binary neutron stars in core-collapsed GCs as well as a critical constraint on the evolution of GCs.