With more than 1000 hours of observation from Feb 2016 to Oct 2019, the Spitzer Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12pc) ultracool dwarf star orbited by seven transiting Earth-sized planets, all well-suited for a detailed atmospheric characterization with the upcoming JWST. In this paper, we present the global results of the project. We analyzed 88 new transits and combined them with 100 previously analyzed transits, for a total of 188 transits observed at 3.6 or 4.5 $\mu$m. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5 $\mu$m to constrain the brightness temperatures of their daysides. We identify several orphan transit-like structures in our Spitzer photometry, but all of them are of low significance. We do not confirm any new transiting planets. We estimate for TRAPPIST-1 transit depth measurements mean noise floors of $\sim$35 and 25 ppm in channels 1 and 2 of Spitzer/IRAC, respectively. most of this noise floor is of instrumental origins and due to the large inter-pixel inhomogeneity of IRAC InSb arrays, and that the much better interpixel homogeneity of JWST instruments should result in noise floors as low as 10ppm, which is low enough to enable the atmospheric characterization of the planets by transit transmission spectroscopy. We construct updated broadband transmission spectra for all seven planets which show consistent transit depths between the two Spitzer channels. We identify and model five distinct high energy flares in the whole dataset, and discuss our results in the context of habitability. Finally, we fail to detect occultation signals of planets b and c at 4.5 $\mu$m, and can only set 3$\sigma$ upper limits on their dayside brightness temperatures (611K for b 586K for c).

We report on the discovery of a transiting Earth-sized (0.95$R_\oplus$) planet around an M3.5 dwarf star at 57\,pc, EPIC~249631677. The planet has a period of $\sim$3.14 days, i.e. ${\sim}\pi$, with an instellation of 7.5\,S$_{\oplus}$. The detection was made using publicly available data from \textit{K2}'s Campaign 15. We observed three additional transits with SPECULOOS Southern and Northern Observatories, and a stellar spectrum from Keck/HIRES, which allowed us to validate the planetary nature of the signal. The confirmed planet is well suited for comparative terrestrial exoplanetology. While exoplanets transiting ultracool dwarfs present the best opportunity for atmospheric studies of terrestrial exoplanets with the \textit{James Webb Space Telescope}, those orbiting mid-M dwarfs within 100\,pc such as EPIC~249631677b will become increasingly accessible with the next generation of observatories (e.g., \textit{HabEx, LUVOIR, OST}).

We report the discovery of two transiting exoplanets from the WASP survey, WASP-150b and WASP-176b. WASP-150b is an eccentric ($e$ = 0.38) hot Jupiter on a 5.6 day orbit around a $V$ = 12.03, F8 main-sequence host. The host star has a mass and radius of 1.4 $\rm M_{\odot}$ and 1.7 $\rm R_{\odot}$ respectively. WASP-150b has a mass and radius of 8.5 $\rm M_J$ and 1.1 $\rm R_J$, leading to a large planetary bulk density of 6.4 $\rm \rho_J$. WASP-150b is found to be $\sim3$ Gyr old, well below its circularisation timescale, supporting the eccentric nature of the planet. WASP-176b is a hot Jupiter planet on a 3.9 day orbit around a $V$ = 12.01, F9 sub-giant host. The host star has a mass and radius of 1.3 $\rm M_{\odot}$ and 1.9 $\rm R_{\odot}$. WASP-176b has a mass and radius of 0.86 $\rm M_J$ and 1.5 $\rm R_J$ respectively, leading to a planetary bulk density of 0.23 $\rm \rho_J$.

SPECULOOS-South, an observatory composed of four independent 1m robotic telescopes, located at ESO Paranal, Chile, started scientific operation in January 2019. This Southern Hemisphere facility operates as part of SPECULOOS, an international network of 1m-class telescopes surveying for transiting terrestrial planets around the nearest and brightest ultra-cool dwarfs. To automatically and efficiently process the observations of SPECULOOS-South, and to deal with the specialised photometric requirements of ultra-cool dwarf targets, we present our automatic pipeline. This pipeline includes an algorithm for automated differential photometry and an extensive correction technique for the effects of telluric water vapour, using ground measurements of the precipitable water vapour. Observing very red targets in the near-infrared can result in photometric systematics in the differential lightcurves, related to the temporally-varying, wavelength-dependent opacity of the Earth's atmosphere. These systematics are sufficient to affect the daily quality of the lightcurves, the longer time-scale variability study of our targets and even mimic transit-like signals. Here we present the implementation and impact of our water vapour correction method. Using the 179 nights and 98 targets observed in the I+z' filter by SPECULOOS-South since January 2019, we show the impressive photometric performance of the facility (with a median precision of ~1.5 mmag for 30-min binning of the raw, non-detrended lightcurves) and assess its detection potential. We compare simultaneous observations with SPECULOOS-South and TESS, to show that we readily achieve high-precision, space-level photometry for bright, ultra-cool dwarfs, highlighting SPECULOOS-South as the first facility of its kind.

