Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBayesian estimation of non-Gaussianity in pulsar timing analysis
249
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic 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.
rate research
Read More
We present a robust approach to incorporating models for the time-variable broadening of the pulse profile due to scattering in the ionized interstellar medium into profile-domain pulsar timing analysis. We use this approach to simultaneously estimate temporal variations in both the dispersion measure (DM) and scattering, together with a model for the pulse profile that includes smooth evolution as a function of frequency, and the pulsars timing model. We show that fixing the scattering timescales when forming time-of-arrival estimates, as has been suggested in the context of traditional pulsar timing analysis, can significantly underestimate the uncertainties in both DM, and the arrival time of the pulse, leading to bias in the timing parameters. We apply our method using a new, publicly available, GPU accelerated code, both to simulations, and observations of the millisecond pulsar PSR J1643$-$1224. This pulsar is known to exhibit significant scattering variability compared to typical millisecond pulsars, and we find including low-frequency ($< 1$ GHz) data without a model for these scattering variations leads to significant periodic structure in the DM, and also biases the astrometric parameters at the $4sigma$ level, for example, changing proper motion in right ascension by $0.50 pm 0.12$. If low frequency observations are to be included when significant scattering variations are present, we conclude it is necessary to not just model those variations, but also to sample the parameters that describe the variations simultaneously with all other parameters in the model, a task for which profile domain pulsar timing is ideally suited.
A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-independent method that parameterises the power at individual frequencies in the signal. We use TempoNest to show that at noise levels representative of current datasets in the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) the linear timing model can underestimate the uncertainties of the timing solution by up to an order of magnitude. We also show how to perform Bayesian model selection between different sets of timing model and stochastic parameters, for example, by demonstrating that in the pulsar B1937+21 both the dispersion measure variations and spin noise in the data are optimally modelled by simple power laws. Finally we show that not including the stochastic parameters simultaneously with the timing model can lead to unpredictable variation in the estimated uncertainties, compromising the robustness of the scientific results extracted from such analysis.
The extremely regular, periodic radio emission from millisecond pulsars makes them useful tools for studying neutron star astrophysics, general relativity, and low-frequency gravitational waves. These studies require that the observed pulse times of arrival be fit to complex timing models that describe numerous effects such as the astrometry of the source, the evolution of the pulsars spin, the presence of a binary companion, and the propagation of the pulses through the interstellar medium. In this paper, we discuss the benefits of using Bayesian inference to obtain pulsar timing solutions. These benefits include the validation of linearized least-squares model fits when they are correct, and the proper characterization of parameter uncertainties when they are not; the incorporation of prior parameter information and of models of correlated noise; and the Bayesian comparison of alternative timing models. We describe our computational setup, which combines the timing models of Tempo2 with the nested-sampling integrator MultiNest. We compare the timing solutions generated using Bayesian inference and linearized least-squares for three pulsars: B1953+29, J2317+1439, and J1640+2224, which demonstrate a variety of the benefits that we posit.
A new Bayesian method for the analysis of folded pulsar timing data is presented that allows for the simultaneous evaluation of evolution in the pulse profile in either frequency or time, along with the timing model and additional stochastic processes such as red spin noise, or dispersion measure variations. We model the pulse profiles using `shapelets - a complete ortho-normal set of basis functions that allow us to recreate any physical profile shape. Any evolution in the profiles can then be described as either an arbitrary number of independent profiles, or using some functional form. We perform simulations to compare this approach with established methods for pulsar timing analysis, and to demonstrate model selection between different evolutionary scenarios using the Bayesian evidence. %s The simplicity of our method allows for many possible extensions, such as including models for correlated noise in the pulse profile, or broadening of the pulse profiles due to scattering. As such, while it is a marked departure from standard pulsar timing analysis methods, it has clear applications for both new and current datasets, such as those from the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA).
We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the latter case we introduce a new parameterisation of the posterior that is extremely efficient in the low signal-to-noise regime and can be readily applied to a wide range of scientific problems. We apply this methodology to a series of simulations, and to between seven and nine yr of observations for PSRs J1713$+$0747, J1744$-$1134, and J1909$-$3744 with frequency coverage that spans 700-3600MHz. We use a smooth model for profile evolution across the full frequency range, and compare smooth and piecewise models for the temporal variations in DM. We find the profile domain framework consistently results in improved timing precision compared to the standard analysis paradigm by as much as 40% for timing parameters. Incorporating smoothness in the DM variations into the model further improves timing precision by as much as 30%. For PSR J1713+0747 we also detect pulse shape variation uncorrelated between epochs, which we attribute to variation intrinsic to the pulsar at a level consistent with previously published analyses. Not accounting for this shape variation biases the measured arrival times at the level of $sim$30ns, the same order of magnitude as the expected shift due to gravitational-waves in the pulsar timing band.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
