TESS exoplanet candidates validated with HARPS archival data. A massive Neptune around GJ143 and two Neptunes around HD23472

We aim at the discovery of new planetary systems by exploiting the transit light curve results from TESS orbital observatory's Sector 1 and 2 observations and validating them with precise Doppler measurements obtained from archival HARPS data. Taking advantage of the reported TESS transit events around GJ143 (TOI 186) and HD23472 (TOI 174) we model their HARPS precise Doppler measurements and derive orbital parameters for these two systems. For the GJ143 system TESS has reported only a single transit, and thus its period is unconstrained from photometry. Our RV analysis of GJ143 reveal the full Keplerian solution of the system, which is consistent with an eccentric planet with a mass almost twice that of Neptune and a period of $P_{\rm b}$ = $35.59_{-0.01}^{+0.01}$ days. Our estimates of the GJ143 b planet are fully consistent with the transit timing from TESS. We confirm the two-planet system around HD23472, which according to our analysis is composed of two Neptune mass planets in a possible 5:3 MMR.


Introduction
The Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015) has begun its planet hunt. As of December 2018, the official TESS data release constitutes a total of 54 days of TESS observations from Sectors 1 and 2 taken between 22 July and 21 September. Since then, a number of TESS planet candidates have been confirmed through ground-based Doppler spectroscopy (Huang et al. 2018;Gandolfi et al. 2018;Wang et al. 2019;Jones et al. 2018), and given this high rate of planet detections, a plethora of new planet discoveries during the 2 yr-mission can easily be predicted. In this paper we report the Doppler validation of two additional TESS systems using archival HARPS spectra. We find an eccentric planet with almost two Neptune masses around the K dwarf GJ 143 (TOI 186). Its period of P b = 35.59 +0.01 −0.01 d matches the single transit-timing event from TESS. We provide evidence of the existence of a twoplanet system around HD 23472 (TOI 174) that is likely composed of two Neptune-mass planets in a period ratio of 5:3. This paper is organized as follows: in Sect. 2 we present estimates of the stellar parameters of GJ 143 and HD 23472, we present the available HARPS data and our approach to modeling the radial velocity (RV) data. We present our results for these Tables with the reprocessed HARPS time series are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5)  two exoplanet systems in Sect. 3, and in Sect. 4 we provide a brief summary and conclusions from this study.

GJ 143 and HD 23472 planet host candidates
We derived stellar parameters using the isochrones package (Morton 2015) together with the MIST stellar evolutionary tracks (Dotter 2016), assuming zero extinction. Input data were the BVJK bands taken from Simbad together with the Gaia DR2 parallaxes (which we corrected for the zero-point offset). The inferred parameters and their precision for the two K-dwarfs GJ 143 (≈16 pc) and HD 23472 (≈39 pc) are listed in Table 1. We neglected systematic uncertainties, but a test using a different isochrone set yielded very similar results.

HARPS data
The HARPS spectra of GJ 143 and HD 23472 were reanalyzed by our team prior to the TESS observations with the ultimate goal to reanalyze all publicly available ESO HARPS spectra in order to derive uniformly processed HARPS-Doppler measurements and study possible small but significant instrumentrelated RV systematics (i.e., Tal-Or et al. 2019). We plan to make these data publicly available as a service to the exoplanet community (Trifonov et al., in prep.). All spectra in our MPIA HARPS-archive are reprocessed with the SpEctrum Radial Velocity AnaLyser (SERVAL, Zechmeister et al. 2018) pipeline, which has been demonstrated to produce more precise RV measurements from HARPS spectra than the standard ESO HARPS-DRS pipeline (e.g., see Trifonov et al. 2018;Kaminski et al. 2018). In addition to the precise Doppler  Zechmeister et al. (2018), H α line index, differential line width (dLW), and chromatic index (CRX). The HARPS RVs clearly exhibit a significant peak at 35.6 d, which has no counterpart in the activity indices.
measurements, we inspected the activity index measurements of the H α line, the differential line width of the spectral lines, and the RV chromaticity (wavelength dependence), for more details see Zechmeister et al. (2018).
GJ 143 has 54 HARPS measurements with a mean RV precision of 0.55 m s −1 taken between November 2003 and December 2009, plus four measurements with a mean RV precision of 1.26 m s −1 taken after the HARPS fiber upgrade (since May 2015, see Lo Curto et al. 2015) in the nights from 18 to 26 December 2016. It is well known that after this intervention, HARPS is effectively a new instrument with a notable RV offset between the data taken before and after the fiber exchange. Lo Curto et al. (2015) have estimated that for a star of spectral type K such as GJ 143, the mutual RV offset is likely about 10−12 m s −1 , which makes the detection of low-amplitude RV signals challenging. We therefore decided to use the full HARPS data set of GJ 143, but as a standard practice in this case, we modeled a mutual RV offset between the pre-and post-upgrade RVs.
HD 23472 has only 14 precise HARPS RV measurements, which were taken between February 2004 and February 2013 (before the fiber upgrade). The mean precision of the RV data is 1.37 m s −1 , but the RV scatter is ≈13.5 m s −1 , which can be taken as an indication that planetary companions are present. Of course, scatter like this is common in stars with no obvious planet detections, and the 14 RVs are not sufficient for a period search. Only after taking into account the a priori knowledge from the TESS transit periods and orbital phases is it possible to model these data to search for the presence of planetary companions.
GJ 143 b Fig. 2. HARPS RV data of GJ 143 phased to the best-fit period of P b = 35.589 d. Blue data points show the pre-fiber upgrade data, while red data points represent post-fiber upgrade data. After fitting for the mutual RV offset, these data sets are consistent with each other. The orbital solution is consistent with the TESS transit event. The RV uncertainties include the estimated RV jitter, which was quadratically added to the error budget.

