Transit Timing Variations Research Articles

The K2-24 planetary system revisited by CHEOPS

The planetary system K2-24 is composed of two transiting low-density Neptunians locked in an almost perfect 2:1 resonance and showing large transit time variations (TTVs), and it is an excellent laboratory to search for signatures of planetary migration. Previous studies performed with K2, Spitzer, and RV data tentatively claimed a significant non-zero eccentricity for one or both planets, possibly high enough to challenge the scenario of pure disk migration through resonant capture. With 13 new CHEOPS light curves (seven of planet b, six of planet c), we carried out a global photometric and dynamical re-analysis by including all the available literature data as well. We obtained the most accurate set of planetary parameters to date for the K2-24 system, including radii and masses at 1% and 5% precision (now essentially limited by the uncertainty on stellar parameters) and non-zero eccentricities eb = 0.0498−0.0018+0.0011, ec = 0.0282−0.0007+0.0003 detected at very high significance for both planets. Such relatively large values imply the need for an additional physical mechanism of eccentricity excitation during or after the migration stage. Also, while the accuracy of the previous TTV model had drifted by up to 0.5 days at the current time, we constrained the orbital solution firmly enough to predict the forthcoming transits for the next ~15 years, thus enabling efficient follow-up with top-level facilities such as JWST or ESPRESSO.

Probing the Possible Causes of the Transit Timing Variation for TrES-2b in the TESS Era

Nowadays, transit timing variations (TTVs) are proving to be a very valuable tool in exoplanetary science to detect exoplanets by observing variations in transit times. To study the TTV of the hot Jupiter TrES-2b, we have combined 64 high-quality transit light curves from all seven sectors of NASA's Transiting Exoplanet Survey Satellite along with 60 best-quality light curves from the ground-based facility Exoplanet Transit Database and 106 midtransit times from the previous works. From the precise transit timing analysis, we have observed a significant improvement in the orbital ephemerides, but we did not detect any short-period TTVs that might result from an additional body. The inability to detect short-term TTVs further motivates us to investigate long-term TTVs, which might be caused by orbital decay, apsidal precession, the Applegate mechanism, and the Rømer effect, and the orbital decay appeared to be a better explanation for the observed TTV with ΔBIC = 4.32. The orbital period of the hot Jupiter TrES-2b appears to be shrinking at a rate of ∼–5.58 ± 1.81 ms yr−1. Assuming this decay is primarily caused by tidal dissipation within the host star, we have subsequently calculated the stellar tidal quality factor value to be ∼9.9 × 103, which is 2–3 orders of magnitude smaller than the theoretically predicted values for other hot-Jupiter systems, and its low value indicates more efficient tidal dissipation within the host star. Additional precise photometric and radial velocity observations are required to pinpoint the cause of the change in the orbital period.

Transit timing variation of K2-237b: hints toward planet disc migration

ABSTRACT Hot Jupiters should initially form at considerable distances from host stars and subsequently migrate towards inner regions, supported directly by transit timing variation (TTV). We report the TTV of K2-237b, using reproduced timings fitted from Kepler K2 and TESS data. The timings span from 2016 to 2021, leading to an observational baseline of 5 yr. The timing evolution presents a significant bias to a constant period scenario. The model evidence is evaluated utilizing the Bayesian Information Criterion (BIC), which favours the scenario of period decay with a $\Delta$BIC of 14.1. The detected TTV induces a period decay rate ($\dot{P}$) of −1.14$\pm 0.28\times 10^{-8}$ d per day ($-$0.36 s/yr). Fitting the spectral energy distribution, we find infrared excess at the significance level of 1.5 $\sigma$ for WISE W1 and W2 bands, and 2 $\sigma$ level for W3 and W4 bands. This potentially reveals the existence of a stellar disc, consisting of hot dust at 800 $\pm$ 300 K, showing a $L_{\mathrm{ dust}}/L_{\ast }$ of 5 $\pm\ 3\times 10^{-3}$. We obtain a stellar age of 1.0$^{+1.4}_{-0.7}\, \times 10^{9}$ yr from isochrone fitting. The properties of K2-237b potentially serve as a direct observational support to the planet disc migration though more observation are needed.

Revisiting WASP-19b’s Orbital Evolution: New TESS Data and Apsidal Precession Analysis

