# TOI-1431b/MASCARA-5b: A Highly Irradiated Ultra-Hot Jupiter Orbiting One of the Hottest & Brightest Known Exoplanet Host Stars

Brett Addison, Emil Knudstrup, Ian Wong, Guillaume Hébrard, Patrick

Dorval, Ignas Snellen, Simon Albrecht, Aaron Bello-Arufe, Jose-Manuel

Almenara, Isabelle Boisse, et al.

► To cite this version:

Brett Addison, Emil Knudstrup, Ian Wong, Guillaume Hébrard, Patrick Dorval, et al. TOI-1431b/MASCARA-5b: A Highly Irradiated Ultra-Hot Jupiter Orbiting One of the Hottest & Brightest Known Exoplanet Host Stars. The Astronomical Journal, 2021, 162 (6), pp.292. 10.3847/1538-3881/ac224e. insu-03458424

HAL Id: insu-03458424

https://insu.hal.science/insu-03458424v1

Submitted on 2 Aug 2025

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

CrossMark

TOI-1431b/MASCARA-5b: A Highly Irradiated Ultrahot Jupiter Orbiting One of the Hottest and Brightest Known Exoplanet Host Stars

Brett C. Addison¹ ID, Emil Knudstrup² ID, Ian Wong³,⁴⁸ ID, Guillaume Hébrard⁴, Patrick Dorval⁵, Ignas Snellen⁵ ID,
Simon Albrecht² ID, Aaron Bello-Arufe⁶ ID, Jose-Manuel Almenara⁷ ID, Isabelle Boisse⁸, Xavier Bonfils⁷, Shweta Dalal⁴,
Olivier D. S. Demangeon⁹,¹⁰,¹¹, Sergio Hoyer⁸ ID, Flavien Kiefer¹², N. C. Santos⁹,¹¹ ID, Grzegorz Nowak¹³,¹⁴ ID,
Rafael Luque¹³,¹⁴ ID, Monika Stangret¹³,¹⁴ ID, Enric Palle¹³,¹⁴ ID, René Tronsgaard⁶ ID, Victoria Antoci²,⁶ ID, Lars A. Buchhave⁶ ID,
Maximilian N. Günther¹⁵,⁴⁹ ID, Tansu Daylan¹⁵,⁵⁰ ID, Felipe Murgas¹³,¹⁴ ID, Hannu Parviainen¹³,¹⁴ ID,
Emma Esparza-Borges¹³,¹⁴ ID, Nicolas Crouzet¹⁶ ID, Norio Narita¹³,¹⁷,¹⁸,¹⁹ ID, Akihiko Fukui¹³,¹⁷ ID, Kiyoe Kawauchi²⁰ ID,
Noriharu Watanabe¹⁹,²¹,²² ID, Markus Rabus²³,²⁴,²⁵ ID, Marshall C. Johnson²⁴ ID, Gilles P. P. L. Otten²,⁸,²⁶ ID, Geert Jan Talens²⁷ ID,
Samuel H. C. Cabot²⁸ ID, Debra A. Fischer²⁸ ID, Frank Grundahl² ID, Mads Fredslund Andersen² ID, Jens Jessen-Hansen²,
Pere Palle¹³,¹⁴ ID, Avi Shporer¹⁵ ID, David R. Ciardi²⁹ ID, Jake T. Clark¹ ID, Robert A. Wittenmyer¹ ID, Duncan J. Wright¹,
Jonathan Horner¹ ID, Karen A. Collins³⁰ ID, Eric L. N. Jensen³¹ ID, John F. Kielkopf³² ID, Richard P. Schwarz³³ ID, Gregor Srdoc³⁴,
Mesut Yilmaz³⁵ ID, Hakan Volkan Senavci³⁵ ID, Brendan Diamond³⁶, Daniel Harbeck²⁴ ID, Thaddeus D. Komacek³⁷,³⁸ ID,
Jeffrey C. Smith³⁹,⁴⁰, Songhu Wang⁴¹ ID, Jason D. Eastman³⁰ ID, Keivan G. Stassun⁴²,⁴³ ID, David W. Latham³⁰ ID,
Roland Vanderspek¹⁵ ID, Sara Seager³,¹⁵,⁴⁴ ID, Joshua N. Winn⁴⁵ ID, Jon M. Jenkins³⁹ ID, Dana R. Louie³⁷,⁴⁶ ID,
Luke G. Bouma⁴⁵ ID, Joseph D. Twicken³⁹,⁴⁰ ID, Alan M. Levine¹⁵ ID, and Brian McLean⁴⁷ ID

¹ University of Southern Queensland, Centre for Astrophysics, USQ Toowoomba, West Street, QLD 4350, Australia; Brett.Addison@usq.edu.au

² Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark

³ Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA

⁴ Institut d'Astrophysique de Paris, UMR7095 CNRS, Université Pierre & Marie Curie, 98bis boulevard Arago, F-75014 Paris, France

⁵ Leiden Observatory, Leiden University, Postbus 9513, 2300 RA Leiden, The Netherlands

⁶ DTU Space, National Space Institute, Technical University of Denmark, Elektrovej 328, DK-2800 Kgs. Lyngby, Denmark

⁷ Université Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France

⁸ Aix Marseille Univ., CNRS, CNES, LAM, Marseille, France

⁹ Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal

¹⁰ Centro de Astrofísica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal

¹¹ Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto, Rua Campo Alegre, 4169-007, Porto, Portugal

¹² LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris, 5 place Jules Janssen, F-92195 Meudon, France

¹³ Instituto de Astrofísica de Canarias (IAC), E-38200 La Laguna, Tenerife, Spain

¹⁴ Dept. Astrofísica, Universidad de La Laguna (ULL), E-38206 La Laguna, Tenerife, Spain

¹⁵ MIT Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA

¹⁶ European Space Agency (ESA), European Space Research and Technology Centre (ESTEC), Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands

¹⁷ Komaba Institute for Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan

¹⁸ JST, PRESTO, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan

¹⁹ Astrobiology Center, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

²⁰ Department of Earth and Planetary Science, The University of Tokyo, Tokyo, Japan

²¹ Department of Astronomical Science, The Graduate University for Advanced Studies, SOKENDAI, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

²² National Astronomical Observatory of Japan, NINS, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

²³ Departamento de Matemática y Física Aplicadas, Universidad Católica de la Santísima Concepción, Alonso de Rivera 2850, Concepción, Chile

²⁴ Las Cumbres Observatory, 6740 Cortona Dr., Ste. 102, Goleta, CA 93117, USA

²⁵ Department of Physics, University of California, Santa Barbara, CA 93106-9530, USA

²⁶ Academia Sinica, Institute of Astronomy and Astrophysics, 11F Astronomy-Mathematics Building, NTU/AS campus, No. 1, Section 4, Roosevelt Rd., Taipei 10617, Taiwan

²⁷ Institut de Recherche sur les Exoplanètes, Département de Physique, Université de Montréal, Montréal, QC H3C 3J7, Canada

²⁸ Department of Astronomy, Yale University, 52 Hillhouse Ave, New Haven, CT 06511, USA

²⁹ NASA Exoplanet Science Institute, Caltech/IPAC, 1200 East California Avenue, Pasadena, CA 91125, USA

³⁰ Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA

³¹ Department of Physics & Astronomy, Swarthmore College, Swarthmore, PA 19081, USA

³² Department of Physics and Astronomy, University of Louisville, Louisville, KY 40292, USA

³³ Patashnick Voorheesville Observatory, Voorheesville, NY 12186, USA

³⁴ Kotizarovci Observatory, Sarsoni 90, 51216 Viskovo, Croatia

³⁵ Ankara University, Department of Astronomy and Space Sciences, TR-06100, Ankara, Turkey

³⁶ Howard Community College, 10901 Little Patuxent Pkwy, Columbia, MD 21044, USA

³⁷ Department of Astronomy, University of Maryland, College Park, MD 20742, USA

³⁸ Department of the Geophysical Sciences, The University of Chicago, Chicago, IL 60637, USA

³⁹ NASA Ames Research Center, Moffett Field, CA 94035, USA

⁴⁰ SETI Institute, 189 Bernardo Ave., Suite 200, Mountain View, CA 94043, USA

⁴¹ Astronomy Department, Indiana University Bloomington, 727 East 3rd Street, Swain West 318, IN 47405, USA

⁴² Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA

⁴³ Department of Physics, Fisk University, Nashville, TN 37208, USA

⁴⁴ Department of Aeronautics and Astronautics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
Received 2021 April 11; revised 2021 August 11; accepted 2021 August 26; published 2021 December 6

Abstract

We present the discovery of a highly irradiated and moderately inflated ultrahot Jupiter, TOI-1431b/MASCARA-5 b (HD 201033b), first detected by NASA's Transiting Exoplanet Survey Satellite mission (TESS) and the Multi-site All-Sky Camera (MASCARA). The signal was established to be of planetary origin through radial velocity measurements obtained using SONG, SOPHIE, FIES, NRES, and EXPRES, which show a reflex motion of $K = 294.1 \pm 1.1$ m s⁻¹. A joint analysis of the TESS and ground-based photometry and radial velocity measurements reveals that TOI-1431b has a mass of $M_p = 3.12 \pm 0.18 M_J$ (990 ± 60 $M_⊕$), an inflated radius of $R_p = 1.49 \pm 0.05 R_J$ (16.7 ± 0.6 $R_⊕$), and an orbital period of $P = 2.650237 \pm 0.000003$ days. Analysis of the spectral energy distribution of the host star reveals that the planet orbits a bright ($V = 8.049$ mag) and young (0.29<sup>+0.32</sup><sub>-0.19</sub> Gyr) Am type star with $T_{\text{eff}} = 7690_{-250}^{+400}$ K, resulting in a highly irradiated planet with an incident flux of $\langle F \rangle = 7.24_{-0.64}^{+0.68} \times 10^9$ erg s⁻¹ cm⁻² (5300<sup>+500</sup><sub>-470</sub> $S_⊕$) and an equilibrium temperature of $T_{\text{eq}} = 2370 \pm 70$ K. TESS photometry also reveals a secondary eclipse with a depth of 127<sup>+4</sup><sub>-5</sub> ppm as well as the full phase curve of the planet's thermal emission in the red-optical. This has allowed us to measure the dayside and nightside temperature of its atmosphere as $T_{\text{day}} = 3004 \pm 64$ K and $T_{\text{night}} = 2583 \pm 63$ K, the second hottest measured nightside temperature. The planet's low day/night temperature contrast (~420 K) suggests very efficient heat transport between the dayside and nightside hemispheres. Given the host star brightness and estimated secondary eclipse depth of ~1000 ppm in the K band, the secondary eclipse is potentially detectable at near-IR wavelengths with ground-based facilities, and the planet is ideal for intensive atmospheric characterization through transmission and emission spectroscopy from space missions such as the James Webb Space Telescope and the Atmospheric Remote-sensing Infrared Exoplanet Large-survey.

*Unified Astronomy Thesaurus concepts: Exoplanets (498); Hot Jupiters (753); Transit photometry (1709); Radial velocity (1332)*

*Supporting material:* data behind figures, machine-readable table

1. Introduction

Exoplanet discoveries over the past 25 years have revealed that the diversity of alien worlds to be far greater than we had ever imagined. Planets more massive than Jupiter were found on orbits measured in mere days (and became known as “hot Jupiters,” e.g., Butler et al. 1998; Henry et al. 2000; Tinney et al. 2001), while others were found moving on highly eccentric orbits that had more in common with the solar system’s cometary bodies than its planets (e.g., Naef et al. 2001; Jones et al. 2006; Wittenmyer et al. 2017, 2019; Bergmann et al. 2021). And more recently, a new class of close-in gas giants has emerged that are hotter than even the coolest main-sequence stars. These superhot gas giants are known as ultrahot Jupiters and have dayside temperatures ≳2200 K (e.g., Parmentier et al. 2018).

In the early years of the exoplanet era, the majority of new exoplanet discoveries were made by radial velocity surveys (e.g., Mayor & Queloz 1995; Cochran et al. 1997; Butler et al. 1999; Pepe et al. 2004; Endl et al. 2008; Vogt et al. 2010) that mainly targeted stars considered to be “Sun-like”: of spectral classes late-F, G, and K. The focus on such stars was, in part, driven by technical necessity—such stars possess ample spectral lines for radial velocity analysis and typically rotate slowly enough that such lines are sufficiently narrow to yield precise radial velocity measurements. As a result, during the first two decades of the exoplanet era, our knowledge of the

diversity of planets was heavily biased toward systems with host stars similar to the Sun.

The situation changed with the launch of the Kepler space observatory, designed to search for exoplanets using the transit technique (Borucki et al. 2010). By surveying a vast number of stars at once, Kepler was able to get a true handle on the frequency of short-period planets around a wide variety of stars. It showed that planets are ubiquitous—a natural byproduct of the star formation process—and that planetary systems have a wide variety of architectures and compositions (e.g., Rogers 2015; Millholland et al. 2016, 2017; Guenther et al. 2017; Cañas et al. 2019; Dorn et al. 2019; Hsu et al. 2019; Bryson et al. 2021).

Although the transit technique is better suited to finding planets around hot, massive stars than the radial velocity method, such discoveries remain challenging due to the inherent difficulty in confirming planet candidates with radial velocity measurements. This is because more massive stars tend to have rapid rotation, resulting in a significant decrease in the obtainable radial velocity precision (e.g., WASP-167b/KELT-13b, KELT-17b, KELT-20b/MASCARA-2b, KELT-19Ab, HAT-P-69b, WASP-189b, HAT-P-70b, KELT-25b, KELT-26b, Zhou et al. 2016; Lund et al. 2017; Temple et al. 2017; Anderson et al. 2018; Talens et al. 2018; Siverd et al. 2018a; Zhou et al. 2019a; Martínez et al. 2020, transiting planets orbiting hot stars exclusively or primarily confirmed through Rossiter–McLaughlin or Doppler tomography observations).

NASA’s Transiting Exoplanet Survey Satellite (TESS: Ricker et al. 2015) was launched in 2018 and is the successor to the hugely successful Kepler and K2 missions. While the Kepler spacecraft could only examine a ~115° patch of sky (~0.3% of the total area of the full sky) at any given time (performing a narrow but deep survey of the sky), TESS is

<footnote n="45"><sup>45</sup> Department of Astrophysical Sciences, Princeton University, Peyton Hall, 4 Ivy Lane, Princeton, NJ 08544, USA</footnote>

<footnote n="46"><sup>46</sup> NASA Goddard Space Flight Center, 8800 Greenbelt Rd, Greenbelt, MD 20771, USA</footnote>

<footnote n="47"><sup>47</sup> Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA</footnote>

<footnote n="48"><sup>48</sup> Pegasi b Fellow.</footnote>

<footnote n="49"><sup>49</sup> Juan Carlos Torres Fellow.</footnote>

<footnote n="50"><sup>50</sup> Kavli Fellow.</footnote>
designed to survey large bands of the sky at a time (96° by 24°). As a result, over the two years of TESS' primary mission, it surveyed the majority of the sky, obtaining high-cadence observations (one target pixel image recorded every two minutes) of more than 200,000 pre-selected main-sequence dwarf stars, and low-cadence observations (one full frame image every 30 minutes) of 20 million bright stars (Ricker et al. 2015). In doing so, it has enabled the search for planets around a wide variety of bright stars (including hot and more massive stars historically avoided in radial velocity surveys)—perfect targets for ground-based follow-up observations. While we are still in the early stages of reaping the full harvest from TESS' remarkable wealth of observations, the spacecraft has already discovered a variety of new exoplanet systems (e.g., Huang et al. 2018; Cañas et al. 2019; Günther et al. 2019; Jones et al. 2019; Newton et al. 2019; Wang et al. 2019; Addison et al. 2021; Brahm et al. 2020; Davis et al. 2020; Gilbert et al. 2020; Jordán et al. 2020; Plavchan et al. 2020; Rodriguez et al. 2020; Teske et al. 2020) with many planet candidates, or TESS Objects of Interests (TOIs), awaiting confirmation (Guerrero et al. 2021).

A great benefit of TESS' discoveries orbiting bright stars is that those stars make ideal targets for follow-up observations, allowing us to better characterize the planets that orbit them. The most obvious route by which such characterization is achieved is through radial velocity observations, which provide measurements of the masses of the planets. Additionally, radial velocity observations also provide vital information on the planet's orbit, including its eccentricity and orbital period (critical if only one or two transits are observed). These observations can also be used to reveal the presence of otherwise undetectable, nontransiting planetary companions (e.g., Mills et al. 2019; Weiss et al. 2020; Sozzetti et al. 2021). When radial velocity observations are made during a transit of a planet across the disk of its host star, they can be used to constrain the inclination of the planet's orbit relative to the plane of the host star's equator through the use of the "Rossiter-McLaughlin effect" (e.g., McLaughlin 1924; Rossiter 1924).

Such observations have yielded startling insights in recent years. Where predictions made based on our knowledge of the solar system would suggest that exoplanets should typically move in roughly the same plane as their host star's equator, Rossiter-McLaughlin observations have revealed that a significant fraction of hot Jupiters move on misaligned orbits. Surprisingly, such observations have even revealed a number of planets moving on polar or retrograde orbits (e.g., Narita et al. 2009; Winn et al. 2009; Hirano et al. 2011; Albrecht et al. 2012; Addison et al. 2013)—a scenario markedly different from the near coplanarity of the solar system (as discussed in Horner et al. 2020 and references therein).

The highly excited orbits of such inclined hot Jupiters are likely linked to their origins. A number of different mechanisms have been proposed to produce hot Jupiters, ranging from the smooth and relatively sedate inward migration of the giant planets through interactions with their host star's protoplanetary disk (e.g., Lin et al. 1996; Ward 1997; Tanaka et al. 2002; Alibert et al. 2005), to planet-planet scattering (where close encounters between two planets throw them both onto highly eccentric orbits, dropping one inwards, toward its host, and flinging the other outwards; see, e.g., Chatterjee et al. 2008; Beaugé & Nesvorný 2012; Li et al. 2019a), and even distant perturbations from binary companions (through the Kozai-Lidov

mechanism; see, e.g., Kozai 1962; Lidov 1962; Nagasawa et al. 2008). Those mechanisms that invoke a period of high orbital eccentricity then require the tidal circularization of the planet's orbit, which, for small pericenter distances, can happen on timescales far shorter than the lifetimes of the stars involved (Nagasawa et al. 2008).

These three methods would naturally produce very different populations of hot Jupiters. Those formed through disk migration would be expected to orbit in essentially the same plane as their host star's equator. Those formed through planet-planet scattering would likely also congregate around such low inclinations, but would be more dispersed/excited, with scattering events able to produce at least moderate orbital inclinations for the planets involved. Planets that migrate as a result of the Kozai-Lidov mechanism, by contrast, would be expected to be found on highly misaligned orbits, because the Kozai-Lidov mechanism causes significant inclination excitation as a byproduct of the eccentricity excitation driving the planetary migration.

Given that more massive stars are statistically more likely to host binary companions (e.g., Perryman et al. 1998; Böhm-Vitense 2007; Zinnecker & Yorke 2007), especially Am stars such as TOI-1431 (see Böhm-Vitense 2006), and also tend to form more massive planets (e.g., Ida & Lin 2005; Bowler et al. 2010; Reffert et al. 2015), it would seem reasonable to expect that such stars might be more likely to host both moderately and highly misaligned planets (since the two mechanisms by which excited hot Jupiters can be produced would be more likely to occur). Indeed, as the catalog of planets for which spin-orbit alignment measurements have been made has grown, this prediction appears to have been borne out. The fraction of highly misaligned planets appears to be a strong function of stellar effective temperature (and to a lesser degree stellar mass), with hotter and more massive stars far more likely to host such planets than those that are cooler and less massive (e.g., Triaud et al. 2010; Winn et al. 2010; Louden et al. 2021).

In this work, we expand the catalog of planets orbiting hot main-sequence stars with the discovery and characterization of TOI-1431 b/MASCARA-5b—an inflated, ultrahot Jupiter orbiting a massive, hot, and metal-peculiar Am star.

In Section 2, we describe the observations and data reduction. The host star properties derived from the analysis of the broadband spectral energy distribution (SED) and spectroscopy are outlined in Section 3. The data analysis, modeling, and results are presented in Section 4. We then discuss the phase curve and secondary eclipse analysis and measurements of the planet's dayside and nightside temperatures in Section 4.1. Lastly, the discussion and conclusions are provided in Sections 5 and 6, respectively.

## 2. Observations and Data Reduction

TOI-1431/MASCARA-5 (HD 201033) was observed with space-based photometry by TESS (Section 2.1) and ground-based photometry from several facilities (Section 2.2). We also obtained follow-up spectroscopic observations from several ground-based facilities (Section 2.3) as well as high-contrast imaging on the Keck II telescope (Section 2.4) to establish the planetary nature of the transit signals detected by TESS and the Multi-site All-Sky Camera (MASCARA). Here we describe the observations collected and the data reduction process.
## 2.1. TESS Photometry

