Search

We study an inflation model with nonminimal derivative coupling that features a coupling between the derivative of the inflaton field and the Einstein tensor. This model naturally amplifies curvature perturbations at small scales via gravitationally enhanced friction, a mechanism critical for the formation of primordial black holes and the associated production of potentially detectable scalar-induced gravitational waves. We derive analytical expressions for the primordial power spectrum, enabling efficient exploration of the model parameter space without requiring computationally intensive numerical solutions of the Mukhanov-Sasaki equation. Using the third data release of the Parkes Pulsar Timing Array (PPTA DR3), we constrain the model parameters characterizing the coupling function: $\phi_c = 3.7^{+0.3}_{-0.5} M_\mathrm{P}$, $\log_{10} \omega_L = 7.1^{+0.6}_{-0.3}$, and $\log_{10} \sigma = -8.3^{+0.3}_{-0.6}$ at 90\% confidence level. Our results demonstrate the growing capability of pulsar timing arrays to probe early Universe physics, complementing traditional cosmic microwave background observations by providing unique constraints on inflationary dynamics at small scales., Comment: 14 pages, 5figures

We conduct the first-ever Pulsar Polarization Array (PPA) analysis to detect the ultralight Axion-Like 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 $\rho_0=0.4\,\textrm{GeV}/\textrm{cm}^3$, the Parkes PPA offers the best global limits on the ALDM Chern-Simons coupling, namely $\lesssim 10^{-13.5}-10^{-12.2}~{\rm GeV}^{-1}$, for the mass range of $10^{-22} - 10^{-21}~{\rm eV}$. The crucial role of pulsar cross-correlation in recognizing the nature of the derived limits is also highlighted., Comment: 6+15 pages, 10 figures, 2 tables, submitted to the journal

This paper demonstrates the impact of state-of-the-art instrumental calibration techniques on the precision of arrival times obtained from 9.6 years 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 (PPTA) 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 percent, and the maximum absolute reduction is 4.5 microseconds. For PSR J0437-4715, METM and MTM reduce the best-fit power-law amplitude (2.7 sigma) and spectral index (1.7 sigma) of the red noise in the arrival time residuals, which can 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., Comment: 30 pages, 9 figures, accepted for publication in The Astrophysical Journal

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 years 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\pm 0.044$ $M_{\odot}$, distance $D=156.96 \pm 0.11$ pc, and orbital inclination angle $i=137.506 \pm 0.016^\circ$, 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., Comment: 13 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

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 $\gamma=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission., Comment: 18 pages, 10 figures. Accepted for publication in ApJL

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 years 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$ years 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., Comment: 15 pages, 6 figures. Accepted for publication in PASA

Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 years. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum $h_c = A(f/1 {\rm yr}^{-1})^{\alpha}$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $\alpha=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $\alpha=-2/3$, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A=2.04^{+0.25}_{-0.22} \times 10^{-15}$. However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on $A$ that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with $\alpha=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2\sigma$). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array., Comment: 19 pages, 10 figures, Accepted for publication in ApJL

In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compared to typical narrow-band methods are, firstly, generation of two-dimensional templates accounting for pulse shape evolution with frequency and, secondly, simultaneous measurements of the pulse time-of-arrival (ToA) and dispersion measure (DM). This is the first time that wideband timing has been applied to a uniform data set collected with a single large-fractional bandwidth receiver, for which such techniques were originally developed. As a result of our study, we present a set of profile evolution models and new timing solutions including initial noise analysis. Precision of our ToA and DM measurements is in the range of 0.005 $-$ 2.08 $\mu$s and (0.043$-$14.24)$\times10^{-4}$ cm$^{-3}$ pc, respectively, with 94% of the pulsars achieving a median ToA uncertainty of less than 1 $\mu$s., Comment: 29 pages, 12 figures, Accepted for publication in ApJ

