Search

Apsidal precession in stellar binaries is the main non-Keplerian dynamical effect impacting the radial-velocities of a binary star system. Its presence can notably hide the presence of orbiting circumbinary planets because many fitting algorithms assume perfectly Keplerian motion. To first order, apsidal precession ($\dot{\omega}$) can be accounted for by adding a linear term to the usual Keplerian model. We include apsidal precession in the kima package, an orbital fitter designed to detect and characterise planets from radial velocity data. In this paper, we detail this and other additions to kima that improve fitting for stellar binaries and circumbinary planets including corrections from general relativity. We then demonstrate that fitting for $\dot{\omega}$ can improve the detection sensitivity to circumbinary exoplanets by up to an order of magnitude in some circumstances, particularly in the case of multi-planetary systems. In addition, we apply the algorithm to several real systems, producing a new measurement of aspidal precession in KOI-126 (a tight triple system), and a detection of $\dot{\omega}$ in the Kepler-16 circumbinary system. Although apsidal precession is detected for Kepler-16, it does not have a large effect on the detection limit or the planetary parameters. We also derive an expression for the precession an outer planet would induce on the inner binary and compare the value this predicts with the one we detect., Comment: 10 pages, 5 figures Re-submitted to MNRAS after reviewer comments

The angle between the rotation and orbital axes of stars in binary systems -- the obliquity -- is an important indicator of how these systems form and evolve but few such measurements exist. We combine the sample of astrometric orbital inclinations from Gaia DR3 with a sample of solar-like oscillators in which rotational inclinations have been measured using asteroseismology. We supplement our sample with one binary whose visual orbit has been determined using speckle interferometry and present the projected spin-orbit alignments in five systems. We find that each system, and the overall sample, is consistent with alignment but there are important caveats. First, the asteroseismic rotational inclinations are fundamentally ambiguous and, second, we can only measure the projected (rather than true) obliquity. If rotational and orbital inclinations are independent and isotropically-distributed, however, the likelihood of drawing our data by chance is less than a few per cent. Though small, our data set argues against uniformly random obliquities in binary systems. We speculate that dozens more measurements could be made using data from NASA's TESS mission, mostly in red giants. ESA's PLATO mission will likely produce hundreds more spin-orbit measurements in systems with main-sequence and subgiant stars., Comment: Accepted for publication in MNRAS. 4 pages, 3 figures, 1 table. Scripts to reproduce Figs 1 & 2 and data in Table 1 are available at https://gitlab.com/warrickball/kepler-gaia-spin-orbit

TOI-1259 consists of a transiting exoplanet orbiting a main sequence star, with a bound outer white dwarf companion. Less than a dozen systems with this architecture are known. We conduct follow-up spectroscopy on the white dwarf TOI-1259B using the Large Binocular Telescope (LBT) to better characterise it. We observe only strong hydrogen lines, making TOI-1259B a DA white dwarf. We see no evidence of heavy element pollution, which would have been evidence of planetary material around the white dwarf. Such pollution is seen in ~ 25 - 50% of white dwarfs, but it is unknown if this rate is higher or lower in TOI-1259-like systems that contain a known planet. Our spectroscopy permits an improved white dwarf age measurement of 4.05 (+1.00 -0.42) Gyrs, which matches gyrochronology of the main sequence star. This is the first of an expanded sample of similar binaries that will allow us to calibrate these dating methods and provide a new perspective on planets in binaries., Comment: 6 pages, accepted to MNRAS, minor revision of first arXiv version

Numerous circumbinary planets have been discovered in surveys of transiting planets. Often, these planets are found to orbit near to the zone of dynamical instability, close to the central binary. The existence of these planets has been explained by hydrodynamical simulations that show that migrating circumbinary planets, embedded in circumbinary discs, halt at the central cavity that is formed by the central binary. Transit surveys are naturally most sensitive to finding circumbinary planets with the shortest orbital periods. The future promise of detecting longer period systems using radial-velocity searches, combined with the anticipated detection of numerous circumbinary planets by ESA's PLATO mission, points to the need to model and understand the formation and evolution of circumbinary planets in a more general sense than has been considered before. With this goal in mind, we present a newly developed global model of circumbinary planet formation that is based on the mercury6 symplectic N-body integrator, combined with a model for the circumbinary disc and prescriptions for a range of processes involved in planet formation such as pebble accretion, gas envelope accretion and migration. Our results show that under reasonable assumptions, the pebble accretion scenario can produce circumbinary systems that are similar to those observed, and in particular is able to produce planets akin to Kepler-16b and Kepler-34b. Comparing our results to other systems, we find that our models also adequately reproduce such systems, including multi-planet systems. Resonances between neighbouring planets are frequently obtained, whilst ejections of planets by the central binary acts as an effective source of free floating planets., Accepted for publication in MNRAS, 23 pages, 16 figures

