ترغب بنشر مسار تعليمي؟ اضغط هنا

New advances in the Gaussian-process approach to pulsar-timing data analysis

104   0   0.0 ( 0 )
 نشر من قبل Rutger van Haasteren
 تاريخ النشر 2014
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

In this work we review the application of the theory of Gaussian processes to the modeling of noise in pulsar-timing data analysis, and we derive various useful and optimized representations for the likelihood expressions that are needed in Bayesian inference on pulsar-timing-array datasets. The resulting viewpoint and formalism lead us to two improved parameter-sampling schemes inspired by Gibbs sampling. The new schemes have vastly lower chain autocorrelation lengths than the Markov Chain Monte Carlo methods currently used in pulsar-timing data analysis, potentially speeding up Bayesian inference by orders of magnitude. The new schemes can be used for a full-noise-model analysis of the large datasets assembled by the International Pulsar Timing Array collaboration, which present a serious computational challenge to existing methods.

قيم البحث

اقرأ أيضاً

127 - Xihao Deng 2014
Gravitational wave burst is a catch-all category for signals whose durations are shorter than the observation period. We apply a method new to gravitational wave data analysis --- Bayesian non-parameterics --- to the problem of gravitational wave det ection, with an emphasis on pulsar timing array observations. In Bayesian non-parametrics, constraints are set on the function space that may be reasonably thought to characterize the range of gravitational-wave signals. This differs from the approaches currently employed or proposed, which focus on introducing parametric signal models or looking for excess power as evidence of the presence of a gravitational wave signal. Our Bayesian nonparametrics analysis method addresses two issues: (1) investigate if a gravitational wave burst is present in the data; (2) infer the sky location of the source and the duration of the burst. Compared with the popular method proposed by Finn & Lommen, our method improves in two aspects: (1) we can estimate the burst duration by adding the prior that the gravitational wave signals are smooth, while Finn & Lommen ignored this important point; (2) we perform a full Bayesian analysis by marginalizing over all possible parameters and provide robust inference on the presence of gravitational waves, while Finn & Lommen chose to optimize over parameters, which would increase false alarm risk and also underestimate the parameter uncertainties.
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 param eters. 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 detection of a stochastic gravitational-wave signal from the superposition of many inspiraling supermassive black holes with pulsar timing arrays (PTAs) is likely to occur within the next decade. With this detection will come the opportunity to l earn about the processes that drive black-hole-binary systems toward merger through their effects on the gravitational-wave spectrum. We use Bayesian methods to investigate the extent to which effects other than gravitational-wave emission can be distinguished using PTA observations. We show that, even in the absence of a detection, it is possible to place interesting constraints on these dynamical effects for conservative predictions of the population of tightly bound supermassive black-hole binaries. For instance, if we assume a relatively weak signal consistent with a low number of bound binaries and a low black-hole-mass to galaxy-mass correlation, we still find that a non-detection by a simulated array, with a sensitivity that should be reached in practice within a few years, disfavors gravitational-wave-dominated evolution with an odds ratio of $sim$30:1. Such a finding would suggest either that all existing astrophysical models for the population of tightly bound binaries are overly optimistic, or else that some dynamical effect other than gravitational-wave emission is actually dominating binary evolution even at the relatively high frequencies/small orbital separations probed by PTAs.
The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the pulsar and the Earth, enabling precise measurements of the pulsar positions, distances and velocities, and the shapes of their orbits. Here we present a timing analysis of 25 pulsars observed as part of the Parkes Pulsar Timing Array (PPTA) project over time spans of up to 24 yr. The data are from the second data release of the PPTA, which we have extended by including legacy data. We make the first detection of Shapiro delay in four Southern pulsars (PSRs J1017$-$7156, J1125$-$6014, J1545$-$4550, and J1732$-$5049), and of parallax in six pulsars. The prominent Shapiro delay of PSR J1125$-$6014 implies a neutron star mass of $M_p = 1.5 pm 0.2 M_odot$ (68% credibility interval). Measurements of both Shapiro delay and relativistic periastron advance in PSR J1600$-$3053 yield a large but uncertain pulsar mass of $M_p = 2.06^{+0.44}_{-0.41}$ M$_odot$ (68% credibility interval). We measure the distance to PSR J1909$-$3744 to a precision of 10 lyr, indicating that for gravitational wave periods over a decade, the pulsar provides a coherent baseline for pulsar timing array experiments.
Folding uncertainty in theoretical models into Bayesian parameter estimation is necessary in order to make reliable inferences. A general means of achieving this is by marginalizing over model uncertainty using a prior distribution constructed using Gaussian process regression (GPR). As an example, we apply this technique to the measurement of chirp mass using (simulated) gravitational-wave signals from binary black holes that could be observed using advanced-era gravitational-wave detectors. Unless properly accounted for, uncertainty in the gravitational-wave templates could be the dominant source of error in studies of these systems. We explain our approach in detail and provide proofs of various features of the method, including the limiting behavior for high signal-to-noise, where systematic model uncertainties dominate over noise errors. We find that the marginalized likelihood constructed via GPR offers a significant improvement in parameter estimation over the standard, uncorrected likelihood both in our simple one-dimensional study, and theoretically in general. We also examine the dependence of the method on the size of training set used in the GPR; on the form of covariance function adopted for the GPR, and on changes to the detector noise power spectral density.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا