The HADES RV Programme with HARPS-N@TNG VIII. Gl15A: A multiple wide planetary system sculpted by binary interaction

We present 20 years of radial velocity (RV) measurements of the M1 dwarf Gl15A, combining 5 years of intensive RV monitoring with the HARPS-N spectrograph with 15 years of archival HIRES/Keck RV data. We carry out an MCMC-based analysis of the RV time series, inclusive of Gaussian Process (GP) approach to the description of stellar activity induced RV variations. Our analysis confirms the Keplerian nature and refines the orbital solution for the 11.44-day period super Earth, Gl15A\,b, reducing its amplitude to $1.68^{+0.17}_{-0.18}$ m s$^{-1}$ ($M \sin i = 3.03^{+0.46}_{-0.44}$ M$_\oplus$), and successfully models a long-term trend in the combined RV dataset in terms of a Keplerian orbit with a period around 7600 days and an amplitude of $2.5^{+1.3}_{-1.0}$ m s$^{-1}$, corresponding to a super-Neptune mass ($M \sin i = 36^{+25}_{-18}$ M$_\oplus$) planetary companion. We also discuss the present orbital configuration of Gl15A planetary system in terms of the possible outcomes of Lidov-Kozai interactions with the wide-separation companion Gl15B in a suite of detailed numerical simulations. In order to improve the results of the dynamical analysis, we derive a new orbital solution for the binary system, combining our RV measurements with astrometric data from the WDS catalogue. The eccentric Lidov-Kozai analysis shows the strong influence of Gl15B on the Gl15A planetary system, which can produce orbits compatible with the observed configuration for initial inclinations of the planetary system between $75^\circ$ and $90^\circ$, and can also enhance the eccentricity of the outer planet well above the observed value, even resulting in orbital instability, for inclinations around $0^\circ$ and $15^\circ - 30^\circ$. The Gl15A system is the multi-planet system closest to Earth, at $3.57$ pc, and hosts the longest period RV sub-jovian mass planet discovered so far.