Issue 
A&A
Volume 642, October 2020



Article Number  A121  
Number of page(s)  20  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/202038922  
Published online  15 October 2020 
Planetary system LHS 1140 revisited with ESPRESSO and TESS^{★}
^{1}
Centro de Astrobiología (CAB, CSICINTA), Departamento de Astrofísica,
ESAC campus 28692 Villanueva de la Cañada,
Madrid, Spain
email: Jorge.Lillo@cab.intacsic.es
^{2}
European Southern Observatory,
Alonso de Cordova 3107,
Vitacura,
Region Metropolitana, Chile
^{3}
Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto,
CAUP, Rua das Estrelas,
4150762
Porto, Portugal
^{4}
Physics Institute, Space Research and Planetary Sciences, Center for Space and Habitability – NCCR PlanetS, University of Bern,
Bern, Switzerland
^{5}
AixMarseille Univ, CNRS, CNES, LAM,
Marseille,
France
^{6}
Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto,
Rua do Campo Alegre,
4169007
Porto,
Portugal
^{7}
Geneva Observatory, University of Geneva,
Chemin des Mailettes 51,
1290
Versoix, Switzerland
^{8}
CFisUC, Department of Physics, University of Coimbra,
3004516
Coimbra,
Portugal
^{9}
IMCCE, Observatoire de Paris, PSL University, CNRS, Sorbonne Universiteé,
77 avenue DenfertRochereau,
75014
Paris,
France
^{10}
CNRS, IPAG, Université Grenoble Alpes,
38000
Grenoble, France
^{11}
Departamento de Matemática y Física Aplicadas, Universidad Católica de la Santísima Concepción,
Alonso de Rivera
2850,
Concepción, Chile
^{12}
LATMOS/CNRS/Sorbonne Université/UVSQ,
11 boulevard d’Alembert,
Guyancourt
78280, France
Received:
14
July
2020
Accepted:
24
September
2020
Context. LHS 1140 is an M dwarf known to host two transiting planets at orbital periods of 3.77 and 24.7 days. They were detected with HARPS and Spitzer. The external planet (LHS 1140 b) is a rocky superEarth that is located in the middle of the habitable zone of this lowmass star. All these properties place this system at the forefront of the habitable exoplanet exploration, and it therefore constitutes a relevant case for further astrobiological studies, including atmospheric observations.
Aims. We further characterize this system by improving the physical and orbital properties of the known planets, search for additional planetarymass components in the system, and explore the possibility of coorbitals.
Methods. We collected 113 new highprecision radial velocity observations with ESPRESSO over a 1.5yr time span with an average photonnoise precision of 1.07 m s^{−1}. We performed an extensive analysis of the HARPS and ESPRESSO datasets and also analyzed them together with the new TESS photometry. We analyzed the Bayesian evidence of several models with different numbers of planets and orbital configurations.
Results. We significantly improve our knowledge of the properties of the known planets LHS 1140 b (P_{b} ~ 24.7 days) and LHS 1140 c (P_{c} ~ 3.77 days). We determine new masses with a precision of 6% for LHS 1140 b (6.48 ± 0.46 M_{⊕}) and 9% for LHS 1140 c (m_{c} = 1.78 ± 0.17 M_{⊕}). This reduces the uncertainties relative to previously published values by half. Although both planets have Earthlike bulk compositions, the internal structure analysis suggests that LHS 1140 b might be ironenriched and LHS 1140 c might be a true Earth twin. In both cases, the water content is compatible to a maximum fraction of 10–12% in mass, which is equivalent to a deep ocean layer of 779 ± 650 km for the habitablezone planet LHS 1140 b. Our results also provide evidence for a new planet candidate in the system (m_{d} = 4.8 ± 1.1M_{⊕}) on a 78.9day orbital period, which is detected through three independent methods. The analysis also allows us to discard other planets above 0.5 M_{⊕} for periods shorter than 10 days and above 2 M_{⊕} for periods up to one year. Finally, our coorbital analysis discards coorbital planets in the tadpole and horseshoe configurations of LHS 1140 b down to 1 M_{⊕} with a 95% confidence level (twice better than with the previous HARPS dataset). Indications for a possible coorbital signal in LHS 1140 c are detected in both radial velocity (alternatively explained by a high eccentricity) and photometric data (alternatively explained by systematics), however.
Conclusions. The new precise measurements of the planet properties of the two transiting planets in LHS 1140 as well as the detection of the planet candidate LHS 1140 d make this system a key target for atmospheric studies of rocky worlds at different stellar irradiations.
Key words: planets and satellites: terrestrial planets / planets and satellites: composition / planets and satellites: interiors / planets and satellites: individual: LHS1140 / techniques: radial velocities / techniques: polarimetric
Full Table B.1 is only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/cat/J/A+A/642/A121
© ESO 2020
1 Introduction
In the past decades, the exploration of exoplanets has moved from the detection scheme to the characterization challenge. The new dedicated ground and spacebased facilities built for this purpose now offer the possibility of fully characterizing the properties of extrasolar planets with exquisite precision (e.g., Pepe et al. 2014; Suárez Mascareño et al. 2020; Damasso et al. 2020). Understanding planet formation and evolution, determining atmospheric properties, and ultimately searching for biosignatures requires a sufficiently large sample of planetary systems. The plethora of properties of the current exoplanet population (more than 4100 known so far, according to the NASA Exoplanet Archive, Akeson et al. 2013) indeed offers a large collection of targets to carry out this endeavor. The combinationof different techniques then becomes critical to provide a complete view of the system, which in turn allows the inference of its history through feeding population synthesis models and allowing subsequent atmospheric characterization campaigns. In particular, the combination of the transits and radial velocity techniques is key to understanding the bulk composition of the planets. Exquisite photometric and radial velocity precision is needed to infer its internal structure, however, and key to setting observational constraints on formation and evolution processes.
Planets around lowmass stars have drawn the attention of the community in the recent years due to the groundbased instrumental capabilities in reaching the rocky domain (e.g., Luque et al. 2019; Zechmeister et al. 2019), also in the temperate region around the star (TRAPPIST1 Gillon et al. 2016). The LHS 1140 planetary system is one of these examples. This M4.5 dwarf is located 10.5 parsec away (Gaia Collaboration 2018) and hosts two known planets, a small shortperiod telluric component (LHS 1140 c, P_{c} ~ 3.77 days) and a temperate rocky superEarth (LHS 1140 b, P_{b} ~ 24.7 days). LHS 1140 b was detected by Dittmann et al. (2017) using MEarth (Nutzman & Charbonneau 2008) photometric time series and HARPS (Mayor et al. 2003) radial velocities. Additional HARPS measurements and observations from the Spitzer Space Telescope by Ment et al. (2019) led to the detection of the inner component in the system, LHS 1140 c, which also transits its host star.
LHS 1140 b lies within the habitable zone of the star and possesses a rocky composition, thus representing a key target for further astrobiological studies. The path toward a habitability analysis involves, among others, the study of the internal structure of the planet (Shahar et al. 2019). Its composition and distribution is a direct consequence of its formation (e.g., the abundances of the different elemental building blocks) and evolution (i.e., the different heating and cooling processes, and impacts throughout the planet history), see Dorn et al. (2018) and references there in. The two planets in LHS 1140 with their different orbital (hence irradiation) properties offer a unique opportunity to understand different evolutionary paths in the same environment.
Added to this, LHS 1140 is also an ideal target for searches of coorbital planets. Coorbital configurations consist of planet pairs trapped in gravitationally stable regions and 1:1 mean motion resonances. The Lagrangian points L_{4} /L_{5} present in the gravitational field of a twobody system (like a star and a planet) are stable points of equilibrium that are located exactly on the same orbit as the planet, but ± 60° ahead and behind it. These gravitational wells have been demonstrated to be very stable once an object is trapped (Laughlin & Chambers 2002). The only condition for stability, once trapped, is that the total mass of the planet pair (planet plus trojan) must be lower than 3.8% of the mass of the star (e.g., Gascheau 1843). This relaxed constraint allows similarmass planets to coorbit together in a longterm stable dance around the star. The two planets can share the same orbital path in different coorbital configurations. In the case of circular or quasicircular orbits, the possibilities are tadpole and horseshoe orbits. In the first case, the coorbital surrounds (or librates around) one of the Lagrangian points in the corotating frame with the planet, as in the case of the Jupiter trojans. In the horseshoe configuration, the coorbital describes a horseshoe shape in the corotating frame, moving from L_{4} to L_{5} and exchanging orbits with the planet. This is the case, for instance, of the Saturnian moons Janus and Epimetheus, which have comparable sizes. In the case of eccentric orbits, other configurations come into play, such as the quasisatellite or the antiLagrange configuration (Giuppone et al. 2010). In general, several studies have shown that the formation and stability of coorbital planets makes these configurations not only possible but probable (Cresswell & Nelson 2008, 2009; Leleu et al. 2019a) and that their dynamical properties could have hidden them in the noise of already existing data (Ford & Holman 2007; Madhusudhan & Winn 2009; Janson 2013; Hippke & Angerhausen 2015). However, despite the ample room for stability estimated from theoretical studies, no coorbital pairs have been found so far. Several attempts have been focused on these configurations and have so far set observational constraints to their existence in different regimes (see, e.g., LilloBox et al. 2018a,b), and one key candidate has already been found (TOI178, Leleu et al. 2019a).
The LHS 1140 planetary system has important properties that make it a key system in which to search for these coorbitals. First, the low mass of the host star allows reaching the subEarth mass domain. The transiting nature of the two known planets allows the application of coorbital detection techniques, thus avoiding the degeneracy with the eccentricity at first order (Leleu et al. 2017). The edgeon orientation of the planetary system also permits the search for coorbitals through the transit technique (Janson 2013; LilloBox et al. 2018b). The slow rotational velocity of the star allows precise radial velocity measurements. Moreover, the multiplanet nature of the system increases the likelihood that it hosts coorbital pairs (Leleu et al. 2019b), especially when the planets are in mean motion resonances (Cresswell & Nelson 2008; Leleu et al. 2019b).
We present the joint analysis of the first observations of this system obtained with the Echelle SPectrograph for Rocky Exoplanetsand Stable Spectroscopic Observations (ESPRESSO) instrument and the precise light curve obtained with the Transiting Exoplanets Survey Satellite (TESS) mission. The data are presented in Sect. 2. In Sect. 3 we explore the new radial velocities in search for other planets in the system. The TESS light curve is also further explored in Sect. 4. In Sect. 5 we explore the possibility of coorbitals to the two known planets in the system from different perspectives and with different techniques. The final joint data analysis including the radial velocity and the light curve is presented in Sect. 6. We discuss the result in Sect. 7 and provide the final conclusions in Sect. 8.
2 Observations
2.1 HARPS
LHS 1140 was intensively observed by with the High Accuracy Radial velocity Planet Searcher (HARPS, Mayor et al. 2003) at the 3.6 m telescope of the European Southern Observatory (ESO) La Silla facilities under the program IDs 191.C0873 and 198.C0838 (PI: X. Bonfils) and 0100.C0884 (PI: N. AstudilloDefru). In total, this radial velocity dataset comprised 293 HARPS radial velocity measurements that were presented in Ment et al. (2019), spanning 783 days between November 2015 (JD = 2 457 349.65) and January 2018 (JD = 2 458 133.54). The averagecadence was one epoch every 1.8 days, but normally, two spectra were obtained every night with separations of a few hours. We performed a night binning of the dataset, as suggested in Dumusque et al. (2012). This simple strategy has proved to be efficient to reduce the effect of shortterm correlated noises on orbital elements (e.g., Hara et al. 2019). In total, 145 individual measurements were available. The corresponding uncertainties per binned data point are distributed around 3.2 ±1.2 m s^{−1}. The periodogram of this dataset is presented in Fig. 1.
Fig. 1
Periodogram of the radial velocity (upper panel) for each dataset individually (ESPRESSO in light green and HARPS in red) and for the joint dataset of both instruments assuming a 26 m s^{−1} offset (see Sect. 6). The activity indicators for the individual datasets are also shown in the lower panels. The periods of the three planets are marked as dotted vertical lines, and the rotation period (P_{rot} = 131 days) and its first harmonic (P_{rot}∕2 = 65 days) are marked as vertical dashed lines. 
2.2 ESPRESSO
ESPRESSO is the new ultrastable highresolution spectrograph of the Very Large Telescope at ESO’s Paranal Observatory(Pepe et al. 2020). This facility has the capability of collecting light from any of the four 8.2 m Unit Telescopes (UTs). The instrument has two arms with a total wavelength coverage of 380–788 nm. ESPRESSO is equipped with a powerful reduction pipeline that performs all steps of the basic reduction and provides highlevel data products, including radial velocities. The aim is a final radial velocity precision down to the 10 cm s^{−1} level.
We obtained 116 spectra of LHS 1140 in three consecutive semesters^{1} (P103, P104, and P105). The total time span of the observations is 404 days; they were taken between October 2018 (JD = 2 458 416.71) and December 2019 (JD = 2 458 820.56), and a mean cadence of one spectrum was obtained every 2.2 days (median of one spectrum per night). Each spectrum was obtained with an exposure time of 1820 s in P102 and 1915 s in P103 and P104, and produced a median signaltonoise ratio (S/N) S∕N = 80 at 700 nm. We processed the whole dataset using version v2.0.0 of the ESPRESSO Data reduction Software^{2} (DRS) pipeline (Pepe et al. 2020). We selected a binary mask corresponding to an M5 spectral type for this M4.5 star to perform the crosscorrelation and obtain the final radial velocities (Baranne et al. 1996; Pepe et al. 2003). The resulting radial velocities have a mean photonnoise uncertainty of 1.07 m s^{−1}. In June 2019, ESPRESSO suffered a major intervention, and the chamber and the vacuum vessel were opened to perform a fiberlink exchange. This introduced a small jump in the radial velocity datasets taken before and after the intervention, and the actual value of this jump is not constant for the different spectral types. The data taken before (labeled ESPRESSOpre in this paper) and after (ESPRESSOpost) should therefore be treated as coming from different instruments. This fiberlink exchange occurred between our first and second observation semester. The mean precision corresponding to the data obtained before the fiberlink exchange is 1.16 m s^{−1}, and the postchange data have an average uncertainty of 1.02 m s^{−1}. This improvement goes in line with the increase in transmission measured after the intervention (Pepe et al. 2020) and a slight increase in exposure time from 1820 s in P102 to 1915 s in P103 and P104. We removed two data points with uncertainties above 3 m s^{−1} because low clouds caused a low S/N (the observations were aborted after half the requested exposure time). The final dataset including the radial velocities, activity indicators, S/N, and exposure time for each spectrum is presented in Table B.1. Figure 2 shows the timeseries radial velocity for both HARPS and ESPRESSO, and Fig. 3 shows the phasefolded curves for the two known planets in the system.
The periodogram of the full dataset (including HARPS and ESPRESSO measurements, see Fig. 1) shows clear signals of the known planets LHS 1140 b and LHS 1140 c. The two signals are detectable independently in the two instrument datasets (ESPRESSO and HARPS), and their significance is boosted when they are combined by only including a 26 m s^{−1} radial velocity offset between the two (see Sect. 3.2). In addition to the signals of these two planets, other signals appear at a significant level. The rotation period of the star at P_{rot} = 131 days is clearly detected in each of the two datasets and in the combined set. Interestingly, both datasets also show a peak at half the rotation period P_{rot}∕2 ~ 65 days, which is notas evident in the HARPS radial velocity dataset. However, the rotation period and its alias are clearly visible in the activity index corresponding to the full width at half maximum (FWHM) of the crosscorrelation function (CCF) in both datasets, see Fig. 1. Added to this, additional signals stand out in the region around 70–100 days. This is further investigated in Sect. 3.2.
Taking advantage of the large number of ESPRESSO spectra, we used ODUSSEAS (AntoniadisKarnavas et al. 2020) to estimate theeffective temperature and metallicity of this M star. This machinelearning tool measures pseudoequivalent widths for more than 4000 spectral lines and compares them to a training dataset composed of HARPS spectra of reference Mdwarfstars. Although the code was originally tested for spectra with resolutions from 48 000 to 115 000, we used the combined ESPRESSO spectrum directly with a higher resolution of 140 000. Although the resolution is slightly higher than the highest resolution grid of the code (115 000), we can safely assume that the method that measures the pseudoequivalent widths provides compatible values when a spectrum with higher resolution is used^{3}. With this code we derived an effective temperature of T_{eff} = 2988 ± 67 K and a metallicity of [Fe/H] = −0.262 ± 0.104 dex. These errorswere estimated by considering on the one hand the machinelearning model error (precision errors of 17 K and 0.03 dex, respectively), and on the other hand the mean absolute errors of the machinelearning models for the training dataset (65 K and 0.10 dex for a resolution of 115 000 for the dataset grid, AntoniadisKarnavas et al. 2020).
Fig. 2
Top: radial velocity of LHS 1140 from the HARPS (red) and ESPRESSO (open for ESPRESSOpre and filled for ESPRESSOpost in green) datasets. The black line shows the median radial velocity model from the joint photometric and radial velocity analysis including the two known Keplerian signals and the GP model (see Sect. 6). The gray shaded regions correspond to the 68.7% (dark gray) and 95% (light gray) confidence intervals of the model. The median GP model is shown as a dashed blue line. Middle: radial velocity dataset after removing the median GP model. The Keplerian model is shown as a solid black line. Bottom: radial velocity residuals of the full model. 
Fig. 3
Phasefolded radial velocity signal of LHS 1140 b (left panels) and LHS 1140 c (right panels). For each planet, upper panel: radial velocity data (colorcoded as in the legend), binned radial velocities (open black symbols) with a bin size corresponding to 15% of the phase, the median model for each planet obtained from the joint fit analysis (see Sect. 6) as the solid black line, and the 68.7 and 95% confidence intervals as shaded dark and light gray regions. Bottom panels: residuals of the median model. 
Fig. 4
TPF of LHS 1140 from the TESS observations in Sector 3 (composed with tpfplotter, Aller et al. 2020). The SPOC pipeline aperture is overplotted with shaded red squares, and the Gaia DR2 catalog is also overlaid with symbol sizes proportional to the magnitude contrast with the target, marked with a white cross. 
2.3 TESS photometry
The Transiting Exoplanet Sky Survey (TESS, Ricker et al. 2014) observed LHS 1140 during Sector 3 in camera 1, from 20 September 2018 (JD = 2 459 115) to 18 October 2018 (JD = 2 459 141), immediately before our ESPRESSO campaign. We used tpfplotter^{4} (Aller et al. 2020) to check for contaminant sources in the automatically selected aperture. This is shown in Fig. 4, where we display all sources from the Gaia DR2 catalog with magnitude contrast up to Δm = 8 mag. In addition to LHS 1140, two additional sources lie inside the TESS aperture with magnitude contrasts in the Gaia passband of 3.8 and 4.0 mag. This imposes an upper limit to the dilution factor of 5.2%.
We used the light curve extracted by the SPOC pipeline. For the purpose of this paper, we use the presearch dataconditioning simple aperture photometry (PSDCSAP) detrending of the data, which has a CDPP of 0.547 parts per thousand (hereafter ppt). In Fig. 5 we show the extracted photometric time series and the phasefolded light curves for the two known planets. The photometric precision of this time series makes it sufficient for an independent analysis and characterization of the two planets. We therefore did not use the previous observations from Spitzer and MEarth presented in Ment et al. (2019).
3 Exploring the radial velocity dataset
3.1 ℓ_{1} periodogram
We first analyzed the radial velocity data with the ℓ_{1} periodogram, as defined in Hara et al. (2017). This tool is designed to search for a representation of the signal as a sum of a small number of sinusoids, where “small” is compared to the number of observations. It has a similar aspect to a regular periodogram, but with far fewer peaks due to aliasing. As in Hara et al. (2020), we computed the ℓ_{1} periodogram of the data with different assumptions on the noise covariance. The covariance models were then ranked by crossvalidation. We considered the ESPRESSO and HARPS data without binning.
Similarly to Haywood et al. (2014), we included in the model two activity indicators per instrument, smoothed with Gaussian kernels with different timescales, as linear predictors. These indicators were chosen because they exhibited significant variations on a short (0–2 days) and longer timescale (order of days and 10–20 days). We included the FWHM, with a smoothing timescale of 20 days and 2 days for HARPS and ESPRESSO, respectively, the pipelinederived line asymmetry for ESPRESSO (10 days) and bisector span for HARPS (1 day). The timescales were chosen after fitting the hyperparameters of the Gaussian kernel. We further added a quadratic trend to the model. We assumed a noise model with a whitenoise component σ_{W}, a calibration error σ_{C}, and a Gaussian component with amplitude σ_{R} and exponential decay τ. We then considered all possible combinations of values for σ_{R} = 0.0,0.5, 1.0, 1.5, 2 m s^{−1}, σ_{W} = 0.5, 1 m s^{−1} σ_{C} = 0.1, 0.5, 0.75 m s^{−1}, and τ = 0.0, 1.0, 3.0, 6.0 d. All these noise models were ranked with cross validation, as in Hara et al. (2020). In Fig. 6 we present the ℓ_{1} periodogram obtained with the noise model with maximum crossvalidation score (corresponding to σ_{W} = 1 m s^{−1}, σ_{R} = 0.5 m s^{−1}, τ = 3 d, and σ_{C} = 0.75 m s^{−1}.).
We then considered the 20% highest ranked model CV_{20}, and computed the number of times that a signal is included in the model (which reaches a falsealarm probability, FAP < 0.05). We find that signals might be included in the models at seven periods. In decreasing order of significance, these are at 24.7, 3.77, 78.8, 62, 424, 129, and 16.9 days (see Table 1). The last two signals are included only in 50% of the CV_{20} models, and16.9 days does not appear in the highest ranked noise model. The other signals appear to be present in the data. The signals at 3.77 and 24.7 days correspond to the known planets LHS 1140 c and LHS 1140 b, respectively. The signals at 130 days and 62 days are likely linked to the rotation period of the star (P_{rot}) and half of it, P_{rot}∕2. The origin of the 433 days signal is less certain because its significance also depends on the way the stellar indicators are included in the model. The signal at 78.8 days (the third most significant in the ℓ1periodogram) might be due to a planet and hence deserves additional attention.
Signals appearing in the ℓ_{1} periodogram.
3.2 Evidence for a planet with a period of ~ 80 days
As reported in Ment et al. (2019) and confirmed in the new ESPRESSO dataset (see Sects. 2.2 and 3.1), the periodogram of the joint dataset shows a power excess between the first and second harmonic of the rotational period of the star; this is between 65 and 131 days. These several signals have a maximum peak at around 80 days. We here explorethe evidence for this third signal by studying different scenarios (different number of planets, from one planet to three planets; and different orbital configurations, eccentric or circular) and comparing the models using Bayesian analysis to unveil the significance of this signal.
We explored each dataset separately (i.e., HARPS and ESPRESSO), together with a final full radial velocity dataset analysis (i.e., HARPS+ESPRESSO). Each planet signal was modeled as a Keplerian function, with independent values for the radial velocity semiamplitude of each planet (K_{i}), orbital period (P_{i}), eccentricity (e_{i}), and argument of the periastron (ω_{i}). Because the clear effect of the stellar activity on the radial velocity is shown in the periodogram and has been pointed out by previous studies, we used Gaussian processes (GPs) to model the correlated noise. We used the quasiperiodic kernel from the george (Ambikasaran et al. 2016) implementation, which is a combination of an exponential decay and a periodic part, (1)
Here η_{1} corresponds to the amplitude of the correlated noise, η_{2} can be interpreted as the timescale of the variations of the stellar features causing the correlated noise (see, e.g., Faria et al. 2016), η_{3} represents the stellar rotation period (P_{rot}), and η_{4} is a balance between the exponential and the periodic components of the kernel. As shown in Suárez Mascareño et al. (2020), we can additionally use the FWHM of the crosscorrelation function as an activity indicator to further constrain the GPs (see the FWHM time series in Fig. 7). In this case, the radial velocities and the FWHM share all hyperparameters of the GP, except for the amplitude, therefore an additional parameter per instrument (η_{1,FWHM,j}) is included. For ESPRESSO, which has the highest resolution of the planet searchers and the most stable instrumental profile, the FWHM is expected to be the best tracker of linedeforming stellar activity. This indicator therefore provides an excellent way to determine the radial velocity variations that are caused by this stellar activity. For every instrument, we also added a jitter value for the radial velocity and the FWHM (σ_{RV,j} and σ_{FWHM,j}, respectively) to account for additional uncorrelated noise and a systemic offset (Δ_{RV,j} and Δ_{FWHM,j}).
We used the emcee^{5} (ForemanMackey et al. 2013) implementation of the Goodman & Weare (2010) affine invariant Markov chain Monte Carlo (MCMC) ensemble sampler to explore the parameter space and sample the posterior distribution of each parameter. For each model, we used a number of walkers equal to four times the number of parameters involved (N_{p}) and 20 000 steps for each walker. Following the recommendations from the george package, we first explored the parameter space with a full run using all the steps and walkers. In a second phase, we resampled the walker positions around the best walker and ran a second iteration with the same number of steps to improve convergence. A final run of the maximum probability region was produced by using 10 000 steps per walker. Convergence was checked by estimating the autocorrelation time (τ) of the chains and ensuring a minimum length of the chain, corresponding to 30 × τ. The final joint chain was then composed of 40 000 × N_{p}, which for the simplest model (no planets) corresponds to 6.8 × 10^{5} steps and for the most complex model (including three Keplerians and three instruments) corresponds to 1.3 × 10^{6} steps.
The priors for each parameter involved in the different models are shown in Table B.2. In brief, we assumed Gaussian priors on the planet periods and midtransit times, centered on the values from Ment et al. (2019), but with a broad width of five times the published uncertainties. The remaining Keplerian parameters were set to broad uniform priors throughout the entire allowed regime. The prior on the GP hyperparameter η_{3} was set to a Gaussian centered on the rotational period of the star, corresponding to 131 ± 5 days (Dittmann et al. 2017). Based on the observed radial velocities, we set the η_{1} and η_{1,FWHM} parameters to a uniform distribution with a maximum amplitude of 100 m s^{−1}. The timescale hyperparameter η_{2} was set to a broad range of values between one and five times the stellar rotation period. The instrument systemic velocities were set up with uniform priors between −13.4 and −13.0 km s^{−1}, and the radial velocity jitter for each instrument was allowed to have values up to 20 m s^{−1} (based on our first model attempts).
In total,15 models (including zero to three planets and all combinations of circular and eccentric orbits for each of them) were tested for each of the three datasets (HARPS, ESPRESSO, and HARPS+ESPRESSO), that is, 45 models were tested in total. In order to statistically compare the different models and datasets, we estimated the Bayesian evidence () of each model by using the perrakis^{6} implementation. Based on this Bayesian evidence, we can estimate the Bayes factor () between two models as the ratio between the evidence of each pair of models as . The model with the strongest evidence has the highest statistical relevance, with considered as a strong evidence for one model against the other.
The three panels in Fig. A.1 display this Bayesian evidence for all 15 models considered in each dataset. We find that the twoplanet model with both planets on circular orbits has the strongest evidence in all datasets. This means that the current data are unable to provide evidence for eccentric architectures of the two known planets. We can only place upper limits of e_{b} < 0.096 and e_{c} < 0.274 (at 95% confidence level) based on the noncircular models. The median and 68.7% confidence interval of the parameters for this twoplanet model are presented in Table. B.2.
Although the evidence in the full dataset is weaker, all threeplanet models in the combined HARPS+ESPRESSO dataset converged to a planet candidate with a period of about 80 days, corresponding to the peak in the ℓ1periodogram described in Sect. 3.1 and already seen in the classical periodogram presented in Sect. 2.2 and Fig. 1. The threeplanet model with all planets on circular orbits (labeled 3p1c2c3c) has the second strongest evidence, even above the eccentric twoplanet cases. The median and 68.7% confidence interval of the parameters for this circular threeplanet model are presented in the last column of Table. B.2, and the phasefolded radial velocity curve of this third planet candidate is shown in Fig. 8. It is important to note that this strong evidence is mainly due to the addition of the exquisite ESPRESSO data, which itself show clear signs for this threeplanet model.
We performed a cumulative Bayesian analysis of the dataset by obtaining the Bayes factor for the two most likely models with two and three planets (i.e., the circular case for all planets, labeled 2p1c2c and 3p1c2c3c in Fig. A.1). We obtained the logevidence using up to N_{i} data points in steps of five data points and starting with the first 20 measurements. The results are presented in Fig. 9 and show that the logevidence of the threeplanet model, although the evidence is weaker than for the twoplanet model, increases progressively when more data points are added. This is an indication that additional data might clearly confirm this signal. The jumps in the Bayes factor correspond to the addition of new data from a different instrument, which correspondingly adds new parameters (and so complexity) to the model. The evolution of the Bayesian evidence shows that about 50 new ESPRESSO measurements would be needed to clearly confirm this candidate signal.
According to our threeplanet model in the combined dataset, the third planet would have a minimum mass of M_{⊕} and an orbital period of days. This describes a world in the rockygaseous frontier beyond the habitable zone. We note here that the detection of its radial velocity signal is at the ~ 3.8σ level (). According to the ephemeris found in this radial velocity solution, the TESS observations would have missed the transit of this small planet (in case of a coplanar orbit with the other known planets in the system, the estimated impact parameter is 0.65), which would have occurred about five days before TESS observations started in this sector. The estimated radius of this planet according to forecaster (Chen & Kipping 2017) is R_{⊕}, which placesthe posterior distribution for the planet radius in the middle of the radius valley. The corresponding transit depth would be about 8 ppt, which is comparable with the transit depth of LHS 1140 b and therefore easily detectable by TESS and groundbased instrumentation. TESS will revisit the system in Sector 30 (22 September 2020 to 21 October 2020, in cycle 3). However, given the derived ephemeris for the candidate planet LHS 1140 d, the probability of a transit within the TESS observations in Sector 30 is unfortunately only 0.03% (see the analysis in Sect. 6).
To verify the results, we performed an independent analysis with kima (Faria et al. 2018) on the HARPS+ESPRESSO dataset. We considered up to three Keplerian signals in addition to a GP model for the correlated noise. Broad and equal priors were assigned to the three sets of orbital parameters (i.e., we did not use information from the transits). For the orbital periods we used a loguniform distribution between 1 and 100 days, for the semiamplitudes a uniform prior between 0 and 20 m s^{−1}, and for the eccentricities a Kumaraswamy distribution with shape parameters a = 0.867 and b = 3.03 (Kipping 2013). The priors for the GP and remaining parameters were otherwise very similar to those in Table B.2.
The results indicate a significant detection of only two signals, corresponding to the orbital periods of the two known transiting planets, and the orbital parameters are fully consistent with those found previously. The probability ratio between the threeplanetand twoplanet models is estimated at 1.15, leading to the marginal detection of a third signal, with an orbital period of 78.9 ± 0.5 days. This is compatible with the results explained above. As before, the data do not constrain the eccentricities in the twoplanet model. Consistent results were also obtained for the GP parameters.
Fig. 5
TESS light curve extracted from the SPOC pipeline and detrended using the PDCSAP flux. Upper panel: complete light curve displaying the transits of the two known transiting planets in the system, LHS 1140 b (red, one transit) and LHS 1140 c (light green, five transits). Lower panel: phasefolded light curve centered on the planet phase for each of the two transiting planets, LHS 1140 b (bottom left) and LHS 1140 c (bottom right). The colored symbols correspond to bins of one fifth of the transit duration (i.e., 30 min for LHS 1140 b and 11 min for LHS 1140c). 
Fig. 6
ℓ_{1} periodogram of the LHS 1140 HARPS and ESPRESSO data with the noise model that has the best crossvalidation score. 
Fig. 7
Dataset of the FWHM of the crosscorrelation function from the HARPS and ESPRESSO data. The systemic offset for each instrument as retrieved from the twoplanet model has been removed. The median GP model is shown as the solid black line, and the 1σ and 3σ confidence intervals are shown as shaded dark and light gray regions, respectively. Bottom panel: residuals after removing the GP model. 
Fig. 8
Phasefolded radial velocity signal of the candidate planet LHS 1140 d in the threeplanet model. Upper panel: radial velocity data (colorcoded as in the legend), binned radial velocities (open black symbols) with a bin size corresponding to 15% of the phase, the median model from the joint threeplanet analysis (see Sect. 6) as the solid black line, and the 68.7 and 95% confidence intervals as shaded dark and light gray regions. Bottom panels: residuals of the median model. 
Fig. 9
Difference in the logevidence of the twoplanet and threeplanet models with circular orbits for a cumulative number of data points. Top panel: Bayes factor against the time span of the data set, and bottom panel: against the number of data points. 
3.3 Limits on additional planets at different periods
Because no additional planets can be confirmed with the current dataset, we explored the sensitivity of the HARPS and ESPRESSO data and sampling by injecting planetary signals (we assumed circular orbits) for different periods and masses ranging from 0.3–1000 days and 0.1–300 M_{⊕}. In total, 10^{4} signals were injected in the GPremoved dataset^{7}. Then we proceed in the same manner as described in Sect. 3.2 (but we now removed the GP part) to recover the injected planet signal. We considered the planet to be detected when the median of the posterior distribution for the radial velocity semiamplitude was more than five times the standard deviation away for zero. This analysis is shown in Fig. 10, where we also add the detection limit assuming the scatter of the radial velocity data from the twoplanet model in Sect. 3.2. We conclude that we can discard planets with masses m > 1 M_{⊕} for periods up to 10 days, m > 2 M_{⊕} for periods below 100 days, and m > 5 M_{⊕} for periods shorter than one year. We note that these limits refer to the detection of the planet signal. The statistical significance of the model with respect to the twoplanet model is not assessed in this step. The signal from the planet candidate LHS 1140 d is therefore above the sensitivity line in Fig. 10 (we can detect it in the data), but it cannot be confirmed with sufficient statistical significance.
Fig. 10
Detection limits of the radial velocity dataset for periods shorter than one year. The color code indicates the significance of the radial velocity semiamplitude parameter. The 5σ contour is shown as a solid line, and the limit corresponding to the scatter of the data is shown as a dashed line. The locations of the two known planets are shown as solid circles, and the location of the third planet candidate is shown as an open circle. 
4 Exploring the TESS dataset
The full analysis of the transit signal of LHS 1140 b and LHS 1140 c is presented in Sect. 6. Here, we explore the possibility of additional signals in the TESS dataset.
We first analyzed the TESS dataset to search for transittiming variations (TTVs) of LHS 1140 c, for which five transits are present in the data. We used allesfitter (Günther & Daylan 2020) to model the photometricdata alone. The priors on the parameters were set to Gaussian distributions according to the values published inMent et al. (2019), with a broad width corresponding to ten times the published uncertainties. We also setuniform priors to the individual TTVs of each transit of ± 3.5 h. The results for the inferred TTVs are shown in Table 2. All values are compatible with zero TTVs, except for the first transit, which occurs slightly earlier. However, the TTV is still compatible with zero within 2σ. Additionally, the TESS transit times are also compatible within one minute with the Spitzer transit time obtained six months before by Ment et al. (2019).
We used the transit leastsquares (TLS) software (tls, see Heller et al. 2019) to search for additional signals in the TESS dataset. We masked out the times of transit corresponding to planets LHS 1140 b and LHS 1140 c and performed a smooth detrending of the light curve using wōtan (Hippke et al. 2019) with a window length of 0.5 and using the biweight method. We then computed the TLS periodogram on this masked and detrended light curve over the period range calculated by the tls software based on the data sampling with a range spanning from 0.6 to10 days. The result is shown in Fig. 11. The highest peak in the TLS periodogram corresponds to a periodicity of 1.26 days at a signal detection efficiency (SDE) of 6. Other intriguing signals are also present: a signal close to the orbital period of LHS 1140 c. Although the signal does not reach the statistical significance level of SDE > 7 required to accept a periodicity as statistically valid (see Heller et al. 2019), the fact that it appears at a similar orbital period as LHS 1140 c is intriguing and is further analyzed in Sect. 5 and in particular in Sect. 5.2.
Fig. 11
TLS periodogram of the TESS dataset after masking the light curve from the transit times corresponding to the two known planets LHS 1140 b and LHS 1140 c. The red line indicates the strongest power peak, and the blue line indicates the period of LHS 1140 c. 
5 Coorbital analysis
The long time span and high precision of the radial velocity dataset obtained for this target allows us to perform a detailed exploration for the possible presence of coorbital planets (exotrojans). We explored this scenario for both planets in the system by applying the technique described in Leleu et al. (2017) (hereafter the αtest), a generalization of the technique proposed by Ford & Gaudi (2006). The αtest has been applied to different planetary systems, for instance, by LilloBox et al. (2018a,b), Armstrong et al. (2020), or ToledoPadrón et al. (2020). In this technique, the radial velocity induced by the coorbital system corresponds the sum of two Keplerians with the same orbital period, with at least one of the components transiting its host star (assumed to be the planet). When we also assume that the mass of the planet pair (m_{p} + m_{t}) is far lower than the stellar mass (M_{⋆}) and that in the nearcircular case (typically, e < 0.1), the radial velocity can be approximated to first order in eccentricity and mass by (2)
where n = 2π∕P_{orb}, γ is the systemic velocity, and is the radial velocity semiamplitude of the coorbital pair. The most relevant parameter in this equation is α, which to first order in eccentricity and trojantoplanet mass ratio is (3)
where ζ is the difference in the mean longitude of the two components of the coorbital pair, and t is the time, with the origin set to the midtransit time of the main planet. This α parameter is therefore zero when no coorbital is present (as m_{t} = 0) and can be qualitatively described as the mass ratio between the trojan and the planet. If the radial velocity data are compatible with a nonzero value for this parameter, we can infer a potential mass imbalance between the planet and the locations of the Lagrangian points. A negative value corresponds to L_{4} and a positive value to L_{5}.
5.1 Coorbital analysis for LHS 1140 b
The main interest in this planet resides on its location (in the middle of the habitable zone of this M dwarf) and its physical properties, including a rocky composition and low eccentricity (if any). These properties make this planet ideal for coorbital searches and add the interest that it might lie in the habitable zone. Dvorak et al. (2004) have demonstrated that rocky trojans that coorbit with planets in the habitable zone can also be habitable.
According to the results presented in Sect. 3.2, we here assumed only planets LHS 1140 b and LHS 1140 c in the system. We therefore added another Keplerian signal to Eq. (2) to account for the radial velocity contribution of LHS 1140 c. We followed the same modeling procedure as in Sect. 3.2. We again tested all possible scenarios corresponding to the twoplanet models in that section (i.e., we assumed circular or eccentric orbits for the two planets).
We find the strongest evidence for the circular model for both planets LHS 1140 b and LHS 1140 c. In this model the inferred α parameter is α = 0.015 ± 0.065. The estimated value is fully compatible with zero, clearly stating that no coorbital is detectable down to our sensitivity limits. Instead, using the posterior distribution of this parameter, we can set limits to the presence of coorbitals to this planet at both Lagrangian points. Assuming the 95% confidence level as an upper limit, we can discard coorbital planets to LHS 1140 b more massive than 1.8 M_{⊕} at L_{4} and 1.0 M_{⊕} at L_{5}. Using only the HARPS measurements, we obtained a broader distribution of , which provided upper limits of 2.1 M_{⊕} at L_{4} and 1.8 M_{⊕} at L_{5}. By adding the 116 ESPRESSO data points, we can decrease this limits by > 30%.
We also inspected the TESS light curve around the location of the Lagrangian points L_{4} /L_{5} of LHS 1140 b. Unfortunately, the L_{4} location falls into a gap in the middle of the TESS sector. The transit of L_{5} region occurs at BTJD = 1404.05, about eight hours before a transit of the inner planet LHS 1140 c. At the exact location of the Lagrangian point passage, no transit is found within the photometric precision. The scatter of this region in the light curve is 2.1 ppt (parts per thousand), which would correspond to a 3σ upper limit on any object corresponding to 3 R_{⊕}.
5.2 Coorbital analysis for LHS 1140 c
We proceeded in the same manner as for LHS 1140 b, but now included the α parameter for planet c. In this case, we obtain a value of . The 95% confidence interval is between α ∈ [−0.27, 0.02]. This means that the 95% upper limit (corresponding to 2σ) is only marginally compatible with zero. A further exploration is therefore valuable. Because of its negative value, the mass imbalance producing this radial velocity signal can be interpreted as an existing body located at L_{4}. Assuming the mass of LHS 1140 c as calculated in Sect. 3.2 and a coplanar orbit, the mass of the trojan^{8} would correspond to 0.26± 0.18 M_{⊕}. Despite this ~2σ signal, the Bayesian evidence of the coorbital model is still lower than the twoplanet model. For the HARPSonly dataset, the evidence is stable throughout the entire dataset. However, in the case of the ESPRESSOonly dataset, the evidence for the coorbital model progressively (although slowly) increases toward positive values as more data points are added. It therefore deserves additional attention.
We also performed an independent test to confirm this signal by following the technique proposed by Ford & Gaudi (2006) to compare the time of conjunction measured from the transit signal (T0_{c,LC}) and from the radial velocity signal (T0_{c,RV}). In the case of a single planet in circular orbit, it is easy to see that the lag between both times is zero, Δt = T0_{c,LC} − T0_{c,RV} = 0. However, in the presence of a coorbital body at one of the Lagrangian points, the gravitational pull of the coorbital will change T0_{c,RV} (this occurs earlier if in L_{4} and later if in L_{5}), while leaving the time of transit (T0_{c,LC}) unchanged. We performed this test by analyzing both datasets independently. We modeled the TESS light curve with two planets and the radial velocity dataset with two Keplerian signals (details of the modeling of the photometric dataset are provided in Sect. 6). As a result, we obtain a time lag between the times of conjunction for LHS 1140 c measured from both techniques of Δt = 2.04 ± 1.27 h. This value is again 1.6σ away from zero and provides consistent results with the αtest, an indication for a coorbital body in L_{4} of LHS 1140 c. We note that this time lag can also be explained by a single planet with e cosω ~ 0.15, which is within the confidence interval of the eccentric twoplanet model we showed in Sect. 3.2.
We further investigated this coorbital candidate by inspecting the TESS light curve in this Lagrangian point. Five transits are observed with TESS of LHS 1140 c, but only four passages of its L_{4} in front of the star because the third Lagrangian point lies inside the midsector gap. We explored the TESS light curve in the region of the Lagrangian point after removing the data points inside the transit of planet LHS 1140 b which occurs close to the L_{4} location of transit 4 of LHS 1140 c. The phasefolded light curve in this regime shows a shallow dimming at the location of the L_{4} Lagrangian point (slightly shifted by 0.006 in phase, or equivalently, ~32 min, but perfectly compatible with a dynamically stable libration orbit), see the middle panel in Fig. 12. Interestingly, the duration of the dimming is perfectly compatible with the transit duration of the planet LHS 1140 c, expected for a coorbital in the coplanar case. The depth of this dimming is at the same level (~ 370 ppm) as the photometric scatter in the outoftransit regions, however. We can compare this with the L_{5} region. The photometric scatter in L_{5} is significantly smaller (~280 ppm), and no signs of dimming are seen at the L_{4} depth. This may indicate that the origin of the dimming is related to the larger systematics around the L_{4} region. However, the different indications (including the radial velocity offset, the duration of the dimming corresponding to the planet transit duration, and its location close to the Lagrangian point) encouraged us to further investigate this signal. We therefore modeled it together with the signal of the two known planets to retrieve the size of the potential coorbital candidate. The result provides a radius for this candidate of 0.44 ± 0.08 R_{⊕}. This result increases the certainty in the coorbital scenario. However, the shallow depth (compatible with the photometric noise in the outoftransit regime around L_{4}) prevents us from definitively confirming the nature of this dimming. If confirmed, it might also explain the peak in the TLS periodogram close to the period of LHS 1140 c described in Sect. 4 (see also Fig. 11). Additional data from the TESS extended mission will shed more light on this signal.
Fig. 12
TESS phasefolded light curve around the Lagrangian points L_{5} (top panel) and L_{4} (middle panel) of LHS 1140 c (see Sect. 5.2). The expected location of the transit is marked by vertical dashed lines. The median transit model inferred from the analysis of the L_{4} dimming is shown as a solid red line (middle panel) and a dashed line in the upper panel to guide the eye on the photometric scatter of the L_{5} region. The 1σ confidence interval is shown in the middle panel as a shaded red region. Bottom panel: residuals of the model for the L_{4} region. All panels show bin data points corresponding to 6.6 min (light blue, ten data points inside the transit duration) and 13.2 min (black, five data points inside the transit duration). 
6 Joint photometric and radial velocity model
Based on the analysis performed in the previous section, we proceeded with a joint simultaneous analysis of the HARPS and ESPRESSO radial velocities and the TESS photometric time series to estimate the physical and orbital parameters of the system. To this end,we modeled the radial velocities as explained in Sect. 3.2, but with a slightly different parameterization that now explicitly included the individual contributions of the planet mass and orbital inclination for each planet (instead of using the radial velocity semiamplitude). We also included the GPs in the modeling and constrained them by simultaneously modeling the FWHM of the CCF of the ESPRESSO and HARPS data (see Sect. 3.2).
The TESS photometry was modeled with the batman code (Kreidberg 2015) to retrieve the transit models. We assumed a quadratic limbdarkening model with Gaussian priors around the values calculated by using the limbdarkening code^{9} (Espinoza & Jordán 2015) for the stellar parameters of LHS 1140 published in Ment et al. (2019). We also added photometric jitter to account for underestimated white noise, a mean level parameter to account for imperfect normalization of the light curve, and a dilution factor to account for the contamination in the TESS light curve due to the additional sources in the aperture (see Fig. 4). The stellar mass and radius were also included in the fit to properly account for their uncertainties by using a Gaussian prior around the published values by Ment et al. (2019).
To sample the posterior distributions, we followed the same principles as in Sect. 3.2. We used an MCMC sampler (emcee) with 112 walkers (four times the number of parameters) and 30 000 steps per walker. We finally removed the first half of each chain and combined all chains (1.68 × 10^{6} steps in total) to compute the final posterior distributions. Table B.3 shows the prior and posterior distributions for each parameter. The final models and confidence intervals for the planet transits are shown in Fig. 5, and the corresponding modeling of the radial velocity dataset is shown in Fig. 2 (full time series) and Fig. 3 (phasefolded dataset). In Fig. 7 we additionally show the joint modeling of the FWHM of the CCF sampling the stellar activity. The variations seen in this parameter mimic the variations seen in the radial velocity dataset, showing the correspondence between both measurements coupled by the activity of the star (see also Suárez Mascareño et al. 2020 for a similar analysis on Proxima b).
The inferred masses are compatible with the minimum masses obtained in the radial velocity analysis (Sect. 3.2). We find masses for the two planets that correspond to m_{b} = 6.38 ± 0.45 M_{⊕} and m_{c} = 1.76 ± 0.17 M_{⊕}, slightly lower than the previously reported values. We also find slightly smaller radii for the two planets with the TESS data (R_{b} = 1.635 ± 0.46 R_{⊕}, R_{c} = 1.169 ± 0.038 R_{⊕}) than the values found by Ment et al. (2019) using Spitzer data ( R_{⊕}, R_{c} = 1.282 ± 0.024 R_{⊕}); they are different by about 1.5σ for LHS 1140 b and by 2σ for LHS 1140 c.
We also performed a threeplanet modeling of the full radial velocity and TESS dataset to extract better constraints on the ephemeris of the third planet candidate. The resulting phasefolded radial velocity curve of the third planet is shown in Fig. 8. We can better constrain the time of conjunction for this planet to .
Fig. 13
Massradius diagram for the known exoplanets with the lowest masses with measured mass precisions better than30% (gray symbols). The new locations of the LHS 1140 planets are shown as filled blue (LHS 1140 b) and orange (LHS 1140 c) circles. Their previous locations using the values from Ment et al. (2019) are shown as open squares with dashed error bars and the same color code. The candidate planet at an orbital period of ~ 78 days is shown as the shaded magenta region. Its radius is not known, therefore its vertical location comes from the simple massradiusrelations from Chen & Kipping (2017). The bulk density lines corresponding to different compositions from Zeng et al. (2019) are shown as solid traces, and the dashed lines correspond to isodensities of 1.33, 5.3, and 10 g cm^{−3} (from top tobottom). 
7 Discussion
7.1 Massradius diagram
We estimated the physical and orbital properties of the two known planets with high precision thanks to the precise radial velocity measurements from ESPRESSO and the highcadence and highprecision photometry from TESS. In Fig. 13 we show a massradius diagram including all planets with masses lower than 10 M_{⊕} and radii smaller than2.2 R_{⊕}. The mass and radius precisions are better than 30%. We include the two confirmed planets in the LHS 1140 system with their derived properties estimated in this work and the previously known properties from Ment et al. (2019). The newly derived properties of both planets place them at the top of the Earthlike bulk density line, while previous measurement suggested a slightly lower density for these planets. The planet radius of LHS 1140 b places it close to the radius gap (Fulton et al. 2017), and the newly derived mass confirms the rocky nature of the planet. LHS 1140 c is one of the few planets with a derived mass lower than 2 M_{⊕} in this diagram. The other two planets are GJ 357 b (Luque et al. 2019; Jenkins et al. 2019) and TRAPPIST1 f (Gillon et al. 2016).
7.2 Transit of the planet candidate LHS 1140 d
The massradius diagram in Fig. 13 also shows the location of the third planet candidate (P_{d}~ 78 days) in the system, with an estimated mass of 3.9 M_{⊕}. The uncertainty on the planet radius is taken from the forecaster estimation. The measured mass and estimated radius do not allow us to discern between rocky or gaseous compositions. Based on its mass, the expected planet radius might be in the range 1.2–2.6 R_{⊕}. This impliesa potential transit depth of 3.5 ppt, which is suitable for groundbased instrumentation. Unfortunately, the precision of the time of conjunction from the radial velocity analysis is only about 3.5 days, which means that the ephemeris is too uncertain toplan groundbased observations. The next transits of this planet candidate will occur on Julian dates 2 459 167.8 ± 5.9 (20201114 ± 5.9 days), 2 459 246.5 ± 6.4 (20210201 ± 6.4 days), 2 459 325.2 ± 6.8 (20210420 ± 6.8 days), 2 459 403.9 ± 7.3 (20210708 ± 7.3 days), and 2 459 482.7 ± 7.6 (20210925 ± 7.6 days). When a circular orbit coplanar with the other two known transiting planets is assumed, the impact parameter of this planet candidate would be b∕R_{⋆} ~ 0.65, and it is therefore expected to transit its host star. LHS 1140 will be reobserved by TESS in its extended mission in Sector 30. However, despite the large uncertainty on the time of midtransit, the TESS observations fall in between two transits of the planet and will therefore most likely not be detected by TESS (see Fig. 14). As a consequence, groundbased observations, and especially highprecision photometry, are needed to confirm this planet and further characterize its properties. Alternatively, LHS 1140 is observable with the Cheops spacebased telescope (Broeg et al. 2013), although the large uncertainty on the ephemeris of the transit would require a large effort from the mission.
7.3 Internal structures of LHS 1140 b and LHS 1140 c
To investigate the interior of the confirmed planets, we used the internal structure model developed by Brugger et al. (2017) for the study of terrestrial planets. The input variables of the interior structure model are the total planetary mass, the core mass fraction (CMF), and the water mass fraction (WMF). In order to explore the parameter space, we performed a complete Bayesian analysis to obtain the probability density distributions of the parameters. The Bayesian analysis was carried out by implementing an MCMC algorithm, following the method proposed by Dorn et al. (2015). The initial values of the three input parameters were randomly drawn from their prior distributions, and we used a Gaussian distribution for the mass and two different uniform distributions for the CMF and the WMF. The uniform distribution for the CMF spans from 0 to 0.65, which is a constraint derived from the Fe/Si ratio in the protoSun (Lodders et al. 2009). Including this upper limit, we assumed that these planets have not undergone dramatic processes during or after their formation, such as mantle evaporation or giant impacts. The WMF ranges from 0 to 0.5, the upper value being derived from the composition of the Saturn moon Titan, which is one of the most hydrated Solar System bodies (Tobie et al. 2006). Along with these constraints, we took into account the masses and the radii obtained in this study for LHS1140 b and c to derive the posterior probability distributions of their CMF and WMF. We assumed a surface pressure and temperatureequal to those prevalent on Earth’s surface. This is valid for LHS 1140 b, but we note that LHS 1140 c is closer to the star, and the WMF estimated in this study therefore only represents an upper limit to its water content (which assumes a layer of liquid water, which is not viable in this case). Furthermore, for LHS 1140 c, we extended the interior structure model to include the modeling of water in supercritical phase, whose implementation is described in Mousis et al. (2020). We considered an atmosphere with a composition of 97% water and 3% carbon dioxide, and a surface pressure at its bottom of 300 bar. The atmospheric mass, thickness, albedo, and surface temperature are provided by a grid generated by the atmospheric model described in Marcq et al. (2017).
Figure 15 shows the results of this internal structure analysis by means of ternary diagrams for both planets. In the case of LHS1140 b, the CMF and WMF are estimated to be 0.49 ± 0.07 and 0.03 ± 0.07, respectively.In the case of LHS1140 c, assuming a liquid water layer, the CMF and WMF are estimated to be 0.45 ± 0.10 and 0.05 ± 0.07, respectively.This implies that the bulk composition of planet b is expected to be at least 10% richer in Fe than the Earth’s bulk composition. Planet c might be compatible with a CMF of 0.35 within its 1σ confidence interval, which is very close to the Earth’s 0.32 value. Furthermore, the WMF of LHS 1140 b implies a range of water content from no liquid water to 100 times more water than Earth. The posterior distribution peaks at WMF = 0.04, around 80 times the water content on Earth (0.0005–0.5% – Sotin et al. 2007). The 1% confidence level corresponds to WMF_{1%} = 0.007, still 1.5 times higher than on Earth. When the supercritical phase of water is considered for LHS 1140 c, the CMF is calculated as 0.59 ± 0.05, and the WMF is constrained as . If we were to assume a CMF similar to that of Earth for LHS 1140 c, its WMF would be lower than 0.05%, which means that the inner planetis very poor in water (see Fig. 16).
Fig. 14
Expected posterior distribution of the midtransit time of planet candidate LHS 1140 d (magenta histograms) and the TESS coverage of Sectors 3 and 30 on the field including the host star (shaded gray region). The median of each posterior distribution is marked by a vertical dotted line, and the corresponding calendar date in Universal Time Coordinate system is annotated. 
Fig. 15
Sampled 2D marginal posterior distribution for the CMF and WMF of LHS1140 b. The color code displays the probability density function (PDF). The solid and dashed black lines represent the isoradius curves for the central value of the radius and its 1σ confidence interval limits, respectively. The gray lines delimit the areas we excluded from our sampling (see text). 
Fig. 16
Massradius diagram for different compositions considering water in supercritical phase under the surface temperatureand pressure conditions of LHS 1140 c. The black dot and error bars indicate the position of LHS 1140 c in the massradius diagram. 
8 Conclusions
We have revisited the planetary system LHS 1140 with its two known Earthlike transiting planets by analyzing an intensive campaign taken with the new ultrastable highresolution spectrograph ESPRESSO and observations from Sector 3 of the TESS mission. We have searched for additional signals in the radial velocity data including HARPS and ESPRESSO measurements. The results show additional evidence for the signal of a third planet candidate at ~ 78 days. Our analysis of the joint radial velocity dataset, including a modeling of the stellar activity using the FWHM of the CCF as an activity indicator, shows positive evidence that this signal is of planetary nature (i.e., not caused by activity). We find a ~ 4σ significance in the semiamplitude of the radial velocity of this signal. This is supported by our statistical analysis with kima and the l1periodogram (see Sect. 3.2). The Bayesian evidence of the twoplanet model against the threeplanet model, however, is still not enough to claim a confirmation, but its evolution when the dataset was increased points toward the planetary scenario. This evolution shows that ~ 50 more ESPRESSO measurements are required to unambiguously confirm the planet signal based on radial velocity data alone. The corresponding planet mass of the candidate planet is 3.9 ± 1.1 M_{⊕}. Unfortunately, no transit signal is present in Sector 3 of the TESS data, as expected from the ephemeris of our joint LC+RV analysis. This ephemeris also suggest that the extended TESS mission will not catch the planet transit when the mission revisits this field in Sector 30 between 22 September 2020 and 21 October 2020. If the system is coplanar (as the two known planet inclinations strongly suggest), the impact parameter of this planet candidate would certainly mean that it transits the host star. The expected total transit duration (assuming circular orbit) is about three hours. Further groundbased photometry to search for this transit would therefore be extremely useful and efficient to confirm this signal. We caution about the relatively large uncertainties of the ephemeris, however, which would imply a dedicated campaign of about five consecutive nights on the same target.
We also explored the possibility of coorbital planets in the two known transiting planets through a dedicated study of the radial velocity data and light curve. For LHS 1140 b, the radial velocity data allowed us to constrain the mass of any tadpole and horseshoe coorbitals up to 2.1 M_{⊕} in L_{4} and 1.8 M_{⊕} in L_{5} (95% confidence levels). Similarly, the TESS data discard nearly coplanar coorbitals around the Lagrangian point regions larger than ~ 3 R_{⊕}. On the other hand, for LHS 1140 c, the radial velocity analysis unveils a 2σ signal in its L_{4} region that would correspond to a subEarth coorbital planet with a 0.26 ± 0.18 M_{⊕}. The Bayesian evidence of this model still does not support this scenario against the null hypothesis, however, and several hundred additional measurements with ESPRESSOlike precision are required to confirm the signal with radial velocity data alone. However, the analysis of the TESS light curve indicates a shallow dimming very close to this Lagrangian point with the same duration as the transit of LHS 1140 c (as expected in the coplanar case). If real, this dimming would correspond to a 0.44± 0.08 R_{⊕} subEarth size body. Unfortunately, the significance of the signal is not strong enough to ensure that this dimming is not caused by systematics of the TESS light curve, and only additional data will shed more light on this coorbital candidate.
Finally, we performed a joint analysis including the radial velocity and light curve datasets assuming the twoplanet scenario to derive precise properties of the known transiting planets. We used the largest evidence model from the radial velocityonly analysis, corresponding to the twoplanet case with circular orbits for both planets (this is favored by the Bayesian evidence with a Bayes factor larger than compared to the eccentric case). The results of this analysis show slightly less massive and smaller planets than previously reported values. For LHS 1140 b, we find a planet mass and radius of m_{b} = 6.48 ± 0.46 M_{⊕} and R_{c} = 1.641 ± 0.048 R_{⊕}, resulting in a bulk density of g cm^{−3}. For LHS 1140 c, we obtain m_{b} = 1.77 ± 0.17 M_{⊕} and R_{c} = 1.185 ± 0.044 R_{⊕}, giving g cm^{−3}. This analysis provides unprecedented mass and density precisions for such small rocky planets (9 and 7% for the outer and inner components, respectively). With these new estimates, both planets lie exactly on the Earthlike bulk density line in the massradius diagram (see Fig. 13). Taking advantage of the precise mass and radius measurements, we performed an internal structure analysis considering liquid water conditions and derived a core mass fraction for LHS 1140 b and LHS 1140 c of 49 ± 7 and 45 ± 10%, respectively. It is remarkable that these values are close to the Earth core mass fraction (32%). Additionally, we found a remarkable water content on LHS 1140 b, with the water mass fraction peaking at 4% (80 times more than on Earth). The posterior distribution for this water mass fraction is truncated at zero, however, thus not allowing a definitive confirmation of a water layer on this planet, but instead providing strong indications that LHS 1140 b is a large true water world. This water mass fraction would imply a huge ocean with a depth of 779 ± 650 km according to our internal structure model. On the other hand, if water in supercritical phase is considered for LHS 1140 c, its maximum WMF is estimated to be 0.06%, which shows that planet c is likely to be dry or very poor in water. Thus we can conclude that the system of LHS 1140 presents a gradient of water mass fraction with irradiation, with an inner dry planet and a wet outer planet. LHS 1140 b and LHS 1140 c might have acquired their water content because they formed beyond the snow line and then migrated inward. In the case of planet c, a waterdominated atmosphere could have undergone atmospheric escape due to the high irradiation it received from its host star, which caused it to loose most of its water content. The analysis performed in Sect. 3.3, however, rejects the presence of additional planets more massive than 1 M_{⊕} in between LHS 1140 c and LHS 1140 b that could form a resonant chain that could explain the inward migration of these planets from orbits beyond the snow line (Delisle 2017). Although several mechanisms can explain the disruption of these chains with hot superEarths (Cossou et al. 2014), additional theoretical studies should focus on reconciling the interior structure of these planets (requiring that they formed beyond the snow line) with the fact that the planets in the system are currently not in resonance (usual outcome of inward migration in multiplanet systems, see, e.g., Delisle 2017), and the system hosts no additional components in orbits between the two known planets that can create a resonant chain (see Sect. 3.3).
With our unprecedentedly precise mass and density measurements for the two transiting rocky worlds, LHS 1140 becomes one of the prime systems for astrobiological studies. The habitablezone planet LHS 1140 b is now known to contain a significant amount of water with a relatively high probability. Previous works have studied the actual state of water layers on the surface of this planet by analyzing the possible eyeball, lobster, or snowball scenarios (Yang et al. 2020). Future atmospheric studies with the James Webb Space Telescope and extremely large telescopes will be able to distinguish among these scenarios. These observations will greatly benefit from the precise measurements presented here. Added to this, the possible presence of a third rocky planet in the system that may be transiting the host star at periods beyond the habitable zone makes this system a key target for understanding atmospheric properties of rocky worlds at different stellar irradiations.
Acknowledgements
This research has been funded by the Spanish State Research Agency (AEI) Projects No.ESP201787676C51R and No. MDM20170737 Unidad de Excelencia “María de Maeztu” – Centro de Astrobiología (INTACSIC). This work was supported by FCT – Fundação para a Ciência e a Tecnologia through national funds and by FEDER through COMPETE2020 – Programa Operacional Competitividade e Internacionalização by these grants: UID/FIS/04434/2019; UIDB/04434/2020; UIDP/04434/2020; PTDC/FISAST/32113/2017 & POCI010145FEDER032113; PTDC/FISAST/28953/2017 & POCI010145FEDER028953. N.A.D. acknowledges the support of FONDECYT project 3180063. X.B. and J.M.A. acknowledge funding from the European Research Council under the ERC Grant Agreement n. 337591ExTrA. A.C. acknowledges support by CFisUC projects (UIDB/04564/2020 and UIDP/04564/2020), ENGAGE SKA (POCI010145FEDER022217), and PHOBOS (POCI010145FEDER029932), funded by COMPETE 2020 and FCT, Portugal. J.P.F. is supported in the form of a work contract funded by national funds through FCT with reference DL57/2016/CP1364/CT0005.
Appendix A: Additional figure
Fig. A.1
Bayesian evidence of different models and datasets (from top to bottom: HARPSonly, ESPRESSOonly, and HARPS+ESPRESSO). Each panel contains 15 different models that include a different number of planets: no planets (GPonly, yellow), oneplanet models (pink), twoplanet models (red), and threeplanet model (green). Each has different orbital configuration assumptions (planets with either circular or eccentric orbits). As an example, model 3p1c2c corresponds to a threeplanet model in which planets 1 (LHS 1140 b) and 2 (LHS 1140 c) have assumed circular orbits, but planet 3 (LHS 1140 d) has free eccentricity. The bars corresponding to each model include this information at the bottom of the bar. At the top of each bar we show the Bayes factor of each model compared to the stronggest evidence model for each dataset. The highest bar corresponds to the strongest evidence model. 
Appendix B: Additional tables
Extracted radial velocity, activity indicators, and spectrum properties for the 113 ESPRESSO data points.
References
 Akeson, R. L., Chen, X., Ciardi, D., et al. 2013, PASP, 125, 989 [NASA ADS] [CrossRef] [Google Scholar]
 Aller, A., LilloBox, J., Jones, D., Miranda, L. F., & Barceló Forteza, S. 2020, A&A, 635, A128 [CrossRef] [EDP Sciences] [Google Scholar]
 Ambikasaran, S., ForemanMackey, D., Greengard, L., Hogg, D. W., & O’Neil, M. 2016, IEEE Trans. Pattern Anal. Mach. Intell., 38, 252 [NASA ADS] [CrossRef] [Google Scholar]
 AntoniadisKarnavas, A., Sousa, S. G., DelgadoMena, E., et al. 2020, A&A, 636, A9 [CrossRef] [EDP Sciences] [Google Scholar]
 Armstrong, D. J., Lopez, T. A., Adibekyan, V., et al. 2020, Nature, 583, 39 [Google Scholar]
 Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Broeg, C., Fortier, A., Ehrenreich, D., et al. 2013, Eur. Phys. J. Web Conf., 47, 3005 [Google Scholar]
 Brugger, B., Mousis, O., Deleuil, M., & Deschamps, F. 2017, ApJ, 850, 93 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, J., & Kipping, D. 2017, ApJ, 834, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Cossou, C., Raymond, S. N., Hersant, F., & Pierens, A. 2014, A&A, 569, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cresswell, P., & Nelson, R. P. 2008, A&A, 482, 677 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cresswell, P., & Nelson, R. P. 2009, A&A, 493, 1141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Damasso, M., Sozzetti, A., Lovis, C., et al. 2020, A&A, 642, A31 [CrossRef] [EDP Sciences] [Google Scholar]
 Delisle, J. B. 2017, A&A, 605, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dittmann, J. A., Irwin, J. M., Charbonneau, D., et al. 2017, Nature, 544, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Dorn, C., Khan, A., Heng, K., et al. 2015, A&A, 577, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dorn, C., Bower, D. J., & Rozel, A. 2018, Assessing the Interior Structure of Terrestrial Exoplanets with Implications for Habitability, 66 [Google Scholar]
 Dumusque, X., Pepe, F., Lovis, C., et al. 2012, Nature, 491, 207 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Dvorak, R., PilatLohinger, E., Schwarz, R., & Freistetter, F. 2004, A&A, 426, L37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Espinoza, N., & Jordán, A. 2015, MNRAS, 450, 1879 [NASA ADS] [CrossRef] [Google Scholar]
 Faria, J. P., Haywood, R. D., Brewer, B. J., et al. 2016, A&A, 588, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Faria, J. P., Santos, N. C., Figueira, P., & Brewer, B. J. 2018, J. Open Source Softw., 3, 487 [NASA ADS] [CrossRef] [Google Scholar]
 Ford, E. B., & Gaudi, B. S. 2006, ApJ, 652, L137 [NASA ADS] [CrossRef] [Google Scholar]
 Ford, E. B., & Holman, M. J. 2007, ApJ, 664, L51 [NASA ADS] [CrossRef] [Google Scholar]
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [NASA ADS] [CrossRef] [Google Scholar]
 Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, AJ, 154, 109 [NASA ADS] [CrossRef] [Google Scholar]
 Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gascheau, G. 1843, C. R. Acad. Sci., 16 [Google Scholar]
 Gillon, M., Jehin, E., Lederer, S. M., et al. 2016, Nature, 533, 221 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Giuppone, C. A., Beaugé, C., Michtchenko, T. A., & FerrazMello, S. 2010, MNRAS, 407, 390 [NASA ADS] [CrossRef] [Google Scholar]
 Goodman, J., & Weare, J. 2010, Comm. Appl. Math. Comp. Sci., 5, 65 [CrossRef] [MathSciNet] [Google Scholar]
 Günther, M. N., & Daylan, T. 2020, AAS J., submitted [arXiv:2003.14371] [Google Scholar]
 Hara, N. C., Boué, G., Laskar, J., & Correia, A. C. M. 2017, MNRAS, 464, 1220 [NASA ADS] [CrossRef] [Google Scholar]
 Hara, N. C., Boué, G., Laskar, J., Delisle, J. B., & Unger, N. 2019, MNRAS, 489, 738 [NASA ADS] [CrossRef] [Google Scholar]
 Hara, N. C., Bouchy, F., Stalport, M., et al. 2020, A&A, 636, L6 [CrossRef] [EDP Sciences] [Google Scholar]
 Haywood, R. D., Collier Cameron, A., Queloz, D., et al. 2014, MNRAS, 443, 2517 [NASA ADS] [CrossRef] [Google Scholar]
 Heller, R., Rodenbeck, K., & Hippke, M. 2019, A&A, 625, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hippke, M., & Angerhausen, D. 2015, ApJ, 811, 1 [Google Scholar]
 Hippke, M., David, T. J., Mulders, G. D., & Heller, R. 2019, AJ, 158, 143 [NASA ADS] [CrossRef] [Google Scholar]
 Janson, M. 2013, ApJ, 774, 156 [Google Scholar]
 Jenkins, J. S., Pozuelos, F. J., Tuomi, M., et al. 2019, MNRAS, 490, 5585 [CrossRef] [Google Scholar]
 Kipping, D. M. 2013, MNRAS, 434, L51 [NASA ADS] [CrossRef] [Google Scholar]
 Kreidberg, L. 2015, PASP, 127, 1161 [NASA ADS] [CrossRef] [Google Scholar]
 Laughlin, G., & Chambers, J. E. 2002, AJ, 124, 592 [Google Scholar]
 Leleu, A., Robutel, P., Correia, A. C. M., & LilloBox, J. 2017, A&A, 599, L7 [CrossRef] [EDP Sciences] [Google Scholar]
 Leleu, A., LilloBox, J., Sestovic, M., et al. 2019a, A&A, 624, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leleu, A., Coleman, G. A. L., & Ataiee, S. 2019b, A&A, 631, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 LilloBox, J., Barrado, D., Figueira, P., et al. 2018a, A&A, 609, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 LilloBox, J., Leleu, A., Parviainen, H., et al. 2018b, A&A, 618, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lodders, K., Palme, H., & Gail, H. P. 2009, in LandoltBörnstein, New Series, Vol. VI/4B, Chap. 4.4, ed. J. E. Trümper (SpringerVerlag), 560 [Google Scholar]
 Luque, R., Pallé, E., Kossakowski, D., et al. 2019, A&A, 628, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Madhusudhan, N., & Winn, J. N. 2009, ApJ, 693, 784 [NASA ADS] [CrossRef] [Google Scholar]
 Marcq, E., Salvador, A., Massol, H., & Davaille, A. 2017, J. Geophys. Res. Planets, 122, 1539 [Google Scholar]
 Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
 Ment, K., Dittmann, J. A., AstudilloDefru, N., et al. 2019, AJ, 157, 32 [NASA ADS] [CrossRef] [Google Scholar]
 Mousis, O., Deleuil, M., Aguichine, A., et al. 2020, ApJ, 896, L22 [CrossRef] [Google Scholar]
 Nutzman, P., & Charbonneau, D. 2008, PASP, 120, 317 [NASA ADS] [CrossRef] [Google Scholar]
 Pepe, F., Bouchy, F., Queloz, D., & Mayor, M. 2003, ASP Conf. Ser., 294, 39 [Google Scholar]
 Pepe, F., Ehrenreich, D., & Meyer, M. R. 2014, Nature, 513, 358 [NASA ADS] [CrossRef] [Google Scholar]
 Pepe, F., Damasso, M., Cristiani, S., et al. 2020, A&A, submitted [Google Scholar]
 Perrakis, K., Ntzoufras, I., & Tsionas, E. G. 2014, Comput. Stat. Data Anal., 77, 54 [CrossRef] [Google Scholar]
 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, SPIE Conf. Ser., 9143, 20 [Google Scholar]
 Shahar, A., Driscoll, P., Weinberger, A., & Cody, G. 2019, Science, 364, 434 [CrossRef] [Google Scholar]
 Sotin, C., Grasset, O., & Mocquet, A. 2007, Icarus, 191, 337 [NASA ADS] [CrossRef] [Google Scholar]
 Suárez Mascareño, A., Faria, J. P., Figueira, P., et al. 2020, A&A, 639, A77 [CrossRef] [EDP Sciences] [Google Scholar]
 Tobie, G., Lunine, J. I., & Sotin, C. 2006, Nature, 440, 61 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 ToledoPadrón, B., Lovis, C., Suárez Mascareño, A., et al. 2020, A&A, 641, A92 [CrossRef] [EDP Sciences] [Google Scholar]
 Yang, J., Ji, W., & Zeng, Y. 2020, Nat. Astron., 4, 58 [CrossRef] [Google Scholar]
 Zechmeister, M., Dreizler, S., Ribas, I., et al. 2019, A&A, 627, A49 [NASA ADS] [EDP Sciences] [Google Scholar]
 Zeng, L., Jacobsen, S. B., Sasselov, D. D., et al. 2019, Proc. Natl. Acad. Sci., 116, 9723 [Google Scholar]
