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

A Bayesian parameter estimation approach to pulsar time-of-arrival analysis

149   0   0.0 ( 0 )
 نشر من قبل Chris Messenger
 تاريخ النشر 2011
  مجال البحث فيزياء
والبحث باللغة English




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

The increasing sensitivities of pulsar timing arrays to ultra-low frequency (nHz) gravitational waves promises to achieve direct gravitational wave detection within the next 5-10 years. While there are many parallel efforts being made in the improvement of telescope sensitivity, the detection of stable millisecond pulsars and the improvement of the timing software, there are reasons to believe that the methods used to accurately determine the time-of-arrival (TOA) of pulses from radio pulsars can be improved upon. More specifically, the determination of the uncertainties on these TOAs, which strongly affect the ability to detect GWs through pulsar timing, may be unreliable. We propose two Bayesian methods for the generation of pulsar TOAs starting from pulsar search-mode data and pre-folded data. These methods are applied to simulated toy-model examples and in this initial work we focus on the issue of uncertainties in the folding period. The final results of our analysis are expressed in the form of posterior probability distributions on the signal parameters (including the TOA) from a single observation.



قيم البحث

اقرأ أيضاً

A large number of binary black holes (BBHs) with longer orbital periods are supposed to exist as progenitors of BBH mergers recently discovered with gravitational wave (GW) detectors. In our previous papers, we proposed to search for such BBHs in tri ple systems through the radial-velocity modulation of the tertiary orbiting star. If the tertiary is a pulsar, high precision and cadence observations of its arrival time enable an unambiguous characterization of the pulsar -- BBH triples located at several kpc, which are inaccessible with the radial velocity of stars. The present paper shows that such inner BBHs can be identified through the short-term R{o}mer delay modulation, on the order of $10$ msec for our fiducial case, a triple consisting of $20~M_odot$ BBH and $1.4~M_odot$ pulsar with $P_mathrm{in}=10$ days and $P_mathrm{out}=100$ days. If the relativistic time delays are measured as well, one can determine basically all the orbital parameters of the triple. For instance, this method is applicable to inner BBHs of down to $sim 1$ hr orbital periods if the orbital period of the tertiary pulsar is around several days. Inner BBHs with $lesssim 1$ hr orbital period emit the GW detectable by future space-based GW missions including LISA, DECIGO, and BBO, and very short inner BBHs with sub-second orbital period can be even probed by the existing ground-based GW detectors. Therefore, our proposed methodology provides a complementary technique to search for inner BBHs in triples, if exist at all, in the near future.
We present a Bayesian parameter-estimation pipeline to measure the properties of inspiralling stellar-mass black hole binaries with LISA. Our strategy (i) is based on the coherent analysis of the three noise-orthogonal LISA data streams, (ii) employs accurate and computationally efficient post-Newtonian waveforms accounting for both spin-precession and orbital eccentricity, and (iii) relies on a nested sampling algorithm for the computation of model evidences and posterior probability density functions of the full 17 parameters describing a binary. We demonstrate the performance of this approach by analyzing the LISA Data Challenge (LDC-1) dataset, consisting of 66 quasi-circular, spin-aligned binaries with signal-to-noise ratios ranging from 3 to 14 and times to merger ranging from 3000 to 2 years. We recover 22 binaries with signal-to-noise ratio higher than 8. Their chirp masses are typically measured to better than $0.02 M_odot$ at $90%$ confidence, while the sky-location accuracy ranges from 1 to 100 square degrees. The mass ratio and the spin parameters can only be constrained for sources that merge during the mission lifetime. In addition, we report on the successful recovery of an eccentric, spin-precessing source at signal-to-noise ratio 15 for which we can measure an eccentricity of $3times 10^{-3}$.
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.
Using the latest numerical simulations of rotating stellar core collapse, we present a Bayesian framework to extract the physical information encoded in noisy gravitational wave signals. We fit Bayesian principal component regression models with know n and unknown signal arrival times to reconstruct gravitational wave signals, and subsequently fit known astrophysical parameters on the posterior means of the principal component coefficients using a linear model. We predict the ratio of rotational kinetic energy to gravitational energy of the inner core at bounce by sampling from the posterior predictive distribution, and find that these predictions are generally very close to the true parameter values, with $90%$ credible intervals $sim 0.04$ and $sim 0.06$ wide for the known and unknown arrival time models respectively. Two supervised machine learning methods are implemented to classify precollapse differential rotation, and we find that these methods discriminate rapidly rotating progenitors particularly well. We also introduce a constrained optimization approach to model selection to find an optimal number of principal components in the signal reconstruction step. Using this approach, we select 14 principal components as the most parsimonious model.
Recent core-collapse supernova (CCSN) simulations have predicted several distinct features in gravitational-wave (GW) spectrograms, including a ramp-up signature due to the g-mode oscillation of the proto-neutron star (PNS) and an excess in the low-f requency domain (100-300 Hz) potentially induced by the standing accretion shock instability (SASI). These predictions motivated us to perform a sophisticated time-frequency analysis (TFA) of the GW signals, aimed at preparation for future observations. By reanalyzing a gravitational waveform obtained in a three-dimensional general-relativistic CCSN simulation, we show that both the spectrogram with an adequate window and the quadratic TFA separate the multimodal GW signatures much more clearly compared with the previous analysis. We find that the observed low-frequency excess during the SASI active phase is divided into two components, a stronger one at 130 Hz and an overtone at 260 Hz, both of which evolve quasi-statically during the simulation time. We also identify a new mode whose frequency varies from 700 to 600 Hz. Furthermore, we develop the quadratic TFA for the Stokes I, Q, U, and V parameters as a new tool to investigate the GW circular polarization. We demonstrate that the polarization states that randomly change with time after bounce are associated with the PNS g-mode oscillation, whereas a slowly changing polarization state in the low-frequency domain is connected to the PNS core oscillation. This study demonstrates the capability of the sophisticated TFA for diagnosing the polarized CCSN GWs in order to explore their complex nature.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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