TATOO: a tidal-chronology standalone tool to estimate the age of massive close-in planetary systems

The presence of a massive close-in planet with an orbital period of a few days or less around a low-mass star can possibly results in a strong variation of the properties of the central star. Indeed, star-planet tidal interactions generate exchanges of angular momentum that can results in tidal spin-up. This effect could then lead to gyrochronological ages biased towards younger ages. Aims. This article provides the community with TATOO, a standalone tool based on tidal-chronology, to estimate the age of a massive close-in planetary system by only using its observed properties: mass of the planet and the star, stellar rotational and planetary orbital periods. I used the numerical code described in Gallet et al. (2018) to create a large multi-parametric grid of synthetic star-planet systems evolution, and 3D interpolation method to provide a fairly precise age estimate, using tidal-chronology technique, of any given planetary systems composed of at least of one massive close-in planet. About half of the planetary systems investigated in this work are subject to tidal spin-up bias. I pointed out that this bias linearly scales with the ratio between rotation to orbital period, making this quantity a useful proxy to rapidly investigate the need to use tidal-chronology. Moreover, while being model dependent, TATOO can also be used even if no rotational departure is present. In that case, it gives results in agreement with the classical gyrochronological analysis. TATOO is a useful tool specifically designed for massive close-in planetary systems that can also be used as a classical gyrochronological tool. For now it is the only publicly available software to estimate the age of massive close-in planetary systems subject to tidal spin-up. In that sense, tidal-chronology can be seen as a first order correction of the impact of tidal interaction on gyrochronology.