https://github.com/exord/bayev. This code is a python implementation by R. Díaz of the formalism explained in Perrakis et al. (2014).
All Tables
Extracted radial velocity, activity indicators, and spectrum properties for the 113 ESPRESSO data points.
All Figures
Fig. 1
Periodogram of the radial velocity (upper panel) for each dataset individually (ESPRESSO in light green and HARPS in red) and for the joint dataset of both instruments assuming a 26 m s^{−1} offset (see Sect. 6). The activity indicators for the individual datasets are also shown in the lower panels. The periods of the three planets are marked as dotted vertical lines, and the rotation period (P_{rot} = 131 days) and its first harmonic (P_{rot}∕2 = 65 days) are marked as vertical dashed lines. 

In the text 
Fig. 2
Top: radial velocity of LHS 1140 from the HARPS (red) and ESPRESSO (open for ESPRESSOpre and filled for ESPRESSOpost in green) datasets. The black line shows the median radial velocity model from the joint photometric and radial velocity analysis including the two known Keplerian signals and the GP model (see Sect. 6). The gray shaded regions correspond to the 68.7% (dark gray) and 95% (light gray) confidence intervals of the model. The median GP model is shown as a dashed blue line. Middle: radial velocity dataset after removing the median GP model. The Keplerian model is shown as a solid black line. Bottom: radial velocity residuals of the full model. 