Due to the ongoing debate of whether WASP19b is experiencing orbital decay or apsidal precession, I reexamine the orbital evolution of the hot Jupiter WASP-19b using new TESS data from sectors 09, 36, 62, and 63. My comprehensive analysis of transit and occultation timing variations compares models of constant period, orbital decay, and apsidal precession. Results show no significant preference between orbital decay and apsidal precession models. I find a period change rate of −1.40 ± 0.54 ms yr−1, slower than some previous studies but consistent with others. My apsidal precession analysis yields an eccentricity of 0.00178 and a precession rate of 0.036 deg day−1, differing from recent literature. These findings underscore the complexity of WASP-19b’s orbital dynamics and the need for continued long-term monitoring.

Characterisation of the warm-Jupiter TOI-1130 system with CHEOPS and a photo-dynamical approach

Context. Among the thousands of exoplanets discovered to date, approximately a few hundred gas giants on short-period orbits are classified as ‘lonely’ and only a few are in a multi-planet system with a smaller companion on a close orbit. The processes that formed multi-planet systems hosting gas giants on close orbits are poorly understood, and only a few examples of this kind of system have been observed and well characterised. Aims. Within the contest of a multi-planet system hosting a gas giant on short orbits, we characterise the TOI-1130 system by measuring masses and orbital parameters. This is a two-transiting planet system with a Jupiter-like planet (c) on a 8.35 days orbit and a Neptune-like planet (b) on an inner (4.07 days) orbit. Both planets show strong anti-correlated transit timing variations (TTVs). Furthermore, radial velocity (RV) analysis showed an additional linear trend, a possible hint of a non-transiting candidate planet on a far outer orbit. Methods. Since 2019, extensive transit and radial velocity observations of the TOI-1130 have been acquired using TESS and various ground-based facilities. We present a new photo-dynamical analysis of all available transit and RV data, with the addition of new CHEOPS and ASTEP+ data, which achieve the best precision to date on the planetary radii and masses and on the timings of each transit. Results. We were able to model interior structure of planet b constraining the presence of a gaseous envelope of H/He, while it was not possible to assess the possible water content. Furthermore, we analysed the resonant state of the two transiting planets, and we found that they lie just outside the resonant region. This could be the result of the tidal evolution that the system underwent. We obtained both masses of the planets with a precision of less than 1.5%, and radii with a precision of about 1% and 3% for planet b and c, respectively.

Precise Transit Photometry Using TESS. II. Revisiting 28 Additional Transiting Systems with Updated Physical Properties

Precise physical properties of the known transiting exoplanets are essential for their precise atmospheric characterization using modern and upcoming instruments. Leveraging the large volume of high-signal-to-noise-ratio photometric follow-up data from TESS, highly precise physical properties can be estimated for these systems, especially for those discovered using ground-based instruments prior to the TESS mission. In this work, I have used the publicly available TESS follow-up data for 28 transiting systems with 10 < V mag < 10.5, with an aim to update their known physical properties. The observed lightcurves have been analyzed by implementing a state-of-the-art critical noise treatment algorithm to effectively reduce both time-correlated and uncorrelated noise components, using sophisticated techniques like wavelet denoising and Gaussian-process regression. Compared with the previous studies, the estimated transit parameters are found to be more precise for most of the targets, including a few cases where a larger space-based instrument like Spitzer, Kepler, or CHEOPS has been used in the previous study. The large volume of transit observations used for each target has also resulted in a more accurate estimation of the physical properties, as this overcomes any error in parameter estimations from bias present in a smaller volume of data. Thus, comparing with the literature values, statistically significant improvements in the known physical properties of several targeted systems have been reported from this work. The large volume of transit-timing information from the analyses was also used to search for transit-timing variation trends in these targets, which has resulted in no significant detection.

Assessing Exoplanetary System Architectures with DYNAMITE Including Observational Upper Limits

The information gathered from observing planetary systems is not limited to the discovery of planets, but also includes the observational upper limits constraining the presence of any additional planets. Incorporating these upper limits into statistical analyses of individual systems can significantly improve our ability to find hidden planets in these systems by narrowing the parameter space in which to search. Here, I include radial velocity (RV), transit, and transit timing variation (TTV) upper limits on additional planets in known multiplanet systems into the Dynamite software package and test their impact on the predicted planets for these systems. The tests are run on systems with previous Dynamite analysis and with updated known planet parameters in the 2–3 yr since the original predictions. I find that the RV limits provide the strongest constraints on additional planets, lowering the likelihood of finding them within orbital periods of ∼10–100 days in the inner planetary systems, as well as truncating the likely planet size (radius and/or mass) distributions toward planets smaller than those currently observed. Transit and TTV limits also provide information on the size and inclination distributions of both the known and predicted planets in the system. Utilizing these limits on a wider range of systems in the near future will help determine which systems might be able to host temperate terrestrial planets and contribute to the search for extraterrestrial biosignatures.

