We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

We carry out a phase curve analysis of the KELT-9 system using photometric observations from NASA's Transiting Exoplanet Survey Satellite. The measured secondary eclipse depth and peak-to-peak atmospheric brightness modulation are $650^{+14}_{-15}$ ppm and $566\pm16$ ppm, respectively. The planet's brightness variation reaches maximum $31\pm5$ minutes before the midpoint of secondary eclipse, indicating a $5\overset{\circ}{.}2\pm0\overset{\circ}{.}9$ eastward shift in the dayside hotspot from the substellar point. We also detect stellar pulsations on KELT-9 with a period of $7.58695\pm0.00091$ hours. The dayside emission of KELT-9b in the TESS bandpass is consistent with a blackbody brightness temperature of $4600\pm100$ K. The corresponding nightside brightness temperature is $3040\pm100$ K, comparable to the dayside temperatures of the hottest known exoplanets. In addition, we detect a significant phase curve signal at the first harmonic of the orbital frequency and a marginal signal at the second harmonic. While the amplitude of the first harmonic component is consistent with the predicted ellipsoidal distortion modulation assuming equilibrium tides, the phase of this photometric variation is shifted relative to the expectation. Placing KELT-9b in the context of other exoplanets with phase curve observations, we find that the elevated nightside temperature and relatively low day-night temperature contrast agree with the predictions of atmospheric models that include H$_{2}$ dissociation and recombination. The nightside temperature of KELT-9b implies an atmospheric composition containing about 50% molecular and 50% atomic hydrogen at 0.1 bar, a nightside emission spectrum that deviates significantly from a blackbody, and a 0.5-2.0 $\mu$m transmission spectrum that is featureless at low resolution.

We present the discovery of three sub-Neptune-sized planets transiting the nearby and bright Sun-like star HD 191939 (TIC 269701147, TOI 1339), a $K_{s}=7.18$ magnitude G8 V dwarf at a distance of only 54 parsecs. We validate the planetary nature of the transit signals by combining five months of data from the Transiting Exoplanet Survey Satellite with follow-up ground-based photometry, archival optical images, radial velocities, and high angular resolution observations. The three sub-Neptunes have similar radii ($R_{b} = 3.42^{+0.11}_{-0.11}\,R_{\oplus}$, $R_{c}=3.23_{-0.11}^{+0.11}\,R_{\oplus}$, and $R_{d}=3.16_{-0.11}^{+0.11}\,R_{\oplus}$) and their orbits are consistent with a stable, circular, and co-planar architecture near mean motion resonances of 1:3 and 3:4 ($P_{b}=8.88$ days, $P_{c}=28.58$ days, and $P_{d}=38.35$ days). The HD~191939 system is an excellent candidate for precise mass determinations of the planets with high-resolution spectroscopy due to the host star's brightness and low chromospheric activity. Moreover, the system's compact and near-resonant nature can provide an independent way to measure planetary masses via transit timing variations while also enabling dynamical and evolutionary studies. Finally, as a promising target for multi-wavelength transmission spectroscopy of all three planets' atmospheres, HD 191939 can offer valuable insight into multiple sub-Neptunes born from a proto-planetary disk that may have resembled that of the early Sun.

Although the $TESS$ primary mission observed the northern and southern ecliptic hemispheres, generally avoiding the ecliptic, and the $Kepler$ space telescope during the $K2$ mission could only observe near the ecliptic, many of the $K2$ fields extend far enough from the ecliptic plane that sections overlap with $TESS$ fields. Using photometric observations from both $K2$ and $TESS$, combined with archival spectroscopic observations, we globally modeled four known planetary systems discovered by $K2$ that were observed in the first year of the primary $TESS$ mission. Specifically, we provide updated ephemerides and system parameters for K2-114 b, K2-167 b, K2-237 b, and K2-261 b. These were some of the first $K2$ planets to be observed by $TESS$ in the first year and include three Jovian sized planets and a sub-Neptune with orbital periods less than 12 days. In each case, the updated ephemeris significantly reduces the uncertainty in prediction of future times of transit, which is valuable for planning observations with the James Webb Space Telescope and other future facilities. The $TESS$ extended mission is expected to observe about half of the $K2$ fields, providing the opportunity to perform this type of analysis on a larger number of systems.