We present the results of two-fluid hydrodynamical simulations of circumbinary discs consisting of gas and dust, with and without embedded planets, to examine the influence of the dust on the structure of the tidally truncated inner cavity and on the parking locations of migrating planets. In this proof-of-concept study, we consider Kepler-16 and -34 analogues, and examine dust fluids with Stokes numbers in the range $10^{-4} \le St \le 10^{-1}$ and dust-to-gas ratios of 0.01 and 1. For the canonical dust-to-gas ratio of 0.01, we find the inclusion of the dust has only a minor effect on the cavity and stopping locations of embedded planets compared to dust-free simulations. However, for the enhanced dust-to-gas ratio of unity, assumed to arise because of significant dust drift and accumulation, we find that the dust can have a dramatic effect by shrinking and circularising the inner cavity, which brings the parking locations of planets closer to the central binary. This work demonstrates the importance of considering both gas and dust in studies of circumbinary discs and planets, and provides a potential means of explaining the orbital properties of circumbinary planets such as Kepler-34b, which have hitherto been difficult to explain using gas-only hydrodynamical simulations., Comment: Accepted for publication in MNRAS, 20 pages, 17 figures

We present the discovery and characterisation of an eclipsing binary identified by the Next Generation Transit Survey in the $\sim$115 Myr old Blanco 1 open cluster. NGTS J0002-29 comprises three M dwarfs: a short-period binary and a companion in a wider orbit. This system is the first well-characterised, low-mass eclipsing binary in Blanco 1. With a low mass ratio, a tertiary companion and binary components that straddle the fully convective boundary, it is an important benchmark system, and one of only two well-characterised, low-mass eclipsing binaries at this age. We simultaneously model light curves from NGTS, TESS, SPECULOOS and SAAO, radial velocities from VLT/UVES and Keck/HIRES, and the system's spectral energy distribution. We find that the binary components travel on circular orbits around their common centre of mass in $P_{\rm orb} = 1.09800524 \pm 0.00000038$ days, and have masses $M_{\rm pri}=0.3978\pm 0.0033$ M$_{\odot}$ and $M_{\rm sec}=0.2245\pm 0.0018$ M$_{\odot}$, radii $R_{\rm pri}=0.4037\pm 0.0048$ R$_{\odot}$ and $R_{\rm sec}=0.2759\pm 0.0055$ R$_{\odot}$, and effective temperatures $T_{\rm pri}=3372\,^{+44}_{-37}$ K and $T_{\rm sec}=3231\,^{+38}_{-31}$ K. We compare these properties to the predictions of seven stellar evolution models, which typically imply an inflated primary. The system joins a list of 19 well-characterised, low-mass, sub-Gyr, stellar-mass eclipsing binaries, which constitute some of the strongest observational tests of stellar evolution theory at low masses and young ages., 21 pages, 9 figures. Accepted for publication in MNRAS

We present the discovery of TOI-2136b, a sub-Neptune planet transiting every 7.85 days a nearby M4.5V-type star, identified through photometric measurements from the TESS mission. The host star is located $33$ pc away with a radius of $R_{\ast} = 0.34\pm0.02\ R_{\odot}$, a mass of $0.34\pm0.02\ M_{\odot}$ and an effective temperature of $\rm 3342\pm100\ K$. We estimate its stellar rotation period to be $75\pm5$ days based on archival long-term photometry. We confirm and characterize the planet based on a series of ground-based multi-wavelength photometry, high-angular-resolution imaging observations, and precise radial velocities from CFHT/SPIRou. Our joint analysis reveals that the planet has a radius of $2.19\pm0.17\ R_{\oplus}$, and a mass measurement of $6.4\pm2.4\ M_{\oplus}$. The mass and radius of TOI2136b is consistent with a broad range of compositions, from water-ice to gas-dominated worlds. TOI-2136b falls close to the radius valley for low-mass stars predicted by the thermally driven atmospheric mass loss models, making it an interesting target for future studies of its interior structure and atmospheric properties., Comment: 19 pages, 15 figures, submitted to MNRAS

