# A multi-technique approach to identifying and/or constraining radial-velocity substellar companions

F. Philipot, A-M Lagrange, Flavien Kiefer, Pascal Rubini, Philippe Delorme, A. Chomez

► To cite this version:

F. Philipot, A-M Lagrange, Flavien Kiefer, Pascal Rubini, Philippe Delorme, et al. A multi-technique approach to identifying and/or constraining radial-velocity substellar companions. Astronomy & Astrophysics - A&A, 2023, 678 (A107), 10.1051/0004-6361/202346612. hal-04465851

HAL Id: hal-04465851

https://hal.science/hal-04465851v1

Submitted on 19 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A multi-technique approach to identifying and/or constraining radial-velocity substellar companions

F. Philipot<fnref n="1" />, A.-M. Lagrange<fnref n="1" />,<fnref n="2" />, F. Kiefer<fnref n="1" />, P. Rubini<fnref n="3" />, P. Delorme<fnref n="2" />, and A. Chomez<fnref n="1" />,<fnref n="2" />

Received 6 April 2023 / Accepted 8 August 2023

## ABSTRACT

**Context.** Although more than one thousand substellar companions have already been detected with the radial velocity (RV) method, many new companions remain to be detected in the public RV archives.

**Aims.** We wish to use the archival data obtained with the ESO/HARPS spectrograph to search for substellar companions.

**Methods.** We used the astronomic acceleration measurements of stars obtained with the HIPPARCOS and Gaia satellites to identify anomalies that could be explained by the presence of a companion. Once hints for a companion were found, we combined the RV data with absolute astrometry data and, when available, relative astrometry data, using a Markov chain Monte Carlo algorithm to determine the orbital parameters and mass of the companion.

**Results.** We find and characterize three new brown dwarfs (GJ 660.1 C, HD 73256 B, and HD 165131 B) and six new planets (HD 75302 b, HD 108202 b, HD 135625 b, HD 185283 b, HIP 10337 b, and HIP 54597 b) with separations between 1 and 6 au and masses between 0.6 and 100 $M_{Jup}$. We also constrain the orbital inclination of ten known substellar companions and determine their true mass. Finally, we identify twelve new stellar companions. This shows that the analysis of proper motion anomalies enables the optimization of the RV search for substellar companions and their characterization.

**Key words.** techniques: radial velocities – techniques: high angular resolution – astrometry – planets and satellites: gaseous planets – brown dwarfs – stars: low-mass

## 1. Introduction

Since the first discovery of a planet orbiting a solar-type star (Mayor & Queloz 1995), various radial velocity (RV) surveys have detected hundreds of extrasolar planets and brown dwarfs (Mayor et al. 2011; Hébrard et al. 2016; Wittenmyer et al. 2020; Rosenthal et al. 2021). In the long term, these surveys aim to estimate the minimum mass and radial distributions of these companions and compare these values with results obtained by different population synthesis models (Mordasini 2018; Emsenhuber et al. 2021).

Despite the many detections, numerous RV data available in the public archives have never been published and could potentially contain the signals of unknown companions. The proper motion acceleration measurements published in Gaia Data Release 3 (hereafter DR3) for more than 1.4 billion stars can be used to select the targets with the highest probability of having a companion. Indeed, the study of the stars' proper motion accelerations has recently allowed the detection in high-contrast imaging (hereafter HCI) of a few substellar companions (e.g., Bonavita et al. 2022a; Swimmer et al. 2022; Franson et al. 2023; Kuzuhara et al. 2022).

Absolute astrometry based on position and proper motion measurements obtained with Gaia (Gaia Collaboration 2020) and HIPPARCOS (Perryman et al. 1997; van Leeuwen 2007) can furthermore remove the indetermination on the companion orbital inclination that is left with the RV method. For

instance, using only the RV data, Kane et al. (2019) report a planet candidate with a minimum mass of about 10 $M_{Jup}$, while Li et al. (2021), combining the RV data with absolute astrometry, found a highly inclined orbit and a true mass of about 20 $M_{Jup}$. Since the first Gaia data release, the combination of these two techniques has allowed us to improve the characterization of tens of substellar companions (Grandjean et al. 2019; Brandt et al. 2019, 2021a,b; Kiefer et al. 2019, 2021; Damasso et al. 2020; Venner et al. 2021; Li et al. 2021; Feng et al. 2022; Philipot et al. 2023).

However, the precise mass of a companion can only be determined if the orbital parameters are properly constrained, which remains difficult when the orbital period (P) is much larger than the RV baseline Lagrange et al. (2023). In the case of poorly constrained periods, the addition of relative astrometry from HCI or interferometry allows us to significantly improve the characterization of the companion. This is, for instance, the case with HD 211847 B, which was first reported using only RV data as a low-mass brown dwarf candidate with very poorly constrained orbital parameters and minimum mass (Sahlmann et al. 2011), then reported as a massive brown dwarf adding absolute astrometry (Feng et al. 2022), and finally identified as a low-mass star by adding relative astrometry data (Philipot et al. 2023).

In this paper, we study the proper motion anomalies of the stars observed with the HARPS/ESO spectrograph, with the aim of detecting new substellar companions. We also combine the available RV data and absolute astrometry with HIPPARCOS

<footnote n="1">¹ LESIA, Observatoire de Paris, Université PSL, CNRS, 5 place Jules Janssen, 92190 Meudon, France
e-mail: florian.philipot@obspm.fr</footnote>

<footnote n="2">² Univ. Grenoble Alpes, CNRS-INSU, Institut de Planétologie et d'Astrophysique de Grenoble (IPAG), UMR 5274, 38041 Grenoble, France</footnote>

<footnote n="3">³ Pixyl S.A., 5 Av. du Grand Sablon, 38700 La Tronche, France</footnote>
and Gaia measurements plus, when available, relative astrometry to further constrain the orbital and mass properties of the newly detected and the already known substellar and stellar companions identified.

In Sect. 2, we describe the selection of our sample. In Sect. 3, we present the RV, astrometric, and HCI data and the algorithm used to perform the orbital fits. We then provide the orbital parameters and masses of the new companions, as well as the update of these parameters for the known companions, in Sects. 4 and 5, respectively. In Sect. 6, we summarize the results obtained for stellar mass companions. Finally, we summarize the results and conclude in Sect. 6.

## 2. Sample

### 2.1. Proper motion anomaly for stars observed by HARPS-South

We first selected the 734 stars observed with the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph mounted on the 3.6 m ESO telescope at La Silla (Mayor et al. 2003) that are listed in the HARPS-RVBANK archive (Trifonov et al. 2020) with more than 30 RV data points spread on a temporal baseline greater than 1000 days. Secondly, we used the astrometric acceleration of stars (the so-called proper motion anomaly, Kervella et al. 2019, hereafter PMa) determined by the difference in proper motion (hereafter PM) between the average HIPPARCOS-Gaia PM (24.5 yr time baseline) and the “instantaneous” PM published in the Gaia DR3 to select systems with the highest probability of hosting a substellar companion. To do so, we selected the stars with a PMa $2\sigma$ larger than the noise. We simulated many samples of Gaia DR3 measurements for single stars of similar *G*-magnitudes and *Gb*–*Gr* colors, out of which we fit a PMa solely based on pure noise. For any of the 734 stars, we obtained PMa distributions of them as if they were single. The percentile of the PMa, measured in Kervella et al. (2022), within these distributions determines by how many $\sigma$ they exceed the noise. To focus our study on the detection of substellar companions, we then excluded systems for which the predicted mass of a companion orbiting at 10 au that could account for the measured PMa value is larger than 200 $M_{Jup}$. We obtained a sample of 116 stars.

### 2.2. Selection of systems with a substellar companion

For the selected sample of stars, we then retrieved all the RV measurements available in the public archives of 12 different spectrographs (see below) and analyzed them to identify the presence of companion(s). For this step, we used the DPASS genetic algorithm<fnref n="1" /> (Lagrange et al. 2019) to detect periodicity in the data and make a rapid estimate of the orbital properties and the minimum mass of the companion(s). We thus find 20 multiple systems (two or more companions) and 80 single systems (one companion).

Then, we selected the systems for which the orbital parameters of the companion are properly determined and the mass is lower than about $80 M_{Jup}$, corresponding to substellar companions. To do so, we combined the RV data with astrometry to determine the orbital inclination of the companion and obtain its

true mass. We used a Markov chain Monte Carlo (MCMC) algorithm to estimate the uncertainties of the different parameters. In the first step, we used the orvarra (Orbits Using Radial Velocity, Absolute, and/or Relative Astrometry) algorithm (Brandt et al. 2021d), which enables the RV data to be combined with the absolute astrometry and the available HCI data. Orvarra converges faster than most other algorithms (Brandt et al. 2021c) and enables us to quickly estimate the orbital parameters and mass of the companion and exclude stellar companions. However, Lagrange et al. (2023) have shown that, when the orbital period is not well covered by the RV data points, MCMC algorithms have difficulties exploring very eccentric solutions. Therefore, we then used our MCMC algorithm (described in Sect. 3.4) based on the MCMC sampling tool used by Lagrange et al. (2023) to enable a better exploration of the possible solutions.

In 27 cases, the coupling of the different methods does not properly constrain the orbit of companions, mainly due to the poor coverage of their periods by RV data. We also report on 25 stellar companions, for which the orbital parameters and the mass can be retrieved in Table C.3. Finally, we obtain a sample of 27 substellar companions for which it is possible to estimate all the orbital parameters and true masses.

## 3. Data available in the sample systems

### 3.1. Radial velocity data

We used the archival data obtained with the HARPS spectrograph available in the European Southern Observatory (ESO) archives<fnref n="2" />. To characterize as well as possible the different systems, we combined the HARPS data with all other publicly available RV data points. We retrieved from the literature the data obtained with the University College London Echelle Spectrograph (UCLES, Diego et al. 1990), the Lick Observatory Hamilton Echelle Spectrometer (Lick, Vogt 1987), the Automated Planet Finder (APF, Vogt et al. 2014), the Coudé Echelle Spectrometer (CES, Enard 1982) long camera (LC) and very long camera (VLC), the Magellan Inamori Kyocera Echelle spectrograph (MIKE, Bernstein et al. 2003), the Carnegie Planet Finder Spectrograph (PFS, Crane et al. 2010), the High-Resolution Spectrograph (HRS, Tull 1998), and the High-Resolution Echelle Spectrometer (HIRES, Vogt et al. 1994). We retrieved from the DACE archive<fnref n="3" /> the data obtained with the CORALIE echelle spectrograph (Queloz et al. 1999) and from the Observatoire de Haute-Provence (OHP) archive<fnref n="4" /> the data obtained with the ELODIE (Baranne et al. 1996) and SOPHIE (Perruchot et al. 2008) spectrographs. For the 53 systems with a characterized companion, the considered RV datasets are listed in Appendix B (Table B.1).

Due to a major optical fiber upgrade in 2015 (Lo Curto et al. 2015), an RV offset may be present between HARPS data obtained before and after this upgrade. We, therefore, separated these data into two different sets as if they were obtained from two different instruments. We labeled the datasets obtained before May 19, 2015 “H03” and the datasets obtained after June 3, 2015 “H15”. As for HARPS, a major upgrade has been performed on the HIRES (Tal-Or et al. 2019) and SOPHIE (Bouchy et al. 2013) spectrographs and two major upgrades have been performed on the CORALIE spectrograph (Ségransan et al. 2010). We labeled the HIRES datasets obtained

<footnote n="1">¹ DPASS uses an evolutionary algorithm to fit RV data with one or more Keplerian orbits. This algorithm converges to the solution corresponding to the smallest root-mean-square (hereafter rms) of the RV residuals.</footnote>

<footnote n="2">² http://archive.eso.org</footnote>

<footnote n="3">³ https://dace.unige.ch</footnote>

<footnote n="4">⁴ http://atlas.obs-hp.fr</footnote>
<table><caption>**Table 1.** Characteristics of the reduced HCI observations.</caption><thead><tr><th>Star</th><th>Date obs.</th><th>Instrument</th><th>Filter</th><th>NDIT×DIT(s)×Nframe (*)</th><th>Seeing (′′)</th><th>Airmass</th><th>τ<sub>0</sub> (ms)<sup>(a)</sup></th><th>Program ID</th></tr></thead><tbody><tr><td>HD 1388</td><td>2006-07-24</td><td>NACO</td><td>DB_Hsdi</td><td>40 × 4 × 1</td><td>1.74</td><td>1.03</td><td>0.0008</td><td>077.C-0293(A)</td></tr><tr><td>HD 78612</td><td>2006-01-04</td><td>NACO</td><td>DB_Hsdi</td><td>25 × 7 × 1</td><td>1.00</td><td>1.034</td><td>0.0037</td><td>076.C-0762(A)</td></tr><tr><td>HD 93351</td><td>2012-06-15</td><td>NACO</td><td>Ks</td><td>200 × 0.15 × 1</td><td>0.81</td><td>1.377</td><td>0.0018</td><td>089.C-0641(A)</td></tr><tr><td>HD 207700</td><td>2012-06-14</td><td>NACO</td><td>ND_L'</td><td>200 × 0.2 × 5</td><td>-1.00</td><td>1.517</td><td>0.0019</td><td>089.C-0641(A)</td></tr><tr><td>HD 21175</td><td>2016-10-04</td><td>SPHERE</td><td>DB_H23</td><td>16 × 8 × 2</td><td>0.75</td><td>1.090</td><td>0.0019</td><td>098.C-0739(A)</td></tr></tbody></table>

**Notes.** (*)NDIT is the number of sub-integration within a DIT. DIT is the detector integration time per frame. τ₀ corresponds to the coherence time.

before and after August 7, 2004 “Hir94” and “Hir04”, respectively, the SOPHIE datasets obtained before and after June 17, 2011 “SOPHIE” and “SOPHIE+”, respectively, and the CORALIE datasets obtained before June 2007, between June 2007 and November 2014, and after November 2014 “C98”, “C07”, and “C14”, respectively.

## 3.2. Relative astrometry data

Whenever possible, we used high-contrast imaging data available in the literature as well as new relative astrometry measurements obtained from observations in the NAOS+CONICA (NACO; Rousset et al. 2003; Lenzen et al. 2003) and the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE; Beuzit et al. 2019) archives. The observing logs for the analyzed observations are summarized in Table 1.

For the NACO observations, even when they are quite shallow, these HCI data are precious in identifying large-separation stellar companions that could mimic the astrometric signature of closer-in substellar objects. Given the low contrast and wide separation of such stellar companions, we carried out a very simple reduction, by just stacking and subtracting one to another raw data cubes acquired at two different offset positions within the same observation, which subtracts the dark and thermal background. In spite of this simple reduction, we detect five high S/N stellar companions (Fig. 1) in the systems HD 1388, HD 78612, HD 93351, and HD 207700. Since we chose to use only non-coronographic observations even when deeper coronographic data did exist, we were also able to measure very accurately the position of the host star and of the high S/N stellar companions, using the moffatian fit of their diffraction-limited cores.

For the SPHERE observations, we detect a massive stellar companion in the system HD 21175. This companion is outside the field of view of the integral field spectrograph (IFS), and hence only the data from the near-infrared imager and spectrograph (IRDIS) are used. All the pre-reduction steps (dark and background subtraction, flat-field correction, bad pixel correction, and centering) were performed on the SPHERE data center (Delorme et al. 2017) using tools provided in the ESO data reduction and handling pipeline (Pavlov et al. 2008). Due to the stellar nature of the companion, its point spread function (PSF) on the residual image provided by angular differential imaging (ADI) based algorithms using the science sequence is completely saturated and prevents any fitting of it. We thus chose to extract the astrometry of the companion using the off-axis unsaturated PSF frames of the primary star taken before and after the coronographic sequence. We created a 4D dataset containing only the PSF frames. We used a simple No-ADI algorithm embedded into the SPHERE speckle calibration tool (SPECIAL) software

<table><caption>**Table 2.** New relative astrometry.</caption><thead><tr><th>Companion</th><th>JD - 2 400 000</th><th>SEP (mas)</th><th>PA (deg)</th></tr></thead><tbody><tr><td>HD 1388 B</td><td>53 940.39</td><td>1490 ± 10</td><td>85.6 ± 1</td></tr><tr><td>HD 21175 B</td><td>57 665.35</td><td>2628.0 ± 2.5</td><td>122.75 ± 0.08</td></tr><tr><td>HD 78612 B</td><td>53 739.29</td><td>613 ± 5</td><td>272.3 ± 1</td></tr><tr><td>HD 93351 B/C</td><td>56 093.02</td><td>1684 ± 30</td><td>76.3 ± 1</td></tr><tr><td>HD 207700 B</td><td>56 092.36</td><td>1090 ± 20</td><td>284.2 ± 1</td></tr></tbody></table>

**Notes.** Since the companion of HD 93351 is binary, we considered the average separation in RA and Dec of the two companions.

(Galicher et al. 2018) installed in the SPHERE data center. The astrometry and photometry were extracted by SPECIAL from the residual stack produced by the No-ADI reduction, taking into account the neutral density filter used during the off-axis calibration frames (ND_3.5).

The separations and position angles (hereafter SEP and PA) found for the new detections can be retrieved in Table 2. For GJ 680, HD 101198, HD 111031, HD 195564, HD 221146, and HIP 113201, we used the separations and position angles published in the literature. All considered relative astrometry values are summarized in Table B.2.

## 3.3. Absolute astrometry data

We used the position and proper motion values from the HIPPARCOS-Gaia Catalog of Accelerations (HGCA) published by Brandt (2021) from HIPPARCOS and Gaia early DR3 (Gaia Collaboration 2016, 2021) measurements to determine the stellar acceleration of each selected target. For each star, the difference between the position measurements obtained by HIPPARCOS and Gaia divided by the mean interval between the two measurements (~25 yr) gives a more accurate tangential proper motion. The corresponding values are provided in Table A.1.

## 3.4. Orbit fitting

We used our MCMC algorithm presented in Philipot et al. (2023) to fit the RV and/or relative astrometry and/or absolute astrometry data, assuming Keplerian orbits. Our tool is based on the emcee 3.0 Python package (Foreman-Mackey et al. 2013). We used the HTOF package (Brandt et al. 2021d) to fit the HIPPARCOS/Gaia intermediate astrometric data. It follows large sections of the orvara code (Brandt et al. 2021d) for the
<figure><img src="image_2.png" /><figcaption>Fig. 1. Images of new stellar companions detected in HCl. The observations considered come from the NACO and SPHERE ESO archives.</figcaption></figure>

likelihood computation. In particular:

$$LLH_{\text{sep,pa}}^2 = -0.5 * \left[ \frac{\Delta_{\text{sep}}^2}{\sigma_{\text{sep}}^2(1 - \rho_{\text{sep,pa}}^2)} + \frac{\Delta_{\text{pa}}^2}{\sigma_{\text{pa}}^2(1 - \rho_{\text{sep,pa}}^2)} + \frac{2\rho_{\text{sep,pa}} * \Delta_{\text{sep}} * \Delta_{\text{pa}}}{\sigma_{\text{sep}} \sigma_{\text{pa}} (1 - \rho_{\text{sep,pa}}^2)} \right] \quad (1)$$

$$LLH_{RV}^2 = -0.5 * \left[ \frac{\Delta_{RV}^2}{\sigma_{RV}^2 + \text{jitter}^2} + \log(\sigma_{\text{err}}^2 + \text{jitter}^2) \right] \quad (2)$$

$$LLH_{HG}^2 = (\mu_{H,o} - \bar{\mu} - p\mu_H)^T C_H^{-1} (\mu_{H,o} - \bar{\mu} - p\mu_H) + (\mu_{HG,o} - \bar{\mu} - p\mu_{HG})^T C_{HG}^{-1} (\mu_{HG,o} - \bar{\mu} - p\mu_{HG}) + (\mu_{G,o} - \bar{\mu} - p\mu_G)^T C_G^{-1} (\mu_{G,o} - \bar{\mu} - p\mu_G), \quad (3)$$

where $\Delta_{\text{sep}}$, $\Delta_{\text{pa}}$, and $\Delta_{\text{RV}}$ are the difference between the measured and predicted separations, the position angle, and the RV. $\rho_{\text{sep,pa}}$ is the correlation coefficient between the separation and the position angle measurement. $\sigma_{\text{sep}}$, $\sigma_{\text{pa}}$, and $\sigma_{\text{err}}$ correspond to the uncertainty on the measured separation, the PA, and the RV, respectively. $\sigma_{\text{RV}}$ is the variance in the RV, and jitter is the stellar jitter measured. For the likelihood of the proper motions

measurements, $\mu_{x,o}$ and $\mu_x$ correspond to the predicted and observed proper motion, p is the parallax, $C_x^{-1}$ is the inverse of the covariance matrix, and $\bar{\mu}$ corresponds to the center-of-mass motion of the system's barycenter.

Unlike orvarra or other similar packages, our code uses modular arithmetic computations for angles, wherever appropriate, to alleviate artifacts at the range boundaries. To minimize multimodality effects, the MCMC code uses a mix of move functions, in particular the Differential Evolution (DE) move function suggested in the Emcee package documentation for these situations (Foreman-Mackey et al. 2013).

We considered ten free parameters for each system: the host star mass, the parallax<fnref n="5" />, the semi-major axis (a), the eccentricity (ecc), the orbital inclination (I), the longitude of the ascending node (Ω), the argument of periastron (ω), the periastron passage time ($P_{\text{time}}$)<fnref n="6" />, the companion mass, and a stellar jitter. It is important to note that in the fitting process we considered $\sqrt{\text{ecc}} \cos \omega$ and $\sqrt{\text{ecc}} \sin \omega$ instead of ecc and ω. In addition, to combine data from different instruments, we added an instrumental offset for each instrument as a free parameter of the model (see above). We considered uniform priors for all fitting parameters, except for the host star mass and the distance of the system, for which

<footnote n="5">⁵ The host star mass and the parallax considered are those reported by Kervella et al. (2022).</footnote>

<footnote n="6">⁶ The periastron passage time is calculated considering an initial time at BJD-2 454 000.</footnote>
we considered Gaussian priors, and for the orbital inclination, for which we considered a sin(I) prior.

When the orbital period of a companion is much shorter than the Gaia DR3 observing window (1038 days), the PMa does not allow us to constrain the orbital inclination. In such cases, we fit Keplerian orbits with our MCMC, considering only the RV data points. We then used GASTON (Gaia Astrometric Noise Simulation to derive Orbit incliNation; (Kiefer et al. 2019, 2021), updated to the improved Gaia DR3 accuracy, to model the DR3 astrometric excess noise (hereafter AEN) with a motion of the system's photocenter induced by the orbiting RV companion at an optimal orbital inclination. This leads to a new estimate of the true mass of the companion. The principle of GASTON is to assume the presence of this companion and (i) simulate Gaia's observations of the system with many possible orbits of the photocenter, (ii) fit out the 5-parameters linear model, as performed by Gaia, that includes position, proper motion, and parallax, (iii) calculate the AEN from the residuals, and (iv) compare the DR3 AEN to the outcomes of the simulations and perform an optimization of the orbital inclination. The GASTON method has been improved since the EDR3 release. It now accounts for the dependence of the noise intrinsic to Gaia's acquisitions and data reductions (instrumental jitter, spacecraft's attitude, and modeling) on magnitude (G-band) and color (Gb - Gr) (Lindegren et al. 2021). It also accounts for a parallax offset due to the orbital motion and produces a corrected parallax.

# 4. New companions

## 4.1. Brown dwarfs

### 4.1.1. GJ 660.1

GJ 660.1 is a binary system composed of a M1V star (GJ 660.1 A) and another low mass star of type M9.5 ± 0.3 (GJ660.1 B, Aganze et al. 2016). The projected separation between the two M dwarfs is estimated at 122 ± 9 au (Schneider et al. 2011). For this system, 59 RV data points were obtained with the HARPS spectrograph between 2009 and 2019 for GJ 660.1 A. A periodicity is clearly present in the periodogram of the data with a peak at about 615 days. Using DPASS, we find the best solution, corresponding to an additional brown dwarf candidate with a minimum mass of 16 $M_{Jup}$ and a period of 620 days (Fig. D.1).

The significant AEN in the Gaia DR3 data for GJ 660.1 of 2.633 mas (83σ larger than noise) is dominated by orbital signals with a period shorter than 100 yr ($a < 16$ au). The effect of GJ660.1 B is thus negligible, and we assume that the AEN is only due to GJ660.1 C. We used Gaston to fit the AEN, considering the RV solution found with our MCMC algorithm (Fig. 2), considering 300 000 iterations (burn-in = 50 000) and 1000 walkers, and varying the orbital inclination and the longitude of the ascending node. We find a semi-major axis of 1.12 ± 0.02 au ($P = 618.6 \pm 0.1$ d), an eccentricity of $0.65^{+0.11}_{-0.08}$, a minimum mass of $59^{+34}_{-15} M_{Jup}$, and an inclination of $87^{+39}_{-43}^\circ$. This leads to a true mass within $78^{+38}_{-16} M_{Jup}$. The posteriors obtained for each free parameter of the MCMC are reported in Table 3 and the corner plot is available in Appendix D (Fig. D.12). Note that, even if the orbital period is properly covered by the RV data, additional measurements near the periastron are necessary to better constrain the mass and eccentricity of GJ660.1 C and define the precise nature of the companion.

Looking at the periodogram of the RV data, corrected for the signal of GJ660.1 C, we observe a significant peak at about 200 days. No correlation between RV variations and bisector

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Fit of the GJ660.1 A RV measurements. The black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.</figcaption></figure>

<table><caption>**Table 3.** Summary of posteriors for new brown dwarfs.</caption><thead><tr><th>Parameter</th><th>GJ660.1 C</th><th>HD 73256 C</th><th>HD 165131 B</th></tr></thead><tbody><tr><td>a (au)</td><td>1.12 ± 0.02</td><td>3.8 ± 0.1</td><td>3.59<sup>+0.06</sup><sub>-0.07</sub></td></tr><tr><td>P (days)</td><td>618.6 ± 0.1</td><td>2690<sup>+60</sup><sub>-102</sub></td><td>2343 ± 2</td></tr><tr><td>Ecc</td><td>0.65<sup>+0.11</sup><sub>-0.08</sub></td><td>0.16 ± 0.07</td><td>0.672 ± 0.003</td></tr><tr><td>I (°)</td><td>87<sup>+39(*)</sup><sub>-43</sub></td><td>29<sup>+5</sup><sub>-3</sub> or 152<sup>+8</sup><sub>-7</sub></td><td>93<sup>+39</sup><sub>-43</sub></td></tr><tr><td>Mass ($M_{Jup}$)</td><td>78<sup>+38(*)</sup><sub>-16</sub></td><td>16 ± 1</td><td>19<sup>+2</sup><sub>-1</sub></td></tr><tr><td>Ω (°)</td><td>112<sup>+122(*)</sup><sub>-123</sub></td><td>303<sup>+5</sup><sub>-6</sub> or 269<sup>+8</sup><sub>-7</sub></td><td>141<sup>+20</sup><sub>-19</sub></td></tr><tr><td>ω (°)</td><td>188<sup>+1</sup><sub>-2</sub></td><td>300<sup>+37</sup><sub>-53</sub></td><td>3.1 ± 0.4</td></tr><tr><td>P<sub>time</sub> (days)</td><td>540<sup>+2</sup><sub>-3</sub></td><td>244<sup>+134</sup><sub>-179</sub></td><td>1588 ± 2</td></tr><tr><td>√<sub>ecc</sub> cos(ω)</td><td>-0.80<sup>+0.06</sup><sub>-0.07</sub></td><td>0.42<sup>+0.08</sup><sub>-0.12</sub></td><td>0.818<sup>+0.002</sup><sub>-0.001</sub></td></tr><tr><td>√<sub>ecc</sub> sin(ω)</td><td>-0.12<sup>+0.02</sup><sub>-0.01</sub></td><td>-0.07<sup>+0.13</sup><sub>-0.14</sub></td><td>0.044 ± 0.006</td></tr><tr><td>Jitter (m s<sup>-1</sup>)</td><td>3.4<sup>+0.6</sup><sub>-0.5</sub></td><td>23<sup>+3</sup><sub>-2</sub></td><td>5.6 ± 0.6</td></tr><tr><td>rms (m s<sup>-1</sup>)</td><td>4.6<sup>+1.6</sup><sub>-0.2</sub></td><td>24.2<sup>+3.8</sup><sub>-1.6</sub></td><td>6.3<sup>+0.2</sup><sub>-0.1</sub></td></tr><tr><td>Instrumental offset (m s<sup>-1</sup>)</td><td>H03 = -50500<sup>+141</sup><sub>-224</sub><br>H15 = -50505<sup>+142</sup><sub>-224</sub></td><td>H03 = 59 ± 6<br>C98 = 74<sup>+13</sup><sub>-16</sub></td><td>H03 = 46204 ± 1<br>H15 = 46218<sup>+1</sup><sub>-2</sub></td></tr><tr><td></td><td></td><td>Hir94 = 101 ± 13<br>Hir04 = 34 ± 9</td><td></td></tr></tbody></table>

