Bayesian estimation of non-Gaussianity in pulsar timing analysis

Abstract

We introduce a method for performing a robust Bayesian analysis of non-Gaussianity present in pulsar timing data, simultaneously with the pulsar timing model, and additional stochastic parameters such as those describing red spin noise and dispersion measure variations. The parameters used to define the presence of non-Gaussianity are zero for Gaussian processes, giving a simple method of defining the strength of non-Gaussian behaviour. We use simulations to show that assuming Gaussian statistics when the noise in the data is drawn from a non-Gaussian distribution can significantly increase the uncertainties associated with the pulsar timing model parameters. We then apply the method to the publicly available 15 year Parkes Pulsar Timing Array data release 1 dataset for the binary pulsar J0437$-$4715. In this analysis we present a significant detection of non-Gaussianity in the uncorrelated non-thermal noise, but we find that it does not yet impact the timing model or stochastic parameter estimates significantly compared to analysis performed assuming Gaussian statistics. The methods presented are, however, shown to be of immediate practical use for current European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) datasets.

arXiv comment17 pages, 9 Figures

References in corpus (18)

• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• TEMPO2, a new pulsar timing package. I: Overview
• Tests of general relativity from timing the double pulsar
• Tempo2, a new pulsar timing package. II: The timing model and precision estimates
• The Parkes Pulsar Timing Array Project
• Gravitational Waves Probe the Coalescence Rate of Massive Black Hole Binaries
• Limits on the Stochastic Gravitational Wave Background from the North American Nanohertz Observatory for Gravitational Waves
• Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data
• Assessing the Role of Spin Noise in the Precision Timing of Millisecond Pulsars
• TempoNest: A Bayesian approach to pulsar timing analysis
• Pulsar timing analysis in the presence of correlated noise
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• Limitations in timing precision due to single-pulse shape variability in millisecond pulsars
• TEMPO2, a new pulsar timing package. III: Gravitational wave simulation
• The optimal schedule for pulsar timing array observations
• A Bayesian estimate of the skewness of the Cosmic Microwave Background
• Bayesian joint estimation of non-Gaussianity and the power spectrum

Cited by in corpus (15)

• Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search
• From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release
• The neutron star mass, distance, and inclination from precision timing of the brilliant millisecond pulsar J0437$-$4715
• The NANOGrav Nine-Year Data Set: Noise Budget for Pulsar Arrival Times on Intraday Timescales
• Wide-band Profile Domain Pulsar Timing Analysis
• Taming outliers in pulsar-timing datasets with hierarchical likelihoods and Hamiltonian sampling
• Generative pulsar timing analysis
• Noise analysis of the Indian Pulsar Timing Array data release I
• Profile Stochasticity in PSR J1909-3744
• Gravitational wave non-linearities and pulsar-timing array angular correlations
• A Detection of Red Noise in PSR J1824$-$2452A and Projections for PSR B1937+21 using NICER X-ray Timing Data
• Structure in Galaxy Distribution. III. Fourier Transforming the Universe
• Improving DM estimates using low-frequency scattering-broadening estimates
• Probing Ultralight Axion-like Dark Matter: A Pulsar Timing Arrays-Pulsar Polarization Arrays Synergy
• Searching for A Generic Gravitational Wave Background via Bayesian Nonparametric Analysis with Pulsar Timing Arrays