We present allesfitter, a public and open-source python software for flexible and robust inference of stars and exoplanets given photometric and radial velocity data. Allesfitter offers a rich selection of orbital and transit/eclipse models, accommodating multiple exoplanets, multi-star systems, transit-timing variations, phase curves, stellar variability, star spots, stellar flares, and various systematic noise models including Gaussian Processes. It features both parameter estimation and Bayesian model selection, allowing to easily run either a Markov Chain Monte Carlo (MCMC) or Nested Sampling fit. For novice users, a graphical user interface allows to specify all input and perform analyses; for python users, all modules can be readily imported in any existing scripts. Allesfitter also produces publication-ready tables, LaTeX commands, and figures. The software is publicly available (https://github.com/MNGuenther/allesfitter), pip-installable (pip install allesfitter) and well documented (www.allesfitter.com). Finally, we demonstrate the software's capabilities on several examples and provide updates to the literature where possible for Pi Mensae, TOI-216, WASP-18, KOI-1003 and GJ 1243.

We report the Transiting Exoplanet Survey Satellite discovery of three small planets transiting one of the nearest and brightest M dwarf hosts to date, TOI-270 (TIC 259377017; K-mag 8.3; 22.5 parsec). The M3V-type star is transited by the super-Earth-sized TOI-270 b (1.247+0.089-0.083 R_earth) and the sub-Neptune-sized exoplanets TOI-270 c (2.42+-0.13 R_earth) and TOI-270 d (2.13+-0.12 R_earth). The planets orbit close to a mean-motion resonant chain, with periods (3.36, 5.66, and 11.38 days) near ratios of small integers (5:3 and 2:1). TOI-270 is a prime target for future studies since: 1) its near-resonance allows detecting transit timing variations for precise mass measurements and dynamical studies; 2) its brightness enables independent radial velocity mass measurements; 3) the outer planets are ideal for atmospheric characterisation via transmission spectroscopy; and 4) the quiet star enables future searches for habitable zone planets. Altogether, very few systems with small, temperate exoplanets are as suitable for such complementary and detailed characterisation as TOI-270.

We perform a study of stellar flares for the 24,809 stars observed with 2 minute cadence during the first two months of the TESS mission. Flares may erode exoplanets' atmospheres and impact their habitability, but might also trigger the genesis of life around small stars. TESS provides a new sample of bright dwarf stars in our galactic neighborhood, collecting data for thousands of M-dwarfs that might host habitable exoplanets. Here, we use an automated search for flares accompanied by visual inspection. Then, our public allesfitter code robustly selects the appropriate model for potentially complex flares via Bayesian evidence. We identify 1228 flaring stars, 673 of which are M-dwarfs. Among 8695 flares in total, the largest superflare increased the stellar brightness by a factor of 16.1. Bolometric flare energies range from 10^31.0 to 10^36.9 erg, with a median of 10^33.1 erg. Furthermore, we study the flare rate and energy as a function of stellar type and rotation period. We solidify past findings that fast rotating M-dwarfs are the most likely to flare, and that their flare amplitude is independent of the rotation period. Finally, we link our results to criteria for prebiotic chemistry, atmospheric loss through coronal mass ejections, and ozone sterilization. Four of our flaring M dwarfs host exoplanet candidates alerted on by TESS, for which we discuss how these effects can impact life. With upcoming TESS data releases, our flare analysis can be expanded to almost all bright small stars, aiding in defining criteria for exoplanet habitability.

