Open Access
Issue
A&A
Volume 674, June 2023
Article Number A132
Number of page(s) 18
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202346096
Published online 15 June 2023

© The Authors 2023

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

Almost 30 yr after the discovery of the first exoplanets (Wolszczan & Frail 1992; Mayor & Queloz 1995), thousands of planetary systems have been discovered, showing an astonishing diversity of orbits, compositions, and host stellar properties. Understanding the origin of this diversity is one of the major goals of modern astrophysics. Young stars are valuable targets that offer the unique opportunity to provide observational constraints on the timescales involved in the formation and migration of planets in close orbits, the physical evolution of the planets, and atmospheric evaporation due to the strong X-ray and UV irradiation of the host star. Planets have been found in young stellar clusters (e.g. Malavolta et al. 2016; Rizzuto et al. 2017), young stellar kinematic groups (e.g. Benatti et al. 2019; Newton et al. 2019; Nardiello et al. 2022) and single young stars (e.g. Plavchan et al. 2020; Carleo et al. 2021). However, the detection and characterisation of planets around young stars via radial velocity (RV) or transit methods are challenging due to the high-activity levels of the host stars, and some previously claimed planets have been retracted after detailed analysis (e.g. Carleo et al. 2018; Damasso et al. 2020). Furthermore, high rotation rates or variability due to pulsations might also hinder the detection of planets around young stars. As a consequence, the sample of known, well-characterised young exoplanets is still rather limited.

On the other hand, M dwarfs have been recognised as promising targets to search for small, rocky planets with the potential to host life (e.g. Dressing & Charbonneau 2013, 2015). While small planets orbiting around M dwarfs seem to be abundant (e.g Howard et al. 2012; Bonfils et al. 2013; Tuomi et al. 2014), results from RV surveys show that, unlike their FGK counterparts, the frequency of gas-giant planets orbiting around low-mass stars is low (Endl et al. 2006; Butler et al. 2006; Cumming et al. 2008; Pinamonti et al. 2022). This result might be understood in the framework of core-accretion models for planet formation. Within this scenario, the low-mass protoplan-etary disc of an M dwarf conspires against the formation of gas-giant planets. However, using data from the Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015), a small number of transiting gas-giant planets orbiting around low-mass stars have been discovered. Some examples include TOI-1728 b (Kanodia et al. 2020), TOI-442 b (Dreizler et al. 2020), and TOI-3757 b (Kanodia et al. 2022). All these planets have close-in orbits and constitute valuable targets to help us in our understanding of planet formation and migration. While it was proposed early on that close-in gas-giant planets should form beyond the snowline and then experience subsequent inward migration (e.g. Lin et al. 1996), the in situ formation of gas-giant planets has been recently revisited (Batygin et al. 2016; Boley et al. 2016; Maldonado et al. 2018). In addition to planetary companions, M dwarfs can also host brown dwarfs in close orbits. Recent works have discussed differences in the period-eccentricity distribution of massive (MP > 45 MJ) and low-mass brown dwarf (MP < 45 MJ) companions, as well as in the chemical composition of their respective host stars (Ma & Ge 2014; Mata Sánchez et al. 2014; Maldonado & Villaver 2017). It has been shown that massive brown dwarfs have longer periods and larger eccentricities and their host stars do not show the so-called metal-rich signature. According to this scenario, massive brown dwarfs should be mainly formed in the same way as low-mass stars (by the fragmentation of a molecular cloud). On the other hand, low-mass brown dwarfs might be formed at high metallicities by core-accretion or by gravitational instability in low-metal environments (Maldonado & Villaver 2017). Furthermore, since these companions are massive and close to their host stars, tidal interactions between the companion and the host star are expected to constrain the evolution of the system (e.g. Mathieu & Mazeh 1988; Mazeh 2008; Barker 2020; Alvarado-Montes 2022).

In this paper, we report the validation and the confirmation of a transiting low-mass companion (Mp ~ 77 MJ) around the very active M0.5 dwarf TOI-5375 by using space-based photometric observations from TESS combined with high-precision RVs from the ground. Given the high activity levels of TOI-5375, this star has a good chance of being young. TESS identified a candidate with a radius of 1.04 ± 0.03 RJ. A companion with such a radius, if confirmed, could be a giant planet, a brown dwarf, or even a very low-mass star. As discussed, a giant planet transiting an active M dwarf is of special interest. Therefore, we included TOI-5375 in the sample of the Global Architecture of Planetary Systems (GAPS) young program (Carleo et al. 2020) in order to perform the validation and the mass determination of the transiting candidate. Our RV monitoring shows that the companion is actually a massive brown dwarf or a very low-mass star. The system represents a remarkable case study as it is observed in the evolutionary phase in which the star is close to synchronising its rotation with the orbital period of the companion.

This paper is organised as follows. The photometric and spectroscopic data analysed in this paper are presented in Sect. 2. We note that the star is very faint for high spectral resolution observations with 4m class telescope, resulting in low signal-to-noise ratio (S/N) spectra. Therefore, the extraction of RVs in this regime is discussed in this section by comparing different methods. Section 3 reviews the stellar properties of TOI-5375. Section 4 describes the modelling of the photometric and RV variations. The results are discussed at length in Sect. 5. Our conclusions follow in Sect. 6.

2 Observations

2.1 TESS photometry

TOI-5375 was observed photometrically by TESS in short-cadence mode (120 s) in sectors 40, 47, 53, and 60. In this work, we use the light curves extracted from all these datasets. We used the data corrected for time-correlated instrumental signatures, that is, the Pre-search Data Conditioning Simple Aperture Photometry (PDCSAP) flux. The corresponding light curve is shown in Fig. 1. We note that the star was also observed in sectors 20 and 26. However, only High Level Science Product (HLSP) data are available from these sectors. Therefore, unless otherwise noted, data from sectors 20 and 26 are not included in the analysis.

The light curve shows a clear modulation that we attribute to the rotation of the star. A search for periodicities by using the generalised Lomb-Scargle (GLS) periodogram (Zechmeister & Kürster 2009) reveals a clear peak at 1.9692 ± 0.0004 d that corresponds to the rotation period of the star, as shown in Fig. 2. Superimposed on the rotational modulation, the light curve shows periodic transit-like signals, due to a possible substellar companion.

2.2 Ground-based photometry

Two transits of TOI-5375 were observed on October 25 and December 19, 2022, by the TASTE programme, a long-term campaign to monitor transiting planets (Nascimbeni et al. 2011). The AFOSC camera, mounted on the Copernico 1.82 m telescope at the Asiago Astrophysical Observatory in northern Italy, was purposely defocused up to 8″ FWHM to avoid saturation and increase the photometric accuracy. Frames were acquired through a Sloan r′ filter and with a constant exposure time of 15 s. A pre-selected set of suitable reference stars was always imaged in the same field of view to perform differential photometry. Unfortunately, during the first transit, the sky transparency was highly variable through the series and thick clouds passed by during the second half of the transit (egress included).

We reduced the Asiago data with the STARSKY code (Nascimbeni et al. 2013), a software pipeline specifically designed to perform differential photometry on defocused images. After a standard correction through bias and flat-field frames, the size of the circular apertures and the weights assigned to each reference star are automatically chosen by the code to minimise the photometric scatter of our target. The time stamps were uniformly converted to the BJD-TDB standard following the prescription by Eastman et al. (2010). The final (unbinned) light curves are shown in Fig. 13.

thumbnail Fig. 1

Original TESS photometry time series. The energetic events discussed in Sect. 4.6 are shown with different symbols according to their Sector. Within each Sector, different colours indicate that the event occurs at a different period (when phase-folded to the stellar rotation period).

thumbnail Fig. 2

GLS periodogram of the TESS photometry dataset. Values corresponding to a False Alarm Probability (FAP) of 10, 1, and 0.1% are shown with horizontal grey lines. The vertical line indicates the period at 1.9692 d.

2.3 Radiai velocity observations

A total of 9 HARPS-N (Cosentino et al. 2012) observations of TOI-5375 were collected in the period between BJD = 2 459 704 d and BJD = 2 459 718 d (May 4-18,2022). HARPS-N spectra cover the wavelength range 383–693 nm with a resolving power of R ~ 115 000. The observations were performed within the framework of the GAPS project (Covino et al. 2013). Data were reduced using the HARPS-N Data Reduction Software PRS; Lovis & Pepe 2007; Dumusque et al. 2021) that takes into account the usual procedures involved in échelle spectra reduction (bias level, flat-fielding, order extraction and wavelength calibration, as well as merging of individual orders). Radial velocities are computed by cross-correlating the spectra of the target star with an optimised binary mask (Baranne et al. 1996; Pepe et al. 2002). However, this procedure might not be optimal for a star like TOI-5375. To start with, TOI-5375 is an M dwarf-type star and the spectra of M dwarfs suffer from heavy blends. Side-lobes appear in the cross-correlation function (CCF) that might affect the RV precision (e.g. Rainer et al. 2020). In addition, TOI-5375 is a very faint star (V ~ 14 mag). This might have an effect on the selection of an appropriate mask to perform the CCF. In the M2 mask, all lines have roughly a similar weight, so weak (not deep enough) lines are lost in the not-well-defined continuum that characterises the spectra of an M dwarf, adding noise to the CCF profile. In contrast, if the K5 mask is used, only deep lines are considered and the CCF profile becomes less noisy. Furthermore, high levels of stellar activity might deform the core of the average line profile, while stellar rotation broadens the lines. Therefore, special care should be taken in selecting the half-window of the CCF in order to include the continuum when fitting the CCF profile (see e.g. Damasso et al. 2020). In addition, since we are dealing with a very red object, the signal in the bluest part of the spectra is very low, so the analysis should be limited to the red part of the spectrum. In order to illustrate these effects, we compare the RVs obtained with the DRS by using the K5 and M2 masks as well as the K5 mask with a half-window of 60 km s−1 (instead of the default value of 20 km s−1). The results from this comparison are given in Table 1, while in Fig. 3 we show the corresponding CCF profiles for one of our observations (we chose the one with the highest S/N, ~13 at spectral order 46, 548–554 nm). It is clear from the figure that the CCF profile obtained with the M2 mask comprises only noise, without any visible signal. The CCF profile obtained with the original K5 mask, although asymmetric has a Gaussian-sharp form but lacks a good continuum for a proper fitting. The enlarged-window K5 mask, on the other hand, provides a better fit to a Gaussian function with a lower error. These calculations have been done by using the DRS version implemented through the YABI interface (Hunter et al. 2012) at the Italian Centre for Astronomical Archives1.