HIP 41378 f is a temperate $9.2\pm0.1 R_{\oplus}$ planet with period of 542.08 days and an extremely low density of $0.09\pm0.02$ g cm$^{-3}$. It transits the bright star HIP 41378 (V=8.93), making it an exciting target for atmospheric characterization including transmission spectroscopy. HIP 41378 was monitored photometrically between the dates of 2019 November 19 and November 28. We detected a transit of HIP 41378 f with NGTS, just the third transit ever detected for this planet, which confirms the orbital period. This is also the first ground-based detection of a transit of HIP 41378 f. Additional ground-based photometry was also obtained and used to constrain the time of the transit. The transit was measured to occur 1.50 hours earlier than predicted. We use an analytic transit timing variation (TTV) model to show the observed TTV can be explained by interactions between HIP 41378 e and HIP 41378 f. Using our TTV model, we predict the epochs of future transits of HIP 41378 f, with derived transit centres of T$_{C,4} = 2459355.087^{+0.031}_{-0.022}$ (May 2021) and T$_{C,5} = 2459897.078^{+0.114}_{-0.060}$ (Nov 2022)., Comment: Accepted for publication in MNRAS Letters. 6 pages, 2 figures

Colour-Magnitude Diagrams provide a convenient way of comparing populations of similar objects. When well populated with precise measurements, they allow quick inferences to be made about the bulk properties of an astronomic object simply from its proximity on a diagram to other objects. We present here a Python toolkit which allows a user to produce colour-magnitude diagrams of transiting exoplanets, comparing planets to populations of ultra-cool dwarfs, of directly imaged exoplanets, to theoretical models of planetary atmospheres, and to other transiting exoplanets. Using a selection of near- and mid-infrared colour-magnitude diagrams, we show how outliers can be identified for further investigation, and how emerging sub-populations can be identified. Additionally, we present evidence that observed differences in the \textit{Spitzer}'s 4.5\mu m flux, between irradiated Jupiters, and field brown dwarfs, might be attributed to phosphine, which is susceptible to photolysis. The presence of phosphine in low irradiation environments may negate the need for thermal inversions to explain eclipse measurements. We speculate that the anomalously low 4.5\mu m flux flux of the nightside of HD 189733b and the daysides of GJ 436b and GJ 3470b might be caused by phosphine absorption. Finally, we use our toolkit to include \textit{Hubble} WFC3 spectra, creating a new photometric band called the `Water band' (\textit{W$_{JH}$}-band) in the process. We show that the colour index [\textit{W$_{JH}$-H}] can be used to constrain the C/O ratio of exoplanets, showing that future observations with \textit{JWST} and \textit{Ariel} will be able to distinguish these populations if they exist, and select members for future follow-up., Comment: Accepted for publication in MNRAS

A dozen short-period detached binaries are known to host transiting circumbinary planets. In all circumbinary systems so far, the planetary and binary orbits are aligned within a couple of degrees. However, the obliquity of the primary star, which is an important tracer of their formation, evolution, and tidal history, has only been measured in one circumbinary system until now. EBLM J0608-59/TOI-1338 is a low-mass eclipsing binary system with a recently discovered circumbinary planet identified by TESS. Here, we perform high-resolution spectroscopy during primary eclipse to measure the projected stellar obliquity of the primary component. The obliquity is low, and thus the primary star is aligned with the binary and planetary orbits with a projected spin-orbit angle $\beta = 2.8 \pm 17.1$ deg. The rotation period of $18.1 \pm 1.6$ days implied by our measurement of $v\sin{i_\star}$ suggests that the primary has not yet pseudo-synchronized with the binary orbit, but is consistent with gyrochronology and weak tidal interaction with the binary companion. Our result, combined with the known coplanarity of the binary and planet orbits, is suggestive of formation from a single disc. Finally, we considered whether the spectrum of the faint secondary star could affect our measurements. We show through simulations that the effect is negligible for our system, but can lead to strong biases in $v\sin{i_\star}$ and $\beta$ for higher flux ratios. We encourage future studies in eclipse spectroscopy test the assumption of a dark secondary for flux ratios $\gtrsim 1$ ppt., Comment: 7 pages, 3 figures. Accepted in MNRAS. Fixed a few typos

A 2014 study of the eclipsing binary star 1SWASPJ011351.29+314909.7 (J0113+31) reported an unexpectedly high effective temperature for the M-dwarf companion to the 0.95-M$_{\odot}$ primary star. The effective temperature inferred from the secondary eclipse depth was $\sim$600 K higher than the value predicted from stellar models. Such an anomalous result questions our understanding of low-mass stars and might indicate a significant uncertainty when inferring properties of exoplanets orbiting them. We seek to measure the effective temperature of the M-dwarf companion using the light curve of J0113+31 recently observed by the Transiting Exoplanet Survey Satellite (TESS). We use the pycheops modelling software to fit a combined transit and eclipse model to the TESS light curve. To calculate the secondary effective temperature, we compare the best-fit eclipse depth to the predicted eclipse depths from theoretical stellar models. We determined the effective temperature of the M dwarf to be ${\rm T}_{\rm eff,2}$ = 3208 $\pm$ 43 K, assuming $\log g_2$ = 5, [Fe/H] = $-0.4$ and no alpha-element enhancement. Varying these assumptions changes ${\rm T}_{\rm eff,2}$ by less than 100 K. These results do not support a large anomaly between observed and theoretical low-mass star temperatures., 5 pages, 3 figures, published in MNRAS