**Notes.** The median and $1\sigma$ confidence interval are given for each parameter. (*) Values obtained with Gaston.

velocity span (BVS) variations that may explain this peak is observed. This signal could therefore correspond to the presence of an additional companion located at 0.5 au ($P = 206$ d) with a minimum mass of 0.3 $M_{Jup}$ (Fig. D.1). It should be noted that, given the mass and semi-major axis of this candidate companion, the AEN induced by it would be negligible compared with that induced by GJ 660.1 C.

### 4.1.2. HD 73256

Eighty-two RV data points were obtained for the G8IV star HD 73256, 26 RV HARPS measurements were obtained between
2003 and 2012, 38 RV CORALIE measurements were obtained between 2001 and 2003, and 18 RV HIRES measurements were obtained between 2004 and 2016. Udry et al. (2003) reported a hot Jupiter with a minimum mass of $1.87 \pm 0.49 M_{Jup}$ and a period of $2.54858 \pm 0.00016$ days, considering the CORALIE data points. Considering all the RV data points available, a periodicity appears in the periodogram of the data with a strong peak around 2.5 days corresponding to the companion identified by Udry et al. (2003). Once the planet signal is removed, a periodicity appears in the periodogram with a peak around 3000 d. Using DPASS, we find the best two-planet solution with inner planet orbital parameters and minimum mass close to those reported by Udry et al. (2003) and a new, outer giant planet candidate with a minimum mass of $10 M_{Jup}$ and a period of 2777 days (Fig. D.2). After removing both planets' signals, a periodicity is again present around 14 days, which corresponds to the rotation period of the host star (Udry et al. 2003).

As the orbital period of HD 73256 b is much smaller than the duration of Gaia DR3, the proper motion induced by the inner planet can be neglected. We thus simultaneously fit all the RV data points available, corrected for the signal of the inner planet, with the absolute astrometry (Fig. 3) with our MCMC, considering 300 000 iterations (burn-in = 50 000) and 1000 walkers. We find a semi-major axis of $3.8 \pm 0.1$ au ($P = 2690^{+60}_{-102}$ d) and an eccentricity of $0.16 \pm 0.07$, and constrained the orbital inclination of HD 73256 c to either $29^{+5}_{-3}^{\circ}$ or $152^{+8}_{-7}^{\circ}$. We find a true mass of $16 \pm 1 M_{Jup}$, corresponding to a low-mass brown dwarf. The posteriors obtained for each free parameter of the MCMC are reported in Table 3 and the corner plots are available in the appendix (Figs. D.13 and D.14).

### 4.1.3. HD 165131

Sixty-nine RV data points were obtained with the HARPS spectrograph between 2006 and 2021 for the G3/5V star HD 165131. A periodicity is clearly present in the periodogram of the data, with a peak at about 2350 days. Using DPASS, we find the best solution, corresponding to a brown dwarf candidate with a minimum mass of $18 M_{Jup}$ and a period of 2344 days (Fig. D.5).

By simultaneously fitting RV HARPS measurements and absolute astrometry with our MCMC (Fig. 4), and considering 300 000 iterations after 50 000 iterations and 1000 walkers, we characterized the orbital parameters of the companion and constrained the orbital inclination to either $70^{+11}_{-9}^{\circ}$ or $109^{+10}_{-12}^{\circ}$. We find a semi-major axis of $3.59^{+0.06}_{-0.07}$ au ($P = 2343 \pm 2$ d), an eccentricity of $0.672 \pm 0.003$ and a true mass of $19^{+2}_{-1} M_{Jup}$. The posteriors obtained for each free parameter of the MCMC are reported in Table 3 and the corner plot is available in Appendix D (Fig. D.18).

## 4.2. New planets

### 4.2.1. HD 75302

Eighty-two RV data points were obtained for this G5V star, HD 75302. 35 RV HARPS measurements were obtained between 2012 and 2018, 26 RV ELODIE measurements were obtained between 1997 and 2005, and 21 RV SOPHIE measurements were obtained between 2008 and 2018. A periodicity appears in the periodogram of the data, with a weak peak around 2000 days (Fig. D.10). Using DPASS, we find the best solution, corresponding to an eccentric planet candidate with a period of 4300 days,

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

**Fig. 3.** Fit of the orbit of HD 73256. Top: fit of the HD 73256 RV measurements after subtracting the signal of the inner planet. Bottom: fit of the HD 73256 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.

a minimum mass of $5.3 M_{Jup}$, and an eccentricity of $0.88$. However, we observe a significant correlation between RV variations and bisector velocity span (BVS) variations (Fig. D.9), indicating a possible stellar activity signal. Correcting the RV data with the linear regression between RVs and the bisector velocity span (Melo et al. 2007), we obtain a new RV dataset, for which a strong peak is present in the periodogram around 3000 days. Considering these data corrected for the stellar activity, we find the best solution for the companion corresponding to a planet that is less massive and eccentric ($m_p \sin i = 2.7 M_{Jup}$ and $e = 0.44$) but that has a similar period (Fig. D.11).

Combining the corrected RV data with absolute astrometry with our MCMC (Fig. 5), considering 300 000 iterations (burn-in = 50 000) and 1000 walkers, we find a well-constrained orbital solution for HD 75302 b, and we constrained its orbital inclination to either $29 \pm 5^{\circ}$ or $150^{+5}_{-6}^{\circ}$. We find a semi-major axis of $5.3^{+0.2}_{-0.1}$ au ($P = 4356^{+173}_{-112}$ d), an eccentricity of $0.39^{+0.10}_{-0.11}$, and a true mass of $5.4^{+0.5}_{-0.4} M_{Jup}$. The posteriors obtained for each free parameter of the MCMC are reported in Table 4 and the corner plot is available in Appendix D (Fig. D.15).
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 4.** Fit of the orbit of HD 165131. Top: fit of the HD 165131 RV measurements. Bottom: fit of the HD 165131 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.</figcaption></figure>

### 4.2.2. HD 108202

Forty-two RV data points were obtained with the HARPS spectrograph between 2005 and 2022 for the K4/5V star HD 108202. A periodicity is present in the periodogram of the data, with a peak of around 3000 days. Using DPASS, we find the best solution, corresponding to a giant planet candidate with a minimum mass of $2.4 M_{\text{Jup}}$ and a period of 2990 days (Fig. D.3).

By simultaneously fitting RV data and absolute astrometry with our MCMC (Fig. 6), and considering 200 000 iterations (burn-in = 50 000) and 1000 walkers, we characterized the orbital parameters of the companion. We find an orbital inclination of either $56_{-13}^{+19\circ}$ or $130_{-24}^{+14\circ}$, a semi-major axis of $3.7 \pm 0.1$ au ($P = 2990_{-51}^{+54}$ d), an eccentricity of $0.52_{-0.03}^{+0.04}$, and a true mass of $3.0_{-0.5}^{+0.7} M_{\text{Jup}}$. The posteriors obtained for each free parameter of the MCMC are reported in Table 4 and the corner plot is available in Appendix D (Fig. D.16).

### 4.2.3. HD 135625

Ninety-nine RV data points were obtained with the HARPS spectrograph between 2004 and 2021 for the G3IV/V star

A&A proofs:

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 5.** Fit of the orbit of HD 75302. Top: fit of the HD 75302 RV measurements. Bottom: fit of the HD 75302 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.</figcaption></figure>

HD 135625. A periodicity is clearly present in the periodogram of the data, with a strong peak at about 4500 days. Using DPASS, we find the best solution, corresponding to a giant planet candidate with a minimum mass of $2 M_{\text{Jup}}$ and a period of 4060 days (Fig. D.4).

By simultaneously fitting RV data and absolute astrometry with our MCMC (Fig. 7), and considering 200 000 iterations (burn-in = 30 000) and 1000 walkers, we characterized the orbital parameters of the companion. Absolute astrometry does not allow us to properly constrain the orbital inclination. We can only exclude the solutions corresponding to a high inclined orbit and we find an orbital inclination of $93_{-43}^{+39\circ}$. We find a semi-major axis of $5.4 \pm 0.1$ au ($P = 4055_{-44}^{+45}$ d), an eccentricity of $0.16_{-0.03}^{+0.04}$, and a true mass of $2.3_{-0.3}^{+0.8} M_{\text{Jup}}$. The posteriors obtained for each free parameter of the MCMC are reported in Table 4 and the corner plot is available in Appendix D (Fig. D.17).
<table><caption>**Table 4.** Summary of posteriors for new planets.</caption><thead><tr><th>Parameter</th><th>HD 75302 b</th><th>HD 108202 b</th><th>HD 135625 b</th><th>HD 185283 b</th><th>HIP 10337 b</th><th>HIP 54597 b</th></tr></thead><tbody><tr><td>a (au)</td><td>5.3<sup>+0.2</sup><sub>-0.1</sub></td><td>3.7 ± 0.1</td><td>5.4 ± 0.1</td><td>4.6<sup>+0.2</sup><sub>-0.1</sub></td><td>5.9 ± 0.1</td><td>4.0 ± 0.1</td></tr><tr><td>P (days)</td><td>4356<sup>+173</sup><sub>-112</sub></td><td>2990<sup>+54</sup><sub>-51</sub></td><td>4055<sup>+45</sup><sub>-44</sub></td><td>4060<sup>+148</sup><sub>-112</sub></td><td>6240<sup>+151</sup><sub>-156</sub></td><td>3274<sup>+47</sup><sub>-46</sub></td></tr><tr><td>Ecc</td><td>0.39<sup>+0.10</sup><sub>-0.11</sub></td><td>0.52<sup>+0.04</sup><sub>-0.03</sub></td><td>0.16<sup>+0.04</sup><sub>-0.03</sub></td><td>0.07<sup>+0.05</sup><sub>-0.04</sub></td><td>0.43 ± 0.05</td><td>0.03<sup>+0.03</sup><sub>-0.02</sub></td></tr><tr><td>I (°)</td><td>29 ± 5 or 150<sup>+5</sup><sub>-6</sub></td><td>56<sup>+19</sup><sub>-13</sub> or 130<sup>+14</sup><sub>-24</sub></td><td>93<sup>+39</sup><sub>-43</sub></td><td>91<sup>+30</sup><sub>-31</sub></td><td>33<sup>+5</sup><sub>-4</sub> or 146<sup>+4</sup><sub>-5</sub></td><td>62<sup>+16</sup><sub>-13</sub> or 116<sup>+13</sup><sub>-14</sub></td></tr><tr><td>Mass (M<sub>Jup</sub>)</td><td>5.4<sup>+0.5</sup><sub>-0.4</sub></td><td>3.0<sup>+0.7</sup><sub>-0.5</sub></td><td>2.3<sup>+0.8</sup><sub>-0.3</sub></td><td>1.3<sup>+0.3</sup><sub>-0.1</sub></td><td>4.8 ± 0.6</td><td>2.4<sup>+0.4</sup><sub>-0.2</sub></td></tr><tr><td>Ω (°)</td><td>211 ± 11 or 246<sup>+11</sup><sub>-10</sub></td><td>308<sup>+10</sup><sub>-9</sub> or 320 ± 8</td><td>141<sup>+20</sup><sub>-19</sub></td><td>238<sup>+20</sup><sub>-21</sub></td><td>49 ± 14 or 261 ± 14</td><td>118<sup>+31</sup><sub>-21</sub> or 204<sup>+20</sup><sub>-23</sub></td></tr><tr><td>ω (°)</td><td>307<sup>+22</sup><sub>-37</sub></td><td>78<sup>+9</sup><sub>-8</sub></td><td>89 ± 12</td><td>281 ± 38</td><td>174 ± 9</td><td>164<sup>+123</sup><sub>-103</sub></td></tr><tr><td>P<sub>time</sub> (days)</td><td>3607<sup>+160</sup><sub>-245</sub></td><td>266<sup>+67</sup><sub>-57</sub></td><td>1946<sup>+125</sup><sub>-126</sub></td><td>136<sup>+378</sup><sub>-399</sub></td><td>6193<sup>+195</sup><sub>-156</sub></td><td>482<sup>+672</sup><sub>-712</sub></td></tr><tr><td>√ecc cos(ω)</td><td>0.38<sup>+0.19</sup><sub>-0.23</sub></td><td>0.14 ± 0.11</td><td>0.01<sup>+0.08</sup><sub>-0.09</sub></td><td>0.16<sup>+0.09</sup><sub>-0.14</sub></td><td>-0.61<sup>+0.08</sup><sub>-0.05</sub></td><td>-0.08<sup>+0.13</sup><sub>-0.28</sub></td></tr><tr><td>√ecc sin(ω)</td><td>-0.46<sup>+0.15</sup><sub>-0.10</sub></td><td>0.70 ± 0.03</td><td>0.39<sup>+0.04</sup><sub>-0.05</sub></td><td>-0.18<sup>+0.16</sup><sub>-0.10</sub></td><td>0.20<sup>+0.10</sup><sub>-0.13</sub></td><td>-0.04<sup>+0.15</sup><sub>-0.12</sub></td></tr><tr><td>Jitter (m s<sup>-1</sup>)</td><td>13.9<sup>+1.5</sup><sub>-1.3</sub></td><td>6.5<sup>+1.0</sup><sub>-0.8</sub></td><td>4.2<sup>+0.4</sup><sub>-0.3</sub></td><td>3.7 ± 0.4</td><td>12 ± 1</td><td>7.3<sup>+0.7</sup><sub>-0.6</sub></td></tr><tr><td>rms (m s<sup>-1</sup>)</td><td>17.0<sup>+1.4</sup><sub>-0.8</sub></td><td>7.1<sup>+1.3</sup><sub>-0.5</sub></td><td>4.7<sup>+0.2</sup><sub>-0.1</sub></td><td>4.5<sup>+0.2</sup><sub>-0.1</sub></td><td>12.1<sup>+0.9</sup><sub>-0.4</sub></td><td>7.2 ± 0.2</td></tr><tr><td>Instrumental offset (m s<sup>-1</sup>)</td><td>H03 = 20 614 ± 6<br>H15 = 20 804<sup>+11</sup><sub>-9</sub></td><td>H03 = -5776<sup>+4</sup><sub>-3</sub><br>H15 = -5766<sup>+2</sup><sub>-3</sub></td><td>H03 = -14 961 ± 1<br>H15 = -14 949 ± 1</td><td>H03 = 51 782 ± 1<br>H15 = 51 792<sup>+1</sup><sub>-2</sub></td><td>H03 = 3203 ± 2<br>H15 = 3201 ± 1</td><td>H03 = 51 754<sup>+1</sup><sub>-2</sub><br>H15 = 51 762<sup>+1</sup><sub>-2</sub></td></tr><tr><td></td><td colspan="6">ELODIE = 12 147 ± 4<br>SOPHIE = 20 524<sup>+5</sup><sub>-6</sub><br>SOPHIE+ = 20 566<sup>+6</sup><sub>-5</sub></td></tr></tbody></table>

**Notes.** The median and $1\sigma$ confidence interval are given for each parameter.

### 4.2.4. HD 185283

Seventy-eight RV data points were obtained with the HARPS spectrograph between 2004 and 2021 for the K3V star HD 185283. A periodicity is clearly present in the periodogram of the data, with a strong peak at about 4500 days. Using DPASS, we find the best solution, corresponding to a giant planet candidate with a minimum mass of 1.5 M<sub>Jup</sub> and a period of 4000 days (Fig. D.6).

By simultaneously fitting RV data and absolute astrometry with our MCMC (Fig. 8), and considering 300 000 iterations (burn-in = 30 000) and 1000 walkers, we characterized the orbital parameters of the companion and constrained the orbital inclination to be of $91^{+30^{\circ}}_{-31^{\circ}}$. We find a semi-major axis of $4.6^{+0.2}_{-0.1}$ au ($P = 4060^{+148}_{-112}$ d), an eccentricity of $0.07^{+0.05}_{-0.04}$, and a true mass of $1.3^{+0.3}_{-0.1}$ M<sub>Jup</sub>. The posteriors obtained for each free parameter of the MCMC are reported in Table 4 and the corner plot is available in Appendix D (Fig. D.19).

### 4.2.5. HIP 10337

One hundred and eight RV data points were obtained for the K7V star HIP 10337. 90 RV HARPS measurements were obtained between 2005 and 2022 and 18 RV HIRES measurements were obtained between 2008 and 2014. A periodicity is present in the periodogram of the data, with a large peak at about 8000 days. Using DPASS, we find the best solution, corresponding to a giant planet candidate with a minimum mass of 2.9 M<sub>Jup</sub> and a period of 6400 days (Fig. D.7).

By simultaneously fitting RV data and absolute astrometry with our MCMC (Fig. 9), and considering 300 000 iterations (burn-in = 30 000) and 1000 walkers, we characterized the orbital parameters of the companion and constrained the orbital inclination to either $33^{+5^{\circ}}_{-4^{\circ}}$ or $146^{+4^{\circ}}_{-5^{\circ}}$. We find a semi-major axis of

$5.9 \pm 0.1$ au ($P = 6240^{+151}_{-156}$ d), an eccentricity of $0.43 \pm 0.05$, and a true mass of $4.8 \pm 0.6$ M<sub>Jup</sub>. The posteriors obtained for each free parameter of the MCMC are reported in Table 4 and the corner plots are available in the appendix (Figs. D.20 and D.21).

### 4.2.6. HIP 54597

Sixty-nine RV data points were obtained with HARPS between 2006 and 2022 for the K5V star HIP 54597. A periodicity is clearly present in the periodogram of the data, with a strong peak at about 3300 days. Using DPASS, we find the best solution, corresponding to a giant planet candidate with a minimum mass of 2.3 M<sub>Jup</sub> and a period of 3270 days (Fig. D.8).

By simultaneously fitting RV data and absolute astrometry with our MCMC (Fig. 10), and considering 300 000 iterations (burn-in = 30 000) and 1000 walkers, we characterized the orbital parameters of the companion and constrained the orbital inclination to either $62^{+16^{\circ}}_{-13^{\circ}}$ or $116^{+13^{\circ}}_{-14^{\circ}}$. We find a semi-major axis of $4.0 \pm 0.1$ au ($P = 3274^{+47}_{-46}$ d), an eccentricity of $0.03^{+0.03}_{-0.02}$, and a true mass of $2.4^{+0.4}_{-0.2}$ M<sub>Jup</sub>. The posteriors obtained for each free parameter of the MCMC are reported in Table 4 and the corner plot is available in Appendix D (Fig. D.22).

### 4.3. Long-term stellar activity

The long-term magnetic activity of the host star can lead to variations in the RVs over periods similar to the orbital periods of the companions detected in our study (2–15 yr, Meunier et al. 2019). In order to estimate the amplitude of these RV variations, Meunier & Lagrange (2019) simulated the impact of the magnetic activity for F-G-K type stars considering a log R'<sub>HK</sub> between –5.1 and –4.6. They show that long-term magnetic activity can induce an RV jitter between 1 m s<sup>-1</sup> for quiet stars and up to 10 m s<sup>-1</sup> for the most active stars. These RV variations
<figure><img src="image_10.png" /><figcaption>**Fig. 6.** Fit of the orbit of HD 108202. Top: fit of the HD 108202 RV measurements. Bottom: fit of the HD 108202 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 7.** Fit of the orbit of HD 135625. Top: fit of the HD 135625 RV measurements. Bottom: fit of the HD 135625 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.</figcaption></figure>

can thus significantly impact the detection of sub-Jupiter planets with a minimum mass typically lower than 0.6 $M_{\text{Jup}}$ (Lagrange et al. 2023). All the newly detected companions have a minimum mass higher than 1.5 $M_{\text{Jup}}$, corresponding to RV amplitudes higher than 40 m s⁻¹, well above the expected activity signal.

Concerning eventual short-term activity, we did not note a significant correlation between BVS and VR for any of these systems (except for HD 75302, see above). This result is consistent with the fact that these companions orbit around mature stars (> typ. 1 Gyr), the stellar jitters of which are generally smaller than those for young systems (Tran et al. 2021).

# 5. Improved parameters of known RV companions

## 5.1. Long-period companions ($P > 1000$ d)

By combining RV and absolute astrometry data, we properly characterized 14 substellar companions that have already been reported in the literature. We significantly improve the characterization and determine the exact nature of ten companions, and

we obtain similar results to those previously reported for four companions: HD 140901 c, HD 169830 c, HD 221420 B, and HIP70849 b.

### 5.1.1. HD 16905

Fifty-two RV data points were obtained with the HARPS spectrograph between 2003 and 2020 for the K3V star HD 16905. By coupling HARPS RV measurements obtained before February 2018 (Julian day: 2 458 160) with absolute astrometry, Feng et al. (2022) report an eccentric, massive giant planet with a mass of $9.1_{-1.4}^{+1.3} M_{\text{Jup}}$, an orbital inclination of $19.7_{-1.9}^{+3.3}^\circ$, a semi-major axis of $6.4_{-0.4}^{+1.6} \text{ au}$, and an eccentricity of $0.67_{-0.07}^{+0.06}$.

By fitting all the available RV data and the absolute astrometry data with our MCMC (Figs. D.23 and D.24), and considering 300 000 iterations (burn-in = 50 000) and 1000 walkers, we find significantly different solutions, with a semi-major axis of $8.8_{-0.3}^{+0.4} \text{ au}$ ($P = 10\,256_{-522}^{+618} \text{ d}$), a mass of $11.3_{-0.7}^{+0.6} M_{\text{Jup}}$, and an orbital inclination of either $43 \pm 3^\circ$ or $136 \pm 3^\circ$.
<figure><img src="image_12.png" /><figcaption>**Fig. 8.** Fit of the orbit of HD 185283. Top: fit of the HD 185283 RV measurements. Bottom: fit of the HD 185283 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and *Gaia* EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.</figcaption></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /><figcaption>**Fig. 9.** Fit of the orbit of HIP 10337. Top: fit of the HIP 10337 RV measurements. Bottom: fit of the HIP 10337 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and *Gaia* EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.</figcaption></figure>

### 5.1.2. HD 28254

Seventy-two RV data points were obtained with the HARPS spectrograph between 2003 and 2021 for the G5V star HD 28254. Naef et al. (2010) reported a very eccentric giant planet candidate with a minimum mass of $1.16_{-0.06}^{+0.10} M_{Jup}$, a semi-major axis of $2.15_{-0.05}^{+0.04}$ au, and an eccentricity of $0.81_{-0.02}^{+0.05}$ based on the RV HARPS data obtained before April 2009 (Julian day: 2 454 931).

The baseline of RV data used by Naef et al. (2010) fully covers the period of HD 28254 b. Combining all HARPS RV measurements available with absolute astrometry with our MCMC (Figs. D.25 and D.26), and considering 600 000 iterations (burn-in = 100 000) and 1000 walkers, we find a solution with a semi-major axis very close to that reported by Naef et al. (2010) but more eccentric: $a = 2.45_{-0.04}^{+0.03}$ au ($P = 1333 \pm 4$ d) and $\text{ecc} = 0.95_{-0.04}^{+0.03}$. We also estimate the orbital inclination to be at either $21_{-11}^{+38.0}$ or $162_{-27}^{+7}$ and the true mass to be between 1.5 and $6.5 M_{Jup}$. Hence, HD 28254 b is a giant planet.

### 5.1.3. HD 62364

Eighty-seven RV data points were obtained with the HARPS spectrograph between 2004 and 2021 for the F7V star HD 62364. By coupling the HARPS RV measurements obtained before February 2018 (Julian day: 2 458 160) with absolute astrometry, Feng et al. (2022) report two companions: a first, $17.4_{-1.7}^{+1.6} M_{Jup}$ brown dwarf at $18.9_{-1.1}^{+1.6}$ au and a second $24.9_{-2.9}^{+2.7} M_{Jup}$ brown dwarf at $36.98_{-2.8}^{+3.0}$ au.

By fitting the RV data considering a one-companion system with DPASS, using three years of additional HARPS RV measurements, we find a significantly different result: a single companion at approximately 6.4 au and no additional signal in the RV residuals. Moreover, when considering a two-companion system, we find an offset of about 78 m s⁻¹ between the HARPS RV measurements obtained before and after the fiber upgrade, while the average offset for a star of this spectral type is smaller than 20 m s⁻¹ (Lo Curto et al. 2015), which is in agreement with the offset of 17 m s⁻¹ found for a one-companion system. We
<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /><figcaption>**Fig. 10.** Fit of the orbit of HIP 54597. Top: fit of the HIP 54597 RV measurements. Bottom: fit of the HIP 54597 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted.</figcaption></figure>

conclude that there is no evidence for the presence of a second companion.

Combining all the RV data points available and the absolute astrometry data with our MCMC, considering a one-companion system (Figs. D.27 and D.28), 300 000 iterations (burn-in = 30 000) and 1000 walkers, we find an eccentric brown dwarf with a mass of $19_{-1}^{+2} M_{\text{Jup}}$, a semi-major axis of $6.4 \pm 0.1$ au ($5 = 5172 \pm 30$ d), an eccentricity of $0.610_{-0.005}^{+0.004}$, and an orbital inclination of either $42 \pm 2^\circ$ or $131 \pm 2^\circ$.

### 5.1.4. HD 89839

Ninety-eight RV data points were obtained with HARPS between 2004 and 2021 for the F7V star HD 89839. Moutou et al. (2011) reported a giant planet candidate based on the RV HARPS measurements obtained before June 2010 (Julian day: 2 455 353). However, this dataset does not fully cover the orbit of the companion, which is consequently poorly characterized. The reported orbital parameters are therefore poorly

constrained, with a semi-major axis of $6.8_{-2.4}^{+3.3}$ au and an eccentricity of $0.32 \pm 0.20$. The available data allow us to reasonably well constrain the amplitude of RV variations, corresponding to a minimum mass of $3.9 \pm 0.4 M_{\text{Jup}}$.

The additional data points obtained between 2010 and 2021 now allow us to properly constrain the orbital parameters of HD 89839 b and, by simultaneously fitting all the RV data points available with the absolute astrometry with our MCMC (Figs. D.29 and D.30), and considering 300 000 iterations (burn-in = 30 000) and 1000 walkers, we constrained the orbital inclination to either $53_{-9}^{+16^\circ}$ or $129_{-15}^{+9^\circ}$. Thus, we find a semi-major axis of $4.9 \pm 0.1$ au ($P = 3448 \pm 26$ d), an eccentricity of $0.20 \pm 0.02$, and a true mass of $5.2 \pm 0.8 M_{\text{Jup}}$. Hence, HD 89839 b is a giant planet.

### 5.1.5. HD 114783

