A massive planet closely orbiting its host star creates tidal forces that distort the typically spherical stellar surface. These distortions, known as ellipsoidal variations, result in changes in the photometric flux emitted by the star, which can be detected within the data from the Kepler Space Telescope. Currently, there exist several models describing such variations and their effect on the photometric flux. By using Bayesian model testing in conjunction with the Bayesian-based exoplanet characterization software package EXONEST, the most probable representation for ellipsoidal variations was determined for synthetic data and the confirmed hot Jupiter exoplanet Kepler-13Ab.The most preferred model for ellipsoidal variations observed in the Kepler-13 light curve was determined to be EVIL-MC. Among the trigonometric models, the Modified Kane & Gelino model provided the best representation of ellipsoidal variations for the Kepler-13 system and may serve as a fast alternative to the more computationally intensive EVIL-MC. The computational feasibility of directly modeling the ellipsoidal variations of a star are examined and future work is outlined. Providing a more accurate model of ellipsoidal variations is expected to result in better planetary mass estimations.

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.