RVs were also computed with the Java-based Template-Enhanced Radial velocity Re-analysis Application (TERRA; Anglada-Escudé & Butler 2012). TERRA measures the RVs by a least-square match of the observed spectrum to co-added high S/N template spectra derived from the same observations. The bluest part of the spectra was excluded from the analysis and only orders with λ > 434 nm were considered. It has been shown that for M dwarf stars the RVs derived by TERRA show a lower root-mean-square (rms) than the DRS and should be preferred (Perger et al. 2017). However, we note that in the case of TOI-5375 the rms of the RVs obtained by using TERRA is ~34 times lower than the rms obtained when using the RVs derived by the DRS. While the CCF RVs show a rms between 9.9 and 17 km s−1 (depending on the mask considered), TERRA RVs show a rms of only 0.5 km s−1. TERRA RVs are not consistent with the RVs variations seen in the spectra. A quick look shows that spectral lines show relative shifts between different nights of the order of several km s−1. A possible explanation for this discrepancy might be related to the low S/N of our spectra. It is likely that TERRA identifies some of the noise as ‘true’ spectral features. In order to check this fact, we made use of the SpEctrum Radial Velocity AnaLyser tool (SERVAL; Zechmeister et al. 2018). The concept behind SERVAL is similar to TERRA, but it constitutes an alternative implementation. Again, we derive RVs that show a very low rms (~1 km s−1). Therefore, we decided to exclude the TERRA and SERVAL RVs from our analysis.

In order to avoid the aforementioned difficulties we decided to derive RVs using a self-written python package that computes the mean line profile using the least-squares deconvolution technique (LSD; Donati et al. 1997) with a stellar mask obtained from the VALD32 database (Piskunov et al. 1995; Ryabchikova et al. 2015) with appropriate stellar parameters for TOI-5375. The mean line profiles are then fitted with a rotational broadening function. Further details can be found in Rainer et al. (2021). An example of one of the derived LSD profiles is shown in Fig. 3 where it can be compared with the CCF profiles obtained by the DRS. The LSD derived RVs have a rms of 12.42 km s−1 and a mean uncertainty of 0.36 km s−1. Figure 4 shows a comparison between the RVs velocities obtained by using different masks and methods. We note that the RVs derived with the enlarged-window K5 mask are similar to the ones derived with the LSD technique. This shows that by performing an appropriate mask and half-window selection it is possible to derived precise RVs using the CCF method even in spectra at low S/N or for the case of fast rotators. As mentioned, RVs based on the template-matching algorithm (TERRA and SERVAL) provide rather flat RV curves when compared with the CCF and LSD values. We note that the CCF and LSD techniques combine the spectral lines to create a mean profile, and thus it should be possible to achieve a reasonable S/N. In fact, the S/N of the CCF increases roughly as N$\sqrt N $ where N is the number of stellar lines used to compute the CCF (typically of the order of several thousands of lines). On the other hand, it seems that the template-matching algorithm, when all the observed spectra have a low S/N, is not able to generate a suitable template to compute accurate RVs. Finally, RVs were also computed from the LSD profiles by using the SpotCCF software (Di Maio et al., in prep.), which takes into account the presence of multiple spots on the stellar surface. It also accounts for the wings of the LSD profile by convolving the rotational profile with a Lorentzian profile. The code was run without considering any spot. We find that the SpotCCF-derived RVs are quite similar to the ones derived from the LSD analysis, thus confirming the reliability of our measurements. The analysis presented in the following refers to the LSD-derived RVs.

Table 1

Comparison between the different methods used to extract the RVs.

thumbnail Fig. 3

Examples of standardised (mean subtracted and division by the standard deviation) CCF profiles obtained by the DRS using different binary masks as well as the LSD profile.

thumbnail Fig. 4

RVs obtained using different masks and methods versus LSD RVs. For TERRA and SERVAL RVs, we subtracted the mean LSD RV for the comparison. Each panel shows the mean value and standard deviation of the difference RVs other methods - LSD RVs. Grey lines show the one-to-one relationship.

3 The host star

TOI-5375 (a,k.a. TIC 71268730, 2MASS J07350822+7124020, 1RXS J073508.9+712415) is a M0.5 dwarf located at a distance of 121.65 ± 0.17 pc from the Sun (Gaia Collaboration 2022). Its main stellar properties are listed in Table 2. Stellar effective temperature, spectral-type and metallicity are determined from a methodology based on the use of a principal component analysis and sparse Bayesian fitting methods3, which make use of the same spectra used in this work for the RV analysis (see Sect. 2). In brief, a relationship between the principal components and the stellar parameters is derived by training the data with the spectra of M dwarfs orbiting around an FGK primary or M dwarfs with effective temperature obtained from interfero-metric measurements of their radii. Further details can be found in Maldonado et al. (2020). We also compute Teff and [Fe/H] by using photometric relationships (Casagrande et al. 2006; Neves et al. 2012). We note that the photometric Teff is significantly lower than the spectroscopic one. This discrepancy has already been noted and discussed (see e.g. Maldonado et al. 2015, and references therein). It is related to the fact that M dwarfs emission deviates from a blackbody for wavelengths beyond ~2 µm. The photometric estimate of [Fe/H] is also lower than the spec-troscopic one. In a recent work, Maldonado et al. (2020) show that the use of photometric estimates for M dwarfs gives a [Fe/H] distribution that is shifted towards lower values when compared with an otherwise similar sample of FGK stars. Therefore, spec-troscopic parameters are preferred, although we caution against the relatively low S/N of our spectra (of the order of 44 in the 604–605 nm region, measured in the coadded spectra). Stellar mass is computed following the relations based on near-infrared photometry by Henry & McCarthy (1993), while the stellar radius is obtained using the mass-radius relationship provided by Maldonado et al. (2015). From the stellar mass and radius values, the surface gravity was estimated. The stellar luminosity is computed by using the bolometric corrections provided by Flower (1996).

The rotational period is computed from the analysis of the TESS light curve (see Sect. 2) while the equatorial velocity is estimated from the rotation period and the stellar radius. The available photometry for TOI-5375 has been used to build and trace the star spectral energy distribution (see Fig. 5) along with PHOENIX/BT-SETTL atmospheric models (Allard et al. 2012).

Table 2

Physical properties of TOI-5375.

thumbnail Fig. 5

Spectral energy distribution of TOI-5375. A PHOENIX/BT-SETTL atmospheric model (Teff = 3885 K, log ɡ = 4.66 cm s−2, [Fe/H] = +0.26 dex) is overplotted. Errors in the observed fluxes are within the symbol’s size.

Age estimates

Stellar age is one of the most difficult parameters to determine in an accurate way, even for main-sequence solar-type stars, as they evolve too slowly to be dated by their position in the Hertzsprung-Russell Diagram. For solar-type stars, it is well known that stellar activity and rotation are related by the stellar dynamo (e.g. Noyes et al. 1984) and both activity and rotation diminish as the star evolves and loses angular momentum. Thus, activity and rotation tracers, such as Ca emission, X-ray, or rotation periods, are often used to estimate stellar ages (see e.g. Mamajek & Hillenbrand 2008). However, in this case, the rotation of TOI-5375 is being accelerated by the tides raised by the companion (as we discuss in Sect. 5). As a consequence, the star can show a high level of chromospheric or X-ray emission despite being rather old. Therefore, age diagnostics based on the rotation period will only provide lower limits to the age of the TOI-5375 system. In the following, we consider as ‘young’ those stars with an age younger than the Hyades, ~750 Myr.

To compute the LX of TOI-5375 we search for an X-ray counterpart in the ROSAT All-Sky Bright Source Catalogue (Voges et al. 1999) and apply the count rate-to-energy flux conversion factor relation given by Fleming et al. (1995). The derived X-ray luminosity is log LX = 29.44 ± 0.22 erg s−1, while the position of TOI-5375 in the log(LX/LBol) versus colour index (B – V) plane is shown in Fig. 6 (left). Bolometric corrections are derived from the (B – V) colour by interpolating in Flower (1996, Table 3). Data from the Pleiades, age ~112 Myr (Stauffer et al. 1994) and Hyades, age ~750 Myr (Stern et al. 1995) are overplotted for comparison. The position of TOI-5375 in this diagram is consistent with the position of the Pleiades.

Figure 6 (centre) shows the rotation period of TOI-5375 as a function of the colour index (BV). For comparison, data from the Pleiades (Prosser et al. 1995) and from the Hyades (Radick et al. 1987) are overplotted. Two gyrochrones (at the ages of the Pleiades and the Hyades) are also overplotted. The figure seems to indicate that TOI-5375 is a young star with a fast rotation period4.