Two hundred and eighty-two RV data points were obtained for the K1V star HD 114783. 91 RV HARPS measurements were obtained between 2004 and 2016, 157 RV HIRES measurements were obtained between 1998 and 2019, and 34 RV HRS measurements were obtained between 2004 and 2007. Firstly, Vogt et al. (2002) reported a giant planet candidate with a minimum mass of $1 M_{\text{Jup}}$ and a period of 501 days. Bryan et al. (2016) report an outer companion with a minimum mass of $0.611_{-0.053}^{+0.056} M_{\text{Jup}}$ and a period of $4319_{-130}^{+151}$ days. Considering the RV HIRES measurements, Rosenthal et al. (2021) report a two planets system at $1.164 \pm 0.016$ au and $4.97_{-0.11}^{+0.12}$ au ($P = 4344_{-14}^{+33}$ d), with a minimum mass of $1.033_{-0.033}^{+0.034} M_{\text{Jup}}$ and $0.660_{-0.047}^{+0.046} M_{\text{Jup}}$, respectively.

As the orbital period of HD 114783 b is much smaller than the duration of Gaia DR3, the proper motion induced by the planet can be neglected. We thus simultaneously fit all the RV data points available, corrected for the signal of the inner planet, with the absolute astrometry with our MCMC (Figs. D.31 and D.32), considering 300 000 iterations (burn-in = 100 000) and 1000 walkers. We find orbital parameters close to those reported by Rosenthal et al. (2021), with a semi-major axis of $5.0 \pm 0.1$ au ($P = 4352_{-76}^{+88}$ d) and an eccentricity of $0.05_{-0.03}^{+0.04}$. We constrained the orbital inclination of HD 114783 c to either $21_{-4}^{+7^\circ}$ or $159_{-6}^{+4^\circ}$. We find a true mass of $2.0 \pm 0.4 M_{\text{Jup}}$. Hence, HD 114783 c is a giant planet.

### 5.1.6. HD 140901

Four hundred and forty-three RV data points were obtained for the G7IV star HD 140901. 52 RV HARPS measurements were obtained between 2004 and 2007, 321 RV AAT measurements were obtained between 1998 and 2014, and 70 RV PFS measurements were obtained between 2011 and 2022. Considering all the RV data points available and the absolute astrometry, Feng et al. (2022) report a massive planet with a semi-major axis of $7.42_{-0.58}^{+0.32}$ au, an inclination of $169.4_{-55.4}^{+2.2}$ $^\circ$, a mass of $6.3_{-5.0}^{+1.3} M_{\text{Jup}}$, and an eccentricity of $0.61_{-0.04}^{+0.03}$. They also report the presence of an inner planet with a minimum mass of $0.049_{-0.008}^{+0.002} M_{\text{Jup}}$ located at $0.085_{-0.004}^{+0.003}$ au. However, the study does not detail the analysis of the system. Although a significant peak is present in the periodogram of the data around ten days, once the signal of the external planet is subtracted the amplitude of the dispersion of the RV measurements is greater than the signal induced by the potential planet. It is not possible to confirm the presence of the planet.
In any case, as the orbital period of the planet candidate is much smaller than the duration of *Gaia* DR3, the proper motion induced by this potential planet is negligible. Thus, by simultaneously fitting the same RV dataset as Feng et al. (2022) and the absolute astrometry with our MCMC (Figs. D.33 and D.34), and by considering 400 000 iterations (burn-in = 100 000) and 1000 walkers, we find a semi-major axis of $11.8_{-2.5}^{+4.1}$ au ($P = 14\,386_{-4415}^{+8099}$ d), an inclination of either $40_{-11}^{+18}^\circ$ or $138_{-21}^{+12}^\circ$, a mass of $1.8 \pm 0.5\ M_{Jup}$, and an eccentricity of $0.77_{-0.07}^{+0.06}$. This solution is within the error bars associated with the values found by Feng et al. (2022) and confirms that HD 140901 b is a giant planet.

### 5.1.7. HD 143361

One hundred and twenty RV data points were obtained for the G6V star HD 143361. 47 RV HARPS measurements were obtained between 2007 and 2019, 44 RV CORALIE measurements were obtained between 2010 and 2014, and 29 RV MIKE measurements were obtained between 1998 and 2005. Jenkins et al. (2017) report a giant planet candidate with a minimum mass of $3.48 \pm 0.24\ M_{Jup}$, a semi-major axis of $1.98 \pm 0.07$ au, and an eccentricity of $0.193 \pm 0.015$, considering the RV MIKE and CORALIE measurements and the RV HARPS measurements obtained before June 2013 (Julian day: 2 456 462).

By simultaneously fitting all the available RV data points with the absolute astrometry with our MCMC (Figs. D.35 and D.36), and considering 300 000 iterations (burn-in = 50 000) and 1000 walkers, we find orbital parameters close to those reported by Jenkins et al. (2017), with a semi-major axis of $2.05 \pm 0.03$ au ($P = 1040 \pm 1$) and an eccentricity of $0.198_{0.006}^{+0.007}$. We constrained the orbital inclination of HD 143361 b to either $50_{-12}^{+19}^\circ$ or $128_{-27}^{+18}^\circ$. Thus, we find a true mass between 3.9 and $6.3\ M_{Jup}$. Hence, HD 143361 b is a giant planet.

### 5.1.8. HD 167677

Forty-four RV data points were obtained with HARPS between 2005 and 2019 for this G5V star HD 167677. Moutou et al. (2011) reported a giant planet candidate with a minimum mass of $1.36 \pm 0.12\ M_{Jup}$, a semi-major axis of $2.90 \pm 0.12$ au, and an eccentricity of $0.17 \pm 0.07$. They considered the RV HARPS measurements obtained before October 2010 (Julian day: 2 455 495), which properly cover the period of the companion.

By simultaneously fitting all the available HARPS RV measurements with the absolute astrometry with our MCMC (Figs. D.37 and D.38), and considering 200 000 iterations (burn-in = 50 000) and 1000 walkers, we find orbital parameters close to those reported by Moutou et al. (2011), with a semi-major axis of $2.9 \pm 0.1$ au ($P = 1804 \pm 73$ d) and an eccentricity of $0.35_{-0.15}^{+0.25}$. We constrained the orbital inclination of HD 167677 b to either $20_{-6}^{+14}^\circ$ or $155_{-12}^{+6}^\circ$. Thus, we find a true mass of $3.7_{-1.1}^{+1.6} M_{Jup}$. Hence, HD 167677 b is a giant planet.

### 5.1.9. HD 169830

Two hundred and seventy-seven RV data points were obtained for the F7V star HD 169830. 86 RV HARPS measurements were obtained between 2004 and 2009, 106 RV CORALIE measurements were obtained between 1999 and 2003, and 85 RV HIRES measurements were obtained between 2000 and 2019. Firstly, Naef et al. (2001) reported a giant planet candidate with a minimum mass of $2.94 \pm 0.12\ M_{Jup}$, a semi-major axis of $0.82$ au, and an eccentricity of $0.35 \pm 0.04$, considering the RV CORALIE measurements obtained before May 2000 (Julian day: 2 451 668).

A second giant planet candidate was then identified by Mayor et al. (2004) and characterized by Rosenthal et al. (2021), considering all available RV data points. They found a minimum mass of $3.51 \pm 0.12\ M_{Jup}$, a semi-major axis of $3.283_{-0.036}^{+0.035}$ au, and an eccentricity of $0.257 \pm 0.019$ for HD 169830 c. Recently, combining all the RV data points with absolute astrometry, Feng et al. (2022) report an orbital inclination of $24.47_{-7.21}^{+12.74}^\circ$ and a true mass of $7.67_{-2.76}^{+1.94} M_{Jup}$ for HD 169830 c$^7$.

By simultaneously fitting the same RV dataset as Feng et al. (2022), corrected for the signal of HD 169830 b, and the absolute astrometry with our MCMC (Figs. D.39 and D.40), and considering 300 000 iterations (burn-in = 50 000) and 1000 walkers, we find a semi-major axis of $3.29_{-0.05}^{+0.04}$ au ($P = 1832 \pm 7$ d), an inclination of either $23_{-6}^{+12}^\circ$ or $158_{-9}^{+5}^\circ$, a mass of $8.9_{-2.6}^{+2.8} M_{Jup}$, and an eccentricity of $0.23_{-0.02}^{+0.03}$. This solution is within the error bars associated with the values found by Feng et al. (2022) and confirms that HD 169830 c is a giant planet.

### 5.1.10. HD 196050

One hundred and thirty-one RV data points were obtained for the G3V star HD 196050 A. 56 RV HARPS measurements were obtained between 2004 and 2021, 31 RV CORALIE measurements were obtained between 1999 and 2003, and 44 RV UCLES measurements were obtained between 1998 and 2005. Based on NACO high contrast imaging observations, Eggenberger et al. (2007) reported a binary companion orbiting at about 500 au from HD 196050 A and composed of an M1.5–M4.5 low-mass star with a mass of $0.29 \pm 0.02\ M_\odot$ and an M2.5–M5.5 low-mass star with a mass of $0.19 \pm 0.02\ M_\odot$. Given the separation and mass of the HD 196050 B/C pair, their movement will induce negligible variations in RV and proper motion. Jones et al. (2002) reported a giant planet candidate with a minimum mass of $2.8 \pm 0.5\ M_{Jup}$, a semi-major axis of $2.4 \pm 0.5$ au, and an eccentricity of $0.19 \pm 0.09$, considering the RV UCLES measurements obtained before June 2002 (Julian day: 2 452 455). Considering the RV CORALIE measurements, Mayor et al. (2004) reported orbital parameters and a minimum mass for HD 196050 Ab close to those reported by Jones et al. (2002).

By simultaneously fitting all the available RV data points with the absolute astrometry with our MCMC (Figs. D.41 and D.42), and considering 300 000 iterations (burn-in = 30 000) and 1000 walkers, we improved the orbital parameters with a semi-major axis of $2.65 \pm 0.03$ au ($P = 1392 \pm 4$ d) and an eccentricity of $0.20 \pm 0.01$. We constrained the orbital inclination of HD 196050 Ab to either $42_{-7}^{+11}^\circ$ or $141_{-12}^{+7}^\circ$. Thus, we find a semi-major axis of $2.65 \pm 0.03$ au, an eccentricity of $0.20 \pm 0.01$, and a true mass of $4.7 \pm 0.8\ M_{Jup}$. Hence, HD 196050 Ab is a giant planet.

### 5.1.11. HD 204961

Two hundred and thirty-eight RV data points were obtained for the M 2/3V star HD 204961. 183 RV HARPS measurements were obtained between 2003 and 2019, 39 RV AAT measurements were obtained between 1998 and 2013, and 16 RV PFS measurements were obtained between 2011 and 2013. Bailey et al. (2009) reported a giant planet candidate with a minimum mass of $0.64 \pm 0.06\ M_{Jup}$ and a period of $3416 \pm 131$ days. An inner super-Earth was reported by Wittenmyer et al. (2014) with a minimum mass of $5.4 \pm 1\ M_{\text{Earth}}$; however, the signal

<footnote>$^7$ As the orbital period of HD 169830 b is much smaller than the duration of *Gaia* DR3, the proper motion induced by this planet is negligible.</footnote>
corresponding to this potential inner planet was then identified as corresponding to the rotation period of the host star by Gorrini et al. (2022).

By simultaneously fitting all the available RV data points with the absolute astrometry with our MCMC (Figs. D.43 and D.44), and considering 300 000 iterations (burn-in = 30 000) and 1000 walkers, we find orbital parameters close to those reported by Gorrini et al. (2022), with a semi-major axis of 3.7 ± 0.1 au ($P = 3853_{-0.47}^{+51}$ d) and an eccentricity of 0.05 ± 0.03. We constrained the orbital inclination of HD 204961 b to either 51 ± 3° or 134 ± 3°. We find a true mass of $0.99_{-0.08}^{+0.09} M_{\text{Jup}}$. Hence, HD 204961 b is a giant planet.

### 5.1.12. HD 216437

One hundred and sixty RV data points were obtained for the G1V star HD 216437. 100 RV HARPS measurements were obtained between 2007 and 2021, 31 RV CORALIE measurements were obtained between 1999 and 2003, and 39 RV UCLES measurements were obtained between 1998 and 2005. Jones et al. (2002) reported a giant planet candidate with a minimum mass of 2.1 ± 0.3 $M_{\text{Jup}}$, a semi-major axis of 2.4 ± 0.5 au, and an eccentricity of 0.33 ± 0.09, considering the RV UCLES measurements obtained before June 2002 (Julian day: 2 452 456). Considering the RV CORALIE measurements, Mayor et al. (2004) reported orbital parameters and minimum mass for HD 216437 b close to those reported by Jones et al. (2002).

By simultaneously fitting all the available RV data points with the absolute astrometry with our MCMC (Figs. D.45 and D.46), and considering 400 000 iterations (burn-in = 100 000) and 1000 walkers, we improved the orbital parameters and constrained the orbital inclination of HD 216437 b to either $36_{-6}^{+8°}$ or $148_{-7}^{+5°}$. Thus, we find a semi-major axis of $2.57 \pm 0.04$ au ($P = 1359 \pm 1$ d), an eccentricity of $0.31 \pm 0.01$, and a true mass of $4.2 \pm 0.7 M_{\text{Jup}}$. Hence, HD 216437 b is a giant planet.

### 5.1.13. HD 221420

Two hundred and three RV data points were obtained for the G1V star HD 221420. 77 RV HARPS measurements were obtained between 2003 and 2018, 88 RV AAT measurements were obtained between 1998 and 2015, and 38 RV PFS measurements were obtained between 2012 and 2021. Kane et al. (2019) report a giant planet candidate with a minimum mass of $9.7_{-1.0}^{+1.1} M_{\text{Jup}}$, a semi-major axis of $18.5 \pm 2.3$ au, and an eccentricity of $0.42_{-0.07}^{+0.05}$, considering the RV AAT measurements. Recently, combining the RV AAT and HARPS measurements with absolute astrometry, Li et al. (2021) report an orbital inclination of either $17.8_{-2.8}^{+2.9°}$ or $162.2_{-2.9}^{+2.8°}$ and a true mass of $20.6_{-1.6}^{+2.0}$ $M_{\text{Jup}}$, corresponding to a brown dwarf. Moreover, they obtain a better constraint on the other orbital parameters and find a semi-major axis of $9.99_{-0.70}^{+0.74}$ au and an eccentricity of $0.162_{-0.030}^{+0.035}$.

By simultaneously fitting all the available RV data points with the absolute astrometry with our MCMC (Figs. D.47 and D.48), and considering 300 000 iterations (burn-in = 50 000) and 1000 walkers, we find a semi-major axis of $10.1 \pm 0.7$ au ($P = 10\,066_{-1125}^{+929}$ d), an inclination of either $14 \pm 1°$ or $162 \pm 2°$, a mass of $22 \pm 2 M_{\text{Jup}}$, and an eccentricity of $0.12_{-0.03}^{+0.04}$. This solution is within the error bars associated with the values found by Li et al. (2021) and confirms that HD 221420 B is a brown dwarf.

### 5.1.14. HIP 70849

Fifty-nine RV data points were obtained with HARPS between 2006 and 2021 for the K7V star HIP 70849.

Ségransan et al. (2011) reported the detection of a very poorly constrained companion with a minimum mass between 3 and 15 $M_{\text{Jup}}$, a semi-major axis between 4.5 and 36 au, and an eccentricity between 0.47 and 0.96 (3σ confidence interval), considering the RV HARPS measurements obtained before June 2010 (Julian day: 2 455 374). In Philipot et al. (2023), we already used our MCMC (Fig. D.49) to combine available HARPS data points with absolute astrometry, considering 300 000 iterations (burn-in = 30 000) and 1000 walkers. We found a semi-major axis of $3.99_{-0.07}^{+0.06}$ au ($P = 3649 \pm 18$ d), an inclination of $96 \pm 16°$, a mass of $4.5_{-0.3}^{+0.4} M_{\text{Jup}}$, and an eccentricity of $0.65_{-0.01}^{+0.02}$, confirming that HIP70849 b is a giant planet. In this study, as we consider the same dataset and the same fitting tools as Philipot et al. (2023), we obtain the same solution.

## 5.2. Short-period companions ($P < 1000$ d)

Four companions (HD 44219 b, HD 100777 b, HD 111998 b, and HD 215497 c) have an orbital period smaller than the duration of Gaia data as delivered in the DR3. It is therefore not possible to precisely determine their inclination and hence their true mass from the variations of proper motion. We therefore fit only the radial velocities with *I* fixed at 90° and *ω* fixed at 0° to obtain the other orbital parameters and the minimum mass of the four companions. We then used GASTON to estimate the maximum orbital inclination and the maximum mass, which can explain the AEN measured by the Gaia DR3.

### 5.2.1. HD 44219

Sixty-five RV data points were obtained with HARPS between 2003 and 2021 for the G3V star HD 44219. Naef et al. (2010) reported a giant planet candidate with a minimum mass of $0.58_{-0.04}^{+0.06} M_{\text{Jup}}$, a semi-major axis of $1.19 \pm 0.02$ au, and an eccentricity of $0.61_{-0.09}^{+0.07}$, considering the HARPS data points obtained before April 2009 (Julian day: 2 454 932).

By fitting all the HARPS data points available with our MCMC, and considering 200 000 iterations (burn-in = 50 000) and 1000 walkers (Fig. D.50), we find orbital parameters close to those reported by Naef et al. (2010), with a semi-major axis of $1.26_{-0.07}^{+0.02}$ au ($P = 472_{-4}^{+3}$ d) and an eccentricity of $0.41_{-0.14}^{+0.12}$. Using Gaston, we could estimate that the orbital inclination of HD 44219 b is between 4° and 176° (confidence interval at 3σ), and thus that the maximum mass of the planet is $11 M_{\text{Jup}}$. Hence, HD 44219 b is a giant planet.

### 5.2.2. HD 100777

Thirty-seven RV data points were obtained with HARPS between 2004 and 2007 for the G8V star HD 100777. Naef et al. (2007) reported a giant planet candidate with a minimum mass of $1.16 \pm 0.03 M_{\text{Jup}}$, a semi-major axis of $1.03 \pm 0.03$ au, and an eccentricity of $0.36 \pm 0.02$, considering all the available HARPS data.

By fitting the same dataset as the one considered by Naef et al. (2007) with our MCMC, and considering 200 000 iterations (burn-in = 50 000) and 1000 walkers (Fig. D.50), we find orbital parameters close to those reported by Naef et al. (2007), with a semi-major axis of $1.06_{-0.01}^{+0.02}$ au ($P = 383 \pm 1$ d) and an eccentricity of $0.38 \pm 0.02$. Using Gaston, we could estimate that the orbital inclination of HD 100777 b is between 4° and 176° (confidence interval at 3σ), and thus that the maximum mass of the planet is $18 M_{\text{Jup}}$. Hence, the use of the AEN measured by Gaia does not allow us to determine with certainty the nature of
the companion. Moreover, it is important to note that the PMa does not seem to agree with that induced by HD 100777 b. This difference could be explained by the presence of an external companion. However, there is no significant signal present in the RV data. This could be due to the short RV temporal baseline (~1200 d).

### 5.2.3. HD 111998

One hundred and fifty RV data points were obtained with the HARPS spectrograph between 2006 and 2014 for HD 111998. As HD 111998 is a fast-rotating F6V star, the HARPS DRS is not able to properly compute the radial velocity of the star. For this reason, Borgniet et al. (2017) used the SAFIR (Spectroscopic data via Analysis of the Fourier Interspectrum Radial velocities) software (Galland et al. 2005) to compute the RV. They report a giant planet candidate with a minimum mass of $4.51 \pm 0.50 M_{Jup}$, a semi-major axis of $1.82 \pm 0.07$ au, and an eccentricity of $0.03 \pm 0.04$.

By fitting the same dataset as the one considered by Borgniet et al. (2017) with our MCMC, and considering 300 000 iterations (burn-in = 100 000) and 1000 walkers (Fig. D.50), we find orbital parameters close to those reported by Borgniet et al. (2017), with a semi-major axis of $1.91^{+0.02}_{-0.03}$ au ($P = 810^{+6}_{-5}$ d) and an eccentricity of $0.13^{+0.05}_{-0.04}$. Using Gaston, we could estimate that the orbital inclination of HD 111998 b is between 23° and 154° (confidence interval at $3\sigma$), and thus that the maximum mass of the planet is $12.8 M_{Jup}$. Hence, HD 111998 is a giant planet. Moreover, it is important to note that, as for HD 100777, the PMa measured for HD 111998 does not seem to agree with that induced by the HD 111998 b. This difference could be due to an additional companion, which could explain the long-term trend reported in Borgniet et al. (2017).

### 5.2.4. HD 215497

One hundred and fifteen RV data points were obtained with HARPS between 2004 and 2009 for the K3V star HD 215497. Lo Curto et al. (2010) reported a super-Earth candidate with a minimum mass of $6.6 M_{Earth}$, a semi-major axis of $0.047$ au, and an eccentricity of $0.16 \pm 0.09$ and a giant planet candidate with a minimum mass of $0.33 M_{Jup}$, a semi-major axis of $1.28$ au, and an eccentricity of $0.49 \pm 0.04$, considering all the available HARPS data points.

HD 215467 b is not massive enough to induce a significant PMa or AEN. We have therefore considered only HD 215497 c. By fitting the HARPS data points, corrected for the signal of HD 215497 b, with our MCMC, and considering 400 000 iterations (burn-in = 100 000) and 1000 walkers (Fig. D.50), we find a semi-major axis of $1.31 \pm 0.02$ au ($P = 566 \pm 4$ d), a minimum mass of $0.35 \pm 0.02 M_{Jup}$, and an eccentricity of $0.47^{+0.03}_{-0.04}$. The solution found by Lo Curto et al. (2010) is within the error bars associated with our values. Using Gaston, we could estimate that the orbital inclination of HD 215467 c is between 4° and 176° (confidence interval at $3\sigma$), and thus that the maximum mass of the planet is $5.0 M_{Jup}$. Hence, HD 215497 c is a giant planet.

## 6. Stellar companions

Among the 116 selected stars, we detect and confirm the nature of 27 stellar companions, including 12 new detections (Table C.3). These companions orbit between 0.56 and 90 au and have masses between 0.1 and $1.1 M_{\odot}$ (uncertainties included).

Even if the addition of absolute astrometry can allow us to lift the $\sin(I)$ indetermination and determine the true mass of a companion, the fit of the variations of the proper motion of a star does not always allow us to choose between a prograde and a retrograde orbit for the companion. For almost all the systems presented in this study, when the orbital period is properly constrained we find very close semi-major axes and masses for the prograde and retrograde solutions. One exception is the system HIP 39470, for which the 32 RV data obtained with HARPS between 2004 and 2019 do not fully cover the orbital period of the companion. Thanks to the significant PMa of the star, adding the absolute astrometry allows us to significantly improve the constraints on the companion's orbital parameters compared to the RV alone. However, considering a prograde and a retrograde solution ($I = 22^{+2}_{-3}$° and $I = 170 \pm 1$° respectively), we find two different solutions, with a semi-major axis of $8.8 \pm 0.2$ or $7.6 \pm 0.1$ au and a mass of $0.138^{+0.007}_{-0.008}$ or $0.147 \pm 0.005 M_{\odot}$, respectively. In the two cases, we find that the new companion HIP 39470 B is a stellar companion. The epoch astrometry of Gaia, expected with the fourth data release (at the end of 2025), or HCI observation, when possible, could constrain the exact orbital inclination of the companions and significantly improve the characterization of HIP 39470 B.

The study of these stellar mass companions shows us the importance of properly determining the orbital inclination of a companion. Indeed, using only RV data, the minimum masses found for three of these companions were in the brown dwarfs' regime. Indeed, we find a minimum mass of about 37 and 52 $M_{Jup}$ for the new companions HD 56380 B and HD 221638 B, and Moutou et al. (2009) reported a minimum mass of about 18 $M_{Jup}$ for HD 131664 B. However, coupling the RV data with absolute astrometry, we find an inclination of either $6.1 \pm 0.1°$ or $173.6 \pm 0.1°$, either $31 \pm 1°$ or $150 \pm 1°$, and either $9.3 \pm 0.3°$ or $170.4 \pm 0.3°$, respectively, corresponding to a true mass of $0.39 \pm 0.02$, $0.132 \pm 0.006$, and $0.109 \pm 0.005 M_{\odot}$ ($410 \pm 19$, $138 \pm 6$, and $114 \pm 5 M_{Jup}$), respectively. These three companions are thus, in reality, low-mass stars. We note that Feng et al. (2021) had already identified HD 131664 B as a stellar companion by coupling RV HARPS measurements with absolute astrometry.

Moreover, in most cases, when the period of the companion is much larger than the available RV baseline, the orbital parameters and the mass of the companion cannot be properly characterized. In these cases, the addition of relative astrometry significantly improves the constraint on the orbit of the companion and, consequently, its mass. Good examples are the cases of GJ680 B and HD 111031 B. Feng et al. (2022) combine RV data points obtained with HARPS between 2004 and 2012 for GJ680, and with HIRES between 1997 and 2016 for HD 111031, and absolute astrometry. They report a semi-major axis of $10.1^{+1.8}_{-1.7}$ and $13.1^{+0.8}_{-1.1}$ au and a true mass of $25.1^{+6.2}_{-11.1}$ and $54.2^{+5.3}_{-6.1} M_{Jup}$, respectively, corresponding to brown dwarfs. However, GJ680 B was detected in HCI in May 2005 and May 2007 with the NACO instrument (Ward-Duong et al. 2015) and HD 111031 B was detected in HCI in February 2018 with the NESSI instrument (Dalba et al. 2021). In this study, we consider the same RV dataset for GJ680 and, for HD 111031, we consider the RV HARPS measurements obtained between 2004 and 2013, the RV HIRES data obtained between 1997 and 2020, the RV APF measurements obtained between 2013 and 2020, and the RV LICK measurements obtained in 2001. Combining these RV data with absolute astrometry and relative astrometry data, we find significantly different results from those published by Feng et al. (2022) for the two companions. We find semi-major
<figure><img src="image_17.png" /><figcaption>**Fig. 11.** Semi-major axis and masses of the different substellar companions. For known companions, the dotted black lines enable the solutions obtained in this study, represented by a star, to be compared with those obtained in previous studies, represented by a square. The orange arrow corresponds to the maximum mass of the planet at $3\sigma$ found with Gaston. The dotted gray horizontal line corresponds to the deuterium burning limit (~13 $M_{Jup}$).</figcaption></figure>

axes of $32_{-6}^{+9}$ and $21.1 \pm 0.6$ au and true masses of $0.178 \pm 0.004$ and $0.129 \pm 0.003 M_{\odot}$ ($186 \pm 4$ and $135 \pm 3 M_{Jup}$), respectively, corresponding to low-mass stars instead of a brown dwarf.

As seen above, the algorithm used for this study, based on the MCMC sampling tool used by Lagrange et al. (2023), allows us to obtain a better exploration of the mass and orbital values compatible with the data. For example, in the case of GJ680 B, the orbit of which is poorly constrained, we obtain, by considering only RV and absolute astrometry data, solutions corresponding to a $1\sigma$ (resp. $3\sigma$) confidence interval for the semi-major axis and the companion mass between 16 and 28 au (resp. 12 and 35 au) and between 70 and 160 $M_{Jup}$ (resp. 50 and 190 $M_{Jup}$), respectively. We notice that the solutions found by adding the HCI observations are well explored by the MCMC, but it is not possible without the imaging to come to any conclusions about the orbital parameters and the mass of the companion.