We report the detection of a Saturn-size exoplanet orbiting HD 332231 (TOI 1456) in light curves from the Transiting Exoplanet Survey Satellite (TESS). HD 332231, an F8 dwarf star with a V-band magnitude of 8.56, was observed by TESS in Sectors 14 and 15. We detect a single-transit event in the Sector 15 presearch data conditioning (PDC) light curve. We obtain spectroscopic follow-up observations of HD 332231 with the Automated Planet Finder, Keck I, and SONG telescopes. The orbital period we infer from the radial velocity (RV) observations leads to the discovery of another transit in Sector 14 that was masked by PDC due to scattered light contamination. A joint analysis of the transit and RV data confirms the planetary nature of HD 332231 b, a Saturn-size ($0.867^{+0.027}_{-0.025} \; R_{\rm J}$), sub-Saturn-mass ($0.244\pm0.021 \; M_{\rm J}$) exoplanet on a 18.71 day circular orbit. The low surface gravity of HD 332231 b and the relatively low stellar flux it receives make it a compelling target for transmission spectroscopy. Also, the stellar obliquity is likely measurable via the Rossiter-McLaughlin effect, an exciting prospect given the 0.14 au orbital separation of HD 332231 b. The spectroscopic observations do not provide substantial evidence for any additional planets in the HD 332231 system, but continued RV monitoring is needed to further characterize this system. We also predict that the frequency and duration of masked data in the PDC light curves for TESS Sectors 14-16 could hide transits of some exoplanets with orbital periods between 10.5 and 17.5 days.

We report the discovery and validation of four extrasolar planets hosted by the nearby, bright, Sun-like (G3V) star HD~108236 using data from the Transiting Exoplanet Survey Satellite (TESS). We present transit photometry, reconnaissance and precise Doppler spectroscopy as well as high-resolution imaging, to validate the planetary nature of the objects transiting HD~108236, also known as the TESS Object of Interest (TOI) 1233. The innermost planet is a possibly-rocky super-Earth with a period of $3.79523_{-0.00044}^{+0.00047}$ days and has a radius of $1.586\pm0.098$ $R_\oplus$. The outer planets are sub-Neptunes, with potential gaseous envelopes, having radii of $2.068_{-0.091}^{+0.10}$ $R_\oplus$, $2.72\pm0.11$ $R_\oplus$, and $3.12_{-0.12}^{+0.13}$ $R_\oplus$ and periods of $6.20370_{-0.00052}^{+0.00064}$ days, $14.17555_{-0.0011}^{+0.00099}$ days, and $19.5917_{-0.0020}^{+0.0022}$ days, respectively. With V and K$_{\rm s}$ magnitudes of 9.2 and 7.6, respectively, the bright host star makes the transiting planets favorable targets for mass measurements and, potentially, for atmospheric characterization via transmission spectroscopy. HD~108236 is the brightest Sun-like star in the visual (V) band known to host four or more transiting exoplanets. The discovered planets span a broad range of planetary radii and equilibrium temperatures, and share a common history of insolation from a Sun-like star ($R_\star = 0.888 \pm 0.017$ R$_\odot$, $T_{\rm eff} = 5730 \pm 50$ K), making HD 108236 an exciting, opportune cosmic laboratory for testing models of planet formation and evolution.

Probabilistic cataloging (PCAT) outperforms traditional cataloging methods on single-band optical data in crowded fields (Portillo et al. 2017). We extend our work to multiple bands, achieving greater sensitivity ($\sim$ 0.4 mag) and greater speed (500x) compared to previous single-band results. We demonstrate the effectiveness of multiband PCAT on mock data, both in terms of recovering accurate posteriors in the catalog space, and in directly deblending sources. When applied to Sloan Digital Sky Survey (SDSS) observations of M2, taking Hubble Space Telescope data as truth, our joint fit on $r$ and $i$ band data goes $\sim0.4$ mag deeper than single-band probabilistic cataloging and has a false discovery rate less than 20\% for F606W$\leq 20$. Compared to DAOPHOT, the two-band SDSS catalog fit goes nearly 1.5 magnitudes deeper using the same data, and maintains a lower false discovery rate down to F606W$\sim 20.5$. Given recent improvements in computational speed, multiband PCAT shows promise in application to large-scale surveys and is a plausible framework for joint analysis of multi-instrument observational data.

