Mid-transit Times Research Articles

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.

The 𝒯ℛ𝒪𝒴 project

Context. Co-orbital objects, also known as trojans, are frequently found in simulations of planetary system formation. In these configurations, a planet shares its orbit with other massive bodies. It is still unclear why there have not been any co-orbitals discovered thus far in exoplanetary systems (exotrojans) or even pairs of planets found in such a 1:1 mean motion resonance. Reconciling observations and theory is an open subject in the field. Aims. The main objective of the 𝒯ℛ𝒪𝒴 project is to conduct an exhaustive search for exotrojans using diverse observational techniques. In this work, we analyze the radial velocity time series informed by transits, focusing the search around low-mass stars. Methods. We employed the α-test method on confirmed planets searching for shifts between spectral and photometric mid-transit times. This technique is sensitive to mass imbalances within the planetary orbit, allowing us to identify non-negligible co-orbital masses. Results. Among the 95 transiting planets examined, we find one robust exotrojan candidate with a significant 3-σ detection. Additionally, 25 exoplanets show compatibility with the presence of exotrojan companions at a 1-σ level, requiring further observations to better constrain their presence. For two of those weak candidates, we find dimmings in their light curves within the predicted Lagrangian region. We established upper limits on the co-orbital masses for either the candidates and null detections. Conclusions. Our analysis reveals that current high-resolution spectrographs effectively rule out co-orbitals more massive than Saturn around low-mass stars. This work points out to dozens of targets that have the potential to better constraint their exotrojan upper mass limit with dedicated radial velocity observations. We also explored the potential of observing the secondary eclipses of the confirmed exoplanets in our sample to enhance the exotrojan search, ultimately leading to a more accurate estimation of the occurrence rate of exotrojans.

Enhancing Exoplanet Ephemerides by Leveraging Professional and Citizen Science Data: A Test Case with WASP-77 A b

We present an updated ephemeris, and physical parameters, for the exoplanet WASP-77 A b. In this effort, we combine 64 ground- and space-based transit observations, 6 space-based eclipse observations, and 32 radial velocity observations to produce this target's most precise orbital solution to date aiding in the planning of James Webb Space Telescope and Ariel observations and atmospheric studies. We report a new orbital period of 1.360029395 ± 5.7 × 10−8 days, a new mid-transit time of 2459957.337860 ± 4.3 × 10−5 Barycentric Julian Date in the Barycentric Dynamical Timescale (BJDTDB) and a new mid-eclipse time of 2459956.658192 ± 6.7 × 10−5 BJDTDB. Furthermore, the methods presented in this study reduce the uncertainties in the planet's mass 1.6654 ± 4.5 × 10−3 M Jup and orbital period 1.360029395 ± 5.7 × 10−8 days by factors of 15.1 and 10.9, respectively. Through a joint fit analysis comparison of transit data taken by space-based and citizen science-led initiatives, our study demonstrates the power of including data collected by citizen scientists compared to a fit of the space-based data alone. Additionally, by including a vast array of citizen science data from ExoClock, Exoplanet Transit Database, and Exoplanet Watch, we can increase our observational baseline and thus acquire better constraints on the forward propagation of our ephemeris than what is achievable with Transiting Exoplanet Survey Satellite data alone.

Looking for timing variations in the transits of 16 exoplanets

ABSTRACT We update the ephemerides of 16 transiting exoplanets using our ground-based observations, new Transiting Exoplanet Survey Satellite data, and previously published observations including those of amateur astronomers. All these light curves were modelled by making use of a set of quantitative criteria with the exofast code to obtain mid-transit times. We searched for statistically significant secular and/or periodic trends in the mid-transit times. We found that the timing data are well modelled by a linear ephemeris for all systems except for XO-2 b, for which we detect an orbital decay with the rate of −12.95 ± 1.85 ms yr−1 that can be confirmed with future observations. We also detect a hint of potential periodic variations in the transit timing variation data of HAT-P-13 b, which also requires confirmation with further precise observations.