Finally, thanks to the HCI observations of HD 93351 obtained with NACO, we find that the companion is actually a stellar binary with a combined mass of $0.70 \pm 0.05 M_{\odot}$. However, the combination of RV, HCI, and absolute astrometry data does not allow us to properly determine the orbital parameters of the binary. We can only obtain an estimate for the semi-major axis, with a solution at $83_{-19}^{+18}$ au.

## 7. Conclusion

Using the PMa of the stars as given by Gaia DR3 and HIPPARCOS, observed with HARPS RV data points spread over at least 1000 days, we identified good candidates to host long-period substellar companions, which allowed us to optimize our

search for substellar companions. We detect at least one companion, stellar or substellar, for 86% of the selected stars, including systems hosting multiple companions as well as those with companion's orbit poorly constrained. Moreover, 23% of the systems are composed of at least one confirmed substellar companion (41% considering companions detected in multiple systems with a minimum mass smaller than $80 M_{Jup}$). This is significantly higher than the detection rates obtained in the case of the large RV surveys (Mayor et al. 2011; Wittenmyer et al. 2020; Rosenthal et al. 2021), which do not use selection criteria to optimize the detection of substellar companions.

Among our sample, we have detected three new brown dwarfs (GJ 660.1 C, HD 73256 C, and HD 165131 B) and six new planets (HD 75302 b, HD 108202 b, HD 135625 b, HD 185283 b, HIP 10337 b, and HIP 54597 b). As already demonstrated in the direct imaging study (Bonavita et al. 2022a, Swimmer et al. 2022, Kuzuhara et al. 2022), we confirm that the study of the PMa enables the selection of targets for the search of substellar companions to be optimized. However, this method remains sensitive only to massive objects (typically $>1 M_{Jup}$) and is not adapted for the detection of low-mass planets (super-Earths and mini-Neptunes).

Combining RV data with absolute astrometry data allows us to properly characterize the orbital parameters of these ten new companions and in particular their orbital inclination. We thus determine their true mass and their exact nature. Moreover, by applying the same method, we are able to improve the characterization and determine precisely the true mass of ten known substellar companions and, using GASTON, we are able to confirm the planetary nature of three known candidate planets (Fig. 11). We determined the properties of 25 stellar companions,
including 12 new detections. We show also that five companions previously reported as brown dwarfs are indeed low-mass stars: GJ 680 B, HD 56380 B, HD 111031 B, HD 131664 B, and HD 221638 B.

The coupling of HIPPARCOS/Gaia measurements with RV data does not always permit the accurate characterization of the orbit of a companion, in particular when its period is not completely covered by the RV data. The addition of HCl data allows for better constraining of the orbit of the companion and, consequently, measuring its mass. However, most of the companions revealed in this study have not been detected or observed in imaging. Thus, neither the orbital parameters nor the exact nature of 27 of the 80 single companions found in our sample can be determined.

In summary, we conclude that studying the PMs of stars allows the search for substellar companions in RV to be optimized so that they can then be properly characterized by combining RV, absolute astrometry, and, when available, relative astrometry data. However, this method is not sensitive to low-mass planets and therefore does not enable Earth-like planets to be searched for.

Acknowledgements. This study was funded by a grant from PSL/OCAV. This project has also received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (COBREX; grant agreement n° 885593). This work presents results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia MultiLateral Agreement (MLA). The Gaia mission website is https://www.cosmos.esa.int/Gaia. The Gaia archive website is https://archives.esac.esa.int/Gaia. This publication makes use of The Data & Analysis Center for Exoplanets (DACE), which is a facility based at the University of Geneva (CH) dedicated to extrasolar planet data visualisation, exchange, and analysis. DACE is a platform of the Swiss National Centre of Competence in Research (NCCR) PlanetS, federating Swiss expertise in exoplanet research. The DACE platform is available at ttps://dace.unige.ch. This work has made use of the SPHERE Data Center, jointly operated by OSUG/IPAG (Grenoble), PYTHEAS/LAM/CeSAM (Marseille), OCA/Lagrange (Nice), Observatoire de Paris/LESIA (Paris), and Observatoire de Lyon (OSUL/CRAL). This research has made use of the SIMBAD database and VizieR catalogue access tool, operated at CDS in Strasbourg, France. It is also based on data retrieved from the SOPHIE archive at Observatoire de Haute-Provence (OHP), available at atlas.-obs-hp.fr/sophie, on spectral data retrieved from the ELODIE archive at Observatoire de Haute-Provence (OHP), and on observations collected at the European Southern Observatory under ESO programmes. We finally thank our anonymous referee, whose feedback helped us to improve the quality of this study.

## References

Adibekyan, V., Santos, N. C., Demangeon, O. D. S., et al. 2021, A&A, 649, A111

Aganze, C., Burgasser, A. J., Faherty, J. K., et al. 2016, AJ, 151, 46

Ammler-von Eiff, M., & Reiners, A. 2012, A&A, 542, A116

Astudillo-Defru, N., Delfosse, X., Bonfils, X., et al. 2017, A&A, 600, A13

Bailey, J., Butler, R. P., Tinney, C. G., et al. 2009, ApJ, 690, 743

Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373

Bernstein, R., Shectman, S. A., Gunnels, S. M., Mochnacki, S., & Athey, A. E. 2003, SPIE Conf. Ser., 4841, 1694

Beuzit, J. L., Vigan, A., Mouillet, D., et al. 2019, A&A, 631, A155

Biller, B. A., Grandjean, A., Messina, S., et al. 2022, A&A, 658, A145

Bonavita, M., Fontanive, C., Gratton, R., et al. 2022a, MNRAS, 513, 5588

Bonavita, M., Gratton, R., Desidera, S., et al. 2022b, A&A, 663, A144

Borgniet, S., Lagrange, A. M., Meunier, N., & Galland, F. 2017, A&A, 599, A57

Boro Saikia, S., Marvin, C. J., Jeffers, S. V., et al. 2018, A&A, 616, A108

Bouchy, F., Díaz, R. F., Hébrard, G., et al. 2013, A&A, 549, A49

Bowler, B. P., Cochran, W. D., Endl, M., et al. 2021, AJ, 161, 106

Brandt, T. D. 2021, ApJS, 254, 42

Brandt, T. D., Dupuy, T. J., & Bowler, B. P. 2019, AJ, 158, 140

Brandt, G. M., Brandt, T. D., Dupuy, T. J., Li, Y., & Michalik, D. 2021a, AJ, 161, 179

Brandt, G. M., Dupuy, T. J., Li, Y., et al. 2021b, AJ, 162, 301

Brandt, T. D., Dupuy, T. J., Li, Y., et al. 2021c, AJ, 162, 186

Brandt, G. M., Michalik, D., Brandt, T. D., et al. 2021d, AJ, 162, 230

Bryan, M. L., Knutson, H. A., Howard, A. W., et al. 2016, ApJ, 821, 89

Butler, R. P., Wright, J. T., Marcy, G. W., et al. 2006, ApJ, 646, 505

Butler, R. P., Vogt, S. S., Laughlin, G., et al. 2017, AJ, 153, 208

Costa Silva, A. R., Delgado Mena, E., & Tsantaki, M. 2020, A&A, 634, A136

Crane, J. D., Shectman, S. A., Butler, R. P., et al. 2010, SPIE Conf. Ser., 7735, 773553

Dalba, P. A., Kane, S. R., Howell, S. B., et al. 2021, AJ, 161, 123

Damasso, M., Sozzetti, A., Lovis, C., et al. 2020, A&A, 642, A31

Delgado Mena, E., Moya, A., Adibekyan, V., et al. 2019, A&A, 624, A78

Delorme, P., Meunier, N., Albert, D., et al. 2017, in SF2A-2017: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, eds. C. Reylé, P. Di Matteo, F. Herpin, et al.

Diego, F., Charalambous, A., Fish, A. C., & Walker, D. D. 1990, SPIE Conf. Ser., 1235, 562

dos Santos, L. A., Meléndez, J., Bedell, M., et al. 2017, MNRAS, 472, 3425

Eggenberger, A., Udry, S., Chauvin, G., et al. 2007, A&A, 474, 273

Ehrenreich, D., Lagrange, A. M., Montagnier, G., et al. 2010, A&A, 523, A73

Emsenhuber, A., Mordasini, C., Burn, R., et al. 2021, A&A, 656, A69

Enard, D. 1982, SPIE Conf. Ser., 331, 232

Feng, F., Butler, R. P., Jones, H. R. A., et al. 2021, MNRAS, 507, 2856

Feng, F., Butler, R. P., Vogt, S. S., et al. 2022, ApJS, 262, 21

Fischer, D. A., Marcy, G. W., & Spronck, J. F. P. 2014, ApJS, 210, 5

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASA, 25, 306

Franson, K., Bowler, B. P., Bonavita, M., et al. 2023, AJ, 165, 39

Gaia Collaboration 2020, VizieR Online Data Catalog: I/350

Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1

Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1

Galicher, R., Boccaletti, A., Mesa, D., et al. 2018, A&A, 615, A92

Galland, F., Lagrange, A. M., Udry, S., et al. 2005, A&A, 443, 337

Ghazaryan, S., Alecian, G., & Hakobyan, A. A. 2019, MNRAS, 487, 5922

Gomes da Silva, J., Santos, N. C., Adibekyan, V., et al. 2021, A&A, 646, A77

Gondoin, P. 2020, A&A, 641, A110

Gorrini, P., Astudillo-Defru, N., Dreizler, S., et al. 2022, A&A, 664, A64

Grandjean, A., Lagrange, A. M., Beust, H., et al. 2019, A&A, 629, C1

Grieves, N., Ge, J., Thomas, N., et al. 2018, MNRAS, 481, 3244

Hartmann, M., & Hatzes, A. P. 2015, A&A, 582, A84

Hébrard, G., Arnold, L., Forveille, T., et al. 2019, A&A, 631, A145

Hojjatpanah, S., Figueira, P., Santos, N. C., et al. 2019, A&A, 629, A80

Hojjatpanah, S., Oshagh, M., Figueira, P., et al. 2020, A&A, 639, A35

Jenkins, J. S., Jones, H. R. A., Tuomi, M., et al. 2017, MNRAS, 466, 443

Jones, H. R. A., Paul Butler, R., Marcy, G. W., et al. 2002, MNRAS, 337, 1170

Kane, S. R., Dalba, P. A., Li, Z., et al. 2019, AJ, 157, 252

Kervella, P., Arenou, F., Mignard, F., & Thévenin, F. 2019, A&A, 623, A72

Kervella, P., Arenou, F., & Thévenin, F. 2022, A&A, 657, A7

Kiefer, F., Hébrard, G., Sahlmann, J., et al. 2019, A&A, 631, A125

Kiefer, F., Hébrard, G., Lecavelier des Etangs, A., et al. 2021, A&A, 645, A7

Kuzuhara, M., Currie, T., Takarada, T., et al. 2022, ApJ, 934, L18

Lagrange, A. M., Meunier, N., Rubini, P., et al. 2019, Nat. Astron., 3, 1135

Lagrange, A. M., Philipot, F., Rubini, P., et al. 2023, A&A, 677, A71

Lenzen, R., Hartung, M., Brandner, W., et al. 2003, SPIE Conf. Ser., 4841, 944

Li, Y., Brandt, T. D., Brandt, G. M., et al. 2021, AJ, 162, 266

Liebing, F., Jeffers, S. V., Reiners, A., & Zechmeister, M. 2021, A&A, 654, A168

Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2

Llorente de Andrés, F., Chavero, C., de la Reza, R., Roca-Fàbrega, S., & Cifuentes, C. 2021, A&A, 654, A137

Lo Curto, G., Mayor, M., Benz, W., et al. 2010, A&A, 512, A48

Lo Curto, G., Pepe, F., Avila, G., et al. 2015, The Messenger, 162, 9

Maldonado, J., Micela, G., Baratella, M., et al. 2020, A&A, 644, A68

Mayor, M., & Queloz, D. 1995, Nature, 378, 355

Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20

Mayor, M., Udry, S., Naef, D., et al. 2004, A&A, 415, 391

Mayor, M., Marmier, M., Lovis, C., et al. 2011, arXiv e-prints [arXiv:1109.2497]

Melo, C., Santos, N. C., Gieren, W., et al. 2007, A&A, 467, 721

Meunier, N., & Lagrange, A. M. 2019, A&A, 628, A125

Meunier, N., Lagrange, A. M., Boulet, T., & Borgniet, S. 2019, A&A, 627, A56

Mordasini, C. 2018, in Handbook of Exoplanets, eds. H. J. Deeg, & J. A. Belmonte (Berlin: Springer), 143

Moutou, C., Mayor, M., Lo Curto, G., et al. 2009, A&A, 496, 513

Moutou, C., Mayor, M., Lo Curto, G., et al. 2011, A&A, 527, A63

Naef, D., Mayor, M., Pepe, F., et al. 2001, A&A, 375, 205

Naef, D., Mayor, M., Benz, W., et al. 2007, A&A, 470, 721

Naef, D., Mayor, M., Lo Curto, G., et al. 2010, A&A, 523, A15
Pavlov, A., Möller-Nilsson, O., Feldt, M., et al. 2008, SPIE Conf. Ser., 7019, 701939

Perruchot, S., Kohler, D., Bouchy, F., et al. 2008, SPIE Conf. Ser., 7014, 70140J

Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49

Philipot, F., Lagrange, A. M., Rubini, P., Kiefer, F., & Chomez, A. 2023, A&A, 670, A65

Queloz, D., Mayor, M., Weber, L., et al. 1999, arXiv e-prints [arXiv:astro-ph/9910223]

Reiners, A., & Zechmeister, M. 2020, ApJS, 247, 11

Rice, M., & Brewer, J. M. 2020, ApJ, 898, 119

Rigaut, F., Salmon, D., Arsenault, R., et al. 1998, PASP, 110, 152

Rosenthal, L. J., Fulton, B. J., Hirsch, L. A., et al. 2021, ApJS, 255, 8

Rousset, G., Lacombe, F., Puget, P., et al. 2003, SPIE Conf. Ser., 4839, 140

Sahlmann, J., Ségransan, D., Queloz, D., et al. 2011, A&A, 525, A95

Santos, N. C., Mayor, M., Bonfils, X., et al. 2011, A&A, 526, A112

Schneider, A., Melis, C., Song, I., & Zuckerman, B. 2011, ApJ, 743, 109

Scott, N. J., Howell, S. B., Horch, E. P., & Everett, M. E. 2018, PASP, 130, 054502

Ségransan, D., Udry, S., Mayor, M., et al. 2010, A&A, 511, A45

Ségransan, D., Mayor, M., Udry, S., et al. 2011, A&A, 535, A54

Service, M., Lu, J. R., Campbell, R., et al. 2016, PASP, 128, 095004

Sousa, S. G., Adibekyan, V., Delgado-Mena, E., et al. 2018, A&A, 620, A58

Stanford-Moore, S. A., Nielsen, E. L., De Rosa, R. J., Macintosh, B., & Czekala, I. 2020, ApJ, 898, 27

Stelzer, B., Marino, A., Micela, G., López-Santiago, J., & Liefke, C. 2013, MNRAS, 431, 2063

Swimmer, N., Currie, T., Steiger, S., et al. 2022, AJ, 164, 152

Tal-Or, L., Trifonov, T., Zucker, S., Mazeh, T., & Zechmeister, M. 2019, MNRAS, 484, L8

Tran, Q. H., Bowler, B. P., Cochran, W. D., et al. 2021, AJ, 161, 173

Trifonov, T., Tal-Or, L., Zechmeister, M., et al. 2020, A&A, 636, A74

Tull, R. G. 1998, SPIE Conf. Ser., 3355, 387

Udry, S., Mayor, M., Clausen, J. V., et al. 2003, A&A, 407, 679

van Leeuwen, F. 2007, A&A, 474, 653

Venner, A., Vanderburg, A., & Pearce, L. A. 2021, AJ, 162, 12

Vogt, S. S. 1987, PASP, 99, 1214

Vogt, S. S., Allen, S. L., Bigelow, B. C., et al. 1994, SPIE Conf. Ser., 2198, 362

Vogt, S. S., Butler, R. P., Marcy, G. W., et al. 2002, ApJ, 568, 352

Vogt, S. S., Radovan, M., Kibrick, R., et al. 2014, PASP, 126, 359

Ward-Duong, K., Patience, J., De Rosa, R. J., et al. 2015, MNRAS, 449, 2618

Wittenmyer, R. A., Endl, M., Cochran, W. D., Levison, H. F., & Henry, G. W. 2009, ApJS, 182, 97

Wittenmyer, R. A., Tuomi, M., Butler, R. P., et al. 2014, ApJ, 791, 114

Wittenmyer, R. A., Wang, S., Horner, J., et al. 2020, MNRAS, 492, 377

Yee, S. W., Petigura, E. A., & von Braun, K. 2017, ApJ, 836, 77

Zechmeister, M., Kürster, M., Endl, M., et al. 2013, A&A, 552, A78
## Appendix A: Target list