We present a systematic phase curve analysis of known transiting systems observed by the Transiting Exoplanet Survey Satellite during Year 1 of the Primary Mission. Using theoretical predictions for the amplitude of the planetary longitudinal atmospheric brightness modulation, stellar ellipsoidal distortion and Doppler boosting, as well as brightness considerations to select targets with likely detectable signals, we applied a uniform data processing and light curve modeling framework to fit the full-orbit phase curves of 19 transiting systems with planet-mass or brown dwarf companions. Statistically significant secondary eclipse depths and/or atmospheric brightness modulation amplitudes were measured for HIP 65A, WASP-72, WASP-100, WASP-111, and WASP-122/KELT-14. For WASP-100b, we found marginal evidence that the brightest region of the atmosphere is shifted eastward by $14\overset{\circ}{.}2\pm4\overset{\circ}{.}6$ away from the substellar point. We found significant ellipsoidal distortion signals in the light curves of HIP 65A, TOI-503, and WASP-30, with TOI-503 also exhibiting Doppler boosting. The measured amplitudes of these signals agree with the predictions of theoretical models. From our light curve fits, we also measured updated and refined orbital ephemerides and transit shape parameters. Combining the optical secondary eclipse depths with previously published Spitzer 3.6 and 4.5 $\mu$m measurements, we derived dayside brightness temperatures and visible-light geometric albedos for a subset of the analyzed systems. We find a tentative correlation between increasing dayside temperature and increasing geometric albedo in the TESS bandpass, suggesting enhanced atmospheric reflectivity and/or additional opacity sources contributing to the visible-light flux for planets with $T_{\mathrm{day}}>2200$ K.

We report the discovery of GJ 1252 b, a planet with a radius of 1.193 $\pm$ 0.074 $R_{\oplus}$ and an orbital period of 0.52 days around an M3-type star (0.381 $\pm$ 0.019 $M_{\odot}$, 0.391 $\pm$ 0.020 $R_{\odot}$) located 20.385 $\pm$ 0.019 pc away. We use TESS data, ground-based photometry and spectroscopy, Gaia astrometry, and high angular resolution imaging to show that the transit signal seen in the TESS data must originate from a transiting planet. We do so by ruling out all false positive scenarios that attempt to explain the transit signal as originating from an eclipsing stellar binary. Precise Doppler monitoring also leads to a tentative mass measurement of 2.09 $\pm$ 0.56 $M_{\oplus}$. The host star proximity, brightness ($V$ = 12.19 mag, $K$ = 7.92 mag), low stellar activity, and the system's short orbital period make this planet an attractive target for detailed characterization, including precise mass measurement, looking for other objects in the system, and planet atmosphere characterization.

We report the discovery of a warm sub-Saturn, TOI-257b (HD 19916b), based on data from NASA's Transiting Exoplanet Survey Satellite (TESS). The transit signal was detected by TESS and confirmed to be of planetary origin based on radial-velocity observations with the Minerva-Australis telescope array. An analysis of the TESS photometry, the Minerva-Australis, FEROS, and HARPS radial velocities, and the asteroseismic data of the stellar oscillations reveals that TOI-257b has a mass of $M_P=0.134^{+0.023}_{-0.022}$$\rm{M_J}$ ($42.6^{+7.3}_{-7.0}$$\rm{M_{\oplus}}$), a radius of $R_P=0.626^{+0.013}_{-0.012}$$\rm{R_J}$ ($7.02^{+0.15}_{-0.13}$$\rm{R_{\oplus}}$), and an orbit with eccentricity $0.242^{+0.040}_{-0.065}$ and period $18.38827\pm0.00072$$\rm{days}$. TOI-257b orbits a bright ($\mathrm{V}=7.570$mag) somewhat evolved late F-type star with $M_*=1.390\pm0.046$$\rm{M_{\odot}}$, $R_*=1.888\pm0.033$$\rm{R_{\odot}}$, $T_{\rm eff}=6075\pm90$$\rm{K}$, and $v\sin{i}=11.3\pm0.5$km/s. Additionally, we statistically validate a second non-transiting sub-Saturn mass planet on a $\sim71$day orbit using the radial velocity data. This system joins the ranks of a small number of exoplanet host stars that have been characterized with asteroseismology. Warm sub-Saturns are rare in the known sample of exoplanets, and thus the discovery of TOI-257b is important in the context of future work studying the formation and migration history of similar planetary systems.

