{We investigate the atmospheric composition of the long period ($P_{\rm orb}=$~10 days), eccentric exo-Saturn WASP-117b. WASP-117b could be in atmospheric temperature and chemistry similar to WASP-107b. In mass and radius WASP-117b is similar to WASP-39b, which allows a comparative study of these planets.} % methods heading (mandatory) { We analyze a near-infrared transmission spectrum of WASP-117b taken with Hubble Space Telescope/WFC3 G141, which was reduced with two independent pipelines to ensure a robust detection of the $3 \sigma$-water spectrum. High resolution measurements were taken with VLT/ESPRESSO in the optical.} % results heading (mandatory) {Using a 1D atmosphere model with isothermal temperature, uniform cloud deck and equilibrium chemistry, the Bayesian evidence of a retrieval analysis of the transmission spectrum indicates a preference for a high atmospheric metallicity ${\rm [Fe/H]}=2.58^{+0.26}_{-0.37}$ and clear skies. The data are also consistent with a lower-metallicity composition ${\rm [Fe/H]}<1.75$ and a cloud deck between $10^{-2.2} - 10^{-5.1}$~bar, but with weaker Bayesian preference. We retrieve a low CH$_4$ abundance of $<10^{-4}$ volume fraction within $1 \sigma$. We cannot constrain the temperature between theoretically imposed limits of 700 and 1000~K. Further observations are needed to confirm quenching of CH$_4$ with $K_{zz}\geq 10^8$~cm$^2$/s. We report indications of Na and K in the VLT/ESPRESSO high resolution spectrum with substantial Bayesian evidence in combination with HST data.}

We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn-mass exoplanet KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 $\mu$m, in addition to a Spitzer 4.5 $\mu$m secondary eclipse. The precise HST transmission spectrum notably reveals a low-amplitude water feature with an unusual shape. We apply a suite of modeling tools to the transmission spectrum to investigate the planet's properties. Based on a retrieval analysis with varying molecular abundances, we find strong evidence for water absorption in the spectrum and tentative evidence for other absorbers (HCN, TiO, and AlO) depending on model assumptions. The retrieved water abundance is generally $\lesssim 0.1\times$ solar (0.001--0.7$\times$ solar over a range of model assumptions), which is several orders of magnitude lower than expected from planet formation models based on the solar system metallicity trend. We also consider chemical equilibrium and self-consistent 1D radiative-convective equilibrium model fits and find that they too prefer low metallicities ($[M/H] \lesssim -2$, consistent with the free retrieval results); however, the equilibrium models fit poorly and should be interpreted with caution. Finally, we find that the dayside flux measured from the Spitzer secondary eclipse is indicative of full heat redistribution from KELT-11b's dayside to nightside, assuming the dayside is not cloudy. Altogether, the unusual results for KELT-11b's composition are suggestive of new challenges on the horizon for atmosphere and formation models in the face of increasingly precise measurements of exoplanet spectra.

In this work we show that modern data-driven machine learning techniques can be successfully applied on lunar surface remote sensing data to learn, in an unsupervised way, sufficiently good representations of the data distribution to enable lunar technosignature and anomaly detection. In particular we train an unsupervised distribution learning neural network model to find the Apollo 15 landing module in a testing dataset, with no dataset specific model or hyperparameter tuning. Sufficiently good unsupervised data density estimation has the promise of enabling myriad useful downstream tasks, including locating lunar resources for future space flight and colonization, finding new impact craters or lunar surface reshaping, and algorithmically deciding the importance of unlabeled samples to send back from power- and bandwidth-constrained missions. We show in this work that such unsupervised learning can be successfully done in the lunar remote sensing and space science contexts.

Exoplanet science is one of the most thriving fields of modern astrophysics. A major goal is the atmospheric characterization of dozens of small, terrestrial exoplanets in order to search for signatures in their atmospheres that indicate biological activity, assess their ability to provide conditions for life as we know it, and investigate their expected atmospheric diversity. None of the currently adopted projects or missions, from ground or in space, can address these goals. In this White Paper we argue that a large space-based mission designed to detect and investigate thermal emission spectra of terrestrial exoplanets in the MIR wavelength range provides unique scientific potential to address these goals and surpasses the capabilities of other approaches. While NASA might be focusing on large missions that aim to detect terrestrial planets in reflected light, ESA has the opportunity to take leadership and spearhead the development of a large MIR exoplanet mission within the scope of the "Voyage 2050" long-term plan establishing Europe at the forefront of exoplanet science for decades to come. Given the ambitious science goals of such a mission, additional international partners might be interested in participating and contributing to a roadmap that, in the long run, leads to a successful implementation. A new, dedicated development program funded by ESA to help reduce development and implementation cost and further push some of the required key technologies would be a first important step in this direction. Ultimately, a large MIR exoplanet imaging mission will be needed to help answer one of mankind's most fundamental questions: "How unique is our Earth?"