BEBOP is a radial-velocity survey that monitors a sample of single-lined eclipsing binaries, in search of circumbinary planets by using high-resolution spectrographs. Here, we describe and test the methods we use to identify planetary signals within the BEBOP data, and establish how we quantify our sensitivity to circumbinary planets by producing detection limits. This process is made easier and more robust by using a diffusive nested sampler. In the process of testing our methods, we notice that contrary to popular wisdom, assuming circular orbits in calculating detection limits for a radial velocity survey provides over-optimistic detection limits by up to $40\%$ in semi-amplitude with implications for all radial-velocity surveys. We perform example analyses using three BEBOP targets from our Southern HARPS survey. We demonstrate for the first time a repeated ability to reach a residual root mean squared scatter of $3~\rm m.s^{-1}$ (after removing the binary signal), and find we are sensitive to circumbinary planets with masses down to that of Neptune and Saturn, for orbital periods up to $1000~\rm days$., Comment: 17 pages, 11 figures

M-dwarfs are the most abundant stars in the galaxy and popular targets for exoplanet searches. However, their intrinsic faintness and complex spectra inhibit precise characterisation. We only know of dozens of M-dwarfs with fundamental parameters of mass, radius and effective temperature characterised to better than a few per cent. Eclipsing binaries remain the most robust means of stellar characterisation. Here we present two targets from the Eclipsing Binary Low Mass (EBLM) survey that were observed with K2: EBLM J0055-00 and EBLM J2217-04. Combined with HARPS and CORALIE spectroscopy, we measure M-dwarf masses with precisions better than 5%, radii better than 3% and effective temperatures on order 1%. However, our fits require invoking a model to derive parameters for the primary star. By investigating three popular models, we determine that the model uncertainty is of similar magnitude to the statistical uncertainty in the model fits. Therefore, whilst these can be considered benchmark M-dwarfs, we caution the community to consider model uncertainty when pushing the limits of precise stellar characterisation., Comment: 13 Pages, MNRAS Submission, Comments welcome

We present ground- and space-based photometric observations of TOI-270 (L231-32), a system of three transiting planets consisting of one super-Earth and two sub-Neptunes discovered by TESS around a bright (K-mag = 8.25) M3V dwarf. The planets orbit near low-order mean-motion resonances (5:3 and 2:1) and are thus expected to exhibit large transit timing variations (TTVs). Following an extensive observing campaign using eight different observatories between 2018 and 2020, we now report a clear detection of TTVs for planets c and d, with amplitudes of ∼10 min and a super-period of ∼3 yr, as well as significantly refined estimates of the radii and mean orbital periods of all three planets. Dynamical modelling of the TTVs alone puts strong constraints on the mass ratio of planets c and d and on their eccentricities. When incorporating recently published constraints from radial velocity observations, we obtain masses of $M_{\mathrm{b}}=1.48\pm 0.18\, M_\oplus$, $M_{\mathrm{c}}=6.20\pm 0.31\, M_\oplus$, and $M_{\mathrm{d}}=4.20\pm 0.16\, M_\oplus$ for planets b, c, and d, respectively. We also detect small but significant eccentricities for all three planets : eb = 0.0167 ± 0.0084, ec = 0.0044 ± 0.0006, and ed = 0.0066 ± 0.0020. Our findings imply an Earth-like rocky composition for the inner planet, and Earth-like cores with an additional He/H2O atmosphere for the outer two. TOI-270 is now one of the best constrained systems of small transiting planets, and it remains an excellent target for atmospheric characterization.

The accuracy of theoretical mass, radius and effective temperature values for M-dwarf stars is an active topic of debate. Differences between observed and theoretical values have raised the possibility that current theoretical stellar structure and evolution models are inaccurate towards the low-mass end of the main sequence. To explore this issue we use the CHEOPS satellite to obtain high-precision light curves of eclipsing binaries with low mass stellar companions. We use these light curves combined with the spectroscopic orbit for the solar-type companion to measure the mass, radius and effective temperature of the M-dwarf star. Here we present the analysis of three eclipsing binaries. We use the pycheops data analysis software to fit the observed transit and eclipse events of each system. Two of our systems were also observed by the TESS satellite -- we similarly analyse these light curves for comparison. We find consistent results between CHEOPS and TESS, presenting three stellar radii and two stellar effective temperature values of low-mass stellar objects. These initial results from our on-going observing programme with CHEOPS show that we can expect to have ~24 new mass, radius and effective temperature measurements for very low mass stars within the next few years., Comment: 12 pages, 8 figures, accepted for publication in MNRAS

