Issue 
A&A
Volume 668, December 2022



Article Number  A110  
Number of page(s)  14  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/202244592  
Published online  13 December 2022 
TESS search for substellar companions through pulsation timing of δ Scuti stars
I. Discovery of companions around Chang 134 and V393 Car
^{1}
Dipartimento di Fisica e Astronomia, Università degli Studi di Padova,
Vicolo dell’Osservatorio 3,
35122
Padova, Italy
email: valentina.vaulato@unige.ch
^{2}
INAF  Osservatorio Astronomico di Padova,
Vicolo dell’Osservatorio 5,
35122
Padova, Italy
email: valerio.nascimbeni@inaf.it
Received:
25
July
2022
Accepted:
2
September
2022
Earlytype mainsequence pulsating stars such as δ Scuti variables are one of the least explored class of targets in the search for exoplanets. Pulsation timing (PT) is an alternative technique to the most effective search methods. It exploits the lighttraveltime effect (LTTE) to infer the presence of additional massive bodies around a pulsating star by measuring a periodic phase modulation of its signal. PT has been extremely successful in discovering and characterizing stellar binaries when it was applied to highprecision light curves over long temporal baselines, such as those delivered by the Kepler mission. In favorable conditions, the sensitivity of PT can reach the planetarymass regime; one such candidate has already been claimed. The advent of TESS, with its nearly fullsky coverage and the availability of fullframe images, opens a great opportunity to expand this field of research. In this work, we present a pilot study aimed to understand the potential of PT as applied to TESS data, which are considerably different from Kepler data in terms of photometric noise, sampling cadence, and temporal baseline. We focused on the most favorable class of δ Scuti, that is, those showing strong pulsations and very simple frequency spectra. After the development of a customized pipeline, we were able to detect candidate companions for two targets in the (sub)stellar mass regime: Chang 134 (43 ± 5 M_{jup}, P ≃ 82 days) and V393 Car (≳100 M_{jup}, P ≳ 700 days). Our results also highlight the limiting factors of this technique and the importance of an accurate absolute time calibration for future missions such as PLATO.
Key words: stars: oscillations / stars: variables: delta Scuti / techniques: photometric / methods: data analysis / planets and satellites: detection / binaries: general
© V. Vaulato et al. 2022
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the SubscribetoOpen model. Subscribe to A&A to support open access publication.
1 Introduction
Since the pioneering detections of 51 Peg b, the first exoplanet orbiting a solartype star (Mayor & Queloz 1995), and HD209458b, the first transiting hot Jupiter (Henry et al. 2000; Charbonneau et al. 2000), contemporary exoplanetary research has mostly focused on latetype (FGKM) mainsequence hosts. This is not just because of the interesting analogies with our own Solar System, but also because these stars are particularly suited for the most fruitful detection techniques in terms of discoveries: radial velocities (RV) and transits. These two techniques are hampered by the presence of stellar activity and pulsations, making the detection and characterization of planets around very young, evolved, or earlytype stars extremely challenging. Along these lines, the somewhat standard road to discovery has become the photometric detection of transitlike signals, followed by an RV monitoring to confirm the planetary nature of the candidate and to measure its mass.
The vast majority of transiting planets known today have been discovered by dedicated spacebased telescopes: starting with the pioneering work of the CoRoT satellite (2006–2013; built by an international consortium led by CNES; Auvergne et al. 2009), which was followed by the Kepler mission by NASA (2009–2014; Borucki et al. 2010). Kepler was stopped because of a technical failure and was then restored for a secondphase observing program called K2 (2014–2018; Howell et al. 2014). Finally, the advent of the NASA Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015), launched in 2018 and currently operating, doubled the number of known candidate exoplanets by starting the first allsky spacebased transit survey. TESS is optimized to search for Earth to Neptunesized planets transiting bright and nearby mainsequence stars with a particular focus on cool, red dwarfs (K and M) because their smaller radii and lower masses make the transit and RV signals much stronger. TESS objects of interest are three magnitudes brighter than those of Kepler on average (Barclay et al. 2018), enabling an easier and more effective spectroscopic followup analysis. Most importantly in this context, and unlike its predecessors, TESS has the capability to download photometric data not only for a limited sample of preselected targets, but also for virtually any other bright star in the sky, although at a longer cadence with respect to the core sample. This is made possible by the availability of fullframe images (FFIs), originally downloaded every 30 min, and since sector 27, at a 10min cadence. FFI photometry opened a vast array of opportunities in many astrophysical fields, including extending the search for planets to a wider range in the stellar parameter space (Montalto et al. 2020; Nardiello et al. 2019).
Several alternative detection techniques have been proposed and implemented to search for exoplanets around stars other than solartype or M dwarfs, or in different Galactic environments. Atype stars, in particular, are one of the least explored class of planethosting stars. Of more than 5000 confirmed planets with known mass and/or radius listed in the NASA Exoplanet Archive, only 22 are hosted by 21 mainsequence stars earlier than the F0V spectral type (T_{eff} ≥ 7220 K; Pecaut & Mamajek 2013). With only four exceptions coming from direct imaging surveys, they were all detected through the transit technique (3 by HATNet, 3 by WASP, 5 by KELT, 2 by MASCARA, 3 by Kepler, 2 by TESS) and were confirmed by means of spectral tomography (Collier Cameron et al. 2010) and/or statistical validation techniques because classical RV measurements are usually not effective on such stars, whose spectral lines are few and Doppler broadened by their very high rotational velocities (v sin i ~ 100 km s^{−1}; Rainer et al. 2021).
Only three of the hosts described above are confirmed or suspected δ Scuti variables, that is, population I pulsators located in the classical Cepheid instability strip (Bowman & Kurtz 2018; Guzik 2021 and references therein): WASP33b (Herrero et al. 2011), HATP57b (Hartman et al. 2015), and WASP167b = KELT13b (Temple et al. 2017). Because ~50% of the stars in this region of the parameter space are expected to show δ Scuti pulsations to some extent (Bowman 2017), some astrophysical reason must explain the lack of planets around δ Scuti stars, and/or a strong observational selection effect is at play. The latter explanation is likely true, as the stellar pulsation pattern is known to limit the efficiency of transitsearch algorithms even after they are filtered in the frequency domain by sophisticated techniques (Ahlers et al. 2019; Hey et al. 2021).
In this context, the pulsation timing technique (PT; Hermes 2018) is a completely independent method for detecting massive companions around stars showing coherent and welldetectable pulsations, such as δ Scuti stars, pulsating white dwarfs, or hot subdwarfs. The underlying physical principle is the socalled lighttraveltime effect (LTTE) that was first discussed in detail by Irwin (1952): Because the speed of light c is finite, the motion of the target along the line of sight will translate into a phase shift of any periodic astrophysical signal originating from the star, including pulsations. The arrival times of the signals can then be compared with those predicted by a linear ephemeris (implying a constant periodicity dP/dt = 0) by computing the socalled O–C diagram (observedcalculated; Sterken 2005). Any deviation from a strict periodicity dP/dt ≠ 0 (assuming that the pulsation themselves are intrinsically stable) would then reveal itself as a nonzero second derivative in the O–C diagram (see Sect. 3 for details). The LTTE signal from a single, massive perturber on a Keplerian orbit is formally analogous to an RV curve having half of the eccentricity e and the argument of the periapsis ω decreased by π/2 (Irwin 1952); if the full shape of the LTTE signal is retrieved, then the orbital period P, the minimum mass M sin i, and other orbital elements of the perturber can be reliably measured. The PT approach shares some similarities with the astrometric method (Black & Scargle 1982) in that it measures the displacement of the target star with respect to an inertial rest frame, but projected along the line of sight rather than on the sky plane. Notably, also the underlying selection effects are similar because both methods are more sensitive to massive bodies at large star–planet separations. PT therefore probes a region of the parameter space that is complementary to the space that is investigated by transits and RVs, which are more sensitive to closein perturbers (Sozzetti 2005).
The PT technique requires photometric data spanning a temporal baseline that is long enough to sample the orbital period of the perturber and a signaltonoise ratio (S/N) high enough to constrain the phase shift of the oscillations to a level comparable to the expected LTTE amplitude^{1}. The main pulsation modes of the star have to be individually identified and measured on the Fourier spectrum of the light curve, implying that continuous time series spanning at least ~10 days are mandatory to achieve the needed resolution in the frequency space (~0.1 cycles per day) when dealing with the harmonic content of typical δ Scuti stars (Murphy et al. 2014). This is the main reason why sparse sampling imposed by the daynight cycle is the main limiting factor of groundbased observations. Targets with an extremely simple harmonic content, such as pseudosinusoidal pulsators, represent an exception, and a few pioneering results were published by analyzing groundbased data (Paparo et al. 1988; Barlow et al. 2011b). Interestingly, the latter claim was subsequently confirmed by an independent RV followup, demonstrating the reliability of this technique (Barlow et al. 2011a).
The full power of the PT technique, however, revealed itself when the nearly uninterrupted 4 yr baseline of Kepler photometry was exploited. Some works (Shibahashi & Kurtz 2012; Murphy et al. 2014, 2018, 2020; Shibahashi et al. 2015; Murphy & Shibahashi 2015), although they focused on the detection and characterization of stellarmass companions, demonstrated that under favorable assumptions, even LTTE signals from planetarymass companions (<13 M_{jup}) are in principle detectable through spacebased photometry (Murphy et al. 2014). Murphy et al. (2016) presented the PT detection of a massive planetary candidate for the first time (11.8 ± 0.7 M_{jup}, P = 840 ± 20 days). This candidate orbits a metalpoor δ Scuti star (KIC 7917485). Together with the discovery by Silvotti et al. (2007) of a giant planet around the hot subdwarf B star V391 Peg, these are the only PT discoveries of planetarymass objects published so far. It should be noted, however, that a followup paper on V391 Peg with new data (Silvotti et al. 2018) failed to fully reproduce the previous claim, probably due to nonlinear interactions between pulsation modes (see also Bowman et al. 2021). This acts as a caution.
While several studies of δ Scuti stars have been published that exploited TESS shortcadence data (Antoci et al. 2019; Hasanzadeh et al. 2021; Southworth 2021 for a review), none of them applied the PT technique. The goal of this paper is to focus on the PT analysis of a small sample of particularly suited targets observed by TESS at short cadence as a pilot study for the exploitation of TESS (and, later on, PLATO) data on a larger scale. In Sect. 2, we describe how our targets were identified and how their relevant stellar parameters were collected. In Sect. 3, we summarize all the steps of our data analysis, starting from the filtering of raw light curves through the harmonic analysis to measure the phase shifts as a function of time. We determine the bestfit orbital solutions through an LTTE model in Sect. 4 and discuss the results in Sect. 5.
2 Target selection and characterization
The vast majority of targets observed by TESS at short and ultrashort photometric cadence (120 and 20 s, respectively) are FGKM dwarfs selected for the core transitsearch survey, drawn from the socalled candidate target list (CTL; Fausnaugh et al. 2021). However, stars of different spectral types, including pulsating variables, are ingested mostly from GI/DDT proposals or from the TESS Asteroseismic Science Operation Center (TASOC^{2}), which provides the TESS Asteroseismic Science Consortium (TASC) with an asteroseismological database of the mission (Schofield et al. 2019).
In order to select our targets, we chose the catalog compiled by Chang et al. (2013) as a starting point. This catalog lists a sample of 1578 δ Scuti stars that were identified at high confidence level by several previous surveys. We crossmatched this catalog with the TESS CTL tables for sectors from 1 to 42 included to obtain only δ Scuti stars that were observed by TESS at 2min cadence. As an additional constraint, we further restricted our sample by (1) requiring a minimum data coverage of at least seven TESS sectors, not necessarily contiguous, to build O–C diagrams with a number of points and a temporal baseline long enough to detect LTTE signals at an approximately 1yr timescale (see Sect. 4); (2) selecting targets brighter than T = 12 in the TESS photometric system, to avoid being limited by photon and background noise; and (3) excluding all the stars that have been identified as binaries in the Chang et al. (2013) catalog. These additional constraints are justified by the nature of this study, which is not focused on a complete sample, but is rather intended as a pilot study to investigate the limiting factors of TESS (especially due the systematic errors) on a very small sample of the most favorable targets.
The final output of this selection process is a short list of 12 targets, which were then individually examined both through a literature search and by inspecting their TESS light curves. In particular, we confirmed that all of them are actually δ Scuti pulsators and carried out a preliminary harmonic analysis on them by computing the generalized Lomb–Scargle periodogram (GLS; Zechmeister & Kürster 2018) of a single sector with the same algorithms as described in Sect. 3. While most targets are classical δ Scuti stars showing smallamplitude (~0.01–0.03 mag) pulsations and a complex pattern of different modes in the frequency domain, two of them stand out as highamplitude pulsators (≥0.1 mag) with a particularly clean periodogram, in which most of the signal is due to a single, welldefined pulsation mode and its harmonics (Figs. 1 and 2): hereafter we refer to them as Chang 134 (from its ID number in the Chang et al. 2013 catalog) and V393 Carinae. This combination of a large pulsation amplitude and a very simple and coherent frequency pattern makes these targets well suited for a precise timing analysis, especially for our pilot study. For this reason, we focused this analysis on them and set aside the other targets (five of which are multimode pulsators, but with modes that are easy to identify) for the next paper of the series. A short review of the available scientific literature is given in the following subsections, and the basic astrophysical parameters of the two targets are reported in Table 1. It is worth mentioning that both our targets could fit the highamplitude δ Scuti subclass (HADS; Breger 2000; Antoci et al. 2019), which is characterized by largeamplitude pulsations, low rotational velocities, and timedomain spectra that are dominated by radial modes. We did not attempt to perform a detailed identification of their pulsation modes, however, because our work is focused on a dynamical technique (where the pulsation pattern is exploited just as a coherent astrophysical clock) rather than on stellar physics.
Fig. 1 TESS photometry of Chang 134 = TYC 91589191. Upper panel: oneday section of the light curve from TESS sector 1. Lower panel: GLS periodogram of the whole light curve including nine TESS sectors, stitched together and filtered as explained in Sect. 3. 
2.1 Chang 134 = TYC 91589191
Chang 134, also known as TYC 91589191 or TIC 431589510, was first identified as a variable star by ASAS (Pojmanski 2002) and classified as a generic δ Scuti star (DSCT variability class). No targeted followup study on it has been published since 2002. An asteroseismological measurement of log(g) and T_{eff} was included in the catalog by Barceló Forteza et al. (2020) through the v_{max} seismic index extracted from the TESS light curves in an automated fashion. All the astrophysical parameters we collected are reported in the second and third columns of Table 1.
The PT technique, just like astrometry or RVs, cannot directly infer the mass of perturbing body, but rather its ratio m_{p}/M_{*} with the stellar mass. The latter is needed at the analysis stage (Sect. 4) to properly interpret our results for the two targets. For Chang 134, an estimate of stellar mass (M_{*} = 1.40 ± 0.16 M_{⊙}) and age was first published in the catalog by Mints & Hekker (2017) using UniDAM models, but without taking advantage of spectroscopic data or Gaia parallaxes, which were not available at that time. Later, Chang 134 was included in the largescale analysis by Queiroz et al. (2020), combining highresolution spectra from the APOGEE2 survey DR16 with broadband photometric data taken from different sources and Gaia DR2 distances. Stellar parameters are derived as the posterior distribution returned by the Bayesian isochronefitting code StarHorse (Queiroz et al. 2018). The stellar mass reported for Chang 134 is M_{*} = 1.38 ± 0.03 M_{⊙}, where the central value is the 50th percentile of the StarHorse posterior distribution, while the symmetrized uncertainty is calculated as the halfdifference between the 84th and the 16th percentiles of the same distribution.
Fig. 2 TESS photometry of V393 Carinae = HD 66260. Upper panel: oneday section of the light curve from TESS sector 1. Lower panel: GLS periodogram of the whole light curve including 14 TESS sectors, stitched together and filtered as explained in Sect. 3. 
Adopted stellar parameters and other basic information for Chang 134 (second and third columns) and V393 Car (fourth and fifth columns).
2.2 V393 Carinae = HD 66260
Unlike Chang 134, the much brighter V393 Car (also known as HD 66260 or TIC 364399376) has been the subject of several targeted studies. It was first discovered as a variable star and classified as a δ Scuti radial pulsator by Helt (1984). Balona & Evers (1999) found in a more detailed analysis that the dominant mode is rather nonradial; the debate of this point is still ongoing. A first claim of change in the period of the main pulsation mode appeared in Garcia et al. (2001), who measured a 7.18min discrepancy with respect to the ephemeris by Helt (1984). Interestingly, further groundbased photometric followup by Axelsen (2014) retrieved the very same period of Helt (1984) within 0.33 s, and found no conclusive evidence of additional overtones. The present work therefore also represents an opportunity to confirm or disprove the claimed period change, and to determine whether it is caused by a secular or an oscillating term.
A measurement of v_{max}, log(g), T_{eff} was listed by for Chang 134 by Barceló Forteza et al. (2020). Unfortunately, V393 Car is missing from the Queiroz et al. (2020) catalog, but is listed in the Anders et al. (2022) catalog, which applied StarHorse on the Gaia EDR3 photometric and astrometric measurements, yielding M_{*} = 2.02 ± 0.10 M_{⊙} for this target. As accurate stellar parameters of V393 Car derived from spectroscopy are lacking in the literature, we attempted to obtain an independent estimate of the stellar mass as a crosscheck by applying the empirical relations derived by Moya et al. (2018) as a function of other stellar parameters. In particular, we exploited the relation (1)
where a = −0.119 ± 0.003, b = 2.14 × 10^{−5} ± 5 × 10^{−7} and c = 0.1837 ± 0.0011 are the coefficients calculated by Moya et al. (2018). This relation is valid within the range of temperatures 4780 ≤ T_{eff} ≤ 10990 K. Adopting T_{eff} from the seismic value given by Barceló Forteza et al. (2020) and L_{*} from the Gaia DR2 catalog^{3}, and propagating the errors, we computed a stellar mass M_{*} = 1.98 ± 0.17 M_{⊙} for V393 Car, which is perfectly consistent with the previous value. We therefore adopt the M_{*} = 2.02 ± 0.10 M_{⊙} estimate from Anders et al. (2022) in the subsequent analysis.
We are aware that a detailed asteroseismic analysis performed on all the available TESS light curves would likely yield much more accurate fundamental parameters for both our targets (Hasanzadeh et al. 2021). However, since an error of about 5% on stellar mass will not be the dominant source of uncertainty in our final parameters (as we discuss in Sect. 4), an analysis like this is beyond the scope of our work.
2.3 TESS light curves
Chang 134 and V393 Car were both observed by TESS in 2min cadence mode. At the time when our analysis started, Chang 134 was observed in 9 nonconsecutive sectors (12361327282936) between July 2018 and March 2021. V393 Car was observed in 14 nonconsecutive sectors (147810112728313435363738) between July 2018 and April 2021. All the photometric data analyzed in this work were processed and extracted from the raw data by the Science Processing Operations Center (SPOC) pipeline (Jenkins et al. 2016) and are publicly available on the Mikulski Archive for Space Telescopes (MAST^{4}). Specifically, we built the light curves for our analysis by extracting the stellar flux from the presearch data conditioning simple aperture photometry column (PDCSAP; Smith et al. 2012; Stumpe et al. 2012) because it results in cleaner time series since systematic longterm trends are removed.
3 Harmonic analysis
Our approach to data analysis is based on the fit of sums of harmonic functions to preselected segments of our light curves in the form (2)
where A_{i} is the amplitude, P_{i} the period, and ϕ_{i} the phase of each component normalized between 0 and 1. The phase of the signal at the dominant pulsation mode ϕ_{0} is then retrieved for each segment, giving the socalled phase shift ϕ_{0}(t) as a function of time. From this quantity, the usual O–C diagram in time units can be computed in a straightforward way by just multiplying the phase shift by the pulsation period P_{0}.
We developed an independent pipeline to perform all the needed tasks from the raw light curves down to the final O–C diagram. Its flow chart (described in the following sections) is conceptually similar to the phase modulation (PM) method described by Murphy et al. (2014), with some substantial differences, in particular, in how the actual fit is performed, that is, through a Markov chain Monte Carlo (MCMC) analysis in our case rather than with a simple leastsquares algorithm. Most of our processing steps were implemented by scripting and modifying routines from the VARTOOLS code version 1.39^{5} released by Hartman & Bakos (2016), to which we refer for a more indepth description of the individual algorithms.
3.1 Light curve preconditioning
We started our process by downloading TESS SPOC light curves of our targets from the MAST archive for each observing sector (23 combined sectors). Each data point was sampled every 2 min as a default for shortcadence observations.
As a first step, we discarded all data points flagged as bad, that is, with a quality factor q ≠ 0 (QUALITY column). We then extracted the time column TIME from each FITS file and converted it into the BJD_{TDB} standard (Eastman et al. 2010). Finally, we extracted the PDCSAP_FLUX column and its associated error PDCSAP_FLUX_ERR and converted them into the magnitude system. The choice of PDCSAP over SAP is justified by our need that systematic errors, especially those manifesting themselves as longterm trends, are corrected for as much as possible to avoid unnecessary noise on our periodograms. As our analysis is focused on retrieving the phase of the pulsation signal, we are not concerned by any small perturbation in the amplitudes that might be introduced by the PDCSAP processing (Cui et al. 2019). In order to filter out the most obvious outliers that survived the q = 0 selection, we carried out an iterative clipping at 15 σ with respect to the mean. This specific threshold was empirically set by confirming that the shape of the main pulsation mode was left unchanged.
The final step of our preconditioning recipe was to split the light curves from each sector into smaller chunks that were individually analyzed later to become single points in the O–C diagram. Working on continuous segments, about 10 days has been shown to be a reasonable compromise between the need of getting 1) a frequency resolution that is high enough to reliably measure the phases of the individual pulsation modes and 2) a time resolution on the O–C diagram as low as possible (Murphy et al. 2014). Because for most TESS sectors the only significant interruption is the oneday central gap to allow the data downlink toward Earth at each perigee of the spacecraft, we chose that gap as a natural boundary and split each sector accordingly into two mostly continuous segments of about 14 days each, to which we refer hereafter as “orbits.” We are therefore left with 18 orbits for Chang 134 and 28 orbits for V393 Car. Each orbit was individually normalized to unit flux (and magnitude zero) before we continued to the next stages.
3.2 Frequency identification
In order to identify the most significant pulsation modes of our targets we made use of the Generalized LombScargle periodogram (GLS, Zechmeister & Kürster 2018), an improved version of the LombScargle periodogram (LS, VanderPlas 2018) which is less sensitive to aliasing and returns more accurate frequencies. Because we searched for LTTE effects, all the frequencies might in principle be slightly varying as a function of time due to the phase shift itself. At this stage, we therefore need to calculate the average value over the full span of our observations for each frequency. This can be accomplished by stitching all the normalized orbits into a single light curve and running GLS on it. We limited our search from 1 to 50 cycles per day, thus including the typical frequency range of δ Scuti modes (Balona & Dziembowski 2011) and being well within the Nyquist limit for our sampling rate (~360 cycles per day). The output are two periodograms for Chang 134 and V393 Car that are visually similar to those shown in Figs. 1 and 2, but with a much higher S/N. In order to improve the estimate of the pulsation frequencies, the light curve is usually fitted by a multisinusoidal function using the periodogram peaks as a starting point (Silvotti et al. 2018). In our case, however, this additional step is not essential given the extremely high S/N of our data.
By identifying the strongest frequencies v_{i} (and corresponding periods P_{i}) as the most prominent peaks in the periodogram, we found one dominant frequency for Chang 124 that was followed by its multiple integer harmonics, sorted hereafter by decreasing GLS power, (3)
and in the same way for V393 Car, (4)
A complete list of all the significant frequencies detected for Chang 134 and V393 Car can be found in Tables B.1 and B.2, respectively.
For both stars, we obtained v_{1} = 2v_{0} and v_{2} = 3v_{0}, that is, the three most prominent frequencies are the second and third harmonics of the main pulsation mode. Chang 134 does not show any significant peak outside the harmonic series of v_{0}. The GLS power of the strongest peak for V393 Car outside the harmonics of v_{0} (at ~ 12.58288 day^{−1}) is three orders of magnitude smaller than the main mode. While irrelevant for our analysis due to its weakness, we note that this minor peak, detected at very high confidence in the TESS data, confirms the presence of an additional pulsation mode at 12.58 or 13.58 cycles per day claimed by Helt (1984) that was never confirmed by subsequent works (Garcia et al. 2001; Axelsen 2014). The very clean window function of the TESS data enables us to state that 13.58 was just a oneday alias of the real peak at 12.58 cycles per day.
3.3 Harmonic fit
After the average frequencies v_{i} = 2π/P_{i} of the most powerful pulsation modes were identified, we proceeded to fit a harmonic series in the form of Eq. (2) to our light curves. As several TESS orbits are affected by longterm systematic errors even after the PDCSAP treatment (on a timescale of typically days to weeks), we added a polynomial term as a function of time to our model to mitigate their impact. After some trial and error, we concluded that a polynomial of order 2 (i.e., a quadratic baseline) is enough to provide us with an accurate fit. The full model M we fit is therefore (5)
where t_{m} is the median time of each orbit, set to minimize the correlations between the fit b_{i} parameters. The P_{i} periods must be fixed to their average values determined in the previous section in order to obtain meaningful phase shifts; for N harmonic components included, the number of free parameters is therefore 3 + 2N. While N might in principle be arbitrarily large, for instance, when all the components that are even barely detectable in the frequency spectrum are included, almost all the information on the phase shift is contained in the three most significant peaks, which sum up to ≫99% of the total power. We therefore adopt N = 3 for our model by including the Pi constants from Eqs. (3) and (4), implying nine free parameters for our fit. We verified this assumption by adding further harmonic terms as a test, which resulted in statistically indistinguishable O–C diagrams at the expense of a much more intensive computation.
We fit the model in Eq. (5) through a differential evolution Monte Carlo Markov chain algorithm (DEMC; Ter Braak 2006), in which multiple MC chains are run in parallel and learn from each other how to converge to a global minimum in the parameter space, rather than running independently as in the classical approach. With respect to the leastsquares techniques used by most previous works to fit harmonic series to light curves, DEMC allowed us to obtain much more reliable error estimates on the final bestfit parameters because the posterior distribution takes any correlation between them into account. In order to increase the efficiency, we first fit only the oscillating term ∑_{i}f_{i} of our model to each orbit with a Levenberg–Marquardt method (Moré 1978), to obtain reasonably good starting points for A_{i} and ϕ_{i}. Then we ran the DEMC code, again orbit by orbit, to fit the full model with uninformative priors centered on the A_{i}, ϕ_{i} values found above, while the b_{i} parameters were initially set to zero. After 1 000 000 steps, the first 10% of the chain was discarded as burnin phase and the posterior distributions of the nine fit parameters were examined to ensure that convergence was reached.
At this stage, when we plotted the residuals of the light curves after the bestfit model had been subtracted, we realized that a nonnegligible number of outliers had escaped the previous filtering steps. These data points are mostly located close to momentum dumps of the spacecraft or at the beginning of each orbit when TESS cameras are not yet thermally stabilized. To deal with them, we discarded all the outliers at more than 5σ and fed the clipped light curves for a second run of the DEMC fit, with the same configuration. The final bestfit parameters, especially the phases ϕ_{i} to which we are most interested in, are statistically consistent with those delivered from the first run, but usually with slightly smaller uncertainties.
Finally, the phase shifts ϕ_{i}, by our definition within the [0,1] interval, were converted into O–C delays through a P_{i} factor. In principle, if the phase shifts were associated with independent pulsation modes and due solely to an LTTE, all the ϕ_{i} computed for the same orbit should result in the same O–C value, within the measurement errors: ϕ_{0}P_{0} = ϕ_{1}P_{1} = ϕ_{2}P_{2}, and computing a weighted average 〈ϕ_{i}P_{i}〉 would be the rigorous way to obtain the overall phase shift. This would also be a great opportunity to independently verify whether the LTTE model is the best explanation because different modes should yield the same O–C (Hermes 2018). In our case, however, all the ϕ_{i} come from a single pulsation mode, and the harmonics do not carry additional information with respect to the fundamental frequency v_{0}, P_{0}, where most of the spectral power lies: σ(ϕ_{0}) ≪ σ(ϕ_{1}) ≪ σ(ϕ_{2}) and thus ϕ_{0}P_{0} ≃ 〈ϕ_{i}P_{i}〉. For this reason, after confirming that the assumption above is true, we defined ϕ_{0}P_{0} as our O–C. We associated the median epoch of the orbit t_{m} as a time stamp for each value as introduced above: (O–C)(t_{m}) = ϕ_{0} × P_{0}. The resulting O–C diagrams for Chang 134 and V393 Car are plotted in Figs. 3 and 4. The formal error bars are 1.5 s on average (range: 1.1–2.6 s) and 0.16 s (range: 0.14–0.19 s), respectively.
Fig. 3 Orbital solution for Chang 134. Upper panel: O–C diagram for the main pulsation mode of Chang 134 in seconds. The bestfit LTTE solution (including a quadratic trend) is overplotted as a continuous cyan line, and the quadratic baseline is plotted as a dashed orange line (see Sect. 3 for details). Lower panel: residuals from the bestfit model in seconds. 
Fig. 4 Quadratic orbital solution for V393 Car. Upper panel: O–C diagram for the main pulsation mode of V393 Car in seconds. The bestfit LTTE solution with a quadratic baseline (O–C)quad (Eq. (7)) is overplotted as a continuous cyan line (see Sect. 3 for details). The quadratic baseline a_{0} + a_{1}t + a_{2}t^{2} is plotted separately as a dashed orange line. The nominal error bars (~0.16 s on average) are smaller than the point size. Lower panel: residuals from the bestfit model in seconds. 
4 LTTE modeling
Even a quick look at both our O–C diagrams (Figs. 3 and 4) reveals that neither can be reasonably fit by a straight line, that is, the resulting reduced is ≫ 1. The presence of a nonzero second derivative implies that the phase of the signal is evolving with time, and we can interpret this modulation by assuming that it is due to an LTTE from unseen companions. Following Kepler et al. (1991), we parameterize our model as (6)
where t_{0} is an arbitrary reference time, ∆P is the difference between the actual average pulsation period P and that estimated in Eqs. (3) and (4), accounts for a linear variation of P as a function of time, and a is the amplitude of an LTTE signal from a perturber on a circular^{6} orbit having an orbital period P_{orb} and with a phase φ.
From a physical point of view, the first term of Eq. (6) does not carry information because it represents just a change in the slope of the O–C and does not imply any variation of the pulsation period. The second term, in the LTTE framework, represents a constant acceleration with respect to the barycenter of the system, possibly caused by a massive, perturbing body whose P_{orb} is much longer than our observing baseline (≈1000 days), such as stellar companions on very wide orbits. We should be aware here that several physical mechanisms other than LTTE can result in a P ≠ 0 term, including nonlinear interactions between different pulsation modes and stellar evolution effects. We return on this point in Sect. 4.2.
Last, the third, oscillating term of Eq. (6) is the LTTE modulation we are most interested in. As the first two terms play the role of a quadratic baseline in our fit, we make it more explicit by rewriting the model with a_{0} = t_{0}, a_{1} = ∆P and , (7)
where the median epoch of each orbit t_{m} has been subtracted from the time variable in order to minimize the correlations between a_{0}, a_{1}, and a_{2}.
In order to confirm the uniqueness of our bestfit solution and the significance of the quadratic coefficient a_{2}, we also performed a second fit to each data set by fixing a_{2} = 0 in our model, that is, by imposing a linear baseline, (8)
According to the Irwin (1952) model, for a circular orbit (eccentricity e = 0) and assuming that the perturbing body is much less massive than the star (m_{p} ≪ M_{*}), the amplitude of the LTTE signal is (9)
where c is the speed of light, a is the orbital semimajor axis of the perturber and i is the inclination of its orbital plane with respect to the sky plane. In other words, when a and P_{orb} are measured from the O–C data, the stellar mass is known (from Table 1) and a is derived from Kepler’s laws, we can constrain the minimum mass of the perturber as m_{p} sin i, with the same degeneracy on i as for the RV and astrometric techniques, (10)
We started the fitting process by computing a GLS periodogram on our O–C diagrams to search for periodic modulations. For both targets, a single, welldefined peak stands out, at about 82 days and 1000 days for Chang 134 and V393 Car, respectively. We then fit our full LTTE model with both a quadratic and linear baseline (parameterized as in Eqs. (7) and (8), respectively) on our O–C diagrams through DEMC, setting uninformative priors on all the five or six free parameters involved: a_{0}, a_{1}, a_{2}, α, P_{orb}, and φ in the (O–C)_{quad} model, and α_{0}, α_{1}, α, P_{orb}, and φ in the (O–C)_{lin} model, and centering the boundaries of P_{orb} on the values previously found by GLS to speed up the convergence. As done at the harmonic analysis stage, after 1000 000 steps, the first 10% of the chain was then discarded as burnin phase, and we extracted the bestfit values and the associated errors of our parameters from their posterior distributions. All these values are reported in Table 2 for the four independent fits. The mass of the perturbing body m_{p} sin i was calculated as well as a derived parameter by propagating the relevant errors in Eq. (10).
4.1 Orbital solution for Chang 134
The bestfit (O–C)_{quad} model for Chang 134 is overplotted on our O–C in the upper panel of Fig. 3 as a solid cyan line, with the quadratic baseline a_{0} + a_{1}t + a_{2}t^{2} plotted as a dashed orange line. The resulting χ^{2} of the residuals (lower panel of Fig. 3) is 56.7, with the reduced χ^{2} being , where DOF = 18 − 6 = 12 are the degrees of freedom of our fit. It is worth noting that the bestfit value for a_{2} is statistically consistent with zero (Table 2, second column), that is, we do not see any evidence of longterm LTTE effects because the derived is also consistent with zero. This is confirmed by the fact that the bestfit linear model (O – C)_{lin} (Table 2, first column) yielded the same results within the error bars for all the fit parameters, with a virtually indistinguishable χ^{2} (56.4 vs. 56.7). In other words, because the (O–C)_{quad} model is not favored by our data, we hereafter adopt the results from the (O–C)_{lin} fit.
The bestfit returned the amplitude of the LTTE term to be α = 6.11 ± 0.65 s (according to Eq. (6)), and after propagating the errors through Eq. (9), we obtain a minimum mass of m_{p} sin i = 43.1 ± 4.7 M_{jup} for the hypothetical companion. This is within the brown dwarf regime (13–80 M_{jup}; Grieves et al. 2021).
4.2 Orbital solution for V393 Car
Unlike the case of Chang 134, the results from V393 Car fits are more difficult to interpret because the two models (O–C)_{quad} and (O–C)_{lin}, while having a similar shape, led us to different orbital solutions. The two bestfit models are overplotted on our O–C in the upper panel of Figs. 4 and 5 as a solid cyan line, respectively, with the linear and quadratic baseline a_{0} + a_{1}t + a_{2}t^{2} plotted as a dashed orange line.
While it is clear that at least an oscillating term is needed to fit the O–C data, the inclusion of a quadratic term yields different parameters for the LTTE sinusoid: P ≃ 1071 days with amplitude α ≃ 107 s for the (O–C)_{lin} fit versus P ≃ 723 days with amplitude α ≃ 48 s for the (O–C)_{quad} fit. The discrepancy of a factor of ~2 in the α parameter translates into a similar ratio for the derived hypothetical companion masses: m_{p} sin i = 175.7 ± 9.2 M_{jup}, or 0.17 ± 0.01 M_{⊙} (linear model) and m_{p} sin i = 102.0 ± 5.3 M_{jup}, or 0.10 ± 0.01 M_{⊙} (quadratic model). We emphasize that for both scenarios, the minimum mass would be consistent with a cool dwarf of spectral type M5V and M6V, respectively (Pecaut & Mamajek 2013).
A technique that is commonly employed for model selection problems like this is the evaluation of metrics such as the Akaike information criterion (AIC; Akaike 1974) or the Bayesian information criterion (BIC; Schwarz 1978). In the case of Gaussian distributions, the latter can be written as BIC = χ^{2} + k ln n, where k is the number of free parameters and n is the number of data points. If applied literally to our problem, the BIC would overwhelmingly favor our (O–C)_{quad} scenario (∆BIC ≫ 10). Unfortunately, this result is distorted by the fact that the χ^{2} of the residuals for both fits (lower panel of Fig. 4 and 5) is exceptionally high due to the extremely low formal error bars with respect to the systematic errors at play (on which we comment in Sect. 5): for the (O–C)_{quad} fit, for instance, we obtained χ^{2} = 8627, with the reduced χ^{2} being , where DOF = 28 − 6 = 22 are the degrees of freedom of our fit. This fact, and the very sparse phase sampling of our signal, prevents us from reaching a firm conclusion on the properties of this candidate companion, which will require further data to be confirmed and better constrained. New TESS observations of V393 Car are planned for Cycle 5, during four consecutive sectors (61–64; from January to May 2023), and again in Sector 68 (August 2023). Since the predictions from the (O–C)_{quad} and (O–C)_{lin} models will diverge by ~12 h by mid2023, we expect that the extension of our analysis to the new TESS light curves will definitely be conclusive about the correct scenario.
As a side note, we could wonder whether the quadratic term we found when we assumed a quadratic baseline (a_{2} = 650 ± 4 ppm s^{−1} as fitted on our O–C plane, translating into yr^{−1}) might be consistent with other nonLTTE mechanisms, as previously mentioned. Significant nonlinear interactions between pulsating modes (Silvotti et al. 2018; Bowman et al. 2021) are not at play here because the spectral power outside the main pulsation mode and its harmonics is essentially negligible. On the other hand, the typical expected from evolutionary effects range from 10^{−9} to 10^{−7} yr^{−1} (Breger & Pamyatnykh 1998; Xue et al. 2022), that is, it is much lower than our fitted value.
The same reasoning can be applied to Chang 134, for which we measure a nonsignificant quadratic term a_{2} = −2.2 ± 7.5 ppm s^{−1} corresponding to yr^{−1}.
Unlike V393 Car, this upper limit is still within the predictions for evolutionary effects. For instance, if we adopt yr^{−1} (as measured on AE UMa by Xue et al. 2022), we obtain a total contribution on the O–C of ~0.35 s over the full baseline of our observations (~1000 days), which is well below our current measurement error.
Fig. 5 Linear orbital solution for V393 Car. Upper panel: same as Fig. 4, but with the bestfit LTTE solution with a linear baseline (O–C)_{lin} (Eq. (8)). The linear baseline a_{0} + a_{1}t is plotted separately as a dashed orange line. Lower panel: residuals from the bestfit model in seconds. 
5 Discussion and conclusions
We applied the PT technique to two δ Scuti stars observed by TESS in a large number of sectors, detecting a periodic modulation on their O–C that is consistent with the presence of companions in both cases: a BD (m_{p} sin i = 43.1 ± 4.7 M_{jup}) on a P ≃ 82 days orbit around Chang 134, and a more massive body (m_{p} sin i ≃ 0.10–0.17 M_{⊙}) around V393 Car, whose orbital parameters are still not firmly constrained by the limited phase sampling of the signal. This double detection on the very first two targets analyzed by our project may appear to be very lucky. Still, the binary fraction of δ Scuti primaries as estimated by a previous PT search on Kepler data is rather high: 15.4 ± 1.4% (Murphy et al. 2018). In addition, unlike the quoted work, we initially restricted our sample in order to target only the most favorable stars in terms of sensitivity to LTTE signals. The overall binary fraction of Atype stars is estimated to be significantly larger than 50% by most authors (Duchêne & Kraus 2013), making our results less surprising. Any further statistical implication is prevented by the small size of our sample and by the heavy selection effects at play. As a side note, however, we mention that Borgniet et al. (2019) constrained the BD frequency within 2–3 au around AF stars to be below 4% (1σ) based on 225 targets observed with SOPHIE and/or HARPS, and that other authors hypothesized the presence of a brown dwarf desert extending to earlytype stars (Murphy et al. 2018). This would make the companion of Chang 134 an uncommon object, worthwhile to be followed up with more data and/or different techniques. The future TESS observations (at least five sectors of Cycle 5 are already scheduled, starting in March 2023) will also possibly help in constraining the eccentricity of the LTTE orbit, forced to zero in our analysis to avoid overfitting given our limited phase coverage and small number of O–C points.
In addition to our results for these two specific stars, our analysis was based on an independent and improved implementation of the PT method (Murphy et al. 2014). It is also intended as a pilot study that is a prelude to a systematic search of LTTE companions around pulsating stars on a much larger scale, exploiting the huge sample of δ Scuti targets that are and will be monitored through the availability of TESS FFIs at 10min cadence. The results for Chang 134, in particular, demonstrates that the detection of bodies in the substellar regime is perfectly feasible from TESS light curves, and that in some particularly favorable configuration, candidates with planetary masses (m_{p} sin i ≲ 13 M_{jup}) might be already within reach. A perturber with the same LTTE amplitude as was found on Chang 134 (α ≃ 7 s), for instance, but on a larger orbit at P ≃ 800 days, would imply a companion with a mass of ≃2 M_{jup} detected at high significance.
On the other hand, our analysis highlighted several limiting factors that are worth discussing further. The first factor, largely selfevident in the case of V393 Car (Sect. 4.2), is related to the sparse sampling of the TESS light curves, an unavoidable consequence of its fixed scanning law (Ricker et al. 2015). Even for targets within or close to the CVZ at ecliptic latitudes β ≳ 78°, that is, targets that are observed in up to 13 contiguous sectors, there is always a large gap every other year, which make the detection of signals at longer periods (P ≳ 1 yr) difficult or ambiguous. Unfortunately, this is also the region in the parameter space where the PT technique is more sensitive to lowmass perturbers, owing to the Eq. (9) relation. As TESS is continuing its mission through extension phases, the accumulated data will mitigate this issue by filling out all or most of the orbital phases on the folded O–C diagram.
A second, subtler limiting factor is the excess scatter in our O–C data, as can be seen from the very high values reported in Table 2 and discussed Sect. 4.2. Our residuals from the bestfit LTTE model have an rms of ~3 s in both fits, regardless of the average nominal error bars, which are about ten times lower on V393 Car than on Chang 134. The most plausible cause for this are systematic errors in the absolute calibration of the time stamp. This issue has been anticipated by the TESS team, and in particular by the TESS Asteroseismic Science Consortium (TASC), who gave accurate timing requirements in the SAC/TESS/0002/6 document^{7}, including RSTASC05: “[…] the time given for each exposure should be accurate over a period of 10 days to better than 1 second.” No requirement is given at longer timescales, however, and the TESS data release notes (DRN^{8}) have reported systematic offsets in the time stamps for several sectors on the order of seconds, that is, consistent with the excess O–C scatter we measured. While some of the reported offsets have been already corrected in subsequent update data releases, others are still waiting to be implemented. A groundbased campaign by von Essen et al. (2020) tried to independently confirm the absolute time calibration of TESS by comparing observations of a sample of eclipsing binaries in common, measuring a global offset of 5.8 ± 2.5 s; the errors of the measurements at individual epochs, however, are too large (≳ 10 s) to confirm or disprove the systematic errors we see. We emphasize that offsets of a few seconds are usually completely negligible when investigating transiting exoplanets, but they are crucial in extending the sensitivity of the PT technique to planetmass companions. An important outcome of a largescale PT analysis of TESS data will allow us also to identify and correct these systematic errors by comparing the timing residuals of tens or hundreds of targets, following a selfcalibration approach.
The same technique as applied in this work, and in general, the expertise gained through this project, will also help the scientific exploitation of PLATO (Rauer et al. 2014), for which the PT analysis could be a compelling case of ancillary science. Unlike TESS, no FFIs will be downloaded from PLATO during the nominal observing phase, and no stars earlier than F5 are to be included in the main target samples, which are focused on transit search around solartype stars (Montalto et al. 2021). Nevertheless, about 8% of the PLATO science data rate will be allocated to the General Observer (GO) program through ESA calls open to the whole astrophysical community. Any δ Scuti star (and, more in general, pulsating star) allocated as GO within a longpointing field (Nascimbeni et al. 2022) would result in a mostly uninterrupted light curve on a baseline of 2–4 yr, with a cadence^{9} and photometric precision far better than what is achieved by TESS, overcoming the sparse phase sampling affecting the latter. Even more important, the time stamps of each PLATO data point will be accurate within 1 s in an absolute sense by formal scientific requirement because of the specifically designed calibration processes that also include the monitoring of a preselected set of detached eclipsing binaries (Nascimbeni et al. 2022).
Acknowledgements
The authors wish to thank the referee, Dr. Roberto Silvotti, for the valuable and fruitful comments and suggestions which significantly improved our manuscript. This research has made use of the SIMBAD database (operated at CDS, Strasbourg, France;Wenger et al. 2000), the VARTOOLS Light Curve Analysis Program (version 1.39 released October 30, 2020, Hartman & Bakos 2016), TOPCAT and STILTS (Taylor 2005, 2006). Valerio Nascimbeni and Giampaolo Piotto recognize support by ASI under program PLATOASI/INAF agreements 2015019R.12018.
Note added in proof
After the acceptance of this paper, we were contacted by Dr. D. Hey, who let us know that he found a similar orbital solution for V393 Car using the maelstrom code (Hey et al. 2020). In particular, his results suggest that the linearbaseline model yields a better fit the timing data.
Appendix A O–C data tables
O–C data points, bestfit models, and residuals for Chang 134.
O–C data points, best−fit models, and residuals for V393 Car.
Appendix B Detected frequencies
Significant GLS frequencies detected for Chang 134.
Significant GLS frequencies detected for V393 Car.
References
 Ahlers, J. P., Barnes, J. W., & Myers, S. A. 2019, AJ, 158, 88 [CrossRef] [Google Scholar]
 Akaike, H. 1974, IEEE Trans. Automatic Control, 19, 716 [Google Scholar]
 Anders, F., Khalatyan, A., Queiroz, A. B. A., et al. 2022, A&A, 658, A91 [Google Scholar]
 Antoci, V., Cunha, M. S., Bowman, D. M., et al. 2019, MNRAS, 490, 4040 [Google Scholar]
 Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 [Google Scholar]
 Axelsen, R. A. 2014, J. Am. Ass. Var. Star Obs., 42, 292 [NASA ADS] [Google Scholar]
 BailerJones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147 [Google Scholar]
 Balona, L. A., & Dziembowski, W. A. 2011, MNRAS, 417, 591 [Google Scholar]
 Balona, L. A., & Evers, E. A. 1999, MNRAS, 302, 349 [Google Scholar]
 Barceló Forteza, S., Moya, A., Barrado, D., et al. 2020, A&A, 638, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barclay, T., Pepper, J., & Quintana, E. V. 2018, ApJS, 239, 2 [Google Scholar]
 Barlow, B. N., Dunlap, B. H., & Clemens, J. C. 2011a, ApJ, 737, L2 [Google Scholar]
 Barlow, B. N., Dunlap, B. H., Clemens, J. C., et al. 2011b, MNRAS, 414, 3434 [Google Scholar]
 Black, D. C., & Scargle, J. D. 1982, ApJ, 263, 854 [NASA ADS] [CrossRef] [Google Scholar]
 Borgniet, S., Lagrange, A. M., Meunier, N., et al. 2019, A&A, 621, A87 [Google Scholar]
 Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [Google Scholar]
 Bowman, D. M. 2017, Amplitude Modulation of Pulsation Modes in Delta Scuti Stars [Google Scholar]
 Bowman, D. M., & Kurtz, D. W. 2018, MNRAS, 476, 3169 [Google Scholar]
 Bowman, D. M., Hermans, J., DaszyńskaDaszkiewicz, J., et al. 2021, MNRAS, 504, 4039 [Google Scholar]
 Breger, M. 2000, in Astronomical Society of the Pacific Conference Series, Delta Scuti and Related Stars, eds. M. Breger, & M. Montgomery, 210, 3 [NASA ADS] [Google Scholar]
 Breger, M., & Pamyatnykh, A. A. 1998, A&A, 332, 958 [NASA ADS] [Google Scholar]
 Chang, S.W., Protopapas, P., Kim, D.W., & Byun, Y.I. 2013, AJ, 145, 132 [NASA ADS] [CrossRef] [Google Scholar]
 Charbonneau, D., Brown, T. M., Latham, D. W., & Mayor, M. 2000, ApJ, 529, L45 [Google Scholar]
 Collier Cameron, A., Guenther, E., Smalley, B., et al. 2010, MNRAS, 407, 507 [Google Scholar]
 Cui, K., Liu, J., Yang, S., et al. 2019, MNRAS, 489, 5513 [Google Scholar]
 Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
 Duchêne, G., & Kraus, A. 2013, ARA&A, 51, 269 [Google Scholar]
 Eastman, J., Siverd, R., & Gaudi, B. S. 2010, PASP, 122, 935 [Google Scholar]
 Fausnaugh, M., Morgan, E., Vanderspek, R., et al. 2021, PASP, 133, 095002 [Google Scholar]
 Gaia Collaboration (Brown, A.G.A., et al.) 2018, A&A, 616, A1 [Google Scholar]
 Gaia Collaboration (Brown, A.G.A., et al.) 2021, A&A, 649, A1 [Google Scholar]
 Gaia Collaboration (Vallenari, A., et al.) 2022, A&A, in press https://doi.org/10.1051/00046361/202243940 [Google Scholar]
 Garcia, J. R., Busso, V. A., & Ferrero Sosa, G. 2001, J. Am. Ass. Var. Star Obs., 30, 44 [NASA ADS] [Google Scholar]
 Grieves, N., Bouchy, F., Lendl, M., et al. 2021, A&A, 652, A127 [Google Scholar]
 Guzik, J. A. 2021, Front. Astron. Space Sci., 8, 55 [NASA ADS] [Google Scholar]
 Hartman, J. D., Bakos, G. Á., Buchhave, L. A., et al. 2015, AJ, 150, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Hartman, J. D., & Bakos, G. Á. 2016, Astron. Comput., 17, 1 [Google Scholar]
 Hasanzadeh, A., Safari, H., & Ghasemi, H. 2021, MNRAS, 505, 1476 [Google Scholar]
 Helt, B. E. 1984, A&As, 56, 457 [NASA ADS] [Google Scholar]
 Henry, G. W., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2000, ApJ, 529, L41 [Google Scholar]
 Hermes, J. J. 2018, Handbook of Exoplanets, 787 [CrossRef] [Google Scholar]
 Herrero, E., Morales, J. C., Ribas, I., & Naves, R. 2011, A&A, 526, A10 [Google Scholar]
 Hey, D. R., Murphy, S. J., ForemanMackey, D., et al. 2020, AJ, 159, 202 [NASA ADS] [CrossRef] [Google Scholar]
 Hey, D. R., Montet, B. T., Pope, B. J. S., Murphy, S. J., & Bedding, T. R. 2021, AJ, 162, 204 [NASA ADS] [CrossRef] [Google Scholar]
 Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [Google Scholar]
 Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398 [Google Scholar]
 Irwin, J. B. 1952, ApJ, 116, 211 [Google Scholar]
 Jenkins, J. M., Twicken, J. D., McCauliff, S., et al. 2016, SPIE Conf. Ser., 9913, 99133E [Google Scholar]
 Kepler, S. O., Winget, D. E., Nather, R. E., et al. 1991, ApJ, 378, L45 [NASA ADS] [CrossRef] [Google Scholar]
 Kim, C. H., Kreiner, J. M., Zakrzewski, B., et al. 2018, ApJS, 235, 41 [Google Scholar]
 Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [Google Scholar]
 Mints, A., & Hekker, S. 2017, A&A, 604, A108 [Google Scholar]
 Montalto, M., Borsato, L., Granata, V., et al. 2020, MNRAS, 498, 1726 [Google Scholar]
 Montalto, M., Piotto, G., Marrese, P. M., et al. 2021, A&A, 653, A98 [Google Scholar]
 Moré, J. J. 1978, in Lecture Notes in Mathematics (Berlin: Springer Verlag), 630, 105 [Google Scholar]
 Moya, A., Zuccarino, F., Chaplin, W. J., & Davies, G. R. 2018, ApJS, 237, 21 [Google Scholar]
 Murphy, S. J., & Shibahashi, H. 2015, MNRAS, 450, 4475 [Google Scholar]
 Murphy, S. J., Bedding, T. R., Shibahashi, H., Kurtz, D. W., & Kjeldsen, H. 2014, MNRAS, 441, 2515 [Google Scholar]
 Murphy, S. J., Bedding, T. R., & Shibahashi, H. 2016, ApJ, 827, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Murphy, S. J., Moe, M., Kurtz, D. W., et al. 2018, MNRAS, 474, 4322 [Google Scholar]
 Murphy, S. J., Barbara, N. H., Hey, D., Bedding, T. R., & Fulcher, B. D. 2020, MNRAS, 493, 5382 [Google Scholar]
 Nardiello, D., Borsato, L., Piotto, G., et al. 2019, MNRAS, 490, 3806 [Google Scholar]
 Nascimbeni, V., Piotto, G., Börner, A., et al. 2022, A&A, 658, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Paparo, M., Szeidl, B., & Mahdy, H. A. 1988, Ap&SS, 149, 73 [NASA ADS] [CrossRef] [Google Scholar]
 Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9 [Google Scholar]
 Pojmanski, G. 2002, Acta Astron., 52, 397 [NASA ADS] [Google Scholar]
 Queiroz, A. B. A., Anders, F., Santiago, B. X., et al. 2018, MNRAS, 476, 2556 [Google Scholar]
 Queiroz, A. B. A., Anders, F., Chiappini, C., et al. 2020, A&A, 638, A76 [Google Scholar]
 Rainer, M., Borsa, F., Pino, L., et al. 2021, A&A, 649, A29 [Google Scholar]
 Rauer, H., Catala, C., Aerts, C., et al. 2014, Exp. Astron., 38, 249 [Google Scholar]
 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telescopes Instrum. Syst., 1, 014003 [Google Scholar]
 Schneider, J. 2005, in Astronomical Society of the Pacific Conference Series, The LightTime Effect in Astrophysics: Causes and cures of the OC diagram, ed. C. Sterken, 335, 191 [NASA ADS] [Google Scholar]
 Schofield, M., Chaplin, W. J., Huber, D., et al. 2019, ApJS, 241, 12 [Google Scholar]
 Schwarz, G. 1978, Ann. Stat., 6, 461 [Google Scholar]
 Shibahashi, H., & Kurtz, D. W. 2012, MNRAS, 422, 738 [Google Scholar]
 Shibahashi, H., Kurtz, D. W., & Murphy, S. J. 2015, MNRAS, 450, 3999 [Google Scholar]
 Silvotti, R., Schuh, S., Janulis, R., et al. 2007, Nature, 449, 189 [Google Scholar]
 Silvotti, R., Schuh, S., Kim, S. L., et al. 2018, A&A, 611, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000 [Google Scholar]
 Southworth, J. 2021, Universe, 7, 369 [Google Scholar]
 Sozzetti, A. 2005, PASP, 117, 1021 [Google Scholar]
 Stassun, K. G., Oelkers, R. J., Paegert, M., et al. 2019, AJ, 158, 138 [Google Scholar]
 Sterken, C. 2005, in Astronomical Society of the Pacific Conference Series, The LightTime Effect in Astrophysics: Causes and cures of the OC diagram, ed. C. Sterken, 335, 3 [NASA ADS] [Google Scholar]
 Stumpe, M. C., Smith, J. C., Van Cleve, J. E., et al. 2012, PASP, 124, 985 [Google Scholar]
 Taylor, M. B. 2005, in Astronomical Society of the Pacific Conference Series, Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert, 347, 29 [NASA ADS] [Google Scholar]
 Taylor, M. B. 2006, in Astronomical Society of the Pacific Conference Series, Astronomical Data Analysis Software and Systems XV, eds. C. Gabriel, C. Arviset, D. Ponz, & S. Enrique, 351, 666 [NASA ADS] [Google Scholar]
 Temple, L. Y., Hellier, C., Albrow, M. D., et al. 2017, MNRAS, 471, 2743 [Google Scholar]
 Ter Braak, C. J. F. 2006, Stat. Comput., 16, 239 [Google Scholar]
 VanderPlas, J. T. 2018, ApJS, 236, 16 [Google Scholar]
 von Essen, C., Lund, M. N., Handberg, R., et al. 2020, AJ, 160, 34 [NASA ADS] [CrossRef] [Google Scholar]
 Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9 [Google Scholar]
 Xue, H.F., Niu, J.S., & Fu, J.N. 2022, Res. Astron. Astrophys., 22, 10 [Google Scholar]
 Zechmeister, M., & Kürster, M. 2018, Astrophysics Source Code Library [record ascl:1807.019] [Google Scholar]