Introduction
Together with rotation, radius, and luminosity, age is a fundamental physical parameter for both stellar and exoplanetary fields, and its determination for a given star or planetary system is currently a key issue in astrophysics. In stellar physics, stellar rotation has been found to be involved in most of the physical mechanisms: magnetic field strength via dynamo effect (Jouve & Brun 2007), orbital evolution in planetary systems (Bolmont & Mathis 2016;Gallet et al. 2018), and stellar angular momentum evolution (Gallet & Bouvier 2013, 2015. Estimating the age of an isolated star or a star that belongs to a planetary system can thus help to constrain its internal structure and chemical species content (e.g. light elements such as Lithium), through the use of stellar evolution models, and to put boundary conditions on planet formation and migration models (Ida & Lin 2008;Mordasini et al. 2009Mordasini et al. , 2012Alibert et al. 2013;Amard et al. 2016;Gallet et al. 2017;Amard et al. 2019). It also has an important role in adding constraints on the star-planet interaction efficiency and characteristic timescale ) that could help improve star-planet magnetic and tidal interaction scenarios (Strugarek et al. 2017;Gallet & Delorme 2019;Benbakoura et al. 2019).
Send offprint requests to: F. Gallet, However, stellar ages cannot be directly measured, except for the Sun, and as a consequence have to be inferred either from theoretical models (e.g. nucleocosmochronology) or by using empirical laws. Gyrochronology (Barnes 2003(Barnes , 2007(Barnes , 2010Angus et al. 2019) and magnetochronology (Vidotto et al. 2014) are two examples of empirical methods that use surface rotation period measurements and magnetic field strength estimation (see Soderblom 2010, for a list of age estimation methods).
Other model-dependant techniques can also be used such as isochrones fitting and main sequence turn-off ages, which both require accurate and precise stellar evolution codes. These techniques suffer from the presence of systematics that are linked to fitting accuracy and uncertain distances (Lebreton et al. 2014). Moreover, they can only be used if there were no interactions between the star and its close environment during the main sequence phase. Indeed, in the case of past and already finished interactions, it is difficult to quantify the impact of these external interactions on the surface stellar rotation by only using the observations. Hence, it is hard to estimate the true and apparent age of the star. In that case, incompatible age estimations provided by different techniques could be found to be due to this external origin (e.g. Bonnefoy et al. 2018). However, for an on-going interaction, for example if a massive planet (M p > 1 M jup ) is seen in the vicinity of the star with a separation less than 0.1 au, these techniques can no longer be used since star-planet tidal interactions have most probably modified the evolution of the surface Article number, page 1 of 11 arXiv:2006.07880v2 [astro-ph.EP] 12 Jul 2020 A&A proofs: manuscript no. 38058corr rotation rate during the system's evolution (Gallet et al. 2018;Qureshi et al. 2018). Consequently, and as predicted in the literature (Maxted et al. 2015;Penev et al. 2018) and shown in Gallet & Delorme (2019), the ages of numerous star-planet systems are currently not known and might appear younger than they are because of tidal spin-up effects. In that case, it is yet possible to use tidal-chronology to get the age of the system from the measurement of the rotation period of the star and the semi-major axis of its orbiting massive planet (see Gallet & Delorme 2019, for a detailed presentation of the technique).
In Gallet & Delorme (2019) we introduced and described the tidal-chronology technique. By applying it to the specific system WASP-43, we pointed out the significant discrepancy between tidal-chronological and gyrochronological age estimation induced by tidal spin-up. In Gallet & Delorme (2019) we used a tailored model grid matching the specific observed properties of WASP-43 to perform this estimation. Since this process is time consuming because of the large initial parameter space of planetary systems to be explored, I decided to provide the community with a standalone tool dedicated to the estimation of the age of planetary systems (composed of one low-mass star and one massive close-in planet). This tool uses only the observation of the surface rotation rate of the host star and the current location of the planet orbiting it. This tool is based on a pre-computed generic grid in which the age of the user's system is interpolated, thus making the age estimation process easy and fast. This paper is separated as follows. In Sect. 2, I recall the basic principle of the tidal-chronology technique and describe the standalone tool TATOO (Tidal-chronology Age TOOl). I then apply this tool to known planetary systems to estimate their age and compare them to classical gyrochronology ones. Finally, I discuss the results in Section 3 and conclude in Section 4.

TATOO: A standalone tool
To properly estimate the age of a given system with the tidalchronology technique (see next section), one needs to construct a dense tailored grid of a planetary system's evolution. This tailored grid should contain the evolution of the stellar rotation period P rot, and planetary orbital period P orb,p of a given system that matches the observed stellar M and planetary M p mass, and this for a large range of initial conditions of stellar rotation P rot,init and planetary orbital period P orb,init . Hence, for the fixed observed properties M ; M p of the system, the tailored grid should explore every possible value of P orb,init for each value of P rot,init . For example, in the case of WASP-43 that we investigated in Gallet & Delorme (2019), we used a 0.71 M stellar model and a 2.052 M jup mass planet for which we explored the initial conditions of the stellar rotation rate and orbital period. This tailored gird is thus only valid for planetary systems that are identical to the WASP-43 ones.
The age of the system is then extracted from this tailored grid by using the observed pair P rot,obs ; P orb,obs . However, this procedure is extremely time consuming and the resulting tailored grid is only valid for the considered planetary system corresponding to the observed couple M ; M p . A generic numerical tool is thus required to provide us with a faster estimate of the age of a given observed system. In this section, I present such a tool named TATOO (https://github.com/GalletFlorian/TATOO) and based on the work of Gallet & Delorme (2019).
To allow a clear reading of the paper, the nomenclature used in this work is explained below: tailored grid, initial grid, pre-compiled age exploration grid: These are three different grids that will be introduced in the text. -P rot, and P orb,p : Without any specific label, they simply refer to the evolution of the stellar rotation and planetary orbital period. -P rot,init and P orb,init : The initial stellar rotation period and planetary orbital period used in the dense initial grid. They are distributed between 1-11 days and 0.1-1.1 P rot,init , respectively. -P rot,i and P orb,i : The stellar rotation period and planetary orbital period of the exploration grid. They are regularly distributed between 2 and 30 days, and 0.2 and 15 days, respectively (see Sect. 2.2), -P rot,obs and P orb,obs : The stellar rotation period and planetary orbital period of the user's favourite observed planetary system.
In the next section, I reiterate the general principles of the tidalchronology technique and how age estimation can be performed using such a technique.

Basic principle of tidal-chronology
In Gallet & Delorme (2019), and following the work of Gallet et al. (2018), we presented the concept of a new age estimation technique based only on the measurement of the surface rotation rate of the star and the location of a massive close-in planet around it: the tidal-chronology. In Gallet & Delorme (2019) we demonstrated that the observed pair composed of the stellar rotation period (P rot,obs ) and planetary orbital period (P orb,obs ) of a given star and massive close-in planet system is only retrieved at a unique age or during a short range of time. We hence highlighted that this characteristic could be used to estimate the age of massive close-in planetary systems. As mentioned above, a dense tailored grid of the evolution planetary systems first needs to be created. From this tailored grid, only systems in which the planet is still orbiting the star (i.e. systems that did not experience planetary engulfment) are selected and the ages at which the observed pair P rot,obs ; P orb,obs is retrieved are extracted. To estimate the validity of this age estimate, in Gallet & Delorme (2019) we introduced the S 2 quantity that needs to be minimized S 2 = (P orb,p − P orb,obs ) 2 σ 2 P orb,obs + (P rot, − P rot,obs ) 2 σ 2 where σ P orb,obs and σ P rot,obs are the errors of the observed P orb,obs and P rot,obs , and P rot, and P orb,p are the stellar rotation and planetary orbital period from the models. This method finally provides the most probable range for the age of the considered planetary system. The model used in this work to create the tailored grid is PROBE (PeRiod and OrBital Evolution code), which is the combination of the stellar angular momentum evolution code JEVOL described in Gallet & Bouvier (2013, 2015 with the modified orbital evolution model used in Bolmont & Mathis (2016). A detailed description of this model can be found in Gallet et al. (2018). To enable an easier reading of the paper, I will recall the main characteristics of the model.

Tidal dissipation
The tidal theory used in this work includes both equilibrium and dynamical tides and it is based on the tidal dissipation formalism described in Bolmont & Mathis (2016) that is parametrized following the simplified model by Ogilvie (2013). As in Bolmont & Mathis (2016) and Gallet et al. (2018), only the frequencyaveraged tidal dissipation is considered and the frequency response of the tides is neglected. Moreover, the dissipation of tidal inertial waves is only treated inside the convective envelope of the star (Mathis 2015;Mathis et al. 2016;Gallet et al. 2017). For now, the dissipation inside the planet's interior is not included as it is still hardly theoretically constrained. However, for Hot Jupiter cases with circular orbit, we can neglect this additional dissipation regarding the short characteristic timescales (10 5 yr) associated to the evolution of both their rotation rate and inclination angle (Leconte et al. 2010;Damiani & Mathis 2018). Finally, I neglected the impact of the magnetic star-planet interactions on the planetary orbit (see Strugarek et al. 2017), while it may have an effect on the rotational evolution of the central star during the early pre-main-sequence (hereafter PMS) phase (Strugarek et al. 2019). The migration timescale of planets during the main-sequence (hereafter MS) phase could be shortened by the inclusion of these additional dissipations, which could slightly change the age estimation given by tidal-chronology.
The impact of the evolution of the planetary orbit on the stellar angular momentum is evaluated using Gallet et al. (2017), who provide an estimation of tidal dissipation effects induced by gravitational interactions between the central star and the orbiting body. I note that this tidal formalism only allows the presence of one planet orbiting the star with a coplanar and circular orbit.

Stellar rotation
As pointed out in the previous section, the evolution of the planetary semi-major axis a strongly depends on the evolution of the surface rotation rate of the star via its impact on the tidal dissipation intensity (Mathis 2015;Mathis et al. 2016;Gallet et al. 2017;Bolmont et al. 2017). As a consequence, in that framework, modelling of the evolution of the stellar rotation should play a crucial role. In PROBE, the rotation rate of the star is modelled using the formalism described in Gallet & Bouvier (2013, 2015.

Stellar model
The evolution of the internal structure and the main physical properties of the star are provided by the stellar evolution code STAREVOL (see Lagarde et al. 2012;Amard et al. 2016Amard et al. , 2019, and references therein). A detailed comparison between STAREVOL and similar stellar evolution codes (including MESA, YREC, and PARSEC, see Choi et al. 2016;Spada et al. 2011;Tognelli et al. 2011, respectively) can be found in Amard et al. (2019). However, in this work I did not explore the possible impact of the choice of stellar model on the outcome of this standalone tool.
Finally, it is worth noting that as in gyrochronology, the tidalchronology technique will provide degenerated solutions for a rotational period below ten days, which corresponds to a system younger than about 100 Myr.

Principle
The initial step to use TATOO is first to model a very large range of generic planetary system evolutions (that will be called the initial grid in the following); then, from this initial grid, precompiled age exploration files are produced at specific pairs P rot,i ; P orb,i for each set of M ; M p masses explored in this work (that will be called the pre-compiled age exploration grid in the following). These specific pairs are regularly spaced in the P rot, ; P orb,p space. TATOO heavily relies on this precompiled exploration grid as it will be used to interpolate the age of a given planetary system by using its observed properties P rot,obs ; P orb,obs . I consider TATOO and PROBE as two different "entities" since TATOO can be used with any other grids as long as they follow the format of the current one.
In contrast to the tailored grid introduced above, the initial grid is designed to be as generic as possible. It is composed of stars between 0.5 and 1.0 M (with P rot,init between 1 and 11 days) and planetary mass from 0.5 to 3.5 M jup (with P orb,init between 0.1 and 1.0 P rot,init ). It is produced using the model presented in Gallet et al. (2018). The range of stellar mass is adopted because the rotational evolution of these stars is well modelled.  The initial rotation period P rot,init is chosen to cover the observed range of rotation in the early-PMS cluster (million-year-old clusters e.g. the Orion Nebulae Cluster and NGC6530). Finally the range of planetary mass between 0.5 to 3.5 is selected so as to ensure a good ratio between a strong enough tidal interaction and sufficiently numerous observed planetary systems in the literature. The upper planetary mass limit can be extended up to the brown dwarf limit (i.e.~13 M jup ) and the stellar mass range down to 0.3 and up to 1.2 M . However, the lower limit of 0.5 M jup is already a strong constraint. Indeed, a massive enough planet is needed in order to have a measurable impact on the rotation of the star through tidal interaction. Section 2.3 describes in more detail the limitations of the present tool ( Fig. A.1 summarizes how the initial grid is created). This grid cannot be directly used to extract the age of a given planetary system, except if the properties of this observed system already match one of the synthetic systems from the initial gird. In any other cases, an interpolation first needs to the realized.
From this initial grid, I then extract the ages of the systems that fulfil the requirement S 2 ≤ 100 (voluntarily larger than the actual maximum critical S 2 value, see Eq. 4). This extraction is done for each stellar mass at the specific pairs P rot,i ; P orb,i and the results are stored in pre-compiled age exploration files as a function of the planet's mass: the pre-compiled age exploration grid.
It is important to understand here that the initial grid and the pre-compiled age exploration grid are two different grids. The initial grid is composed of the temporal evolution of planetary systems, with given initial conditions, which are produced with the PROBE code. The pre-compiled age exploration grid is composed of the ages found from the initial grid at the specific pairs P rot,i ; P orb,i for a given M star, but for the whole range of planetary masses M p . The interpolation in planetary mass to extract the age of the actual observed system is done in the pre-compiled age exploration grid.
The relation between the orbital period and the semi-major axis is given by but I note that even if the planetary distances are expressed in P orb,p units, the model itself works with the semi-major axis, which is independent of the properties of the system. Indeed a given value of P orb,p will correspond to different physical distances (a) depending on the masses of the star and planet of the system. In this work P rot,i ranges from 2 to 30 days and P orb,i from 0.36-0.51 days (semi-major axis = 0.01 au) to 10.67-15.11 days (semi-major axis = 0.095 au) are considered. I explored this range of P rot,i for the exploration grid because most of the stars have a surface rotation period within 2-30 days between the PMS and the end of the MS phase (see Gallet & Bouvier 2013, 2015. For P orb,i , I investigated planets with a semi-major axis a from 0.095 au down to 0.01 au because below 0.01 au the planet is most probably engulfed by its host star, and beyond 0.05 au no tidal interaction and migration are expected given the planetary masses considered here. Finally, TATOO follows these steps: 1 1. It attains the observed periods and masses ( P rot,obs ; P orb,obs ; M p ; M ) together with error bars (σ P rot,obs and σ P orb,obs ) of the user's planetary system. A draw is then done on P rot,obs and P orb,obs given σ P rot,obs and σ P orb,obs . These random periods are used as the input P rot,obs and P orb,obs . 2. Then TATOO finds the four synthetic pre-compiled age exploration files (defined by the values of P rot,i and P orb,i ) that encompass the observed P rot,obs ; P orb,obs pair. 3. For each of these couples, linear interpolation on the age as a function of planetary mass is performed, which ends up with four numerical relations of the form: Age = θ 1 ×M p +θ 0 (with θ 1 and θ 0 two numerical constants) and four evaluations of the linearity via the Pearson coefficient Pe c . 4. TATOO then finds the age of these four couples at the observed planetary mass, which creates a 3D map composed of age; P rot,i ; P orb,i for the four nearest points. To assess the robustness of the final age estimation process, I introduced a quality number into the rest of this paper that depends on the value Pe c . 5. Finally, TATOO uses a 3D interpolation method (Python-SciPy griddata routine) to provide the age of the requested system with the observed P rot,obs ; P orb,obs ; M p ; M quadruplet. 6. Procedures 1 to 5 are repeated 100 times for the two stellar masses that border the observed stellar mass M and the final age is given using a linear interpolation between the two explored masses. The error on this final age estimate is given from the standard deviation of these 100 age evaluations. Figure 1 describes the four nearest points selection method. In this framework, the correct definition of the S 2 quantity is finally where σ P orb,p ∼ 10 −5 days (error on a: σ a = 0.00018 au) and σ P rot, = 0.4 days in this instance.
Given that the expression of S 2 above is expressed in "standard deviation units", we can directly use it to define a critical S 2 value. Since S 2 is here a sum, we can consider an "ellipse" of semi-major axis 3σ defined by 1 = P orb,p − P orb,i 3σ P orb,p 2 + P rot, − P rot,i 3σ P rot, 2 , which then produces a threshold S 2 lim = 9. In Table A.1 I explored the impact of the choice of S 2 lim on the age estimation.

Limitation
In TATOO the age of the system is interpolated from the planetary mass after the interpolations in stellar mass, stellar rotation period, and planetary orbital period. The main problem with this tool is the assumption that there is a linear correlation between the age of the system and the mass of the planet. To investigate the linearity of this relation, the Pearson correlation coefficient Pe c between the mass of the planet M p and the age of the system, for a given P rot,i ; P orb,i ; M triplet, is extracted from each synthetic pre-compiled age exploration file. The Pearson correlation coefficient is a quantity that estimates the linear correlation between two variables; the closer this value is to one, the higher the linear correlation between these two variables is. Figure 2 shows that most of the Pearson correlation coefficients, extracted from the whole pre-compiled exploration grid, are above 0.5 and close to 1. More than 64% of the Pearson coefficient estimates calculated from the exploration files are above the moderate positive limit of 0.5, and the median of these estimates is 0.6 (average = 0.58). This suggests an actual linear correlation between M p and the age of the system. To fully assess the validity of the tidal-chronology technique, age estimates provided by TATOO should be compared to fundamentally independent age predictions such as ages based on asteroseismic data (see Bellinger et al. 2019, and all the work of Dr Bellinger related to asteroseismic age). Unfortunately, this test will have to wait for the TESS (Transiting Exoplanet Survey Satellite) and PLATO (PLAnetary Transits and Oscillations of Stars). After that please use only the abbreviation. An abbreviation or acronym may be introduced if it is used more than three times in the Abstract and more than five times in the main text (otherwise spell out). missions, since for now only a few asteroseismic age estimations exist the host stars of massive close-in planets. Indeed, most of these estimates are for red giant branch stars or for Kepler targets that are composed of either a planet that is not sufficiently massive or one whose semi-major axis is too large. Moreover, most of the targets listed in Table 2 are located in the southern hemisphere, while the Kepler observations and the asteroseismic data are for targets in the northern hemisphere. Hence, PLATO and the future TESS observations are fundamental to constrain and test the tidal-chronology technique. It is also worth noting that sometimes only the minimum mass of the planet can be estimated, this is the case for instance in a radial velocity survey. In those cases, only a lower age limit could be estimated via TATOO.
Finally, it is important to be aware that even if the principle of tidal-chronology is a physical reality (see Gallet & Delorme 2019), the outcome of this technique is only valid in the framework of the considered hypothesis regarding the geometry of the system (coplanar and circular orbit) as well as in relation to how the tidal interaction is treated. Moreover, because of the nature of the initial grid delivered with this version, TATOO currently only works for planetary systems composed of a star between 0.5 and 1.0 M around which orbits a planet between 0.5 and 3.5 M jup .

Simple gyrochronological mode
I first investigate whether TATOO can be used as a simple gyrochronological tool. For this, the ages of the observed systems listed in Table 2 are extracted using the same procedure as above but by considering a higher P orb,obs for the planet so as to remove the impact of this latter on the surface rotation of the star. These ages, hereafter named tidal-gyrochronological ages, are thus the counterpart of the empirical gyrochronological ages but are estimated using the full Gallet et al. (2018) period and orbital evolution numerical model. These tidal-gyrochronological ages (Age tidal−gyro ) are listed in Table 1. For the gyrochronological ages, the tool and calibration provided by Angus et al. (2015) were used. They are calibrated using the Kepler asteroseismic targets with the (B-V)-mass relation that is given by the YREC (The Yale Rotating Stellar Evolution Code) isochrones extracted at 600 Myr for solar metallicity stars (Z = 0.01757, see An et al. 2007) 2 . It is worth noticing that gyrochronology analysis can only be applied to systems with an actual age older than about 100-200 Myr for a 1.0 M star and 500-600 Myr for a 0.6 M star (Barnes 2010;Delorme et al. 2011). The gyrochronological age estimates from Angus et al. (2015) are on average 10% smaller than the gyrochronological ages from Delorme et al. (2011) used in Gallet & Delorme (2019). Each gyrochronological age estimate is the median of 100 evaluations using a random draw in the range of the observed σ P rot, of the actual stellar rotation period as the input rotation period; the error is given by the standard deviation of these 100 age estimations. Notes. All of these systems are from Maxted et al. (2015) and Penev et al. (2018), except NGTS-10 that is from McCormac et al. (2020).
In Table 1, the absolute deviation between Age gyro and Age tidal−gyro is on average 15%. From these estimates, 34% have an absolute deviation smaller than 10%, 53% have an absolute deviation smaller than 15%, 70% have an absolute deviation smaller than 20%, and 85% have an absolute deviation smaller than 25%. The agreement between Age gyro and Age tidal−gyro is better by using the Delorme et al. (2011) calibration. In that case, 70% of the sample have an absolute deviation lower than 10%, and 92% lower than 25%. Table 1 shows that there is a good agreement between the two techniques, gyrochronology and tidal-gyrochronology, which is interesting since these two methods are fundamentally opposed. Indeed, gyrochronology is based on observations (empirical method) while tidalchronology is strongly based on a theoretical approach.
The apparent agreement between these two age estimate paradigms validates the theoretical age determination employed by TATOO and shows that the model correctly works. It is even more interesting as empirical relations, such as gyrochronology, should in the long term be replaced by a theoretical-based model in which the physics is understood.

Application to several known systems
In this section, TATOO is used to estimate the age of several massive close-in planetary systems. These systems mostly come from Maxted et al. (2015), in which they compared isochrones to gyrochronological ages for various star-planets systems, and from Penev et al. (2018). In Maxted et al. (2015), they already Fig. 2: Three-dimensional plot between P orb,i , P rot,i , and the Pearson coefficient Pe c . The two additional two-dimensional plots display the density map of P orb,i against Pe c (left) and P rot,i against Pe c (right).
highlighted the possible impact of tidal spin-up induced by the presence of a massive planet to explain the large discrepancies between isochrone and gyrochronological ages.
The properties of the planetary systems considered in this article are displayed in Table 2. This table also lists their estimated tidal-chronology and gyrochronology ages. Systems for which the tidal-age is greater than the gyrochronological age by 20% are considered to have been impacted by tidal spin-up. It seems that a bit less than half of these systems have a tidal age estimate close to the gyrochronological one. The quality factor of each age estimation is defined as the sum between a value that depends on the value of Pe c and the ratio between the age estimate and its standard deviation σ Age : The greatness of the age estimation increases with the value of the quality factor. A quality factor above 60 indicates an intermediate quality while a quality factor of 90 is already a good estimate. Table 3 shows the values of the precision of the age estimates extracted from TATOO as a function of the stellar and planetary mass, and on the input errors σ P rot,obs and σ P orb,obs . With an error below 3% on P rot,obs and P orb,obs the precision of TATOO is lower than 10%. For a given value of σ P orb,obs and σ P rot,obs , the precision of TATOO increases for increasing rotation period and orbital period. Table 3: Median precision of TATOO given the initial input conditions.

Discussion
The apparent existence of the linear relation between the age and the planetary mass of the system suggests that for some systems the heavier the planetary mass is, the older the system. It means that for two given planetary systems that are subject to tidal interaction, if they have the same set of observed properties then their age is proportional to the mass of their orbiting planet. Indeed, since the orbital evolution timescale T evolves as a 8 (see Bolmont & Mathis 2016), then for a given P orb,obs and stellar properties, if two massive planets are still observed orbiting close to their central star, then it means that the more massive planet has started its evolution at a greater P orb,init than the less massive one, and that therefore this system is the oldest one. For example, say that we want to get the age a system in which the star is a 0.7 M with a rotation of 16 days and that there is a planet of unknown mass orbiting this latter at 0.02 au. At such a small star-planet distance, we expect that the planet will be impacted by tidal interaction that will most probably induce inward migration. Since more massive planets evolve faster because the characteristic timescale T decreases for increasing planetary mass, if two systems have the same properties (same stellar mass M , same stellar rotation period P rot, , and same planetary orbital period P orb,p ) then the age of the system increases for increasing planetary mass. This effect is depicted in Fig. 3, which shows the estimated age of a system composed of a 0.7 M with a rotation period of 16 days as a function of the mass of the planet that is located at 0.02 au around the star. The colour gradient depicts the value of the S 2 .
In the sample explored in this article, half of the systems experienced tidal-spin up during their evolution. These systems are the ones, on average, with the highest ratio P rot, /P orb,p and R co /a, where R co is the corotation radius defined as which is the distance at which the orbital period of the planet equals the rotation period of the star. For the tidal spin-up cases the averaged P rot, /P orb,p and R co /a ratios are around 9.2 and 1.24, respectively. For the non tidal spin-up cases the averaged P rot, /P orb,p and R co /a ratios are around 6.0 and 0.98, respectively. A ratio of R co /a > 1.0 refers to a planet inside the corotation radius, which may indicate that this latter has had an impact on the surface rotation rate of its host's star. In all of the explored cases, P orb,p is smaller than P rot, . Figure 4 shows the deviation between the age estimated with tidal-chronology and the one using gyrochronology as a function of the ratio P rot, /P orb,p . The Fig. 4: Deviation between tidal-chronology age and gyrochronology age as a function of the ratio P rot, /P orb,p . The dashed line is a linear fit considering the full dataset. The dashdotted line is also a linear fit but omitting the WASP-4 and Qatar-1 systems from the dataset. The red dots are the systems for which tidal-spin up is suspected while the blue dots are systems without tidal spin-up. The dashed line corresponds to a linear fit to the full data while the dash-dotted line is a linear fit to the same data but without the WASP-4 and Qatar-1 systems that are located in the lower-right corner of the plot (with a ratio~16% and a deviation lower about 7 and 30%, respectively). The size of the points increases for increasing quality factor extracted from Table 2. quantity R 2 is the coefficient of determination defined by with y i the data, f i the linear fit to the data, andȳ the data's averaged value. A perfect linear fit gives R 2 = 1. This figure shows that the deviation, which is linked to tidal spin-up, scales linearly with the ratio P rot, /P orb,p : the higher this latter is, the stronger the deviation will be. The mass of the systems, both stellar and planetary, does not seems to affect this linearity. This behaviour can be explained by the fact that an increasing P rot, /P orb,p corresponds to an increasing R co /a quantity that is linked to systems in which the planet is located increasingly close inside the corotation radius. Therefore, a high P rot, /P orb,p ratio system is expected to have experienced stronger tidal interaction and thus more powerful tidal spin-up effect. From the linear fit in Fig. 4 we can infer a limit around P rot, /P orb,p = 7 (corresponding to ã 1.01 R co ) above which the deviation is larger than 20% and therefore tidal-chronology should be used. This trend could then suggest that either the orbital and rotational period of the WASP-be placed in the tidal spin-up region for which gyrochronological ages are biased. In the case of WASP-4, it seems that this system is highly sensitive to the S 2 lim value. Indeed, Table A.1 shows that the deviation for this system could reaches 50% with S 2 lim = 18. However, the Qatar-1 system appears to be stable regarding the choice of S 2 lim , which could point towards a incorrectly estimated rotational period (assuming that the orbital period is correct). Covino et al. (2013) reported a rotational period of 24±6 days for Qatar-1, which is very close to the 23.7±0.5 days found by Penev et al. (2018) but with a higher deviation. This produces a tidal-age of 3728 Myr corresponding to a deviation of about 40%. Given the orbital period of Qatar-1 b of 1.42 days, the ratio P rot, /P orb,p could range from 12.7 to 21, placing the Qatar-1 system in better agreement with the linear model.
Finally TATOO, like gyrochronology, is also not well suited for measuring the age of binary systems in which an external torque due to the stellar companion could have been applied to the surface rotation of the primary body. This is probably the case in this work for the TrES-2 system, which is known to be a binary system. For this system the tidal-gyrochronological age differs by -35.6% from the simple gyrochronological estimation. For this kind of system, the estimated age using tidal-chronology is most probably biased depending on the strength of the binary interaction. Removing this system from the sample increases the linearity of the fit up to R 2 = 0.5 and reduces the tidal spin-up limit down to 6.74 corresponding to a~0.98 R co .

Conclusion
While the age of stars is a fundamental parameter that can strongly constrain on-going stellar and planetary models, it is currently the most uncertain one. It is especially true in the case of stars that are or were subject to tidal spin-up induced by the presence of a massive close-in companion. In those specific cases, gyrochronology analysis cannot be used without having a bias towards a younger age estimate. This apparent young age will then lead to an incorrect estimation of the system properties.
To solve this problem, this paper provides the astrophysical community with a standalone tool, based on the initial work of Gallet & Delorme (2019), which can be used to rapidly and easily get an estimate of the age of a given massive close-in planetary system. This method only requires information on the stellar and planetary fundamental parameters: M , M p , P rot,obs , P orb,obs . This interconnectivity between stellar and planetary characteristics clearly highlights the need for the stellar and planetary communities to systematically provide all the available properties of the star and planet of their systems of interest.
While the most recent techniques only consider isolated stars (Angus et al. 2019;Bellinger 2019), TATOO is currently the only tool that proposes to estimate the age of massive closein planetary systems by taking into account the possible impact of the star-planet tidal interaction on the rotational evolution of the stellar surface. This tool can also be used even if no rotational departure is present. In that case, it gives results in agreement with the classical gyrochronological analysis. However, as with isochrone fitting, TATOO is model dependent and relies on the physics adopted in the numerical model used to produce the initial grid (here PROBE from Gallet et al. 2018). Hence, this model dependence should be taken into account when choosing to use TATOO as an age estimation tool. In this work I found that half of the investigated massive close-in systems for which TATOO can be applied, are subject to tidal spin-up and therefore have biased gyrochronological age estimations. I point out that this bias strongly scales with the ratio P rot, /P orb,p , which therefore could be used to easily predict the correctness of the gyrochronological age estimates. According to the investigated systems, a ratio of P rot, /P orb,p 7 is the indication that tidalchronology should be preferred to a simple gyrochronological tool. Indeed, above this limit the deviation between the two techniques starts to be greater than 20%. A list of such systems is given in Table 4. Moreover, correct ages can be retrieved, even without tidal spin-up effect, by using TATOO. In that case TATOO provides age estimates in good agreement with the gyrochronology technique. This proves that tidal-chronology can be seen as a first-order correction of the impact of tidal interaction on gyrochronology.
Finally, in the framework of the next generation of space missions such as PLATO (Rauer et al. 2014) and TESS (Ricker et al. 2015), TATOO will be a valuable age estimation tool as it will be crucial to have the best information on the planetary systems chosen as targets, including the estimation of their age, to fully interpret the data. The recent work of Amard & Matt (2020) also shows that metallicity could have a non-negligible effect for example on the validity of the current gyrochronology calibrations. This suggests the need for a more detailed analysis of the stability of tidal-chronology regarding the metallicity of the planetary systems investigated. This should be done in future works and is currently beyond the aim of this present one.