The star TOI-1431 (HD 201033; TIC 375506058, Stassun et al. 2019) was observed in Sectors 15 (on Camera 2 and CCD chip number 4) and 16 (on Camera 2 and CCD chip number 3) by TESS in 2 minute cadence mode nearly continuously between 2019 August 15 and 2019 October 7. The photometric data were processed by the TESS Science Processing Operations Center (SPOC) pipeline (see Jenkins et al. 2016 for a description of the SPOC pipeline), resulting in two versions of the light curves: simple aperture photometry (SAP, see Twicken et al. 2010; Morris et al. 2020) and multiscale maximum a posteriori (MAP) presearch data conditioning (PDC, see Smith et al. 2012; Stumpe et al. 2012, 2014). Multiscale MAP PDC divides each light curve into three wavelet bandpasses and then performs a MAP-like correction in each using separate cotrending basis vectors derived within each bandpass. A fit to neighboring targets is used as the Bayesian prior during the a posteriori fit. We downloaded both versions of the light curves from NASA's Mikulski Archive for Space Telescopes (MAST). Additionally, we produced a third version of the TESS light curve, single-scale MAP PDC (where no band-splitting is performed, Smith et al. 2012), for the analysis described in Section 4.

The transiting planet candidate was promoted to TOI status (see TOI catalog, Guerrero et al. 2021) and designated TOI-1431.01 by the TESS Science Office based on model fit results and a passing grade on all diagnostic tests in the SPOC Data Validation (DV) report (Twicken et al. 2018; Li et al. 2019b) for Sector 15. Twenty transits (10 per sector) were observed by TESS. Each transit had a duration of ~2.5 hr, a depth of ~6000 parts per million (ppm), and recurred with a period of 2.65 days. The first detected transit occurred on BJD$_{TDB}$ 2,458,712 in Sector 15. The raw SAP and multiscale MAP PDC light curves are shown in Figure 1.

Before detrending and fitting the SAP and the two versions of the PDC light curves, we first removed all quality-flagged data. The light curves were then split into a total of 10 segments, with each segment split at the spacecraft momentum dumps<fnref n="51" /> at intervals of 4.25 and 5.83 days in Sectors 15 and 16, respectively, to minimize the offset effects on the light curves during the detrending and fitting analysis. The transits were then masked and a 5σ median filter was applied to the out-of-transit data to remove remaining outliers in each light-curve segment. Lastly, the light-curve segments were normalized with the mean of the out-of-transit flux before performing the detrending and fitting analysis as described in Section 4.

## 2.2. Ground-based Transit Photometry

Transit signals of TOI-1431b were first detected by the ground-based MASCARA telescope (Talens et al. 2017) at the Roque de los Muchachos Observatory, La Palma, starting 2015 February. Transits were observed by MASCARA between 2015 February and 2018 March and were discovered before TESS observed them. Figure 2 shows the phase-folded light curve of TOI-1431b from MASCARA with the transit signal clearly detected at the same period of ~2.65 days and approximate transit depth as measured by TESS. The MASCARA observations were taken with exposure times of

6.4 s without using any filters and achieved a photometric precision of ~5000 ppm per 5 m binned observation.

Both MASCARA and TESS have large on-sky pixel sizes of $1'$ pixel$^{-1}$ (Talens et al. 2017) and $0.35$ pixel$^{-1}$ (Ricker et al. 2015), respectively. As such, further photometric follow-up is required to rule out potential false positive scenarios such as nearby eclipsing binaries, where an eclipsing binary pair falls on or very close to the same pixel as the target star in the MASCARA and TESS images.

To confirm the transit signals detected by TESS and MASCARA were coming from TOI-1431 and to check for any changes in transit depth between multiple photometric bands (chromaticity), we acquired additional photometric follow-up observations from several ground-based facilities through the TESS Follow-up Observing Program Working Group as summarized in Table 3 in the Appendix. We used the TESS Transit Finder, which is a customized version of the Tapir software package (Jensen 2013), to schedule the observations. These facilities include the 0.36 m CDK14 telescope at Howard Community College, the MuSCAT2 imager on the 1.52 m Carlos Sanchez Telescope (TCS) at the Teide Observatory (OT), the 0.8 m telescope at the Ankara University Kreiken Observatory (AUKR), the 0.3 m telescope at the Kotizarovci Observatory (SCT) in Rijeka, Croatia, the 0.6 m University of Louisville Manner Telescope (ULMT) at Mt Lemmon, and the Las Cumbres Observatory Global Telescope (LCOGT; Brown et al. 2013) 1 m node at the McDonald Observatory. For the final global analysis with Allesfitter, we included only the high-quality and high-precision (rms ≲3 ppt) ground-based observations with full transit coverage that are free of significant systematics to provide precise transit depth measurements and accurate transit ephemeris. As such, we have not included the light curve taken on 2020 May 24 from MuSCAT2 (but we do include the 2020 May 16 observation) due to significant systematics (see Section 2.2.2), ULMT taken on 2020 September 28 due to systematics and moderately low precision photometry obtained (~3.4 ppt), nor the photometry from MASCARA due to its relatively low photometric precision of ~5 ppt.

### 2.2.1. CDK14

One full transit of TOI-1431b was observed on 2019 December 24–25 with the CDK14 telescope at Howard Community College. The observations were taken using alternating 10 s exposures in the Sloan $g'$ and $z_s$ filters, starting at 23:02 UT on December 24 at an airmass of 1.2 and finishing at 02:56 UT on December 25 at an airmass of 2.4. We reduced the data and extracted the light curves for each filter using the AstroImageJ (AIJ) software package (Collins et al. 2017) with 7 pixel apertures. Figure 3 shows the transit signatures in the two bands with fitted transit depths of $5.1_{-0.7}^{+0.6}$ ppt and $5.1_{-0.6}^{+0.5}$ ppt for the Sloan $g'$ and $z_s$ filters, respectively. These transit depths are consistent with the SPOC pipeline transit depth of $5.68 \pm 0.03$ ppt to within ~1σ.

### 2.2.2. MuSCAT2

Two full transits of TOI-1431b were observed with MuSCAT2, one on 2020 May 16 and a second one on 2020 May 24, using simultaneous multicolor photometry in $g'$, $r'$, $i'$, and $z_s$ bands (Narita et al. 2019). For the joint fitting, we only used the transit observed on 2020 May 16 since the second

<footnote n="51">⁵¹ As provided in the data release notes found at https://archive.stsci.edu/tess/tess_drn.html.</footnote>
<figure><img src="image_3.png" /></figure>