TOI-1408: Discovery and Photodynamical Modeling of a Small Inner Companion to a Hot Jupiter Revealed by Transit Timing Variations

We report the discovery and characterization of a small planet, TOI-1408 c, on a 2.2 day orbit located interior to a previously known hot Jupiter, TOI-1408 b (P = 4.42 days, M = 1.86 ± 0.02 M Jup, R = 2.4 ± 0.5 R Jup) that exhibits grazing transits. The two planets are near 2:1 period commensurability, resulting in significant transit timing variations (TTVs) for both planets and transit duration variations for the inner planet. The TTV amplitude for TOI-1408 c is 15% of the planet’s orbital period, marking the largest TTV amplitude relative to the orbital period measured to date. Photodynamical modeling of ground-based radial velocity (RV) observations and transit light curves obtained with the Transiting Exoplanet Survey Satellite and ground-based facilities leads to an inner planet radius of 2.22 ± 0.06 R ⊕ and mass of 7.6 ± 0.2 M ⊕ that locates the planet into the sub-Neptune regime. The proximity to the 2:1 period commensurability leads to the libration of the resonant argument of the inner planet. The RV measurements support the existence of a third body with an orbital period of several thousand days. This discovery places the system among the rare systems featuring a hot Jupiter accompanied by an inner low-mass planet.

Photo-dynamical characterisation of the TOI-178 resonant chain

Context. The TOI-178 system consists of a nearby, late-K-dwarf with six transiting planets in the super-Earth to mini-Neptune regime, with radii ranging from to 2.9 R⊕ and orbital periods between 1.9 and 20.7 days. All the planets, but the innermost one, form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision event has taken place since the formation and migration of the planets in the protoplanetary disc, thereby providing important anchors for planet formation models. Aims. We aim to improve the characterisation of the architecture of this key system and, in particular, the masses and radii of its planets. In addition, since this system is one of the few resonant chains that can be characterised by both photometry and radial velocities, we propose to use it as a test bench for the robustness of the planetary mass determination with each technique. Methods. We performed a global analysis of all the available photometry from CHEOPS, TESS and NGTS, and radial velocity from ESPRESSO, using a photo-dynamical modelling of the light curve. We also tried different sets of priors on the masses and eccentricity, as well as different stellar activity models, to study their effects on the masses estimated by transit-timing variations (TTVs) and radial velocities (RVs). Results. We demonstrate how stellar activity prevents a robust mass estimation for the three outer planets using radial velocity data alone. We also show that our joint photo-dynamical and radial velocity analysis has resulted in a robust mass determination for planets c to 𝑔, with precision of ~ 12% for the mass of planet c, and better than 10% for planets d to 𝑔. The new precisions on the radii range from 2 to 3%. The understanding of this synergy between photometric and radial velocity measurements will be valuable for the PLATO mission. We also show that TOI-178 is indeed currently locked in the resonant configuration, librating around an equilibrium of the chain.

Repelling Planet Pairs by Ping-pong Scattering

The Kepler mission reveals a peculiar trough-peak feature in the orbital spacing of close-in planets near mean-motion resonances: a deficit and an excess that are, respectively, a couple of percent interior to and wide of the resonances. This feature has received two main classes of explanations: one involving eccentricity damping and the other scattering with small bodies. Here, we point out a few issues with the damping scenario and study the scattering scenario in more detail. We elucidate why scattering small bodies tends to repel two planets. As the small bodies random-walk in energy and angular momentum space, they tend to absorb fractionally more energy than angular momentum. This, which we call “ping-pong repulsion,” transports angular momentum from the inner to the outer planet and pushes the two planets apart. Such a process, even if ubiquitous, leaves identifiable marks only near first-order resonances: diverging pairs jump across the resonance quickly and produce the mean-motion resonance asymmetry. To explain the observed positions of the trough-peaks, a total scattering mass of order a few percent of the planet masses is required. Moreover, if this mass is dominated by a handful of Mercury-sized bodies, one can also explain the planet eccentricities as inferred from transit-time variations. Last, we suggest how these conditions may have naturally arisen during the late stages of planet formation.

TOI-2015 b: A Warm Neptune with Transit Timing Variations Orbiting an Active Mid-type M Dwarf

