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.

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.

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.

In this paper we review the concepts of Bayesian evidence and Bayes factors, also known as log odds ratios, and their application to model selection. The theory is presented along with a discussion of analytic, approximate and numerical techniques. Specific attention is paid to the Laplace approximation, variational Bayes, importance sampling, thermodynamic integration, and nested sampling and its recent variants. Analogies to statistical physics, from which many of these techniques originate, are discussed in order to provide readers with deeper insights that may lead to new techniques. The utility of Bayesian model testing in the domain sciences is demonstrated by presenting four specific practical examples considered within the context of signal processing in the areas of signal detection, sensor characterization, scientific model selection and molecular force characterization.

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.

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.

The study of exoplanets (planets orbiting other stars) is revolutionizing the way we view our universe. High-precision photometric data provided by the Kepler Space Telescope (Kepler) enables not only the detection of such planets, but also their characterization. This presents a unique opportunity to apply Bayesian methods to better characterize the multitude of previously confirmed exoplanets. This paper focuses on applying the EXONEST algorithm to characterize the transiting short-period-hot-Jupiter, HAT-P-7b. EXONEST evaluates a suite of exoplanet photometric models by applying Bayesian Model Selection, which is implemented with the MultiNest algorithm. These models take into account planetary effects, such as reflected light and thermal emissions, as well as the effect of the planetary motion on the host star, such as Doppler beaming, or boosting, of light from the reflex motion of the host star, and photometric variations due to the planet-induced ellipsoidal shape of the host star. By calculating model evidences, one can determine which model best describes the observed data, thus identifying which effects dominate the planetary system. Presented are parameter estimates and model evidences for HAT-P-7b.

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.