**Figure 1.** The TESS light curves of TOI-1431 from Sectors 15 (left) and 16 (right). The top panels show the simple aperture photometry and the bottom panels show the preresearch data conditioning versions of the light curves. The circles plotted in blue represent data that have a quality flag greater than 0 (but not equal to 128) and the red squares are data with quality flag equal to 128 (flagged as manual exclude, i.e., the cadence was excluded because of an anomaly). The TESS quality flags are described in the TESS Data Products Description Document (see https://archive.stsci.edu/missions/tess/doc/EXP-TESS-ARC-ICD-TM-0014-Rev-F.pdf.

<figure><img src="image_4.png" /></figure>

transit observation was affected by systematics due to weather, resulting in different transit depths between the filter bands and residuals larger than 1.0 ppt. On the night of 2020 May 16, defocused observations were taken from 01:50 UT to 05:22 UT, with exposure times of 8 s for all channels, except for the *r*′, where an exposure time of 5 s was used for autoguiding. We reduced the data and extracted the light curves using the standard MuSCAT2 pipeline (Parviainen et al. 2019), recovering the expected transit signal as shown in Figure 3. Figure 14 in the Appendix shows the transit light curve in all four bands with fits to the transits and systematics, fits to the transits with systematics removed, and the resulting residuals.

### 2.2.3. AUKR

A full transit of TOI-1431b was observed with the 0.8 m Prof. Dr. Berahitdin Albayrak Telescope (T80) at the AUKR on 2020 June 16 in the Sloan *z*<sub>s</sub> band. Observations started approximately 30 minutes before transit ingress at 21:15 UT and lasted until 01:13 UT the following day, with exposure times of 40 s, and the target was observed at an airmass between 1.0 and 1.3. We extracted the photometry using the AIJ software package with 12 pixel (8′′16) apertures, and the resulting fit to the light curve is shown in Figure 4. The measured transit depth of 4.5<sup>+0.5</sup><sub>-0.6</sub> ppt is moderately shallower than the transit depths found from the TESS and other ground-based light curves, likely due to unaccounted for systematics.

### 2.2.4. SCT

We used the 0.3 m telescope at the Kotizarovci Observatory near Viskovo, Croatia, to observe a full transit of TOI-1431b on 2020 August 8. The transit was observed in the TESS-like band with 20 s exposures continuously for nearly 6.5 hr starting at 20:23 UT and ending at 02:45 UT in low airmass that ranged from 1′′1 at the start to 1′′3 at the end of the observations. The photometry was processed using the AIJ software package with an aperture of 11 pixels (13′′2) and Figure 4 shows the resulting fit to the light curve. We measure a transit depth of 5.6<sup>+0.2</sup><sub>-0.3</sub> ppt, consistent with the depth measured from the TESS light curves.

### 2.2.5. ULMT

On 2020 September 20, we observed a full transit of TOI-1431b from the 0.6 m ULMT at Mt Lemmon using the Sloan *z*′ filter. We set the exposure length to 32 s and started the observations at 04:29 UT. The observing was interrupted between 06:57 UT and 07:31 UT, around the start of transit ingress, due to a guiding error and then resumed afterwards until 10:26 UT. The target was observed at an airmass of 1.1 at the start and airmass of 1.8 at the end of the observations. The images were bias- and dark-corrected and the light curve was extracted using an aperture of 13 pixels (5′′1) with the AIJ software package. We measure a transit depth of 5.8 ± 0.3 ppt
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Photometry of TOI-1431 taken from the Multi-site All-Sky Camera (MASCARA) telescope at the Roque de los Muchachos Observatory, La Palma. The top panels are the box least-squares periodograms from the MASCARA light curves showing the strongest peak period at 2.65023 days. The middle panels show the phase-folded light curve from MASCARA, full phase on the left and zoomed in on the transit on the right. The blue points represent binned data in 0.003 phase steps, the black points are the original observations in bins of 50 images (320 s), and the red line is the best-fit box least-squares model (see Kovács et al. 2002) used in detecting the transit. The bottom panels show the photometry phase-folded at half and double the period (bottom left and right, respectively) found at the strongest peak in the periodogram. No obvious dip features are seen in the half-period phase-folded plot but two dip features are apparent in the double-period phase-folded plot at phase 0 and 0.5, as expected for the transiting planet with a period of 2.65023 days. The MASCARA light curves were not used in the global modeling in Allesfitter.</figcaption></figure>

(see Figure 4), in excellent agreement with the depth measured from the TESS photometry.

### 2.2.6. LCOGT

A full transit of TOI-1431b was observed on 2020 October 14 from the LCOGT 1.0 m network node at McDonald Observatory near Fort Davis, Texas (LCO-McD), in PAN-STARRS Y band. The SINISTRO camera images were calibrated using the standard LCOGT BANZAI pipeline (McCully et al. 2018), and the photometry was extracted using the AIJ software package with 20 pixel (7/8) apertures. Figure 4 shows a clear transit detection with an individually fitted transit depth of 5.7 ± 0.3 ppt, consistent with the SPOC pipeline transit depth of 5.68 ± 0.03 ppt to within $1\sigma$.

## 2.3. Spectroscopic Observations

We obtained high-resolution (HR) spectroscopic observations of TOI-1431 with SONG, SOPHIE, FIES, NRES (TLV and ELP), and EXPRESS to establish the planetary nature of the TESS transiting candidate and measure its mass. Here we describe the observations from each spectrograph and list the derived radial velocities (RVs) in Table 4 in the Appendix.

### 2.3.1. SONG

High-resolution spectroscopic observations of TOI-1431 were obtained using the robotic Stellar Observations Network Group (SONG) 1 m Hertzsprung telescope at the Teide Observatory in Tenerife (Andersen et al. 2014; Fredslund Andersen et al. 2019). In total, 19 spectra were taken with SONG between 2020 March 6 and 2020 June 29. SONG is equipped with an HR coudé echelle spectrograph with wavelength coverage between 4400 and 6900 Å across 51 spectral orders. The observations were taken using the iodine cell to allow for precise wavelength calibration. For each observation, we took a 2400 s exposure in slit mode 6 (slit width 1/22), providing a spectral resolution of $R=90,000$, a median count per pixel at 5560 Å of 1795.2 ADU (signal-to-noise ratio (S/N) of at least 50 per resolution element for all observations), and median radial velocity precision of 29 m s⁻¹.

Before each target observation, calibration frames (including bias frames, flat fields, and ThAr spectra) were taken on the same night and applied to target images. The target observations were then reduced and the spectra extracted with a pipeline written in Python. This pipeline uses a C++ implementation of the optimal extraction method from Ritter et al. (2014). We also acquired observations of the bright fast-rotating star HR 5191 that was used as an intrinsic stellar template for determining the spectral-line-spread function of
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

**Figure 3.** The ground-based transit light curves of TOI-1431b/MASCARA-5 b as a function of time from the observations taken with the two CDK14 photometric bands and the four MuSCAT2 photometric bands. The normalized photometric measurements are represented by the gray points while the blue points are the binned photometry. The 15 red lines are the randomly drawn transit light-curve models from the nested sampling posterior.

(The data used to create this figure are available.)

the spectrograph and to deconvolve a high S/N spectrum of
TOI-1431 taken without the iodine cell. We then analyzed the
extracted spectra using the iSONG pipeline (Antoci et al. 2013)
to produce radial velocities (given in Table 4 and shown in
Figure 5), following the procedure of Grundahl et al. (2017).

2.3.2. SOPHIE

TOI-1431 was observed with the fiber-fed SOPHIE HR echelle spectrograph on the 1.93 m telescope at the Haute-Provence Observatory (Perruchot et al. 2008, 2011) between 2019 December 18 and 2020 January 12. We collected a total
<figure><img src="image_12.png" /><figcaption>**Figure 4.** Same as Figure 3, but showing the observations taken with AUKR, SCT, ULMT, and LCO.
(The data used to create this figure are available.)</figcaption></figure>

of eight spectra in the HR mode (R = 75,000 at 5500 Å). The second SOPHIE aperture allowed us to check that there was no background sky pollution.

Thorium-argon calibration spectra were taken at the start of each night of the observations to determine the wavelength calibration zero-point, as well as Fabry-Perot étalon spectra regularly during the nights to monitor any possible instrumental drift. The spectra cover the wavelength range from 3872 to 6943 Å across 41 spectral orders, of which 39 were extracted and used for computing the radial velocities. The exposure times ranged from 214 to 490 s, resulting in a uniform S/N for each observation of 50 per resolution element at 5500 Å and median Doppler velocity precision of 7 m s⁻¹.

The spectroscopic data were reduced and then extracted using the standard SOPHIE pipeline (Bouchy et al. 2009), which performs optimal order extraction, cosmic-ray rejection, and wavelength calibration. Radial velocities were computed using the cross-correlation technique and a numerical binary mask corresponding to a star of F0 spectral type (e.g., see Courcol et al. 2015), the closest available binary mask to TOI-1431's spectral type. Table 4 lists the radial velocities from SOPHIE and they are shown in Figure 5. The width of the cross-correlation function (CCF) allows us to measure the rotational velocity of the star as $v \sin i = 6.0 \pm 0.2$ km s⁻¹.

### 2.3.3. FIES

We acquired 52 spectra of TOI-1431 using the Fibre-fed Echelle Spectrograph (FIES; Frandsen & Lindberg 1999; Telting et al. 2014) at the 2.56 m Nordic Optical Telescope (NOT; Djupvik & Andersen 2010) of Roque de los Muchachos Observatory (La Palma, Spain). The observations were carried out between 2019 November 27 and 2020 June 20 (UT). We used the FIES HR mode, which provides a resolving power of R = 67,000 in the spectral range 3760–8220 Å. We traced the RV drift of the instrument by acquiring long-exposure ThAr spectra ($T_{exp} \approx 90$ s) immediately before and after each science exposure. The exposure time was set to 300–900 s, depending on the sky conditions and scheduling constraints.

The data reduction follows the steps of Buchhave et al. (2010) and includes optimal extraction of the spectrum and interpolation of the ThAr wavelength solutions. The S/N per resolution element at 5500 Å of the extracted spectra is in the range 42–127. We cross-correlate each observation order-by-order in velocity space against the strongest exposure of the star (BJD = 2,458,995.703476). Finally, the radial velocity of each observation is measured from the peak of the summed CCF and transformed to the barycentric frame. The uncertainties for the
<figure><img src="image_13.png" /></figure>

**Figure 5.** Radial velocity measurements of TOI-1431b/MASCARA-5 b as a function of orbital phase. The radial velocities from SONG, NRES (ELP and TLV), FIES, SOPHIE, and EXPRESS are represented by the cyan colored circles, purple triangles pointed down, black triangles pointed up, green stars, blue squares, and gold triangles pointed to the left, respectively. A random selection of 20 radial velocity curves drawn from the posterior of the nested sample modeling are plotted as the red lines. The bottom panel shows the residuals between the data and the best-fit (posterior median) model.

derived radial velocities are 7.2–19.1 m s⁻¹ with a mean value of 12.2 m s⁻¹.

### 2.3.4. NRES

We triggered observations of TOI-1431 on the Network of Robotic Echelle Spectrographs (NRES) (Siverd et al. 2018b), operated by the LCOGT (Brown et al. 2013). NRES consists of four similar spectrographs with two in the Northern Hemisphere and two in the Southern. The spectrographs cover the wavelength range from 380 to 860 nm with a resolution of R ~ 53,000. They are fiber-fed, where one fiber observes the science target and the second fiber is used for the wavelength calibration source. In the case of NRES the calibration source is a ThAr lamp.

Due to the target's Northern declination, TOI-1431 was accessible and observed by the units at the Wise Observatory (TLV) and the McDonald Observatory (ELP). We obtained a total of 34 spectra. For spectral classification we used the standard NRES pipeline in combination with SpecMatch-Synth code.<fnref n="52" /> In order to obtain high-precision RV measurements, we processed the spectra using an adapted version of the CERES pipeline (Brahm et al. 2017).

### 2.3.5. EXPRESS

We observed TOI-1431 with the Extreme Precision Spectrometer (EXPRES) (Jurgenson et al. 2016), which was recently commissioned at the 4.3 m Lowell Discovery Telescope (Levine et al. 2012). EXPRESS is a vacuum-stabilized and fiber-fed optical spectrograph, operating between 380 and 780 nm at a resolution of R ~ 137,500. It achieves a single-measurement radial velocity precision of approximately 30 cm s⁻¹ at S/N of 250 per pixel (Petersburg et al. 2020). The wavelength solution is determined from both a ThAr lamp and a Menlo Systems laser frequency comb. Recent performance benchmarks and science results are shown by Blackman et al. (2020) and Brewer et al. (2020), and the radial velocity pipeline is discussed by Petersburg et al. (2020). For TOI-1431, we acquired 15 observations with EXPRESS and achieved a mean radial velocity precision of 2.1 m s⁻¹.

## 2.4. High-contrast Imaging Observations

We obtained high-angular-resolution imaging of TOI-1431 using the NIRC2 instrument on the Keck II. The observations were taken on the night of 2020 September 9 UTC under clear skies, very good seeing conditions of ~0.4, and at an airmass of 1.25. The resulting 5σ detection sensitivity and adaptive optics (AO) image from the NIRC2 observations are plotted in

<footnote n="52">⁵² https://github.com/petigura/specmatch-syn</footnote>
<figure><img src="image_14.png" /><figcaption>**Figure 6.** The 5σ detection sensitivity and AO image from Keck II NIRC2 high-angular-resolution imaging of TOI-1431 in the Brγ band. The orientation of the AO image has north pointed up and east to the left. No stars were detected within 3" of TOI-1431.</figcaption></figure>

Figure 6. No nearby companions or background sources were detected within 3" of TOI-1431.

### 3. Host Star Properties

TOI-1431/MASCARA-5 (HD 201033) is a bright ($V=8.049$) Am (chemically peculiar, kA5mF2) type star (Figueras et al. 1991; Renson et al. 1991; Renson & Manfroid 2009) with a Gaia DR2 parallax of $6.686 \pm 0.046$ mas (Gaia Collaboration et al. 2018) and distance of $149.6_{-1.0}^{+1.1}$ pc. The star has a radius of $1.92 \pm 0.07 R_\odot$, mass of $1.90_{-0.08}^{+0.10} M_\odot$, surface gravity of $\log g = 4.15 \pm 0.04$ dex, effective temperature of $7690_{-250}^{+400}$ K, metallicity of [Fe/H] = $0.43_{-0.28}^{+0.20}$ dex, and luminosity of $11.7_{-1.2}^{+2.3} L_\odot$, derived from an analysis of the broadband SED and the Yonsei Yale (YY) stellar evolutionary models (Yi et al. 2001) using EXOFASTv2 (Eastman et al. 2019) as discussed in Section 3.1. We also derive independent values for the effective temperature and metallicity from the analysis of SONG, SOPHIE, and FIES spectra in Section 3.2.

#### 3.1. SED Analysis

We performed a SED analysis of TOI-1431 using EXOFASTv2 to derive stellar parameters independent of those from spectroscopy. We used broadband photometry obtained from the Tycho (Høg et al. 2000), Two Micron All Sky Survey (2MASS, Cutri et al. 2003), WISE (Cutri et al. 2013), and Gaia (Gaia Collaboration et al. 2018) catalogs for the fitting of the SED (Figure 7). Gaussian priors were placed on the parallax from Gaia DR2, adding 82 µas to correct for the systematic offset found by Stassun & Torres (2018) and adding the 33 µas uncertainty in their offset in quadrature to the Gaia-reported uncertainty. We applied an upper limit on the V-band extinction of 3.4905 from the Schlafly & Finkbeiner (2011) dust maps at the location of TOI-1431 as well as an upper limit on the dilution effects of nearby stars of 0.00484 as reported in the TESS SPOC light curve file. The SED analysis ran until convergence, i.e., the Gelman-Rubin statistic ($R_z$) and the number of independent chain draws ($T_z$) were less than 1.01

<figure><img src="image_15.png" /><figcaption>**Figure 7.** The SED for TOI-1431. The red symbols are the broadband photometric measurements used in the SED analysis (provided in Table 1) with the horizontal uncertainty bars representing the effective width of the passband. The blue symbols are the model fluxes from the best-fit Kurucz atmosphere model.</figcaption></figure>

and greater than 1000, respectively (see Eastman et al. 2019). The YY stellar evolutionary models were used instead of the default MESA isochrones and stellar tracks (MIST) to determine the stellar radius ($R_\star$), mass ($M_\star$), effective temperature ($T_{\text{eff}}$), luminosity ($L_\star$), metallicity ([Fe/H]), surface gravity ($\log g$), and age of TOI-1431, which are given in Table 1. We found that using the MIST models resulted in many of the Markov Chain Monte Carlo (MCMC) steps reaching the [Fe/H] grid's upper limit of 0.5 dex while the MCMC chains never reached the [Fe/H] upper limit of 0.78 dex using the YY tracks. The YY evolutionary models are better suited for hot metal-rich stars such as TOI-1431.

From the best-fit Kurucz stellar atmosphere model from the SED and the best-fitting YY stellar evolutionary model (Figure 8), we find that TOI-1431 is a very hot and metal-rich A-type star with $R_\star = 1.92 \pm 0.07 R_\odot$, $M_\star = 1.90_{-0.08}^{+0.10} M_\odot$, $T_{\text{eff}} = 7690_{-250}^{+400}$ K, $\log g = 4.15 \pm 0.04$ (where g is in units of cm s⁻²), and [Fe/H] = $0.43_{-0.28}^{+0.20}$ dex. The effective temperature found for this star makes it one of the hottest known exoplanet hosts. The only known host stars hotter than TOI-1431 are KELT-9 ($T_{\text{eff}} = 10,170 \pm 450$ K, Gaudi et al. 2017), WASP-178 ($T_{\text{eff}} = 9360 \pm 150$ K, Hellier et al. 2019), KELT-20 ($T_{\text{eff}} = 8720_{-260}^{+250}$ K, Lund et al. 2017), Kepler-1115 ($T_{\text{eff}} = 8480_{-220}^{+300}$ K, Morton et al. 2016), HAT-P-70 ($T_{\text{eff}} = 8450_{-690}^{+540}$ K, Zhou et al. 2019a), HATS-70 ($T_{\text{eff}} = 7930_{-820}^{+630}$ K, Zhou et al. 2019b), and MASCARA-4 ($T_{\text{eff}} = 7800 \pm 200$ K, Dorval et al. 2020). These stellar parameters (in particular $T_{\text{eff}}$, $\log g$, and metallicity), however, are in strong disagreement with the spectroscopic parameters derived from the analysis of the SONG, SOPHIE, and FIES data. The source of this disagreement is likely due to the “anomalous luminosity effect” observed in Am stars (e.g., Bolton 1971) such as TOI-1431, which we discuss further in Section 5.1. The $\log g$ derived from the SED fit is in good agreement with the value obtained from the stellar density ($\log g \sim 4.17$ dex) from fitting the transits, and we have therefore adopted this value along with the values for $T_{\text{eff}}$, $R_\star$, and $M_\star$ as the preferred values for the star.
<table><caption>**Table 1**

Stellar Parameters for HD 201033 as Used in This Work</caption><thead><tr><th>Parameter</th><th>Value</th><th>Source</th></tr></thead><tbody><tr><td>R.A. (hh:mm:ss)</td><td>21:04:48.89</td><td>(1) Gaia DR2</td></tr><tr><td>Decl. (dd:mm:ss)</td><td>55:35:16.88</td><td>(1) Gaia DR2</td></tr><tr><td>μ<sub>α</sub> (mas yr<sup>-1</sup>)</td><td>11.74 ± 0.08</td><td>(1) Gaia DR2</td></tr><tr><td>μ<sub>δ</sub> (mas yr<sup>-1</sup>)</td><td>23.87 ± 0.06</td><td>(1) Gaia DR2</td></tr><tr><td>Parallax (mas)</td><td>6.686 ± 0.046</td><td>(1) Gaia DR2</td></tr><tr><td>Luminosity (L<sub>☉</sub>)</td><td>10.69 ± 0.10<sup>b</sup></td><td>(1) Gaia DR2</td></tr><tr><td>A<sub>V</sub> (mag)</td><td>1.099 (≤3.491)<sup>c</sup></td><td>(2) Schlegel dust maps</td></tr><tr><td>Distance (pc)</td><td>149.6<sup>+1.1</sup><sub>-1.0</sub></td><td>(1) Gaia DR2</td></tr><tr><td>Spectral type</td><td>Am C (kA5mF2)</td><td>(3, 4, 5) A-type star catalogs</td></tr><tr><td colspan="3"><strong>Broadband magnitudes:</strong></td></tr><tr><td>B<sub>T</sub> (mag)</td><td>8.368 ± 0.016</td><td>(6) Tycho</td></tr><tr><td>V<sub>T</sub> (mag)</td><td>8.049 ± 0.011</td><td>(6) Tycho</td></tr><tr><td>TESS (mag)</td><td>7.798 ± 0.006</td><td>(7) TESS TIC v8</td></tr><tr><td>J (mag)</td><td>7.541 ± 0.030</td><td>(8) 2MASS</td></tr><tr><td>H (mag)</td><td>7.452 ± 0.040</td><td>(8) 2MASS</td></tr><tr><td>K<sub>s</sub> (mag)</td><td>7.439 ± 0.030</td><td>(8) 2MASS</td></tr><tr><td>WISE1 (mag)</td><td>7.335 ± 0.034</td><td>(9) WISE</td></tr><tr><td>WISE2 (mag)</td><td>7.430 ± 0.019</td><td>(9) WISE</td></tr><tr><td>WISE3 (mag)</td><td>7.446 ± 0.015</td><td>(9) WISE</td></tr><tr><td>WISE4 (mag)</td><td>7.371 ± 0.096</td><td>(9) WISE</td></tr><tr><td>Gaia (mag)</td><td>7.976 ± 0.001</td><td>(1) Gaia DR2</td></tr><tr><td>GaiaBP (mag)</td><td>8.123 ± 0.001</td><td>(1) Gaia DR2</td></tr><tr><td>GaiaRP (mag)</td><td>7.772 ± 0.001</td><td>(1) Gaia DR2</td></tr><tr><td colspan="3"><strong>Stellar properties from SED and YY tracks:</strong></td></tr><tr><td>T<sub>eff</sub> (K)</td><td>7690<sup>+400a</sup><sub>-250</sub></td><td>(10) EXOFASTv2; this paper</td></tr><tr><td>log g (dex)</td><td>4.15 ± 0.04<sup>b</sup></td><td>(10) EXOFASTv2; this paper</td></tr><tr><td>[Fe/H] (dex)</td><td>0.43<sup>+0.20</sup><sub>-0.28</sub></td><td>(10) EXOFASTv2; this paper</td></tr><tr><td>R<sub>*</sub> (R<sub>☉</sub>)</td><td>1.92 ± 0.07<sup>a b</sup></td><td>(10) EXOFASTv2; this paper</td></tr><tr><td>M<sub>*</sub> (M<sub>☉</sub>)</td><td>1.90<sup>+0.10a</sup><sub>-0.08</sub></td><td>(10) EXOFASTv2; this paper</td></tr><tr><td>ρ<sub>*</sub> (g cm<sup>-3</sup>)</td><td>0.38 ± 0.05</td><td>(10) EXOFASTv2; this paper</td></tr><tr><td>L<sub>*</sub> (L<sub>☉</sub>)</td><td>11.7<sup>+2.3</sup><sub>-1.2</sub></td><td>(10) EXOFASTv2; this paper</td></tr><tr><td>Age (Gyr)</td><td>0.29<sup>+0.32b</sup><sub>-0.19</sub></td><td>(10) EXOFASTv2; this paper</td></tr><tr><td colspan="3"><strong>Spectroscopic properties from SONG spectra:</strong></td></tr><tr><td>T<sub>eff</sub> (K)</td><td>6764 ± 120</td><td>(11, 12) iSpec; this paper</td></tr><tr><td>log g (dex)</td><td>2.76 ± 0.26</td><td>(11, 12) iSpec; this paper</td></tr><tr><td>[M/H] (dex)</td><td>-0.15 ± 0.10</td><td>(11, 12) iSpec; this paper</td></tr><tr><td>[α/H] (dex)</td><td>-0.27 ± 0.10</td><td>(11, 12) iSpec; this paper</td></tr><tr><td>ν sin i (km s<sup>-1</sup>)</td><td>7.1 ± 1.0</td><td>(11, 12) iSpec; this paper</td></tr><tr><td colspan="3"><strong>Spectroscopic properties from SOPHIE spectra:</strong></td></tr><tr><td>T<sub>eff</sub> (K)</td><td>6950 ± 60</td><td>(13) FASMA; this paper</td></tr><tr><td>log g (dex)</td><td>4.72 ± 0.08</td><td>(13) FASMA; this paper</td></tr><tr><td>[Fe/H] (dex)</td><td>0.09 ± 0.03<sup>b</sup></td><td>(13) FASMA; this paper</td></tr><tr><td>ν sin i (km s<sup>-1</sup>)</td><td>6.0 ± 0.2<sup>b</sup></td><td>CCF; this paper</td></tr><tr><td colspan="3"><strong>Spectroscopic properties from FIES spectra:</strong></td></tr><tr><td>T<sub>eff</sub> (K)</td><td>6910 ± 50</td><td>(14) SPC; this paper</td></tr><tr><td>log g (dex)</td><td>3.29 ± 0.10</td><td>(14) SPC; this paper</td></tr><tr><td>[M/H] (dex)</td><td>0.03 ± 0.08</td><td>(14) SPC; this paper</td></tr><tr><td>ν sin i (km s<sup>-1</sup>)</td><td>9.2 ± 0.5</td><td>(14) SPC; this paper</td></tr></tbody></table>

Notes. 1. Gaia Collaboration et al. (2018); 2. Schlaflay & Finkbeiner (2011); 3. Figueras et al. (1991); 4. Renson et al. (1991); 5. Renson & Manfroid (2009); 6. Høg et al. (2000); 7. Stassun et al. (2019); 8. Cutri et al. (2003); 9. Cutri et al. (2013); 10. Eastman et al. (2019); 11. Blanco-Cuaresma et al. (2014); 12. Blanco-Cuaresma (2019); 13. Tsantaki et al. (2018); 14. Buchhave et al. (2012).

a Priors used in the Allesfitter analysis.

b Preferred solution for the stellar and spectroscopic parameters.

c Upper limit on the V-band extinction from Schlegel dust maps.

<figure><img src="image_16.png" /><figcaption>**Figure 8.** The YY mass track for TOI-1431. The black line represents the mass track interpolated at the model values for $M_*$, [Fe/H], and age. The black circle is the best-fit model value for $T_{eff}$ and log g while the red asterisk is the best-fit model age along the track. The black circle and red asterisk almost perfectly overlap, indicating excellent consistency among all the components of the global model.</figcaption></figure>

## 3.2. Spectroscopy

We used the ISPEC spectral analysis tool (Blanco-Cuaresma et al. 2014; Blanco-Cuaresma 2019) to derive the effective temperature, surface gravity, rotational velocity, alpha elemental abundance ([α/H]), and overall metallicity ([M/H]) of TOI-1431 from the reduced and stacked SONG spectra. The ISPEC synthetic grid was configured to incorporate a MARCS atmospheric model (Gustafsson et al. 2008) and the SPECTRUM (Gray & Corbally 1994) radiative transfer code. [M/H] and [α/H] were derived using version 5.0 of Gaia-ESO Survey’s (GES) line list (Heiter et al. 2015) normalized by solar values obtained by Asplund et al. (2009). Initial values for $T_{eff}$, log g, and [M/H] as determined by the SED analysis with EXOFASTv2 (see Section 3.1) were used to construct the synthetic spectral fit. From the ISPEC analysis of the SONG spectra, we find that TOI-1431 has $T_{eff} = 6764 \pm 120$ K, $log g = 2.76 \pm 0.26$ dex, [M/H] = $-0.15 \pm 0.10$ dex, [α/H] = $-0.27 \pm 0.10$ dex, and $ν sin i = 7.1 \pm 1.0$ km s$^{-1}$. These results would seem to suggest that the star is quite evolved, is slightly metal-poor, and is deficient in α elements. However, $T_{eff}$ and log g, in particular, are in very strong disagreement with the results obtained from the SED analysis as a consequence of the nature of the host being an Am star (see discussion in Section 5.1).

The atmospheric stellar parameters $T_{eff}$, log g, and [Fe/H] were derived from the SOPHIE spectra for TOI-1431 using the FASMA spectral synthesis package (Tsantaki et al. 2018). We provided initial stellar parameters for $T_{eff}$, log g, and [Fe/H] from the SED analysis with EXOFASTv2 as a starting point for the spectral synthesis in FASMA. The synthetic spectra were created using the radiative transfer code MOOG (Sneden 1973), and the grid of model atmospheres was generated by the ATLAS-APOGEE (Kurucz 1993), MARCS (Gustafsson et al. 2008), and ATLAS9 (Mészáros et al. 2012) models. The line list used to determine the stellar parameters comes from the Vienna Atomic Line Database (Piskunov et al.
1995; Ryabchikova et al. 2015) and includes all of the iron lines in the regions 5399–5619 Å and 6347–6790 Å as well as all the atomic and molecular lines within intervals of ±2 Å of the iron lines. From this analysis with the SOPHIE spectra, we find the star has $T_{\text{eff}} = 6950 \pm 60$ K, $\log g = 4.72 \pm 0.08$ dex, and $[\text{Fe}/\text{H}] = 0.09 \pm 0.03$ dex, consistent with a late A-type main-sequence star that is modestly metal-rich. We have adopted the [Fe/H] value derived here as the preferred value for the star’s iron abundance.

The stellar parameters for this star were also derived from the FIES spectra using the stellar parameter classification (SPC) technique (Torres et al. 2012; Buchhave et al. 2012, 2014). To derive $T_{\text{eff}}$, $\log g$, $v \sin i$, and [M/H] for this star using SPC, we cross-correlated the reduced and the stacked FIES spectra against a library of synthetic spectra calculated from the Kurucz model atmospheres (Castelli & Kurucz 2003). The synthetic spectra cover the wavelength range 5050–5360 Å; five of the FIES echelle orders span that region and could be used for the spectral analysis. The best-fit values were determined by finding the maximum of the cross-correlation coefficient as a function of the stellar parameters. The results of the SPC analysis on the FIES spectra give $T_{\text{eff}} = 6910 \pm 50$ K, $\log g = 3.29 \pm 0.10$ dex, $[\text{M/H}] = 0.03 \pm 0.08$ dex, and $v \sin i = 9.2 \pm 0.5$ km s⁻¹, which suggest an A-type subgiant star with solar metallicity. It should be noted that SPC was not designed for hot stars such as TOI-1431 and that surface gravities are often poorly constrained by spectral analyses (e.g., see Sozzetti et al. 2007; Winn et al. 2008).

The spectral analyses carried out on the SONG, SOPHIE, and FIES spectra give significantly different results for the host star properties, and in particular for the surface gravity. However, the stellar density (and surface gravity) determined from the SED and transit light-curve analysis are in good agreement with each other. Winn et al. (2008) has demonstrated that stellar densities can be more accurately determined through transit light-curve analysis and should be preferred over those derived through spectroscopy if there are strong disagreements between them. In this light, we therefore have adopted the $T_{\text{eff}}$, $\log g$, $R_*$, and $M_*$ from the SED analysis as the preferred values for the star.

# 4. Modelling and Results

To determine the system parameters for TOI-1431 and its planet, including the planet's thermal emissions and day/night temperature contrast, we used Allesfitter (Günther & Daylan 2019, 2021) to perform a joint analysis of the TESS light-curve segments (from the multiscale MAP PDC, single-scale MAP PDC, and SAP photometry), the photometric ground-based light curves, and the radial velocity measurements. The ground-based light curves included in the analysis are from CDK14 (Sloan *g*' and *z*<sub>s</sub> bands) taken on 2019 December 24, MuSCAT2 (four Sloan bands *g*', *r*', *i*', and *z*<sub>s</sub>) obtained on 2020 May 16, AUKR (Sloan *z*<sub>s</sub> band) observed on 2020 June 16, SCT (TESS-equivalent band) acquired on 2020 August 8, ULMT in the Sloan *z*'<sub>s</sub> band collected on 2020 September 20, and LCOGT captured on 2020 October 14. We used the radial velocity measurements from the SONG, SOPHIE, FIES, NRES, and EXPRESS instruments in the joint analysis.

## 4.1. Joint Transit, Phase-curve, and Radial Velocity Analysis

For the joint analysis, Gaussian priors were placed on the stellar parameters $R_*$, $M_*$, and $T_{\text{eff}}$ from the SED and YY isochrones fitting using EXOFASTv2. We also applied a Gaussian prior on the dilution parameter ($D_0$) from the TESS SPOC light-curve parameter CROWDSAP, after transforming the CROWDSAP value to conform to the definition of the dilution parameter used in Allesfitter ($D_0 = 1 - \text{CROWDSAP}$). For the other physical model parameters, we used uniform priors with reasonable boundaries and starting values. These parameters include the planet-to-star radius ratio ($R_p/R_*$), the ratio of the sum of the planet and star radii to the semimajor axis ($((R_* + R_p)/a_p)$), cosine of the inclination angle ($\cos i_p$), mid-transit time ($T_{0,b}$), orbital period ($P_p$), the radial velocity semiamplitude ($K$), eccentricity ($\sqrt{e_p^2 \cos^2 \omega_p}$ and $\sqrt{e_p^2 \sin^2 \omega_p}$), surface brightness ratio ($J_p$), set of quadratic limb-darkening coefficients ($q_1$ and $q_2$) and flux error scaling ($\ln \sigma_{F_{\text{inst}}}$) for each of the light curves, and a radial velocity baseline offset ($\Delta RV_{\text{inst}}$) and jitter term ($\log \sigma_{RV_{\text{inst}}}$) for each radial velocity instrument.

We experimented with several different detrending techniques in our analysis of each of the three versions of the TESS photometry, including (1) baseline flux offset, (2) rspline detrending prior to the light-curve fitting with a baseline flux offset, (3) hybrid cubic spline detrending run simultaneously with the light-curve fitting, and (4) a red noise detrending model using a Matérn 3/2 Gaussian process (GP) kernel run simultaneously with the light-curve fitting. For the GP parameters, we applied uniform priors with reasonable boundaries and starting values on the flux offset ($GP_{\text{offset}}$, between -0.01 and 0.01), characteristic amplitude ($\ln \sigma_{\text{flux}}$, between -15 and 0), and characteristic length scale ($\ln \rho_{\text{flux}}$, between 2.07 and 12). These parameters were also coupled across the TESS light-curve segments so that a single GP model is sampled at every sampling step. The lower boundary on the characteristic length-scale term was set to three times the planet's orbital period to ensure that the GP preserves any phase-curve modulations while still removing most of the remaining systematics present in the light-curve segments.

Close inspection of the TESS light curves reveals periodic photometric modulations that are synchronized with the orbital period of the planet. These phase-curve variations are driven by the changing viewing angle of the planet and the planet-star gravitational interaction. There are three main components to the phase curve that can be detected in visible-light photometry (see detailed review in Shporer 2017): (1) the atmospheric brightness modulation, which arises due to the variations in atmospheric temperature, thermal emission, and reflected starlight across the surface of a tidally locked planet, (2) ellipsoidal distortion of the host star due to the tidal bulge raised by the orbiting planet, and (3) beaming, i.e., modulations in the star's flux within the bandpass due primarily to the radial velocity-driven Doppler shifting of the stellar spectrum. Setting the zero-point of orbital phase at mid-transit, these three phase-curve modulations appear, to the leading order, as the cosine of the fundamental frequency, the cosine of the first harmonic, and the sine of the fundamental.

In this work, we also fit for these phase-curve terms by including free parameters for the respective full amplitudes: $A_{\text{atmospheric}}$, $A_{\text{ellipsoidal}}$, and $A_{\text{beaming}}$. Similar methodologies have been used in extensive previous phase-curve analyses of TESS-observed exoplanet systems (e.g., Shporer et al. 2019;
<table><caption>**Table 2**
Priors Used in the Analysis</caption><thead><tr><th>Parameter</th><th>Description</th><th>Prior<sup>a</sup></th><th>Value</th></tr></thead><tbody><tr><td colspan="4" style="text-align: center;">Fitted parameters</td></tr><tr><td>R<sub>b</sub>/R<sub>*</sub></td><td>Radius ratio</td><td>U(0.082; 0.010; 0.300)</td><td>0.07955<sup>+0.00063</sup><sub>-0.00053</sub></td></tr><tr><td>(R<sub>*</sub> + R<sub>b</sub>)/a</td><td>Radii sum to semimajor axis</td><td>U(0.21; 0.05; 0.50)</td><td>0.2096 ± 0.0017</td></tr><tr><td>cos i<sub>b</sub></td><td>Cosine of inclination</td><td>U(0.169; 0; 1)</td><td>0.1715 ± 0.0022</td></tr><tr><td>T<sub>0</sub></td><td>Mid-transit time (2,458,700 – BJD<sub>TDB</sub>)</td><td>U(39.177; 38.177; 40.177)</td><td>39.17737 ± 0.00007</td></tr><tr><td>P<sub>b</sub></td><td>Orbital period (d)</td><td>U(2.65; 2.35; 2.95)</td><td>2.650237 ± 0.000003</td></tr><tr><td>√e cos ω</td><td>Eccentricity term</td><td>U(0.0; -0.3; 0.3)</td><td>-0.009<sup>+0.023</sup><sub>-0.021</sub></td></tr><tr><td>√e sin ω</td><td>Eccentricity term</td><td>U(0.0; -0.3; 0.3)</td><td>0.036<sup>+0.033</sup><sub>-0.041</sub></td></tr><tr><td>K<sub>b</sub></td><td>RV semiamplitude (km s<sup>-1</sup>)</td><td>U(0.300; 0; 1.0)</td><td>0.2941 ± 0.0011</td></tr><tr><td>D<sub>0</sub><sub>TESS</sub></td><td>Dilution</td><td>N(0.004904; 0.000588)</td><td>0.00486 ± 0.00051</td></tr><tr><td>q<sub>1</sub><sub>TESS</sub></td><td>Transformed limb darkening</td><td>U(0.40; 0; 1)</td><td>0.22 ± 0.05</td></tr><tr><td>q<sub>2</sub><sub>TESS</sub></td><td>Transformed limb darkening</td><td>U(0.21; 0; 1)</td><td>0.18<sup>+0.13</sup><sub>-0.09</sub></td></tr><tr><td>q<sub>1</sub><sub>LCO Y</sub></td><td>Transformed limb darkening</td><td>U(0.24; 0; 1)</td><td>0.19<sup>+0.13</sup><sub>-0.10</sub></td></tr><tr><td>q<sub>2</sub><sub>LCO Y</sub></td><td>Transformed limb darkening</td><td>U(0.44; 0; 1)</td><td>0.50 ± 0.29</td></tr><tr><td>q<sub>1</sub><sub>MuSCAT2 z'</sub></td><td>Transformed limb darkening</td><td>U(0.35; 0; 1)</td><td>0.26<sup>+0.12</sup><sub>-0.10</sub></td></tr><tr><td>q<sub>2</sub><sub>MuSCAT2 z'</sub></td><td>Transformed limb darkening</td><td>U(0.42; 0; 1)</td><td>0.52 ± 0.28</td></tr><tr><td>q<sub>1</sub><sub>MuSCAT2 r'</sub></td><td>Transformed limb darkening</td><td>U(0.50; 0; 1)</td><td>0.41<sup>+0.13</sup><sub>-0.12</sub></td></tr><tr><td>q<sub>2</sub><sub>MuSCAT2 r'</sub></td><td>Transformed limb darkening</td><td>U(0.47; 0; 1)</td><td>0.51 ± 0.25</td></tr><tr><td>q<sub>1</sub><sub>MuSCAT2 i'</sub></td><td>Transformed limb darkening</td><td>U(0.26; 0; 1)</td><td>0.19<sup>+0.09</sup><sub>-0.08</sub></td></tr><tr><td>q<sub>2</sub><sub>MuSCAT2 i'</sub></td><td>Transformed limb darkening</td><td>U(0.53; 0; 1)</td><td>0.58<sup>+0.24</sup><sub>-0.28</sub></td></tr><tr><td>q<sub>1</sub><sub>MuSCAT2 g'</sub></td><td>Transformed limb darkening</td><td>U(0.62; 0; 1)</td><td>0.54<sup>+0.13</sup><sub>-0.11</sub></td></tr><tr><td>q<sub>2</sub><sub>MuSCAT2 g'</sub></td><td>Transformed limb darkening</td><td>U(0.61; 0; 1)</td><td>0.63<sup>+0.19</sup><sub>-0.21</sub></td></tr><tr><td>q<sub>1</sub><sub>ULMT i'</sub></td><td>Transformed limb darkening</td><td>U(0.26; 0; 1)</td><td>0.14<sup>+0.18</sup><sub>-0.10</sub></td></tr><tr><td>q<sub>2</sub><sub>ULMT i'</sub></td><td>Transformed limb darkening</td><td>U(0.53; 0; 1)</td><td>0.35<sup>+0.32</sup><sub>-0.23</sub></td></tr><tr><td>q<sub>1</sub><sub>CDK14 z'</sub></td><td>Transformed limb darkening</td><td>U(0.35; 0; 1)</td><td>0.54 ± 0.28</td></tr><tr><td>q<sub>2</sub><sub>CDK14 z'</sub></td><td>Transformed limb darkening</td><td>U(0.42; 0; 1)</td><td>0.51 ± 0.30</td></tr><tr><td>q<sub>1</sub><sub>CDK14 g'</sub></td><td>Transformed limb darkening</td><td>U(0.62; 0; 1)</td><td>0.45<sup>+0.29</sup><sub>-0.27</sub></td></tr><tr><td>q<sub>2</sub><sub>CDK14 g'</sub></td><td>Transformed limb darkening</td><td>U(0.61; 0; 1)</td><td>0.63<sup>+0.23</sup><sub>-0.31</sub></td></tr><tr><td>q<sub>1</sub><sub>AUKR z'</sub></td><td>Transformed limb darkening</td><td>U(0.35; 0; 1)</td><td>0.83<sup>+0.11</sup><sub>-0.18</sub></td></tr><tr><td>q<sub>2</sub><sub>AUKR z'</sub></td><td>Transformed limb darkening</td><td>U(0.42; 0; 1)</td><td>0.60<sup>+0.24</sup><sub>-0.29</sub></td></tr><tr><td>q<sub>1</sub><sub>SCTTESS</sub></td><td>Transformed limb darkening</td><td>U(0.40; 0; 1)</td><td>0.22<sup>+0.15</sup><sub>-0.11</sub></td></tr><tr><td>q<sub>2</sub><sub>SCTTESS</sub></td><td>Transformed limb darkening</td><td>U(0.21; 0; 1)</td><td>0.57<sup>+0.26</sup><sub>-0.30</sub></td></tr><tr><td>J<sub>p</sub></td><td>Surface brightness ratio</td><td>U(0.0074; 0; 0.1)</td><td>0.0076 ± 0.0023</td></tr><tr><td>A<sub>beaming</sub></td><td>Phase curve term (ppt)</td><td>U(0.0099; 0; 1)</td><td>0.0010<sup>+0.0012</sup><sub>-0.0007</sub></td></tr><tr><td>A<sub>atmospheric</sub></td><td>Phase curve term (ppt)</td><td>U(0.0976; 0; 1)</td><td>0.0785<sup>+0.0073</sup><sub>-0.0077</sub></td></tr><tr><td>A<sub>ellipsoidal</sub></td><td>Phase curve term (ppt)</td><td>U(0.0409; 0; 1)</td><td>0.0291 ± 0.0064</td></tr><tr><td>ln σ<sub>TESS</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-7.545; -10; -5)</td><td>-7.779<sup>+0.003</sup><sub>-0.004</sub></td></tr><tr><td>ln σ<sub>LCO Y</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-6.53; -10; -3)</td><td>-6.528 ± 0.040</td></tr><tr><td>ln σ<sub>MuSCAT2 z'</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-6.92; -10; -3)</td><td>-6.926 ± 0.041</td></tr><tr><td>ln σ<sub>MuSCAT2 r'</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-7.00; -10; -3)</td><td>-7.007 ± 0.041</td></tr><tr><td>ln σ<sub>MuSCAT2 i'</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-7.02; -10; -3)</td><td>-7.019 ± 0.041</td></tr><tr><td>ln σ<sub>MuSCAT2 g'</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-7.05; -10; -3)</td><td>-7.045 ± 0.042</td></tr><tr><td>ln σ<sub>ULMT i'</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-6.5; -10; -3)</td><td>-5.826 ± 0.030</td></tr><tr><td>ln σ<sub>CDK14 z'</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-5.9; -10; -3)</td><td>-5.330 ± 0.043</td></tr><tr><td>ln σ<sub>CDK14 g'</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-6.1; -10; -3)</td><td>-5.358 ± 0.044</td></tr><tr><td>ln σ<sub>AUKR z'</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-5.8; -10; -3)</td><td>5.808 ± 0.031</td></tr><tr><td>ln σ<sub>SCTTESS</sub></td><td>Flux error scaling (ln rel. flux)</td><td>U(-6.5; -10; -3)</td><td>-5.845 ± 0.023</td></tr><tr><td>log σ<sub>RVSONG</sub></td><td>RV jitter (km s<sup>-1</sup>)</td><td>U(-7.35; -10; -1)</td><td>-7.2 ± 1.6</td></tr><tr><td>log σ<sub>RVELP</sub></td><td>RV jitter (km s<sup>-1</sup>)</td><td>U(-2.93; -10; -1)</td><td>-2.97<sup>+0.34</sup><sub>-0.31</sub></td></tr><tr><td>log σ<sub>RVTLV</sub></td><td>RV jitter (km s<sup>-1</sup>)</td><td>U(-2.79; -10; -1)</td><td>-2.79<sup>+0.14</sup><sub>-0.14</sub></td></tr><tr><td>log σ<sub>RVFIES</sub></td><td>RV jitter (km s<sup>-1</sup>)</td><td>U(-4.43; -10; -1)</td><td>-7.3<sup>+1.3</sup><sub>-1.5</sub></td></tr><tr><td>log σ<sub>RVSOPHIE</sub></td><td>RV jitter (km s<sup>-1</sup>)</td><td>U(-5.09; -10; -1)</td><td>-5.54<sup>+0.73</sup><sub>-1.3</sub></td></tr><tr><td>log σ<sub>RVEXPRES</sub></td><td>RV jitter (km s<sup>-1</sup>)</td><td>U(-2.50; -10; -1)</td><td>-7.15<sup>+0.70</sup><sub>-1.1</sub></td></tr><tr><td>ΔRV<sub>SONG</sub></td><td>RV offset (km s<sup>-1</sup>)</td><td>U(-25.15; -26.0; -24.4)</td><td>-25.154 ± 0.006</td></tr><tr><td>ΔRV<sub>ELP</sub></td><td>RV offset (km s<sup>-1</sup>)</td><td>U(-25.57; -26.5; -24.8)</td><td>-25.579 ± 0.019</td></tr></tbody></table>
<table><caption>**Table 2**
(Continued)</caption><thead><tr><th>Parameter</th><th>Description</th><th>Prior<sup>a</sup></th><th>Value</th></tr></thead><tbody><tr><td>ΔRV<sub>TLV</sub></td><td>RV offset (km s<sup>-1</sup>)</td><td>U(-25.56; -26.0; -25.0)</td><td>-25.570 ± 0.011</td></tr><tr><td>ΔRV<sub>FIES</sub></td><td>RV offset (km s<sup>-1</sup>)</td><td>U(-0.28; -1.2; 0.2)</td><td>-0.274 ± 0.002</td></tr><tr><td>ΔRV<sub>SOPHIE</sub></td><td>RV offset (km s<sup>-1</sup>)</td><td>U(-25.24; -26.0; -24.4)</td><td>-25.249 ± 0.003</td></tr><tr><td>ΔRV<sub>EXPRES</sub></td><td>RV offset (km s<sup>-1</sup>)</td><td>U(0.0; -1.0; 1.0)</td><td>0.0051 ± 0.0007</td></tr><tr><td colspan="4" style="text-align: center;">Derived parameters</td></tr><tr><td>R<sub>*</sub>/a</td><td>Host radius to semimajor axis</td><td></td><td>0.194 ± 0.002</td></tr><tr><td>a/R<sub>*</sub></td><td>Semimajor axis to host radius</td><td></td><td>5.15 ± 0.05</td></tr><tr><td>R<sub>b</sub></td><td>Planet radius (R<sub>⊕</sub>)</td><td></td><td>16.7 ± 0.6</td></tr><tr><td>R<sub>b</sub></td><td>Planet radius (R<sub>J</sub>)</td><td></td><td>1.49 ± 0.05</td></tr><tr><td>a</td><td>Semimajor axis (au)</td><td></td><td>0.046 ± 0.002</td></tr><tr><td>i</td><td>Inclination (deg)</td><td></td><td>80.13 ± 0.13</td></tr><tr><td>b<sub>tra</sub></td><td>Impact parameter</td><td></td><td>0.881 ± 0.004</td></tr><tr><td>e</td><td>Eccentricity</td><td></td><td>0.0022<sup>+0.0030</sup><sub>-0.0016</sub></td></tr><tr><td>ω</td><td>Argument of periastron (deg)</td><td></td><td>108<sup>+92</sup><sub>-31</sub></td></tr><tr><td>q<sub>b</sub></td><td>Mass ratio</td><td></td><td>0.00157 ± 0.00006</td></tr><tr><td>M<sub>b</sub></td><td>Companion mass (M<sub>⊕</sub>)</td><td></td><td>990 ± 60</td></tr><tr><td>M<sub>b</sub></td><td>Companion mass (M<sub>J</sub>)</td><td></td><td>3.12 ± 0.18</td></tr><tr><td>T<sub>tot</sub></td><td>Total transit duration (h)</td><td></td><td>2.489 ± 0.009</td></tr><tr><td>T<sub>full</sub></td><td>Full-transit duration (h)</td><td></td><td>1.057 ± 0.034</td></tr><tr><td>T<sub>0;occ</sub></td><td>Epoch occultation (2,458,700 - BJD<sub>TDB</sub>)</td><td></td><td>40.502 ± 0.002</td></tr><tr><td>b<sub>occ</sub></td><td>Impact parameter occultation</td><td></td><td>0.885<sup>+0.006</sup><sub>-0.005</sub></td></tr><tr><td>ρ<sub>*</sub></td><td>Host density from orbit (cgs)</td><td></td><td>0.37 ± 0.01</td></tr><tr><td>ρ<sub>b</sub></td><td>Companion density (cgs)</td><td></td><td>1.17<sup>+0.18</sup><sub>-0.16</sub></td></tr><tr><td>g<sub>b</sub></td><td>Companion surface gravity (cgs)</td><td></td><td>3430<sup>+50</sup><sub>-60</sub></td></tr><tr><td>T<sub>eq</sub><sup>b</sup></td><td>Equilibrium temperature (K)</td><td></td><td>2370 ± 70</td></tr><tr><td>H<sub>b</sub></td><td>Companion atmospheric scale height (km)</td><td></td><td>230 ± 30</td></tr><tr><td>u<sub>1;TESS</sub></td><td>Limb darkening</td><td></td><td>0.13<sup>+0.13</sup><sub>-0.08</sub></td></tr><tr><td>u<sub>2;TESS</sub></td><td>Limb darkening</td><td></td><td>0.33<sup>+0.10</sup><sub>-0.12</sub></td></tr><tr><td>u<sub>1;LCO Y</sub></td><td>Limb darkening</td><td></td><td>0.40<sup>+0.27</sup><sub>-0.24</sub></td></tr><tr><td>u<sub>2;LCO Y</sub></td><td>Limb darkening</td><td></td><td>0.00<sup>+0.25</sup><sub>-0.22</sub></td></tr><tr><td>u<sub>1;MuSCAT2 z'</sub></td><td>Limb darkening</td><td></td><td>0.51 ± 0.26</td></tr><tr><td>u<sub>2;MuSCAT2 z'</sub></td><td>Limb darkening</td><td></td><td>-0.02<sup>+0.29</sup><sub>-0.24</sub></td></tr><tr><td>u<sub>1;MuSCAT2 r'</sub></td><td>Limb darkening</td><td></td><td>0.64<sup>+0.28</sup><sub>-0.31</sub></td></tr><tr><td>u<sub>2;MuSCAT2 r'</sub></td><td>Limb darkening</td><td></td><td>-0.01<sup>+0.33</sup><sub>-0.28</sub></td></tr><tr><td>u<sub>1;MuSCAT2 i'</sub></td><td>Limb darkening</td><td></td><td>0.49 ± 0.24</td></tr><tr><td>u<sub>2;MuSCAT2 i'</sub></td><td>Limb darkening</td><td></td><td>-0.07<sup>+0.24</sup><sub>-0.20</sub></td></tr><tr><td>u<sub>1;MuSCAT2 g'</sub></td><td>Limb darkening</td><td></td><td>0.91<sup>+0.25</sup><sub>-0.28</sub></td></tr><tr><td>u<sub>2;MuSCAT2 g'</sub></td><td>Limb darkening</td><td></td><td>-0.18<sup>+0.31</sup><sub>-0.26</sub></td></tr><tr><td>u<sub>1;ULMT i'</sub></td><td>Limb darkening</td><td></td><td>0.23<sup>+0.25</sup><sub>-0.16</sub></td></tr><tr><td>u<sub>2;ULMT i'</sub></td><td>Limb darkening</td><td></td><td>0.09<sup>+0.23</sup><sub>-0.19</sub></td></tr><tr><td>u<sub>1;CDK14 z'</sub></td><td>Limb darkening</td><td></td><td>0.67<sup>+0.50</sup><sub>-0.40</sub></td></tr><tr><td>u<sub>2;CDK14 z'</sub></td><td>Limb darkening</td><td></td><td>-0.02<sup>+0.42</sup><sub>-0.39</sub></td></tr><tr><td>u<sub>1;CDK14 g'</sub></td><td>Limb darkening</td><td></td><td>0.76<sup>+0.46</sup><sub>-0.43</sub></td></tr><tr><td>u<sub>2;CDK14 g'</sub></td><td>Limb darkening</td><td></td><td>-0.15<sup>+0.35</sup><sub>-0.33</sub></td></tr><tr><td>u<sub>1;AUKR z'</sub></td><td>Limb darkening</td><td></td><td>1.04<sup>+0.44</sup><sub>-0.50</sub></td></tr><tr><td>u<sub>2;AUKR z'</sub></td><td>Limb darkening</td><td></td><td>-0.17<sup>+0.51</sup><sub>-0.42</sub></td></tr><tr><td>u<sub>1;SCT(TESS)</sub></td><td>Limb darkening</td><td></td><td>0.48<sup>+0.27</sup><sub>-0.24</sub></td></tr><tr><td>u<sub>2;SCT(TESS)</sub></td><td>Limb darkening</td><td></td><td>-0.05<sup>+0.29</sup><sub>-0.21</sub></td></tr><tr><td>δ<sub>TESS<sub>tr;undil</sub></sub></td><td>Transit depth (undil.) (ppt)</td><td></td><td>5.825 ± 0.007</td></tr><tr><td>δ<sub>TESS<sub>tr;dil</sub></sub></td><td>Transit depth (dil.) (ppt)</td><td></td><td>5.797 ± 0.007</td></tr><tr><td>δ<sub>TESS<sub>occ;undil</sub></sub></td><td>Occultation depth (undil.) (ppt)</td><td></td><td>0.127<sup>+0.004</sup><sub>-0.005</sub></td></tr><tr><td>δ<sub>TESS<sub>occ;dil</sub></sub></td><td>Occultation depth (dil.) (ppt)</td><td></td><td>0.127<sup>+0.004</sup><sub>-0.005</sub></td></tr><tr><td>F<sub>TESS<sub>night;undil</sub></sub></td><td>Nightside flux (undil.) (ppt)</td><td></td><td>0.049<sup>+0.004</sup><sub>-0.005</sub></td></tr><tr><td>F<sub>TESS<sub>night;dil</sub></sub></td><td>Nightside flux (dil.) (ppt)</td><td></td><td>0.049<sup>+0.004</sup><sub>-0.005</sub></td></tr><tr><td>δ<sub>LCO Y</sub></td><td>Transit depth (ppt)</td><td></td><td>5.67<sup>+0.25</sup><sub>-0.27</sub></td></tr><tr><td>δ<sub>MuSCAT2 z'</sub></td><td>Transit depth (ppt)</td><td></td><td>5.54<sup>+0.22</sup><sub>-0.24</sub></td></tr><tr><td>δ<sub>MuSCAT2 r'</sub></td><td>Transit depth (ppt)</td><td></td><td>5.30<sup>+0.22</sup><sub>-0.27</sub></td></tr></tbody></table>
<table><caption>**Table 2**
(Continued)</caption><thead><tr><th>Parameter</th><th>Description</th><th>Prior<sup>a</sup></th><th>Value</th></tr></thead><tbody><tr><td>δ<sub>MUSCAT2 i'</sub></td><td>Transit depth (ppt)</td><td></td><td>5.62<sup>+0.21</sup><sub>-0.23</sub></td></tr><tr><td>δ<sub>MUSCAT2 g'</sub></td><td>Transit depth (ppt)</td><td></td><td>4.90 ± 0.26</td></tr><tr><td>δ<sub>ULMT i'</sub></td><td>Transit depth (ppt)</td><td></td><td>5.84<sup>+0.26</sup><sub>-0.31</sub></td></tr><tr><td>δ<sub>CDK14 z'</sub></td><td>Transit depth (ppt)</td><td></td><td>5.13<sup>+0.47</sup><sub>-0.61</sub></td></tr><tr><td>δ<sub>CDK14 g'</sub></td><td>Transit depth (ppt)</td><td></td><td>5.09<sup>+0.55</sup><sub>-0.68</sub></td></tr><tr><td>δ<sub>AUKR z'</sub></td><td>Transit depth (ppt)</td><td></td><td>4.51<sup>+0.53</sup><sub>-0.56</sub></td></tr><tr><td>δ<sub>SCTTESS</sub></td><td>Transit depth (ppt)</td><td></td><td>5.58<sup>+0.24</sup><sub>-0.27</sub></td></tr></tbody></table>