<table><caption>Table A.1: Proper motion values from HGCA of the 116 selected stars.</caption><thead><tr><th>Star</th><th>σ PMa<sup>a</sup></th><th>M<sub>b</sub> min at 10au<br>(M<sub>hip</sub>)</th><th>Parallax<sup>b</sup><br>(pc)</th><th>μ<sup>a</sup><sub>Hip</sub><br>(mas/yr)</th><th>μ<sup>b</sup><sub>Hip</sub><br>(mas/yr)</th><th>μ<sup>a</sup><sub>EDR3</sub><br>(mas/yr)</th><th>μ<sup>b</sup><sub>EDR3</sub><br>(mas/yr)</th><th>μ<sup>a</sup><sub>Hip-EDR3</sub><br>(mas/yr)</th><th>μ<sup>b</sup><sub>Hip-EDR3</sub><br>(mas/yr)</th><th>No. of<br>companion(s)</th></tr></thead><tbody><tr><td>GJ660.1</td><td>> 4</td><td>11.8</td><td>43.297 ± 0.397</td><td>192.695 ± 0.525</td><td>-704.419 ± 0.354</td><td>184.702 ± 4.626</td><td>-693.716 ± 3.139</td><td>190.922 ± 0.129</td><td>-702.635 ± 0.080</td><td>1</td></tr><tr><td>GJ676</td><td>> 4</td><td>3.4</td><td>62.609 ± 0.034</td><td>-258.759 ± 0.046</td><td>-185.119 ± 0.034</td><td>-260.087 ± 1.371</td><td>-185.157 ± 0.831</td><td>-258.838 ± 0.048</td><td>-185.880 ± 0.033</td><td>4</td></tr><tr><td>GJ680</td><td>> 4</td><td>16.4</td><td>103.345 ± 0.024</td><td>74.069 ± 0.033</td><td>470.166 ± 0.020</td><td>83.508 ± 3.523</td><td>454.994 ± 2.580</td><td>77.747 ± 0.082</td><td>464.975 ± 0.057</td><td>1</td></tr><tr><td>HD 142</td><td>> 4</td><td>11.7</td><td>38.213 ± 0.041</td><td>575.099 ± 0.024</td><td>-40.874 ± 0.033</td><td>575.297 ± 0.447</td><td>-39.223 ± 0.469</td><td>575.093 ± 0.012</td><td>-39.845 ± 0.013</td><td>2</td></tr><tr><td>HD 457</td><td>> 4</td><td>6.1</td><td>18.102 ± 0.025</td><td>113.112 ± 0.025</td><td>-18.873 ± 0.028</td><td>112.462 ± 0.592</td><td>-19.548 ± 0.678</td><td>112.791 ± 0.018</td><td>-18.826 ± 0.018</td><td>1*</td></tr><tr><td>HD 1388</td><td>> 4</td><td>28.5</td><td>37.176 ± 0.029</td><td>401.642 ± 0.043</td><td>-0.238 ± 0.024</td><td>396.154 ± 0.824</td><td>-0.481 ± 0.457</td><td>398.610 ± 0.021</td><td>-0.437 ± 0.016</td><td>1</td></tr><tr><td>HD7449</td><td>> 4</td><td>37.4</td><td>25.9132 ± 0.0287</td><td>-164.544 ± 0.053</td><td>-134.382 ± 0.039</td><td>-161.045 ± 0.621</td><td>-138.692 ± 0.464</td><td>-163.019 ± 0.022</td><td>-136.855 ± 0.017</td><td>2</td></tr><tr><td>HD11195</td><td>3.4</td><td>4.2</td><td>17.4998 ± 0.0184</td><td>-6.488 ± 0.028</td><td>14.994 ± 0.025</td><td>-6.833 ± 0.508</td><td>15.066 ± 0.502</td><td>-6.512 ± 0.018</td><td>14.763 ± 0.018</td><td>0</td></tr><tr><td>HD11397</td><td>> 4</td><td>11.9</td><td>18.4075 ± 0.0196</td><td>132.287 ± 0.027</td><td>-352.861 ± 0.020</td><td>126.848 ± 1.338</td><td>-363.072 ± 1.007</td><td>131.014 ± 0.045</td><td>-358.017 ± 0.033</td><td>1*</td></tr><tr><td>HD13350</td><td>2.0</td><td>0.1</td><td>54.6825 ± 0.0354</td><td>310.909 ± 0.064</td><td>1049.061 ± 0.082</td><td>311.667 ± 0.397</td><td>1050.236 ± 0.444</td><td>310.720 ± 0.013</td><td>1049.658 ± 0.012</td><td>0</td></tr><tr><td>HD13357</td><td>2.0</td><td>0.1</td><td>12.1502 ± 0.0157</td><td>-30.155 ± 0.022</td><td>-19.200 ± 0.020</td><td>-30.826 ± 1.055</td><td>-195.625 ± 0.919</td><td>-30.082 ± 0.029</td><td>-195.058 ± 0.031</td><td>0</td></tr><tr><td>HD13612</td><td>> 4</td><td>17.3</td><td>26.1678 ± 0.0274</td><td>366.258 ± 0.043</td><td>-67.275 ± 0.030</td><td>376.253 ± 9.264</td><td>-72.567 ± 6.630</td><td>367.067 ± 0.684</td><td>-66.783 ± 0.411</td><td>1*</td></tr><tr><td>HD15337</td><td>> 4</td><td>13.2</td><td>22.2922 ± 0.0160</td><td>-73.581 ± 0.021</td><td>-211.935 ± 0.020</td><td>-72.014 ± 1.023</td><td>-210.732 ± 0.920</td><td>-72.829 ± 0.034</td><td>-211.474 ± 0.031</td><td>1*</td></tr><tr><td>HD16160</td><td>> 4</td><td>75.3</td><td>35.8970 ± 0.1225</td><td>114.381 ± 0.209</td><td>-199.348 ± 0.131</td><td>99.098 ± 2.307</td><td>-187.790 ± 1.360</td><td>105.221 ± 0.079</td><td>-190.355 ± 0.046</td><td>2</td></tr><tr><td>HD16548</td><td>> 4</td><td>108.2</td><td>18.9061 ± 0.0226</td><td>261.625 ± 0.042</td><td>-62.568 ± 0.028</td><td>273.099 ± 0.730</td><td>-51.351 ± 0.795</td><td>266.536 ± 0.025</td><td>-61.001 ± 0.023</td><td>1</td></tr><tr><td>HD16905</td><td>> 4</td><td>5.1</td><td>25.1176 ± 0.0123</td><td>36.720 ± 0.024</td><td>-184.705 ± 0.018</td><td>37.464 ± 0.966</td><td>-183.811 ± 0.898</td><td>37.153 ± 0.029</td><td>-184.468 ± 0.029</td><td>1</td></tr><tr><td>HD19467</td><td>3.0</td><td>2.3</td><td>22.0344 ± 0.0177</td><td>36.061 ± 0.026</td><td>-80.180 ± 0.020</td><td>36.537 ± 1.713</td><td>-78.766 ± 1.067</td><td>36.267 ± 0.045</td><td>-80.069 ± 0.028</td><td>1*</td></tr><tr><td>HD19641</td><td>> 4</td><td>20</td><td>19.0259 ± 0.0216</td><td>-25.024 ± 0.024</td><td>114.795 ± 0.030</td><td>-25.973 ± 0.673</td><td>117.150 ± 0.749</td><td>-25.519 ± 0.020</td><td>115.886 ± 0.021</td><td>1*</td></tr><tr><td>HD20010</td><td>> 4</td><td>43</td><td>71.4337 ± 0.1320</td><td>359.565 ± 0.112</td><td>619.186 ± 0.161</td><td>371.142 ± 0.427</td><td>611.897 ± 0.550</td><td>364.729 ± 0.015</td><td>616.214 ± 0.017</td><td>1*</td></tr><tr><td>HD21019</td><td>2.8</td><td>0.1</td><td>26.7913 ± 0.0912</td><td>1.295 ± 0.246</td><td>-219.292 ± 0.170</td><td>2.035 ± 0.751</td><td>-219.220 ± 0.511</td><td>1.799 ± 0.022</td><td>-219.483 ± 0.014</td><td>1*</td></tr><tr><td>HD21175</td><td>> 4</td><td>59</td><td>56.5570 ± 0.0270</td><td>63.985 ± 0.033</td><td>35.175 ± 0.049</td><td>40.646 ± 0.780</td><td>42.888 ± 0.778</td><td>53.036 ± 0.026</td><td>39.153 ± 0.025</td><td>1</td></tr><tr><td>HD24633</td><td>> 4</td><td>32.4</td><td>18.7655 ± 0.0884</td><td>4.290 ± 0.147</td><td>64.156 ± 0.111</td><td>4.088 ± 1.494</td><td>64.266 ± 1.419</td><td>5.886 ± 0.042</td><td>65.529 ± 0.035</td><td>1*</td></tr><tr><td>HD27894</td><td>> 4</td><td>7</td><td>22.8888 ± 0.0121</td><td>182.473 ± 0.017</td><td>270.012 ± 0.023</td><td>182.047 ± 0.915</td><td>272.210 ± 1.075</td><td>182.032 ± 0.030</td><td>270.376 ± 0.033</td><td>0</td></tr><tr><td>HD28254</td><td>2.2</td><td>0.1</td><td>15.2528 ± 0.0183</td><td>181.376 ± 0.027</td><td>-164.130 ± 0.026</td><td>182.195 ± 1.058</td><td>-163.440 ± 0.488</td><td>181.367 ± 0.032</td><td>-164.082 ± 0.019</td><td>1</td></tr><tr><td>HD38459</td><td>> 4</td><td>2.9</td><td>26.5541 ± 0.0114</td><td>-19.556 ± 0.021</td><td>-39.926 ± 0.021</td><td>-20.643 ± 0.894</td><td>-40.680 ± 0.962</td><td>-19.706 ± 0.030</td><td>-40.138 ± 0.030</td><td>0</td></tr><tr><td>HD39091</td><td>2.7</td><td>2.3</td><td>18.2833 ± 0.0216</td><td>74.236 ± 0.026</td><td>-214.258 ± 0.023</td><td>75.863 ± 1.357</td><td>-213.220 ± 0.847</td><td>74.058 ± 0.042</td><td>-214.357 ± 0.024</td><td>0</td></tr><tr><td>HD40397</td><td>> 4</td><td>20</td><td>41.5366 ± 0.0210</td><td>74.849 ± 0.029</td><td>-200.354 ± 0.027</td><td>70.996 ± 0.789</td><td>-203.208 ± 0.541</td><td>73.753 ± 0.024</td><td>-202.096 ± 0.017</td><td>1</td></tr><tr><td>HD42581</td><td>> 4</td><td>7.7</td><td>138.3400 ± 0.3177</td><td>177.585 ± 0.589</td><td>1477.306 ± 0.340</td><td>1807.221 ± 1.454</td><td>1443.037 ± 0.713</td><td>1792.759 ± 0.033</td><td>1450.768 ± 0.022</td><td>0</td></tr><tr><td>HD42659</td><td>> 4</td><td>42</td><td>7.7597 ± 0.0308</td><td>-0.118 ± 0.033</td><td>-18.771 ± 0.036</td><td>0.465 ± 0.617</td><td>-17.998 ± 0.513</td><td>0.284 ± 0.017</td><td>-18.558 ± 0.016</td><td>1</td></tr><tr><td>HD42936</td><td>> 4</td><td>0.9</td><td>17.8589 ± 0.0131</td><td>323.791 ± 0.013</td><td>-281.239 ± 0.017</td><td>324.151 ± 0.515</td><td>-281.307 ± 0.860</td><td>323.830 ± 0.017</td><td>-280.993 ± 0.026</td><td>2</td></tr><tr><td>HD43197</td><td>3.7</td><td>5.7</td><td>16.0171 ± 0.0106</td><td>147.648 ± 0.013</td><td>16.600 ± 0.015</td><td>148.028 ± 1.432</td><td>13.923 ± 1.615</td><td>147.878 ± 0.044</td><td>16.381 ± 0.052</td><td>2</td></tr><tr><td>HD44219</td><td>2.2</td><td>0.1</td><td>31.2191 ± 0.0240</td><td>-8.694 ± 0.036</td><td>-260.642 ± 0.035</td><td>-7.903 ± 0.734</td><td>-260.163 ± 0.791</td><td>-8.445 ± 0.022</td><td>-260.600 ± 0.017</td><td>1</td></tr><tr><td>HD46569</td><td>> 4</td><td>87</td><td>27.1217 ± 0.0476</td><td>102.014 ± 0.074</td><td>104.577 ± 0.074</td><td>102.415 ± 0.418</td><td>104.925 ± 0.423</td><td>107.769 ± 0.014</td><td>103.583 ± 0.014</td><td>1</td></tr><tr><td>HD49095</td><td>> 4</td><td>74</td><td>41.7838 ± 0.0270</td><td>-217.138 ± 0.038</td><td>-317.283 ± 0.038</td><td>-204.905 ± 0.450</td><td>-304.276 ± 0.461</td><td>-212.357 ± 0.013</td><td>-310.692 ± 0.015</td><td>0</td></tr><tr><td>HD56380</td><td>> 4</td><td>162</td><td>18.7261 ± 0.2264</td><td>48.032 ± 0.304</td><td>113.332 ± 0.388</td><td>66.521 ± 2.491</td><td>105.182 ± 2.474</td><td>57.329 ± 0.080</td><td>104.799 ± 0.083</td><td>1</td></tr><tr><td>HD61383</td><td>> 4</td><td>136</td><td>18.2219 ± 0.0241</td><td>226.754 ± 0.034</td><td>-126.880 ± 0.026</td><td>228.806 ± 0.875</td><td>122.982 ± 0.509</td><td>230.681 ± 0.025</td><td>-122.299 ± 0.015</td><td>0</td></tr><tr><td>HD61986</td><td>> 4</td><td>5.7</td><td>32.9777 ± 0.0239</td><td>-148.238 ± 0.034</td><td>60.209 ± 0.028</td><td>-147.092 ± 0.945</td><td>59.709 ± 0.581</td><td>-148.613 ± 0.028</td><td>59.841 ± 0.014</td><td>0</td></tr><tr><td>HD62364</td><td>> 4</td><td>16</td><td>18.8811 ± 0.0229</td><td>-123.632 ± 0.045</td><td>-8.469 ± 0.039</td><td>-124.392 ± 0.536</td><td>-9.721 ± 0.541</td><td>-123.989 ± 0.018</td><td>-9.190 ± 0.017</td><td>1</td></tr><tr><td>HD63454</td><td>> 4</td><td>2.0</td><td>23.3044 ± 0.0283</td><td>114.008 ± 0.043</td><td>72.063 ± 0.037</td><td>112.003 ± 1.660</td><td>73.320 ± 1.260</td><td>114.225 ± 0.059</td><td>-72.085 ± 0.044</td><td>2*</td></tr><tr><td>HD65277</td><td>> 4</td><td>5.0</td><td>56.6412 ± 0.0254</td><td>-157.142 ± 0.032</td><td>7.580 ± 0.022</td><td>-155.772 ± 1.537</td><td>3.088 ± 1.089</td><td>-158.025 ± 0.043</td><td>6.923 ± 0.032</td><td>1*</td></tr><tr><td>HD67200</td><td>2.2</td><td>0.1</td><td>17.6470 ± 0.0197</td><td>-88.625 ± 0.028</td><td>-39.919 ± 0.024</td><td>-89.442 ± 0.818</td><td>-39.482 ± 0.872</td><td>-88.698 ± 0.028</td><td>-39.832 ± 0.025</td><td>1*</td></tr><tr><td>HD68161</td><td>> 4</td><td>152</td><td>3.7118 ± 0.0994</td><td>-9.131 ± 0.143</td><td>-1.733 ± 0.154</td><td>-9.992 ± 0.443</td><td>-1.443 ± 0.379</td><td>-9.854 ± 0.014</td><td>-1.243 ± 0.013</td><td>0</td></tr><tr><td>HD72659</td><td>> 4</td><td>11</td><td>19.2582 ± 0.0271</td><td>-112.313 ± 0.036</td><td>-96.386 ± 0.024</td><td>-112.566 ± 0.689</td><td>-98.187 ± 0.601</td><td>-112.625 ± 0.024</td><td>-97.048 ± 0.016</td><td>2</td></tr><tr><td>HD73256</td><td>> 4</td><td>14</td><td>27.2441 ± 0.0217</td><td>-182.193 ± 0.024</td><td>67.373 ± 0.029</td><td>-180.349 ± 0.430</td><td>66.238 ± 0.542</td><td>-181.090 ± 0.013</td><td>67.182 ± 0.019</td><td>2</td></tr><tr><td>HD74957</td><td>> 4</td><td>27</td><td>20.5730 ± 0.0153</td><td>64.486 ± 0.025</td><td>78.146 ± 0.022</td><td>59.966 ± 0.706</td><td>76.040 ± 0.669</td><td>62.862 ± 0.023</td><td>77.328 ± 0.022</td><td>1*</td></tr><tr><td>HD75302</td><td>> 4</td><td>5.7</td><td>34.6094 ± 0.0239</td><td>27.283 ± 0.029</td><td>112.918 ± 0.034</td><td>27.126 ± 0.665</td><td>112.245 ± 0.590</td><td>27.068 ± 0.021</td><td>112.434 ± 0.020</td><td>1</td></tr><tr><td>HD78612</td><td>> 4</td><td>101</td><td>24.4502 ± 0.0210</td><td>-368.061 ± 0.028</td><td>-108.564 ± 0.021</td><td>-355.978 ± 0.782</td><td>-108.083 ± 0.621</td><td>-362.803 ± 0.024</td><td>-109.072 ± 0.020</td><td>1</td></tr><tr><td>HD85390</td><td>2.4</td><td>0.1</td><td>65.5889 ± 0.0342</td><td>-415.589 ± 0.053</td><td>-213.991 ± 0.032</td><td>-415.899 ± 0.630</td><td>-214.685 ± 0.548</td><td>-415.721 ± 0.019</td><td>-213.898 ± 0.014</td><td>2</td></tr><tr><td>HD88072</td><td>2.4</td><td>2.2</td><td>29.9081 ± 0.0377</td><td>-262.971 ± 0.080</td><td>-3.649 ± 0.093</td><td>-263.271 ± 0.511</td><td>-3.995 ± 0.392</td><td>-263.333 ± 0.018</td><td>-3.496 ± 0.014</td><td>1*</td></tr><tr><td>HD88218</td><td>> 4</td><td>30</td><td>31.5137 ± 0.0307</td><td>-445.632 ± 0.036</td><td>18.774 ± 0.040</td><td>-443.157 ± 0.447</td><td>16.019 ± 0.521</td><td>-444.675 ± 0.014</td><td>16.690 ± 0.016</td><td>1*</td></tr><tr><td>HD89839</td><td>3.3</td><td>2.4</td><td>21.9392 ± 0.0275</td><td>-39.173 ± 0.037</td><td>-120.068 ± 0.031</td><td>-38.414 ± 1.088</td><td>-119.640 ± 0.458</td><td>-39.236 ± 0.033</td><td>-120.267 ± 0.018</td><td>1</td></tr><tr><td>HD93351</td><td>> 4</td><td>10</td><td>17.9129 ± 0.0168</td><td>88.176 ± 0.021</td><td>-142.442 ± 0.021</td><td>87.116 ± 1.324</td><td>-143.170 ± 0.942</td><td>87.655 ± 0.045</td><td>-142.596 ± 0.029</td><td>2</td></tr><tr><td>HD96116</td><td>> 4</td><td>128</td><td>17.5165 ± 0.0164</td><td>18.249 ± 0.022</td><td>-39.264 ± 0.019</td><td>30.843 ± 0.823</td><td>-35.226 ± 0.686</td><td>23.825 ± 0.026</td><td>-38.921 ± 0.025</td><td>1</td></tr><tr><td>HD97037</td><td>> 4</td><td>15</td><td>30.8043 ± 0.0212</td><td>-238.340 ± 0.039</td><td>-173.313 ± 0.025</td><td>-240.783 ± 0.558</td><td>-171.223 ± 0.559</td><td>-239.437 ± 0.020</td><td>-172.792 ± 0.017</td><td>1*</td></tr><tr><td>HD100777</td><td>> 4</td><td>7.8</td><td>20.1590 ± 0.0285</td><td>-11.535 ± 0.038</td><td>35.586 ± 0.031</td><td>-11.575 ± 0.990</td><td>34.643 ± 0.852</td><td>-11.268 ± 0.033</td><td>35.034 ± 0.030</td><td>1</td></tr><tr><td>HD101198</td><td>> 4</td><td>42</td><td>37.2450 ± 0.0730</td><td>91.849 ± 0.104</td><td>125.335 ± 0.069</td><td>99.625 ± 0.453</td><td>125.765 ± 0.337</td><td>95.408 ± 0.018</td><td>126.203 ± 0.012</td><td>1</td></tr><tr><td>HD103743</td><td>2.2</td><td>0.1</td><td>28.2336 ± 0.0231</td><td>-171.869 ± 0.038</td><td>-8.144 ± 0.026</td><td>-171.494 ± 7.626</td><td>-11.607 ± 4.054</td><td>-171.894 ± 0.382</td><td>-8.165 ± 0.241</td><td>0</td></tr><tr><td>HD104731</td><td>> 4</td><td>33</td><td>40.4345 ± 0.0884</td><td>318.554 ± 0.079</td><td>-111.120 ± 0.076</td><td>323.541 ± 0.384</td><td>-111.825 ± 0.371</td><td>321.144 ± 0.012</td><td>-111.265 ± 0.011</td><td>0</td></tr><tr><td>HD106906</td><td>3.5</td><td>9.2</td><td>9.7673 ± 0.0185</td><td>-39.066 ± 0.020</td><td>-12.698 ± 0.022</td><td>-38.695 ± 0.617</td><td>-12.303 ± 0.651</td><td>-39.222 ± 0.019</td><td>-12.684 ± 0.021</td><td>2</td></tr><tr><td>HD107094</td><td>> 4</td><td>19</td><td>18.6573 ± 0.0234</td><td>-3.315 ± 0.025</td><td>-45.172 ± 0.022</td><td>-8.639 ± 0.741</td><td>-44.534 ± 0.643</td><td>-2.702 ± 0.027</td><td>-46.030 ± 0.022</td><td>1*</td></tr><tr><td>HD108063</td><td>2.0</td><td>0.1</td><td>19.0734 ± 0.0231</td><td>-135.319 ± 0.027</td><td>-22.538 ± 0.018</td><td>-135.004 ± 0.451</td><td>-22.583 ± 0.432</td><td>-135.170 ± 0.014</td><td>-22.539 ± 0.014</td><td>1*</td></tr><tr><td>HD108202</td><td>2.3</td><td>0.4</td><td>25.8379 ± 0.0177</td><td>-309.967 ± 0.032</td><td>-77.177 ± 0.022</td><td>-311.265 ± 1.621</td><td>-77.031 ± 1.142</td><td>-309.702 ± 0.061</td><td>-77.415 ± 0.045</td><td>1*</td></tr><tr><td>HD110537</td><td>> 4</td><td>4.3</td><td>21.2496 ± 0.1143</td><td>-0.716 ± 0.182</td><td>-382.755 ± 0.176</td><td>0.561 ± 0.896</td><td>-383.939 ± 0.940</td><td>-1.046 ± 0.028</td><td>-382.604 ± 0.028</td><td>1*</td></tr><tr><td>HD111031</td><td>> 4</td><td>39</td><td>32.0310 ± 0.0219</td><td>-279.627 ± 0.034</td><td>46.791 ± 0.026</td><td>-282.073 ± 0.633</td><td>54.801 ± 0.454</td><td>-281.564 ± 0.022</td><td>49.571 ± 0.016</td><td>1</td></tr><tr><td>HD111232</td><td>> 4</td><td>5.5</td><td>24.4620 ± 0.0455</td><td>108.284 ± 0.077</td><td>-50.040 ± 0.081</td><td>108.246 ± 0.657</td><td>-49.514 ± 0.569</td><td>108.551 ± 0.025</td><td>-49.692 ± 0.021</td><td>2</td></tr><tr><td>HD111998</td><td>2.4</td><td>1.8</td><td>19.7872 ± 0.0205</td><td>-191.118 ± 0.019</td><td>-65.020 ± 0.021</td><td>-190.718 ± 0.513</td><td>-65.308 ± 0.649</td><td>-191.179 ± 0.017</td><td>-65.216 ± 0.021</td><td>1</td></tr><tr><td>HD114174</td><td>> 4</td><td>37</td><td>37.8677 ± 0.0243</td><td>85.533 ± 0.037</td><td>-680.324 ± 0.037</td><td>84.368 ± 0.759</td><td>-670.201 ± 0.603</td><td>84.993 ± 0.024</td><td>-676.775 ± 0.018</td><td>1*</td></tr><tr><td>HD114330</td><td>2.0</td><td>0.1</td><td>11.1834 ± 0.4110</td><td>-35.252 ± 0.752</td><td>-33.302 ± 0.552</td><td>-35.629 ± 0.990</td><td>-32.322 ± 0.829</td><td>-36.673 ± 0.049</td><td>-32.368 ± 0.038</td><td>0</td></tr><tr><td>HD114783</td><td>2.9</td><td>1.1</td><td>47.5529 ± 0.0291</td><td>-138.362 ± 0.047</td><td>10.284 ± 0.030</td><td>-138.541 ± 0.817</td><td>9.834 ± 0.851</td><td>-138.488 ± 0.028</td><td>10.508 ± 0.025</td><td>2</td></tr><tr><td>HD125276</td><td>> 4</td><td>11</td><td>55.5900 ± 0.0531</td><td>-359.037 ± 0.062</td><td>364.426 ± 0.056</td><td>-356.395 ± 0.887</td><td>366.796 ± 0.737</td><td>-357.805 ± 0.026</td><td>365.366 ± 0.019</td><td>0</td></tr><tr><td>HD128571</td><td>> 4</td><td>6.1</td><td>17.0083 ± 0.0184</td><td>-71.750 ± 0.017</td><td>19.160 ± 0.025</td><td>-71.718 ± 0.503</td><td>19.211 ± 0.687</td><td>-71.476 ± 0.016</td><td>18.924 ± 0.021</td><td>1*</td></tr></tbody></table>