We report on four new transiting hot Jupiters discovered by the WASP-South survey. WASP-178b transits a V = 9.9, A1V star with Teff = 9350 +/- 150 K, the second-hottest transit host known. It has a highly bloated radius of 1.81 +/- 0.09 Rjup, in line with the known correlation between high irradiation and large size. With an estimated temperature of 2470 +/- 60 K, the planet is one of the best targets for studying ultra-hot Jupiters that is visible from the Southern hemisphere. The three host stars WASP-184, WASP-185 and WASP-192 are all post-main-sequence G0 stars of ages 4-8 Gyr. The larger stellar radii (1.3-1.7 Msun) mean that the transits are relatively shallow (0.7-0.9%) even though the planets have moderately inflated radii of 1.2-1.3 Rjup. WASP-185b has an eccentric orbit (e = 0.24) and a relatively long orbital period of 9.4 d. A star that is 4.6 arcsec from WASP-185 and 4.4 mag fainter might be physically associated., Version as Accepted for MNRAS (RV table added)

We report the discovery and characterisation of WASP-180Ab, a hot Jupiter confirmed by the detection of its Doppler shadow and by measuring its mass using radial velocities. We find the 0.9 $\pm$ 0.1 $M_{\rm Jup}$, 1.24 $\pm$ 0.04 $R_{\rm Jup}$ planet to be in a misaligned, retrograde orbit around an F7 star with $T_{\rm eff}$ = 6500K and a moderate rotation speed of vsini = 19.9 km s$^{-1}$. The host star is the primary of a $V$ = 10.7 binary, where a secondary separated by 5$''$ ($\sim$1200 AU) contributes $\sim$30% of the light. WASP-180Ab therefore adds to a small sample of transiting hot Jupiters known in binary systems. A 4.6-day modulation seen in the WASP data is likely to be the rotational modulation of the companion star, WASP-180B., Comment: 9 pages, 6 figures, accepted by MNRAS

NS acknowledges the support of the National Priorities Research Program (NPRP) grant #X-019-1-006 from the Qatar National Research Fund (a member of Qatar Foundation). ACC acknowledges support from the Science & Technology Facilities Council (STFC) consolidated grant number ST/R000824/1. SuperWASP, the Northern hemisphere WASP observatory, has been observing the skies from La Palma since 2004. In that time, more than 50 planets have been discovered with data contributions from SuperWASP. In the process of validating planets, many false-positive candidates have also been identified. The TESS telescope is set to begin observations of the northern sky in 2019. Similar to the WASP survey, the TESS pixel size is relatively large (13arcsec for WASP and 21arcsec for TESS), making it susceptible to many blended signals and false detections caused principally by grazing and blended stellar eclipsing binary systems. In order to reduce duplication of effort on targets, we present a catalogue of 1041 Northern hemisphere SuperWASP targets that have been rejected as planetary transits through follow-up observation. Publisher PDF

Plate tectonics is a geophysical process currently unique to Earth, has an important role in regulating the Earth's climate, and may be better understood by identifying rocky planets outside our solar system with tectonic activity. The key criterion for whether or not plate tectonics may occur on a terrestrial planet is if the stress on a planet's lithosphere from mantle convection may overcome the lithosphere's yield stress. Although many rocky exoplanets closely orbiting their host stars have been detected, all studies to date of plate tectonics on exoplanets have neglected tidal stresses in the planet's lithosphere. Modeling a rocky exoplanet as a constant density, homogeneous, incompressible sphere, we show the tidal stress from the host star acting on close-in planets may become comparable to the stress on the lithosphere from mantle convection. Tidal stress of this magnitude may aid mantle convection stress in subduction of plates, or drive the subduction of plates without the need for mantle convective stresses. We also show that tidal stresses from planet-planet interactions are unlikely to be significant for plate tectonics, but may be strong enough to trigger Earthquakes. Our work may imply planets orbiting close to their host stars are more likely to experience plate tectonics, with implications for exoplanetary geophysics and habitability. We produce a list of detected rocky exoplanets under the most intense stresses. Atmospheric and topographic observations may confirm our predictions in the near future. Investigations of planets with significant tidal stress can not only lead to observable parameters linked to the presence of active plate tectonics, but may also be used as a tool to test theories on the main driving force behind tectonic activity.