**Notes.** Median values and 68% confidence intervals of the fitted and derived parameters for TOI-1431b from the joint nested sampling Allesfitter analysis of the TESS single-scale MAP corrected PDC photometry, ground-based light curves, and radial velocity measurement.

<sup>a</sup> $\mathcal{N}(\mu; \sigma)$ is a normal distribution with mean $\mu$ and width $\sigma$, $\mathcal{U}(s; a; b)$ is a uniform prior with a starting value $s$ and lower and upper limits of $a$ and $b$, respectively.

<sup>b</sup> Calculated assuming zero Bond albedo.

Wong et al. 2020a, 2020b, 2020c; Daylan et al. 2021). The full list of priors used in the analysis is given in Table 2.

### 4.1.1. Multiscale MAP PDC Photometry

We started the joint analysis using the TESS multiscale MAP corrected PDC light curve (after removing all quality-flagged data and outliers as described in Section 2.1) since most of the instrumental systematics are removed by the SPOC PDC pipeline from the use of cotrending basis vectors and flux corrections due to crowding from other stars (e.g., see Stumpe et al. 2014). This results in a cleaner light curve than the SAP version while still preserving astrophysical signals on time-scales shorter than the orbital period of TOI-1431b, important for phase-curve and secondary eclipse analysis.