Notes: $^a$: Significance of the PMa. $^b$: Gaia DR3. $\mu_{Hip}$ corresponds to the proper motion obtained by HIPPARCOS. $^*$: Orbital parameters and mass of the companion poorly constrained. $\mu_{EDR3}$ corresponds to the proper motion obtained by Gaia EDR3. $\mu_{Hip-EDR3}$ corresponds to the proper motion obtained by the HIPPARCOS-Gaia EDR3 positional difference.
<table><caption>Table A.2: Table A.1 continued.</caption><thead><tr><th>Star</th><th>σ PMa<sup>a</sup></th><th>M<sub>b</sub> at 10au<br>(M<sub>10au</sub>)</th><th>Parallax<sup>b</sup><br>(mas)</th><th>μ<sup>2</sup><sub>HIP</sub><br>(mas/yr)</th><th>μ<sup>2</sup><sub>HIP</sub><br>(mas/yr)</th><th>μ<sup>2</sup><sub>EDR3</sub><br>(mas/yr)</th><th>μ<sup>2</sup><sub>EDR3</sub><br>(mas/yr)</th><th>μ<sup>a</sup><sub>HIP-EDR3</sub><br>(mas/yr)</th><th>μ<sup>b</sup><sub>HIP-EDR3</sub><br>(mas/yr)</th><th>No. of<br>companion(s)</th></tr></thead><tbody><tr><td>HD129191</td><td>&gt; 4</td><td>11</td><td>18.4622 ± 0.0242</td><td>26.520 ± 0.042</td><td>-152.832 ± 0.035</td><td>26.278 ± 1.394</td><td>-152.293 ± 1.540</td><td>26.058 ± 0.039</td><td>-152.433 ± 0.038</td><td>1*</td></tr><tr><td>HD131664</td><td>&gt; 4</td><td>99</td><td>19.1360 ± 0.0829</td><td>8.008 ± 0.105</td><td>22.284 ± 0.126</td><td>14.979 ± 0.710</td><td>28.685 ± 0.740</td><td>11.370 ± 0.022</td><td>25.899 ± 0.023</td><td>1</td></tr><tr><td>HD131977</td><td>&gt; 4</td><td>3.7</td><td>169.8843 ± 0.0653</td><td>1031.472 ± 0.093</td><td>-1723.619 ± 0.075</td><td>1035.702 ± 1.343</td><td>-1725.446 ± 1.007</td><td>1032.968 ± 0.051</td><td>-1724.781 ± 0.033</td><td>0</td></tr><tr><td>HD135625</td><td>2.0</td><td>0.1</td><td>18.9018 ± 0.0191</td><td>72.385 ± 0.026</td><td>-16.574 ± 0.027</td><td>72.305 ± 0.974</td><td>-17.487 ± 0.728</td><td>72.283 ± 0.023</td><td>-16.433 ± 0.016</td><td>1</td></tr><tr><td>HD140901</td><td>2.0</td><td>0.1</td><td>26.2048 ± 0.0179</td><td>-4.040 ± 0.019</td><td>-176.046 ± 0.025</td><td>-3.664 ± 0.508</td><td>-176.463 ± 0.808</td><td>-4.053 ± 0.015</td><td>-176.132 ± 0.022</td><td>1</td></tr><tr><td>HD142070</td><td>&gt; 4</td><td>81</td><td>5.5771 ± 0.0307</td><td>-14.315 ± 0.043</td><td>-24.244 ± 0.033</td><td>-15.266 ± 1.073</td><td>-25.030 ± 1.016</td><td>-13.895 ± 0.031</td><td>-23.580 ± 0.023</td><td>0</td></tr><tr><td>HD143361</td><td>3.0</td><td>3.0</td><td>14.5456 ± 0.0202</td><td>-156.561 ± 0.026</td><td>-120.231 ± 0.019</td><td>-157.333 ± 1.362</td><td>-119.387 ± 1.039</td><td>-156.630 ± 0.040</td><td>-120.369 ± 0.027</td><td>2</td></tr><tr><td>HD150139</td><td>2.3</td><td>0.1</td><td>20.5407 ± 0.0251</td><td>139.201 ± 0.044</td><td>-143.615 ± 0.037</td><td>138.903 ± 0.807</td><td>-142.494 ± 0.594</td><td>139.101 ± 0.031</td><td>-143.648 ± 0.023</td><td>1*</td></tr><tr><td>HD156079</td><td>&gt; 4</td><td>4.4</td><td>17.5648 ± 0.0239</td><td>-69.179 ± 0.031</td><td>-165.076 ± 0.023</td><td>-69.736 ± 0.909</td><td>-164.554 ± 0.561</td><td>-68.923 ± 0.028</td><td>-165.238 ± 0.015</td><td>0</td></tr><tr><td>HD156274</td><td>&gt; 4</td><td>9.9</td><td>113.7513 ± 0.0725</td><td>1029.609 ± 0.117</td><td>106.935 ± 0.101</td><td>1036.484 ± 1.039</td><td>109.036 ± 0.520</td><td>1033.069 ± 0.039</td><td>107.932 ± 0.021</td><td>1*</td></tr><tr><td>HD165131</td><td>&gt; 4</td><td>15</td><td>17.2691 ± 0.0317</td><td>-61.918 ± 0.043</td><td>-124.587 ± 0.030</td><td>-59.570 ± 1.126</td><td>-125.994 ± 0.717</td><td>-62.476 ± 0.036</td><td>-124.074 ± 0.023</td><td>1</td></tr><tr><td>HD166724</td><td>3.2</td><td>1.6</td><td>37.4641 ± 0.0253</td><td>-43.140 ± 0.040</td><td>72.728 ± 0.037</td><td>-43.593 ± 0.393</td><td>73.244 ± 0.361</td><td>-43.414 ± 0.013</td><td>72.873 ± 0.012</td><td>1*</td></tr><tr><td>HD167677</td><td>3.2</td><td>3.0</td><td>16.7894 ± 0.0166</td><td>33.127 ± 0.020</td><td>12.219 ± 0.017</td><td>33.530 ± 1.050</td><td>12.868 ± 0.828</td><td>33.235 ± 0.035</td><td>12.392 ± 0.027</td><td>1</td></tr><tr><td>HD169830</td><td>2.0</td><td>0.1</td><td>27.1461 ± 0.1469</td><td>-0.341 ± 0.221</td><td>16.103 ± 0.158</td><td>-0.721 ± 0.881</td><td>15.434 ± 0.537</td><td>0.064 ± 0.023</td><td>16.374 ± 0.014</td><td>2</td></tr><tr><td>HD171587</td><td>2.0</td><td>0.1</td><td>28.5668 ± 0.0217</td><td>126.024 ± 0.025</td><td>-193.681 ± 0.025</td><td>126.315 ± 0.938</td><td>-194.794 ± 0.733</td><td>125.979 ± 0.034</td><td>-193.715 ± 0.023</td><td>0</td></tr><tr><td>HD174429</td><td>3.1</td><td>2.1</td><td>21.1621 ± 0.0223</td><td>16.272 ± 0.025</td><td>-85.519 ± 0.024</td><td>17.353 ± 1.217</td><td>-83.766 ± 0.830</td><td>16.081 ± 0.035</td><td>-85.590 ± 0.025</td><td>0</td></tr><tr><td>HD179346</td><td>&gt; 4</td><td>53</td><td>18.0318 ± 0.0942</td><td>-7.800 ± 0.126</td><td>-53.982 ± 0.116</td><td>-9.922 ± 0.877</td><td>-52.875 ± 0.720</td><td>-4.943 ± 0.030</td><td>-53.287 ± 0.021</td><td>0</td></tr><tr><td>HD181433</td><td>&gt; 4</td><td>2.1</td><td>37.0511 ± 0.0211</td><td>-230.723 ± 0.018</td><td>235.806 ± 0.025</td><td>-230.886 ± 0.841</td><td>235.382 ± 0.996</td><td>-230.958 ± 0.030</td><td>235.640 ± 0.029</td><td>3</td></tr><tr><td>HD183414</td><td>2.2</td><td>0.2</td><td>36.1158 ± 0.0264</td><td>1.440 ± 0.037</td><td>-174.330 ± 0.042</td><td>2.701 ± 0.701</td><td>-174.656 ± 0.464</td><td>1.268 ± 0.023</td><td>-174.389 ± 0.015</td><td>0</td></tr><tr><td>HD185283</td><td>2.3</td><td>0.8</td><td>31.9046 ± 0.0217</td><td>-163.887 ± 0.022</td><td>-53.650 ± 0.017</td><td>-162.481 ± 1.218</td><td>-53.718 ± 1.109</td><td>-163.982 ± 0.033</td><td>-53.707 ± 0.025</td><td>1</td></tr><tr><td>HD191797</td><td>&gt; 4</td><td>58</td><td>22.7336 ± 0.0272</td><td>-8.443 ± 0.038</td><td>51.173 ± 0.033</td><td>-12.621 ± 2.324</td><td>54.340 ± 1.456</td><td>-5.479 ± 0.051</td><td>48.047 ± 0.033</td><td>1</td></tr><tr><td>HD195145</td><td>&gt; 4</td><td>3.1</td><td>16.6547 ± 0.0305</td><td>4.132 ± 0.052</td><td>-4.200 ± 0.050</td><td>1.989 ± 0.769</td><td>-4.698 ± 0.481</td><td>4.050 ± 0.021</td><td>-4.430 ± 0.012</td><td>0</td></tr><tr><td>HD195564</td><td>&gt; 4</td><td>21</td><td>40.4241 ± 0.0497</td><td>309.736 ± 0.068</td><td>109.910 ± 0.051</td><td>307.694 ± 0.651</td><td>107.216 ± 0.518</td><td>308.366 ± 0.021</td><td>108.386 ± 0.012</td><td>1</td></tr><tr><td>HD196050</td><td>&gt; 4</td><td>6.1</td><td>22.2309 ± 0.0234</td><td>-219.316 ± 0.041</td><td>-179.612 ± 0.039</td><td>-219.148 ± 0.675</td><td>-180.370 ± 0.547</td><td>-219.413 ± 0.026</td><td>-179.984 ± 0.023</td><td>1</td></tr><tr><td>HD196385</td><td>&gt; 4</td><td>24</td><td>20.9103 ± 0.0239</td><td>73.036 ± 0.031</td><td>-16.412 ± 0.024</td><td>73.293 ± 0.772</td><td>-15.056 ± 0.445</td><td>73.633 ± 0.022</td><td>-15.383 ± 0.012</td><td>1*</td></tr><tr><td>HD202206</td><td>3.2</td><td>2.2</td><td>22.2420 ± 0.0740</td><td>352.541 ± 0.107</td><td>-31.590 ± 0.064</td><td>350.543 ± 0.977</td><td>-32.314 ± 0.509</td><td>352.754 ± 0.031</td><td>-31.388 ± 0.022</td><td>2</td></tr><tr><td>HD204313</td><td>&gt; 4</td><td>11</td><td>20.7705 ± 0.0343</td><td>42.751 ± 0.054</td><td>-270.443 ± 0.034</td><td>42.551 ± 0.747</td><td>-271.125 ± 0.483</td><td>42.234 ± 0.027</td><td>-270.727 ± 0.018</td><td>3</td></tr><tr><td>HD204961</td><td>&gt; 4</td><td>0.7</td><td>201.3252 ± 0.0237</td><td>-45.917 ± 0.032</td><td>-816.875 ± 0.025</td><td>-45.677 ± 1.028</td><td>-817.318 ± 0.656</td><td>-46.049 ± 0.030</td><td>-816.302 ± 0.025</td><td>1</td></tr><tr><td>HD207700</td><td>&gt; 4</td><td>11</td><td>26.0873 ± 0.0172</td><td>-196.563 ± 0.025</td><td>-437.609 ± 0.030</td><td>-195.605 ± 0.674</td><td>-438.053 ± 0.573</td><td>-195.753 ± 0.019</td><td>-437.915 ± 0.018</td><td>1</td></tr><tr><td>HD207869</td><td>3.1</td><td>1.6</td><td>19.5789 ± 0.0118</td><td>-156.358 ± 0.014</td><td>-248.956 ± 0.017</td><td>-157.311 ± 0.951</td><td>-249.431 ± 0.807</td><td>-156.228 ± 0.032</td><td>-248.970 ± 0.030</td><td>0</td></tr><tr><td>HD210797</td><td>&gt; 4</td><td>13</td><td>15.4128 ± 0.0210</td><td>258.946 ± 0.027</td><td>-122.008 ± 0.024</td><td>258.354 ± 1.303</td><td>-118.366 ± 1.076</td><td>258.487 ± 0.037</td><td>-121.500 ± 0.036</td><td>1*</td></tr><tr><td>HD212036</td><td>&gt; 4</td><td>69</td><td>19.1167 ± 0.0803</td><td>21.555 ± 0.105</td><td>-234.978 ± 0.116</td><td>29.075 ± 1.182</td><td>-224.715 ± 1.505</td><td>22.904 ± 0.036</td><td>-231.443 ± 0.033</td><td>1</td></tr><tr><td>HD215257</td><td>&gt; 4</td><td>6.0</td><td>23.6251 ± 0.0251</td><td>150.456 ± 0.035</td><td>333.617 ± 0.033</td><td>150.934 ± 0.889</td><td>331.705 ± 0.653</td><td>150.716 ± 0.026</td><td>333.292 ± 0.020</td><td>1*</td></tr><tr><td>HD215497</td><td>2.0</td><td>0.1</td><td>24.6569 ± 0.0152</td><td>-54.666 ± 0.016</td><td>-61.147 ± 0.019</td><td>-53.881 ± 0.939</td><td>-60.575 ± 0.933</td><td>-54.665 ± 0.031</td><td>-61.222 ± 0.033</td><td>2</td></tr><tr><td>HD216437</td><td>2.8</td><td>1.7</td><td>26.9269 ± 0.0154</td><td>52.792 ± 0.025</td><td>215.141 ± 0.027</td><td>51.858 ± 0.767</td><td>214.317 ± 0.865</td><td>52.863 ± 0.026</td><td>214.857 ± 0.028</td><td>1</td></tr><tr><td>HD218396</td><td>2.1</td><td>3.5</td><td>28.2075 ± 0.0230</td><td>-179.085 ± 0.037</td><td>-6.549 ± 0.023</td><td>-196.276 ± 10.121</td><td>-9.146 ± 5.480</td><td>-179.244 ± 0.344</td><td>-6.934 ± 0.226</td><td>0</td></tr><tr><td>HD219834</td><td>&gt; 4</td><td>0.1</td><td>173.5740 ± 0.0170</td><td>-135.692 ± 0.015</td><td>-719.178 ± 0.023</td><td>-137.136 ± 0.575</td><td>-713.779 ± 0.834</td><td>-136.426 ± 0.018</td><td>-715.027 ± 0.028</td><td>0</td></tr><tr><td>HD221146</td><td>&gt; 4</td><td>29</td><td>27.5348 ± 0.0247</td><td>-10.513 ± 0.033</td><td>-24.287 ± 0.024</td><td>-11.763 ± 1.093</td><td>-28.467 ± 0.843</td><td>-11.268 ± 0.030</td><td>-25.794 ± 0.022</td><td>1</td></tr><tr><td>HD221420</td><td>&gt; 4</td><td>16.3</td><td>32.1023 ± 0.0325</td><td>16.306 ± 0.056</td><td>0.736 ± 0.058</td><td>15.885 ± 0.437</td><td>1.647 ± 0.385</td><td>14.973 ± 0.014</td><td>0.455 ± 0.012</td><td>1</td></tr><tr><td>HD221638</td><td>&gt; 4</td><td>40</td><td>19.1599 ± 0.0269</td><td>56.624 ± 0.041</td><td>72.839 ± 0.042</td><td>52.210 ± 0.623</td><td>77.925 ± 0.609</td><td>55.789 ± 0.019</td><td>74.527 ± 0.019</td><td>1</td></tr><tr><td>HD222237</td><td>&gt; 4</td><td>3.7</td><td>87.3724 ± 0.0187</td><td>143.736 ± 0.030</td><td>-736.907 ± 0.032</td><td>141.756 ± 0.524</td><td>-736.300 ± 0.507</td><td>142.849 ± 0.018</td><td>-736.651 ± 0.017</td><td>1*</td></tr><tr><td>HD223238</td><td>2.5</td><td>0.1</td><td>17.8497 ± 0.0171</td><td>-1.916 ± 0.016</td><td>-202.994 ± 0.018</td><td>-2.808 ± 1.100</td><td>-202.424 ± 1.016</td><td>-1.814 ± 0.037</td><td>-203.140 ± 0.035</td><td>0</td></tr><tr><td>HIP9603</td><td>&gt; 4</td><td>109</td><td>18.8225 ± 0.0866</td><td>-111.175 ± 0.143</td><td>257.497 ± 0.147</td><td>-110.184 ± 0.933</td><td>245.727 ± 0.990</td><td>-119.933 ± 0.033</td><td>255.468 ± 0.033</td><td>1</td></tr><tr><td>HIP10337</td><td>&gt; 4</td><td>1.9</td><td>42.1225 ± 0.0167</td><td>375.454 ± 0.019</td><td>55.718 ± 0.020</td><td>377.499 ± 1.622</td><td>55.290 ± 1.164</td><td>375.597 ± 0.051</td><td>55.414 ± 0.040</td><td>1</td></tr><tr><td>HIP39470</td><td>&gt; 4</td><td>91</td><td>22.6284 ± 0.0325</td><td>16.231 ± 0.044</td><td>-82.449 ± 0.047</td><td>9.228 ± 1.155</td><td>-82.499 ± 1.131</td><td>18.669 ± 0.039</td><td>-76.357 ± 0.041</td><td>1</td></tr><tr><td>HIP54597</td><td>2.3</td><td>0.1</td><td>19.3611 ± 0.0183</td><td>40.676 ± 0.019</td><td>53.190 ± 0.025</td><td>41.256 ± 0.462</td><td>54.595 ± 0.690</td><td>40.793 ± 0.016</td><td>53.317 ± 0.022</td><td>1</td></tr><tr><td>HIP70849</td><td>&gt; 4</td><td>4.6</td><td>41.4618 ± 0.0175</td><td>-44.051 ± 0.023</td><td>-201.577 ± 0.028</td><td>-47.003 ± 2.099</td><td>-203.284 ± 1.890</td><td>-44.428 ± 0.060</td><td>-202.047 ± 0.044</td><td>1</td></tr><tr><td>HIP113201</td><td>&gt; 4</td><td>132</td><td>68.7356 ± 0.0272</td><td>341.147 ± 0.037</td><td>-219.108 ± 0.034</td><td>343.238 ± 2.100</td><td>-218.067 ± 3.265</td><td>341.100 ± 0.055</td><td>-219.083 ± 0.063</td><td>1</td></tr></tbody></table>
<table><caption>Table A.3: Characteristics of the 53 stars with characterized companions.</caption><thead><tr><th>Name</th><th>HIP</th><th>Mass<sup>1</sup><br>(M<sub>☉</sub>)</th><th>Distance<sup>1</sup><br>(pc)</th><th>Spectral<sup>2</sup><br>type</th><th>V<sup>2</sup><br>(mag)</th><th>logR'<sub>HK</sub></th><th>Vsin(I)<br>(km s<sup>-1</sup>)</th><th>Age<br>(Gyr)</th><th>χ<sup>2</sup> 22</th></tr></thead><tbody><tr><td>GJ660.1</td><td>84212</td><td>0.43 ± 0.01</td><td>23.10 ± 0.21</td><td>M1V</td><td>11.62</td><td>-5.30<sup>7</sup></td><td>-</td><td>> 2<sup>23</sup></td><td>69</td></tr><tr><td>GJ680</td><td>86057</td><td>0.47 ± 0.01</td><td>9.676 ± 0.003</td><td>M3V</td><td>10.13</td><td>-5.16<sup>7</sup></td><td>3.0<sup>12</sup></td><td>8.1 ± 4.2<sup>24</sup></td><td>8601</td></tr><tr><td>HD1388</td><td>1444</td><td>1.18 ± 0.06</td><td>26.90 ± 0.02</td><td>G0V</td><td>6.50</td><td>-4.98<sup>6</sup></td><td>3.1<sup>13</sup></td><td>4.6 ± 1.7<sup>25</sup></td><td>4141</td></tr><tr><td>HD16160</td><td>12114</td><td>0.80 ± 0.04</td><td>7.23 ± 0.02</td><td>K3V</td><td>5.60<sup>3</sup></td><td>-4.87<sup>3</sup></td><td>4.7<sup>13</sup></td><td>4.7 ± 3.9<sup>13</sup></td><td>6078</td></tr><tr><td>HD16548</td><td>12350</td><td>1.31 ± 0.07</td><td>52.81 ± 0.07</td><td>G5V</td><td>6.98</td><td>-5.13<sup>3</sup></td><td>5.2<sup>13</sup></td><td>4.9 ± 0.3<sup>25</sup></td><td>11460</td></tr><tr><td>HD16905</td><td>12436</td><td>0.85 ± 0.04</td><td>39.80 ± 0.02</td><td>K3V</td><td>9.46</td><td>-4.72<sup>3</sup></td><td>1.5<sup>14</sup></td><td>6.1 ± 4.7<sup>25</sup></td><td>184</td></tr><tr><td>HD21175</td><td>15799</td><td>0.86 ± 0.04</td><td>17.68 ± 0.01</td><td>K0V</td><td>6.90</td><td>-4.74<sup>6</sup></td><td>1.8<sup>14</sup></td><td>3.7<sup>+3.3 6</sup><sub>-1.7</sub></td><td>69960</td></tr><tr><td>HD28254</td><td>20606</td><td>1.10 ± 0.05</td><td>55.16 ± 0.08</td><td>G5V</td><td>7.68</td><td>-5.17<sup>3</sup></td><td>2.3<sup>5</sup></td><td>6.0 ± 1.3<sup>25</sup></td><td>10.3</td></tr><tr><td>HD42659</td><td>29365</td><td>2.20 ± 0.11</td><td>128.5 ± 0.6</td><td>A3</td><td>6.75</td><td>-</td><td>19.0<sup>15</sup></td><td>-</td><td>15</td></tr><tr><td>HD44219</td><td>30114</td><td>1.20 ± 0.06</td><td>52.84 ± 0.06</td><td>G3V</td><td>7.70</td><td>-5.04<sup>3</sup></td><td>2.2<sup>13</sup></td><td>8.4 ± 0.9<sup>25</sup></td><td>23</td></tr><tr><td>HD46569</td><td>31079</td><td>1.50 ± 0.07</td><td>36.87 ± 0.07</td><td>F8V</td><td>5.57</td><td>-4.93<sup>6</sup></td><td>5.9<sup>16</sup></td><td>7.6<sup>+3.26 6</sup><sub>-3.0</sub></td><td>6232</td></tr><tr><td>HD56380</td><td>34961</td><td>0.98 ± 0.05</td><td>53.1 ± 0.3</td><td>G8V</td><td>9.20</td><td>-5.00<sup>3</sup></td><td>0.5<sup>14</sup></td><td>12.2 ± 1.1<sup>25</sup></td><td>3697</td></tr><tr><td>HD61383</td><td>37233</td><td>1.17 ± 0.06</td><td>53.05 ± 0.08</td><td>G3V</td><td>7.57</td><td>-5.02<sup>3</sup></td><td>2.7<sup>13</sup></td><td>11.4 ± 0.2<sup>13</sup></td><td>36110</td></tr><tr><td>HD62364</td><td>36941</td><td>1.30 ± 0.06</td><td>52.95 ± 0.08</td><td>F7V</td><td>7.31</td><td>-4.95<sup>3</sup></td><td>3.6<sup>13</sup></td><td>3.9 ± 0.7<sup>25</sup></td><td>364</td></tr><tr><td>HD73256</td><td>42214</td><td>1.00 ± 0.05</td><td>36.68 ± 0.04</td><td>G8IV</td><td>8.06</td><td>-4.46<sup>5</sup></td><td>3.50<sup>13</sup></td><td>2.2 ± 2.1<sup>25</sup></td><td>2040</td></tr><tr><td>HD75302</td><td>43297</td><td>1.06 ± 0.05</td><td>30.29 ± 0.03</td><td>G5V</td><td>7.46<sup>4</sup></td><td>-4.62<sup>6</sup></td><td>0.46<sup>17</sup></td><td>2.1<sup>+2.6 6</sup><sub>-1.1</sub></td><td>299</td></tr><tr><td>HD78612</td><td>44896</td><td>1.20 ± 0.06</td><td>40.85 ± 0.04</td><td>G3V</td><td>7.15</td><td>-4.98<sup>3</sup></td><td>1.6<sup>13</sup></td><td>10.2 ± 0.8<sup>25</sup></td><td>20060</td></tr><tr><td>HD89839</td><td>50653</td><td>1.30 ± 0.06</td><td>50.10 ± 0.07</td><td>F7V</td><td>7.64</td><td>-4.94<sup>3</sup></td><td>4.2<sup>14</sup></td><td>1.6 ± 0.9<sup>25</sup></td><td>56</td></tr><tr><td>HD93351</td><td>52720</td><td>0.98 ± 0.05</td><td>55.72 ± 0.05</td><td>G6V</td><td>9.14</td><td>-4.98<sup>3</sup></td><td>-</td><td>4.3 ± 3.9<sup>13</sup></td><td>143</td></tr><tr><td>HD96116</td><td>54102</td><td>1.08 ± 0.05</td><td>57.05 ± 0.06</td><td>G3V</td><td>8.69</td><td>-4.69<sup>3</sup></td><td>1.7<sup>13</sup></td><td>2.1 ± 1.8<sup>25</sup></td><td>28020</td></tr><tr><td>HD100777</td><td>56572</td><td>1.08 ± 0.05</td><td>49.52 ± 0.08</td><td>G8V</td><td>8.42</td><td>-5.09<sup>3</sup></td><td>1.8<sup>13</sup></td><td>5.7 ± 3.3<sup>25</sup></td><td>173</td></tr><tr><td>HD101198</td><td>56802</td><td>1.33 ± 0.07</td><td>26.82 ± 0.03</td><td>F6.5V</td><td>5.48</td><td>-4.87<sup>6</sup></td><td>5<sup>18</sup></td><td>6.2<sup>+3.26 6</sup><sub>-2.6</sub></td><td>1529</td></tr><tr><td>HD107094</td><td>60051</td><td>0.98 ± 0.05</td><td>53.54 ± 0.07</td><td>G5V</td><td>9.13</td><td>-4.83<sup>3</sup></td><td>1.8<sup>13</sup></td><td>8.6 ± 4.0<sup>25</sup></td><td>1225</td></tr><tr><td>HD108202</td><td>60648</td><td>0.74 ± 0.04</td><td>38.65 ± 0.03</td><td>K4/5V</td><td>10.22</td><td>-4.50<sup>3</sup></td><td>-</td><td>6.5 ± 4.0<sup>28</sup></td><td>28</td></tr><tr><td>HD111031</td><td>62345</td><td>1.17 ± 0.06</td><td>31.19 ± 0.02</td><td>G5V</td><td>6.87</td><td>-5.09<sup>3</sup></td><td>1.9<sup>13</sup></td><td>4.1 ± 1.6<sup>25</sup></td><td>8191</td></tr><tr><td>HD111998</td><td>62875</td><td>1.40 ± 0.07</td><td>33.39 ± 0.05</td><td>F6V</td><td>6.11</td><td>-4.58<sup>6</sup></td><td>20.0<sup>18</sup></td><td>1.7<sup>+3.4 6</sup><sub>-0.9</sub></td><td>27</td></tr><tr><td>HD114783</td><td>64457</td><td>0.90 ± 0.04</td><td>21.01 ± 0.02</td><td>K1V</td><td>7.55</td><td>-4.97<sup>3</sup></td><td>2.0<sup>13</sup></td><td>6.0 ± 4.9<sup>25</sup></td><td>33</td></tr><tr><td>HD131664</td><td>73408</td><td>1.16 ± 0.06</td><td>52.2 ± 0.3</td><td>G3V</td><td>8.13</td><td>-4.85<sup>5</sup></td><td>3.0<sup>13</sup></td><td>1.2 ± 1.0<sup>25</sup></td><td>1823</td></tr><tr><td>HD135625</td><td>74856</td><td>1.28 ± 0.06</td><td>56.58 ± 0.07</td><td>G3IV/V</td><td>7.71</td><td>-4.99<sup>6</sup></td><td>4.5<sup>13</sup></td><td>4.4 ± 0.9<sup>25</sup></td><td>7</td></tr><tr><td>HD140901</td><td>77358</td><td>1.06 ± 0.05</td><td>15.24 ± 0.01</td><td>G7IV</td><td>6.01</td><td>-4.79<sup>6</sup></td><td>2.0<sup>13</sup></td><td>3.2 ± 2.8<sup>25</sup></td><td>18</td></tr><tr><td>HD143361</td><td>78521</td><td>1.06 ± 0.05</td><td>68.62 ± 0.10</td><td>G6V</td><td>9.16</td><td>-5.01<sup>5</sup></td><td>0.68<sup>5</sup></td><td>6.3 ± 3.7<sup>29</sup></td><td>21</td></tr><tr><td>HD165131</td><td>88595</td><td>1.10 ± 0.06</td><td>57.8 ± 0.1</td><td>G3/5V</td><td>8.41</td><td>-4.98<sup>3</sup></td><td>2.3<sup>19</sup></td><td>2.7 ± 2.0<sup>25</sup></td><td>264</td></tr><tr><td>HD167677</td><td>89583</td><td>1.00 ± 0.05</td><td>54.59 ± 0.07</td><td>G5V</td><td>8.92</td><td>-4.98<sup>3</sup></td><td>0.7<sup>19</sup></td><td>11.5 ± 1.7<sup>25</sup></td><td>22</td></tr><tr><td>HD169830</td><td>90485</td><td>1.40 ± 0.07</td><td>36.8 ± 0.2</td><td>F7V</td><td>5.90</td><td>-4.95<sup>5</sup></td><td>3.4<sup>13</sup></td><td>2.2 ± 0.2<sup>25</sup></td><td>12</td></tr><tr><td>HD179346</td><td>94370</td><td>1.16 ± 0.06</td><td>55.4 ± 0.3</td><td>G0V</td><td>8.05</td><td>-4.84<sup>3</sup></td><td>3.8<sup>13</sup></td><td>0.79 ± 0.68<sup>25</sup></td><td>496</td></tr><tr><td>HD185283</td><td>96861</td><td>0.80 ± 0.04</td><td>31.32 ± 0.02</td><td>K3V</td><td>8.95</td><td>-4.94<sup>3</sup></td><td>0.9<sup>14</sup></td><td>6.0 ± 4.7<sup>25</sup></td><td>9</td></tr></tbody></table>

Notes: ¹Kervella et al. (2022). ²Trifonov et al. (2020). ³Gomes da Silva et al. (2021). ⁴Reiners & Zechmeister (2020). ⁵Liebing et al. (2021). ⁶Stanford-Moore et al. (2020). ⁷Astudillo-Defru et al. (2017). ⁸Yee et al. (2017). ⁹Gondoin (2020). ¹⁰Boro Saikia et al. (2018). ¹¹Hojjatpanah et al. (2020). ¹²Stelzer et al. (2013). ¹³Llorente de Andrés et al. (2021). ¹⁴Hojjatpanah et al. (2020).

¹⁵Ghazaryan et al. (2019). ¹⁶Ammler-von Eiff & Reiners (2012). ¹⁷Rice & Brewer (2020). ¹⁸Borgniet et al. (2017). ¹⁹Costa Silva et al. (2020). ²⁰Hojjatpanah et al. (2019). ²¹Sousa et al. (2018). ²²Brandt (2021). ²³Schneider et al. (2011). ²⁴Maldonado et al. (2020).

²⁵Delgado Mena et al. (2019). ²⁶Ségransan et al. (2011). ²⁷Biller et al. (2022). ²⁸Grieves et al. (2018). ²⁹Adibekyan et al. (2021).
<table><caption>Table A.4: Table A.3 continued.</caption><thead><tr><th>Name</th><th>HIP</th><th>Mass<sup>1</sup><br>(M<sub>☉</sub>)</th><th>Distance<sup>1</sup><br>(pc)</th><th>Spectral<sup>2</sup><br>type</th><th>V<sup>2</sup><br>(mag)</th><th>logR'<sub>HK</sub></th><th>Vsin(I)<br>(km s<sup>-1</sup>)</th><th>Age<br>(Gyr)</th><th>χ<sup>2</sup> <sup>22</sup></th></tr></thead><tbody><tr><td>HD191797</td><td>99695</td><td>0.85 ± 0.04</td><td>43.94 ± 0.06</td><td>K0V</td><td>9.50</td><td>-4.37<sup>3</sup></td><td>3.6<sup>13</sup></td><td>4.6 ± 4.3<sup>25</sup></td><td>6461</td></tr><tr><td>HD195564</td><td>101345</td><td>1.24 ± 0.06</td><td>24.71 ± 0.04</td><td>G2.5IV</td><td>5.65</td><td>-5.18<sup>3</sup></td><td>1.9<sup>13</sup></td><td>7.4 ± 0.3<sup>25</sup></td><td>1054</td></tr><tr><td>HD196050</td><td>101806</td><td>1.27 ± 0.06</td><td>50.49 ± 0.05</td><td>G3V</td><td>7.49</td><td>-4.98<sup>5</sup></td><td>3.3<sup>13</sup></td><td>4.3 ± 0.9<sup>25</sup></td><td>655</td></tr><tr><td>HD204961</td><td>106440</td><td>0.45 ± 0.01</td><td>4.966 ± 0.001</td><td>M2/3V</td><td>8.67</td><td>-5.18<sup>7</sup></td><td>2.1<sup>20</sup></td><td>8.1 ± 4.2<sup>24</sup></td><td>278</td></tr><tr><td>HD207700</td><td>108158</td><td>1.16 ± 0.06</td><td>38.31 ± 0.03</td><td>G4V</td><td>7.43</td><td>-5.10<sup>5</sup></td><td>1.2<sup>13</sup></td><td>9.9 ± 1.6<sup>25</sup></td><td>655</td></tr><tr><td>HD212036</td><td>110440</td><td>1.08 ± 0.05</td><td>52.2 ± 0.3</td><td>G5V</td><td>8.47</td><td>-4.86<sup>6</sup></td><td>1.3<sup>14</sup></td><td>5.7 ± 3.2<sup>25</sup></td><td>894</td></tr><tr><td>HD215497</td><td>112441</td><td>0.92 ± 0.05</td><td>40.53 ± 0.02</td><td>K3V</td><td>8.95</td><td>-5.07<sup>3</sup></td><td>-</td><td>5.8 ± 4.9<sup>29</sup></td><td>4</td></tr><tr><td>HD216437</td><td>113137</td><td>1.22 ± 0.06</td><td>26.68 ± 0.02</td><td>G1V</td><td>6.05</td><td>-5.04<sup>9</sup></td><td>3.7<sup>13</sup></td><td>4.9 ± 1.0<sup>13</sup></td><td>46</td></tr><tr><td>HD221146</td><td>115951</td><td>1.24 ± 0.06</td><td>36.27 ± 0.04</td><td>G1V</td><td>6.89</td><td>-4.98<sup>5</sup></td><td>2.9<sup>13</sup></td><td>6.4 ± 1.1<sup>25</sup></td><td>2220</td></tr><tr><td>HD221420</td><td>116250</td><td>1.34 ± 0.07</td><td>31.14 ± 0.04</td><td>G1V</td><td>6.89</td><td>-5.28<sup>3</sup></td><td>0.7<sup>13</sup></td><td>3.7 ± 0.5<sup>25</sup></td><td>541</td></tr><tr><td>HD221638</td><td>116350</td><td>1.20 ± 0.06</td><td>52.15 ± 0.09</td><td>F6V</td><td>7.55</td><td>-4.86<sup>3</sup></td><td>7.2<sup>13</sup></td><td>2.2 ± 1.4<sup>25</sup></td><td>1586</td></tr><tr><td>HIP9603</td><td>9603</td><td>0.60 ± 0.01</td><td>27.8 ± 0.1</td><td>K7V</td><td>11.10</td><td>-</td><td>-</td><td>-</td><td>6339</td></tr><tr><td>HIP10337</td><td>10337</td><td>0.69 ± 0.01</td><td>23.725 ± 0.010</td><td>K7V</td><td>9.83</td><td>-4.43<sup>10</sup></td><td>3.2<sup>5</sup></td><td>4.0<sup>+9.5<sup>8</sup><sub>-2.8</sub></sup></td><td>62</td></tr><tr><td>HIP39470</td><td>39470</td><td>0.69 ± 0.01</td><td>44.16 ± 0.07</td><td>M0</td><td>10.83</td><td>-</td><td>-</td><td>5.1 ± 4.6<sup>25</sup></td><td>10660</td></tr><tr><td>HIP54597</td><td>54597</td><td>0.80 ± 0.04</td><td>39.51 ± 0.02</td><td>K5V</td><td>9.82</td><td>-5.04<sup>11</sup></td><td>0.5<sup>14</sup></td><td>4.3 ± 4.0<sup>13</sup></td><td>9</td></tr><tr><td>HIP70849</td><td>70849</td><td>0.65 ± 0.01</td><td>24.10 ± 0.01</td><td>K7V</td><td>10.38</td><td>-4.48<sup>10</sup></td><td>0.3<sup>21</sup></td><td>1 - 5<sup>26</sup></td><td>119</td></tr><tr><td>HIP113201</td><td>113201</td><td>0.53 ± 0.01</td><td>23.5 ± 0.2</td><td>M4V</td><td>11.51</td><td>-4.42<sup>7</sup></td><td>-</td><td>1.2 ± 0.1<sup>27</sup></td><td>4869</td></tr></tbody></table>
## Appendix B: Data considered for orbital fits