In the text 
Fig. 3
Phasefolded radial velocity signal of LHS 1140 b (left panels) and LHS 1140 c (right panels). For each planet, upper panel: radial velocity data (colorcoded as in the legend), binned radial velocities (open black symbols) with a bin size corresponding to 15% of the phase, the median model for each planet obtained from the joint fit analysis (see Sect. 6) as the solid black line, and the 68.7 and 95% confidence intervals as shaded dark and light gray regions. Bottom panels: residuals of the median model. 

In the text 
Fig. 4
TPF of LHS 1140 from the TESS observations in Sector 3 (composed with tpfplotter, Aller et al. 2020). The SPOC pipeline aperture is overplotted with shaded red squares, and the Gaia DR2 catalog is also overlaid with symbol sizes proportional to the magnitude contrast with the target, marked with a white cross. 

In the text 
Fig. 5
TESS light curve extracted from the SPOC pipeline and detrended using the PDCSAP flux. Upper panel: complete light curve displaying the transits of the two known transiting planets in the system, LHS 1140 b (red, one transit) and LHS 1140 c (light green, five transits). Lower panel: phasefolded light curve centered on the planet phase for each of the two transiting planets, LHS 1140 b (bottom left) and LHS 1140 c (bottom right). The colored symbols correspond to bins of one fifth of the transit duration (i.e., 30 min for LHS 1140 b and 11 min for LHS 1140c). 

