Rapid parameter estimation of a two-component neutron star model with spin wandering using a Kalman filter

Abstract

The classic, two-component, crust-superfluid model of a neutron star can be formulated as a noise-driven, linear dynamical system, in which the angular velocities of the crust and superfluid are tracked using a Kalman filter applied to electromagnetic pulse timing data and gravitational-wave data, when available. Here it is shown how to combine the marginal likelihood of the Kalman filter and nested sampling to estimate full posterior distributions of the six model parameters, extending previous analyses based on a maximum-likelihood approach. The method is tested across an astrophysically plausible parameter domain using Monte Carlo simulations. It recovers the injected parameters to ≲10 pe..