We report the discovery of TOI-677 b, first identified as a candidate in light curves obtained within Sectors 9 and 10 of the Transiting Exoplanet Survey Satellite (TESS) mission and confirmed with radial velocities. TOI-677 b has a mass of M_p = 1.236$^{+0.069}_{-0.067}$ M_J, a radius of R_p = 1.170 +- 0.03 R_J,and orbits its bright host star (V=9.8 mag) with an orbital period of 11.23660 +- 0.00011 d, on an eccentric orbit with e = 0.435 +- 0.024. The host star has a mass of M_* = 1.181 +- 0.058 M_sun, a radius of R_* = 1.28 +- 0.03 R_sun, an age of 2.92$^{+0.80}_{-0.73}$ Gyr and solar metallicity, properties consistent with a main sequence late F star with T_eff = 6295 +- 77 K. We find evidence in the radial velocity measurements of a secondary long term signal which could be due to an outer companion. The TOI-677 b system is a well suited target for Rossiter-Mclaughlin observations that can constrain migration mechanisms of close-in giant planets.

The exoplanet HD 118203 b, orbiting a bright (V = 8.05) host star, was discovered using the radial velocity method by da Silva et al. (2006), but was not previously known to transit. TESS photometry has revealed that this planet transits its host star. Five planetary transits were observed by TESS, allowing us to measure the radius of the planet to be $1.133 \pm 0.031 R_J$, and to calculate the planet mass to be $2.173 \pm 0.078 M_J$. The host star is slightly evolved with an effective temperature of $T_{\rm eff} = 5692 \pm 83$ K and a surface gravity of ${\rm log}(g) = 3.891 \pm 0.019$. With an orbital period of $6.134980 \pm 0.000038$ days and an eccentricity of $0.316 \pm 0.021$, the planet occupies a transitional regime between circularized hot Jupiters and more dynamically active planets at longer orbital periods. The host star is among the ten brightest known to have transiting giant planets, providing opportunities for both planetary atmospheric and asteroseismic studies.

We report the discovery of the 1.008-day, ultra-short period (USP) super-Earth HD 213885b (TOI-141b) orbiting the bright ($V=7.9$) star HD 213885 (TOI-141, TIC 403224672), detected using photometry from the recently launched TESS mission. Using FEROS, HARPS and CORALIE radial-velocities, we measure a precise mass of $8.8\pm0.6$ $M_\oplus$ for this $1.74 \pm 0.05$ $R_\oplus$ exoplanet, which provides enough information to constrain its bulk composition, which is similar to Earth's but enriched in iron. The radius, mass and stellar irradiation of HD 213885b are, given our data, very similar to 55 Cancri e, making this exoplanet a good target to perform comparative exoplanetology of short period, highly irradiated super-Earths. Our precise radial-velocities reveal an additional $4.78$-day signal which we interpret as arising from a second, non-transiting planet in the system, HD 213885c (TOI-141c), whose minimum mass of $19.95\pm 1.4$ $M_\oplus$ makes it consistent with being a Neptune-mass exoplanet. The HD 213885 system is very interesting from the perspective of future atmospheric characterization, being the second brightest star to host an ultra-short period transiting super-Earth (with the brightest star being, in fact, 55 Cancri). Prospects for characterization with present and future observatories are discussed.

We study the red-optical photometry of the ultra-hot Jupiter WASP-121 b as observed by the Transiting Exoplanet Survey Satellite (TESS) and model its atmosphere through a radiative transfer simulation. Given its short orbital period of 1.275 days, inflated state and bright host star, WASP-121 b is exceptionally favorable for detailed atmospheric characterization. Towards this purpose, we use allesfitter to characterize its full red-optical phase curve, including the planetary phase modulation and the secondary eclipse. We measure the day and nightside brightness temperatures in the TESS passband as 2940+38-41 K and 2190+294-146 K, respectively, and find no phase shift between the brightest and substellar points. This is consistent with an inefficient heat recirculation on the planet. We then perform an atmospheric retrieval analysis to infer the dayside atmospheric properties of WASP-121 b such as its bulk composition, albedo and heat recirculation. We confirm the temperature inversion in the atmosphere and suggest H-, TiO and VO as potential causes of the inversion, absorbing heat at optical wavelengths at low pressures. Future HST and JWST observations of WASP-121 b will benefit from its first full phase curve measured by TESS.