In the text 
Fig. 6
ℓ_{1} periodogram of the LHS 1140 HARPS and ESPRESSO data with the noise model that has the best crossvalidation score. 

In the text 
Fig. 7
Dataset of the FWHM of the crosscorrelation function from the HARPS and ESPRESSO data. The systemic offset for each instrument as retrieved from the twoplanet model has been removed. The median GP model is shown as the solid black line, and the 1σ and 3σ confidence intervals are shown as shaded dark and light gray regions, respectively. Bottom panel: residuals after removing the GP model. 

In the text 
Fig. 8
Phasefolded radial velocity signal of the candidate planet LHS 1140 d in the threeplanet model. Upper panel: radial velocity data (colorcoded as in the legend), binned radial velocities (open black symbols) with a bin size corresponding to 15% of the phase, the median model from the joint threeplanet analysis (see Sect. 6) as the solid black line, and the 68.7 and 95% confidence intervals as shaded dark and light gray regions. Bottom panels: residuals of the median model. 

In the text 
Fig. 9
Difference in the logevidence of the twoplanet and threeplanet models with circular orbits for a cumulative number of data points. Top panel: Bayes factor against the time span of the data set, and bottom panel: against the number of data points. 

In the text 
Fig. 10
Detection limits of the radial velocity dataset for periods shorter than one year. The color code indicates the significance of the radial velocity semiamplitude parameter. The 5σ contour is shown as a solid line, and the limit corresponding to the scatter of the data is shown as a dashed line. The locations of the two known planets are shown as solid circles, and the location of the third planet candidate is shown as an open circle. 