We report the discovery of a close-in (P orb = 3.349 days) warm Neptune with clear transit timing variations (TTVs) orbiting the nearby (d = 47.3 pc) active M4 star, TOI-2015. We characterize the planet's properties using Transiting Exoplanet Survey Satellite (TESS) photometry, precise near-infrared radial velocities (RVs) with the Habitable-zone Planet Finder Spectrograph, ground-based photometry, and high-contrast imaging. A joint photometry and RV fit yields a radius Rp=3.37−0.20+0.15R⊕ , mass mp=16.4−4.1+4.1M⊕ , and density ρp=2.32−0.37+0.38gcm−3 for TOI-2015 b, suggesting a likely volatile-rich planet. The young, active host star has a rotation period of P rot = 8.7 ± 0.9 days and associated rotation-based age estimate of 1.1 ± 0.1 Gyr. Though no other transiting planets are seen in the TESS data, the system shows clear TTVs of super-period Psup≈430days and amplitude ∼100 minutes. After considering multiple likely period-ratio models, we show an outer planet candidate near a 2:1 resonance can explain the observed TTVs while offering a dynamically stable solution. However, other possible two-planet solutions—including 3:2 and 4:3 resonances—cannot be conclusively excluded without further observations. Assuming a 2:1 resonance in the joint TTV-RV modeling suggests a mass of mb=13.3−4.5+4.7M⊕ for TOI-2015 b and mc=6.8−2.3+3.5M⊕ for the outer candidate. Additional transit and RV observations will be beneficial to explicitly identify the resonance and further characterize the properties of the system.

HD152843 b &amp; c: the masses and orbital periods of a sub-Neptune and a superpuff Neptune

ABSTRACT We present the characterization of the two transiting planets around HD 152843 (TOI 2319, TIC 349488688) using an intensive campaign of HARPS-N radial velocities, and two sectors of TESS data. These data reveal a unique and fascinating system: HD 152843 b and c have near equal masses of around 9 $M_{\hbox{$\oplus $}}$ but differing radii of $3.05 \pm 0.11$$R_{\hbox{$\oplus $}}$ and $5.94 _{ - 0.16 } ^ { + 0.18 }$$R_{\hbox{$\oplus $}}$, respectively, and orbital periods of $11.62071 _{ - 0.000106 } ^ { + 9.6e-05 }$ and $19.502104 _{ - 8.5e-05 } ^ { + 7.4e-05 }$ d. This indicates that HD 152843 c is in the lowest fifth-percentile in density of the known exoplanet population, and has the longest orbital period among these low-density planets. Further, HD 152843 c’s radius places it in the ‘Saturn valley’, the observed lack of planets larger than Neptune, but smaller than Saturn. The orbital periods of these planets indicate they are near a $5:3$ mean motion resonance, indicating the possibility of transit timing variations, and hints at the possibility of interaction with a third planet at some point in the evolution of this system. Further, the brightness of the host star and the low density of HD 152843 c make it a key target for atmospheric characterization.

TOI-2447 b / NGTS-29 b: a 69-day Saturn around a Solar analogue

ABSTRACT Discovering transiting exoplanets with relatively long orbital periods (&amp;gt;10 d) is crucial to facilitate the study of cool exoplanet atmospheres (Teq &amp;lt; 700 K) and to understand exoplanet formation and inward migration further out than typical transiting exoplanets. In order to discover these longer period transiting exoplanets, long-term photometric, and radial velocity campaigns are required. We report the discovery of TOI-2447 b (=NGTS-29 b), a Saturn-mass transiting exoplanet orbiting a bright (T = 10.0) Solar-type star (Teff = 5730 K). TOI-2447 b was identified as a transiting exoplanet candidate from a single transit event of 1.3 per cent depth and 7.29 h duration in TESS Sector 31 and a prior transit event from 2017 in NGTS data. Four further transit events were observed with NGTS photometry which revealed an orbital period of P = 69.34 d. The transit events establish a radius for TOI-2447 b of $0.865 \pm 0.010\, \rm R_{\rm J}$, while radial velocity measurements give a mass of $0.386 \pm 0.025\, \rm M_{\rm J}$. The equilibrium temperature of the planet is 414 K, making it much cooler than the majority of TESS planet discoveries. We also detect a transit signal in NGTS data not caused by TOI-2447 b, along with transit timing variations and evidence for a ∼150 d signal in radial velocity measurements. It is likely that the system hosts additional planets, but further photometry and radial velocity campaigns will be needed to determine their parameters with confidence. TOI-2447 b/NGTS-29 b joins a small but growing population of cool giants that will provide crucial insights into giant planet composition and formation mechanisms.