The analysis was carried out using an MCMC process that ran until convergence (all chains are >42 times their autocorrelation lengths) for the four different detrending techniques mentioned in Section 4.1. Figure 9 is the resulting phase-folded and detrended PDC light curve of TOI-1431 from the hybrid cubic spline detrending, overplotted with 20 light-curve models drawn from the MCMC posteriors. It is apparent in the residuals to the best-fit transit model that there are two significant dip features in the light curve (~200 ppm depth), one starting ~1 hr before transit ingress and continuing until just after the start of ingress and another starting just prior to the end of egress with a duration of ~1 hr. The residuals from the fit using the other three detrending techniques also show similar dip features. As additional checks, we ran the Allesfitter analysis without any detrending separately on the Sector 15 and 16 PDC light curves and observed similar dip signatures, ruling out the possibility that the light-curve detrending was responsible for these features. These signatures are also present in the DV light curve provided in the DV Report Summary available on exo.MAST.

Further investigations revealed the source of the dip features to be a wavelet artifact caused by the multiscale MAP correction applied by PDC (see Section 2.1). Normally, multiscale MAP performs better than single-scale MAP for removing long-period systematics while preserving signals at transit scales. However, there is a mechanism where multiscale MAP can slightly corrupt a very deep transit in a light curve with little long-period systematics, such as we see in TOI-1431. The single-scale MAP corrected light curve was found to not contain the dip features and, as such, we chose to report the results of the analysis performed on that light curve.

<figure><img src="image_17.png" /></figure>

**Figure 9.** Phase-folded and detrended TESS multiscale MAP PDC light curve of TOI-1431 (Sectors 15 and 16) from the simultaneous hybrid cubic spline detrending model, focusing on the transit. Systematics and other quality-flagged data (as shown in Figure 1) that are the result of, for example, spacecraft pointing anomalies have been removed prior to performing the analysis on the light curve. The blue points are binned at a cadence of 15 minutes and the red solid lines are 20 light-curve models drawn from the posteriors of the Markov Chain Monte Carlo analysis in Allesfitter. Dip features starting ~1 hr before and after transit ingress and egress, respectively, are apparent in the light curve and were introduced during the PDC multiscale MAP cotrending process. These dip features are not present in the single-scale MAP corrected PDC and SAP versions of the light curve (see Figure 10).

### 4.1.2. Single-scale MAP PDC Photometry

The joint analysis with the precleaned (as described in Section 2.1) single-scale MAP corrected PDC light-curve segments in Allesfitter follows the same procedure as with the multiscale MAP corrected PDC light curves. We ran MCMC fits until convergence for the four different detrending techniques. We then used the best-fit parameter values found from the MCMC as starting values in subsequent nested sampling fits to estimate the Bayesian evidence for use in a robust comparison of the four detrending models. We used the default Allesfitter settings for the nested sampling fits (see
<figure><img src="image_18.png" /><figcaption>**Figure 10.** Phase-folded and detrended light curve of TOI-1431 from the simultaneous hybrid cubic spline detrending model, similar to Figure 9, but with the single-scale MAP corrected PDC data instead. The colored points represent the individual phased light-curve segments across Sectors 15 and 16 that have been binned at a cadence of 15 minutes. The red solid lines are 20 light-curve models drawn from the posteriors of the nested sampling analysis in Allesfitter.</figcaption></figure>

(The data used to create this figure are available.)

Günther & Daylan 2021) that ran until they reached the convergence criterion threshold of $\Delta \ln Z \le 0.01$.

The results from the nested sampling fits are in good agreement with the MCMC results to within $1\sigma$. Comparing the Bayesian evidence of the nested sampling models, we find that the hybrid cubic spline detrending run simultaneously with the light-curve fitting is decisively favored over the second most favored model, the red noise GP detrending model, with a Bayes model comparison factor of $\ln R = 119$. We also applied the hybrid cubic spline detrending to each of the ground-based transit light curves used in the joint fit.

Figure 10 shows the resulting phase-folded and detrended single-scale MAP PDC light curve of TOI-1431 from the simultaneous cubic spline detrending nested sampling model. Unlike the multiscale MAP PDC light curve shown in Figure 9, there is no evidence of the strange dip features near the transit, motivating our decision to use the results derived from the single-scale MAP version of light-curve analysis.

Both the MCMC and nested sampling give consistent results, and we report both the fitted and derived parameter results of the nested sampling model in Table 2. In Figures 3, 4, and 5, we show the ground-based transit light curves and radial velocity measurements with the corresponding models drawn from the nested sampling posterior, respectively.

Our preferred solution from the joint analysis using the single-scale MAP PDC light curve shows that TOI-1431 hosts an inflated and highly irradiated Jovian planet with a radius of $R_p = 1.49 \pm 0.05 R_j$ (16.7 $\pm 0.6 R_⊕$), mass of $3.12 \pm 0.18 M_J$, high equilibrium temperature of $T_{eq} = 2370 \pm 70$ K (assuming zero Bond albedo), and a circular orbit with a period of $P = 2.650237 \pm 0.000003$ days. The mean density of the planet

is $1.17^{+0.18}_{-0.16}$ g cm$^{-3}$, slightly less than but consistent with the density of Jupiter ($1.326$ g cm$^{-3}$) and similar to other Jovian worlds. The large number and high quality of radial velocity measurements (as shown in Figure 5), coupled with the host star's low $v \sin i$ of $6.0 \pm 0.2$ km s$^{-1}$, have allowed us to measure the radial velocity semiamplitude of the orbit ($294.1 \pm 1.1$ m s$^{-1}$) with unusually high precision for a planet orbiting such a hot star ($7690^{+400}_{-250}$ K). In fact, this is the highest precision orbit measured for any planet that has a host star hotter than 6600 K. The largest contributor to the uncertainty in the planet mass ($17.4\sigma$ detection) is from the uncertainty on the host star mass.

There is also a clear signature of the secondary eclipse (occultation, i.e., when the planet passes behind the star) with a depth of $\delta_{occ;dil} = 0.127^{+0.004}_{-0.005}$ ppt. We also robustly detect the longitudinal modulation of the planet's light in the TESS bandpass ($A_{atmospheric} = 0.078^{+0.007}_{-0.008}$ ppt) from the full phase curve, as shown in Figure 11, providing us with a rare opportunity to measure the planet's day/night brightness temperature contrast (see Section 4.1).

From the light curve, we see a very weak beaming signal with a measured amplitude of $A_{beaming} = 0.0010^{+0.0012}_{-0.0006}$ ppt and a marginal detection of ellipsoidal modulation with an amplitude of $A_{ellipsoidal} = 0.029 \pm 0.006$ ppt. The predicted full amplitudes of both of these phase-curve components can be derived from theory and depend on the planet-star mass ratio $q_p = M_p/M_*$ (see, for example, the review in Shporer 2017):

$$A_{\text{ellipsoidal}} = 2\alpha_{\text{ellip}}q_p\left(\frac{R_*}{a}\right)^3 \sin^2 i_p, \quad (1)$$

$$A_{\text{beaming}} = 2 \left[ \frac{2\pi G q_p^2 M_p \sin^3 i_p}{P_p c^3 (1 + q_p)^2} \right]^{1/3} \left\langle \frac{xe^x}{e^x - 1} \right\rangle_{\text{TESS}}. \quad (2)$$

Here, we have assumed for simplicity that the stellar spectrum is a blackbody, and $x \equiv hc/k\lambda T_*$; the angled brackets indicate averaging over the TESS bandpass. The prefactor $\alpha_{ellip}$ in the ellipsoidal distortion is an expression that includes the limb- and gravity-darkening coefficients for the host star.

To calculate the predicted amplitudes, we use the stellar parameters derived from the SED fit (Table 1) and the best-fit system parameters from the joint analysis. For the limb- and gravity-darkening coefficients, we take tabulated values from Claret (2017) and interpolate to the measured stellar parameters. Uncertainties are propagated to the estimates using Monte Carlo sampling. We obtain predicted amplitudes of $A_{ellipsoidal} = 0.028 \pm 0.004$ ppt and $A_{beaming} = 0.0052 \pm 0.0004$ ppt. These theoretical values are in excellent agreement with our measured values to within the $1\sigma$ level.

Our phase-curve analysis did not consider a possible phase shift in the time of maximum brightness, which would correspond to an offset in the dayside hotspot. Previous studies of ultrahot Jupiters have indicated very small (if any) phase-curve offsets (see, for example, the case of KELT-9 b; Wong et al. 2020b). Any offset would manifest itself in our phase-curve fits through its effect on the measured beaming amplitude. We have shown that the measured and predicted values of $A_{beaming}$ are consistent to within $1\sigma$. Therefore, we conclude that there is no evidence for a significant phase offset in the planet's atmospheric brightness modulation.
<figure><img src="image_19.png" /><figcaption>Figure 11. The phase-folded TESS single-scale MAP PDC light curve for TOI-1431b, zoomed in to show the phase curve and secondary eclipse. The gray points are the detrended photometry, the dark gray points are the binned photometry with a cadence of P/300 (~12.7 minutes), and the black stars are the binned photometry with a cadence of P/150 (~25.5 minutes). Included in the phase-curve model fit are the ellipsoidal (blue dashed-dotted line), the atmospheric (orange dashed line), and the beaming (green dotted line) modulations. The combined median model fit is plotted as the solid red line. The dark gray dashed line represents no emission from TOI-1431b. The bottom panel shows the residuals to the combined median model.</figcaption></figure>

### 4.1.3. SAP Photometry

We also performed the analysis of the precleaned SAP light-curve segments following the same procedures as done in Sections 4.1.1 and 4.1.2. The SAP light curve (see Figure 15 in the Appendix) also shows no evidence of the dip features near the transit.

We do find some discrepancies in the results between the joint fit using the SAP light curve and the single-scale MAP PDC light curve that are likely due to spacecraft point excursions (see Figure 16 in the Appendix).

However, most of our derived parameters were consistent within $1\sigma$. We adopt the results from the joint analysis using the single-scale MAP PDC light curve based on Bayesian model comparison for the discussion of this paper.

### 4.2. The Atmosphere

We have used the TESS single-scale MAP corrected PDC photometry and Allesfitter to characterize the red-optical (6000–9500 Å) phase curve of TOI-1431b, which includes contributions from beaming, ellipsoidal modulation, and atmospheric brightness modulation, as well as the secondary eclipse (see Section 4.1.2). We also fit the phase curve and secondary eclipse using the TESS SAP photometry (see Figure 18 in the Appendix) as a comparison to our preferred solution using the single-scale MAP PDC photometry.

From the measured secondary eclipse depth $\delta_{occ}$ and nightside planetary flux $\delta_{night}$, we can derive the dayside and nightside brightness temperatures of TOI-1431b using the relationship between the planet-star flux ratio $\delta$, planetary emission $F_p(\lambda, T_p)$ (assumed to be a blackbody), and stellar

spectrum $F_*(\lambda, T_{eff})$:

$$ \delta = \left( \frac{R_p}{R_*} \right)^2 \frac{\int F_p(\lambda, T_p) \tau(\lambda) d\lambda}{\int F_*(\lambda, T_{eff}) \tau(\lambda) d\lambda} + A_g \left( \frac{R_p}{a} \right)^2. \quad (3) $$

Any light reflected by the dayside atmosphere contributes to the measured flux ratio through the geometric albedo $A_g$; this second term is zero when computing the nightside temperature. The planetary and stellar spectra are integrated over the TESS bandpass, which has a transmission function (in energy units) represented by $\tau(\lambda)$.

We model the star’s spectrum using PHOENIX models (Husser et al. 2013). To properly interpolate the stellar models and propagate the uncertainties in stellar parameters to the estimates of planetary temperature, we follow the methodology described in Wong et al. (2020b) and calculate the integrated stellar flux within the TESS bandpass for a grid of PHOENIX models, before fitting a cubic polynomial in ($T_{eff}$, log g, [Fe/H]) to the full set of values. We then use a standard MCMC routine to compute the posterior distribution of the planetary brightness temperature $T_p$, with Gaussian priors on the stellar and system parameters, as well as the measured dayside or nightside flux ($\delta_{occ}$ or $\delta_{night}$; see Table 2).

For the dayside brightness temperature, we obtain $T_{day} = 3004 \pm 64$ K when assuming zero geometric albedo; across an albedo range of 0–0.2, we find dayside brightness temperatures spanning 2700–3100 K. We measure a very high nightside brightness temperature of $T_{night} = 2583 \pm 63$ K. The extremely hot dayside is expected to preclude the formation of condensates, resulting in a cloud-free dayside hemisphere (e.g., Wakeford et al. 2017). Previous analyses of other ultrahot Jupiters combining secondary eclipses at optical and thermal
<figure><img src="image_20.png" /><figcaption>**Figure 12.** Plots of nightside brightness temperature $T_{\text{night}}$, Bond albedo $A_B$, and day–night heat recirculation efficiency $\epsilon$ as a function of equilibrium temperature $T_{\text{eq}}$ for a population of hot and ultrahot Jupiters. TOI-1431b is indicated in blue; the black points are taken/derived from the results of Bell et al. (2021). KELT-16b has been omitted due to its very large nightside temperature uncertainty (>400 K) and poorly constrained $\epsilon$ value. With its high nightside temperature, low day/night temperature contrast, and strongly negative Bond albedo, TOI-1431b is a notable outlier among gas giants with similar levels of stellar irradiation, and is very similar to HAT-P-7b.</figcaption></figure>

infrared wavelengths confirm the predictions of low reflectiv-
ity, yielding near-zero geometric albedos (e.g., WASP-18b and
WASP-33b; Shporer et al. 2019; von Essen et al. 2020).

We use the measured dayside and nightside brightness temperatures to jointly constrain the Bond albedo $A_B$ and the efficiency of heat recirculation from the dayside to nightside of the planet $\epsilon$ (0 for no recirculation; 1 for full recirculation), following the thermal balance formalism outlined in Cowan & Agol (2011) and the methodology described in Wong et al. (2020b). We find $A_B = -0.82_{-0.38}^{+0.30}$ and $\epsilon = 0.76 \pm 0.05$. The very efficient day-night heat recirculation is borne out by the low day/night temperature contrast (~420 K). However, the unusual negative Bond albedo suggests that the overall thermal emission level across the entire planet cannot be accounted for by the energy input from stellar insolation alone.

By placing TOI-1431b in context, we can appreciate the
exceptional atmospheric properties of this planet. Figure 12
shows the nightside brightness temperatures, Bond albedos,
and recirculation efficiencies for a sample of hot and ultrahot
Jupiters, taken and/or analogously derived from the brightness
temperatures listed in the comprehensive Spitzer 4.5 µm phase-
curve analysis by Bell et al. (2021). From the plots, it is evident
that TOI-1431b has a much higher nightside temperature and
day-night heat recirculation efficiency than most other hot
Jupiters with comparable equilibrium temperatures. In part-
icular, this planet has a higher heat recirculation efficiency and
significantly lower Bond albedo than KELT-9 b ($\epsilon = 0.39 \pm$
0.05 and $A_B = 0.19_{-0.11}^{+0.12}$; Wong et al. 2020b), the hottest
exoplanet ever discovered (Gaudi et al. 2017), with dayside and
nightside brightness temperatures of $T_{\text{day}} = 4450_{-210}^{+220}$ K and
$T_{\text{night}} = 3290 \pm 170$ K (Bell et al. 2021), respectively. Mean-
while, only HAT-P-7 b has a similar strongly (>2σ) negative
inferred Bond albedo; indeed, HAT-P-7 b is very similar to
TOI-1431b in all respects. Broadly speaking, KELT-9 b, HAT-
P-7 b, and TOI-1431b are outliers amid the weak overall trends
across gas giants with $T_{\text{eq}} > 1200$ K, i.e., increasing nightside
temperature, decreasing Bond albedo, and decreasing day-
night heat recirculation efficiency with increasing equilibrium
temperature. In the case of KELT-9 b in particular, given its
extreme dayside temperature, its anomalously high heat
recirculation efficiency could be explained from the moderation
of the day/night temperature contrast due to the transport of
thermally dissociated atomic hydrogen from the dayside to the
nightside and subsequent recombination into molecular

hydrogen, which releases a significant amount of heat (e.g.,
Bell & Cowan 2018).

The unexpected negative Bond albedo may indicate the
limitations of the simple thermal balance considerations
underpinning our estimates of A<sub>B</sub> and ϵ. In particular, gradients
in chemical composition between the dayside and nightside
hemispheres can entail drastically different atmospheric
opacities, meaning that the pressure levels probed by our
broadband photometric measurements may vary significantly
across the planet's surface. Another explanation for a negative
Bond albedo is additional thermal emission from residual heat
of formation. This scenario would increase the emitted flux at
all longitudes, raising the measured brightness temperatures on
both the dayside and the nightside hemispheres.

This explanation may be especially applicable to TOI-1431b: the stellar age inferred from SED modeling is $0.29_{-0.19}^{+0.32}$ Gyr, making the system among the youngest hosting giant planets hitherto discovered. However, the age of the planet is sufficient for it to have deflated through the initial Kelvin–Helmholtz contraction phase and for its atmosphere to have cooled and reached the equilibrium temperature. The timescale over which this initial cooling and deflation occurred is ~1 Myr (see Equation (17) from Ginzburg & Sari 2015), the photospheric cooling time to reach the regime where stellar irradiation acts to slow cooling. The current level of incident stellar flux is sufficient to slow the planet’s interior cooling and allow it to remain hot and inflated at its present-day radius (see Figure 17 in the Appendix and Komacek et al. 2020; Mol Louis & Miguel 2020). We do not require deposited heating in the deep interior to explain the present-day radius, however; a weak conversion of <0.05% of the incident stellar flux to deposited heating is allowed. Importantly, we find that the effects of irradiation slowing cooling can only explain the present-day radius if TOI-1431b arrived at its current orbit within ~1 Myr after formation. This may imply either that TOI-1431b formed in situ or that a nearby stellar or massive planetary companion rapidly scattered the planet inward very soon after formation. The high-obliquity orbit as measured by the Rossiter–McLaughlin effect (Stangret et al. 2021, see Section 5.2) suggests that the planet experienced just such a scattering event early in its history.