TOI-6883.01: A Single-transit Planet Candidate Detected from TESS

A new candidate exoplanet, proposed in ExoFOP by authors, it was promoted from TIC 393818343 to TOI-6883.01 at coordinates R.A.(J2000)20:41:10.01 decl.(J2000) + 3:38:17.87 in the Delphinus constellation and its distance is (93.73 ± 0.35) pc from Earth. The target star is a Sun-type having G0 class according to Skiff spectral classification and it has photometric magnitude V(Johnson) = 9.5 mag. We analyzed the transit light curve using Transit Exoplanet Survey Satellite’s ExoMAST tool, and evaluated the Lightcurve Analysis Tool for Transiting Exoplanets-report by excluding false positives, studying the background, satellite and centroid motion, and analyzing the pixels in sector 55 to define the source goodness. The candidate exoplanet transit has a duration of TD = 3.957 hr, Depth = 0.1196% and mid-transit time occurs at TBJD = 2811.24. Only one event was observed, so the orbital period cannot be evaluated.

The Advantages of Global Photometric Models in Fitting Transit Variations

Estimation of planetary orbital and physical parameters from light-curve data relies heavily on the accurate interpretation of Transit Timing Variation (TTV) measurements. In this letter, we review the process of TTV measurement and compare two fitting paradigms—one that relies on making transit-by-transit timing estimates and then fitting a TTV model to the observed timings, and one that relies on fitting a global flux model to the entire light-curve data set simultaneously. The latter method is achieved either by solving for the underlying planetary motion (often referred to as “photodynamics”), or by using an approximate or empirical shape of the TTV signal. We show that, across a large range of the transit S/N regime, the probability distribution function of the mid-transit time significantly deviates from a Gaussian, even if the flux errors do distribute normally. Treating the timing uncertainties as if they are distributed normally leads, in such a case, to a wrong interpretation of the TTV measurements. We illustrate these points using numerical experiments and conclude that a fitting process that relies on a global flux fitting, rather than the derived TTVs, should be preferred.

Search for the wide-orbit massive companion of XO-7b in the follow-up radial-velocity and transit-timing data: no significant clues

ABSTRACT XO-7b is a hot Jupiter transiting a V = 10.52 mag G0V-type star. The planetary system is interesting because the linear slope in the discovery radial-velocity (RV) data indicated a wide-orbit massive companion. In 2020 we started an RV campaign for the system with the main scientific goal to follow-up this linear slope, and to put constraints on the orbital period of the companion. Furthermore, we aimed at refining the system parameters and we wanted to probe transit time variations (TTVs) of XO-7b in order to search for long-term dynamical signs of the companion of XO-7b in the observed-minus-calculated (O-C) data of mid-transit times. Apart from the discovery RVs, we obtained and analysed 20 follow-up RV observations and TESS photometric data. The previously observed significant linear RV slope was not confirmed with the follow-up RV data, where we detected only a marginal linear slope with the opposite trend. If the announced companion really exists, the most convincing explanation is that both RV data sets were collected near its quadrature position. Based on the RVs we estimated the minimum orbital period, which is Porb, min, 3 ≳ 7900 ± 1660 d, and the ‘minimum’ minimum mass of the companion, which is (M3sin i)min = 16.7 ± 3.5 MJup. We did not find significant evidence of the companion of XO-7b in the O-C data set of mid-transit times. We can again conclude that if the announced companion really exists, this is in agreement with previous results that distant companions of exoplanets are only known by RV solutions.

Investigation on transit observations of the WASP-10 and WASP-11 systems