RV modeling
To derive the best-fit model of the RV data, we adopted a maximum likelihood estimator (MLE) scheme coupled with a Nelder-Mead algorithm (Nelder & Mead 1965), which optimizes the likelihood function (− ln L) of a Keplerian or an N-body model. We modeled the standard the RV curve parameters such as semi-amplitude K, period P, eccentricity e, and arguments of periastron ω. Instead of the time of periastron passage t p , we modeled the planetary mean anomalies M for a given epoch. This is a more convenient parameter in the case of an N-body model, which yields osculating orbital parameters of a given epoch. In our study the GJ 143, and HD 23472 mean anomalies were defined at the first HARPS observational epoch for each data set. We assumed edge-on systems, that is, the orbital inclination i was fixed to 90 • , while the longitude of the node Ω was fixed to 0 • . In addition to the fitted parameters, we followed the method of Baluev (2009) and added the unknown RV variance (RV jitter) as an additional parameter. We analyzed the parameter distributions and estimated parameter uncertainties by coupling our MLE fitting code to a L7, page 2 of 5 T. Trifonov et al.: HARPS archival data confirmation of TESS targets  (14) by itself is not sufficient to discover the planets, but taking advantage of the a priory knowledge of the TESS transit ephemerids, it enables constructing a two-planet model, which is in excellent agreement with the RV data (see Table 2 for details). As in Fig. 2, the RV uncertainties include the RV jitter. . More details about our fitting tools, including a GUI interface, can be found online 1 (Trifonov et al., in prep.).

GJ 143
The official TESS Sectors 1 and 2 data release has not provided a transit period of the GJ 143 candidate companion because only one single but significant transit (transit depth = 929 ± 56 ppm, R b = 2.56 ± 0.32 R ⊕ ) of GJ 143 has been observed between 22 July to 21 September for this target. From the transit alone, we 1 https://github.com/3fon3fonov/trifon therefore have no prior knowledge of the period, but the RV data together with the mid-time of the detected transit might in principle reveal the planetary parameters. Figure 1 shows a generalized Lomb-Scargle periodogram (GLS; Zechmeister & Kürster 2009) of the SERVAL time-series products of the GJ 143 HARPS spectra. The RV measurements exhibit a strong peak at 35.62 d with a significant GLS power of 0.6721 that indicates our planetary candidate. There is no significant power at this period in the HARPS-SERVAL activity indicators. It is worth noting that the measurements of the differential line-width activity indicator display significant excess power at very small frequencies and at a period of ≈400 d, the latter confirmed by the H α index. This may mean that GJ 143 is a somewhat active star.
Because the indication of a planet-induced signal is strong, we applied a full Keplerian model to the RV data. However, L7, page 3 of 5 A&A 622, L7 (2019) while we were sampling around the best-fit model, we rejected all the confident MCMC samples whose configurations would not lead to a transit event at BJD = 2458350.312 with a tolerance of ±2 h. We converged to a best-fit solution consistent with a planet with a mass of 30.6 M ⊕ (almost twice the mass of Neptune) orbiting with a period of P b = 35.6 d on a moderately eccentric orbit with e b = 0.33. Orbital parameter estimates and MCMC uncertainties for GJ 143 b are provided in Table 2. The best-fit solution and its confidence region are fully consistent with the transit event reported by TESS.

