MatBYIB: A Matlab-based code for Bayesian inference of extreme mass-ratio inspiral binary with arbitrary eccentricity

Abstract

Accurate parameter estimation(PE) of gravitational waves(GW) is essential for GW data analysis. In extreme mass-ratio inspiral binary(EMRI) systems, orbital eccentricity is a critical parameter for PE. However, current software for for PE of GW often neglects the direct estimation of orbital eccentricity. To fill this gap, we have developed the MatBYIB, a MATLAB-based software package for PE of GW with arbitrary eccentricity. The MatBYIB employs the Analytical Kludge (AK) waveform as a computationally efficient signal generator and computes parameter uncertainties via the Fisher Information Matrix (FIM) and the Markov Chain Monte Carlo (MCMC). For Bayesian inference, we implement the Metropolis-Hastings (M-H) algorithm to derive posterior distributions. To guarantee convergence, the Gelman-Rubin convergence criterion (the Potential Scale Reduction Factor R) is used to determine sampling adequacy, with MatBYIB dynamically increasing the sample size until R < 1.05 for all parameters. Our results demonstrate strong agreement between FIM- based predictions and full MCMC sampling. This program is user-friendly and allows for estimation of gravitational wave parameters with arbitrary eccentricity on standard personal computers. Code availability:The implementation is open-source at https://github.com/GenliangLi/MatBYIB.

11 pages, 6 figures

References in corpus (37)

• emcee: The MCMC Hammer
• Advanced Virgo: a 2nd generation interferometric gravitational wave detector
• Advanced LIGO
• Exchange Monte Carlo Method and Application to Spin Glass Simulations
• dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
• TianQin: a space-borne gravitational wave detector
• Gravitational Waves from Mergin Compact Binaries: How Accurately Can One Extract the Binary's Parameters from the Inspiral Waveform?
• Effective one-body approach to general relativistic two-body dynamics
• Frequency-domain gravitational waves from non-precessing black-hole binaries. II. A phenomenological model for the advanced detector era
• Dark matter annihilation at the galactic center
• General state space Markov chains and MCMC algorithms
• A simple model of complete precessing black-hole-binary gravitational waveforms
• The construction and use of LISA sensitivity curves
• LISA Capture Sources: Approximate Waveforms, Signal-to-Noise Ratios, and Parameter Estimation Accuracy
• Angular Resolution of the LISA Gravitational Wave Detector
• Effective-one-body model for black-hole binaries with generic mass ratios and spins
• Gravitational radiation from inspiralling compact binaries completed at the third post-Newtonian order
• Dynamic temperature selection for parallel-tempering in Markov chain Monte Carlo simulations
• Dimensional regularization of the gravitational interaction of point masses
• GW190521: orbital eccentricity and signatures of dynamical formation in a binary black hole merger signal
• PyCBC Inference: A Python-based parameter estimation toolkit for compact binary coalescence signals
• "Kludge" gravitational waveforms for a test-body orbiting a Kerr black hole
• Gravitational-Wave Inspiral of Compact Binary Systems to 7/2 Post-Newtonian Order
• Gravitational Waves from a Compact Star in a Circular, Inspiral Orbit, in the Equatorial Plane of a Massive, Spinning Black Hole, as Observed by LISA
• FastEMRIWaveforms: New tools for millihertz gravitational-wave data analysis
• Frequency domain reduced order models for gravitational waves from aligned-spin compact binaries
• A new probe of dark matter properties: gravitational waves from an intermediate mass black hole embedded in a dark matter mini-spike
• Distinguishing Between Formation Channels for Binary Black Holes with LISA
• eLISA eccentricity measurements as tracers of binary black hole formation
• The Fast and the Fiducial: Augmented kludge waveforms for detecting extreme-mass-ratio inspirals
• Constraining stellar binary black hole formation scenarios with eLISA eccentricity measurements
• Science with the TianQin observatory: Preliminary result on extreme-mass-ratio inspirals
• SWIGLAL: Python and Octave interfaces to the LALSuite gravitational-wave data analysis libraries
• First Constraints on Nuclear Coupling of Axionlike Particles from the Binary Neutron Star Gravitational Wave Event GW170817
• Probing Dark Matter Spikes via Gravitational Waves of Extreme Mass Ratio Inspirals
• Testing General Relativity with Gravitational Waves: An Overview
• EMRI_MC: A GPU-based Python code for Bayesian inference of EMRI waveforms