<table><caption>Table B.1: RV data considered for orbital fits.</caption><thead><tr><th>Star</th><th>Instrument</th><th>ΔT (days)</th><th>N<sub>obs</sub></th><th>Reference</th><th>Star</th><th>Instrument</th><th>ΔT (days)</th><th>N<sub>obs</sub></th><th>Reference</th></tr></thead><tbody><tr><td>GJ660.1</td><td>HARPS</td><td>3649</td><td>59</td><td>ESO archive</td><td>HD140901</td><td>HARPS</td><td>1293</td><td>52</td><td>ESO archive</td></tr><tr><td>GJ680</td><td>HARPS</td><td>3008</td><td>39</td><td>ESO archive</td><td>HD140901</td><td>HARPS</td><td>1293</td><td>52</td><td>ESO archive</td></tr><tr><td>HD1388</td><td>HARPS</td><td>5082</td><td>212</td><td>ESO archive</td><td>HD140901</td><td>AAT</td><td>5878</td><td>321</td><td>Feng et al. (2022)</td></tr><tr><td>HD1388</td><td>HIRES</td><td>8357</td><td>53</td><td>Rosenthal et al. (2021)</td><td>HD140901</td><td>PFS</td><td>3994</td><td>70</td><td>Feng et al. (2022)</td></tr><tr><td>HD16160</td><td>HARPS</td><td>1776</td><td>229</td><td>ESO archive</td><td>HD143361</td><td>HARPS</td><td>4667</td><td>59</td><td>ESO archive</td></tr><tr><td>HD16160</td><td>HIRES</td><td>4858</td><td>146</td><td>Rosenthal et al. (2021)</td><td>HD143361</td><td>CORALIE</td><td>2012</td><td>45</td><td>Jenkins et al. (2017)</td></tr><tr><td>HD16160</td><td>APF</td><td>2216</td><td>252</td><td>Rosenthal et al. (2021)</td><td>HD143361</td><td>MIKE</td><td>2155</td><td>17</td><td>Jenkins et al. (2017)</td></tr><tr><td>HD16160</td><td>LICK</td><td>8153</td><td>81</td><td>Rosenthal et al. (2021)</td><td>HD165131</td><td>HARPS</td><td>5596</td><td>69</td><td>ESO archive</td></tr><tr><td>HD16548</td><td>HARPS</td><td>6333</td><td>52</td><td>ESO archive</td><td>HD167677</td><td>HARPS</td><td>5205</td><td>44</td><td>ESO archive</td></tr><tr><td>HD16548</td><td>LICK</td><td>1349</td><td>9</td><td>Fischer et al. (2014)</td><td>HD169830</td><td>HARPS</td><td>1920</td><td>86</td><td>ESO archive</td></tr><tr><td>HD16905</td><td>HARPS</td><td>6263</td><td>52</td><td>ESO archive</td><td>HD169830</td><td>CORALIE</td><td>1506</td><td>106</td><td>Mayor et al. (2004)*</td></tr><tr><td>HD21175</td><td>HARPS</td><td>1147</td><td>61</td><td>ESO archive</td><td>HD169830</td><td>HIRES</td><td>6959</td><td>85</td><td>Rosenthal et al. (2021)</td></tr><tr><td>HD28254</td><td>HARPS</td><td>6338</td><td>72</td><td>ESO archive</td><td>HD179346</td><td>HARPS</td><td>5861</td><td>44</td><td>ESO archive</td></tr><tr><td>HD42659</td><td>HARPS</td><td>1535</td><td>49</td><td>Hartmann & Hatzes (2015)</td><td>HD185283</td><td>HARPS</td><td>4827</td><td>78</td><td>ESO archive</td></tr><tr><td>HD44219</td><td>HARPS</td><td>6353</td><td>65</td><td>ESO archive</td><td>HD191797</td><td>HARPS</td><td>5523</td><td>42</td><td>ESO archive</td></tr><tr><td>HD46569</td><td>HARPS</td><td>1297</td><td>63</td><td>ESO archive</td><td>HD195564</td><td>HARPS</td><td>4363</td><td>95</td><td>ESO archive</td></tr><tr><td>HD46569</td><td>LC</td><td>1925</td><td>53</td><td>Zechmeister et al. (2013)</td><td>HD195564</td><td>HIRES</td><td>7936</td><td>56</td><td>Rosenthal et al. (2021)</td></tr><tr><td>HD46569</td><td>VLC</td><td>2040</td><td>54</td><td>Zechmeister et al. (2013)</td><td>HD195564</td><td>APF</td><td>1856</td><td>299</td><td>Rosenthal et al. (2021)</td></tr><tr><td>HD56380</td><td>HARPS</td><td>6287</td><td>96</td><td>ESO archive</td><td>HD196050</td><td>HARPS</td><td>6243</td><td>56</td><td>ESO archive</td></tr><tr><td>HD61383</td><td>HARPS</td><td>6293</td><td>63</td><td>ESO archive</td><td>HD196050</td><td>CORALIE</td><td>1256</td><td>31</td><td>Mayor et al. (2004)*</td></tr><tr><td>HD62364</td><td>HARPS</td><td>6233</td><td>87</td><td>ESO archive</td><td>HD196050</td><td>UCLES</td><td>2513</td><td>44</td><td>Butler et al. (2006)*</td></tr><tr><td>HD73256</td><td>HARPS</td><td>2986</td><td>26</td><td>ESO archive</td><td>HD204961</td><td>HARPS</td><td>5864</td><td>183</td><td>ESO archive</td></tr><tr><td>HD73256</td><td>CORALIE</td><td>755</td><td>38</td><td>?*</td><td>HD204961</td><td>UCLES</td><td>3519</td><td>32</td><td>Bailey et al. (2009)*</td></tr><tr><td>HD73256</td><td>HIRES</td><td>4700</td><td>18</td><td>Butler et al. (2017)*</td><td>HD207700</td><td>HARPS</td><td>3260</td><td>36</td><td>ESO archive</td></tr><tr><td>HD75302</td><td>HARPS</td><td>2160</td><td>35</td><td>ESO archive</td><td>HD212036</td><td>HARPS</td><td>5152</td><td>51</td><td>ESO archive</td></tr><tr><td>HD75302</td><td>ELODIE</td><td>3218</td><td>26</td><td>OHP archive</td><td>HD215497</td><td>HARPS</td><td>1842</td><td>99</td><td>?</td></tr><tr><td>HD75302</td><td>SOPHIE</td><td>3227</td><td>20</td><td>OHP archive</td><td>HD216437</td><td>HARPS</td><td>5031</td><td>100</td><td>ESO archive</td></tr><tr><td>HD78612</td><td>HARPS</td><td>3073</td><td>31</td><td>ESO archive</td><td>HD216437</td><td>CORALIE</td><td>1405</td><td>21</td><td>Mayor et al. (2004)*</td></tr><tr><td>HD89839</td><td>HARPS</td><td>6191</td><td>98</td><td>ESO archive</td><td>HD216437</td><td>UCLES</td><td>2801</td><td>39</td><td>Butler et al. (2006)*</td></tr><tr><td>HD93351</td><td>HARPS</td><td>6537</td><td>40</td><td>ESO archive</td><td>HD221146</td><td>HARPS</td><td>3612</td><td>69</td><td>ESO archive</td></tr><tr><td>HD96116</td><td>HARPS</td><td>6296</td><td>40</td><td>ESO archive</td><td>HD221146</td><td>HRS</td><td>6957</td><td>50</td><td>Bowler et al. (2021)</td></tr><tr><td>HD100777</td><td>HARPS</td><td>1231</td><td>37</td><td>ESO archive</td><td>HD221420</td><td>HARPS</td><td>5493</td><td>77</td><td>ESO archive</td></tr><tr><td>HD101198</td><td>HARPS</td><td>1988</td><td>91</td><td>Borgniet et al. (2017)</td><td>HD221420</td><td>AAT</td><td>6515</td><td>88</td><td>Feng et al. (2022)</td></tr><tr><td>HD107094</td><td>HARPS</td><td>6566</td><td>98</td><td>ESO archive</td><td>HD221420</td><td>PFS</td><td>3323</td><td>38</td><td>Feng et al. (2022)</td></tr><tr><td>HD108202</td><td>HARPS</td><td>6088</td><td>42</td><td>ESO archive</td><td>HD221638</td><td>HARPS</td><td>5670</td><td>57</td><td>ESO archive</td></tr><tr><td>HD111031</td><td>HARPS</td><td>3337</td><td>40</td><td>ESO archive</td><td>HIP9603</td><td>HARPS</td><td>5790</td><td>39</td><td>ESO archive</td></tr><tr><td>HD111998</td><td>HARPS</td><td>2992</td><td>150</td><td>Borgniet et al. (2017)</td><td>HIP10337</td><td>HARPS</td><td>6435</td><td>90</td><td>ESO archive</td></tr><tr><td>HD114783</td><td>HARPS</td><td>4415</td><td>91</td><td>ESO archive</td><td>HIP10337</td><td>HIRES</td><td>2161</td><td>18</td><td>Butler et al. (2017)*</td></tr><tr><td>HD114783</td><td>HIRES</td><td>7849</td><td>157</td><td>Rosenthal et al. (2021)</td><td>HIP39470</td><td>HARPS</td><td>5554</td><td>32</td><td>ESO archive</td></tr><tr><td>HD114783</td><td>HRS</td><td>826</td><td>34</td><td>Wittenmyer et al. (2009)</td><td>HIP54597</td><td>HARPS</td><td>5836</td><td>69</td><td>ESO archive</td></tr><tr><td>HD131664</td><td>HARPS</td><td>6100</td><td>65</td><td>ESO archive</td><td>HIP70849</td><td>HARPS</td><td>5521</td><td>59</td><td>ESO archive</td></tr><tr><td>HD135625</td><td>HARPS</td><td>6292</td><td>99</td><td>ESO archive</td><td>HIP113201</td><td>HARPS</td><td>3259</td><td>66</td><td>ESO archive</td></tr></tbody></table>

Notes: *Data retrieved from the DACE archives.
<table><caption>Table B.2: Relative astrometry considered for orbital fits.</caption><thead><tr><th>Companion</th><th>JD - 2400000</th><th>SEP (mas)</th><th>PA (deg)</th><th>Instrument</th><th>Reference</th></tr></thead><tbody><tr><td>GJ680 B</td><td>53491.03</td><td>3760 ± 50</td><td>325.4 ± 0.3</td><td>NACO</td><td>Ward-Duong et al. (2015)</td></tr><tr><td>GJ680 B</td><td>54235.64</td><td>3940 ± 50</td><td>323.7 ± 0.3</td><td>NACO</td><td>Ward-Duong et al. (2015)</td></tr><tr><td>HD1388 B</td><td>57999.50</td><td>1845.9 ± 1.3</td><td>85.8 ± 0.6</td><td>NESSI<sup>a</sup></td><td>Dalba et al. (2021)</td></tr><tr><td>HD1388 B</td><td>53940.39</td><td>1490 ± 10</td><td>85.6 ± 1</td><td>NACO</td><td>This paper</td></tr><tr><td>HD21175 B</td><td>57665.35</td><td>2628.0 ± 2.5</td><td>122.75 ± 0.08</td><td>SPHERE</td><td>This paper</td></tr><tr><td>HD78612 B</td><td>53739.29</td><td>613 ± 5</td><td>272.3 ± 1</td><td>NACO</td><td>This paper</td></tr><tr><td>HD93351 B/C</td><td>56093.02</td><td>1684 ± 30</td><td>76.3 ± 1</td><td>NACO</td><td>This paper</td></tr><tr><td>HD101198 B</td><td>54421.70</td><td>891.2 ± 13.6</td><td>250.36 ± 1.16</td><td>PUEO<sup>b</sup></td><td>Ehrenreich et al. (2010)</td></tr><tr><td>HD101198 B</td><td>54948.30</td><td>926.4 ± 6.3</td><td>251.23 ± 0.55</td><td>NACO</td><td>Ehrenreich et al. (2010)</td></tr><tr><td>HD101198 B</td><td>54948.30</td><td>927.0 ± 10.6</td><td>250.93 ± 0.86</td><td>NACO</td><td>Ehrenreich et al. (2010)</td></tr><tr><td>HD111031 B</td><td>58151.50</td><td>1055.8 ± 1.3</td><td>121.3 ± 0.6</td><td>NESSI</td><td>Dalba et al. (2021)</td></tr><tr><td>HD195564 B</td><td>58006.70</td><td>1113.0 ± 1.3</td><td>49.4 ± 0.6</td><td>NESSI</td><td>Dalba et al. (2021)</td></tr><tr><td>HD195564 B</td><td>58006.70</td><td>1104.9 ± 1.3</td><td>49.8 ± 0.6</td><td>NESSI</td><td>Dalba et al. (2021)</td></tr><tr><td>HD207700 B</td><td>56092.36</td><td>1090 ± 20</td><td>284.2 ± 1</td><td>NACO</td><td>This paper</td></tr><tr><td>HD221146 B</td><td>53262.30</td><td>1723.0 ± 5.0</td><td>28.6 ± 0.1</td><td>NIRC2<sup>c</sup></td><td>Bowler et al. (2021)</td></tr><tr><td>HD221146 B</td><td>53571.14</td><td>1720.0 ± 5.0</td><td>28.5 ± 0.2</td><td>NIRC2</td><td>Bowler et al. (2021)</td></tr><tr><td>HD221146 B</td><td>58041.14</td><td>1623.0 ± 5.0</td><td>25.10 ± 0.14</td><td>NIRC2</td><td>Bowler et al. (2021)</td></tr><tr><td>HD221146 B</td><td>58480.47</td><td>1600.0 ± 5.0</td><td>24.33 ± 0.12</td><td>NIRC2</td><td>Bowler et al. (2021)</td></tr><tr><td>HD221146 B</td><td>58480.47</td><td>1604.0 ± 5.0</td><td>24.26 ± 0.14</td><td>NIRC2</td><td>Bowler et al. (2021)</td></tr><tr><td>HD221146 B</td><td>58676.11</td><td>1592.0 ± 5.0</td><td>24.4 ± 0.2</td><td>NIRC2</td><td>Bowler et al. (2021)</td></tr><tr><td>HIP113201 B</td><td>57649.13</td><td>161.21 ± 0.38</td><td>182.36 ± 0.13</td><td>SPHERE</td><td>Bonavita et al. (2022b)</td></tr><tr><td>HIP113201 B</td><td>57676.03</td><td>167.80 ± 0.53</td><td>180.90 ± 0.53</td><td>SPHERE</td><td>Bonavita et al. (2022b)</td></tr><tr><td>HIP113201 B</td><td>57709.01</td><td>172.61 ± 1.36</td><td>177.72 ± 0.13</td><td>SPHERE</td><td>Bonavita et al. (2022b)</td></tr><tr><td>HIP113201 B</td><td>57906.44</td><td>198.44 ± 1.76</td><td>162.72 ± 0.51</td><td>SPHERE</td><td>Bonavita et al. (2022b)</td></tr><tr><td>HIP113201 B</td><td>58621.39</td><td>285.49 ± 1.29</td><td>123.59 ± 0.26</td><td>SPHERE</td><td>Bonavita et al. (2022b)</td></tr></tbody></table>

Notes: <sup>a</sup> Scott et al. (2018). <sup>b</sup> Rigaut et al. (1998). <sup>c</sup> Service et al. (2016).
<table><caption>Table C.1: Summary of posteriors obtained with our MCMC algorithm for known planets.</caption><thead><tr><th>Parameter</th><th>HD16905 b</th><th>HD28254 b</th><th>HD44219 b*</th><th>HD62364 B</th><th>HD89839 b</th><th>HD100777 b*</th><th>HD111998 b*</th><th>HD114783 c</th><th>HD140901 b</th></tr></thead><tbody><tr><td>a (au)</td><td>8.8<sup>+0.4</sup><sub>-0.3</sub></td><td>2.45<sup>+0.03</sup><sub>-0.04</sub></td><td>1.26<sup>+0.02</sup><sub>-0.07</sub></td><td>6.4 ± 0.1</td><td>4.9 ± 0.1</td><td>1.06<sup>+0.02</sup><sub>-0.01</sub></td><td>1.91<sup>+0.02</sup><sub>-0.03</sub></td><td>5.0 ± 0.1</td><td>11.8<sup>+2.5</sup><sub>-4.1</sub></td></tr><tr><td>P (days)</td><td>10256<sup>+618</sup><sub>-522</sub></td><td>1333 ± 4</td><td>472<sup>+3</sup><sub>-4</sub></td><td>5172 ± 30</td><td>3448 ± 26</td><td>383 ± 1</td><td>810<sup>+6</sup><sub>-5</sub></td><td>4352<sup>+88</sup><sub>-76</sub></td><td>14386<sup>+8099</sup><sub>-4415</sub></td></tr><tr><td>Ecc</td><td>0.68<sup>+0.02</sup><sub>-0.01</sub></td><td>0.95<sup>+0.03</sup><sub>-0.04</sub></td><td>0.41<sup>+0.12</sup><sub>-0.14</sub></td><td>0.610<sup>+0.004</sup><sub>-0.005</sub></td><td>0.20 ± 0.02</td><td>0.38 ± 0.02</td><td>0.13<sup>+0.05</sup><sub>-0.04</sub></td><td>0.05<sup>+0.04</sup><sub>-0.03</sub></td><td>0.77<sup>+0.06</sup><sub>-0.07</sub></td></tr><tr><td>I (°)</td><td>43 ± 3 or 136 ± 3</td><td>21<sup>+38</sup><sub>-11</sub> or 162<sup>+7</sup><sub>-27</sub></td><td>fix to 90</td><td>42 ± 2 or 131 ± 2</td><td>53<sup>+16</sup><sub>-9</sub> or 129<sup>+9</sup><sub>-15</sub></td><td>fix to 90</td><td>fix to 90</td><td>21<sup>+7</sup><sub>-4</sub> or 159<sup>+4</sup><sub>-6</sub></td><td>40<sup>+18</sup><sub>-11</sub> or 138<sup>+12</sup><sub>-21</sub></td></tr><tr><td>Mass (M<sub>lup</sub>)</td><td>11.3<sup>+0.6</sup><sub>-0.7</sub></td><td>3.8<sup>+3.0</sup><sub>-2.0</sub></td><td>0.58<sup>+0.09</sup><sub>-0.08</sub> (&lt; 11.0)</td><td>19<sup>+2</sup><sub>-1</sub></td><td>5.2 ± 0.8</td><td>1.23 ± 0.05 (&lt; 18.0)</td><td>4.8 ± 0.2 (&lt; 12.8)</td><td>1.9<sup>+0.5</sup><sub>-0.4</sub></td><td>1.8 ± 0.5</td></tr><tr><td>Ω (°)</td><td>353<sup>+4</sup><sub>-8</sub> or 125 ± 5</td><td>248<sup>+53</sup><sub>-33</sub> or 163<sup>+38</sup><sub>-36</sub></td><td>fix to 0</td><td>136 ± 3 or 273 ± 5</td><td>195<sup>+9</sup><sub>-10</sub> or 175 ± 9</td><td>fix to 0</td><td>fix to 0</td><td>307<sup>+14</sup><sub>-11</sub> or 355<sup>+11</sup><sub>-9</sub></td><td>153<sup>+14</sup><sub>-16</sub> or 81<sup>+22</sup><sub>-17</sub></td></tr><tr><td>ω (°)</td><td>248 ± 1</td><td>289<sup>+20</sup><sub>-10</sub></td><td>0 ± 1</td><td>0 ± 1</td><td>171<sup>+4</sup><sub>-5</sub></td><td>204<sup>+4</sup><sub>-3</sub></td><td>190<sup>+22</sup><sub>-17</sub></td><td>153<sup>+46</sup><sub>-48</sub></td><td>299<sup>+8</sup><sub>-9</sub></td></tr><tr><td>P<sub>time</sub></td><td>1856 ± 8</td><td>246<sup>+9</sup><sub>-10</sub></td><td>145<sup>+23</sup><sub>-34</sub></td><td>147<sup>+42</sup><sub>-27</sub></td><td>1788<sup>+45</sup><sub>-47</sub></td><td>224 ± 3</td><td>521<sup>+45</sup><sub>-36</sub></td><td>1618<sup>+583</sup><sub>-555</sub></td><td>1182<sup>+47</sup><sub>-46</sub></td></tr><tr><td>√eccos(ω)</td><td>-0.31 ± 0.02</td><td>0.32<sup>+0.30</sup><sub>-0.29</sub></td><td>-0.33<sup>+0.24</sup><sub>-0.17</sub></td><td>0.778 ± 0.003</td><td>-0.44<sup>+0.2</sup><sub>-0.3</sub></td><td>-0.56 ± 0.02</td><td>-0.34<sup>+0.09</sup><sub>-0.07</sub></td><td>-0.17<sup>+0.17</sup><sub>-0.10</sub></td><td>0.42<sup>+0.12</sup><sub>-0.13</sub></td></tr><tr><td>√ecsin(ω)</td><td>-0.766 ± 0.007</td><td>-0.90<sup>+0.14</sup><sub>-0.06</sub></td><td>0.38<sup>+0.16</sup><sub>-0.26</sub></td><td>0.01<sup>+0.01</sup><sub>-0.02</sub></td><td>0.07<sup>+0.04</sup><sub>-0.03</sub></td><td>-0.25<sup>+0.04</sup><sub>-0.03</sub></td><td>-0.06<sup>+0.09</sup><sub>-0.11</sub></td><td>0.08<sup>+0.12</sup><sub>-0.13</sub></td><td>-0.76<sup>+0.07</sup><sub>-0.06</sub></td></tr><tr><td>Jitter (m s<sup>-1</sup>)</td><td>2.9<sup>+0.5</sup><sub>-0.3</sub></td><td>2.9<sup>+0.4</sup><sub>-0.3</sub></td><td>9.7<sup>+1.0</sup><sub>-0.9</sub></td><td>4.0<sup>+0.5</sup><sub>-0.4</sub></td><td>4.3 ± 0.4</td><td>2.4 ± 0.4</td><td>18.7 ± 1.2</td><td>3.2<sup>+0.1</sup><sub>-0.2</sub></td><td>9.2 ± 0.3</td></tr><tr><td>rms (m s<sup>-1</sup>)</td><td>3.6<sup>+0.8</sup><sub>-0.2</sub></td><td>5.8<sup>+0.3</sup><sub>-0.1</sub></td><td>10.8 ± 0.3</td><td>4.9<sup>+0.5</sup><sub>-0.2</sub></td><td>5.7<sup>+0.2</sup><sub>-0.1</sub></td><td>2.6 ± 0.02</td><td>19.4<sup>+0.6</sup><sub>-0.3</sub></td><td>3.9 ± 0.1</td><td>9.4<sup>+0.2</sup><sub>-0.1</sub></td></tr><tr><td>Instrumental offset (m s<sup>-1</sup>)</td><td>H03 = 64804 ± 2<br>H15 = 64817 ± 1</td><td>H03 = -9324 ± 1<br>H15 = -9308<sup>+2</sup><sub>-1</sub></td><td>H03 = -24155<sup>+1</sup><sub>-5</sub><br>H15 = -24211<sup>+6</sup><sub>-5</sub></td><td>H03 = 35510 ± 1<br>H15 = 35525 ± 2</td><td>H03 = 31758 ± 1<br>H15 = 31770<sup>+2</sup><sub>-1</sub></td><td>H03 = 1250 ± 1</td><td>H03 = -3 ± 3</td><td>H03 = 0 ± 1<br>H15 = 0<sup>+2</sup><sub>-1</sub></td><td>H03 = -6701 ± 2<br>AAT = -3 ± 1<br>HRS = 0 ± 1<br>Hir94 = 0 ± 1<br>Hir04 = 0 ± 1</td></tr></tbody></table>

**Notes:** The median and $1\sigma$ confidence interval are given for each parameter. * corresponds to the solutions found considering only the RV data points. The reported masses are therefore the minimum masses. In these cases, the maximum value of the true mass of the companions, estimated with Gaston, with a confidence index of $3\sigma$, are noted between parentheses.