Machine learning is now used in many areas of astrophysics, from detecting exoplanets in Kepler transit signals to removing telescope systematics. Recent work demonstrated the potential of using machine learning algorithms for atmospheric retrieval by implementing a random forest to perform retrievals in seconds that are consistent with the traditional, computationally-expensive nested-sampling retrieval method. We expand upon their approach by presenting a new machine learning model, \texttt{plan-net}, based on an ensemble of Bayesian neural networks that yields more accurate inferences than the random forest for the same data set of synthetic transmission spectra. We demonstrate that an ensemble provides greater accuracy and more robust uncertainties than a single model. In addition to being the first to use Bayesian neural networks for atmospheric retrieval, we also introduce a new loss function for Bayesian neural networks that learns correlations between the model outputs. Importantly, we show that designing machine learning models to explicitly incorporate domain-specific knowledge both improves performance and provides additional insight by inferring the covariance of the retrieved atmospheric parameters. We apply \texttt{plan-net} to the Hubble Space Telescope Wide Field Camera 3 transmission spectrum for WASP-12b and retrieve an isothermal temperature and water abundance consistent with the literature. We highlight that our method is flexible and can be expanded to higher-resolution spectra and a larger number of atmospheric parameters.

This paper outlines some of the possible advancements for the technosignatures searches using the new methods currently rapidly developing in computer science, such as machine learning and deep learning. It also showcases a couple of case studies of large research programs where such methods have been already successfully implemented with notable results. We consider that the availability of data from all sky, all the time observations paired with the latest developments in computational capabilities and algorithms currently used in artificial intelligence, including automation, will spur an unprecedented development of the technosignatures search efforts.

Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$.

Accurately and rapidly classifying exoplanet candidates from transit surveys is a goal of growing importance as the data rates from space-based survey missions increases. This is especially true for NASA's TESS mission which generates thousands of new candidates each month. Here we created the first deep learning model capable of classifying TESS planet candidates. We adapted the neural network model of Ansdell et al. (2018) to TESS data. We then trained and tested this updated model on 4 sectors of high-fidelity, pixel-level simulations data created using the Lilith simulator and processed using the full TESS SPOC pipeline. We find our model performs very well on our simulated data, with 97% average precision and 92% accuracy on planets in the 2-class model. This accuracy is also boosted by another ~4% if planets found at the wrong periods are included. We also performed 3- and 4-class classification of planets, blended & target eclipsing binaries, and non-astrophysical false positives, which have slightly lower average precision and planet accuracies, but are useful for follow-up decisions. When applied to real TESS data, 61% of TCEs coincident with currently published TOIs are recovered as planets, 4% more are suggested to be EBs, and we propose a further 200 TCEs as planet candidates.

Despite the existence of co-orbital bodies in the solar system, and the prediction of the formation of co-orbital planets by planetary system formation models, no co-orbital exoplanets (also called trojans) have been detected thus far. Here we study the signature of co-orbital exoplanets in transit surveys when two planet candidates in the system orbit the star with similar periods. Such pair of candidates could be discarded as false positives because they are not Hill-stable. However, horseshoe or long libration period tadpole co-orbital configurations can explain such period similarity. This degeneracy can be solved by considering the Transit Timing Variations (TTVs) of each planet. We then focus on the three planet candidates system TOI-178: the two outer candidates of that system have similar orbital period and had an angular separation near $\pi/3$ during the TESS observation of sector 2. Based on the announced orbits, the long-term stability of the system requires the two close-period planets to be co-orbitals. Our independent detrending and transit search recover and slightly favour the three orbits close to a 3:2:2 resonant chain found by the TESS pipeline, although we cannot exclude an alias that would put the system close to a 4:3:2 configuration. We then analyse in more detail the co-orbital scenario. We show that despite the influence of an inner planet just outside the 2:3 mean-motion resonance, this potential co-orbital system can be stable on the Giga-year time-scale for a variety of planetary masses, either on a trojan or a horseshoe orbit. We predict that large TTVs should arise in such configuration with a period of several hundred days. We then show how the mass of each planet can be retrieved from these TTVs.