The TRAPPIST-1 planetary system is a favorable target for the atmospheric characterization of temperate earth-sized exoplanets by means of transmission spectroscopy with the forthcoming James Webb Space Telescope (JWST). A possible obstacle to this technique could come from the photospheric heterogeneity of the host star that could affect planetary signatures in the transit transmission spectra. To constrain further this possibility, we gathered an extensive photometric data set of 25 TRAPPIST-1 transits observed in the near-IR J band (1.2 $\mu$m) with the UKIRT and the AAT, and in the NB2090 band (2.1 $\mu$m) with the VLT during the period 2015-2018. In our analysis of these data, we used a special strategy aiming to ensure uniformity in our measurements and robustness in our conclusions. We reach a photometric precision of $\sim0.003$ (RMS of the residuals), and we detect no significant temporal variations of transit depths of TRAPPIST-1 b, c, e, and g over the period of three years. The few transit depths measured for planets d and f hint towards some level of variability, but more measurements will be required for confirmation. Our depth measurements for planets b and c disagree with the stellar contamination spectra originating from the possible existence of bright spots of temperature 4500 K. We report updated transmission spectra for the six inner planets of the system which are globally flat for planets b and g and some structures are seen for planets c, d, e, and f., Comment: accepted for publication in MNRAS

We present the discovery of 3 transiting planets from the WASP survey, two hot-Jupiters: WASP-177b (~0.5 M_Jup, ~1.6 R_Jup) in a 3.07-d orbit of a V = 12.6 K2 star, WASP-183b (~0.5 M_Jup, ~1.5 R_Jup) in a 4.11-d orbit of a V = 12.8 G9/K0 star; and one hot-Saturn planet WASP-181b (~0.3 M_Jup, ~1.2 R_Jup) in a 4.52-d orbit of a V = 12.9 G2 star. Each planet is close to the upper bound of mass-radius space and has a scaled semi-major axis, a/R_star, between 9.6 and 12.1. These lie in the transition between systems that tend to be in orbits that are well aligned with their host-star's spin and those that show a higher dispersion., 11 pages, 9 figures, 5 tables

We report the discovery of a transiting, temperate, Neptune-sized exoplanet orbiting the nearby ($d$ = 27.5 pc), M3V star TOI-1231 (NLTT 24399, L 248-27, 2MASS J10265947-5228099). The planet was detected using photometric data from the Transiting Exoplanet Survey Satellite and followed up with observations from the Las Cumbres Observatory and the Antarctica Search for Transiting ExoPlanets program. Combining the photometric data sets, we find that the newly discovered planet has a radius of 3.65$^{+0.16}_{-0.15}$ R$_{\oplus}$, and an orbital period of 24.246 days. Radial velocity measurements obtained with the Planet Finder Spectrograph on the Magellan Clay telescope confirm the existence of the planet and lead to a mass measurement of 15.5$\pm$3.3 M$_{\oplus}$. With an equilibrium temperature of just 330K TOI-1231 b is one of the coolest small planets accessible for atmospheric studies thus far, and its host star's bright NIR brightness (J=8.88, K$_{s}$=8.07) make it an exciting target for HST and JWST. Future atmospheric observations would enable the first comparative planetology efforts in the 250-350 K temperature regime via comparisons with K2-18 b. Furthermore, TOI-1231's high systemic radial velocity (70.5 k\ms) may allow for the detection of low-velocity hydrogen atoms escaping the planet by Doppler shifting the H I Ly-alpha stellar emission away from the geocoronal and ISM absorption features., 20 pages, 9 figures. Accepted for publication in The Astronomical Journal

Warm Jupiters -- defined here as planets larger than 6 Earth radii with orbital periods of 8--200 days -- are a key missing piece in our understanding of how planetary systems form and evolve. It is currently debated whether Warm Jupiters form in situ, undergo disk or high eccentricity tidal migration, or have a mixture of origin channels. These different classes of origin channels lead to different expectations for Warm Jupiters' properties, which are currently difficult to evaluate due to the small sample size. We take advantage of the \TESS survey and systematically search for Warm Jupiter candidates around main-sequence host stars brighter than the \TESS-band magnitude of 12 in the Full-Frame Images in Year 1 of the \TESS Prime Mission data. We introduce a catalog of 55 Warm Jupiter candidates, including 19 candidates that were not originally released as \TESS Objects of Interest (TOIs) by the \TESS team. We fit their \TESS light curves, characterize their eccentricities and transit-timing variations (TTVs), and prioritize a list for ground-based follow-up and \TESS Extended Mission observations. Using hierarchical Bayesian modeling, we find the preliminary eccentricity distributions of our Warm-Jupiter-candidate catalog using a Beta distribution, a Rayleigh distribution, and a two-component Gaussian distribution as the functional forms of the eccentricity distribution. Additional follow-up observations will be required to clean the sample of false positives for a full statistical study, derive the orbital solutions to break the eccentricity degeneracy, and provide mass measurements., 30 pages, 8 figures, 6 tables. submitted to ApJS, revised in response to referee report