New ephemerides and detection of transit-timing variations in the K2-138 system using high-precision CHEOPS photometry

Multi-planet systems are a perfect laboratory for constraining planetary formation models. A few of these systems present planets that come very close to mean motion resonance, potentially leading to significant transit-timing variations (TTVs) due to their gravitational interactions. Of these systems represents a excellent laboratory for studying the dynamics of its six small planets (with radii ranging between $ as the five innermost planets are in a near 3:2 resonant chain. In this work, we aim to constrain the orbital properties of the six planets in the system by monitoring their transits with CHaracterising ExOPlanets Satellite ( We also seek to use this new data to lead a TTV study on this system. We obtained twelve light curves of the system with transits of planets $d$, $e$, $f,$ and $g$. With these data, we were able to update the ephemerides of the transits for these planets and search for timing transit variations. With our measurements, we reduced the uncertainties in the orbital periods of the studied planets, typically by an order of magnitude. This allowed us to correct for large deviations, on the order of hours, in the transit times predicted by previous studies. This is key to enabling future reliable observations of the planetary transits in the system. We also highlight the presence of potential TTVs ranging from 10 minutes to as many as 60 minutes for planet $d$.

HD 110067 c has an aligned orbit

Planetary systems in mean motion resonances hold a special place among the planetary population. They allow us to study planet formation in great detail as dissipative processes are thought to have played an important role in their existence. Additionally, planetary masses in bright resonant systems can be independently measured via both radial velocities and transit timing variations. In principle, they also allow us to quickly determine the inclination of all planets in the system since, for the system to be stable, they are likely all in coplanar orbits. To describe the full dynamical state of the system, we also need the stellar obliquity, which provides the orbital alignment of a planet with respect to the spin of its host star and can be measured thanks to the Rossiter–McLaughlin effect. It was recently discovered that HD 110067 harbors a system of six sub-Neptunes in resonant chain orbits. We here analyze an ESPRESSO high-resolution spectroscopic time series of HD 110067 during the transit of planet c. We find the orbit of HD 110067 c to be well aligned, with a sky-projected obliquity of λ =6+24-26 deg. This result indicates that the current architecture of the system was reached through convergent migration without any major disruptive events. Finally, we report transit-timing variation in this system as we find a significant offset of 19 ± 4 min in the center of the transit compared to the published ephemeris.

Updated Catalog of Kepler Planet Candidates: Focus on Accuracy and Orbital Periods

We present a new catalog of Kepler planet candidates that prioritizes accuracy of planetary dispositions and properties over uniformity. This catalog contains 4376 transiting planet candidates, including 1791 residing within 709 multiplanet systems, and provides the best parameters available for a large sample of Kepler planet candidates. We also provide a second set of stellar and planetary properties for transiting candidates that are uniformly derived for use in occurrence rate studies. Estimates of orbital periods have been improved, but as in previous catalogs, our tabulated values for period uncertainties do not fully account for transit timing variations (TTVs). We show that many planets are likely to have TTVs with long periodicities caused by various processes, including orbital precession, and that such TTVs imply that ephemerides of Kepler planets are not as accurate on multidecadal timescales as predicted by the small formal errors (typically 1 part in 106 and rarely >10−5) in the planets’ measured mean orbital periods during the Kepler epoch. Analysis of normalized transit durations implies that eccentricities of planets are anticorrelated with the number of companion transiting planets. Our primary catalog lists all known Kepler planet candidates that orbit and transit only one star; for completeness, we also provide an abbreviated listing of the properties of the two dozen nontransiting planets that have been identified around stars that host transiting planets discovered by Kepler.

Characterisation of the TOI-421 planetary system using CHEOPS, TESS, and archival radial velocity data

Context. The TOI-421 planetary system contains two sub-Neptune-type planets (Pb ~ 5.2 days, Teq,b ~ 900 K, and Pc ~ 16.1 days, Teq,c ~ 650 K) and is a prime target to study the formation and evolution of planets and their atmospheres. The inner planet is especially interesting as the existence of a hydrogen-dominated atmosphere at its orbital separation cannot be explained by current formation models without previous orbital migration. Aims. We aim to improve the system parameters to further use them to model the interior structure and simulate the atmospheric evolution of both planets, to finally gain insights into their formation and evolution. We also investigate the possibility of detecting transit timing variations (TTVs). Methods. We jointly analysed photometric data of three TESS sectors and six CHEOPS visits as well as 156 radial velocity data points to retrieve improved planetary parameters. We also searched for TTVs and modelled the interior structure of the planets. Finally, we simulated the evolution of the primordial H-He atmospheres of the planets using two different modelling frameworks. Results. We determine the planetary radii and masses of TOI-421 b and c to be Rb = 2.64 ± 0.08 R⊕, Mb = 6.7 ± 0.6 M⊕, Rc = 5.09 ± 0.07 R⊕, and Mc = 14.1 ± 1.4 M⊕. Using these results we retrieved average planetary densities of ρb = 0.37 ± 0.05ρ⊕ and ρc = 0.107 ± 0.012 ρ⊕. We do not detect any statistically significant TTV signals. Assuming the presence of a hydrogen-dominated atmosphere, the interior structure modelling results in both planets having extensive envelopes. While the modelling of the atmospheric evolution predicts for TOI-421 b to have lost any primordial atmosphere that it could have accreted at its current orbital position, TOI-421 c could have started out with an initial atmospheric mass fraction somewhere between 10 and 35%. Conclusions. We conclude that the low observed mean density of TOI-421 b can only be explained by either a bias in the measured planetary parameters (e.g. driven by high-altitude clouds) and/or in the context of orbital migration. We also find that the results of atmospheric evolution models are strongly dependent on the employed planetary structure model.

The Six-planet Resonant Chain of HD 110067

HD 110067 is the brightest star known to have six transiting planets. Each adjacent pair of planets has a period ratio that is nearly equal to a ratio of small integers, suggesting the planets are in a chain of mean-motion resonances, but the limited time span of the available data has prevented firm conclusions. Here, we show that the requirement of long-term dynamical stability implies that all six planets are very likely to form a resonant chain. Dynamical simulations of nonresonant systems with initial conditions compatible with the available data almost always suffer an instability within 25 Myr (∼0.3% of the system’s age). Assuming the system is in resonance, we place upper limits on the planets’ eccentricities and lower limits on the masses of the planets that have not yet been measured. We also predict the characteristics of transit timing variations and the values of the three-body libration centers.

Frequency-Dependent Variability of Pulse Wave Transit Time: Pilot Study.

The dynamics of the pulse wave (PW) associated with the PW transit time variability (PWTTV) determines the peripheral pulse rate variability, which is used as a surrogate for heart rate variability (HRV). The aim of the work is to analyze the frequency-dependent dynamics of PWTTV and to identify the possible frequency-phase modulation of PW velocity oscillations on the transit from the heart to the soft tissues of the distal parts of the upper extremities. RR-interval recordings and synchronous records of photoplethysmograms of 12 conditionally healthy subjects from the PhysioNet open database were used in this work. Using the Hilbert-Huang transform 3 spectral components of PWTTV and HRV were identified. It was shown that the amplitudes of PWTTV oscillations were many times (up to 8.4 times) smaller than the amplitudes of HRV, and the peaks of PWTTV spectral components were shifted towards higher frequencies than those of HRV. Functional relations between PWTTV and HRV, which can determine the phase modulation of periodic changes in the PW propagation velocity, were revealed.

Star-spot activity, orbital obliquity, transmission spectrum, physical properties, and transit time variations of the HATS-2 planetary system

Aims. Our aim in this paper is to refine the orbital and physical parameters of the HATS-2 planetary system and study transit timing variations and atmospheric composition thanks to transit observations that span more than 10 yr and that were collected using different instruments and pass-band filters. We also investigate the orbital alignment of the system by studying the anomalies in the transit light curves induced by starspots on the photosphere of the parent star. Methods. We analysed new transit events from both ground-based telescopes and NASA’s TESS mission. Anomalies were detected in most of the light curves and modelled as starspots occulted by the planet during transit events. We fitted the clean and symmetric light curves with the JKTEBOP code and those affected by anomalies with the PRISM+GEMC codes to simultaneously model the photometric parameters of the transits and the position, size, and contrast of each starspot. Results. We found consistency between the values we found for the physical and orbital parameters and those from the discovery paper and ATLAS9 stellar atmospherical models. We identified different sets of consecutive starspot-crossing events that temporally occurred in less than five days. Under the hypothesis that we are dealing with the same starspots, occulted twice by the planet during two consecutive transits, we estimated the rotational period of the parent star and, in turn the projected and the true orbital obliquity of the planet. We find that the system is well aligned. We identified the possible presence of transit timing variations in the system, which can be caused by tidal orbital decay, and we derived a low-resolution transmission spectrum.