Over the past decade, the study of extrasolar planets has evolved rapidly from plain detection and identification to comprehensive categorization and characterization of exoplanet systems and their atmospheres. Atmospheric retrieval, the inverse modeling technique used to determine an exoplanetary atmosphere's temperature structure and composition from an observed spectrum, is both time-consuming and compute-intensive, requiring complex algorithms that compare thousands to millions of atmospheric models to the observational data to find the most probable values and associated uncertainties for each model parameter. For rocky, terrestrial planets, the retrieved atmospheric composition can give insight into the surface fluxes of gaseous species necessary to maintain the stability of that atmosphere, which may in turn provide insight into the geological and/or biological processes active on the planet. These atmospheres contain many molecules, some of them biosignatures, spectral fingerprints indicative of biological activity, which will become observable with the next generation of telescopes. Runtimes of traditional retrieval models scale with the number of model parameters, so as more molecular species are considered, runtimes can become prohibitively long. Recent advances in machine learning (ML) and computer vision offer new ways to reduce the time to perform a retrieval by orders of magnitude, given a sufficient data set to train with. Here we present an ML-based retrieval framework called Intelligent exoplaNet Atmospheric RetrievAl (INARA) that consists of a Bayesian deep learning model for retrieval and a data set of 3,000,000 synthetic rocky exoplanetary spectra generated using the NASA Planetary Spectrum Generator. Our work represents the first ML retrieval model for rocky, terrestrial exoplanets and the first synthetic data set of terrestrial spectra generated at this scale.

Space-based missions such as Kepler, and soon TESS, provide large datasets that must be analyzed efficiently and systematically. Recent work by Shallue & Vanderburg (2018) successfully used state-of-the-art deep learning models to automatically classify Kepler transit signals as either exoplanets or false positives; our application of their model yielded 95.8% accuracy and 95.5% average precision. Here we expand upon that work by including additional scientific domain knowledge into the network architecture and input representations to significantly increase overall model performance to 97.5% accuracy and 98.0% average precision. Notably, we achieve 15-20% gains in recall for the lowest signal-to-noise transits that can correspond to rocky planets in the habitable zone. We input into the network centroid time-series information derived from Kepler data plus key stellar parameters taken from the Kepler DR25 catalogue. We also implement data augmentation techniques to alleviate model over-fitting. These improvements allow us to drastically reduce the size of the model, while still maintaining improved performance; smaller models are better for generalization, for example from Kepler to TESS data. This work illustrates the importance of including expert domain knowledge in even state-of-the-art deep learning models when applying them to scientific research problems that seek to identify weak signals in noisy data. This classification tool will be especially useful for upcoming space-based photometry missions focused on finding small planets, such as TESS and PLATO.

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

Exoplanet hunting efforts have revealed the prevalence of exotic worlds with diverse properties, including Earth-sized bodies, which has fueled our endeavor to search for life beyond the Solar System. Accumulating experiences in astrophysical, chemical, and climatological characterization of uninhabitable planets are paving the way to characterization of potentially habitable planets. In this paper, we review our possibilities and limitations in characterizing temperate terrestrial planets with future observational capabilities through 2030s and beyond, as a basis of a broad range of discussions on how to advance "astrobiology" with exoplanets. We discuss the observability of not only the proposed biosignature candidates themselves, but also of more general planetary properties that provide circumstantial evidence, since the evaluation of any biosignature candidate relies on their context. Characterization of temperate Earth-size planets in the coming years will focus on those around nearby late-type stars. JWST and later 30 meter-class ground-based telescopes will empower their chemical investigations. Spectroscopic studies of potentially habitable planets around solar-type stars will likely require a designated spacecraft mission for direct imaging, leveraging technologies that are already being developed and tested as part of the WFIRST mission. Successful initial characterization of a few nearby targets will be an important touchstone toward a more detailed scrutiny and a larger survey that are envisioned beyond 2030. The broad outlook this paper presents may help develop new observational techniques to detect relevant features as well as frameworks to diagnose planets based on the observables.

Transmission spectroscopy during planetary transits is expected to be a major source of information on the atmospheres of small (approximately Earth-sized) exoplanets in the next two decades. This technique, however, is intrinsically affected by stellar spectral contamination caused by the fact that stellar photo- and chromospheres are not perfectly homogeneous. Such stellar contamination will often reach or exceed the signal introduced by the planetary spectral features. Finding effective methods to correct stellar contamination -- or at least to quantify its possible range -- for the most important exoplanets is a necessary step for our understanding of exoplanet atmospheres. This will require significantly deepening our understanding of stellar heterogeneity, which is currently limited by the available data.

We combine and calibrate publicly available data for Boyajian's star including photometry from ASAS (SN, V, I), Kepler, Gaia, SuperWASP, and citizen scientist observations (AAVSO, HAO and Burke-Gaffney). Precise (mmag) photometry covers the years 2006-2017. We show that the year-long flux variations with an amplitude of ~4% can not be explained with cyclical symmetric or asymmetric models with periods shorter than ten years. If the dips are transits, their period must exceed ten years, or their structure must evolve significantly during each 4-year long cycle.