We report the discovery of TOI-172 b from the Transiting Exoplanet Survey Satellite (TESS) mission, a massive hot Jupiter transiting a slightly evolved G-star with a 9.48-day orbital period. This is the first planet to be confirmed from analysis of only the TESS full frame images, because the host star was not chosen as a two minute cadence target. From a global analysis of the TESS photometry and follow-up observations carried out by the TESS Follow-up Observing Program Working Group, TOI-172 (TIC 29857954) is a slightly evolved star with an effective temperature of $T_{\rm eff}$ =$5645\pm50$ K, a mass of $M_{\star}$ = $1.128^{+0.065}_{-0.061}$ $M_{\odot}$, radius of $R_{\star}$ = $1.777^{+0.047}_{-0.044}$ $R_{\odot}$, a surface gravity of $\log$ $g_{\star}$ = $3.993^{+0.027}_{-0.028}$, and an age of $7.4^{+1.6}_{-1.5}$ Gyr. Its planetary companion (TOI-172 b) has a radius of $R_{\rm P}$ = $0.965^{+0.032}_{-0.029}$ $R_{\rm J}$, a mass of $M_{\rm P}$ = $5.42^{+0.22}_{-0.20}$ $M_{\rm J}$, and is on an eccentric orbit ($e = 0.3806^{+0.0093}_{-0.0090}$). TOI-172 b is one of the few known massive giant planets on a highly eccentric short-period orbit. Future study of the atmosphere of this planet and its system architecture offer opportunities to understand the formation and evolution of similar systems.

The Transiting Exoplanet Survey Satellite (TESS) recently observed 18 transits of the hot Jupiter WASP-4b. The sequence of transits occurred 81.6 $\pm$ 11.7 seconds earlier than had been predicted, based on data stretching back to 2007. This is unlikely to be the result of a clock error, because TESS observations of other hot Jupiters (WASP-6b, 18b, and 46b) are compatible with a constant period, ruling out an 81.6-second offset at the 6.4$\sigma$ level. The 1.3-day orbital period of WASP-4b appears to be decreasing at a rate of $\dot{P} = -12.6 \pm 1.2$ milliseconds per year. The apparent period change might be caused by tidal orbital decay or apsidal precession, although both interpretations have shortcomings. The gravitational influence of a third body is another possibility, though at present there is minimal evidence for such a body. Further observations are needed to confirm and understand the timing variation.

NASA's Transiting Exoplanet Survey Satellite (TESS) presents us with an unprecedented volume of space-based photometric observations that must be analyzed in an efficient and unbiased manner. With at least $\sim1,000,000$ new light curves generated every month from full frame images alone, automated planet candidate identification has become an attractive alternative to human vetting. Here we present a deep learning model capable of performing triage and vetting on TESS candidates. Our model is modified from an existing neural network designed to automatically classify Kepler candidates, and is the first neural network to be trained and tested on real TESS data. In triage mode, our model can distinguish transit-like signals (planet candidates and eclipsing binaries) from stellar variability and instrumental noise with an average precision (the weighted mean of precisions over all classification thresholds) of 97.0% and an accuracy of 97.4%. In vetting mode, the model is trained to identify only planet candidates with the help of newly added scientific domain knowledge, and achieves an average precision of 69.3% and an accuracy of 97.8%. We apply our model on new data from Sector 6, and present 288 new signals that received the highest scores in triage and vetting and were also identified as planet candidates by human vetters. We also provide a homogeneously classified set of TESS candidates suitable for future training.

Warm, large exoplanets with 10-100 day orbital periods pose a major challenge to our understanding of how planetary systems form and evolve. Although high eccentricity tidal migration has been invoked to explain their proximity to their host stars, a handful reside in or near orbital resonance with nearby planets, suggesting a gentler history of in situ formation or disk migration. Here we confirm and characterize a pair of warm, large exoplanets discovered by the TESS Mission orbiting K-dwarf TOI-216. Our analysis includes additional transits and transit exclusion windows observed via ground-based follow-up. We find two families of solutions, one corresponding to a sub-Saturn-mass planet accompanied by a Neptune-mass planet and the other to a Jupiter in resonance with a sub-Saturn-mass planet. We prefer the second solution based on the orbital period ratio, the planet radii, the lower free eccentricities, and libration of the 2:1 resonant argument, but cannot rule out the first. The free eccentricities and mutual inclination are compatible with stirring by other, undetected planets in the system, particularly for the second solution. We discuss prospects for better constraints on the planets' properties and orbits through follow-up, including transits observed from the ground.

