An official website of the United States government
Abstract Over the past few decades, the measurement precision of some pulsar timing experiments has advanced from ∼10 μ s to ∼10 ns, revealing many subtle phenomena. Such high precision demands both careful data handling and sophisticated timing models to avoid systematic error. To achieve these goals, we present PINT ( P INT I s N ot T empo3 ), a high-precision Python pulsar timing data analysis package, which is hosted on GitHub and available on the Python Package Index (PyPI) as pint-pulsar . PINT is well tested, validated, object oriented, and modular, enabling interactive data analysis and providing an extensible and flexible development platform for timing applications. It utilizes well-debugged public Python packages (e.g., the N um P y and A stropy libraries) and modern software development schemes (e.g., version control and efficient development with git and GitHub) and a continually expanding test suite for improved reliability, accuracy, and reproducibility. PINT is developed and implemented without referring to, copying, or transcribing the code from other traditional pulsar timing software packages (e.g., Tempo / Tempo2 ) and therefore provides a robust tool for cross-checking timing analyses and simulating pulse arrival times. In this paper, we describe the design, use, and validation of PINT , and we compare timing results between it and Tempo and Tempo2 .
more »
« less
PINT: Maximum-likelihood Estimation of Pulsar Timing Noise Parameters
Abstract PINTis a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework withinPINTto characterize the single-pulsar noise processes present in pulsar timing data sets. This framework enables parameter estimation for both uncorrelated and correlated noise processes, as well as model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated data sets as well as a real data set of PSR B1855+09. We also describe the new features implemented inPINTsince it was first described in the literature.
more »
« less
- Award ID(s):
- 2020265
- PAR ID:
- 10579513
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- Astrophysical Journal
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 971
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 150
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Context.High-precision pulsar timing is highly dependent on the precise and accurate modelling of any effects that can potentially impact the data. In particular, effects that contain stochastic elements contribute to some level of corruption and complexity in the analysis of pulsar-timing data. It has been shown that commonly used solar wind models do not accurately account for variability in the amplitude of the solar wind on both short and long timescales. Aims.In this study, we test and validate a new, cutting-edge solar wind modelling method included in theenterprisesoftware suite (widely used for pulsar noise analysis) through extended simulations. We use it to investigate temporal variability in LOFAR data. Our model testing scheme in itself provides an invaluable asset for pulsar timing array (PTA) experiments. Since, improperly accounting for the solar wind signature in pulsar data can induce false-positive signals, it is of fundamental importance to include in any such investigations. Methods.We employed a Bayesian approach utilising a continuously varying Gaussian process to model the solar wind. It uses a spherical approximation that modulates the electron density. This method, which we refer to as a solar wind Gaussian process (SWGP), has been integrated into existing noise analysis software, specificallyenterprise. Our Validation of this model was performed through simulations. We then conduct noise analysis on eight pulsars from the LOFAR dataset, with most pulsars having a time span of ∼11 years encompassing one full solar activity cycle. Furthermore, we derived the electron densities from the dispersion measure values obtained by the SWGP model. Results.Our analysis reveals a strong correlation between the electron density at 1 AU and the ecliptic latitude (ELAT) of the pulsar. Pulsars with |ELAT|< 3° exhibit significantly higher average electron densities. Furthermore, we observed distinct temporal patterns in the electron densities in different pulsars. In particular, pulsars within |ELAT|< 3° exhibit similar temporal variations, while the electron densities of those outside this range correlate with the solar activity cycle. Notably, some pulsars exhibit sensitivity to the solar wind up to 45° away from the Sun in LOFAR data. Conclusions.The continuous variability in electron density offered in this model represents a substantial improvement over previous models, that assume a single value for piece-wise bins of time. This advancement holds promise for solar wind modelling in future International Pulsar Timing Array (IPTA) data combinations.more » « less
-
ABSTRACT The ionized interstellar medium disperses pulsar radio signals, resulting in a stochastic time-variable delay known as the dispersion measure (DM) noise. In the wideband paradigm of pulsar timing, we measure a DM together with a time of arrival from a pulsar observation to handle frequency-dependent profile evolution, interstellar scintillation, and radio frequency interference more robustly, and to reduce data volumes. In this paper, we derive a method to incorporate arbitrary models of DM variation, including Gaussian process models, in pulsar timing and noise analysis and pulsar timing array analysis. This generalizes the existing method for handling DM noise in wideband data sets.more » « less
-
Abstract Astrometry, the precise measurement of star motions, offers an alternative avenue to investigate low-frequency gravitational waves through the spatial deflection of photons, complementing pulsar timing arrays reliant on timing residuals. Upcoming data from Gaia, Theia, and Roman can not only cross-check pulsar timing array findings but also explore the uncharted frequency range bridging pulsar timing arrays and LISA. We present an analytical framework to evaluate the feasibility of detecting a gravitational wave background, considering measurement noise and the intrinsic variability of the stochastic background. Furthermore, we highlight astrometry's crucial role in uncovering key properties of the gravitational wave background, such as spectral index and chirality, employing information-matrix analysis. Finally, we simulate the emergence of quadrupolar correlations, commonly referred to as the generalized Hellings-Downs curves.more » « less
-
ABSTRACT Time-domain data sets of many varieties can be prone to statistical outliers that result from instrumental or astrophysical anomalies. These can impair searches for signals within the time series and lead to biased parameter estimation. Versatile outlier mitigation methods tuned toward multimessenger time-domain searches for supermassive binary black holes have yet to be fully explored. In an effort to perform robust outlier isolation with low computational costs, we propose a Gibbs sampling scheme. This provides structural simplicity to outlier modelling and isolation, as it requires minimal modifications to adapt to time-domain modelling scenarios with pulsar-timing array or photometric data. We robustly diagnose outliers present in simulated pulsar-timing data sets, and then further apply our methods to pulsar J1909−3744 from the NANOGrav 9-year Data set. We also explore the periodic binary-AGN candidate PG1302−102 using data sets from the Catalina Real-time Transient Survey, All-Sky Automated Survey for Supernovae, and the Lincoln Near-Earth Asteroid Research. We present our findings and outline future work that could improve outlier modelling and isolation for multimessenger time-domain searches.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.3847/1538-4357/ad59f7