Future spectrally resolved emission spectra with the James
Webb Space Telescope (JWST) can enable detailed analyses of
the atmospheric composition and temperature–pressure profiles
<figure><img src="image_21.png" /><figcaption>**Figure 13.** Comparison of TOI-1431b and its host star with other hot and ultrahot Jupiter systems. Top left: the atmospheric scale heights of hot Jupiters as a function of host star brightness in $J_{mag}$. The points are colored based on the distance from the host star, with TOI-1431b plotted as a triangle. This planet transits one of the brightest stars, making it an excellent candidate for future follow-up atmospheric studies. Top right: planet mass versus radius of hot Jupiters. The red point is the location of TOI-1431b. Planets located above the dark gray solid line are predicted to have no solid core, while those plotted below the light gray solid line are predicted to have a solid core of at least 100 $M_{⊕}$ (Fortney et al. 2007). The dotted lines are isodensity lines, with the density value given in the top left, in cgs units. This plot shows that TOI-1431b is moderately inflated compared to other Jovians of similar mass. Bottom left: incident stellar irradiation versus planet radius for all known Jovians and sub-Jovians that have measured masses and radii with the planet masses indicated by the color scale on the left. TOI-1431b is plotted as the colored triangle on the right and it is clear that this planet is one of the most highly irradiated planets. Bottom right: plot of stellar effective temperature versus surface gravity for planet-hosting stars. TOI-1431 is plotted as the red point and is seen to be one of the hottest (top 1%) known planet-hosting stars.</figcaption></figure>

across the planet's surface, providing a full picture of the thermal energy budget for TOI-1431b.

From the derived equilibrium temperature and surface gravity of the planet and assuming an H/He atmosphere with a mean molecular mass of $\mu = 2.3$ amu, the calculated atmospheric scale height (from $H_p = kT_{eq}/(\mu g_p)$) is $H_p = 220 \pm 30$ km. While this is not the largest scale height among the population of hot Jupiters, this planet orbits one of the brightest host stars ($J_{mag} = 7.541 \pm 0.030$ and $K_{mag} = 7.439 \pm 0.030$, see Figure 13), making it a potential target for future atmospheric characterization through transmission spectroscopy.

The transmission spectroscopy metric (TSM, see Kempton et al. 2018), used to assess the suitability of transmission

spectroscopy observations with JWST, is ~110 for TOI-1431b. Planets with TSM values greater than 90 (for Jovians and sub-Jovians), such as for this planet, are considered suitable for these observations with JWST. However, transmission spectroscopy carried out by Stangret et al. (2021) from two HARPS-N and one EXPRES transit observations finds no absorption signatures in the planet's atmosphere, likely due to its high surface gravity. Additionally, given the detection of the phase curve and secondary eclipse in the TESS red-optical band photometry, this will be an excellent target for emission spectroscopy with JWST. In particular, phase-curve observations carried out with the Near Infrared Imager and Slitless Spectrograph (NIRISS) over wavelengths between 0.6 and
5.0 µm should provide a high-precision global temperature map of this planet's atmosphere (Parmentier & Crossfield 2018). Measurements of the abundances of molecular species as a function of longitude (chemical mapping) could also be probed through phase-resolved spectroscopic observations taken with JWST's NIRCam and NIRSpec instruments.

Previous atmospheric modeling of the extreme end-member ultrahot Jupiter KELT-9 b can serve as a guide for exploring what future intensive emission spectroscopy might reveal for TOI-1431b. The comparison is particularly apropos due to the similar planetary surface gravities of the two planets —log g<sub>p</sub> = 3.30 and 3.54 (in cgs units) for KELT-9 b and TOI-1431b, respectively. In Wong et al. (2020b), radiative transfer calculations of the dayside emission spectrum of KELT-9 b produced largely featureless spectra across the visible and near-infrared wavelength range. In contrast, model nightside emission spectra spanning 2500–3500 K showed large near-infrared absorption features due to H<sub>2</sub>O and CO, as well as excess continuum opacity in the visible and near-infrared through ~1.7 µm due to dissociated H<sup>-</sup>; this latter spectral feature is a unique property of the hottest ultrahot Jupiters and has been studied in detail by numerous earlier theoretical works (e.g., Arcangeli et al. 2018; Lothringer et al. 2018; Parmentier et al. 2018). Measuring the detailed shape of this H<sup>-</sup> feature can provide the abundance of H<sup>-</sup> and the temperature–pressure profile across the nightside, which will in turn probe the effect of hydrogen dissociation on the global atmospheric dynamics and thermal energy budget.

TOI-1431b will also be a great target for detailed atmo-
spheric characterization by the European Space Agency's
Atmospheric Remote-sensing Infrared Exoplanet Large-survey
(ARIEL, Tinetti et al. 2018) telescope. ARIEL will operate in
the infrared with a spectral range of 1.25–7.8 µm as well as
multiple narrowband photometry in the optical. As such, it will
be well suited to potentially measure this planet's equilibrium
chemistry, trace gases, and vertical and horizontal thermal
structures, and to detect clouds and cloud composition
(assuming the atmosphere is not cloud-free) through transmis-
sion and emission spectroscopy and phase-curve observations.

5. Discussion

TOI-1431b joins the growing list of ~16 hot and ultrahot Jovians with measured full phase curves and secondary eclipses (e.g., see Bell et al. 2021). With dayside and nightside brightness temperatures of T<sub>day</sub> = 3004 ± 64 K and T<sub>night</sub> = 2583 ± 63 K, respectively, TOI-1431b is one of the hottest known exoplanets.

We have measured the mass of TOI-1431b to be M<sub>p</sub> = 3.12 ± 0.18 M<sub>J</sub> from radial velocity observations obtained from the SONG, SOPHIE, FIES, NRES TLV and ELP, and EXPRES spectrographs. When combined with the planet radius of R<sub>p</sub> = 1.49 ± 0.05 R<sub>J</sub>, this mass gives a bulk density of ρ<sub>p</sub> = 1.17<sup>+0.18</sup><sub>-0.16</sub> g cm<sup>-3</sup>, similar to Jupiter (1.326 g cm<sup>-3</sup>). The measured radial velocity semiamplitude of K = 294.1 ± 1.1 m s<sup>-1</sup> is the most precise K measurement for a planet-hosting star hotter than 6600 K. A large contributing factor for this is the unusually low projected stellar rotational velocity of v sin i < 10 km s<sup>-1</sup>.

Figure 13 shows a plot of the planet radius as a function of
planet mass as well as a plot of planet radius as a function of
stellar irradiation for all known Jovians and sub-Jovians with
measured masses and radii. The radius of TOI-1431b is inflated

compared to other exoplanets of similar mass as shown in
Figure 13, but is not so inflated when accounting for the very
high insolation flux (5300<sup>+500</sup><sub>-470</sub> S<sub>⊕</sub>), indicating that the planet
inflation is likely the result of stellar irradiation and youthful
age (0.29<sup>+0.32</sup><sub>-0.19</sub> Gyr).

5.1. Metal-peculiar Am Star

TOI-1431b orbits a bright (V<sub>T</sub> ~ 8.0 and J ~ 7.5) and very
hot (T<sub>eff</sub> ~ 7700 K) host star. In fact, TOI-1431 is one of the
hottest planet-hosting stars as shown in Figure 13. The star is
also classified as being a nonmagnetic metallic-line chemically
peculiar Am star (see Conti 1970; Figueras et al. 1991;
Renson et al. 1991; Renson & Manfroid 2009). Am stars
typically have slow rotation rates compared to typical "field" A
stars of similar effective temperatures (A star v sin i range
between ~100 and ~150 km s<sup>-1</sup>; see, for example, Abt &
Morrell 1995), which is certainly the case for TOI-1431
(assuming the star is not being observed nearly pole-on). The
slow rotation of these stars is believed to be due to a close
orbiting (2.5 &le; P &le; 100 days) stellar companion (at least 64%
of Am stars are found to be in binary systems; see, e.g., Böhm-
Vitense 2006; Carquillat & Prieur 2007) that raises tidal forces
on the primary star, causing the primary to lose angular
momentum and spin down (Michaud et al. 1983).

This in turn results in the onset of gravitational settling and radiative levitation due to the lack of mixing in the shallow convective envelopes of slowly rotating A stars that produces a chemical peculiarity observed as a photospheric overabundance of iron-peak elements (such as Cr, Mn, Fe, and Ni and specifically of Ba, Y, and Sr) but depleted abundances of light elements such as Ca, Sc, and Mg (e.g., see Preston 1974; Michaud et al. 1983; Charbonneau 1993; Böhm-Vitense 2006; Xiang et al. 2020).

The spectroscopic analysis of TOI-1431 from spectra
collected from SONG reveals that the star has an overall
metallicity slightly lower than the solar abundance ([M/H] = -0.15 ± 0.10 dex) and even less so in alpha elements
(such as O, Mg, Si, S, Ca, and Ti) with [α/H] = -0.27 ± 0.10
dex, while analysis of the SOPHIE spectra gives a somewhat
higher than solar abundance of iron ([Fe/H] = 0.09 ± 0.03).
This seems to suggest that the star has an overabundance of
iron-group elements and an underabundance of light elements
(compared with normal A stars) that is characteristic of Am
stars. However, a more detailed analysis of the stellar spectra is
required to confirm the Am star classification.

Another interesting feature of Am stars is the "anomalous luminosity effect", which causes deviations in the strengths of specific spectral lines from what is expected from normal main-sequence A stars at a given effective temperature (Bolton 1971). This means that, depending on the spectral lines used when fitting the stellar spectrum model, one can obtain very different values for a star's effective temperature and surface gravity, hence the overall stellar classification. This appears to be the case for TOI-1431 (especially for log g); the analysis of SONG spectra gives T<sub>eff</sub> = 6764 ± 120 K and log g = 2.76 ± 0.26 dex, the FIES spectral analysis gives T<sub>eff</sub> = 6910 ± 50 K and log g = 3.29 ± 0.10 dex, and the SOPHIE spectral analysis gives T<sub>eff</sub> = 6950 ± 60 K and log g = 4.72 ± 0.08 dex.

Planets orbiting main-sequence Am stars appear to be quite
rare—only five that transit have been discovered to date: TOI-
1431b, WASP-33 b (Collier Cameron et al. 2010b), KELT-17 b
(Zhou et al. 2016; Saffe et al. 2020, 2021), KELT-19A b
(Siverd et al. 2018c), and WASP-178 b (Hellier et al. 2019)/ KELT-26 b (Rodríguez Martínez et al. 2020). This raises some interesting questions. First, what is the occurrence rate of planets around Am stars? Radial velocity surveys have historically avoided hot and early type stars ($T_{\text{eff}} > 6500$ K and mid-F and earlier), and instead have targeted "solar analogs" (i.e., late-F, G, and K type stars, see Cumming et al. 1999; Tinney et al. 2001; Pepe et al. 2004; Valenti & Fischer 2005; Wright et al. 2008) in the search for planets. As such, planets orbiting hot main-sequence A stars like TOI-1431 have gone nearly undiscovered, until now thanks to TESS indiscriminately surveying the sky for transiting planets orbiting bright stars. Planets have been discovered orbiting former A stars, i.e., stars that have evolved off the main sequence to become giant and subgiant stars, in radial velocity surveys targeting such stars (e.g., Johnson et al. 2006; Trifonov et al. 2014; Reffert et al. 2015; Jones et al. 2016; Luhn et al. 2019; Wittenmyer et al. 2020). Recent results from surveys of evolved stars indicate that the occurrence rate of giant planets around giant stars is ~10% (e.g., Wittenmyer et al. 2020; Wolthoff et al., under review), indicating that giant planets are relatively common around hot stars and could also be common around Am stars as well.

If giant planet formation around A-type stars is indeed a relatively common occurrence, could the migration of a Jovian planet to a close-in orbit around its host (as presumed to be the case for TOI-1431b) play a significant role in its tidal spin-down and, as such, contribute to its nature as being an Am star? For TOI-1431b, we do not see any conclusive evidence for this given that the system does not appear to be tidally synchronized ($P_p \approx P_{\text{rot}}$) based on the host star's rotation period of between ~10 and ~16 days (from spectroscopic measurements of $v \sin i$) and the planet's orbital period of ~2.65 days. However, the one caveat is that the orientation of the stellar spin axis ($I_\star$) is unknown and the star could perhaps be spinning more rapidly if we are observing it nearly pole-on (leaving open the possibility of $P_p \approx P_{\text{rot}}$). If indeed giant planets can contribute to the nature of their host's being an Am star, for Am stars without evidence of a stellar binary, do most or all host a hot sub-stellar or massive planetary companion? Radial velocity and transit searches targeting Am stars could resolve this question as well as provide valuable insights into the formation and evolution of Am stars, a process that is not fully understood.

## 5.2. Projected Spin-Orbit Angle

We attempted to measure the projected spin-orbit alignment ($\lambda$) of this system, using the planetary shadow technique (e.g., see Collier Cameron et al. 2010a; Johnson et al. 2014; Zhou et al. 2016), by acquiring in-transit HR spectroscopic observations with the FIES instrument. We observed a transit of TOI-1431b on the night of 2020 May 23 using FIES, obtaining a total of 30 exposures, each with an exposure time of 300 s and with airmass decreasing from 1.88 to 1.16 throughout the observation. The wavelength-calibrated 2D reduced spectra were extracted as outlined in Section 2.3.3. We then followed a procedure similar to that of Hoeijmakers et al. (2020) to attempt to extract the absorption line deformations during the transit of the planet. This involved first cross-correlating a continuum-normalized model template of the stellar spectrum of TOI-1431 with the telluric-corrected FIES spectra (in the stellar rest frame) to produce in-transit CCFs.

The in-transit CCFs were then divided by the mean out-of-transit CCF to retrieve the Doppler shadow map.

From our analysis, the planetary shadow is not clearly detected from the FIES data, likely because of the slow rotation of the star. However, observations taken with HARPS-N and EXPRES of three transits of TOI-1431b do successfully detect the Rossiter-McLaughlin effect, revealing that the planet is on a retrograde orbit with $\lambda = -155.3_{-11.3}^{+16.9} (\text{Stangret et al. 2021})$. The results of Stangret et al. (2021) suggest that TOI-1431b likely experienced high-eccentricity migration in the past that produced its high-obliquity orbit and then later the orbit was tidally circularized to the period of ~2.65 days we observe today (for dissenting views of the formation of close-in gas giant planets in situ via the core-accretion process and spin-orbit misalignments due to processes unrelated to planet migration, see, e.g., Batygin et al. 2016; Hasegawa et al. 2019; Louden et al. 2021; Hjorth et al. 2021). This result follows the general pattern observed by other studies that the hottest stars tend to host planets on misaligned orbits (e.g., Winn et al. 2009; Albrecht et al. 2012; Louden et al. 2021).

## 6. Conclusions

We have presented the discovery of the transiting ultrahot Jupiter TOI-1431b. This planet orbits one of the hottest ($T_{\text{eff}} = 7690_{-250}^{+400}$ K) and brightest ($V_T \sim 8.0$) of the known host stars with a period of just $P = 2.650237 \pm 0.000003$ days, resulting in it receiving a high amount of insolation flux and being moderately inflated. A joint analysis of the TESS light curve, ground-based light curves from MuSCAT2 and LCOGT, and radial velocities from SONG, SOPHIE, FIES, NRES, and EXPRES instruments results in a planet radius of $R_p = 1.49 \pm 0.05 R_J$ ($16.7 \pm 0.6 R_\oplus$) and a planet mass of $3.12 \pm 0.18 M_J$, corresponding to a bulk density of $1.17_{-0.16}^{+0.18} \text{ g cm}^{-3}$.

The planet's phase curve and secondary eclipse have been detected from the TESS photometry, providing us with the exciting opportunity to measure the planet's dayside and nightside temperatures as $T_{\text{day}} = 3004 \pm 64$ K and $T_{\text{night}} = 2583 \pm 63$ K, respectively, when assuming zero dayside geometric albedo. Among the population of hot/ultrahot Jupiters, TOI-1431b has the second highest measured nightside temperature and day-night heat recirculation efficiency.

It is also an excellent candidate for future follow-up observations with JWST and ARIEL to measure its transmission and emission spectra as well as to obtain a high-precision global temperature and cloud map of this planet's atmosphere. Furthermore, the discovery and characterization of planets orbiting Am stars, for which few planets have been found, provide good opportunities to probe the tidal interactions between Jovian planets and hot host stars and the potential mechanisms responsible for the creation and evolution of Am stars.

Based on observations made with the Hertzsprung SONG telescope operated on the Spanish Observatorio del Teide on the island of Tenerife by the Aarhus and Copenhagen Universities and by the Instituto de Astrofísica de Canarias.

Funding for the Stellar Astrophysics Centre is provided by The Danish National Research Foundation (grant agreement no.: DNRF106).

Based on observations made with the Nordic Optical Telescope, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos,
La Palma, Spain, of the Instituto de Astrofísica de Canarias under programs 59-210, 59-503, and 61-804.

Based on observations collected with the SOPHIE spectrograph on the 1.93 m telescope at the Observatoire de Haute-Provence (CNRS), France.

This article is based on observations made with the MuSCAT2 instrument, developed by ABC, at Telescopio Carlos Sánchez operated on the island of Tenerife by the IAC in the Spanish Observatorio del Teide.

This work makes use of observations from the LCOGT network. LCOGT telescope time was granted by NOIRLab through the Mid-Scale Innovations Program (MSIP). MSIP is funded by NSF.

B.A. is supported by Australian Research Council Discovery Grant DP180100972. The SONG telescope network is partially supported by the Australian Research Council Linkage, Infrastructure, Equipment and Facilities grant LE190100036.

V.A. was supported by a research grant (00028173) from VILLUM FONDEN.

I.W. is supported by a Heising-Simons 51 Pegasi b postdoctoral fellowship.

Funding for the TESS mission is provided by NASA's Science Mission directorate. We acknowledge the use of public TESS Alert data from pipelines at the TESS Science Office and at the TESS Science Processing Operations Center. This research has made use of the Exoplanet Follow-up Observation Program website, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. Resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center for the production of the SPOC data products. This paper includes data collected by the TESS mission, which are publicly available from the Mikulski Archive for Space Telescopes (MAST).

T.D. acknowledges support from MIT's Kavli Institute as a Kavli postdoctoral fellow.

M.Y. and H.V.S. acknowledge the support by the research fund of Ankara University (BAP) through the project 18A0759001.

This work was supported by Fundação para a Ciência e a Tecnologia (FCT) and Fundo Europeu de Desenvolvimento

Regional (FEDER) via COMPETE2020 through the research grants UIDB/04434/2020, UIDP/04434/2020, PTDC/FIS-AST/32113/2017 and POCI-01-0145-FEDER-032113, PTDC/FIS-AST/28953/2017 and POCI-01-0145-FEDER-028953.

O.D.S.D. is supported in the form of work contract (DL 57/2016/CP1364/CT0004) funded by FCT.

This work was supported by FCT—Fundação para a Ciência e a Tecnologia through national funds and by FEDER through COMPETE2020—Programa Operacional Competitividade e Internacionalização by these grants: UID/FIS/04434/2019; UIDB/04434/2020; UIDP/04434/2020; PTDC/FIS-AST/32113/2017 and POCI-01-0145-FEDER-032113; PTDC/FIS-AST/28953/2017 and POCI-01-0145-FEDER-028953; PTDC/FIS-AST/28987/2017 and POCI-01-0145-FEDER-028987.

S.H. acknowledges CNES funding through the grant 837319.

G.H. and I.B. received funding from the French Programme National de Planétologie (PNP) of CNRS (INSU).

This work is partly supported by JSPS KAKENHI Grant Nos. JP18H01265 and JP18H05439, and JST PRESTO grant No. JPMJPR1775, and a University Research Support Grant from the National Astronomical Observatory of Japan (NAOJ).

M.N.G. acknowledges support from MIT's Kavli Institute as a Juan Carlos Torres Fellow.

T.D.K. acknowledges support from the 51 Pegasi b Fellowship in Planetary Astronomy sponsored by the Heising-Simons Foundation.

D.R.L. is supported by an appointment to the NASA Postdoctoral Program at NASA Goddard Space Flight Center, administered by Universities Space Research Association under contract with NASA.

**Facilities:** TESS, SONG, SOPHIE, NRES ELP, NRES TLV, FIES, EXPRES, LCOGT, MuSCAT2, MASCARA.

**Software:** AstroImageJ (Collins et al. 2017), TAPIR (Jensen 2013), Allesfitter (Günther & Daylan 2019, 2021), EXOFASTv2 (Eastman et al. 2013, 2019; Eastman 2017), Astropy (Astropy Collaboration et al. 2018), Matplotlib (Hunter 2007), Wotan (Hippke et al. 2019).

## Appendix
<figure><img src="image_22.png" /><figcaption>**Figure 14.** A transit observed on 2020 May 16 by MuSCAT2 using simultaneous multicolor photometry in g', r', i', and z_s bands. The top panel shows the distribution of residuals to the fit of the light curve in each band. The second panel from the top shows the fit to the transit plus systematics. The third panel shows the fit to the transit with the systematics removed. The three blue vertical bars in the second, third, and fourth panels represent the predicted ingress, mid-transit, and egress times along with their $1\sigma$ uncertainties based on the ephemeris used at the time. The horizontal dotted blue lines in the third panel show the predicted transit depth. The bottom panel shows the resulting residuals to the fit of the light curves.</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>**Figure 15.** Phase-folded and detrended light curve of TOI-1431 from the simultaneous hybrid cubic spline detrending model using the SAP data. The colored points represent the individual phased light-curve segments across Sectors 15 and 16 that have been binned at a cadence of 15 minutes. The red solid lines are 20 light-curve models drawn from the posteriors of the nested sampling analysis in Allesfitter.</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>**Figure 16.** TESS Sector 15 photometry of TOI-1431 from the SAP (blue points) and single-scale MAP PDC (red points) light curves. The SAP light curve shows “whisker”-like flux dips that are caused by brief spacecraft pointing excursions, some of which overlap with the transit events. The correction applied to the single-scale MAP PDC light curve from the use of the cotrending basis vectors has removed these features.</figcaption></figure>