In the text 
Fig. 11
TLS periodogram of the TESS dataset after masking the light curve from the transit times corresponding to the two known planets LHS 1140 b and LHS 1140 c. The red line indicates the strongest power peak, and the blue line indicates the period of LHS 1140 c. 

In the text 
Fig. 12
TESS phasefolded light curve around the Lagrangian points L_{5} (top panel) and L_{4} (middle panel) of LHS 1140 c (see Sect. 5.2). The expected location of the transit is marked by vertical dashed lines. The median transit model inferred from the analysis of the L_{4} dimming is shown as a solid red line (middle panel) and a dashed line in the upper panel to guide the eye on the photometric scatter of the L_{5} region. The 1σ confidence interval is shown in the middle panel as a shaded red region. Bottom panel: residuals of the model for the L_{4} region. All panels show bin data points corresponding to 6.6 min (light blue, ten data points inside the transit duration) and 13.2 min (black, five data points inside the transit duration). 

In the text 
Fig. 13
Massradius diagram for the known exoplanets with the lowest masses with measured mass precisions better than30% (gray symbols). The new locations of the LHS 1140 planets are shown as filled blue (LHS 1140 b) and orange (LHS 1140 c) circles. Their previous locations using the values from Ment et al. (2019) are shown as open squares with dashed error bars and the same color code. The candidate planet at an orbital period of ~ 78 days is shown as the shaded magenta region. Its radius is not known, therefore its vertical location comes from the simple massradiusrelations from Chen & Kipping (2017). The bulk density lines corresponding to different compositions from Zeng et al. (2019) are shown as solid traces, and the dashed lines correspond to isodensities of 1.33, 5.3, and 10 g cm^{−3} (from top tobottom). 