TOI-5375 shows a clear emission in the Ca II H & К lines. Balmer lines (from Hα to Hє), as well as other known activity indexes such as the Na I D1, D2 doublet and the He I D3 line are also shown in emission. To quantify the stellar chro-mospheric activity we make use of the log(RHK)$\left( {{{R'}_{{\rm{HK}}}}} \right)$, defined as the ratio of the chromospheric emission in the cores of the broad Ca II H & К absorption lines to the total bolometric emission (e.g. Noyes et al. 1984). We start by computing the so-called S index, which is measured in the Mount Wilson scale by the HARPS-N DRS implemented through the YABI interface. The S index contains both the contribution of the photosphere and the chromosphere. In order to account for the photospheric correction we used the calibration based on the colour index (BV) developed by Suárez Mascareño et al. (2015) who updated and extended the original calibration from Noyes et al. (1984) to the low-mass regime up to a (B – V) value of 1.90. The derived log(RHK)$\left( {{{R'}_{{\rm{HK}}}}} \right)$ for TOI-5375 is −3.06 ± 0.10. This value can be compared with the activity level of stars in clusters with well-known ages. The comparison is shown in Fig. 6 (right). Data from the Sco-Cen complex (age ~ 5 Myr), Pleiades, Hyades, and M67 (age ~ 4000 Myr) are from Mamajek & Hillenbrand (2008). TOI-5375 is located in the range of very active stars with a mean log(RHK)$\left( {{{R'}_{{\rm{HK}}}}} \right)$ higher than other stars in young stellar associations and clusters.

As an additional age diagnostic we plot the position of TOI-5375 in a colour-magnitude diagram (see Fig. 7). The ZAMS position as well as several pre-main sequence isochrones (computed at metallicity +0.20 dex) from Siess et al. (2000) are plotted. We use the Mv versus (B–V) (left panel) and Mv versus (VKS)5 plots (right panel), as it has been shown that young rotating stars in the Pleiades tend to show blue excesses, that is, bluer (B–V) colours, than other relatively old stars (Stauffer et al. 2003). Indeed, we note that the position of the star in the diagram, when considering the (B – V) colour, suggests an age between 30 and 200 Myr (we note the large uncertainty in В magnitude), whilst when considering the (VKS) colour, the position of TOI-5375 is more consistent with a younger, ∼30 Myr age.

Membership in stellar associations and young moving groups has been used as a methodology to identify young stars and to assign ages, especially after the release of the HIPPARCOS data (see e.g. Maldonado et al. 2010, and references therein). The Galactic spatial-velocity components (U, V, W) of TOI-5375 are computed from the stellar RV from our own orbital solution, together with Gaia parallaxes and proper motions (Gaia Collaboration 2022). We follow the procedure described in Montes et al. (2001) and Maldonado et al. (2010), where the original algorithm (Johnson & Soderblom 1987) is adapted to epoch J2000 in the International Celestial Reference System (ICRS) following Sect. 1.5 of the HIPPARCOS and Tycho catalogues (ESA 1997). We made use of the full covariance matrix to compute the uncertainties. In this way, we account for the possible correlation between the astrometric parameters. The corresponding (U, V, W) values are provided in Table 2, while the position of TOI-5375 in the (U, V) plane is shown in Fig. 8. From the analysis of its velocity components we see that TOI-5375 is located outside the region of the young stars so we conclude that it is a field thin disc star. This conclusion is supported by the use of Bayesian methods (BANYAN; Gagné et al. 2018)6 that confirm that TOI-5375 is, with high-probability (99.9%), a field star. The high-metallicity content of the star might be an additional indication that the star is not young. This is because typically, the metallicity of young stars in the solar neighbourhood is close-to-solar or slightly sub-solar (Spina et al. 2017).

To summarise, age estimates based on the rotation of the star seem to indicate that TOI-5375 is a young star. However, other age estimates based on the stellar metal content or the kinematics do not support this statement. The disagreement is due to the fact the stellar rotation period is being accelerated by the companion, as we shall discuss in Sect. 5. Therefore, the only conclusion that we can draw from this analysis is that TOI-5375 is a field star that should have an age older than the one of stars in young stellar associations like the Pleiades or the Hyades.

thumbnail Fig. 6

Age diagnostics for TOI-5375. Left: fractional X-ray luminosity log(Lx/LBol) versus colour index (BV). Centre: rotation period as a function of the (BV) colour. Two gyrochrones (at the ages of the Pleiades and the Hyades) have been overplotted for comparison. Right: log(RHK)$\left( {{{R'}_{{\rm{HK}}}}} \right)$ versus colour index (BV). Dashed lines are the limits for very active, active, inactive, and very inactive stars, according to Henry et al. (1996).

thumbnail Fig. 7

Colour-magnitude diagram. Left: Mv as a function of (BV) colour. Right: Mv as a function of (VKs) colour.

thumbnail Fig. 8

(U, V) plane showing the position of the young disc population (dashed line) as defined by Eggen (1984, 1989) as well as the known kinematic groups (large crosses and asterisks) in the solar neighbourhood (see Maldonado et al. 2010; Maldonado 2012, and references therein). The dotted line represents the velocity ellipsoid determined by Francis & Anderson (2009). The position of TOI-5375 is shown with a black star (uncertainties are within the symbol size). All velocities are in the Solar System reference frame.

4 Analysis

4.1 Validation and vetting tests

As mentioned, the TESS photometric light curve of TOI-5375 shows clear periodic dips compatible with the transits of a substellar companion orbiting around this active, M dwarf. Because of the low spatial resolution of the TESS cameras (≈21 arcsecpixel−1) some of the objects initially identified as a substellar companion might indeed be a false positive (FP). We, therefore, perform an analysis to identify possible nearby contaminating stars that might be the origin of blended eclipsing binaries (BEB). The analysis is based on the use of the photometric data provided by Gaia and is described at length in Mantovan et al. (2022).

The analysis shows that, besides TOI-5375, there is a star at 33.5 arcsec (Gaia DR2 1110586939683894528) separation from TOI-5375 that might reproduce the transit signal (see Fig. 9). To be sure that the transits are not due to contamination by this neighbour, we perform the analysis of in- or out-of-transit centroid check following the procedure described in Nardiello et al. (2020) and Nardiello (2020). Data from sectors 20 and 26 are included in this analysis. It should be noticed that the in- or out-of-transit centroid analysis is performed on the full frame images and therefore it is independent from the analysis of the light curve. The results are shown in Fig. 10. It can be seen that, within the errors, the in- or -out-of-transit mean centroids are in agreement with the position of TOI-5375.

Finally, in order to rule out the possibility that the transiting candidate is an FP we use the VESPA7 software (Morton 2012, 2015). We obtain a 100% probability of having a Keplerian transiting companion around TOI-5375, while the probability of a FP is of the order of ~7 × 10−10. All these analyses confirm that TOI-5375 b is fully and statistically vetted and only requires reconnaissance spectroscopy to be promoted to a ‘statistically validated’ companion.

thumbnail Fig. 9

Gaia stars identified in the TESS field. The position of TOI-5375 is shown with the letter ‘ť.

thumbnail Fig. 10

In- or out-oſ-transit centroid tests calculated for each sector. Black dots are background stars.

4.2 Photometric time series analysis

In order to characterise the properties of TOI-5375 b, the photometric data were modelled by using a Bayesian framework based on a Monte Carlo sampling of the parameter space The photometric light curve contains the contribution of the proposed planetary signal as well as the effect of the stellar activity of the host star. The stellar signal is modelled with a Gaussian process (GP) analysis The covariance (or kernel) function adopted is a simple harmonic oscillator (SHO) kernel that is defined as follows k(i,j)=σGPexp(ΔtτGP)  [ cos  (qProtΔt)+ProtqτGPsin  (qProtΔt) ]           q=2π(2πτGPProt)1,$\matrix{ {k\left( {i,j} \right) = {\sigma _{{\rm{GP}}}}\,\exp \,\left( {{{ - {\rm{\Delta }}t} \over {{\tau _{{\rm{GP}}}}}}} \right)\,\,\left[ {\cos \,\,\left( {{q \over {{P_{{\rm{rot}}}}}}{\rm{\Delta }}t} \right) + {{{P_{{\rm{rot}}}}} \over {q{\tau _{{\rm{GP}}}}}}\sin \,\,\left( {{q \over {{P_{{\rm{rot}}}}}}{\rm{\Delta }}t} \right)} \right]} \hfill \cr {\,\,\,\,\,\,\,\,\,\,\,q = 2\pi \,\sqrt {\left( {2\pi {{{\tau _{{\rm{GP}}}}} \over {{P_{{\rm{rot}}}}}}} \right) - 1} ,} \hfill \cr } $

where k(i, j) is the ij element of the covariance matrix, Δt = titj, where ti and tj are two times from the photometric dataset, σGP is the amplitude of the covariance, Prot is the rotation period, and τGP is the timescale of the exponential decay. In order to obtain the parameters corresponding to the best transit model we use the batman package (Kreidberg 2015) with zero eccentricity and a linear limb darkening law. This hypothesis is justifiable thanks to the strong tidal interaction that tends to quickly damp any initial eccentricity. Of course, this is assuming that there are no other bodies in the system that could re-excite the eccentricity. Our full model also contains an extra white noise (jitter) term (σj) per each observed sector. Before fitting, each sector data were normalised by their median value. The (hyper)-parameter space was sampled by using the publicly available nested sampler and Bayesian inference tool MULTINEST V3.1Ø, (e.g. Feroz et al. 2019) through the PYMULTINEST (Buchner et al. 2014) wrapper, set to run with 5000 live points and a sampling efficiency 0.3. Priors of the model are provided in Table 3.

The resulting posterior distributions are shown in Fig. 11 while Fig. 12 shows the best transit model fits all transits folded in the light curve after the subtraction of the stellar activity signal (top) as well as the residuals of the light curve once the companion transit model is subtracted (bottom). The best-fit parameters (median values) are listed in Table 3. As seen, the results are consistent with the presence of a companion orbiting around TOI-5375 with a period of 1.72 d. This value is close to, but significantly different from, the rotational period of the star. Therefore, it is likely that the star is close to synchronising its rotation with the orbital period of the companion (further discussion will be provided in Sect. 5). This justifies our choice of a circular orbit.

We performed a search for additional companions by performing a GLS periodogram of the best-fit residuals8. The rms of the residuals is very low and no significant additional signals were found. The ‘best’ period is found at ~ 1.97 d (that is, at the stellar rotation period) but with a Scargle power of ~3 × 10−6, which implies a FAP of 100% (we note that the analytical FAP level of 10% corresponds to a power of 0.0014). This result confirms that our treatment of the stellar activity by means of a SHO kernel does effectively remove it from the light curve. It also shows that if there were more companions in the system, they would be difficult to reveal (at least with the data at hand).

Table 3

Best-fit values obtained for the star-planet model.

thumbnail Fig. 11

Posterior distributions of the Bayesian fit of the photometric time series.

thumbnail Fig. 12

Photometric modelling of TOI-5375. Top: folded transits in the light curve after the subtraction of the stellar activity signal. The green line shows the best model fit of the transits. Bottom: residuals of the light curve after the subtraction of the companion transit model.

4.3 Ground-based photometry analysis

We also fit the ground-based photometry dataset obtained from the Asiago observatory. These observations confirm the transit of a companion around TOI-5375 and are additional proof that TOI-5375 b is real and not a false positive due to an BEB. We followed a similar analysis as described in Sect. 4.4 (see below), but we did not consider the dilution factor (which defines the total flux that falls into the photometric aperture from contaminants divided by the flux contribution of the target star) and a second-order polynomial trend. We treated the stellar limb-darkening (LD) contribution by setting the Sloan r′ filter in PyLDTk. From the analysis of the ground-based photometry, we derive a companion-to-stellar-radius ratio, Rb/R, of 0.1850 ± 0.0012, which is larger than the value derived from the analysis of the TESS dataset (see Fig. 13). This difference can be explained by several reasons. To start with, this effect can be due to the fact that the TESS and ground-based observations are taken at different wavelengths (atmospheres absorb larger at some wavelengths and less at others).

Other explanations can be due to stellar activity. A visual inspection of the individual TESS transits reveals that some of them show ‘bumps’ during the transit that can be related to the presence of stellar spots. In some cases, there are also negative bumps, likely due to the occultation of faculae. We also note that the TESS and the Asiago datasets are not simultaneous. To clarify this point we analyse each individual TESS transit as well as the ones observed from the Asiago Observatory. In order to analyse the space and the ground-based photometry in a homogeneous way, we use the TESS SAP light curve using a dilution factor to take into account the possible contamination by nearby stars. In addition, the transit depth is fitted simultaneously to both TESS and Asiago datasets.

Figure 14 shows the histogram of the derived Rb/R for each individual transit. Excluding one clear outlier, the values of Rb/R obtained from the analysis of the individual TESS transits show a rms of 0.053 covering the range between 0.151 and 0.172. The apparent variation in the companion radius due to the effects of starspots can be written as: ΔRR12Aλ(fif0),${{{\rm{\Delta }}R} \over R} \simeq - {1 \over 2}{A_\lambda }\,\left( {{f_i} - {f_0}} \right),$(1)

where Aλ is the contrast, which is a function of the bandpass, ƒ0 is the filling factor due to un-occulted spots, and ƒi is the filling factor of spots occulted during the transit. Following Ballerini et al. (2012), for an early-M dwarf like TOI-5375, and taking into account that the TESS passband covers the range from 600 to 1000 nm and is centred on the Cousins I-band, we consider Aλ = 0.924. Then, in the extreme case of ƒi = 0, we derive a value of ƒ0 of the order of 0.285, which is a reasonable value for a very active star like this case. Therefore, we conclude that the spread in Rb/R from one transit to another can be explained by the stellar activity of the host star. On the other hand, the difference in transit depths between the TESS and the Asiago observations is likely due to a combination of activity and the effect of observing at two different bandpass.

Another possibility could be due to an additional companion (of redder colour) spatially unresolved in Gaia observations. It should be noted that very close binaries (P < 3 days) have very often tertiary companions (Tokovinin et al. 2006). Also, an unresolved companion might explain the slight offset in (V-K) colour with respect to the main sequence. However, this hypothesis is quite speculative, and we do not have indications of an additional component in our spectra.

4.4 Independent joint analysis: TESS and ground-based photometry

In order to have an independent check of our results we perform an additional test by analysing the TESS light curve using the PyORBIT software9 (Malavolta 2016; Malavolta et al. 2018) finding similar results. In this case, we analyse the ground-based photometry simultaneously with the TESS data. Furthermore, to test whether or not the difference in transit depth between the Asiago and the TESS data could be due to a possible contamination of the TESS aperture, we used the Simple Aperture Photometry (SAP) fluxes over Presearch Data Conditioning (PDC) ones. In this way, we conserve the contribution of stellar contamination from neighbour stars and later perform our correction for stellar dilution. We followed Mantovan et al. (2022) to calculate the fraction of light of the stellar PSF which falls into a given photometric aperture (k parameter) for each star, and we then got the dilution factor (and its associated error).

We selected from the SAP light curves each transit event and an out-of-transit portion twice as long as the transit duration (almost two hours, both before the ingress and after the egress). To perform the selections, we created a mask that flags each transit and cuts portions of the light curves accordingly. For each TESS transit, we fitted the following: the central time of transit (T0), the planetary-to-stellar radius ratio, the impact parameter, b, the stellar density (ρ, in solar units), the quadratic LD law with the parameters q1 and q2 introduced in Kipping (2013) – where u1 and u2 have been estimated using PyLDTk10 (Husser et al. 2013; Parviainen et al. 2015) and considering the TESS filter, the dilution factor, a second-order polynomial trend (with c0 as the intercept, c1 as the linear coefficient, and c2 as the quadratic coefficient), and a jitter term (different for each sector, to take into account possible TESS systematics and short-term stellar activity noise) to be added in quadrature to the errors of TESS light curves. We imposed a Gaussian prior on the stellar density, LD coefficients, and dilution factor, while we kept fixed the period, eccentricity, and argument of pericentre.

We performed a global optimisation of the parameters by running a differential evolution algorithm, pyDE (Storn & Price 1997) and then a Bayesian analysis of each selected light curve around each transit. To perform the latter, we used the affme-invariant ensemble sampler (Goodman & Weare 2010) for Markov chain Monte Carlo (MCMC) implemented within the emcee package (Foreman-Mackey et al. 2013). We modelled each transit with batman (Kreidberg 2015). We fit each transit with 4ndim walkers (with ndim being the model's dimensionality) for 20 000 generations with pyDE and then with 300 000 steps with emcee – where we applied a thinning factor of 200 to reduce the effect of the chain auto-correlation. We discarded the first 25 000 steps (burn-in) after checking the convergence of the chains with the Gelman–Rubin (GR) statistics (Gelman & Rubin 1992, threshold value Ř = 1.01). From this analysis, we derive a companion-to-stellar-radius ratio, Rb/R, of 0.1699 ± 0.0011, similar to the value obtained in the previous analysis (Rb/R, of 0.165 ± 0.002). The small difference is likely due to the fact that this time, we analyse the ground-based photometry simultaneously with the TESS data. The best-fitting parameters and the corresponding plots are shown in Appendix A.

thumbnail Fig. 13

Photometric observations from Asiago, folded on the transit model fit period of TOI-5375 b. Left: October 25, 2022, observations. Right: December 19, 2022, observations. The green line shows the best fit to the Asiago dataset, while for comparison purposes we also show in red the best fit derived for the TESS light curves. The residuals of the fit to the Asiago ground-based photometry are shown in the bottom panel.

thumbnail Fig. 14

Histogram of the derived Rb/R values from each individual TESS and Asiago transit.

4.5 Search for TTVs and other variability effects

Finally, even though the residuals of the fit of the TESS photometry do not reveal additional signals, a search for transit timing variations (TTVs) due to a possible additional companion orbiting around TOI-5375 was performed (see e.g. Borsato et al. 2019). To this purpose we analyse each individual TESS transit as well as the second transit observed from the Asiago Observatory (the first transit is not included, since, as we have seen, its sampling is not complete). Once the transit times were derived, we fit them with a straight line to obtain a linear ephemeris: T0=2459604.837650(±3×107)${T_0} = 2459604.837650\left( { \pm 3 \times 10 - 7} \right)$(2) +N×1.721557(±2×106)BJDTDB,$ + N \times 1.721557\left( { \pm 2 \times {{10}^{ - 6}}} \right)\,{\rm{BJ}}{{\rm{D}}_{{\rm{TDB}}}},$(3)

where N is the epoch of the transit. The corresponding fit is shown in Fig. 15. A search for periodicities in the residuals was performed by means of a GLS periodogram. The periodogram identifies a significant peak at 12.72 ± 0.03 d with a FAP lower than 1%. The semi-amplitude of the signal is ~1 min. As said, some of the TESS transits are affected by stellar activity. Indeed, according to the simulations performed by Oshagh et al. (2013, Fig. 8), for a transit with Rb/R ~ 0.15, TTVs of the order of one minute can be explained by the presence of stellar spots with a filling factor of ~0.6%. Given the very active nature of TOI-5375, our conclusion is that the ~ 1 min amplitude TTV signal seen in this star is most likely related to stellar activity. We note that the synodic period between the stellar rotation and the orbital period is 13.59 days, which is close to the apparent TTV period of 12.7 days. This reinforces the interpretation in terms of surface activity features of the apparent TTVs.

Notes.(‡) 0: Accept null hypothesis; 1: Reject null hypothesis.

Finally, we estimate the photometric flux variations due to other variability effects. In particular, we consider the Doppler beaming (the reflex motion of the star induced by its planetary companion produces Doppler variations in its photometric flux), the ellipsoidal light variations due to tidal effects on the star from the companion as well as the reflected light in the companion surface. Following Loeb & Gaudi (2003) we derive flux variations, ΔF/F0, of the order of 2.8 × 10−4, 9.4 × 10−5, and 2.4 × 10−4 for Doppler beaming, ellipsodial light variations, and reflected light, respectively. The calculation of Doppler beaming was done in the TESS band (λ ~ 7850 Å), while for the reflected light effect a geometric albedo of 2/3 was considered. All these effects are extremely small compared with the intrinsic variability of the star (the TESS light curve shows a rms of 26.7, 54.35, 68.68, and 43.96 e s−1 in Sectors 40, 47, 53 and 60, respectively).

Table 4

Results of the K-S tests performed in this work.

thumbnail Fig. 15

Residuals of the timings of mid-transits versus a linear ephemeris.

4.6 High energy events in the TESS light curve

Many energetic events or flares seem to be present in the TESS light curve (see Fig. 1). This is in agreement with the high rotation level of this star but it might also be evidence of a possible interaction between TOI-5375 and its companion. In order to test this hypothesis we use the methodology described by Colombo et al. (2022) to identify the impulsive events and derive their amplitudes and the amount of energy emitted. We then phase-fold the times at which the events take place with the companion period (1.72 d), the rotation period (1.97 d) as well as the synodic period of the TOI-5375 system. The synodic period is defined as 1/psynodic = 1/Pcompanion 1/Protation, that is, Psynodic =13.59 d.

Figure 16 shows the number of events (blue empty histogram) and the released energy (orange filled histogram) phase-folded with the companion period (left) and with the synodic period (centre), and the stellar period (right). In order to test whether the distribution of the released energy is or not uniform we made use of a Kolmogorov-Smirnov (KS) test11. The results are shown in Table 4. In all cases, we test the null hypothesis that the phase-folded released energy distribution is drawn from a uniform distribution. The results from the K-S test show that the distribution of events is far from uniform, which supports the star-companion interaction scenario. When phased to the stellar rotation period, four main peaks of events can be seen at phases ~0.1, 0.4, 0.7 and ~0.9. The first peak is driven by an energetic event in Sector 53 (see Fig. 1, green square). The peak at phase ~0.4 corresponds to an energetic event in Sector 47 (shown as a purple triangle in Fig. 1), which is followed by many more, less energetic, events. The peak at phase ~0.7 is driven by two high-energetic events in Sectors 60 and 53 (marked with a red diamond and a red circle in Fig. 1), again followed by several less energetic events. Finally, it can be seen a fourth peak close to a phase 0.9 corresponding to a high-energy event in Sector 40 (see Fig. 1, red circle). On the other hand, if the distribution of events is folded with respect to the period of the orbital companion (Fig. 16, left) or the synodic period (Fig. 16, centre), the impulsive events are mainly at phase values between 0.0 and 0.6 and are driven by the five energetic events already discussed.

4.7 RV time series analysis

Table 5 provides the RVs time series for TOI-5375 analysed in this work, which, as discussed earlier, are those obtained by the LSD method. The RV data show a rms of 12.42 km s−1, which is approximately 35 times the mean error of the measured RVs (0.36 km s−1). The GLS periodogram (see Fig. 17) reveals a double peak structure with significant frequencies at f1 = 0.4202 ± 0.0008 d−1 (P = 2.380 ± 0.004 d), and f2 = 0.5810 ± 0.0006 d−1 (P = 1.721 ± 0.002 d). For frequencies values larger than 0.5 d−1 the periodogram is clearly a mirror copy of the one between 0 and 0.5 d−1 so both peaks should be related. The significance of these frequencies was tested by using a bootstrapping analysis to compute the FAPs values (see e.g. Endl et al. 2001). The signal at 1.72 d agrees well with the periodicity found for the transit of TOI-5375 b in the light curve analysis. On the other hand, the period at 2.38 d is likely due to the observing sampling as our observations have a gap of ~6 d (∣f1 - f2∣ ~ 1/6 d−1).

We applied the same Bayesian framework as for the photometric curve to the RV data but, in this case, the model contains only a Keplerian signal and the number of live points has been set to 2000. Given the small number of RVs points, we do not perform a simultaneous analysis of the RV with the photometry.

Instead, we take advantage of the photometry fit and set tight priors in the RV fitting on all parameters but in Kb, where a wide prior was considered (see Table 3). The resulting posterior distributions are shown in Fig. 18 while Fig. 19 shows the best-fit planet model (top) and the corresponding RV residuals (bottom). The best-fit parameters are also listed in Table 3. No significant signals are found in the analysis of the RV residuals. We note that the resiudals show a larger dispersion during the first observing season. A possible explanation for this could be related to the high variability and activity level of TOI-5375.

For the sake of completeness, we also measured the main activity indexes of chromospheric activity, namely the Ca II H & К index, Hα, and Na I D1, D2 (see e.g. Maldonado et al. 2019). All indexes show a symmetric periodogram around 0.5 d−1 with several peaks, all with a significance below the 10% analytical FAP value. Only the Hα index shows a dominant peak, close to the rotation period at ~ 1.99 d. The time series of the different activity indexes as well as their corresponding periodograms are shown in Appendix B.

thumbnail Fig. 16

Histogram of impulsive events (blue empty line) and amount of energy emitted (orange filled histogram) phase-folded to the companion orbital period (left), the synodic period (centre) and the stellar period (right).

thumbnail Fig. 17

GLS periodogram of the RV dataset. Values corresponding to a FAP of 10, 1, and 0.1% are shown with horizontal grey lines. The vertical line indicates the period of TOI-5375 b at 1.72 d.

Table 5

RVs (with their corresponding uncertainties) of TOI-5375 used in this work.

thumbnail Fig. 18

Posterior distributions of the Bayesian fit of the RVs.

thumbnail Fig. 19

Best-fit planetary model for TOI-5375 b (top), and RV residuals (bottom).

thumbnail Fig. 20

Known substellar companions and low-mass stars (radius versus mass) around late-K and M dwarfs (as listed at http://exoplanet.eu/). Different colours and symbol sizes indicate the effective temperature of the host star. Models for substellar objects and low-mass stars with solar metallicity from Baraffe et al. (2015) are over-plotted with grey lines.

5 Discussion

From the best values of Kb and pb derived in the previous section, we derived a mass for TOI-5375 b of 77 ± 8 MJ and a semi-major axis a = 0.0251 ± 0.0012 au. On the other hand, from the analysis of the light curve, we obtain a radius of 0.99 ± 0.16 RJ. TOI-5375 b is therefore an object in the boundary between the brown dwarfs and the very-low-mass stars. This is a remarkable result, as massive substellar companions orbiting around M dwarfs are rare (see Sect. 1). Figure 20 shows the position of TOI-5375 b in the planetary radius versus planetary mass diagram. For comparison purposes, the location of the known substellar companions around late-K and M dwarf stars (with known values of planetary radius and mass) are shown. In order to discuss our results in a broader context, we also show the location of several known brown dwarfs and low-mass stars orbiting around a late-K or M dwarf primary (Mireles et al. 2020; Carmichael 2023).

The figure shows two regimes of the mass-radius relationship. Up to planetary masses of the order of ~0.2 MJ, the planetary radius steadily increases with the planetary mass. For higher masses, the radius remains constant and even shows a slight decrease with increasing planetary mass. TOI-5375 b is clearly located in the region of the diagram where Rp is roughly constant (~1 RJ) with Mp. The two regimes of the mass-radius relationship are related to the boundary between objects with a dominating hydrogen and helium composition and objects that consist of different composition and different mass-radius relationships (see e.g. Bashi et al. 2017, and references therein).

The TOI-5375 system constitutes an interesting case study that can help us shed light on the tidal interactions between close companions and their host stars. With an orbital period for TOI-5375 b of 1.72 d, a rotational period for the host star of 1.97 d, the system is likely being observed in the short time frame in which the star is close to synchronising its rotation with the orbital period of the companion. Tidally driven orbital circularisation and synchronisation is a well-known phenomenon in binary systems containing low-mass stars (e.g. Zahn & Bouchet 1989), whilst the evidence for similar phenomena in systems harbouring short period hot Jupiter planets has been recently discussed (Birkby et al. 2014; Maciejewski et al. 2016; Patra et al. 2017; Wilkins et al. 2017; Bouma et al. 2019; Yee et al. 2020).

Taking into account the properties of the TOI-5375 system, tidal dissipation should be mainly driven by inertial waves (Ogilvie & Lin 2007). Under these circumstances, and following the model described in Barker (2020), we derive a tidal quality factor due to inertial waves of QIW=3.9×105${{Q'}_{{\rm{IW}}}} = 3.9 \times {10^5}$. We used the formulae of Leconte et al. (2010) which take into account the rotation of the star in the evolution of the system parameters to derive the characteristic timescales (the e-folding time τ of a given parameter q is defined as 1τ=1qdqdt${\textstyle{1 \over \tau }} = {\textstyle{1 \over q}}{\textstyle{{{\rm{d}}q} \over {{\rm{d}}t}}}$) for the TOI-5375 system. We derive an estimated timescale for tidal synchronisation of the stellar spin of only ~3.36 Myr and a circularisation timescale of ~6.76 Myr. Such a synchronisation timescale is much shorter than the angular momentum loss timescale of the stellar magnetised wind (of the order of 100 Myr given the fast rotation of the star), thus supporting the statement of Sect. 3 that we cannot use gyrochronology to estimate the age of this star. As a matter of fact, the rotation period of TOI-5375 is expected to be longer than about 10 days from the central plot in Fig. 6, given that it is likely to belong to the thin disc stars, that is, it is probably older than the Hyades. Conversely, it shows a rotation period as short as 1.97 days that we interpret as the consequence of the tidal acceleration given an estimated present synchronisation timescale much shorter than the expected magnetic braking timescale. The characteristic timescale for the semi-major axis decay is ~115 Myr, while the obliquity of the orbit has a decay timescale of ~0.72 Myr. The total angular momentum of the system is about 1.61 times larger than the minimum value required to reach a final stable tidal equilibrium with the star and the low-mass companion fully synchronised on a circular orbit. After that stage, the orbit will slowly decay as a result of the angular momentum losses in the magnetised stellar wind of the star (see Damiani & Lanza 2015).

Given the mass ratio between the star and its companion, M / Mb ~ 8.7, it is likely that the TOI-5375 system was formed as a low-mass binary system. One possibility is that TOI-5375 and its companion formed from the same molecular cloud that fragmented into two objects. Another possibility is that TOI-5375 formed a protoplanetary disc susceptible to gravitational instabilities. A relationship between the disc mass and the stellar mass of the form MdiscM1.2${M_{{\rm{disc}}}} \propto M_ \star ^{1.2}$ has been suggested (Alibert et al. 2011). Thus, the relatively low-mass protoplanetary disc around an Μ dwarf star like TOI-5375 avoided the formation of two stellar objects of equal mass. Disc fragmentation has been shown to successfully explain the formation of relatively close binaries, and due to the interactions with their natal disc, significant migration can also occur (e.g. Kratter 2011). Simulations show that discs with a moderate size and mass (100 au, 0.25 M) might fragment around late-K/early-M dwarfs forming mainly brown dwarfs, but also low-mass hydrogen-burning stars in generally close-in orbits (Stamatellos et al. 2011). Therefore, another piece of evidence supporting the disc fragmentation might come from the fact that TOI-5375 and its companion are very close and the star is to align its stellar spin with the binary orbit. However, we should note that the current spin-orbit alignment of the star may not be primordial, but it could have been established recently because the present timescale for the obliquity damping by tides is only 0.72 Myr (see above).

thumbnail Fig. 21

Simulated observed spectrum by the Ariel mission of TOI-5375 b generated from a PHOENIX/BT-SETTL atmospheric model. The grey line shows the spectra degraded to a resolution of R = 50, while the blue circles indicate the expected fluxes at Ariel wavelengths. Some characteristic spectral features of a brown dwarf are indicated (see e.g. Burrows et al. 2002; Burgasser et al. 2006; Sorahana & Yamamura 2012).

6 Conclusions

In this work, we report the discovery and confirmation of a transiting companion in the boundary between the brown dwarfs and the very-low-mass stars, TOI-5375 b, orbiting an active M0.5 dwarf with a period of 1.72 d. We make use of the available TESS photometric data and high-resolution HARPS-N spectra to model the signals present in the data. We derive a rotational period of 1.97 d for the host star and conclude that the orbit of the companion decays quickly and that the system should be in the phase in which the star is very close to synchronising its rotation with the orbital period of the companion.

TOI-5375 b, with a mass of 77 ± 8 MJ and a radius of 0.99 ± 0.16 RJ, joins the population of massive substellar or very-low-mass stellar companions in close-in orbits around low-mass stars that is currently emerging. Further studies of the properties of these objects will help us explore our models of low-mass companion formation and migration as well as to study how the tidal interactions with their host stars affect their evolution.

We can measure the change in the orbit of TOI-5375 b by computing the associated shift in transit arrival time (transits should occur earlier as the companion’s orbit decay). We obtain a shift in the transit time of 20.53 seconds in a timespan of 10 yr. In our analysis, we combined data from four TESS sectors containing a total of 53 transits. The precision in the estimated transit-mid time (from the simultaneous analysis of all transits) is of the order of 5 × 10−5 d or ~4.32 s, while the relative error in the transit's depth is of the order of ~0.2%. Therefore, while TOI-5375 is a good candidate to observe TTVs due to the effects of tidal waves, a long-term monitoring observation from ground-based facilities of this star would be needed (unfortunately the star is too faint to be observed by the CHEOPS satellite).

The TOI-5375 system can also be an interesting target for ongoing and future space missions aimed at the characterisation of the atmospheres of substellar companions like the James Webb Space Telescope (Gardner et al. 2006) or Ariel (Tinetti et al. 2018). We computed the transmission spectroscopy metric (TSM) as well as the emission spectroscopy metric (ESM) following the procedure defined in Kempton et al. (2018). The derived TSM value, 0.94, is significantly low, but we caution that TSM values were defined for small, rocky planets, so it is unclear whether or not they are still valid for a massive brown dwarf like TOI-5375 b or not. On the other hand, we obtain an ESM value of 125, which suggests that TOI-5375 b is indeed a good target for atmospheric characterisation.

In order to test what we might expect in the future, we generate a synthetic spectrum of TOI-5375 b as observed by Ariel. For this purpose, we made use of a PHOENIX/BT-SETTL atmospheric model (Allard et al. 2012) with parameters Teff = 1000 K, log ɡ = 4.5 cm s−2, and [Fe/H] = +0.00 dex. We convolved the original model with a Gaussian filter in order to match the expected resolution for Ariel Tier 2 observations (R ~ 50 at 1.95 − 3.95 µm) and computed the expected fluxes at Ariel wavelengths. The synthetic spectrum is shown in Fig. 21. The comparison of the spectral properties of TOI-5375 b with other brown dwarfs and stars will help us discover whether their atmospheric properties are equal or deviating, which may be a signspot for different formation scenarios.

We acknowledge that while this manuscript was under the referee review process, an independent analysis of the TOI-5375 system has been published (Lambert et al. 2023). These authors found similar results to ours. Nevertheless, we have intentionally maintained our analysis completely independent of these authors.

Acknowledgements

This paper includes data collected by the TESS mission. Funding for the TESS mission is provided by the NASA’s Science Mission Directorate. Based on observations collected at Copernico 1.82 m telescope (Asiago Mount Ekar, Italy) INAF – Osservatorio Astronomico di Padova. This work has been supported by the PRIN-INAF 2019 “Planetary systems at young ages (PLATEA)” and ASI-INAF agreement n.2018-16-HH.0. J. M., A. P., S. C. and C. D. M. acknowledge support from the Accordo Attuativo ASI-INAF n. 2021-5-HH.0, Partecipazione alla fase B2/C della missione Ariel (ref. G. Micela). We sincerely appreciate the careful reading of the manuscript and the constructive comments of an anonymous referee.

Appendix A Results from the joint analysis: TESS and ground-based photometry

Table A.1 provides the best-fit values obtained for the simultaneous fit of the TESS SAP light curve and the Asiago observations, as described in Sect. 4.4. The best transit model fit is shown in Figure A.1.

thumbnail Fig. A.1

Photometric modelling of TOI-5375 (TESS and Asiago photometry). Top: Folded transits in the light curve after the subtraction of the stellar activity signal. The black line shows the best model fit of the transits. Bottom: Residuals of the light curve after the subtraction of the companion transit model.

Table A.1

Best-fit values obtained for the simultaneous fit of the TESS SAP light curve and the Asiago observations.

Appendix В Activity indexes

Figure B.1 shows the temporal variation of the different activity indexes measured in this work, while the periodogram analysis is given in Fig. B.2. Table B.1 provides the measured values with their corresponding uncertainties.

thumbnail Fig. B.1

From top to bottom: Ca II H & K index, Hα index, and Na I D1 D2 index time series for TOI-5375.

Table B.1

Activity indexes measured for TOI-5375.

thumbnail Fig. B.2

From top to bottom: GLS periodogram of the Ca II H & K, Нα, and Na I activity indexes. The vertical line indicates the period of TOI-5375 b at 1.72 d. Values corresponding to a FAP of 10%, 1%, and 0.1% are shown with horizontal grey lines.

References

  1. Alibert, Y., Mordasini, C., & Benz, W. 2011, A&A, 526, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Allard, F., Homeier, D., & Freytag, B. 2012, Philos. Trans. R. Soc. A, 370, 2765 [NASA ADS] [CrossRef] [Google Scholar]
  3. Alvarado-Montes, J. A. 2022, MNRAS, 517, 2831 [NASA ADS] [CrossRef] [Google Scholar]
  4. Anglada-Escudé, G., & Butler, R. P. 2012, ApJS, 200, 15 [Google Scholar]
  5. Ballerini, P., Micela, G., Lanza, A. F., & Pagano, I. 2012, A&A, 539, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Barker, A. J. 2020, MNRAS, 498, 2270 [Google Scholar]
  9. Bashi, D., Helled, R., Zucker, S., & Mordasini, C. 2017, A&A, 604, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Batygin, K., Bodenheimer, P. H., & Laughlin, G. P. 2016, ApJ, 829, 114 [Google Scholar]
  11. Benatti, S., Nardiello, D., Malavolta, L., et al. 2019, A&A, 630, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bilir, S., Ak, S., Karaali, S., et al. 2008, MNRAS, 384, 1178 [Google Scholar]
  13. Birkby, J. L., Cappetta, M., Cruz, P., et al. 2014, MNRAS, 440, 1470 [NASA ADS] [CrossRef] [Google Scholar]
  14. Boley, A. C., Granados Contreras, A. P., & Gladman, B. 2016, ApJ, 817, L17 [Google Scholar]
  15. Bonfils, X., Delfosse, X., Udry, S., et al. 2013, A&A, 549, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Borsato, L., Malavolta, L., Piotto, G., et al. 2019, MNRAS, 484, 3233 [NASA ADS] [CrossRef] [Google Scholar]
  17. Bouma, L. G., Winn, J. N., Baxter, C., et al. 2019, AJ, 157, 217 [NASA ADS] [CrossRef] [Google Scholar]
  18. Buchner, J., Georgakakis, A., Nandra, K., et al. 2014, A&A, 564, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Burgasser, A. J., Geballe, T. R., Leggett, S. K., Kirkpatrick, J. D., & Golimowski, D. A. 2006, ApJ, 637, 1067 [Google Scholar]
  20. Burrows, A., Burgasser, A. J., Kirkpatrick, J. D., et al. 2002, ApJ, 573, 394 [NASA ADS] [CrossRef] [Google Scholar]
  21. Butler, R. P., Johnson, J. A., Marcy, G. W., et al. 2006, PASP, 118, 1685 [Google Scholar]
  22. Carleo, I., Benatti, S., Lanza, A. F., et al. 2018, A&A, 613, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Carleo, I., Malavolta, L., Lanza, A. F., et al. 2020, A&A, 638, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Carleo, I., Desidera, S., Nardiello, D., et al. 2021, A&A, 645, A71 [EDP Sciences] [Google Scholar]
  25. Carmichael, T. W. 2023, MNRAS, 519, 5177 [NASA ADS] [CrossRef] [Google Scholar]
  26. Casagrande, L., Portinari, L., & Flynn, C. 2006, MNRAS, 373, 13 [Google Scholar]
  27. Colombo, S., Petralia, A., & Micela, G. 2022, A&A, 661, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Cosentino, R., Lovis, C., Pepe, F., et al. 2012, Proc. SPIE, 8446, 84461V [Google Scholar]
  29. Covino, E., Esposito, M., Barbieri, M., et al. 2013, A&A, 554, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Cumming, A., Butler, R. P., Marcy, G. W., et al. 2008, PASP, 120, 531 [CrossRef] [Google Scholar]
  31. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
  32. Cutri, R. M., Wright, E. L., Conrow, T., et al. 2021, VizieR Online Data Catalog: II/328 [Google Scholar]
  33. Damasso, M., Lanza, A. F., Benatti, S., et al. 2020, A&A, 642, A133 [EDP Sciences] [Google Scholar]
  34. Damiani, C., & Lanza, A. F. 2015, A&A, 574, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Donati, J. F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658 [Google Scholar]
  36. Dreizler, S., Crossfield, I. J. M., Kossakowski, D., et al. 2020, A&A, 644, A127 [EDP Sciences] [Google Scholar]
  37. Dressing, C. D., & Charbonneau, D. 2013, ApJ, 767, 95 [Google Scholar]
  38. Dressing, C. D., & Charbonneau, D. 2015, ApJ, 807, 45 [Google Scholar]
  39. Dumusque, X., Cretignier, M., Sosnowska, D., et al. 2021, A&A, 648, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Eastman, J., Siverd, R., & Gaudi, B. S. 2010, PASP, 122, 935 [Google Scholar]
  41. Eggen, O. J. 1984, AJ, 89, 1358 [NASA ADS] [CrossRef] [Google Scholar]
  42. Eggen, O. J. 1989, PASP, 101, 366 [NASA ADS] [CrossRef] [Google Scholar]
  43. Endl, M., Kürster, M., Els, S., Hatzes, A. P., & Cochran, W. D. 2001, A&A, 374, 675 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Endl, M., Cochran, W. D., Kürster, M., et al. 2006, ApJ, 649, 436 [Google Scholar]
  45. ESA 1997, The HIPPARCOS and Tycho catalogues. Astrometric and photometric star catalogues derived from the ESA Hipparcos Space Astrometry Mission (Noordwijk, The Netherlands: ESA Publications Division) [Google Scholar]
  46. Feroz, F., Hobson, M. P., Cameron, E., & Pettitt, A. N. 2019, Open J. Astrophys., 2, 10 [Google Scholar]
  47. Fleming, T. A., Schmitt, J. H. M. M., & Giampapa, M. S. 1995, ApJ, 450, 401 [NASA ADS] [CrossRef] [Google Scholar]
  48. Flower, P. J. 1996, ApJ, 469, 355 [NASA ADS] [CrossRef] [Google Scholar]
  49. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  50. Francis, C., & Anderson, E. 2009, New A, 14, 615 [NASA ADS] [CrossRef] [Google Scholar]
  51. Gagné, J., Mamajek, E. E., Malo, L., et al. 2018, ApJ, 856, 23 [Google Scholar]
  52. Gaia Collaboration 2022, VizieR Online Data Catalog: I/355 [Google Scholar]
  53. Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, Space Sci. Rev., 123, 485 [Google Scholar]
  54. Gelman, A., & Rubin, D. B. 1992, Statis. Sci., 7, 457 [NASA ADS] [Google Scholar]
  55. Goodman, J., & Weare, J. 2010, Commun. Appl. Math. Comput. Sci., 5, 65 [Google Scholar]
  56. Henry, T. J., & McCarthy, J., Donald W. 1993, AJ, 106, 773 [NASA ADS] [CrossRef] [Google Scholar]
  57. Henry, T. J., Soderblom, D. R., Donahue, R. A., & Baliunas, S. L. 1996, AJ, 111, 439 [Google Scholar]
  58. Howard, A. W., Marcy, G. W., Bryson, S. T., et al. 2012, ApJS, 201, 15 [Google Scholar]
  59. Hunter, A., Macgregor, A., Szabo, T., Wellington, C., & Bellgard, M. 2012, Source Code Biol. Med., 7, 1 [Google Scholar]
  60. Husser, T. O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [Google Scholar]
  62. Kanodia, S., Cañas, C. I., Stefansson, G., et al. 2020, ApJ, 899, 29 [Google Scholar]
  63. Kanodia, S., Libby-Roberts, J., Cañas, C. I., et al. 2022, AJ, 164, 81 [NASA ADS] [CrossRef] [Google Scholar]
  64. Kempton, E. M. R., Bean, J. L., Louie, D. R., et al. 2018, PASP, 130, 114401 [Google Scholar]
  65. Kipping, D. M. 2013, MNRAS, 435, 2152 [Google Scholar]
  66. Kratter, K. M. 2011, ASP Conf. Ser., 447, 47 [NASA ADS] [Google Scholar]
  67. Kreidberg, L. 2015, PASP, 127, 1161 [Google Scholar]
  68. Lambert, M., Bender, C. F., Kanodia, S., et al. 2023, AJ, 165, 218 [NASA ADS] [CrossRef] [Google Scholar]
  69. Leconte, J., Chabrier, G., Baraffe, I., & Levrard, B. 2010, A&A, 516, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Lin, D. N. C., Bodenheimer, P., & Richardson, D. C. 1996, Nature, 380, 606 [Google Scholar]
  71. Loeb, A., & Gaudi, B. S. 2003, ApJ, 588, L117 [Google Scholar]
  72. Lovis, C., & Pepe, F. 2007, A&A, 468, 1115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Ma, B., & Ge, J. 2014, MNRAS, 439, 2781 [Google Scholar]
  74. Maciejewski, G., Dimitrov, D., Fernández, M., et al. 2016, A&A, 588, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Malavolta, L. 2016, Astrophysics Source Code Library [record ascl:1612.008] [Google Scholar]
  76. Malavolta, L., Nascimbeni, V., Piotto, G., et al. 2016, A&A, 588, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Malavolta, L., Mayo, A. W., Louden, T., et al. 2018, AJ, 155, 107 [NASA ADS] [CrossRef] [Google Scholar]
  78. Maldonado, J. 2012, Ph.D. Thesis, Universidad Autónoma de Madrid, Spain [Google Scholar]
  79. Maldonado, J., & Villaver, E. 2017, A&A, 602, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Maldonado, J., Martínez-Arnáiz, R. M., Eiroa, C., Montes, D., & Montesinos, B. 2010, A&A, 521, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Maldonado, J., Affer, L., Micela, G., et al. 2015, A&A, 577, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Maldonado, J., Villaver, E., & Eiroa, C. 2018, A&A, 612, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Maldonado, J., Phillips, D. F., Dumusque, X., et al. 2019, A&A, 627, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Maldonado, J., Micela, G., Baratella, M., et al. 2020, A&A, 644, A68 [EDP Sciences] [Google Scholar]
  85. Mamajek, E. E., & Hillenbrand, L. A. 2008, ApJ, 687, 1264 [Google Scholar]
  86. Mantovan, G., Montalto, M., Piotto, G., et al. 2022, MNRAS, 516, 4432 [NASA ADS] [CrossRef] [Google Scholar]
  87. Mata Sánchez, D., González Hernández, J. I., Israelian, G., et al. 2014, A&A, 566, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Mathieu, R. D., & Mazeh, T. 1988, ApJ, 326, 256 [NASA ADS] [CrossRef] [Google Scholar]
  89. Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [Google Scholar]
  90. Mazeh, T. 2008, EAS Pub. Ser., 29, 1 [Google Scholar]
  91. Mireles, I., Shporer, A., Grieves, N., et al. 2020, AJ, 160, 133 [Google Scholar]
  92. Montes, D., López-Santiago, J., Gálvez, M. C., et al. 2001, MNRAS, 328, 45 [NASA ADS] [CrossRef] [Google Scholar]
  93. Morton, T. D. 2012, ApJ, 761, 6 [NASA ADS] [CrossRef] [Google Scholar]
  94. Morton, T. D. 2015, Astrophysics Source Code Library [record ascl:1503.011] [Google Scholar]
  95. Nardiello, D. 2020, MNRAS, 498, 5972 [Google Scholar]
  96. Nardiello, D., Piotto, G., Deleuil, M., et al. 2020, MNRAS, 495, 4924 [Google Scholar]
  97. Nardiello, D., Malavolta, L., Desidera, S., et al. 2022, A&A, 664, A163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  98. Nascimbeni, V., Piotto, G., Bedin, L. R., & Damasso, M. 2011, A&A, 527, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Nascimbeni, V., Cunial, A., Murabito, S., et al. 2013, A&A, 549, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Neves, V., Bonfils, X., Santos, N. C., et al. 2012, A&A, 538, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  101. Newton, E. R., Mann, A. W., Tofflemire, B. M., et al. 2019, ApJ, 880, L17 [Google Scholar]
  102. Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H. 1984, ApJ, 279, 763 [Google Scholar]
  103. Ogilvie, G. I., & Lin, D. N. C. 2007, ApJ, 661, 1180 [Google Scholar]
  104. Oshagh, M., Santos, N. C., Boisse, I., et al. 2013, A&A, 556, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Paegert, M., Stassun, K. G., Collins, K. A., et al. 2022, VizieR Online Data Catalog: IV/39 [Google Scholar]
  106. Parviainen, H., Aigrain, S., Thatte, N., et al. 2015, MNRAS, 453, 3875 [NASA ADS] [Google Scholar]
  107. Patra, K. C., Winn, J. N., Holman, M. J., et al. 2017, AJ, 154, 4 [Google Scholar]
  108. Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  109. Perger, M., García-Piquer, A., Ribas, I., et al. 2017, A&A, 598, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  110. Pinamonti, M., Sozzetti, A., Maldonado, J., et al. 2022, A&A, 664, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  111. Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525 [Google Scholar]
  112. Plavchan, P., Barclay, T., Gagné, J., et al. 2020, Nature, 582, 497 [Google Scholar]
  113. Prosser, C. F., Shetrone, M. D., Dasgupta, A., et al. 1995, PASP, 107, 211 [NASA ADS] [CrossRef] [Google Scholar]
  114. Radick, R. R., Thompson, D. T., Lockwood, G. W., Duncan, D. K., & Baggett, W. E. 1987, ApJ, 321, 459 [NASA ADS] [CrossRef] [Google Scholar]
  115. Rainer, M., Borsa, F., & Affer, L. 2020, Exp. Astron., 49, 73 [Google Scholar]
  116. Rainer, M., Borsa, F., Pino, L., et al. 2021, A&A, 649, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  117. Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telesc. Instrum. Syst., 1, 014003 [Google Scholar]
  118. Rizzuto, A. C., Mann, A. W., Vanderburg, A., Kraus, A. L., & Covey, K. R. 2017, AJ, 154, 224 [Google Scholar]
  119. Ryabchikova, T., Piskunov, N., Kurucz, R. L., et al. 2015, Phys. Scr., 90, 054005 [Google Scholar]
  120. Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 [Google Scholar]
  121. Sorahana, S., & Yamamura, I. 2012, ApJ, 760, 151 [NASA ADS] [CrossRef] [Google Scholar]
  122. Spina, L., Randich, S., Magrini, L., et al. 2017, A&A, 601, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  123. Stamatellos, D., Maury, A., Whitworth, A., & André, P. 2011, MNRAS, 413, 1787 [Google Scholar]
  124. Stauffer, J. R., Caillault, J. P., Gagne, M., Prosser, C. F., & Hartmann, L. W. 1994, ApJS, 91, 625 [Google Scholar]
  125. Stauffer, J. R., Jones, B. F., Backman, D., et al. 2003, AJ, 126, 833 [NASA ADS] [CrossRef] [Google Scholar]
  126. Stern, R. A., Schmitt, J. H. M. M., & Kahabka, P. T. 1995, ApJ, 448, 683 [Google Scholar]
  127. Storn, R., & Price, K. 1997, J. Global Optimiz., 11, 341 [Google Scholar]
  128. Suárez Mascareño, A., Rebolo, R., González Hernández, J. I., & Esposito, M. 2015, MNRAS, 452, 2745 [Google Scholar]
  129. Tinetti, G., Drossart, P., Eccleston, P., et al. 2018, Exp. Astron., 46, 135 [NASA ADS] [CrossRef] [Google Scholar]
  130. Tokovinin, A., Thomas, S., Sterzik, M., & Udry, S. 2006, A&A, 450, 681 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  131. Tuomi, M., Jones, H. R. A., Barnes, J. R., Anglada-Escudé, G., & Jenkins, J. S. 2014, MNRAS, 441, 1545 [Google Scholar]
  132. Voges, W., Aschenbach, B., Boller, T., et al. 1999, A&A, 349, 389 [NASA ADS] [Google Scholar]
  133. Wilkins, A. N., Delrez, L., Barker, A. J., et al. 2017, ApJ, 836, L24 [NASA ADS] [CrossRef] [Google Scholar]
  134. Wolszczan, A., & Frail, D. A. 1992, Nature, 355, 145 [CrossRef] [Google Scholar]
  135. Yee, S. W., Winn, J. N., Knutson, H. A., et al. 2020, ApJ, 888, L5 [Google Scholar]
  136. Zahn, J. P., & Bouchet, L. 1989, A&A, 223, 112 [NASA ADS] [Google Scholar]
  137. Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [CrossRef] [EDP Sciences] [Google Scholar]
  138. Zechmeister, M., Reiners, A., Amado, P. J., et al. 2018, A&A, 609, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

4

The Hyades star with a rotation period equal to 3.66 d is VB 190. Radick et al. (1987) discussed two possibilities to explain its peculiar behaviour. First, the star lies somewhat above the cluster main sequence, suggesting that it might be a binary. Then, its rapid rotation might be due to tidal coupling with the companion. Second, about one-third of the M-type Hyades stars are known to have projected rotational velocities larger than 10 km s−1. VB 190 may be a member of this population.

5

Original B, V, R, I magnitudes from the isochrones have been transformed into (VKS) following the relationship by Bilir et al. (2008).

8

We note that two additional candidates have been identified by TESS, one using data from Sectors 40 and 47, and a second candidate in Sector 60. A careful look at their corresponding TESS reports shows that they are likely FP due to the procedure used to flatten the light curve used by TESS. This might happen when the star shows a high level of variability and when, as it is our case, flares are present in the light curve.

11

Computed with the Python scipy ks_1samp routine.

All Tables

Table 1

Comparison between the different methods used to extract the RVs.

Table 2

Physical properties of TOI-5375.

Table 3

Best-fit values obtained for the star-planet model.

Table 4

Results of the K-S tests performed in this work.

Table 5

RVs (with their corresponding uncertainties) of TOI-5375 used in this work.

Table A.1

Best-fit values obtained for the simultaneous fit of the TESS SAP light curve and the Asiago observations.

Table B.1

Activity indexes measured for TOI-5375.

All Figures

thumbnail Fig. 1

Original TESS photometry time series. The energetic events discussed in Sect. 4.6 are shown with different symbols according to their Sector. Within each Sector, different colours indicate that the event occurs at a different period (when phase-folded to the stellar rotation period).

In the text
thumbnail Fig. 2

GLS periodogram of the TESS photometry dataset. Values corresponding to a False Alarm Probability (FAP) of 10, 1, and 0.1% are shown with horizontal grey lines. The vertical line indicates the period at 1.9692 d.

In the text
thumbnail Fig. 3

Examples of standardised (mean subtracted and division by the standard deviation) CCF profiles obtained by the DRS using different binary masks as well as the LSD profile.

In the text
thumbnail Fig. 4

RVs obtained using different masks and methods versus LSD RVs. For TERRA and SERVAL RVs, we subtracted the mean LSD RV for the comparison. Each panel shows the mean value and standard deviation of the difference RVs other methods - LSD RVs. Grey lines show the one-to-one relationship.

In the text
thumbnail Fig. 5

Spectral energy distribution of TOI-5375. A PHOENIX/BT-SETTL atmospheric model (Teff = 3885 K, log ɡ = 4.66 cm s−2, [Fe/H] = +0.26 dex) is overplotted. Errors in the observed fluxes are within the symbol’s size.

In the text
thumbnail Fig. 6

Age diagnostics for TOI-5375. Left: fractional X-ray luminosity log(Lx/LBol) versus colour index (BV). Centre: rotation period as a function of the (BV) colour. Two gyrochrones (at the ages of the Pleiades and the Hyades) have been overplotted for comparison. Right: log(RHK)$\left( {{{R'}_{{\rm{HK}}}}} \right)$ versus colour index (BV). Dashed lines are the limits for very active, active, inactive, and very inactive stars, according to Henry et al. (1996).

In the text
thumbnail Fig. 7

Colour-magnitude diagram. Left: Mv as a function of (BV) colour. Right: Mv as a function of (VKs) colour.

In the text
thumbnail Fig. 8

(U, V) plane showing the position of the young disc population (dashed line) as defined by Eggen (1984, 1989) as well as the known kinematic groups (large crosses and asterisks) in the solar neighbourhood (see Maldonado et al. 2010; Maldonado 2012, and references therein). The dotted line represents the velocity ellipsoid determined by Francis & Anderson (2009). The position of TOI-5375 is shown with a black star (uncertainties are within the symbol size). All velocities are in the Solar System reference frame.

In the text
thumbnail Fig. 9

Gaia stars identified in the TESS field. The position of TOI-5375 is shown with the letter ‘ť.

In the text
thumbnail Fig. 10

In- or out-oſ-transit centroid tests calculated for each sector. Black dots are background stars.

In the text
thumbnail Fig. 11

Posterior distributions of the Bayesian fit of the photometric time series.

In the text
thumbnail Fig. 12

Photometric modelling of TOI-5375. Top: folded transits in the light curve after the subtraction of the stellar activity signal. The green line shows the best model fit of the transits. Bottom: residuals of the light curve after the subtraction of the companion transit model.

In the text
thumbnail Fig. 13

Photometric observations from Asiago, folded on the transit model fit period of TOI-5375 b. Left: October 25, 2022, observations. Right: December 19, 2022, observations. The green line shows the best fit to the Asiago dataset, while for comparison purposes we also show in red the best fit derived for the TESS light curves. The residuals of the fit to the Asiago ground-based photometry are shown in the bottom panel.

In the text
thumbnail Fig. 14

Histogram of the derived Rb/R values from each individual TESS and Asiago transit.

In the text
thumbnail Fig. 15

Residuals of the timings of mid-transits versus a linear ephemeris.

In the text
thumbnail Fig. 16

Histogram of impulsive events (blue empty line) and amount of energy emitted (orange filled histogram) phase-folded to the companion orbital period (left), the synodic period (centre) and the stellar period (right).

In the text
thumbnail Fig. 17

GLS periodogram of the RV dataset. Values corresponding to a FAP of 10, 1, and 0.1% are shown with horizontal grey lines. The vertical line indicates the period of TOI-5375 b at 1.72 d.

In the text
thumbnail Fig. 18

Posterior distributions of the Bayesian fit of the RVs.

In the text
thumbnail Fig. 19

Best-fit planetary model for TOI-5375 b (top), and RV residuals (bottom).

In the text
thumbnail Fig. 20

Known substellar companions and low-mass stars (radius versus mass) around late-K and M dwarfs (as listed at http://exoplanet.eu/). Different colours and symbol sizes indicate the effective temperature of the host star. Models for substellar objects and low-mass stars with solar metallicity from Baraffe et al. (2015) are over-plotted with grey lines.

In the text
thumbnail Fig. 21

Simulated observed spectrum by the Ariel mission of TOI-5375 b generated from a PHOENIX/BT-SETTL atmospheric model. The grey line shows the spectra degraded to a resolution of R = 50, while the blue circles indicate the expected fluxes at Ariel wavelengths. Some characteristic spectral features of a brown dwarf are indicated (see e.g. Burrows et al. 2002; Burgasser et al. 2006; Sorahana & Yamamura 2012).

In the text
thumbnail Fig. A.1

Photometric modelling of TOI-5375 (TESS and Asiago photometry). Top: Folded transits in the light curve after the subtraction of the stellar activity signal. The black line shows the best model fit of the transits. Bottom: Residuals of the light curve after the subtraction of the companion transit model.

In the text
thumbnail Fig. B.1

From top to bottom: Ca II H & K index, Hα index, and Na I D1 D2 index time series for TOI-5375.

In the text
thumbnail Fig. B.2

From top to bottom: GLS periodogram of the Ca II H & K, Нα, and Na I activity indexes. The vertical line indicates the period of TOI-5375 b at 1.72 d. Values corresponding to a FAP of 10%, 1%, and 0.1% are shown with horizontal grey lines.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.