The ability to make independent detections of the signatures of exoplanets with complementary telescopes and instruments brings a new potential for robust identification of exoplanets and precision characterization. We introduce PEXO, a package for Precise EXOplanetology to facilitate the efficient modeling of timing, astrometry, and radial velocity data, which will benefit not only exoplanet science but also various astrophysical studies in general. PEXO is general enough to account for binary motion and stellar reflex motions induced by planetary companions and is precise enough to treat various relativistic effects both in the solar system and in the target system. We also model the post-Newtonian barycentric motion for future tests of general relativity in extrasolar systems. We benchmark PEXO with the pulsar timing package TEMPO2 and find that PEXO produces numerically similar results with timing precision of about 1 ns, space-based astrometry to a precision of 1 $\mu$as, and radial velocity of 1 $\mu$m/s and improves on TEMPO2 for decade-long timing data of nearby targets, due to its consideration of third-order terms of Roemer delay. PEXO is able to avoid the bias introduced by decoupling the target system and the solar system and to account for the atmospheric effects which set a practical limit for ground-based radial velocities close to 1 cm/s. Considering the various caveats in barycentric correction and ancillary data required to realize cm/s modeling, we recommend the preservation of original observational data. The PEXO modeling package is available at GitHub (https://github.com/phillippro/pexo).

The Alpha Centauri system is the primary target for planet search as it is the closest star system composed of a solar twin Alpha Centauri A, a K-dwarf Alpha Centauri B and an M-dwarf Proxima Centauri, which has a confirmed planet in the temperate zone. Alpha Centauri A & B were monitored intensively with the HARPS spectrograph for over 10 years, providing high-precision radial velocity measurements. In this work we study the available data to better understand the stellar activity and other contaminating signals. We highlight the importance of telluric contamination and its impact on the radial velocity measurements. Our suggested procedures lead to discarding about 5% of HARPS data, providing a dataset with an RMS improved by a factor of 2. We compile and quantify the behaviour of 345 spectral lines with a wide range of line shapes and sensitivity to activity.