New Doppler velocity measurements of GJ 876

The nearby M4V dwarf hosts four known planets, the outer three of which are trapped in a Laplace mean-motion resonance. A thorough characterization of the complex resonant perturbations exhibited by the orbiting planets, and the chaotic dynamics therein, is key to a complete picture of the system's formation and evolutionary history.

Here we present a reanalysis of the system using 6 yr of new radial velocity (RV) data from four instruments. These new data augment and more than double the size of the decades-long collection of existing velocity measurements. We provide updated estimates of the system parameters by employing a computationally efficient Wisdom-Holman N-body symplectic integrator, coupled with a Gaussian process (GP) regression model to account for correlated stellar noise.