In the study, we present the results of new transit observations of two hot Jupiter WASP-10 b and WASP-11 b using the T100 telescope in Turkey. By analyzing the new transit light curves with the MCMC method through EXOFASTv2, we updated the physical parameters of two hot-Jupiter systems. The updated parameters are in agreement with the previously determined values within the uncertainties except eccentricity and impact factor. Starspot activity, noticed in previous studies, was not clearly detectable due to our highly scattering transit light curves of both systems. We determined an updated linear ephemeris and orbital period for the two systems by combining our new mid-transit times with the ones in the literature. There is no significant evidence for the orbital period variation of the systems in the most-updated O−C diagram. In addition, we could not find any strong signals for both two systems in the period analysis.

Searching for candidates of orbital decays among transit exoplanets

Transit observations have become an important technique to probe exoplanets. Therefore, there are many projects carrying on organized observations of transit events, which make a huge amount of light-curve and transit timing data available. We consider this as an excellent opportunity to search for possible orbital decays of exoplanets from this big number of mid-transit times through data-model fitting with both fixed-orbit and orbit-decay models. In order to perform this task, we collect mid-transit-time data from several sources and construct the most complete database up to date. Among 144 hot Jupiters in our study, HAT-P-51b, HAT-P-53b, TrES-5b, WASP-12b are classified as the orbit-decay cases. Thus, in addition to reconfirming WASP-12b as an orbit-decay planet, our results indicate that HAT-P-51b, HAT-P-53b, TrES-5b are potential orbit-decay candidates.

Investigating the Variation of Selected Kepler Objects Mid-Transit Times

We derived minima times from the transit curves of star-planet systems, Kepler-412, Kepler-422, Kepler-427 ve Kepler-435 observed by Kepler space telescope, using the Kwee - van Woerden method and Gauss function fitting. We examined the O-C diagram of each system separately, modelled them with linear and quadratic functions. We obtained that linear models give best fit for O-C distributions. We presented updated light elements of systems and concluded that O-C diagrams of systems can best be represented by the linear model.

Discovering planets with PLATO: Comparison of algorithms for stellar activity filtering

Context. To date, stellar activity is one of the main limitations in detecting small exoplanets via the transit photometry technique. Since this activity is enhanced in young stars, traditional filtering algorithms may severely underperform in attempting to detect such exoplanets, with shallow transits often obscured by the photometric modulation of the light curve. Aims. This paper aims to compare the relative performances of four algorithms developed by independent research groups specifically for the filtering of activity in the light curves of young active stars, prior to the search for planetary transit signals: Notch and LOCoR (N&L), Young Stars Detrending (YSD), K2 Systematics Correction (K2SC), and VARLET. Our comparison also includes the two best-performing algorithms implemented in the Wōtan package: Tukey’s biweight and Huber spline algorithms. Methods. For this purpose, we performed a series of injection-retrieval tests of planetary transits of different types, from Jupiter down to Earth-sized planets, moving both on circular and eccentric orbits. These experiments were carried out over a set of 100 realistically simulated light curves of both quiet and active solar-like stars (i.e., F and G types) that will be observed by the ESA Planetary Transits and Oscillations of stars (PLATO) space telescope, starting 2026. Results. From the experiments for transit detections, we found that N&L is the best choice in many cases, since it misses the lowest number of transits. However, this algorithm is shown to underperform when the planetary orbital period closely matches the stellar rotation period, especially in the case of small planets for which the biweight and VARLET algorithms work better. Moreover, for light curves with a large number of data-points, the combined results of two algorithms, YSD and Huber spline, yield the highest recovery percentage. Filtering algorithms allow us to obtain a very precise estimate of the orbital period and the mid-transit time of the detected planets, while the planet-to-star radius is underestimated most of the time, especially in cases of grazing transits or eccentric orbits. A refined filtering that takes into account the presence of the planet is thus compulsory for proper planetary characterization analyses.

A super-Earth and a mini-Neptune near the 2:1 MMR straddling the radius valley around the nearby mid-M dwarf TOI-2096