Composed of ultralight bosons, fuzzy dark matter provides an intriguing solution to challenges that the standard cold dark matter model encounters on sub-galactic scales. The ultralight dark matter with mass $m\sim10^{-23} \rm{eV}$ will induce a periodic oscillation in gravitational potentials with a frequency in the nanohertz band, leading to observable effects in the arrival times of radio pulses from pulsars. Unlike scalar dark matter, pulsar timing signals induced by the vector dark matter are dependent on the oscillation direction of the vector fields. In this work, we search for ultralight vector dark matter in the mass range of $[2\times 10^{-24}, 2\times 10^{-22}]{\rm{eV}}$ through its gravitational effect in the Parkes Pulsar Timing Array (PPTA) second data release. Since no statistically significant detection is made, we place $95\%$ upper limits on the local dark matter density as $\rho_{\rm{\tiny{VF}}} \lesssim 5{\rm{GeV/cm^{3}}}$ for $m\lesssim 10^{-23}{\rm{eV}}$. As no preferred direction is found for the vector dark matter, these constraints are comparable to those given by the scalar dark matter search with an earlier 12-year data set of PPTA., Comment: 9 pages, 3 figures, 2 tables; accepted for publication as a Letter in Phys. Rev. D

Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pulsar Timing Array data indicate that such a detection may arise spuriously. In this paper, we use simulated pulsar timing array datasets to further test the robustness of the inference methods for spectral and spatial correlations from a GWB. Expanding on our previous results, we find strong support (Bayes factors exceeding $10^5$) for the presence of a common-spectrum noise process in datasets where no common process is present, under a wide range of timing noise prescriptions per pulsar. We show that these results are highly sensitive to the choice of Bayesian priors on timing noise parameters, with priors that more closely match the injected distributions of timing noise parameters resulting in diminished support for a common-spectrum noise process. These results emphasize shortcomings in current methods for inferring the presence of a common-spectrum process, and imply that the detection of a common process is not a reliable precursor to detection of the GWB. Future searches for the nanohertz GWB should remain focussed on detecting spatial correlations, and make use of more tailored specifications for a common-spectrum noise process., Comment: 12 pages, 12 figures. Accepted for publication in MNRAS

Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA) collaboration that the same models used in recent analyses resulted in strong evidence for a common-spectrum process in simulations where none is present. In this work, we introduce a methodology to distinguish pulsar power spectra with the same amplitude from noise power spectra of similar but distinct amplitudes. The former is the signature of a spatially uncorrelated pulsar term of a nanohertz gravitational-wave background, whereas the latter could represent ensemble pulsar noise properties. We test the methodology on simulated data sets. We find that the reported common process in PPTA pulsars is indeed consistent with the spectral feature of a pulsar term. We recommend this methodology as one of the validity tests that the real astrophysical and cosmological backgrounds should pass, as well as for inferences about the spatially uncorrelated component of the background., Comment: 9 pages, 2 figures

The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations and predictions from simulations of small-scale structures in the cold dark matter paradigm, while remaining consistent with other observations. Because of its large de Broglie wavelength and large local occupation number within galaxies, ultralight dark matter behaves like a coherently oscillating background field with an oscillating frequency dependent on its mass. If the dark matter particle is a spin-1 dark photon, such as the $U(1)_B$ or $U(1)_{B-L}$ gauge boson, it can induce an external oscillating force and lead to displacements of test masses. Such an effect would be observable in the form of periodic variations in the arrival times of radio pulses from highly stable millisecond pulsars. In this study, we search for evidence of ultralight dark photon dark matter (DPDM) using 14-year high-precision observations of 26 pulsars collected with the Parkes Pulsar Timing Array. While no statistically significant signal is found, we place constraints on coupling constants for the $U(1)_B$ and $U(1)_{B-L}$ DPDM. Compared with other experiments, the limits on the dimensionless coupling constant $\epsilon$ achieved in our study are improved by up to two orders of magnitude when the dark photon mass is smaller than $3\times10^{-22}$~eV ($10^{-22}$~eV) for the $U(1)_{B}$ ($U(1)_{B-L}$) scenario., Comment: 12 pages, 4 figures

A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave background with the latest data set of the Parkes Pulsar Timing Array. We find no evidence for a Hellings-Downs spatial correlation as expected for a stochastic gravitational wave background. Therefore, we present constraints on first-order phase transition model parameters. Our analysis shows that pulsar timing is particularly sensitive to the low-temperature ($T \sim 1 - 100$ MeV) phase transition with a duration $(\beta/H_*)^{-1}\sim 10^{-2}-10^{-1}$ and therefore can be used to constrain the dark and QCD phase transitions., Comment: 9 pages, 4 figures, 2 tables

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 years for the bright millisecond pulsar, PSR J0437-4715. These variations are the signature of the relative transverse motions of the Earth, pulsar, and 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 \pm 0.4$ pc and $D_2 = 124 \pm 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 \pm 0.3^\circ$, and longitude of ascending node, $\Omega=206.3\pm0.4^\circ$. Using scintillation arcs for precise astrometry and orbital dynamics can be superior to modelling 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., Comment: 19 pages, 10 figures. Accepted for publication in ApJ

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities., Comment: Fixed a typo in Fig. 5; Published in Phys. Rev. D as Editor's Suggestion article, 17 pages, 7 figures, 4 tables