As a comparison, the LTTE effect induced by the orbital motion of Jupiter around the Sun, as seen edgeon, is ~10 s (Schneider 2005), with aperiod of ~12 yr. See also Sect. 5.
Unfortunately, the stellar parameters reported for V393 Car by Gaia DR3 (Gaia Collaboration 2022) are highly discrepant with each other and associated with unreasonably small errors, possibly because the pipeline is unable to deal with a pulsating star. We therefore adopted the DR2 values for consistency.
Throughout this paper we assume a circular orbit as a simplifying hypothesis for our LTTE model, although a significant fraction of exoplanetary systems and binary stars are known to have nonzero eccentricity (Kim et al. 2018). This is justified by our need to avoid two additional free parameters (e and ω) and is further discussed in Sect. 5.
The nominal cadence of the “normal” PLATO cameras will be 25 s, i.e., a factor of about 5 shorter than the TESS standard cadence (Rauer et al. 2014).
All Tables
Adopted stellar parameters and other basic information for Chang 134 (second and third columns) and V393 Car (fourth and fifth columns).
All Figures
Fig. 1 TESS photometry of Chang 134 = TYC 91589191. Upper panel: oneday section of the light curve from TESS sector 1. Lower panel: GLS periodogram of the whole light curve including nine TESS sectors, stitched together and filtered as explained in Sect. 3. 