Context.Several planetary formation models have been proposed to explain the observed abundance and variety of compositions of super-Earths and mini-Neptunes. In this context, multitransiting systems orbiting low-mass stars whose planets are close to the radius valley are benchmark systems, which help to elucidate which formation model dominates.Aims.We report the discovery, validation, and initial characterization of one such system, TOI-2096 (TIC 142748283), a two-planet system composed of a super-Earth and a mini-Neptune hosted by a mid-type M dwarf located 48 pc away.Methods.We characterized the host star by combining optical spectra, analyzing its broadband spectral energy distribution, and using evolutionary models for low-mass stars. Then, we derived the planetary properties by modeling the photometric data from TESS and ground-based facilities. In addition, we used archival data, high-resolution imaging, and statistical validation to support our planetary interpretation.Results.We found that the stellar properties of TOI-2096 correspond to a dwarf star of spectral type M4±0.5. It harbors a super-Earth (R= 1.24 ± 0.07R⊕) and a mini-Neptune (R= 1.90 ± 0.09R⊕) in likely slightly eccentric orbits with orbital periods of 3.12 d and 6.39 d, respectively. These orbital periods are close to the first-order 2:1 mean-motion resonance (MMR), a configuration that may lead to measurable transit timing variations (TTVs). We computed the expected TTVs amplitude for each planet and found that they might be measurable with high-precision photometry delivering mid-transit times with accuracies of ≲2 min. Moreover, we conclude that measuring the planetary masses via radial velocities (RVs) could also be possible. Lastly, we found that these planets are among the best in their class to conduct atmospheric studies using the NIRSpec/Prism onboard theJames WebbSpace Telescope (JWST).Conclusions.The properties of this system make it a suitable candidate for further studies, particularly for mass determination using RVs and/or TTVs, decreasing the scarcity of systems that can be used to test planetary formation models around low-mass stars.

Radial velocity confirmation of a hot super-Neptune discovered by TESS with a warm Saturn–mass companion

ABSTRACT We report the discovery and confirmation of the planetary system TOI-1288. This late G dwarf harbours two planets: TOI-1288 b and TOI-1288 c. We combine TESS space-borne and ground-based transit photometry with HARPS-N and HIRES high-precision Doppler measurements, which we use to constrain the masses of both planets in the system and the radius of planet b. TOI-1288 b has a period of $2.699835^{+0.000004}_{-0.000003}$ d, a radius of 5.24 ± 0.09 R⊕, and a mass of 42 ± 3 M⊕, making this planet a hot transiting super-Neptune situated right in the Neptunian desert. This desert refers to a paucity of Neptune-sized planets on short period orbits. Our 2.4-yr-long Doppler monitoring of TOI-1288 revealed the presence of a Saturn–mass planet on a moderately eccentric orbit ($0.13^{+0.07}_{-0.09}$) with a minimum mass of 84 ± 7 M⊕ and a period of $443^{+11}_{-13}$ d. The five sectors worth of TESS data do not cover our expected mid-transit time for TOI-1288 c, and we do not detect a transit for this planet in these sectors.

A 16 hr Transit of Kepler-167 e Observed by the Ground-based Unistellar Telescope Network

More than 5000 exoplanets have been confirmed and among them almost 4000 were discovered by the transit method. However, few transiting exoplanets have an orbital period greater than 100 days. Here we report a transit detection of Kepler-167 e, a “Jupiter analog” exoplanet orbiting a K4 star with a period of 1071 days, using the Unistellar ground-based telescope network. From 2021 November 18 to 20, citizen astronomers located in nine different countries gathered 43 observations, covering the 16 hr long transit. Using a nested sampling approach to combine and fit the observations, we detected the midtransit time to be UTC 2021 November 19 17:20:51 with a 1σ uncertainty of 9.8 minutes, making it the longest-period planet to ever have its transit detected from the ground. This is the fourth transit detection of Kepler-167 e, but the first made from the ground. This timing measurement refines the orbit and keeps the ephemeris up to date without requiring space telescopes. Observations like this demonstrate the capabilities of coordinated networks of small telescopes to identify and characterize planets with long orbital periods.