<table><caption>Table C.2: Table C.1 continued.</caption><thead><tr><th>Parameter</th><th>HD143361 b</th><th>HD167677 b</th><th>HD169830 c</th><th>HD196050 b</th><th>HD204961 b</th><th>HD215497 c*</th><th>HD216437 b</th><th>HD221420 B</th><th>HIP70849 b</th></tr></thead><tbody><tr><td>a (au)</td><td>2.05 ± 0.03</td><td>2.9 ± 0.1</td><td>3.29<sup>+0.04</sup><sub>-0.05</sub></td><td>2.65 ± 0.03</td><td>3.7 ± 0.1</td><td>1.31 ± 0.02</td><td>2.57 ± 0.04</td><td>10.1 ± 0.7</td><td>3.99<sup>+0.06</sup><sub>-0.07</sub></td></tr><tr><td>P (days)</td><td>1040 ± 1</td><td>1804 ± 73</td><td>1832 ± 7</td><td>1392 ± 4</td><td>3853<sup>+51</sup><sub>-47</sub></td><td>566 ± 4</td><td>1359 ± 1</td><td>10066<sup>+929</sup><sub>-1125</sub></td><td>3649 ± 18</td></tr><tr><td>Ecc</td><td>0.198<sup>+0.007</sup><sub>-0.006</sub></td><td>0.35<sup>+0.25</sup><sub>-0.15</sub></td><td>0.23<sup>+0.03</sup><sub>-0.02</sub></td><td>0.20 ± 0.01</td><td>0.05 ± 0.03</td><td>0.47<sup>+0.03</sup><sub>-0.04</sub></td><td>0.31 ± 0.01</td><td>0.12<sup>+0.04</sup><sub>-0.03</sub></td><td>0.65<sup>+0.02</sup><sub>-0.01</sub></td></tr><tr><td>I (°)</td><td>50<sup>+19</sup><sub>-12</sub> or 128<sup>+18</sup><sub>-27</sub></td><td>20<sup>+14</sup><sub>-6</sub> or 155<sup>+6</sup><sub>-12</sub></td><td>23<sup>+12</sup><sub>-6</sub> or 158<sup>+5</sup><sub>-12</sub></td><td>42<sup>+11</sup><sub>-7</sub> or 141<sup>+7</sup><sub>-12</sub></td><td>51 ± 3 or 134 ± 3</td><td>fix to 90</td><td>36<sup>+8</sup><sub>-6</sub> or 148<sup>+5</sup><sub>-7</sub></td><td>14 ± 1 or 162 ± 2</td><td>96 ± 16</td></tr><tr><td>Mass (M<sub>lup</sub>)</td><td>4.9<sup>+1.4</sup><sub>-0.9</sub></td><td>3.7<sup>+1.6</sup><sub>-1.1</sub></td><td>8.9<sup>+2.8</sup><sub>-2.6</sub></td><td>4.7 ± 0.8</td><td>0.99<sup>+0.09</sup><sub>-0.08</sub></td><td>0.35 ± 0.02 (&lt; 5.0)</td><td>4.2 ± 0.7</td><td>22 ± 2</td><td>4.5<sup>+0.4</sup><sub>-0.3</sub></td></tr><tr><td>Ω (°)</td><td>352<sup>+30</sup><sub>-18</sub> or 3.91<sup>+23</sup><sub>-16</sub></td><td>236 ± 24</td><td>158<sup>+24</sup><sub>-29</sub> or 291<sup>+25</sup><sub>-28</sub></td><td>5<sup>+13</sup><sub>-8</sub> or 18<sup>+8</sup><sub>-7</sub></td><td>61<sup>+17</sup><sub>-13</sub> or 265<sup>+12</sup><sub>-15</sub></td><td>fix to 0</td><td>129 ± 8 or 101<sup>+7</sup><sub>-6</sub></td><td>272 ± 3 or 241 ± 3</td><td>35 ± 6</td></tr><tr><td>ω (°)</td><td>242 ± 2</td><td>283 ± 25</td><td>305 ± 10</td><td>178 ± 6</td><td>207<sup>+22</sup><sub>-31</sub></td><td>42 ± 4</td><td>65 ± 2</td><td>217<sup>+24</sup><sub>-10</sub></td><td>182 ± 1</td></tr><tr><td>P<sub>time</sub> (days)</td><td>725<sup>+5</sup><sub>-6</sub></td><td>488<sup>+188</sup><sub>-150</sub></td><td>614<sup>+40</sup><sub>-44</sub></td><td>962<sup>+23</sup><sub>-26</sub></td><td>2696<sup>+54</sup><sub>-58</sub></td><td>431<sup>+7</sup><sub>-6</sub></td><td>667 ± 8</td><td>2624<sup>+53</sup><sub>-52</sub></td><td>2719<sup>+14</sup><sub>-13</sub></td></tr><tr><td>√eccos(ω)</td><td>-0.21<sup>+0.02</sup><sub>-0.01</sub></td><td>0.09<sup>+0.30</sup><sub>-0.23</sub></td><td>0.36<sup>+0.06</sup><sub>-0.07</sub></td><td>-0.44<sup>+0.01</sup><sub>-0.02</sub></td><td>-0.16<sup>+0.14</sup><sub>-0.08</sub></td><td>0.51 ± 0.04</td><td>0.23 ± 0.02</td><td>-0.28<sup>+0.07</sup><sub>-0.07</sub></td><td>-0.81 ± 0.01</td></tr><tr><td>√ecsin(ω)</td><td>-0.39 ± 0.01</td><td>-0.56<sup>+0.15</sup><sub>-0.12</sub></td><td>-0.35 ± 0.05</td><td>0.02<sup>+0.04</sup><sub>-0.05</sub></td><td>-0.13<sup>+0.11</sup><sub>-0.07</sub></td><td>0.46<sup>+0.03</sup><sub>-0.04</sub></td><td>0.50<sup>+0.02</sup><sub>-0.01</sub></td><td>-0.21<sup>+0.03</sup><sub>-0.05</sub></td><td>-0.03<sup>+0.01</sup><sub>-0.02</sub></td></tr><tr><td>Jitter (m s<sup>-1</sup>)</td><td>2.1<sup>+0.3</sup><sub>-0.2</sub></td><td>9 ± 1</td><td>10.4<sup>+0.6</sup><sub>-0.5</sub></td><td>3.9<sup>+0.5</sup><sub>-0.4</sub></td><td>2.3 ± 0.1</td><td>1.5 ± 0.2</td><td>2.4 ± 0.2</td><td>3.4 ± 0.2</td><td>6.7<sup>+0.8</sup><sub>-0.7</sub></td></tr><tr><td>rms (m s<sup>-1</sup>)</td><td>8.0 ± 0.1</td><td>9.3<sup>+2.2</sup><sub>-0.7</sub></td><td>10.7<sup>+0.2</sup><sub>-0.1</sub></td><td>6.7 ± 0.1</td><td>3.7<sup>+0.8</sup><sub>-0.3</sub></td><td>2.1 ± 0.1</td><td>4.3 ± 0.1</td><td>4.5<sup>+2.7</sup><sub>-0.5</sub></td><td>6.9<sup>+0.3</sup><sub>-0.2</sub></td></tr><tr><td>Instrumental offset (m s<sup>-1</sup>)</td><td>H03 = 444.4<sup>+0.4</sup><sub>-0.5</sub><br>H15 = 434<sup>+1</sup><sub>-2</sub><br>C07 = 3 ± 2</td><td>H03 = -57241<sup>+5</sup><sub>-3</sub><br>H15 = -57232 ± 4<br>Hir94 = -2 ± 3</td><td>H03 = -6 ± 2<br>C98 = -1<sup>+2</sup><sub>-1</sub><br>Hir04 = -2<sup>+2</sup><sub>-1</sub></td><td>H03 = 61417 ± 1<br>H15 = -1<sup>+2</sup><sub>-1</sub><br>C98 = 61348<sup>+2</sup><sub>-1</sub><br>Hir04 = -2<sup>+2</sup><sub>-1</sub></td><td>H03 = 13352 ± 1<br>H15 = 13350 ± 1<br>UCLES = 0 ± 1</td><td>H03 = 49307 ± 1<br>H15 = -2209 ± 1<br>UCLES = 0 ± 1</td><td>H03 = -2227 ± 1<br>H15 = 26591<sup>+4</sup><sub>-5</sub><br>AAT = 43<sup>+4</sup><sub>-5</sub><br>C98 = -2276 ± 1</td><td>H03 = 53 ± 1<br>H15 = 65<sup>+2</sup><sub>-1</sub><br>PFS = 25<sup>+4</sup><sub>-6</sub></td></tr></tbody></table>
<table><caption>Table C.3: Orbital parameters and mass of single stellar companions.</caption><thead><tr><th>Star</th><th>a (au)</th><th>Mass (M<sub>☉</sub>)</th><th>Ecc</th><th>I (°)</th><th>Ω (°)</th><th>P<sub>time</sub> (years)</th><th>Methods</th><th>Discovery reference</th></tr></thead><tbody><tr><td>Gl680</td><td>32<sup>+9</sup><sub>-6</sub></td><td>0.178 ± 0.004</td><td>0.75<sup>+0.11</sup><sub>-0.13</sub></td><td>167<sup>+4</sup><sub>-10</sub></td><td>171<sup>+18</sup><sub>-14</sub></td><td>60<sup>+22</sup><sub>-20</sub></td><td>RV+HG+HCI</td><td>Ward-Duong et al. (2015)</td></tr><tr><td>HD1388</td><td>77<sup>+13</sup><sub>-10</sub></td><td>0.474<sup>+0.006</sup><sub>-0.005</sub></td><td>0.43<sup>+0.07</sup><sub>-0.05</sub></td><td>91 ± 2</td><td>86 ± 1</td><td>283<sup>+10</sup><sub>-13</sub></td><td>RV+HG+HCI</td><td>Dalba et al. (2021)</td></tr><tr><td>HD16160</td><td>17.1 ± 0.3</td><td>0.099 ± 0.004</td><td>0.635 ± 0.003</td><td>47 ± 1</td><td>307 ± 1</td><td>133.5 ± 0.3</td><td>RV+HG</td><td>Feng et al. (2021)</td></tr><tr><td>HD16548</td><td>30<sup>+11</sup><sub>-8</sub></td><td>0.51<sup>+0.007</sup><sub>-0.06</sub></td><td>0.66 ± 0.06</td><td>24<sup>+5</sup><sub>-4</sub></td><td>125<sup>+13</sup><sub>-10</sub></td><td>4<sup>+9</sup><sub>-10</sub></td><td>105<sup>+73</sup><sub>-43</sub></td><td>This paper</td></tr><tr><td>HD21175</td><td>43<sup>+6</sup><sub>-4</sub></td><td>1.05 ± 0.05</td><td>0.23 ± 0.13</td><td>27<sup>+9</sup><sub>-6</sub></td><td>93<sup>+248</sup><sub>-61</sub></td><td>284<sup>+247</sup><sub>-9</sub></td><td>134<sup>+60</sup><sub>-24</sub></td><td>This paper</td></tr><tr><td>HD42659</td><td>0.561<sup>+0.007</sup><sub>-0.008</sub></td><td>0.50 ± 0.02</td><td>0.15<sup>+0.02</sup><sub>-0.01</sub></td><td>fix to 90</td><td>fix to 0</td><td>227<sup>+9</sup><sub>-8</sub></td><td>RV</td><td>Hartmann & Hatzes (2015)</td></tr><tr><td>HD46569</td><td>10.4 ± 0.2</td><td>0.29 ± 0.01</td><td>0.580 ± 0.004</td><td>45 ± 1 or 138 ± 1</td><td>252 ± 3 or 313 ± 3</td><td>278.7 ± 0.3</td><td>RV+HG</td><td>Zechmeister et al. (2013)</td></tr><tr><td>HD56380</td><td>4.79 ± 0.07</td><td>0.40 ± 0.01</td><td>0.650 ± 0.001</td><td>6.1 ± 0.1</td><td>314 ± 1</td><td>24.0 ± 0.1</td><td>RV+HG</td><td>This paper</td></tr><tr><td>HD61383</td><td>12.7 ± 0.2</td><td>0.227<sup>+0.007</sup><sub>-0.006</sub></td><td>0.33 ± 0.02</td><td>87 ± 4</td><td>222 ± 2</td><td>169<sup>+2</sup><sub>-1</sub></td><td>RV+HG</td><td>This paper</td></tr><tr><td>HD78612</td><td>18.2<sup>+5.5</sup><sub>-5.6</sub></td><td>0.29 ± 0.01</td><td>0.42<sup>+0.13</sup><sub>-0.17</sub></td><td>48<sup>+11</sup><sub>-9</sub></td><td>95.9 ± 0.8</td><td>11<sup>+22</sup><sub>-19</sub></td><td>RV+HG+HCI</td><td>This paper</td></tr><tr><td>HD96116</td><td>23<sup>+5</sup><sub>-3</sub></td><td>0.41<sup>+0.08</sup><sub>-0.05</sub></td><td>0.40<sup>+0.07</sup><sub>-0.08</sub></td><td>55 ± 10</td><td>86<sup>+2</sup><sub>-3</sub></td><td>67<sup>+22</sup><sub>-16</sub></td><td>RV+HG</td><td>dos Santos et al. (2017)</td></tr><tr><td>HD101198</td><td>26<sup>+2</sup><sub>-1</sub></td><td>0.50 ± 0.02</td><td>0.64<sup>+0.09</sup><sub>-0.12</sub></td><td>108 ± 3</td><td>227 ± 5</td><td>147 ± 6</td><td>RV+HG+HCI</td><td>Borgniet et al. (2017)</td></tr><tr><td>HD107094</td><td>54<sup>+8</sup><sub>-9</sub></td><td>0.28 ± 0.1</td><td>0.81 ± 0.03</td><td>95.4 ± 0.4</td><td>278<sup>+5</sup><sub>-6</sub></td><td>257 ± 1</td><td>RV+HG</td><td>Santos et al. (2011)</td></tr><tr><td>HD111031</td><td>21.1 ± 0.6</td><td>0.129 ± 0.003</td><td>0.71 ± 0.01</td><td>137<sup>+4</sup><sub>-3</sub></td><td>325<sup>+2</sup><sub>-1</sub></td><td>42<sup>+2</sup><sub>-3</sub></td><td>RV+HG+HCI</td><td>Feng et al. (2022)</td></tr><tr><td>HD131664</td><td>3.4 ± 0.1</td><td>0.69 ± 0.01</td><td>0.624 ± 0.002</td><td>67 ± 1</td><td>150<sup>+2</sup><sub>-1</sub></td><td>151.7<sup>+0.5</sup><sub>-0.4</sub></td><td>RV+HG</td><td>Moutou et al. (2009)</td></tr><tr><td>HD179346</td><td>9.4 ± 0.2</td><td>0.192 ± 0.007</td><td>0.78 ± 0.01</td><td>91 ± 6</td><td>315 ± 1</td><td>126.3<sup>+0.2</sup><sub>-0.1</sub></td><td>RV+HG</td><td>This paper</td></tr><tr><td>HD191797</td><td>5.4 ± 0.1</td><td>0.139 ± 0.005</td><td>0.94<sup>+0.02</sup><sub>-0.07</sub></td><td>61<sup>+7</sup><sub>-9</sub></td><td>263<sup>+24</sup><sub>-11</sub></td><td>270 ± 10</td><td>RV+HG</td><td>This paper</td></tr><tr><td>HD195564</td><td>73<sup>+17</sup><sub>-12</sub></td><td>0.42<sup>+0.02</sup><sub>-0.01</sub></td><td>0.94<sup>+0.02</sup><sub>-0.01</sub></td><td>154<sup>+10</sup><sub>-16</sub></td><td>79<sup>+24</sup><sub>-30</sub></td><td>318<sup>+75</sup><sub>-66</sub></td><td>RV+HG+HCI</td><td>Dalba et al. (2021)</td></tr><tr><td>HD207700</td><td>40<sup>+7</sup><sub>-6</sub></td><td>0.13 ± 0.01</td><td>0.28<sup>+0.29</sup><sub>-0.19</sub></td><td>104<sup>+4</sup><sub>-7</sub></td><td>35<sup>+4</sup><sub>-7</sub></td><td>181.1<sup>+0.06</sup><sub>-0.05</sub></td><td>RV+HG</td><td>This paper</td></tr><tr><td>HD212036</td><td>4.1 ± 0.1</td><td>0.164<sup>+0.006</sup><sub>-0.005</sub></td><td>0.51 ± 0.01</td><td>142<sup>+12</sup><sub>-15</sub></td><td>350<sup>+18</sup><sub>-69</sub></td><td>126<sup>+10</sup><sub>-66</sub></td><td>RV+HG+HCI</td><td>This paper</td></tr><tr><td>HD221146</td><td>38.0<sup>+3.3</sup><sub>-1.3</sub></td><td>0.633<sup>+0.009</sup><sub>-0.008</sub></td><td>0.89<sup>+0.02</sup><sub>-0.07</sub></td><td>142<sup>+12</sup><sub>-15</sub></td><td>124<sup>+2</sup><sub>-2</sub></td><td>196.7<sup>+0.9</sup><sub>-0.8</sub></td><td>RV+HG</td><td>Bowler et al. (2021)</td></tr><tr><td>HD221638</td><td>7.2 ± 0.2</td><td>0.108 ± 0.005</td><td>0.68 ± 0.02</td><td>149 ± 1</td><td>71 ± 5</td><td>16.71 ± 0.06</td><td>RV+HG</td><td>This paper</td></tr><tr><td>HIP9603</td><td>15 ± 3</td><td>0.21 ± 0.02</td><td>0.63 ± 0.05</td><td>45 ± 1</td><td>260 ± 2 or 142 ± 2</td><td>128 ± 3</td><td>RV+HG</td><td>This paper</td></tr><tr><td>HIP39470</td><td>8.8 ± 0.2 or 7.6 ± 0.1</td><td>0.138<sup>+0.007</sup><sub>-0.008</sub></td><td>0.33<sup>+0.05</sup><sub>-0.06</sub></td><td>22<sup>+2</sup><sub>-3</sub> or 170 ± 1</td><td>201 ± 2</td><td>151 ± 3 or 150 ± 4</td><td>RV+HG</td><td>This paper</td></tr><tr><td>HIP113201</td><td>9.7 ± 0.2</td><td>0.118 ± 0.002</td><td>0.69 ± 0.01</td><td>144.8 ± 0.4</td><td>317 ± 1</td><td>8.96 ± 0.02</td><td>RV+HG+HCI</td><td>Biller et al. (2022)</td></tr></tbody></table>

Notes: HG corresponds to (HIPARCOS/Gaia) absolute astrometry data. For HD42659 B, the reported mass is the minimum mass.
<figure><img src="image_18.png" /><figcaption>Fig. C.1: Semi-major axis and masses of the different stellar companions. For known companions, when previous characterization is available, the dotted lines allow the solutions obtained in this study to be compared with those obtained in previous studies. The dotted gray line corresponds to the hydrogen burning limit (~80 $M_{\text{Jup}}$).</figcaption></figure>
# Appendix D: additional figures

## Appendix D.1. new companions

### Appendix D.1.1. DPASS results

<figure><img src="image_19.png" /><figcaption>Fig. D.1: Periodograms of GJ660.1. *Top left:* Periodogram of the GJ660.1 A RV measurements. *Top right:* Periodogram of the RV residuals. *Bottom left:* Best fit of the GJ660.1 A RV measurements obtained with DPASS. *Bottom left:* Best fit of the GJ660.1 A RV measurements, corrected for the signal of GJ660.1 C, obtained with DPASS.</figcaption></figure>
## Appendix D.1.2. HD75302

<figure><img src="image_20.png" /></figure>

Fig. D.2: Periodograms of HD73256. Top left: Periodogram of the HD73256 RV measurements. Top right: Periodogram of the HD73256 RV measurements, corrected for the signal of HD73256 b. Bottom left: Periodogram of the RV residuals. Bottom right: Fit of the HD73256 RV measurements obtained with DPASS, considering a two-companion orbital fit with a period of 2.5 d and 2785 d, respectively.

<figure><img src="image_21.png" /></figure>

Fig. D.3: Periodograms of HD108202. Left: Periodogram of the HD108202 RV measurements. Middle: Periodogram of the RV residuals. Right: Best fit of the HD108202 RV measurements obtained with DPASS.

## Appendix D.1.3. Corner plots
<figure><img src="image_22.png" /><figcaption>Fig. D.4: Periodograms of HD135625. Left: Periodogram of the HD135625 RV measurements. Middle: Periodogram of the RV residuals. Right: Best fit of the HD135625 RV measurements obtained with DPASS.</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>Fig. D.5: Periodograms of HD165131. Left: Periodogram of the HD165131 RV measurements. Middle: Periodogram of the RV residuals. Right: Best fit of the HD165131 RV measurements obtained with DPASS.</figcaption></figure>

<figure><img src="image_24.png" /><figcaption>Fig. D.6: Periodograms of HD185283. Left: Periodogram of the HD185283 RV measurements. Middle: Periodogram of the RV residuals. Right: Best fit of the HD185283 RV measurements obtained with DPASS.</figcaption></figure>

<figure><img src="image_25.png" /><figcaption>Fig. D.7: Periodograms of HIP10337. Left: Periodogram of the HIP10337 RV measurements. Middle: Periodogram of the RV residuals. Right: Best fit of the HIP10337 RV measurements obtained with DPASS.</figcaption></figure>
<figure><img src="image_26.png" /><figcaption>Fig. D.8: Periodograms of HIP54597. Left: Periodogram of the HIP54597 RV measurements. Middle: Periodogram of the RV residuals. Right: Best fit of the HIP54597 RV measurements obtained with DPASS.</figcaption></figure>

## Appendix D.2. Plots for known substellar companions

### Appendix D.2.1. Long-period companions

#### Appendix D.2.2. short-period companions
<figure><img src="image_27.png" /></figure>

Fig. D.9: Radial velocity vs. bisector velocity span of HD75302 before (left) and after (right) RV correction of the star's short-term activity. The black line corresponds to the best linear fit. The RV and BVS data, from top to bottom, correspond to the data obtained with the instruments H03, H15, ELODIE, SOPHIE, and SOPHIE+.
<figure><img src="image_28.png" /><figcaption>Fig. D.10: Periodograms of HD75302. *Top left:* Periodogram of the HD75302 RV measurements. *Top right:* Periodogram of the HD75302 RV measurements after correction of the star’s short-term activity. *Bottom left:* Periodogram of the RV corrected residuals. *Bottom right:* Sampling window periodogram of the RV measurements.</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>Fig. D.11: Fit of the HD75302 RV measurements before (left) and after (right) correction of the star’s short-term activity.</figcaption></figure>
<figure><img src="image_30.png" /><figcaption>Fig. D.12: Corner plot of the fitted parameters of GJ660.1 C. V0 and V1 correspond to H03 and H15 dataset.</figcaption></figure>
<figure><img src="image_31.png" /><figcaption>Fig. D.13: Corner plot of the fitted parameters of HD73256 C, considering a prior on *I* between 0 and 90°. V0, V1, V2, and V3 correspond to the C98, Hir94, Hir04, and H03 datasets, respectively.</figcaption></figure>
<figure><img src="image_32.png" /><figcaption>Fig. D.14: Corner plot of the fitted parameters of HD73256 C, considering a prior on *I* between 90 and 180°. V0, V1, V2, and V3 correspond to the C98, Hir94, Hir04, and H03 datasets, respectively.</figcaption></figure>
<figure><img src="image_33.png" /><figcaption>Fig. D.15: Corner plot of the fitted parameters of HD75302 b. V0, V1, V2, V3 and V4 correspond to the H03, H15, ELODIE, SOPHIE, and SOPHIE+ datasets, respectively.</figcaption></figure>
<figure><img src="image_34.png" /><figcaption>Fig. D.16: Corner plot of the fitted parameters of HD108202 b. An offset of 5.7 km/s is added to V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_35.png" /><figcaption>Fig. D.17: Corner plot of the fitted parameters of HD135625 b. An offset of 14.9 km/s is added to V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_36.png" /><figcaption>Fig. D.18: Corner plot of the fitted parameters of HD165131 B. An offset of 46.2 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_37.png" /><figcaption>Fig. D.19: Corner plot of the fitted parameters of HD185283 b. An offset of 51.7 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_38.png" /><figcaption>Fig. D.20: Corner plot of the fitted parameters of HIP10337 b, considering a prior on *I* between 0 and 90°. An offset of 3.2 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability. V2 corresponds to the Hir04 dataset.</figcaption></figure>
<figure><img src="image_39.png" /><figcaption>Fig. D.21: Corner plot of the fitted parameters of HIP10337 b, considering a prior on *I* between 90 and 180°. An offset of 3.2 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability. V2 corresponds to the Hir04 dataset.</figcaption></figure>
<figure><img src="image_40.png" /><figcaption>Fig. D.22: Corner plot of the fitted parameters of HIP54597 b. An offset of 51.75 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_41.png" /></figure>

Fig. D.23: Fit of the orbit of HD16905. Top left: Fit of the HD16905 RV measurements. Top right: Fit of the HD16905 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD16905 b, considering a prior on *I* between 0 and 90°. An offset of 64.8 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.
<figure><img src="image_42.png" /><figcaption>Fig. D.24: Corner plot of the fitted parameters of HD16905 b, considering a prior on *I* between 90 and 180°. An offset of 64.8 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_43.png" /><figcaption>Fig. D.25: Fit of the orbit of HD28254. Top left: Fit of the HD28254 RV measurements. Top right: Fit of the HD28254 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD28254 b, considering a prior on *i* between 0 and 90°. An offset of 9.3 km/s is added to V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_44.png" /><figcaption>Fig. D.26: Corner plot of the fitted parameters of HD28254 b, considering a prior on *I* between 90 and 180°. An offset of 9.3 km/s is added to V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_45.png" /><figcaption>Fig. D.27: Fit of the orbit of HD62364. Top left: Fit of the HD62364 RV measurements. Top right: Fit of the HD62364 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD62364 B, considering a prior on *I* between 0 and 90°. An offset of 35.5 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_46.png" /><figcaption>Fig. D.28: Corner plot of the fitted parameters of HD62364 B, considering a prior on *I* between 90 and 180°. An offset of 35.5 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_47.png" /></figure>

Fig. D.29: Fit of the orbit of HD89839. Top left: Fit of the HD89839 RV measurements. Top right: Fit of the HD89839 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the Hipparcos and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD89839 b, considering a prior on *I* between 0 and 90°. An offset of 31.7 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.
<figure><img src="image_48.png" /><figcaption>Fig. D.30: Corner plot of the fitted parameters of HD89839 b, considering a prior on *I* between 90 and 180°. An offset of 31.7 km/s is subtracted from V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_49.png" /></figure>

Fig. D.31: Fit of the orbit of HD114783. Top left: Fit of the HD114783 RV measurements after subtracting the signal of the inner planet. Top right: Fit of the HD114783 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD114783 c, considering a prior on *i* between 0 and 90°. V0, V1, V2, V3, and V4 correspond to H03, H15, HRS, Hir94, and Hir04, respectively.
<figure><img src="image_50.png" /><figcaption>Fig. D.32: Corner plot of the fitted parameters of HD114783 c, considering a prior on *I* between 90 and 180°. V0, V1, V2, V3, and V4 correspond to H03, H15, HRS, Hir94, and Hir04, respectively.</figcaption></figure>
<figure><img src="image_51.png" /><figcaption>Fig. D.33: Fit of the orbit of HD140901. Top left: Fit of the HD140901 RV measurements. Top right: Fit of the HD140901 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the Hipparcos and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD140901 b, considering a prior on *i* between 0 and 90°. V0, V1, V2, and V3 correspond to H03, H15, C07, and MIKE, respectively.</figcaption></figure>
<figure><img src="image_52.png" /><figcaption>Fig. D.34: Corner plot of the fitted parameters of HD140901 b, considering a prior on *I* between 90 and 180°. V0, V1, V2, and V3 correspond to H03, H15, C07, and MIKE, respectively.</figcaption></figure>
<figure><img src="image_53.png" /><figcaption>Fig. D.35: Fit of the orbit of HD143361. Top left: Fit of the HD143361 RV measurements. Top right: Fit of the HD143361 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD143361 b, considering a prior on *I* between 0 and 90°. V0, V1, V2, and V3 correspond to H03, H15, C07, and MIKE, respectively.</figcaption></figure>
<figure><img src="image_54.png" /><figcaption>Fig. D.36: Corner plot of the fitted parameters of HD143361 b, considering a prior on *I* between 90 and 180°. V0, V1, V2, and V3 correspond to H03, H15, C07, and MIKE, respectively.</figcaption></figure>
<figure><img src="image_55.png" /></figure>

Fig. D.37: Fit of the orbit of HD167677. Top left: Fit of the HD167677 RV measurements. Top right: Fit of the HD167677 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the Hipparcos and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD167677 b, considering a prior on *I* between 0 and 90°. An offset of 57.2 km/s is added to V0 (H03) and V1 (H15) to improve readability.
<figure><img src="image_56.png" /><figcaption>Fig. D.38: Corner plot of the fitted parameters of HD167677 b, considering a prior on *I* between 90 and 180°. An offset of 57.2 km/s is added to V0 (H03) and V1 (H15) to improve readability.</figcaption></figure>
<figure><img src="image_57.png" /><figcaption>Fig. D.39: Fit of the orbit of HD169830. Top left: Fit of the HD169830 RV measurements after subtracting the signal of the inner planet. Top right: Fit of the HD169830 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD169830 c, considering a prior on *I* between 0 and 90°. V0, V1, V2, and V3 correspond to C98, Hir94, Hir04, and H03, respectively.</figcaption></figure>
<figure><img src="image_58.png" /><figcaption>Fig. D.40: Corner plot of the fitted parameters of HD169830 c, considering a prior on *I* between 90 and 180°. V0, V1, V2, and V3 correspond to C98, Hir94, Hir04, and H03, respectively.</figcaption></figure>
<figure><img src="image_59.png" /><figcaption>Fig. D.41: Fit of the orbit of HD196050. Top left: Fit of the HD196050 RV measurements. Top right: Fit of the HD196050 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD196050 b, considering a prior on *i* between 0 and 90°. An offset of 61.4 km/s is subtracted from V0 and V1 to improve readability. V0, V1, V2, and V3 correspond to H03, H15, C98, and UCLES, respectively.</figcaption></figure>
<figure><img src="image_60.png" /><figcaption>Fig. D.42: Corner plot of the fitted parameters of HD196050 b, considering a prior on *I* between 90 and 180°. An offset of 61.4 km/s is subtracted from V0 and V1 to improve readability. V0, V1, V2, and V3 correspond to H03, H15, C98, and UCLES, respectively.</figcaption></figure>
<figure><img src="image_61.png" /><figcaption>Fig. D.43: Fit of the orbit of HD204961. Top left: Fit of the HD204961 RV measurements. Top right: Fit of the HD204961 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the Hipparcos and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD204961 b. An offset of 13.0 km/s is subtracted from V0 and V1 to improve readability, considering a prior on *I* between 0 and 90°. V0, V1, and V2 correspond to H03, H15, and UCLES, respectively.</figcaption></figure>
<figure><img src="image_62.png" /><figcaption>Fig. D.44: Corner plot of the fitted parameters of HD204961 b, considering a prior on *I* between 90 and 180°. An offset of 13.0 km/s is subtracted from V0 and V1 to improve readability. V0, V1, and V2 correspond to H03, H15, and UCLES, respectively.</figcaption></figure>
<figure><img src="image_63.png" /><figcaption>Fig. D.45: Fit of the orbit of HD216437. Top left: Fit of the HD216437 RV measurements. Top right: Fit of the HD216437 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD216437 b, considering a prior on *i* between 0 and 90°. An offset of 2.2 km/s is added to V0 and V1 to improve readability. V0, V1, V2, and V3 correspond to H03, H15, UCLES, and C98, respectively.</figcaption></figure>
<figure><img src="image_64.png" /><figcaption>Fig. D.46: Corner plot of the fitted parameters of HD216437 b, considering a prior on *I* between 90 and 180°. An offset of 2.2 km/s is added to V0 and V1 to improve readability. V0, V1, V2, and V3 correspond to H03, H15, UCLES, and C98, respectively.</figcaption></figure>
<figure><img src="image_65.png" /><figcaption>Fig. D.47: Fit of the orbit of HD221420. Top left: Fit of the HD221420 RV measurements. Top right: Fit of the HD221420 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HD221420 B, considering a prior on *I* between 0 and 90°. An offset of 26.5 km/s is subtracted from V0 and V1 to improve readability. V0, V1, V2, and V3 correspond to H03, H15, AAT, and PFS, respectively.</figcaption></figure>
<figure><img src="image_66.png" /><figcaption>Fig. D.48: Corner plot of the fitted parameters of HD221420 B, considering a prior on *I* between 90 and 180°. An offset of 26.5 km/s is subtracted from V0 and V1 to improve readability. V0, V1, V2, and V3 correspond to H03, H15, AAT, and PFS, respectively.</figcaption></figure>
<figure><img src="image_67.png" /><figcaption>Fig. D.49: Fit of the orbit of HIP70849. Top left: Fit of the HIP70849 RV measurements. Top right: Fit of the HIP70849 astrometric acceleration in right ascension (left) and declination (right). The black points correspond to the HIPPARCOS and Gaia EDR3 proper motion measurements. In each plot, the black curve corresponds to the best fit. The color bar indicates the log-likelihood corresponding to the different fits plotted. Bottom: Corner plot of the fitted parameters of HIP70849 b. V0 and V1 correspond to H03 and H15, respectively.</figcaption></figure>
<figure><img src="image_68.png" /><figcaption>Fig. D.50: Fit of RV data points of the stars hosting short-period planets with MCMC. The stars shown are, from left to right and top to bottom, HD44219, HD100777, HD111998, and HD215497.</figcaption></figure>