For the first time in human history, we will soon be able to apply the scientific method to the question "Are We Alone?" The rapid advance of exoplanet discovery, planetary systems science, and telescope technology will soon allow scientists to search for life beyond our Solar System through direct observation of extrasolar planets. This endeavor will occur alongside searches for habitable environments and signs of life within our Solar System. While the searches are thematically related and will inform each other, they will require separate observational techniques. The search for life on exoplanets holds potential through the great diversity of worlds to be explored beyond our Solar System. However, there are also unique challenges related to the relatively limited data this search will obtain on any individual world. This white paper reviews the scientific community's ability to use data from future telescopes to search for life on exoplanets. This material summarizes products from the Exoplanet Biosignatures Workshop Without Walls (EBWWW). The EBWWW was constituted by a series of online and in person activities, with participation from the international exoplanet and astrobiology communities, to assess state of the science and future research needs for the remote detection of life on planets outside our Solar System.

The fields of astronomy and astrophysics are currently engaged in an unprecedented era of discovery as recent missions have revealed thousands of exoplanets orbiting other stars. While the Kepler Space Telescope mission has enabled most of these exoplanets to be detected by identifying transiting events, exoplanets often exhibit additional photometric effects that can be used to improve the characterization of exoplanets. The EXONEST Exoplanetary Explorer is a Bayesian exoplanet inference engine based on nested sampling and originally designed to analyze archived Kepler Space Telescope and CoRoT (Convection Rotation et Transits plan\'etaires) exoplanet mission data. We discuss the EXONEST software package and describe how it accommodates plug-and-play models of exoplanet-associated photometric effects for the purpose of exoplanet detection, characterization and scientific hypothesis testing. The current suite of models allows for both circular and eccentric orbits in conjunction with photometric effects, such as the primary transit and secondary eclipse, reflected light, thermal emissions, ellipsoidal variations, Doppler beaming and superrotation. We discuss our new efforts to expand the capabilities of the software to include more subtle photometric effects involving reflected and refracted light. We discuss the EXONEST inference engine design and introduce our plans to port the current MATLAB-based EXONEST software package over to the next generation Exoplanetary Explorer, which will be a Python-based open source project with the capability to employ third-party plug-and-play models of exoplanet-related photometric effects.

We present near infrared high-precision photometry for eight transiting hot Jupiters observed during their predicted secondary eclipses. Our observations were carried out using the staring mode of the WIRCam instrument on the Canada-France-Hawaii Telescope (CFHT). We present the observing strategies and data reduction methods which delivered time series photometry with statistical photometric precisionas low as 0.11%. We performed a Bayesian analysis to model the eclipse parameters and systematics simultaneously. The measured planet-to-star flux ratios allowed us to constrain the thermal emission from the day side of these hot Jupiters, as we derived the planet brightness temperatures. Our results combined with previously observed eclipses reveal an excess in the brightness temperatures relative to the blackbody prediction for the equilibrium temperatures of the planets for a wide range of heat redistribution factors. We find a trend that this excess appears to be larger for planets with lower equilibrium temperatures. This may imply some additional sources of radiation, such as reflected light from the host star and/or thermal emission from residual internal heat from the formation of the planet.

When the Stratospheric Observatory for Infrared Astronomy (SOFIA) was conceived and its first science cases defined, exoplanets had not been detected. Later studies, however, showed that optical and near-infrared photometric and spectrophotometric follow-up observations during planetary transits and eclipses are feasible with SOFIA's instrumentation, in particular with the HIPO-FLITECAM and FPI+ optical and near infrared (NIR) instruments. Additionally, the airborne-based platform SOFIA has a number of unique advantages when compared to other ground- and space-based observatories in this field of research. Here we will outline these theoretical advantages, present some sample science cases and the results of two observations from SOFIA's first five observation cycles -- an observation of the Hot Jupiter HD 189733b with HIPO and an observation of the Super-Earth GJ 1214b with FLIPO and FPI+. Based on these early products available to this science case, we evaluate SOFIA's potential and future perspectives in the field of optical and infrared exoplanet spectrophotometry in the stratosphere.

Secondary eclipses and phase curves reveal information about the reflectivity and heat distribution in exoplanet atmospheres. The phase curve is composed of a combination of reflected, and thermally emitted light from the planet, and for circular orbits the phase curve peaks during the secondary eclipse or at an orbital phase of 0.5. Physical mechanisms have been discovered which shift the phase curve maximum of tidally locked close in planets to the right, or left, of the secondary eclipse. These mechanisms include cloud formations, and atmospheric superrotation, both of which serve to shift the thermally bright hot-spot, or highly reflective bright spot, of the atmosphere away from the sub-stellar point. Here we present a methodology for optimizing observing strategies for both secondary eclipses and phase curves with the goal of maximizing the information gained about the planetary atmosphere while minimizing the (assumed) continuous observation time. We show that we can increase the duty cycle of observations aimed at the measurements of phase curve characteristics (Amplitude, Phase offset) by up to $50\%$ for future platforms such as CHEOPS and JWST. We apply this methodology to the test cases of the Spitzer phase curve of 55-Cancri-e, which displays an eastward shift in its phase curve maximum, as well as model-generated observations of an ultra-short period planet observed with CHEOPS.

The benchmark exoplanet GJ 1214b is one of the best studied transiting planets in the transition zone between rocky Earth-sized planets and gas or ice giants. This class of super-Earth/mini-Neptune planets is unknown in our Solar System, yet is one of the most frequently detected classes of exoplanets. Understanding the transition from rocky to gaseous planets is a crucial step in the exploration of extrasolar planetary systems, in particular with regard to the potential habitability of this class of planets. GJ 1214b has already been studied in detail from various platforms at many different wavelengths. Our airborne observations with SOFIA add information in the Paschen-alpha cont. 1.9 micron infrared wavelength band, which is not accessible by any other current ground- or space-based instrument due to telluric absorption or limited spectral coverage. We used FLIPO and FPI+ on SOFIA to comprehensively analyse the transmission signal of the possible water-world GJ 1214b through photometric observations during transit in three optical and one infrared channels. We present four simultaneous light curves and corresponding transit depths in three optical and one infrared channel, which we compare to previous observations and state-of-the-art synthetic atmospheric models of GJ 1214b. The final precision in transit depth is between 1.5 and 2.5 times the theoretical photon noise limit, not sensitive enough to constrain the theoretical models any better than previous observations. This is the first exoplanet observation with SOFIA that uses its full set of instruments available to exoplanet spectrophotometry. Therefore we use these results to evaluate SOFIA's potential in this field and suggest future improvements.

The F3 main sequence star KIC 8462852 (Boyajian's Star) showed deep (up to 20%) day-long brightness dips of unknown cause during the 4 years of the Kepler mission. A 0.164 mag (16%) dimming between 1890 and 1990 was claimed, based on the analysis of photographic plates from the Harvard Observatory. We have gathered an independent set of historic plates from Sonneberg Observatory, Germany, covering the years 1934-1995. With 861 magnitudes in B, and 397 magnitudes in V, we find the star to be of constant brightness within 0.03 mag per century (3%). Consistent outcomes are found using by-eye estimates of the best 119 plates. Results are supported by data from Sternberg Observatory, Moscow, which show the star as constant between 1895 and 1995. The previously claimed century-long dimming is inconsistent with our results at the 5-sigma level, however the recently reported modest dimming of 3% in the Kepler data is not inconsistent with our data. We find no periodicities or shorter trends within our limits of 5% per 5-year bin, but note a possible dimming event on 24 Oct 1978.

We announce the discovery of KELT-16b, a highly irradiated, ultra-short period hot Jupiter transiting the relatively bright ($V = 11.7$) star TYC 2688-1839-1. A global analysis of the system shows KELT-16 to be an F7V star with $T_\textrm{eff} = 6236\pm54$ K, $\log{g_\star} = 4.253_{-0.036}^{+0.031}$, [Fe/H] = -0.002$_{-0.085}^{+0.086}$, $M_\star = 1.211_{-0.046}^{+0.043} M_\odot$, and $R_\star = 1.360_{-0.053}^{+0.064} R_\odot$. The planet is a relatively high mass inflated gas giant with $M_\textrm{P} = 2.75_{-0.15}^{+0.16} M_\textrm{J}$, $R_\textrm{P} = 1.415_{-0.067}^{+0.084} R_\textrm{J}$, density $\rho_\textrm{P} = 1.20\pm0.18$ g cm$^{-3}$, surface gravity $\log{g_\textrm{P}} = 3.530_{-0.049}^{+0.042}$, and $T_\textrm{eq} = 2453_{-47}^{+55}$ K. The best-fitting linear ephemeris is $T_\textrm{C} = 2457247.24791\pm0.00019$ BJD$_{tdb}$ and $P = 0.9689951 \pm 0.0000024$ d. KELT-16b joins WASP-18b, -19b, -43b, -103b, and HATS-18b as the only giant transiting planets with $P < 1$ day. Its ultra-short period and high irradiation make it a benchmark target for atmospheric studies by HST, Spitzer, and eventually JWST. For example, as a hotter, higher mass analog of WASP-43b, KELT-16b may feature an atmospheric temperature-pressure inversion and day-to-night temperature swing extreme enough for TiO to rain out at the terminator. KELT-16b could also join WASP-43b in extending tests of the observed mass-metallicity relation of the Solar System gas giants to higher masses. KELT-16b currently orbits at a mere $\sim$ 1.7 Roche radii from its host star, and could be tidally disrupted in as little as a few $\times 10^{5}$ years (for a stellar tidal quality factor of $Q_*' = 10^5$). Finally, the likely existence of a widely separated bound stellar companion in the KELT-16 system makes it possible that Kozai-Lidov oscillations played a role in driving KELT-16b inward to its current precarious orbit.

Distance estimates from Gaia parallax and expected luminosities are compared for KIC 8462852. Gaia DR1 yields a parallax of $2.55\pm0.31$mas, that is a distance of $391.4\substack{+53.6 \\ -42.0}$pc, or $391.4\substack{+122.1 \\ -75.2}$pc including systematic uncertainty. The distance estimate based on the absolute magnitude of an F3V star and measured reddening is $\sim454\pm35$pc. Both estimates agree within $<1\sigma$, which only excludes some of the most extreme theorized scenarios for KIC 8462852. Future Gaia data releases will determine the distance to within 1% and thus allow for the determination a precise absolute luminosity.

We present new ways to identify single and multiple moons around extrasolar planets using planetary transit timing variations (TTVs) and transit duration variations (TDVs). For planets with one moon, measurements from successive transits exhibit a hitherto undescribed pattern in the TTV-TDV diagram, originating from the stroboscopic sampling of the planet's orbit around the planet-moon barycenter. This pattern is fully determined and analytically predictable after three consecutive transits. The more measurements become available, the more the TTV-TDV diagram approaches an ellipse. For planets with multiple moons in orbital mean motion resonance (MMR), like the Galilean moons, the pattern is much more complex and addressed numerically in this report. Exomoons in MMR can also form closed, predictable TTV-TDV figures if the drift of the moons' pericenters is sufficiently slow. We find that MMR exomoons produce loops in the TTV-TDV diagram and that the number of these loops is equal to the order of the MMR, or the largest integer in the MMR ratio. We use a Bayesian model and Monte Carlo simulations to test the discoverability of exomoons using TTV-TDV diagrams with current and near-future technology. In a blind test, two of us (BP, DA) successfully retrieved a large moon from simulated TTV-TDV by co-authors MH and RH, which resembled data from a known Kepler planet candidate. Single exomoons with a 10% moon-to-planet mass ratio, like to Pluto-Charon binary, can be detectable in the archival data of the Kepler primary mission. Multi-exomoon systems, however, require either larger telescopes or brighter target stars. Complementary detection methods invoking a moon's own photometric transit or its orbital sampling effect can be used for validation or falsification. A combination of TESS, CHEOPS, and PLATO data would offer a compelling opportunity for exomoon discoveries around bright stars.

The century-long photometric record of the DASCH project provides a unique window into the variability of stars normally considered to be photometrically inactive. In this paper, we look for long-term trends in the brightness of F stars, with particular attention to KIC 8462852,an F3 main sequence star that has been identified as significant short-term variability according to Kepler observations. Although a simple search for variability suggests long-term dimming of a number of F stars, we find that such trends are artifacts of the 'Menzel Gap' in the DASCH data. That includes the behavior of KIC 8462852, which we believe is consistent with constant flux over the full duration of observations. We do, however, present a selection of F stars that do have significant photometric trends, even after systematics are taken into account.

We present a statistical analysis of the accuracy of the digitized magnitudes of photometric plates on the time scale of decades. In our examination of archival Johnson B photometry from the Harvard DASCH archive, we find a median RMS scatter of lightcurves of order 0.15mag over the range B~9-17 for all calibrations. Slight underlying systematics (trends or flux discontinuities) are on a level of ~<=0.2mag per century (1889-1990) for the majority of constant stars. These historic data can be unambiguously used for processes that happen on scales of magnitudes, and need to be carefully examined in cases approaching the noise floor. The characterization of these limits in photometric stability may guide future studies in their use of plate archives. We explain these limitations for the example case of KIC8462852, which has been claimed to dim by 0.16mag per century, and show that this trend cannot be considered as significant.

The James Webb Space Telescope will revolutionize transiting exoplanet atmospheric science due to its capability for continuous, long-duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful future transiting exoplanet characterization programs. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed "community targets") that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations.(Abridged)

Presented here is an independent re-analysis of the Kepler light curve of Kepler-91 (KIC 8219268). Using the EXONEST software package, which provides both Bayesian parameter estimation and Bayesian model testing, we were able to re-confirm the planetary nature of Kepler-91b. In addition to the primary and secondary eclipses of Kepler-91b, a third dimming event appears to occur approximately $60^o$ away (in phase) from the secondary eclipse, leading to the hypothesis that a Trojan planet may be located at the L4 or L5 Lagrange points. Here, we present a comprehensive investigation of four possibilities to explain the observed dimming event using all available photometric data from the Kepler Space Telescope, recently obtained radial velocity measurements, and N-body simulations. We find that the photometric model describing Kepler-91b and a Trojan planet is highly favored over the model involving Kepler-91b alone. However, it predicts an unphysically high temperature for the Trojan companion, leading to the conclusion that the extra dimming event is likely a false-postive.

Planets emit thermal radiation and reflect incident light that they recieve from their host stars. As a planet orbits it's host star the photometric variations associated with these two effects produce very similar phase curves. If observed through only a single bandpass this leads to a degeneracy between certain planetary parameters that hinder the precise characterization of such planets. However, observing the same planet through two different bandpasses gives one much more information about the planet. Here, we develop a Bayesian methodology for combining photometry from both \emph{Kepler} and the Transiting Exoplanet Survey Satellite (TESS). In addition, we demonstrate via simulations that one can disentangle the reflected and thermally emitted light from the atmosphere of a hot-Jupiter as well as more precisely constrain both the geometric albedo and dayside temperature of the planet. This methodology can further be employed using various combinations of photometry from the James Webb Space Telescope (JWST), the Characterizing ExOplanet Satellite (CHEOPS), or the PLATO mission.

We present a comprehensive study of phase curves and secondary eclipses in the Kepler data set using all data from 16 quarters that were available in 2013-2014. Our sample consists of 20 confirmed planets with R_p > 4 R_e ,P < 10d, V_mag < 15. Here we derive their temperatures and albedos, with an eye towards constraining models for the formation and evolution of such planets. Where there was overlap our results confirm parameters derived by previous studies, whereas we present new results for Kepler 1b-8b, 12b-15b, 17b, 40b, 41b, 43b, 44b, 76b, 77b, and 412b derived in a consistent manner. We also present lightcurve analyses for Kepler 91b and Kepler 74b, which both show extra dimmings at times other than from the expected primary and secondary eclipses. Corrected for thermal emission we find most of the massive planets from our sample to be low in albedo (<0.1) with a few having higher albedo (>0.1).

Trojans are small bodies in planetary Lagrangian points. In our solar system, Jupiter has the largest number of such companions. Their existence is assumed for exoplanetary systems as well, but none has been found so far. We present an analysis by super-stacking $\sim4\times10^4$ Kepler planets with a total of $\sim9\times10^5$ transits, searching for an average trojan transit dip. Our result gives an upper limit to the average Trojan transiting area (per planet) corresponding to one body of radius $<460$km at $2\sigma$ confidence. We find a significant Trojan-like signal in a sub-sample for planets with more (or larger) Trojans for periods $>$60 days. Our tentative results can and should be checked with improved data from future missions like PLATO2.0, and can guide planetary formation theories.

Time-series transit photometry from the Kepler space telescope has allowed for the discovery of thousands of exoplanets. We explore the potential of yet improved future missions such as PLATO 2.0 in detecting solar system analogues. We use real-world solar data and end-to-end simulations to explore the stellar and instrumental noise properties. By injecting and retrieving planets, rings and moons of our own solar system, we show that the discovery of Venus- and Earth-analogues transiting G-dwarfs like our Sun is feasible at high S/N after collecting 6yrs of data, but Mars and Mercury will be difficult to detect due to stellar noise. In the best cases, Saturn's rings and Jupiter's moons will be detectable even in single transit observations. Through the high number (>1 bn) of observed stars by PLATO 2.0, it will become possible to detect thousands of single-transit events by cold gas giants, analogue to our Jupiter, Saturn, Uranus and Neptune. Our own solar system aside, we also show, through signal injection and retrieval, that PLATO 2.0-class photometry will allow for the secure detection of exomoons transiting quiet M-dwarfs. This is the first study analyzing in-depth the potential of future missions, and the ultimate limits of photometry, using realistic case examples.

Here we report on the first successful exoplanet transit observation with the Stratospheric Observatory for Infrared Astronomy (SOFIA). We observed a single transit of the hot Jupiter HD 189733 b, obtaining two simultaneous primary transit lightcurves in the B and z' bands as a demonstration of SOFIA's capability to perform absolute transit photometry. We present a detailed description of our data reduction, in particular the correlation of photometric systematics with various in-flight parameters unique to the airborne observing environment. The derived transit depths at B and z' wavelengths confirm a previously reported slope in the optical transmission spectrum of HD 189733 b. Our results give new insights to the current discussion about the source of this Rayleigh scattering in the upper atmosphere and the question of fixed limb darkening coefficients in fitting routines.

[Abridged] We have only been able to comprehensively characterize the atmospheres of a handful of transiting planets, because most orbit faint stars. TESS will discover transiting planets orbiting the brightest stars, enabling, in principle, an atmospheric survey of 10^2 to 10^3 bright hot Jupiters and warm sub-Neptunes. Uniform observations of such a statistically significant sample would provide leverage to understand---and learn from---the diversity of short-period planets. We argue that the best way to maximize the scientific returns of TESS is with a follow-up space mission consisting of a ~1 m telescope with an optical--NIR spectrograph: it could measure molecular absorption for non-terrestrial planets, as well as eclipses and phase variations for the hottest jovians. Such a mission could observe up to 10^3 transits per year, thus enabling it to survey a large fraction of the bright (J<11) TESS planets. JWST could be used to perform detailed atmospheric characterization of the most interesting transiting targets (transit, eclipse, and---when possible---phase-resolved spectroscopy). TESS is also expected to discover a few temperate terrestrial planets transiting nearby M-Dwarfs. Characterizing these worlds will be time-intensive: JWST will need months to provide tantalizing constraints on the presence of an atmosphere, planetary rotational state, clouds, and greenhouse gases. Future flagship missions should be designed to provide better constraints on the habitability of M-Dwarf temperate terrestrial planets.

EXONEST is an algorithm dedicated to detecting and characterizing the photometric signatures of exoplanets, which include reflection and thermal emission, Doppler boosting, and ellipsoidal variations. Using Bayesian Inference, we can test between competing models that describe the data as well as estimate model parameters. We demonstrate this approach by testing circular versus eccentric planetary orbital models, as well as testing for the presence or absence of four photometric effects. In addition to using Bayesian Model Selection, a unique aspect of EXONEST is the capability to distinguish between reflective and thermal contributions to the light curve. A case-study is presented using Kepler data recorded from the transiting planet KOI-13b. By considering only the non-transiting portions of the light curve, we demonstrate that it is possible to estimate the photometrically-relevant model parameters of KOI-13b. Furthermore, Bayesian model testing confirms that the orbit of KOI-13b has a detectable eccentricity.

We introduce a new thermochemical kinetics and photochemical model. We use high-temperature bidirectional reaction rates for important H, C, O and N reactions (most importantly for CH$_4$ to CO interconversion), allowing us to attain thermochemical equilibrium, deep in an atmosphere, purely kinetically. This allows ab initio chemical modeling of an entire atmosphere, from deep-atmosphere thermochemical equilibrium to the photochemically dominated regime. We use our model to explore the atmospheric chemistry of cooler ($T_{eff} < 10^3$ K) extrasolar giant planets. In particular, we choose to model the nearby hot Neptune GJ436b, the only planet in this temperature regime for which spectroscopic measurements and estimates of chemical abundances now exist. Recent {\it Spitzer} measurements with retrieval have shown that methane is driven strongly out of equilibrium and is deeply depleted on the dayside of GJ 436b, whereas quenched carbon monoxide is abundant. This is surprising because GJ 436b is cooler than many of the heavily irradiated hot Jovians and thermally favorable for CH$_4$, and thus requires an efficient mechanism for destroying it. We include realistic estimates of ultraviolet flux from the parent dM star GJ 436, to bound the direct photolysis and photosensitized depletion of CH$_4$. While our models indicate fairly rich disequilibrium conditions are likely in cooler exoplanets over a range of planetary metallicities, we are unable to generate the conditions for substantial CH$_4$ destruction. One possibility is an anomalous source of abundant H atoms between 0.01-1 bars (which attack CH$_4$), but we cannot as yet identify an efficient means to produce these hot atoms.

Detection of molecules using infrared spectroscopy probes the conditions and compositions of exoplanet atmospheres. Water (H2O), methane (CH4), carbon dioxide (CO2), and carbon monoxide (CO) have been detected in two hot Jupiters. These previous results relied on space-based telescopes that do not provide spectroscopic capability in the 2.4 - 5.2 micron spectral region. Here we report ground-based observations of the dayside emission spectrum for HD 189733b between 2.0-2.4 micron and 3.1-4.1 micron, where we find a bright emission feature. Where overlap with space-based instruments exists, our results are in excellent agreement with previous measurements. A feature at ~3.25 micron is unexpected and difficult to explain with models that assume local thermodynamic equilibrium (LTE) conditions at the 1 bar to 1 x 10-6 bar pressures typically sampled by infrared measurements. The most likely explanation for this feature is that it arises from non-LTE emission from CH4, similar to what is seen in the atmospheres of planets in our own Solar System. These results suggest that non-LTE effects may need to be considered when interpreting measurements of strongly irradiated exoplanets.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

Using the NICMOS instrument on the Hubble Space Telescope, we have measured the dayside spectrum of HD 209458b between 1.5--2.5 microns. The emergent spectrum is dominated by features due to the presence of methane (CH4) and water vapor (H2O), with smaller contributions from carbon dioxide (CO2). Combining this near-infrared spectrum with existing mid-infrared measurements shows the existence of a temperature inversion and confirms the interpretation of previous photometry measurements. We find a family of plausible solutions for the molecular abundance and detailed temperature profile. Observationally resolving the ambiguity between abundance and temperature requires either (1) improved wavelength coverage or spectral resolution of the dayside emission spectrum, or (2) a transmission spectrum where abundance determinations are less sensitive to the temperature structure.

We present first results of an exploratory study to use integral field spectroscopy to observe extrasolar planets. We focus on transiting "Hot Jupiters" and emphasize the importance of observing strategy and exact timing. We demonstrate how integral field spectroscopy compares with other spectroscopic techniques currently applied. We have tested our concept with a time series observation of HD209458b obtained with SINFONI at the VLT during a superior conjunction.