In the text 
Fig. 2 TESS photometry of V393 Carinae = HD 66260. Upper panel: oneday section of the light curve from TESS sector 1. Lower panel: GLS periodogram of the whole light curve including 14 TESS sectors, stitched together and filtered as explained in Sect. 3. 

In the text 
Fig. 3 Orbital solution for Chang 134. Upper panel: O–C diagram for the main pulsation mode of Chang 134 in seconds. The bestfit LTTE solution (including a quadratic trend) is overplotted as a continuous cyan line, and the quadratic baseline is plotted as a dashed orange line (see Sect. 3 for details). Lower panel: residuals from the bestfit model in seconds. 

In the text 
Fig. 4 Quadratic orbital solution for V393 Car. Upper panel: O–C diagram for the main pulsation mode of V393 Car in seconds. The bestfit LTTE solution with a quadratic baseline (O–C)quad (Eq. (7)) is overplotted as a continuous cyan line (see Sect. 3 for details). The quadratic baseline a_{0} + a_{1}t + a_{2}t^{2} is plotted separately as a dashed orange line. The nominal error bars (~0.16 s on average) are smaller than the point size. Lower panel: residuals from the bestfit model in seconds. 

In the text 
Fig. 5 Linear orbital solution for V393 Car. Upper panel: same as Fig. 4, but with the bestfit LTTE solution with a linear baseline (O–C)_{lin} (Eq. (8)). The linear baseline a_{0} + a_{1}t is plotted separately as a dashed orange line. Lower panel: residuals from the bestfit model in seconds. 

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