<figure><img src="image_25.png" /><figcaption>**Figure 17.** Radius evolution of TOI-1431b for various assumptions about the presence of irradiation and deep deposited heating. The curve labeled Thorngren & Fortney (2018) includes irradiation and applies deposited heating in the deep interior with a strength that varies with equilibrium temperature, as in Equation (34) of Thorngren & Fortney (2018). The models labeled $\gamma = 5 \times 10^{-4}$ and $\gamma = 1 \times 10^{-4}$ respectively assume that 0.05% and 0.01% of the incident stellar power is converted to heat deposited deep in the planetary interior and include irradiation. The model labeled “irradiation” considers only irradiation with no deposited heating, and the model labeled “no heating” does not include irradiation or deposited heating. Evolution models are conducted with MESA (Paxton et al. 2011, 2013, 2015, 2018, 2019) using the model setup described in Komacek et al. (2020) and stellar evolution tracks from Choi et al. (2016) and Dotter (2016). The present-day radius of TOI-1431b can be fit by planetary evolution models including only irradiation slowing cooling, but a small amount of deep deposited heating corresponding to <0.05% of the incident stellar power converted to heat is allowed.</figcaption></figure>
<figure><img src="image_26.png" /><figcaption>Figure 18. The phase-folded TESS SAP light curve for TOI-1431b, zoomed in to show the phase curve and secondary eclipse, similar to Figure 11.</figcaption></figure>

<table><caption>**Table 3**
A Summary of the Ground-based Transit Follow-up Observations Taken of TOI-1431</caption><thead><tr><th>Telescope</th><th>Camera</th><th>Filter</th><th>UT Date</th><th>Coverage</th><th>Precision (ppt)</th><th>Joint Fit</th><th>Comments</th></tr></thead><tbody><tr><td>CDK14</td><td>STXL-6303E</td><td>z'</td><td>2019-12-24</td><td>Full</td><td>2.62</td><td>Yes</td><td>Minimal systematics</td></tr><tr><td>CDK14</td><td>STXL-6303E</td><td>g'</td><td>2019-12-24</td><td>Full</td><td>2.32</td><td>Yes</td><td>Minimal systematics</td></tr><tr><td>TCS</td><td>MuSCAT2</td><td>g'</td><td>2020-05-16</td><td>Full</td><td>0.93</td><td>Yes</td><td>Minimal systematics</td></tr><tr><td>TCS</td><td>MuSCAT2</td><td>r'</td><td>2020-05-16</td><td>Full</td><td>0.97</td><td>Yes</td><td>Minimal systematics</td></tr><tr><td>TCS</td><td>MuSCAT2</td><td>i'</td><td>2020-05-16</td><td>Full</td><td>0.94</td><td>Yes</td><td>Minimal systematics</td></tr><tr><td>TCS</td><td>MuSCAT2</td><td>z<sub>s</sub></td><td>2020-05-16</td><td>Full</td><td>1.01</td><td>Yes</td><td>Minimal systematics</td></tr><tr><td>TCS</td><td>MuSCAT2</td><td>g'</td><td>2020-05-23</td><td>Full</td><td>1.23</td><td>No</td><td>Limited pre-ingress data and systematics</td></tr><tr><td>TCS</td><td>MuSCAT2</td><td>r'</td><td>2020-05-23</td><td>Full</td><td>1.07</td><td>No</td><td>Limited pre-ingress data and systematics</td></tr><tr><td>TCS</td><td>MuSCAT2</td><td>i'</td><td>2020-05-23</td><td>Full</td><td>1.08</td><td>No</td><td>Limited pre-ingress data and systematics</td></tr><tr><td>TCS</td><td>MuSCAT2</td><td>z<sub>s</sub></td><td>2020-05-23</td><td>Full</td><td>1.09</td><td>No</td><td>Limited pre-ingress data and systematics</td></tr><tr><td>AUKR</td><td>ALTA U47</td><td>z'</td><td>2020-06-16</td><td>Full</td><td>3.00</td><td>Yes</td><td>Some systematics</td></tr><tr><td>LCOGT-McD</td><td>Sinistro</td><td>PANSTARRS Y</td><td>2020-07-13</td><td>Ingress</td><td>0.48</td><td>No</td><td>Ingress only</td></tr><tr><td>SCT</td><td>ST7XME</td><td>TESS band</td><td>2020-08-08</td><td>Full</td><td>1.49</td><td>Yes</td><td>Minimal systematics</td></tr><tr><td>ULMT</td><td>STX 16803</td><td>i'</td><td>2020-09-20</td><td>Full</td><td>1.51</td><td>Yes</td><td>Some systematics</td></tr><tr><td>ULMT</td><td>STX 16803</td><td>i'</td><td>2020-09-28</td><td>Full</td><td>3.41</td><td>No</td><td>Significant systematics</td></tr><tr><td>LCOGT-McD</td><td>Sinistro</td><td>PANSTARRS Y</td><td>2020-10-14</td><td>Full</td><td>1.18</td><td>Yes</td><td>Minimal systematics</td></tr></tbody></table>
<table><caption>**Table 4**
Radial Velocity Measurements</caption><thead><tr><th>Date<sup>a</sup><br>BJD – 2,400,000</th><th>RV<sup>b</sup><br>(m s<sup>-1</sup>)</th><th>σ<sub>RV</sub><br>(m s<sup>-1</sup>)</th><th>Instrument</th></tr></thead><tbody><tr><td>58,805.564199</td><td>-25,696.6</td><td>14.5</td><td>ELP</td></tr><tr><td>58,806.565002</td><td>-25,761.0</td><td>15.3</td><td>ELP</td></tr><tr><td>58,807.566340</td><td>-25,346.6</td><td>13.3</td><td>ELP</td></tr><tr><td>58,808.570467</td><td>-25,830.7</td><td>29.9</td><td>ELP</td></tr><tr><td>58,814.663425</td><td>-25,557.9</td><td>26.2</td><td>ELP</td></tr><tr><td>58,815.307593</td><td>9.2</td><td>14.2</td><td>FIES</td></tr><tr><td>58,818.618086</td><td>-25,457.5</td><td>18.0</td><td>ELP</td></tr><tr><td>58,836.224800</td><td>-25,039.0</td><td>7.0</td><td>SOPHIE</td></tr><tr><td>58,838.252200</td><td>-25,441.0</td><td>7.0</td><td>SOPHIE</td></tr><tr><td>58,840.399100</td><td>-25,532.0</td><td>45.0</td><td>SOPHIE</td></tr><tr><td>58,841.235500</td><td>-25,222.0</td><td>7.0</td><td>SOPHIE</td></tr><tr><td>58,841.352100</td><td>-25,152.0</td><td>6.0</td><td>SOPHIE</td></tr><tr><td>58,859.232100</td><td>-25,528.0</td><td>7.0</td><td>SOPHIE</td></tr><tr><td>58,860.267700</td><td>-24,987.0</td><td>6.0</td><td>SOPHIE</td></tr><tr><td>58,861.259600</td><td>-25,374.0</td><td>6.0</td><td>SOPHIE</td></tr><tr><td>58,914.768952</td><td>-25,465.4</td><td>53.1</td><td>SONG</td></tr><tr><td>58,919.763375</td><td>-25,434.9</td><td>81.7</td><td>SONG</td></tr><tr><td>58,921.764647</td><td>-24,969.6</td><td>74.1</td><td>SONG</td></tr><tr><td>58,942.562272</td><td>-25,283.8</td><td>19.2</td><td>TLV</td></tr><tr><td>58,943.567282</td><td>-25,788.4</td><td>15.4</td><td>TLV</td></tr><tr><td>58,948.564452</td><td>-25,514.9</td><td>16.4</td><td>TLV</td></tr><tr><td>58,952.584719</td><td>-25,528.9</td><td>16.3</td><td>TLV</td></tr><tr><td>58,955.532830</td><td>-25,299.7</td><td>28.7</td><td>TLV</td></tr><tr><td>58,956.668739</td><td>-25,303.1</td><td>139.4</td><td>SONG</td></tr><tr><td>58,961.645589</td><td>-25,067.3</td><td>56.1</td><td>SONG</td></tr><tr><td>58,965.721756</td><td>-25,219.0</td><td>44.1</td><td>SONG</td></tr><tr><td>58,980.495174</td><td>-25,900.9</td><td>19.9</td><td>TLV</td></tr><tr><td>58,982.464359</td><td>-25,256.1</td><td>27.5</td><td>TLV</td></tr><tr><td>58,983.637401</td><td>-25,477.9</td><td>64.2</td><td>SONG</td></tr><tr><td>58,987.442635</td><td>-25,265.9</td><td>22.1</td><td>TLV</td></tr><tr><td>58,988.470752</td><td>-25,659.4</td><td>27.9</td><td>TLV</td></tr><tr><td>58,988.619270</td><td>-25,374.0</td><td>26.6</td><td>SONG</td></tr><tr><td>58,989.466532</td><td>-25,588.9</td><td>42.2</td><td>TLV</td></tr><tr><td>58,990.449283</td><td>-25,278.6</td><td>20.1</td><td>TLV</td></tr><tr><td>58,990.567802</td><td>-24.1</td><td>7.9</td><td>FIES</td></tr><tr><td>58,991.456242</td><td>-25,896.5</td><td>26.1</td><td>TLV</td></tr><tr><td>58,991.562807</td><td>-581.4</td><td>15.6</td><td>FIES</td></tr><tr><td>58,991.594498</td><td>-565.9</td><td>12.9</td><td>FIES</td></tr><tr><td>58,991.642527</td><td>-25,435.4</td><td>18.0</td><td>SONG</td></tr><tr><td>58,991.698593</td><td>-565.4</td><td>11.0</td><td>FIES</td></tr><tr><td>58,992.606914</td><td>-69.5</td><td>12.5</td><td>FIES</td></tr><tr><td>58,992.615274</td><td>-75.2</td><td>18.9</td><td>FIES</td></tr><tr><td>58,993.658735</td><td>-318.3</td><td>11.6</td><td>FIES</td></tr><tr><td>58,993.662679</td><td>-339.5</td><td>15.6</td><td>FIES</td></tr><tr><td>58,993.666687</td><td>-314.2</td><td>14.1</td><td>FIES</td></tr><tr><td>58,993.670878</td><td>-329.4</td><td>18.3</td><td>FIES</td></tr><tr><td>58,993.674862</td><td>-315.8</td><td>11.2</td><td>FIES</td></tr><tr><td>58,993.678937</td><td>-333.8</td><td>13.4</td><td>FIES</td></tr><tr><td>58,993.686045</td><td>-334.8</td><td>9.8</td><td>FIES</td></tr><tr><td>58,993.690092</td><td>-329.0</td><td>13.9</td><td>FIES</td></tr><tr><td>58,993.694137</td><td>-337.2</td><td>13.8</td><td>FIES</td></tr><tr><td>58,993.698146</td><td>-345.5</td><td>9.7</td><td>FIES</td></tr><tr><td>58,994.571980</td><td>-500.4</td><td>9.3</td><td>FIES</td></tr><tr><td>58,994.621049</td><td>-471.6</td><td>8.0</td><td>FIES</td></tr><tr><td>58,994.704411</td><td>-429.4</td><td>11.7</td><td>FIES</td></tr><tr><td>58,995.443468</td><td>-25,324.0</td><td>28.7</td><td>TLV</td></tr><tr><td>58,995.570851</td><td>15.5</td><td>11.5</td><td>FIES</td></tr><tr><td>58,995.620876</td><td>13.7</td><td>7.2</td><td>FIES</td></tr><tr><td>58,995.703476</td><td>0.0</td><td>11.9</td><td>FIES</td></tr><tr><td>58,996.502720</td><td>-25,732.7</td><td>26.1</td><td>TLV</td></tr><tr><td>58,997.422153</td><td>-25,702.1</td><td>38.1</td><td>TLV</td></tr><tr><td>58,998.422132</td><td>-25,287.5</td><td>23.5</td><td>TLV</td></tr><tr><td>58,999.424288</td><td>-25,873.9</td><td>19.0</td><td>TLV</td></tr></tbody></table>

<table><caption>**Table 4**
(Continued)</caption><thead><tr><th>Date<sup>a</sup><br>BJD – 2,400,000</th><th>RV<sup>b</sup><br>(m s<sup>-1</sup>)</th><th>σ<sub>RV</sub><br>(m s<sup>-1</sup>)</th><th>Instrument</th></tr></thead><tbody><tr><td>58,999.624293</td><td>-25,458.0</td><td>76.4</td><td>SONG</td></tr><tr><td>59,001.471298</td><td>-25,606.2</td><td>20.1</td><td>TLV</td></tr><tr><td>59,001.641856</td><td>-25,224.1</td><td>28.8</td><td>SONG</td></tr><tr><td>59,002.428683</td><td>-25,802.8</td><td>21.2</td><td>TLV</td></tr><tr><td>59,002.787631</td><td>-55.3</td><td>2.5</td><td>EXPRES</td></tr><tr><td>59,003.954076</td><td>167.7</td><td>1.8</td><td>EXPRES</td></tr><tr><td>59,004.452546</td><td>-25,721.2</td><td>16.9</td><td>TLV</td></tr><tr><td>59,004.602433</td><td>-25,375.3</td><td>24.6</td><td>SONG</td></tr><tr><td>59,004.953876</td><td>-279.8</td><td>2.0</td><td>EXPRES</td></tr><tr><td>59,005.450710</td><td>-25,560.8</td><td>18.8</td><td>TLV</td></tr><tr><td>59,006.382260</td><td>-25,320.8</td><td>19.6</td><td>TLV</td></tr><tr><td>59,007.421885</td><td>-25,873.3</td><td>17.2</td><td>TLV</td></tr><tr><td>59,007.604415</td><td>-25,444.7</td><td>39.9</td><td>SONG</td></tr><tr><td>59,008.456417</td><td>-25,235.3</td><td>34.6</td><td>TLV</td></tr><tr><td>59,009.487869</td><td>-25,570.9</td><td>18.7</td><td>TLV</td></tr><tr><td>59,009.566599</td><td>-25,187.8</td><td>24.2</td><td>SONG</td></tr><tr><td>59,010.502769</td><td>-25,751.8</td><td>19.5</td><td>TLV</td></tr><tr><td>59,010.624528</td><td>-25,295.0</td><td>24.9</td><td>SONG</td></tr><tr><td>59,012.426456</td><td>-25,777.6</td><td>16.5</td><td>TLV</td></tr><tr><td>59,012.640398</td><td>-25,419.6</td><td>29.3</td><td>SONG</td></tr><tr><td>59,012.955625</td><td>-271.6</td><td>1.9</td><td>EXPRES</td></tr><tr><td>59,018.520808</td><td>-447.3</td><td>17.6</td><td>FIES</td></tr><tr><td>59,018.597262</td><td>-406.4</td><td>15.7</td><td>FIES</td></tr><tr><td>59,018.693433</td><td>-330.5</td><td>8.5</td><td>FIES</td></tr><tr><td>59,019.513133</td><td>34.7</td><td>8.4</td><td>FIES</td></tr><tr><td>59,019.594250</td><td>2.8</td><td>7.5</td><td>FIES</td></tr><tr><td>59,019.694408</td><td>-25.8</td><td>9.3</td><td>FIES</td></tr><tr><td>59,020.514030</td><td>-516.7</td><td>10.5</td><td>FIES</td></tr><tr><td>59,020.590422</td><td>-546.9</td><td>9.0</td><td>FIES</td></tr><tr><td>59,020.688977</td><td>-552.9</td><td>12.6</td><td>FIES</td></tr><tr><td>59,023.610494</td><td>-25,418.5</td><td>21.3</td><td>SONG</td></tr><tr><td>59,027.569361</td><td>-24,874.0</td><td>24.4</td><td>SONG</td></tr><tr><td>59,029.510059</td><td>-25,065.6</td><td>20.3</td><td>SONG</td></tr><tr><td>59,033.957319</td><td>-287.7</td><td>2.2</td><td>EXPRES</td></tr><tr><td>59,059.915400</td><td>-36.2</td><td>2.1</td><td>EXPRES</td></tr><tr><td>59,059.919605</td><td>-44.0</td><td>2.0</td><td>EXPRES</td></tr><tr><td>59,059.923762</td><td>-38.9</td><td>2.0</td><td>EXPRES</td></tr><tr><td>59,059.927997</td><td>-45.9</td><td>2.2</td><td>EXPRES</td></tr><tr><td>59,059.932407</td><td>-51.7</td><td>2.2</td><td>EXPRES</td></tr><tr><td>59,059.936606</td><td>-50.3</td><td>2.3</td><td>EXPRES</td></tr><tr><td>59,059.940835</td><td>-50.7</td><td>2.2</td><td>EXPRES</td></tr><tr><td>59,059.945055</td><td>-58.2</td><td>2.1</td><td>EXPRES</td></tr><tr><td>59,059.949247</td><td>-59.5</td><td>1.9</td><td>EXPRES</td></tr><tr><td>59,059.953490</td><td>-62.4</td><td>2.3</td><td>EXPRES</td></tr></tbody></table>

**Notes.**

<sup>a</sup> The dates for each observation are reported as BJD at the UTC time at the midpoint of the exposure.

<sup>b</sup> FIES and EXPRES radial velocities are given at an arbitrary zero-point.

(This table is available in machine-readable form.)

**ORCID iDs**

Brett C. Addison https://orcid.org/0000-0003-3216-0626

Emil Knudstrup https://orcid.org/0000-0001-7880-594X

Ian Wong https://orcid.org/0000-0001-9665-8429

Ignas Snellen https://orcid.org/0000-0003-1624-3667

Simon Albrecht https://orcid.org/0000-0003-1762-8235

Aaron Bello-Arufe https://orcid.org/0000-0003-3355-1223

Jose-Manuel Almenara https://orcid.org/0000-0003-3208-9815
Sergio Hoyer https://orcid.org/0000-0003-3477-2466

N. C. Santos https://orcid.org/0000-0003-4422-2919

Grzegorz Nowak https://orcid.org/0000-0002-7031-7754

Rafael Luque https://orcid.org/0000-0002-4671-2957

Monika Stangret https://orcid.org/0000-0002-1812-8024

Enric Palle https://orcid.org/0000-0003-0987-1593

René Tronsgaard https://orcid.org/0000-0003-1001-0707

Victoria Antoci https://orcid.org/0000-0002-0865-3650

Lars A. Buchhave https://orcid.org/0000-0003-1605-5666

Maximilian N. Günther https://orcid.org/0000-0002-3164-9086

Tansu Daylan https://orcid.org/0000-0002-6939-9211

Felipe Murgas https://orcid.org/0000-0001-9087-1245

Hannu Parviainen https://orcid.org/0000-0001-5519-1391

Emma Esparza-Borges https://orcid.org/0000-0002-2341-3233

Nicolas Crouzet https://orcid.org/0000-0001-7866-8738

Norio Narita https://orcid.org/0000-0001-8511-2981

Akihiko Fukui https://orcid.org/0000-0002-4909-5763

Kiyoe Kawauchi https://orcid.org/0000-0003-1205-5108

Noriharu Watanabe https://orcid.org/0000-0002-7522-8195

Markus Rabus https://orcid.org/0000-0003-2935-7196

Marshall C. Johnson https://orcid.org/0000-0002-5099-8185

Gilles P. P. L. Otten https://orcid.org/0000-0002-6717-1977

Geert Jan Talens https://orcid.org/0000-0003-4787-2335

Samuel H. C. Cabot https://orcid.org/0000-0001-9749-6150

Debra A. Fischer https://orcid.org/0000-0003-2221-0861

Frank Grundahl https://orcid.org/0000-0002-8736-1639

Mads Fredslund Andersen https://orcid.org/0000-0002-9194-8520

Pere Pallé https://orcid.org/0000-0003-3803-4823

Avi Shporer https://orcid.org/0000-0002-1836-3120

David R. Ciardi https://orcid.org/0000-0002-5741-3047

Jake T. Clark https://orcid.org/0000-0003-3964-4658

Robert A. Wittenmyer https://orcid.org/0000-0001-9957-9304

Jonathan Horner https://orcid.org/0000-0002-1160-7970

Karen A. Collins https://orcid.org/0000-0001-6588-9574

Eric L. N. Jensen https://orcid.org/0000-0002-4625-7333

John F. Kielkopf https://orcid.org/0000-0003-0497-2651

Richard P. Schwarz https://orcid.org/0000-0001-8227-1020

Mesut Yilmaz https://orcid.org/0000-0002-3276-0704

Hakan Volkan Senavci https://orcid.org/0000-0002-8961-277X

Daniel Harbeck https://orcid.org/0000-0002-8590-007X

Thaddeus D. Komacek https://orcid.org/0000-0002-9258-5311

Songhu Wang https://orcid.org/0000-0002-7846-6981

Jason D. Eastman https://orcid.org/0000-0003-3773-5142

Keivan G. Stassun https://orcid.org/0000-0002-3481-9052

David W. Latham https://orcid.org/0000-0001-9911-7388

Roland Vanderspek https://orcid.org/0000-0001-6763-6562