TOI-216 hosts a pair of warm, large exoplanets discovered by the TESS Mission. These planets were found to be in or near the 2:1 resonance, and both of them exhibit transit timing variations (TTVs). Precise characterization of the planets' masses and radii, orbital properties, and resonant behavior can test theories for the origins of planets orbiting close to their stars. Previous characterization of the system using the first six sectors of TESS data suffered from a degeneracy between planet mass and orbital eccentricity. Radial velocity measurements using HARPS, FEROS, and PFS break that degeneracy, and an expanded TTV baseline from TESS and an ongoing ground-based transit observing campaign increase the precision of the mass and eccentricity measurements. We determine that TOI-216c is a warm Jupiter, TOI-216b is an eccentric warm Neptune, and that they librate in the 2:1 resonance with a moderate libration amplitude of 60 +/- 2 degrees; small but significant free eccentricity of 0.0222 +0.0005/-0.0003 for TOI-216b; and small but significant mutual inclination of 1.2-3.9 degrees (95% confidence interval). The libration amplitude, free eccentricity, and mutual inclination imply a disturbance of TOI-216b before or after resonance capture, perhaps by an undetected third planet., Comment: AJ accepted

Leleu, A., et al., Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at the possible presence of a near 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152to 2.87Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02to 0.177times the Earth's density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 (H = 8.76 mag, J = 9.37 mag, V = 11.95 mag) allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes., The authors acknowledge support from the Swiss NCCR PlanetS and the Swiss National Science Foundation. Y.A. and M.J.H. acknowledge the support of the Swiss National Fund under grant 200020_172746. A.C.C. and T.W. acknowledge support from STFC consolidated grant number ST/M001296/1. This work was granted access to the HPC resources of MesoPSL financed by the Region Ile de France and the project Equip@Meso (reference ANR-10-EQPX-29-01) of the programme Investissements d’Avenir supervised by the Agence Nationale pour la Recherche. SH acknowledges CNES funding through the grant 837319. Based on data collected under the NGTS project at the ESO La Silla Paranal Observatory. The NGTS facility is operated by the consortium institutes with support from the UK Science and Technology Facilities Council (STFC) project ST/M001962/1. The Belgian participation to CHEOPS has been supported by the Belgian Federal Science Policy Office (BELSPO) in the framework of the PRODEX Program, and by the University of Liège through an ARC grant for Concerted Research Actions financed by the Wallonia-Brussels Federation. V.A. acknowledges the support from FCT through Investigador FCT contract nr. IF/00650/2015/CP1273/CT0001. We acknowledge support from the Spanish Ministry of Science and Innovation and the European Regional Development Fund through grants ESP2016-80435-C2-1-R, ESP2016-80435-C2-2-R, PGC2018-098153-B-C33, PGC2018-098153-B-C31, ESP2017-87676-C5-1-R, MDM-2017-0737 Unidad de Excelencia “María de Maeztu”- Centro de Astrobiología (INTA-CSIC), as well as the support of the Generalitat de Catalunya/CERCA programme. The MOC activities have been supported by the ESA contract No. 4000124370. S.C.C.B. acknowledges support from FCT through FCT contracts nr. IF/01312/2014/CP1215/CT0004. X.B., S.C., D.G., M.F. and J.L. acknowledge their role as ESA-appointed CHEOPS science team members. ABr was supported by the SNSA. A.C. acknowledges support by CFisUC projects (UIDB/04564/2020 and UIDP/04564/2020), GRAVITY (PTDC/FIS-AST/7002/2020), ENGAGE SKA (POCI-01-0145-FEDER-022217), and PHOBOS (POCI-01-0145-FEDER-029932), funded by COMPETE 2020 and FCT, Portugal. This work was supported by FCT - Fundaçãopara a Ciência e a Tecnologia through national funds and by FEDER through COMPETE2020 - Programa OperacionalCompetitividade 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. O.D.S.D. is supported in the form of work contract (DL 57/2016/CP1364/CT0004) funded by national funds through FCT. B.-O.D. acknowledges support from the Swiss National Science Foundation (PP00P2-190080). M.F. and C.M.P. gratefully acknowledgethe support of the Swedish National Space Agency (DNR 65/19, 174/18). D.G. gratefully acknowledges financial support from the CRT foundation under Grant No. 2018.2323 “Gaseousor rocky? Unveiling the nature of small worlds”. E.G. gratefully acknowledges support from the David and Claudia Harding Foundation in the form of a WintonExoplanet Fellowship. M.G. is an F.R.S.-FNRS Senior Research Associate. J.I.G.H. acknowledges financial support from Spanish Ministry of Science and Innovation (MICINN) under the 2013 Ramón y Cajal programme RYC-2013-14875. J.I.G.H., A.S.M., R.R., and C.A.P. acknowledge financial support from the Spanish MICINN AYA2017-86389-P. A.S.M. acknowledges financial support from the Spanish Ministry of Science and Innovation (MICINN) under the 2019 Juan de la Cierva Programme. MNG ackowledges support from the MIT Kavli Institute as a Juan Carlos Torres Fellow. J.H. acknowledges the support of the Swiss National Fund under grant 200020_172746. KGI is the ESA CHEOPS Project Scientist and is responsible for the ESA CHEOPS Guest Observers Programme.She does not participate in, or contribute to, the definition of the Guaranteed Time Programme of the CHEOPS mission through which observations described in this paper have been taken, nor to any aspect of target selection forthe programme. J.S.J. acknowledges support by FONDECYT grant 1201371, and partial support from CONICYT project Basal AFB-170002. A.J. acknowledges support from ANID - Millennium Science Initiative - ICN12_009 and from FONDECYT project 1171208. P.M. acknowledges support from STFC research grant number ST/M001040/1. N.J.N is supported by the contract and exploratory project IF/00852/2015, and projects UID/FIS/04434/2019, PTDC/FIS-OUT/29048/2017, COMPETE2020: POCI-01-0145-FEDER-028987 & FCT: PTDC/FIS-AST/28987/2017. N.J.N is supported by the contract and exploratory project IF/00852/2015, and project PTDC/FIS-OUT/29048/2017. This work was also partially supported by a grant from the Simons Foundation (PI Queloz, grant number 327127). Acknowledges support from the Spanish Ministry of Science and Innovation and the European Regional Development Fund through grant PGC2018-098153-B- C33, as well as the support of the Generalitat de Catalunya/CERCA programme. S.G.S. acknowledges support from FCT through FCT contract nr. CEECIND/00826/2018 and POPH/FSE (EC). This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This project has been supported by the Hungarian National Research, Development and Innovation Office (NKFIH) grants GINOP-2.3.2-15-2016-00003, K-119517, K-125015, and the City of Szombathely under Agreement No. 67.177-21/2016. This research received funding from the MERAC foundation, from the European Research Council under the European Union’s Horizon 2020 research and innovation programme (grant agreement no 803193/ BEBOP, and from the Science and Technology Facilites Council (STFC, grant no ST/S00193X/1). V.V.G. is an F.R.S-FNRS Research Associate. J.I.V. acknowledges support of CONICYT-PFCHA/Doctorado Nacional-21191829. M. R. Z. O. acknowledges financial support from projects AYA2016-79425-C3-2-P and PID2019-109522GB-C51 from the Spanish Ministry of Science, Innovation and Universities.

