A new pulsar timing model for scalar-tensor gravity with applications to PSR J2222-0137 and pulsar-black hole binaries

Context. Scalar-tensor gravity (STG) theories are well-motivated alternatives to general relativity (GR). One class of STG theories, the Damour-Esposito-Farese (DEF) gravity, has a massless scalar field with two arbitrary coupling parameters. We are interested in this theory because, despite its simplicity, it predicts a wealth of different phenomena, such as dipolar gravitational wave emission and spontaneous scalarization of neutron stars (NSs). These phenomena of DEF gravity can be tested by timing binary radio pulsars. Aims. We aim to develop a new binary pulsar timing model DDSTG to enable more precise tests of STG theories based on a minimal set of binary parameters. The expressions for post-Keplerian (PK) parameters in DEF gravity are self-consistently incorporated into the model. The new technique takes into account all possible correlations between PK parameters naturally. Methods. Grids of physical parameters of NSs are calculated in the framework of DEF gravity for a set of 11 equations of state. The automatic Differentiation (AutoDiff) technique is employed, which aids in the calculation of gravitational form factors of NSs with higher precision than in previous works. The pulsar timing program TEMPO is selected as a framework for the realization of the DDSTG model. The implemented model is applicable to any type of pulsar companions. Results. We apply the DDSTG model to the most recently published observational data for PSR J2222-0137. The obtained limits on DEF gravity parameters for this system confirm and improve previous results. New limits are also the most reliable because DEF gravity is directly fitted to the data. We argue that future observations of PSR J2222-0137 can significantly improve the limits and that PSR-BH systems have the potential to place the tightest limits in certain areas of the DEF gravity parameter space.