Mass, radius, and age are three of the most fundamental parameters for celestial objects, enabling studies of the evolution and internal physics of stars, brown dwarfs, and planets. Brown dwarfs are hydrogen-rich objects that are unable to sustain core fusion reactions but are supported from collapse by electron degeneracy pressure. As they age, brown dwarfs cool, reducing their radius and luminosity. Young exoplanets follow a similar behaviour. Brown dwarf evolutionary models are relied upon to infer the masses, radii and ages of these objects. Similar models are used to infer the mass and radius of directly imaged exoplanets. Unfortunately, only sparse empirical mass, radius and age measurements are currently available, and the models remain mostly unvalidated. Double-line eclipsing binaries provide the most direct route for the absolute determination of the masses and radii of stars. Here, we report the SPECULOOS discovery of 2M1510A, a nearby, eclipsing, double-line brown dwarf binary, with a widely-separated tertiary brown dwarf companion. We also find that the system is a member of the $45\pm5$ Myr-old moving group, Argus. The system's age matches those of currently known directly-imaged exoplanets. 2M1510A provides an opportunity to benchmark evolutionary models of brown dwarfs and young planets. We find that widely-used evolutionary models do reproduce the mass, radius and age of the binary components remarkably well, but overestimate the luminosity by up to 0.65 magnitudes, which could result in underestimated photometric masses for directly-imaged exoplanets and young field brown dwarfs by 20 to 35%.

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.

We present the discovery of four new giant planets from \WASP, three hot Jupiters and one bloated sub-Saturn mass planet; WASP-169b, WASP-171b, WASP-175b and WASP-182b. Besides the discovery photometry from \wasp\ we use radial velocity measurements from CORALIE and HARPS as well as follow-up photometry from EulerCam, TRAPPIST-North and -South and SPECULOOS. WASP-169b is a low density Jupiter ($M=0.561 \pm 0.061~\mathrm{M_{Jup}}, R=1.304^{+0.150}_{-0.073} ~\mathrm{R_{Jup}}$) orbiting a V=12.17 F8 sub-giant in a 5.611~day orbit. WASP-171b is a typical hot Jupiter ($M=1.084 \pm 0.094~\mathrm{M_{Jup}}, R=0.98^{+0.07}_{-0.04} ~\mathrm{R_{Jup}}$, $P=3.82~\mathrm{days}$) around a V=13.05 G0 star. We find a linear drift in the radial velocities of WASP-171 spanning 3.5 years, indicating the possibility of an additional outer planet or stellar companion. WASP-175b is an inflated hot Jupiter ($M=0.99 \pm 0.13~\mathrm{M_{Jup}}, R=1.208 \pm 0.081 ~\mathrm{R_{Jup}}$, $P=3.07~\mathrm{days}$) around a V=12.04 F7 star, which possibly is part of a binary system with a star 7.9\arcsec\ away. WASP-182b is a bloated sub-Saturn mass planet ($M=0.148 \pm 0.011~\mathrm{M_{Jup}}, R=0.850\pm 0.030 ~\mathrm{R_{Jup}}$) around a metal rich V=11.98 G5 star ([Fe/H]$=0.27 \pm 0.11$). With a orbital period of $P=3.377~\mathrm{days}$, it sits right in the apex of the sub-Jovian desert, bordering the upper- and lower edge of the desert in both the mass-period and radius-period plane. WASP-169b, WASP-175b and WASP-182b are promising targets for atmospheric characterisation through transmission spectroscopy, with expected transmission signals of 121, 150 and 264 ppm respectively.

We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet radii ranging from 0.8REarth to 1.6REarth. All three planets have short orbital periods, ranging from 2.25 to 7.45 days with the outer pair just wide of a 2:1 period resonance. Diagnostic tests produced by the TESS Data Validation Report and the vetting package DAVE rule out common false positive sources. These analyses, along with dedicated follow-up and the multiplicity of the system, lend confidence that the observed signals are caused by planets transiting L 98-59 and are not associated with other sources in the field. The L 98-59 system is interesting for a number of reasons: the host star is bright (V = 11.7 mag, K = 7.1 mag) and the planets are prime targets for further follow-up observations including precision radial-velocity mass measurements and future transit spectroscopy with the James Webb Space Telescope; the near resonant configuration makes the system a laboratory to study planetary system dynamical evolution; and three planets of relatively similar size in the same system present an opportunity to study terrestrial planets where other variables (age, metallicity, etc.) can be held constant. L 98-59 will be observed in 4 more TESS sectors, which will provide a wealth of information on the three currently known planets and have the potential to reveal additional planets in the system.

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.