Revisiting the Transit Timing Variations in the TrES-3 and Qatar-1 Systems with TESS Data

We present and analyze 58 transit light curves of TrES-3b and 98 transit light curves of Qatar-1b, observed by the Transiting Exoplanet Survey Satellite, plus two transit light curves of Qatar-1b, observed by us, using a ground-based 1.23 m telescope. These light curves are combined with the best-quality light curves taken from the Exoplanet Transit Database and the literature. The precisely determined midtransit times from these light curves enable us to obtain the refined orbital ephemerides, with improved precision, for both hot Jupiters. From the timing analysis, we find indications of the presence of transit timing variations (TTVs) in both systems. Since the observed TTVs are unlikely to be short-term and periodic, the possibility of additional planets in orbits close to TrES-3b and Qatar-1b is ruled out. The possible causes of long-term TTVs, such as orbital decay, apsidal precession, the Applegate mechanism, and line-of-sight acceleration, are also examined. However, none of these possibilities are found to explain the observed TTV of TrES-3b. In contrast to this, line-of-sight acceleration appears to be a plausible explanation for the observed TTV of Qatar-1b. In order to confirm these findings, further high-precision transit and radial velocity observations of both systems would be worthwhile.

The study on transmission spectrum and TTV behaviour of the hot Jupiter WASP-12b

ABSTRACT By utilizing the 1m and 2.4m telescopes of Yunnan Observatories, we obtained 13 new photometric transit light curves of WASP-12 exoplanetary system from 2009 to 2016. The data were analysed by using the Transit Analysis Package with the aid of Markov chain Monte Carlo technique. Based on the physical parameters from previous work, we have obtained the variation of radius with wavelength and mid-transit times of WASP-12b. Compared with a set of theoretical transmission spectra, its radius variation is almost consistent with a flatten and featureless model along with a slope in the bluer band, except the weak Na signal. The flatten and featureless transmission spectrum implies a possible evidence for the absence of TiO and VO in the upper atmosphere of WASP-12b. The slope in the bluer band may result from the Rayleigh scattering caused by H2/He and there would be hazes or hydride in the upper atmosphere. Combined new mid-transit times with all of the available timing data in literature, we have found that there is a significant change in the planetary orbit, the residuals of the timing prefer the orbital decay model to the apsidal precession model and the planetary interaction-induced TTV model. We updated the decay rate as −37.14 ± 1.31 ms yr−1 and derived a slightly lower modified quality factor $Q^{^{\prime }}_{*}$ for the constant-phase lag model. Furthermore, the current observations cannot rule out the planetary interaction-induced TTV model conclusively.

Variable and Supersonic Winds in the Atmosphere of an Ultrahot Giant Planet

Abstract Hot Jupiters (HJs) receive intense irradiation from their stellar hosts. The resulting extreme environments in their atmospheres allow us to study the conditions that drive planetary atmospheric dynamics, e.g., global-scale winds. General circulation models predict day-to-nightside winds and equatorial jets with speeds of the order of a few km s−1. To test these models, we apply high-resolution transmission spectroscopy using the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI) spectrograph on the Large Binocular Telescope to study the atmosphere of KELT-9 b, an ultrahot Jupiter and currently the hottest known planet. We measure ∼10 km s−1 day-to-nightside winds traced by Fe ii features in the planet’s atmosphere. This is at odds with previous literature (including data taken with PEPSI), which report no significant day-to-nightside winds on KELT-9 b. We identify the cause of this discrepancy as due to an inaccurate ephemeris for KELT-9 b in previous literature. We update the ephemeris, which shifts the midtransit time by up to 10 minutes for previous data sets, resulting in consistent detections of blueshifts in all the data sets analyzed here. Furthermore, a comparison with archival data sets from the High-accuracy Radial velocity Planet Searcher for the Northern hemisphere suggests a temporal wind variability of ∼5–8 km s−1 over timescales between weeks to years. Temporal variability of atmospheric dynamics on HJs is a phenomenon anticipated by certain general circulation models that has not been observed over these timescales until now. However, such large variability as we measure on KELT-9 b challenges general circulation models, which predict much lower amplitudes of wind variability over timescales between days to weeks.