We present early time light curves of 18 Type Ia supernovae observed in the first six sectors of TESS data. Nine of these supernovae were discovered by ASAS-SN, four by ATLAS, four by ZTF, and one by \textit{Gaia}. For eight of these objects with sufficient dynamic range ($>$3.0 mag from detection to peak), we fit power law models and search for signatures of companion stars. Most of our sources are consistent with fireball models where the flux increases $\propto t^2$, while two display a flatter rise with flux $\propto t$. We do not find any obvious evidence for additional structure such as multiple power law components in the early rising light curves. Assuming a favorable viewing angle, we place conservative upper limits on the radii of any companion stars of less than 25 R$_\odot$ for six supernovae and less than 4 R$_\odot$ for 4 supernovae. If such systems are commonplace, the odds of non-detection in this sample are 38\% for companions $\gtrsim$25\,R$\odot$ and 52\% for companions $\gtrsim$4\,R$\odot$. Based on the examples in this work, we expect that TESS will either detect emission from a dwarf companion during its initial two year survey, or place stringent limits on single-degenerate progenitor models.

We present a visible-light full orbital phase curve of the transiting planet WASP-18b measured by the TESS Mission. The phase curve includes the transit, secondary eclipse, and sinusoidal modulations across the orbital phase shaped by the planet's atmospheric characteristics and the star-planet gravitational interaction. We measure the beaming (Doppler boosting) and tidal ellipsoidal distortion phase modulations and show that the amplitudes of both agree with theoretical expectations. We find that the light from the planet's day-side hemisphere occulted during secondary eclipse, with a relative brightness of $341_{-18}^{+17}$ ppm, is dominated by thermal emission, leading to an upper limit on the geometric albedo in the TESS band of 0.048 (2$\sigma$). We also detect the phase modulation due to the planet's atmosphere longitudinal brightness distribution. We find that its maximum is well-aligned with the sub-stellar point, to within 2.9 deg (2$\sigma$). We do not detect light from the planet's night-side hemisphere, with an upper limit of 43 ppm (2$\sigma$), which is 13 % of the day-side brightness. The low albedo, lack of atmospheric phase shift, and inefficient heat distribution from the day to night hemispheres that we deduce from our analysis are consistent with theoretical expectations and similar findings for other strongly irradiated gas giant planets. This work demonstrates the potential of TESS data for studying full orbital phase curves of transiting systems. Finally, we complement our study by looking for transit timing variations (TTVs) in the TESS data and combined with previously published transit times, although we do not find a statistically significant TTV signal.

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.

Strong lensing is a sensitive probe of the small-scale density fluctuations in the Universe. We implement a novel approach to modeling strongly lensed systems using probabilistic cataloging, which is a transdimensional, hierarchical, and Bayesian framework to sample from a metamodel (union of models with different dimensionality) consistent with observed photon count maps. Probabilistic cataloging allows us to robustly characterize modeling covariances within and across lens models with different numbers of subhalos. Unlike traditional cataloging of subhalos, it does not require model subhalos to improve the goodness of fit above the detection threshold. Instead, it allows the exploitation of all information contained in the photon count maps, for instance, when constraining the subhalo mass function. We further show that, by not including these small subhalos in the lens model, fixed-dimensional inference methods can significantly mismodel the data. Using a simulated Hubble Space Telescope (HST) dataset, we show that the subhalo mass function can be probed even when many subhalos in the sample catalogs are individually below the detection threshold and would be absent in a traditional catalog. With the planned Wide Field Infrared Space Telescope (WFIRST), simultaneous probabilistic cataloging of dark subhalos in high-resolution, deep strong lens images has the potential to constrain the subhalo mass function at even lower masses.

The DECam Plane Survey is a five-band optical and near-infrared survey of the southern Galactic plane with the Dark Energy Camera at Cerro Tololo. The survey is designed to reach past the main-sequence turn-off at the distance of the Galactic center through a reddening E(B-V) of 1.5 mag. Typical single-exposure depths are 23.7, 22.8, 22.3, 21.9, and 21.0 mag in the grizY bands, with seeing around 1 arcsecond. The footprint covers the Galactic plane with |b| < 4 degrees, 5 degrees > l > -120 degrees. The survey pipeline simultaneously solves for the positions and fluxes of tens of thousands of sources in each image, delivering positions and fluxes of roughly two billion stars with better than 10 mmag precision. Most of these objects are highly reddened and deep in the Galactic disk, probing the structure and properties of the Milky Way and its interstellar medium. The full survey is publicly available.