In the text 
Fig. 14
Expected posterior distribution of the midtransit time of planet candidate LHS 1140 d (magenta histograms) and the TESS coverage of Sectors 3 and 30 on the field including the host star (shaded gray region). The median of each posterior distribution is marked by a vertical dotted line, and the corresponding calendar date in Universal Time Coordinate system is annotated. 

In the text 
Fig. 15
Sampled 2D marginal posterior distribution for the CMF and WMF of LHS1140 b. The color code displays the probability density function (PDF). The solid and dashed black lines represent the isoradius curves for the central value of the radius and its 1σ confidence interval limits, respectively. The gray lines delimit the areas we excluded from our sampling (see text). 

In the text 
Fig. 16
Massradius diagram for different compositions considering water in supercritical phase under the surface temperatureand pressure conditions of LHS 1140 c. The black dot and error bars indicate the position of LHS 1140 c in the massradius diagram. 

In the text 
Fig. A.1
Bayesian evidence of different models and datasets (from top to bottom: HARPSonly, ESPRESSOonly, and HARPS+ESPRESSO). Each panel contains 15 different models that include a different number of planets: no planets (GPonly, yellow), oneplanet models (pink), twoplanet models (red), and threeplanet model (green). Each has different orbital configuration assumptions (planets with either circular or eccentric orbits). As an example, model 3p1c2c corresponds to a threeplanet model in which planets 1 (LHS 1140 b) and 2 (LHS 1140 c) have assumed circular orbits, but planet 3 (LHS 1140 d) has free eccentricity. The bars corresponding to each model include this information at the bottom of the bar. At the top of each bar we show the Bayes factor of each model compared to the stronggest evidence model for each dataset. The highest bar corresponds to the strongest evidence model. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.