Measurements of the physical properties of stars at the lower end of the main sequence are scarce. In this context we report masses, radii and surface gravities of ten very-low-mass stars in eclipsing binary systems, with orbital periods of the order of several days. The objects probe the stellar mass-radius relation in the fully convective regime, $M_\star \lesssim 0.35$ M$_\odot$, down to the hydrogen burning mass-limit, $M_{\mathrm{HB}} \sim 0.07$ M$_\odot$. The stars were detected by the WASP survey for transiting extra-solar planets, as low-mass, eclipsing companions orbiting more massive, F- and G-type host stars. We use eclipse observations of the host stars (TRAPPIST, Leonhard Euler, SPECULOOS telescopes), and radial velocities of the host stars (CORALIE spectrograph), to determine physical properties of the low-mass companions. Companion surface gravities are derived from the eclipse and orbital parameters of each system. Spectroscopic measurements of the host star effective temperature and metallicity are used to infer the host star mass and age from stellar evolution models. Masses and radii of the low-mass companions are then derived from the eclipse and orbital parameters of each system. The objects are compared to stellar evolution models for low-mass stars, to test for an effect of the stellar metallicity and orbital period on the radius of low-mass stars in close binary systems. Measurements are in good agreement with stellar models; an inflation of the radii of low-mass stars with respect to model predictions is limited to 1.6 $\pm$ 1.2% in the fully convective regime. The sample of ten objects indicates a scaling of the radius of low-mass stars with the host star metallicity. No correlation between stellar radii and orbital periods of the binary systems is determined. A combined analysis with comparable objects from the literature is consistent with this result.

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.

We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance.

We present the discovery by the WASP-South transit survey of three new transiting hot Jupiters, WASP-161 b, WASP-163 b and WASP-170 b. Follow-up radial velocities obtained with the Euler/CORALIE spectrograph and high-precision transit light curves obtained with the TRAPPIST-North, TRAPPIST-South, SPECULOOS-South, NITES, and Euler telescopes have enabled us to determine the masses and radii for these transiting exoplanets. WASP-161\,b completes an orbit around its $V=11.1$ F6V-type host star in 5.406 days, and has a mass and radius of $2.5\pm 0.2$$M_{Jup}$ and $1.14\pm 0.06$ $R_{Jup}$ respectively. WASP-163\,b has an orbital period of 1.609 days, a mass of $1.9\pm0.2$ $M_{Jup}$, and a radius of $1.2\pm0.1$ $R_{Jup}$. Its host star is a $V=12.5$ G8-type dwarf. WASP-170\,b is on a 2.344 days orbit around a G1V-type star of magnitude $V=12.8$. It has a mass of $1.7\pm0.2$ $M_{Jup}$ and a radius of $1.14\pm0.09$ $R_{Jup}$. Given their irradiations ($\sim10^9$ erg.s$^{-1}$.cm$^{-2}$) and masses, the three new planets sizes are in good agreement with classical structure models of irradiated giant planets.

The TRAPPIST-1 planetary system represents an exceptional opportunity for the atmospheric characterization of temperate terrestrial exoplanets with the upcoming James Webb Space Telescope (JWST). Assessing the potential impact of stellar contamination on the planets' transit transmission spectra is an essential precursor step to this characterization. Planetary transits themselves can be used to scan the stellar photosphere and to constrain its heterogeneity through transit depth variations in time and wavelength. In this context, we present our analysis of 169 transits observed in the optical from space with K2 and from the ground with the SPECULOOS and Liverpool telescopes. Combining our measured transit depths with literature results gathered in the mid/near-IR with Spitzer/IRAC and HST/WFC3, we construct the broadband transmission spectra of the TRAPPIST-1 planets over the 0.8-4.5 $\mu$m spectral range. While planets b, d, and f spectra show some structures at the 200-300ppm level, the four others are globally flat. Even if we cannot discard their instrumental origins, two scenarios seem to be favored by the data: a stellar photosphere dominated by a few high-latitude giant (cold) spots, or, alternatively, by a few small and hot (3500-4000K) faculae. In both cases, the stellar contamination of the transit transmission spectra is expected to be less dramatic than predicted in recent papers. Nevertheless, based on our results, stellar contamination can still be of comparable or greater order than planetary atmospheric signals at certain wavelengths. Understanding and correcting the effects of stellar heterogeneity therefore appears essential to prepare the exploration of TRAPPIST-1's with JWST.