Cataloging is challenging in crowded fields because sources are extremely covariant with their neighbors and blending makes even the number of sources ambiguous. We present the first optical probabilistic catalog, cataloging a crowded (~0.1 sources per pixel brighter than 22nd magnitude in F606W) Sloan Digital Sky Survey r band image from M2. Probabilistic cataloging returns an ensemble of catalogs inferred from the image and thus can capture source-source covariance and deblending ambiguities. By comparing to a traditional catalog of the same image and a Hubble Space Telescope catalog of the same region, we show that our catalog ensemble better recovers sources from the image. It goes more than a magnitude deeper than the traditional catalog while having a lower false discovery rate brighter than 20th magnitude. We also present an algorithm for reducing this catalog ensemble to a condensed catalog that is similar to a traditional catalog, except it explicitly marginalizes over source-source covariances and nuisance parameters. We show that this condensed catalog has a similar completeness and false discovery rate to the catalog ensemble. Future telescopes will be more sensitive, and thus more of their images will be crowded. Probabilistic cataloging performs better than existing software in crowded fields and so should be considered when creating photometric pipelines in the Large Synoptic Space Telescope era.

Detection of point sources in images is a fundamental operation in astrophysics, and is crucial for constraining population models of the underlying point sources or characterizing the background emission. Standard techniques fall short in the crowded-field limit, losing sensitivity to faint sources and failing to track their covariance with close neighbors. We construct a Bayesian framework to perform inference of faint or overlapping point sources. The method involves probabilistic cataloging, where samples are taken from the posterior probability distribution of catalogs consistent with an observed photon count map. In order to validate our method we sample random catalogs of the gamma-ray sky in the direction of the North Galactic Pole (NGP) by binning the data in energy and Point Spread Function (PSF) classes. Using three energy bins spanning $0.3 - 1$, $1 - 3$ and $3 - 10$ GeV, we identify $270\substack{+30 \\ -10}$ point sources inside a $40^\circ \times 40^\circ$ region around the NGP above our point-source inclusion limit of $3 \times 10^{-11}$/cm$^2$/s/sr/GeV at the $1-3$ GeV energy bin. Modeling the flux distribution as a power law, we infer the slope to be $-1.92\substack{+0.07 \\ -0.05}$ and estimate the contribution of point sources to the total emission as $18\substack{+2 \\ -2}$\%. These uncertainties in the flux distribution are fully marginalized over the number as well as the spatial and spectral properties of the unresolved point sources. This marginalization allows a robust test of whether the apparently isotropic emission in an image is due to unresolved point sources or of truly diffuse origin.

Past studies have identified a spatially extended excess of $\sim$1-3 GeV gamma rays from the region surrounding the Galactic Center, consistent with the emission expected from annihilating dark matter. We revisit and scrutinize this signal with the intention of further constraining its characteristics and origin. By applying cuts to the \textit{Fermi} event parameter CTBCORE, we suppress the tails of the point spread function and generate high resolution gamma-ray maps, enabling us to more easily separate the various gamma-ray components. Within these maps, we find the GeV excess to be robust and highly statistically significant, with a spectrum, angular distribution, and overall normalization that is in good agreement with that predicted by simple annihilating dark matter models. For example, the signal is very well fit by a 36-51 GeV dark matter particle annihilating to $b\bar{b}$ with an annihilation cross section of $\sigma v = (1-3)\times 10^{-26}$ cm$^3$/s (normalized to a local dark matter density of 0.4 GeV/cm$^3$). Furthermore, we confirm that the angular distribution of the excess is approximately spherically symmetric and centered around the dynamical center of the Milky Way (within $\sim$$0.05^{\circ}$ of Sgr A$^*$), showing no sign of elongation along the Galactic Plane. The signal is observed to extend to at least $\simeq10^{\circ}$ from the Galactic Center, disfavoring the possibility that this emission originates from millisecond pulsars.