Homogeneous transit timing analyses of 10 exoplanet systems

ABSTRACT We study the transit timings of 10 exoplanets in order to investigate potential transit timing variations in them. We model their available ground-based light curves, some presented here and others taken from the literature, and homogeneously measure the mid-transit times. We statistically compare our results with published values and find that the measurement errors agree. However, in terms of recovering the possible frequencies, homogeneous sets can be found to be more useful, of which no statistically relevant example has been found for the planets in our study. We corrected the ephemeris information of all 10 planets we studied and provide these most precise light elements as references for future transit observations with space-borne and ground-based instruments. We found no evidence for secular or periodic changes in the orbital periods of the planets in our sample, including the ultra-short period WASP-103 b, whose orbit is expected to decay on an observable time-scale. Therefore, we derive the lower limits for the reduced tidal quality factors (Q$^{\prime }_{\star }$) for the host stars based on best-fitting quadratic functions to their timing data. We also present a global model of all available data for WASP-74 b, which has a Gaia parallax-based distance value ∼25 per cent larger than the published value.

Nauyaca: a New Tool to Determine Planetary Masses and Orbital Elements through Transit Timing Analysis

Abstract The transit timing variations method is currently the most successful method to determine dynamical masses and orbital elements for Earth-sized transiting planets. Precise mass determination is fundamental to restrict planetary densities and thus infer planetary compositions. In this work, we present Nauyaca, a Python package dedicated to finding planetary masses and orbital elements through the fitting of observed midtransit times from an N-body approach. The fitting strategy consists of performing a sequence of minimization algorithms (optimizers) that are used to identify high probability regions in the parameter space. These results from optimizers are used for initialization of a Markov chain Monte Carlo method, using an adaptive Parallel-Tempering algorithm. A set of runs are performed in order to obtain posterior distributions of planetary masses and orbital elements. In order to test the tool, we created a mock catalog of synthetic planetary systems with different numbers of planets where all of them transit. We calculate their midtransit times to give them as an input to Nauyaca, testing statistically its efficiency in recovering the planetary parameters from the catalog. For the recovered planets, we find typical dispersions around the real values of ∼1–14 M ⊕ for masses, between 10–110 s for periods, and between ∼0.01–0.03 for eccentricities. We also investigate the effects of the signal-to-noise ratio and number of transits on the correct determination of the planetary parameters. Finally, we suggest choices of the parameters that govern the tool for the usage with real planets, according to the complexity of the problem and computational facilities.

Probing Transit Timing Variations of three hot Jupiters: HATP-36b, HATP-56b, and WASP-52b

ABSTRACT We report the results of new transit observations for the three hot Jupiter-like planets, HATP-36b, HATP-56b, and WASP-52b, respectively. Transit timing variations (TTVs) are presented for these systems based on observations that span the period 2016–2020. The data were collected with the 0.6-m telescope at Adiyaman University (ADYU60, Turkey) and the 1.0 m telescope at TÜBİTAK National Observatory (TUG, Turkey). Global fits were performed to the combined light curves for each system along with the corresponding radial velocity (RV) data taken from the literature. The extracted parameters (for all three systems) are found to be consistent with the values from previous studies. Through fits to the combined mid-transit times data from our observations and the data available in the literature, an updated linear ephemeris is obtained for each system. Although a number of potential outliers are noted in the respective O-C diagrams, the majority of the data are consistent within the 3σ confidence level implying a lack of convincing evidence for the existence of additional objects in the systems studied.