HD 23472
The two transiting exoplanet candidates around HD 23472 inferred from TESS observations have periods of P b = 17.6800 ± 0.0015 d and P c = 29.8102 ± 0.0047 d, respectively. The durations of both transits are ≈3 ± 0.3 h, with transit depths of 663 ± 48 ppm for the inner and 627 ± 61 ppm for the outer planet. This leads to estimated planetary radii consistent with sub-Neptune-size objects, R b = 1.872 ± 1.321 R ⊕ and R c = 2.149 ± 0.345 R ⊕ , respectively. To measure the most likely masses of the planetary candidates, we modeled the available HARPS data using a two-planet circular Keplerian model, and used the transit information as a prior. In this case, the only unconstrained parameters in this model are the RV semi-amplitudes K b,c (related to the planetary masses) because the transit epoch and periods are well established by the photometry. However, we decided to include the planetary periods and phases in the modeling. Our reason for this is the following: Over the ≈14 years of RV measurements, small deviations of the oscillating orbit may accumulate into strong deviations from an unperturbed sine-like curve in some phases where RV data are obtained. Moreover, although they are precise, the number of TESS transits is limited (only two for HD 23472 c), and thus we do not know the scale of the gravitational interaction in the system, which may have a strong impact on the orbital perturbation and transit timing variations (TTVs). We note that taking the exact period estimates from TESS leads to a very poor fit, and allowing the planetary periods to be adjusted slightly in the modeling is therefore well justified. Our best-fit circular model is shown in Figs. 3 and 4, while Table 2 summarizes its parameters and MCMC uncertainties. We find very consistent estimates of HD 23471 b, whose period is very similar to that found by TESS, whereas for HD 23471 c, the period deviates by about 4 h, which is still within the MCMC-derived errors, however.
For completeness, we also tested a full two-planet dynamical model, which led to a marginally better fit. This fit suggests a somewhat eccentric inner planet with e b = 0.2, and the dynamical integration of this configuration revealed a much larger oscillation of the period ratio of P c /P b from 1.64 to 1.70 with a mean period of 1.67, hence in a potential 5:3 mean motion resonance (MMR). Figure 5 shows the N-body evolution of the P c /P b of the two-planet circular and the full dynamical models within the baseline of the observations (using SyMBA; Duncan et al. 1998). Clearly, both configurations are dynamically active, which shows how challenging the multiplanet modeling of sparse archival data could be because of the unknown magnitude of the dynamical perturbations in the system. The full N-body model is more realistic and may to some extent explain the small deviation of the period of HD 23472 c, but when compared with the TESS results, modeling these sparse 14 HARPS RVs with such a complex model is not justified. We therefore limit our dynamical analysis of the HD 23472 system to an initially circu-  Table 2. Panel b: evolution of a full dynamical model to the HARPS data, which is consistent with a moderately eccentric inner planet at start (e bint. = 0.2) and thus stronger orbital oscillations. The mean period ratio in both cases is slightly above 5:3, which may indicate a resonant system. lar solution (i.e., fixed e b,c = 0, ω b,c = 0), which is already a good model.
We caution that more RV data may be needed to set firm constraints to the orbital parameters and dynamical masses of the HD 23472 system. This can be effectively done by modeling the transit photometry and RVs simultaneously when more data are collected. Until then, the ambiguity in the models continues, although the magnitude of the HD 23472 b and c planetary masses is already determined in the range of a few Earth masses to more likely Neptune-mass planets. This was our scientific goal with this paper.

Conclusions
We reported the confirmation of two new planetary systems discovered with TESS. The GJ 143 system is an interesting case because based on the orbital period of the GJ 143 companion, TESS was able to detect only one single transit, and thus no reported period was available. However, the RV data of GJ 143 together with the transit constraints from TESS revealed an eccentric massive Neptune with 30.6 +2.7 −2.6 M ⊕ and a period of P b = 35.59 +0.01 −0.01 d. We provided strong arguments for a twoplanet system around HD 23472, which is likely composed of two Neptune-mass planets with 17.9 +1.4 −14.0 M ⊕ and 17.2 +1.1 −13.8 M ⊕ , respectively. The two planets are consistent with a period ratio of 5:3, which may indicate a second-order MMR system. GJ 143 would benefit from more photometric data to fully constrain the transit period and the system configuration in general. Similarly, HD 23472 would benefit from more RV measurements, which together with the transit photometry could lead to a more conclusive dynamical analysis, revealing the exact planetary masses and possible resonant motion involved.
We validated the existence of these two TESS planetary system candidates using HARPS archival data taken prior to the TESS observations. This shows the importance of archival Doppler data in the TESS era.