We report the first discovery of a transiting circumbinary planet detected from a single sector of TESS data. During Sector 21, the planet TIC 172900988b transited the primary star and then 5 days later it transited the secondary star. The binary is itself eclipsing, with a period of P = 19.7 days and an eccentricity of e = 0.45. Archival data from ASAS-SN, Evryscope, KELT, and SuperWASP reveal a prominent apsidal motion of the binary orbit, caused by the dynamical interactions between the binary and the planet. A comprehensive photodynamical analysis of the TESS, archival and follow-up data yields stellar masses and radii of M1 = 1.2384 +/- 0.0007 MSun and R1 = 1.3827 +/- 0.0016 RSun for the primary and M2 = 1.2019 +/- 0.0007 MSun and R2 = 1.3124 +/- 0.0012 RSun for the secondary. The radius of the planet is R3 = 11.25 +/- 0.44 REarth (1.004 +/- 0.039 RJup). The planet's mass and orbital properties are not uniquely determined - there are six solutions with nearly equal likelihood. Specifically, we find that the planet's mass is in the range of 824 < M3 < 981 MEarth (2.65 < M3 < 3.09 MJup), its orbital period could be 188.8, 190.4, 194.0, 199.0, 200.4, or 204.1 days, and the eccentricity is between 0.02 and 0.09. At a V = 10.141 mag, the system is accessible for high-resolution spectroscopic observations, e.g. Rossiter-McLaughlin effect and transit spectroscopy., 57 pages, 30 figures, 25 tables; Accepted AJ