One of the most powerful methods used to estimate sky-projected spin-orbit angles of exoplanetary systems is through a spectroscopic transit observation known as the Rossiter-McLaughlin (RM) effect. So far mostly single RM observations have been used to estimate the spin-orbit angle, and thus there have been no studies regarding the variation of estimated spin-orbit angle from transit to transit. Stellar activity can alter the shape of photometric transit light curves and in a similar way they can deform the RM signal. In this paper we discuss several RM observations, obtained using the HARPS spectrograph, of known transiting planets that all transit extremely active stars, and by analyzing them individually we assess the variation in the estimated spin-orbit angle. Our results reveal that the estimated spin-orbit angle can vary significantly (up to 42 degrees) from transit to transit, due to variation in the configuration of stellar active regions over different nights. This finding is almost two times larger than the expected variation predicted from simulations. We could not identify any meaningful correlation between the variation of estimated spin-orbit angles and the stellar magnetic activity indicators. We also investigated two possible approaches to mitigate the stellar activity influence on RM observations. The first strategy was based on obtaining several RM observations and folding them to reduce the stellar activity noise. Our results demonstrated that this is a feasible and robust way to overcome this issue. The second approach is based on acquiring simultaneous high-precision short-cadence photometric transit light curves using TRAPPIST/SPECULOOS telescopes, which provide more information about the stellar active region's properties and allow a better RM modeling.

We report the discovery of NGTS-4b, a sub-Neptune-sized planet transiting a 13th magnitude K-dwarf in a 1.34d orbit. NGTS-4b has a mass M=$20.6\pm3.0$M_E and radius R=$3.18\pm0.26$R_E, which places it well within the so-called "Neptunian Desert". The mean density of the planet ($3.45\pm0.95$g/cm^3) is consistent with a composition of 100% H$_2$O or a rocky core with a volatile envelope. NGTS-4b is likely to suffer significant mass loss due to relatively strong EUV/X-ray irradiation. Its survival in the Neptunian desert may be due to an unusually high core mass, or it may have avoided the most intense X-ray irradiation by migrating after the initial activity of its host star had subsided. With a transit depth of $0.13\pm0.02$%, NGTS-4b represents the shallowest transiting system ever discovered from the ground, and is the smallest planet discovered in a wide-field ground-based photometric survey.

We apply the transit light curve self-contamination technique of Morris et al. (2018) to search for the effect of stellar activity on the transits of the ultracool dwarf TRAPPIST-1 with 2018 Spitzer photometry. The self-contamination method fits the transit light curves of planets orbiting spotted stars, allowing the host star to be a source of contaminating positive or negative flux which influences the transit depths but not the ingress/egress durations. We find that none of the planets show statistically significant evidence for self-contamination by bright or dark regions of the stellar photosphere. However, we show that small-scale magnetic activity, analogous in size to the smallest sunspots, could still be lurking in the transit photometry undetected.

We present here SPECULOOS, a new exoplanet transit search based on a network of 1m-class robotic telescopes targeting the $\sim$1200 ultracool (spectral type M7 and later) dwarfs bright enough in the infrared ($K$-mag $\leq 12.5$) to possibly enable the atmospheric characterization of temperate terrestrial planets with next-generation facilities like the $\textit{James Webb Space Telescope}$. The ultimate goals of the project are to reveal the frequency of temperate terrestrial planets around the lowest-mass stars and brown dwarfs, to probe the diversity of their bulk compositions, atmospheres and surface conditions, and to assess their potential habitability.

We present the discovery of NGTS-3Ab, a hot Jupiter found transiting the primary star of an unresolved binary system. We develop a joint analysis of multi-colour photometry, centroids, radial velocity (RV) cross-correlation function (CCF) profiles and their bisector inverse slopes (BIS) to disentangle this three-body system. Data from the Next Generation Transit Survey (NGTS), SPECULOOS and HARPS are analysed and modelled with our new blendfitter software. We find that the binary consists of NGTS-3A (G6V-dwarf) and NGTS-3B (K1V-dwarf) at <1 arcsec separation. NGTS-3Ab orbits every P = 1.675 days. The planet radius and mass are R_planet = 1.48+-0.37 R_J and M_planet = 2.38+-0.26 M_J, suggesting it is potentially inflated. We emphasise that only combining all the information from multi-colour photometry, centroids and RV CCF profiles can resolve systems like NGTS-3. Such systems cannot be disentangled from single-colour photometry and RV measurements alone. Importantly, the presence of a BIS correlation indicates a blend scenario, but is not sufficient to determine which star is orbited by the third body. Moreover, even if no BIS correlation is detected, a blend scenario cannot be ruled out without further information. The choice of methodology for calculating the BIS can influence the measured significance of its correlation. The presented findings are crucial to consider for wide-field transit surveys, which require wide CCD pixels (>5 arcsec) and are prone to contamination by blended objects. With TESS on the horizon, it is pivotal for the candidate vetting to incorporate all available follow-up information from multi-colour photometry and RV CCF profiles.