Sara Seager https://orcid.org/0000-0002-6892-6948

Joshua N. Winn https://orcid.org/0000-0002-4265-047X

Jon M. Jenkins https://orcid.org/0000-0002-4715-9460

Dana R. Louie https://orcid.org/0000-0002-2457-272X

Luke G. Bouma https://orcid.org/0000-0002-0514-5538

Joseph D. Twicken https://orcid.org/0000-0002-6778-7552

Alan M. Levine https://orcid.org/0000-0001-8172-0453

Brian McLean https://orcid.org/0000-0002-8058-643X

References

Abt, H. A., & Morrell, N. I. 1995, ApJS, 99, 135

Addison, B. C., Tinney, C. G., Wright, D. J., et al. 2013, ApJL, 774, L9

Addison, B. C., Wright, D. J., Nicholson, B. A., et al. 2021, MNRAS, 502, 3704

Albrecht, S., Winn, J. N., Butler, R. P., et al. 2012, ApJ, 744, 189

Alibert, Y., Mordasini, C., Benz, W., & Winisdoerffer, C. 2005, A&A, 434, 343

Andersen, M. F., Grundahl, F., Christensen-Dalsgaard, J., et al. 2014, RMxAA, 45, 83

Anderson, D. R., Temple, L. Y., Nielsen, L. D., et al. 2018, arXiv:1809.04897

Antoci, V., Handler, G., Grundahl, F., et al. 2013, MNRAS, 435, 1563

Arcangeli, J., Désert, J.-M., Line, M. R., et al. 2018, ApJL, 855, L30

Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481

Astropy Collaboration, Price-Whelan, A. M., Sipőcz, B. M., et al. 2018, AJ, 156, 123

Batygin, K., Bodenheimer, P. H., & Laughlin, G. P. 2016, ApJ, 829, 114

Beaugé, C., & Nesvorný, D. 2012, ApJ, 751, 119

Bell, T. J., & Cowan, N. B. 2018, ApJL, 857, L20

Bell, T. J., Dang, L., Cowan, N. B., et al. 2021, MNRAS, 504, 3316

Bergmann, C., Jones, M. I., Zhao, J., et al. 2021, PASA, 38, 19

Blackman, R. T., Fischer, D. A., Jurgenson, C. A., et al. 2020, AJ, 159, 238

Blanco-Cuaresma, S. 2019, MNRAS, 486, 2075

Blanco-Cuaresma, S., Soubiran, C., Heiter, U., & Jofré, P. 2014, A&A, 569, A111

Böhm-Vitense, E. 2006, PASP, 118, 419

Böhm-Vitense, E. 2007, AJ, 133, 1903

Bolton, C. T. 1971, A&A, 14, 233

Borucki, W. J., Koch, D., Basri, G., et al. 2010, Sci, 327, 977

Bouchy, F., Hébrard, G., Udry, S., et al. 2009, A&A, 505, 853

Bowler, B. P., Johnson, J. A., Marcy, G. W., et al. 2010, ApJ, 709, 396

Brahm, R., Jordán, A., & Espinoza, N. 2017, PASP, 129, 034002

Brahm, R., Nielsen, L. D., Wittenmyer, R. A., et al. 2020, AJ, 160, 235

Brewer, J. M., Fischer, D. A., Blackman, R. T., et al. 2020, AJ, 160, 67

Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031

Bryson, S., Kunimoto, M., Kopparapu, R. K., et al. 2021, AJ, 161, 36

Buchhave, L. A., Bizzarro, M., Latham, D. W., et al. 2014, Natur, 509, 593

Buchhave, L. A., Bakos, G. Á., Hartman, J. D., et al. 2010, ApJ, 720, 1118

Buchhave, L. A., Latham, D. W., Johansen, A., et al. 2012, Natur, 486, 375

Butler, R. P., Marcy, G. W., Fischer, D. A., et al. 1999, ApJ, 526, 916

Butler, R. P., Marcy, G. W., Vogt, S. S., & Apps, K. 1998, PASP, 110, 1389

Cañas, C. I., Wang, S., Mahadevan, S., et al. 2019, ApJL, 870, L17

Carquillat, J. M., & Prieur, J. L. 2007, MNRAS, 380, 1064

Castelli, F., & Kurucz, R. L. 2003, in IAU Symp. 210, Modelling of Stellar Atmospheres, ed. N. Piskunov, W. W. Weiss, & D. F. Gray (Cambridge: Cambridge Univ. Press), A20

Charbonneau, P. 1993, in ASP Conf. Ser. 44, IAU Coll. 138: Peculiar versus Normal Phenomena in A-type and Related Stars, ed. M. M. Dworetsky, F. Castelli, & R. Faraggiana (San Francisco, CA: ASP), 474

Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, ApJ, 686, 580

Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102

Claret, A. 2017, A&A, 600, A30

Cochran, W. D., Hatzes, A. P., Butler, R. P., & Marcy, G. W. 1997, ApJ, 483, 457

Collier Cameron, A., Bruce, V. A., Miller, G. R. M., Triaud, A. H. M. J., & Queloz, D. 2010a, MNRAS, 403, 151

Collier Cameron, A., Guenther, E., Smalley, B., et al. 2010b, MNRAS, 407, 507

Collins, K. A., Kielkopf, J. F., Stassun, K. G., & Hessman, F. V. 2017, AJ, 153, 77

Conti, P. S. 1970, PASP, 82, 781

Courcol, B., Bouchy, F., Pepe, F., et al. 2015, A&A, 581, A38

Cowan, N. B., & Agol, E. 2011, ApJ, 729, 54

Cumming, A., Marcy, G. W., & Butler, R. P. 1999, ApJ, 526, 890

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, yCat, 2246, 0

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2013, Explanatory Supplement to the AllWISE Data Release Products

Davis, A. B., Wang, S., Jones, M., et al. 2020, AJ, 160, 229

Daylan, T., Günther, M. N., Mikal-Evans, T., et al. 2021, AJ, 161, 131

Djupvik, A. A., & Andersen, J. 2010, ASSP, 14, 211

Dorn, C., Harrison, J. H. D., Bonsor, A., & Hands, T. O. 2019, MNRAS, 484, 712

Dorval, P., Talens, G. J. J., Otten, G. P. P. L., et al. 2020, A&A, 635, A60

Dotter, A. 2016, ApJS, 222, 8
Eastman, J. 2017, EXOFASTv2: Generalized publication-quality exoplanet modeling code, Astrophysics Source Code Library, http://ascl.net/ 1710.003

Eastman, J., Gaudí, B. S., & Agol, E. 2013, PASP, 125, 83

Eastman, J. D., Rodriguez, J. E., Agol, E., et al. 2019, arXiv:1907.09480

Endl, M., Cochran, W. D., Wittenmyer, R. A., & Boss, A. P. 2008, ApJ, 673, 1165

Figueras, F., Torra, J., & Jordi, C. 1991, A&AS, 87, 319

Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661

Frandsen, S., & Lindberg, B. 1999, in Astrophysics with the NOT, ed. H. Karttunen & V. Piiriola (Piikkio: University of Turku, Tuorla Observatory), 71

Fredslund Andersen, M., Handberg, R., Weiss, E., et al. 2019, PASP, 131, 045003

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018a, A&A, 616, A1

Gaudí, B. S., Stassun, K. G., Collins, K. A., et al. 2017, Natur, 546, 514

Gilbert, E. A., Barclay, T., Schlieder, J. E., et al. 2020, AJ, 160, 116

Ginzburg, S., & Sari, R. 2015, ApJ, 803, 111

Gray, R. O., & Corbally, C. J. 1994, AJ, 107, 742

Grundahl, F., Fredslund Andersen, M., Christensen-Dalsgaard, J., et al. 2017, ApJ, 836, 142

Guenther, E. W., Barragán, O., Dai, F., et al. 2017, A&A, 608, A93

Guerrero, N. M., Seager, S., Huang, C. X., et al. 2021, ApJS, 254, 39

Günther, M. N., & Daylan, T. 2019, Allesfitter: Flexible Star and Exoplanet Inference From Photometry and Radial Velocity, Astrophysics Source Code Library, http://ascl.net/1903.003

Günther, M. N., & Daylan, T. 2021, ApJS, 254, 13

Günther, M. N., Pozuelos, F. J., Dittmann, J. A., et al. 2019, NatAs, 3, 1099

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

Hasegawa, Y., Yu, T. Y. M., & Hansen, B. M. S. 2019, A&A, 629, L1

Heiter, U., Lind, K., Asplund, M., et al. 2015, PhyS, 90, 054010

Hellier, C., Anderson, D. R., Barkaoui, K., et al. 2019, MNRAS, 490, 1479

Henry, G. W., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2000, ApJL, 529, L41

Hippke, M., David, T. J., Mulders, G. D., & Heller, R. 2019, AJ, 158, 143

Hirano, T., Narita, N., Sato, B., et al. 2011, PASJ, 63, L57

Hjorth, M., Albrecht, S., Hirano, T., et al. 2021, PNAS, 118, 2017418118

Hoejmakers, H. J., Seidel, J. V., Pino, L., et al. 2020, A&A, 641, A123

Høg, E., Fabricius, C., Makarov, V. V., et al. 2009, A&A, 355, L27

Horner, J., Kane, S. R., Marshall, J. P., et al. 2020, PASP, 132, 102001

Hsu, D. C., Ford, E. B., Ragozzine, D., & Ashby, K. 2019, AJ, 158, 109

Huang, C. X., Burt, J., Vanderburg, A., et al. 2018, ApJL, 868, L39

Hunter, J. D. 2007, CSE, 9, 90

Husser, T. O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6

Ida, S., & Lin, D. N. C. 2005, ApJ, 626, 1045

Jenkins, J. M., Twicken, J. D., McCauliff, S., et al. 2016, Proc. SPIE, 9913, 99133E

Jensen, E. 2013, Tapir: A web interface for transit/eclipse observability, Astrophysics Source Code Library, http://ascl.net/1306.007

Johnson, J. A., Marcy, G. W., Fischer, D. A., et al. 2006, ApJ, 652, 1724

Johnson, M. C., Cochran, W. D., Albrecht, S., et al. 2014, ApJ, 790, 30

Jones, H. R. A., Butler, R. P., Tinney, C. G., et al. 2006, MNRAS, 369, 249

Jones, M. I., Brahm, R., Espinoza, N., et al. 2019, A&A, 625, A16

Jones, M. I., Jenkins, J. S., Brahm, R., et al. 2016, A&A, 590, A38

Jordán, A., Brahm, R., Espinoza, N., et al. 2020, AJ, 159, 145

Jurgenson, C., Fischer, D., McCracken, T., et al. 2016, Proc. SPIE, 9908, 99086T

Kempton, E. M. R., Bean, J. L., Louie, D. R., et al. 2018, PASP, 130, 114401

Komacek, T. D., Thorngren, D. P., Lopez, E. D., & Ginzburg, S. 2020, ApJ, 893, 36

Kovács, G., Zucker, S., & Mazeh, T. 2002, A&A, 391, 369

Kozai, Y. 1962, AJ, 67, 591

Kurucz, R. L. 1993, SYNTHE spectrum synthesis programs and line data (None) (Cambridge, MA: Smithsonian Astrophysical Observatory)

Levine, S. E., Bida, T. A., Chylek, T., et al. 2012, Proc. SPIE, 8444, 844419

Li, D., Mustill, A. J., & Davies, M. B. 2019a, MNRAS, 488, 1366

Li, J., Tenenbaum, P., Twicken, J. D., et al. 2019b, PASP, 131, 024506

Lidov, M. L. 1962, P&SS, 9, 719

Lin, D. N. C., Bodenheimer, P., & Richardson, D. C. 1996, Natur, 380, 606

Lothringer, J. D., Barman, T., & Koskinen, T. 2018, ApJ, 866, 27

Louden, E. M., Winn, J. N., Petigura, E. A., et al. 2021, AJ, 161, 68

Luhn, J. K., Schliestien, F. A., Wright, J. T., et al. 2009, AJ, 157, 149

Lund, M. B., Rodriguez, J. E., Zhou, G., et al. 2017, AJ, 154, 194

Martínez, R. R., Gaudí, B. S., Rodriguez, J. E., et al. 2020, AJ, 160, 111

Mayor, M., & Queloz, D. 1995, Natur, 378, 355

McCully, C., Volgenau, N. H., Harbeck, D.-R., et al. 2018, Proc. SPIE, 10707, 107070K

McLaughlin, D. B. 1924, ApJ, 60, 22

Mészáros, S., Allende Prieto, C., Edvardsson, B., et al. 2012, AJ, 144, 120

Michaud, G., Tarasick, D., Charland, Y., & Pelletier, C. 1983, ApJ, 269, 239

Millholland, S., Wang, S., & Laughlin, G. 2016, ApJL, 823, L7

Millholland, S., Wang, S., & Laughlin, G. 2017, ApJL, 849, L33

Mills, S. M., Howard, A. W., Weiss, L. M., et al. 2019, AJ, 157, 145

Mol Louis, M., & Miguel, Y. 2020, MNRAS, 495, 294

Morris, R. L., Twicken, J. D., Smith, J. C., et al. 2020, Kepler Data Processing Handbook (Moffett Field, CA: NASA Ames Research Center), 87

Morton, T. D., Bryson, S. T., Coughlin, J. L., et al. 2016, ApJ, 822, 86

Naef, D., Latham, D. W., Mayor, M., et al. 2001, A&A, 375, L27

Nagasawa, M., Ida, S., & Bessho, T. 2008, ApJ, 678, 498

Narita, N., Fukui, A., Kusakabe, N., et al. 2019, JATIS, 5, 015001

Narita, N., Sato, B., Hirano, T., & Tamura, M. 2009, PASJ, 61, L35

Newton, E. R., Mann, A. W., Tofflemire, B. M., et al. 2019, ApJL, 880, L17

Parmentier, V., & Crossfield, I. J. M. 2018, Exoplanet Phase Curves: Observations and Theory (Berlin: Springer), 116

Parmentier, V., Line, M. R., Bean, J. L., et al. 2018, A&A, 617, A110

Parviainen, H., Tingley, B., Deeg, H. J., et al. 2019, A&A, 630, A89

Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3

Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4

Paxton, B., Marchant, P., Schwab, J., et al. 2015, ApJS, 220, 15

Paxton, B., Schwab, J., Bauer, E., et al. 2018, ApJS, 234, 34

Paxton, B., Smolec, R., Schwab, J., et al. 2019, ApJS, 243, 10

Pepe, F., Mayor, M., Queloz, D., et al. 2004, A&A, 423, 385

Perruchot, S., Bouchy, F., Chazelas, B., et al. 2011, Proc. SPIE, 8151, 815115

Perruchot, S., Kohler, D., Bouchy, F., et al. 2008, Proc. SPIE, 7010, 70140J

Perryman, M. A. C., Brown, A. G. A., Lebreton, Y., et al. 1998, A&A, 331, 81

Petersburg, R. R., Ong, J. M. J., Zhao, L. L., et al. 2020, AJ, 159, 187

Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525

Plavchan, P., Barclay, T., Gagné, J., et al. 2020, Natur, 582, 497

Preston, G. W. 1974, ARA&A, 12, 257

Reffert, S., Bergmann, C., Quirrenbach, A., Trifonov, T., & Künstler, A. 2015, A&A, 574, A116

Renson, P., Gerbaldi, M., & Catalano, F. A. 1991, A&AS, 89, 429

Renson, P., & Manfroid, J. 2009, A&A, 498, 961

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, JATIS, 1, 014003

Ritter, A., Hyde, E. A., & Parker, Q. A. 2014, PASP, 126, 170

Rodriguez, J. E., Vanderburg, A., Zieba, S., et al. 2020, AJ, 160, 117

Rodríguez Martínez, R., Gaudí, B. S., Rodriguez, J. E., et al. 2020, AJ, 160, 111

Rogers, L. A. 2015, ApJ, 801, 41

Rossiter, R. A. 1924, ApJ, 60, 15

Ryabchikova, T., Piskunov, N., Kurucz, R. L., et al. 2015, PhyS, 90, 054005

Saffe, C., Miquelarena, P., Alacoria, J., et al. 2020, A&A, 641, A145

Saffe, C., Miquelarena, P., Alacoria, J., et al. 2021, A&A, 647, 49

Schlaflay, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103

Spórker, A. 2017, PASP, 129, 072001

Spórker, A., Wong, I., Huang, C. X., et al. 2019, AJ, 157, 178

Siverd, R. J., Brown, T. M., Barnes, S., et al. 2018b, Proc. SPIE, 10702, 107026C

Siverd, R. J., Brown, T. M., Barnes, S., et al. 2018c, AJ, 155, 35

Siverd, R. J., Collins, K. A., Zhou, G., et al. 2018a, AJ, 155, 35

Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000

Sneden, C. 1973, ApJ, 184, 839

Sozzetti, A., Damasso, M., Bonomo, A. S., et al. 2021, A&A, 648, 75

Sozzetti, A., Torres, G., Charbonneau, D., et al. 2007, ApJ, 664, 1190

Stangret, M., Pallé, E., Casasayas-Barris, N., et al. 2021, A&A, 654, 73

Stassun, K. G., Oelkers, R. J., Paegert, M., et al. 2019, AJ, 158, 138

Stassun, K. G., & Torres, G. 2018, ApJ, 862, 61

Stumpe, M. C., Smith, J. C., Catanzarite, J. H., et al. 2014, PASP, 126, 100

Stumpe, M. C., Smith, J. C., Van Cleve, J. E., et al. 2012, PASP, 124, 985

Talens, G. J. J., Justesen, A. B., Albrecht, S., et al. 2018, A&A, 612, A57

Talens, G. J. J., Spronck, J. F. P., Lesage, A. L., et al. 2017, A&A, 601, A11

Tanaka, H., Takeuchi, T., & Ward, W. R. 2002, ApJ, 565, 1257

Telting, J. H., Avila, G., Buchhave, L., et al. 2014, AN, 335, 41

Temple, L. Y., Hellier, C., Albrow, M. D., et al. 2017, MNRAS, 471, 2743

Teske, J., Díaz, M. R., Luque, R., et al. 2020, AJ, 160, 96

Thorngren, D., & Fortney, J. 2018, AJ, 155, 214

Tinetti, G., Drossart, P., Eccleston, P., et al. 2018, ExA, 46, 135

Tinney, C. G., Butler, R. P., Marcy, G. W., et al. 2001, ApJ, 551, 507

Torres, G., Fischer, D. A., Sozzetti, A., et al. 2012, ApJ, 757, 161
Triaud, A. H. M. J., Collier Cameron, A., Queloz, D., et al. 2010, A&A, 524, A25

Trifonov, T., Reffert, S., Tan, X., Lee, M. H., & Quirrenbach, A. 2014, A&A, 568, A64

Tsantaki, M., Andreasen, D. T., Teixeira, G. D. C., et al. 2018, MNRAS, 473, 5066

Twicken, J. D., Catanzarite, J. H., Clarke, B. D., et al. 2018, PASP, 130, 064502

Twicken, J. D., Clarke, B. D., Bryson, S. T., et al. 2010, Proc. SPIE, 7740, 774023

Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141

Vogt, S. S., Wittenmyer, R. A., Butler, R. P., et al. 2010, ApJ, 708, 1366

von Essen, C., Mallonn, M., Borre, C. C., et al. 2020, A&A, 639, A34

Wakeford, H. R., Visscher, C., Lewis, N. K., et al. 2017, MNRAS, 464, 4247

Wang, S., Jones, M., Shporer, A., et al. 2019, AJ, 157, 51

Ward, W. R. 1997, Icar, 126, 261

Weiss, L. M., Fabrycky, D. C., Agol, E., et al. 2020, AJ, 159, 242

Winn, J. N., Fabrycky, D., Albrecht, S., & Johnson, J. A. 2010, ApJL, 718, L145

Winn, J. N., Holman, M. J., Torres, G., et al. 2008, ApJ, 683, 1076

Winn, J. N., Johnson, J. A., Albrecht, S., et al. 2009, ApJL, 703, L99

Wittenmyer, R. A., Butler, R. P., Horner, J., et al. 2020, MNRAS, 491, 5248

Wittenmyer, R. A., Clark, J. T., Zhao, J., et al. 2019, MNRAS, 484, 5859

Wittenmyer, R. A., Jones, M. I., Horner, J., et al. 2017, AJ, 154, 274

Wong, I., Benneke, B., Shporer, A., et al. 2020a, AJ, 159, 104

Wong, I., Shporer, A., Daylan, T., et al. 2020c, AJ, 160, 155

Wong, I., Shporer, A., Kitzmann, D., et al. 2020b, AJ, 160, 88

Wright, J. T., Marcy, G. W., Butler, R. P., et al. 2008, ApJL, 683, L63

Xiang, M.-S., Rix, H.-W., Ting, Y.-S., et al. 2020, ApJ, 898, 28

Yi, S., Demarque, P., Kim, Y.-C., et al. 2001, ApJS, 136, 417

Zhou, G., Bakos, G. Á., Bayliss, D., et al. 2019b, AJ, 157, 31

Zhou, G., Huang, C. X., Bakos, G. Á., et al. 2019a, AJ, 158, 141

Zhou, G., Rodriguez, J. E., Collins, K. A., et al. 2016, AJ, 152, 136

Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481