# EXONEST: The Bayesian Exoplanetary Explorer

###### Abstract

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étaires) 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.

## 1 Introduction

We are currently enjoying an unprecedented era of exploration and discovery. July of 2015 saw the New Horizons probe’s fly-by of Pluto and Charon, which marked the end of mankind’s initial exploration of the solar system. We are now beginning to explore the neighboring star systems by discovering and characterizing their planets (exoplanets). As stated in the 2014 NASA Strategic Plan, “We are navigating a voyage of unprecedented scope and ambition: seeking to discover and study planets orbiting around other stars and to explore whether they could harbor life” NASA (2014). To date, we have identified 4496 exoplanet candidates of which 3502 are confirmed planets in over 1648 star systems Caltech (2017).

One of the most successful missions to date is the Kepler Space Telescope (Kepler), which was designed to monitor the light intensity (photometry) from approximately 150,000 stars in the constellations of Cygnus and Lyra Borucki et al. (2010). While the characterization of exoplanets in terms of the chemistry and the detection of life will require a serious commitment to spectroscopic studies, photometry will continue to play an important role in exoplanet detection and characterization. Future surveys like the Transiting Exoplanet Survey Satellite (TESS) Ricker et al. (2010); Placek et al. (2016) will pave the way for higher precision photometry missions like the James Webb Space Telescope (JWST) Beichman et al. (2014); Placek et al. (2017), the CHaracterising ExOPlanets Satellite (CHEOPS) ESA (2016a); CHEOPS Science Team - University of Bern (2017); Broeg et al. (2013); Placek et al. (2017), the PLAnetary Transits and Oscillations (PLATO) mission ESA (2016b); ESA Study Team and PLATO Science Study Team (2010); Hippke and Angerhausen (2015) and the Wide Field Infrared Survey Telescope (WFIRST) NASA (2016). Increasingly accurate representations of photometric effects will become more important with these higher precision missions.

In this paper, we summarize the EXONEST software package for detecting and characterizing exoplanets Placek et al. (2014); Placek (2014), as well as introduce our new efforts toward more careful modeling of reflected light, refracted light and atmospheric effects and their incorporation into the next generation of the Exoplanetary Explorer software package. We are currently focused on porting and expanding the MATLAB-based EXONEST software package into the next generation Exoplanetary Explorer software package, which will accommodate multiple planet systems, three-body and multi-body orbital mechanics, as well as more subtle photometric effects. The Exoplanetary Explorer will be a Python-based open source project with the capability to employ third-party plug-and-play models of exoplanet-related photometric effects.

## 2 Exonest

The EXONEST software package (Figure 1) was originally developed as a part of Ben Placek’s Ph.D. thesis Placek et al. (2014); Placek (2014). Written in MATLAB, EXONEST consists of a core set of stellar and planetary models, a nested sampling-based Bayesian inference engine that can utilize the original nested sampling algorithm Skilling (2004, 2006); Sivia and Skilling (2006), the MultiNest variant Feroz et al. (2009, 2011, 2013) and the Metropolis–Hastings Markov chain Monte Carlo (MCMC) sampling algorithm Metropolis et al. (1953); Hastings (1970). EXONEST can also accept a set of plug-and-play models, which consist of orbital models, instrument likelihood functions, transit models, a set of photometric models and any additional user-defined models. This is currently implemented by modifying the code by calling one or more subroutines. In the future, these options will be selected, or unselected, via a run file and/or a graphical user interface.

EXONEST incorporates an efficient orbit integrator, as well as detailed models for four photometric effects: reflected light, thermal emissions, Doppler boosting or beaming and ellipsoidal variations of the host star Placek et al. (2014); Placek (2014). Both likelihood functions and priors on the parameter values can be specified.

In the following sections, we will summarize the orbital and photometric models employed by EXONEST and our current work to extend and improve them. We will then conclude by discussing the Bayesian inference engine and the need for improvement.

## 3 Orbital Models

EXONEST incorporates an efficient orbit integrator that can accommodate both circular and elliptical orbits Brown (2009), known as Keplerian orbits, producing both the star-planet distance and the orbital phase as a function of time. The current implementation of EXONEST is designed to accommodate one exoplanet per system so that only one period of the orbit need be integrated. The data are phase-wrapped to decrease computation time.

The basic anatomy of a Keplerian orbit is illustrated in Figure 2. The star is situated at the focus of the elliptical orbit. The periastron (point of closest approach) may be rotated by an angle , which is called the argument of the periastron, and the position of the planet along the orbit is designated by the angle from periastron, , which is called the true anomaly. The length of the semimajor axis, , is used to characterize the star-planet distance, , which varies as a function of the true anomaly according to:

(1) |

where is the eccentricity of the orbit.

The upcoming re-write of the EXONEST code will accommodate multiple exoplanets in different orbits. This is accomplished by integrating a single period of each exoplanet orbit and labeling the time series of the data in terms of the orbital phase of each orbit.

A full three-body orbital model has been implemented for EXONEST. The model can handle cases where the planets do or do not orbit in a common plane. Despite its utility, the model is computationally expensive as the initial conditions of each planet must be fully parameterized, thus significantly increasing the number of parameters to be estimated. In addition, orbits in three-body situations are, in general, not periodic. Unlike two-body periodic elliptical (Keplerian) orbits, which need only be integrated over a single period to make predictions for the entire dataset, three-body orbits must be integrated over the entire length of the dataset, dramatically increasing the computations necessary to generate the predicted photometric flux.

## 4 Photometric Effects

In this section, we consider the main photometric effects that are detectable with instruments exhibiting precision ranging from 10 parts per million to 30 parts per million. These photometric effects include both planetary components and stellar components. We will begin with the planetary components by discussing transits, which produce relatively large photometric signals. This is followed by a discussion of reflected light, thermal emissions, refracted light, and atmospheric effects. We then conclude with the stellar components, which consist of Doppler boosting and ellipsoidal variations.

Figure 3 illustrates the recorded data and light curve fit from the exoplanet HAT-P-7b Borucki et al. (2009). The primary and secondary eclipses are clearly seen at orbital phase zero days and 1.1 days, respectively. In Figure 3B, planetary flux variations, such as reflected light, are evident in the sinusoidally-varying flux and in the presence of the secondary eclipse as the planet passes behind the host star. The fact that the flux during the secondary eclipse dips below the baseline (red line) indicates that the planet is also emitting a great deal of thermal radiation, which is hidden as the planet passes behind the host star.

### 4.1 Transits and Eclipses

The most prominent photometric effect is known as the transit or primary transit, which occurs when the planet passes in front of the host star blocking a fraction of the starlight Seager and Mallěn-Ornelas (2003). The primary transit of the exoplanet HAT-P-7b can be seen in Figure 3 at an orbital phase of zero days and again at 2.2 days, which is its orbital period. The bottoms of the primary transits are typically rounded because the light emitted by the star is not uniform across the stellar disk, but is darker near the limbs as one is looking through the cooler upper regions of the photosphere Mandel and Agol (2002). This effect, known as limb darkening, is illustrated by the coloration of the stars on the left side of Figure 4A,B. The transit depth is proportional to the ratio of cross-sectional areas of the planet and star:

where is the radius of the planet and is the radius of the star.

Transits are relatively rare since the planet must pass between the stellar disk and the observer. Borucki and Summers found the probability that a distant observer could witness an exoplanetary transit to be Borucki and Summers (1984). This expression was found by considering the solid angle subtended by the exoplanet’s shadow in a circular orbit. Here, we consider eccentric orbits with eccentricity , and more-or-less follow the derivation by Barnes Barnes (2007), taking care to rigorously treat the probability distribution associated with the true anomaly, the argument of the periastron and the inclination.

To find the probability of a transit, we first consider circular orbits and work out the prior probability of observing a planet at a given angular position defined by the true anomaly and the inclination, . It is expected that there is no preferred orientation for the inclination of a planet’s orbital plane. Moreover, for a circular orbit, the planet orbits at a constant speed so that the prior probability of the true anomaly is also constant. As a result, the joint prior probability is constant, , where it should be noted that in our expressions for the probability, , we adopt the convention of Jaynes Jaynes (2003) and Sivia and Skilling Sivia and Skilling (2006) in which all probabilities are written as being conditional on all of the prior information about the problem:

(2) |

so that:

(3) |

In the case of eccentric orbits, one must also consider the argument of the periastron since:

(4) |

and the speed of the planet changes as it orbits so that the probability density of finding the planet at a true anomaly of varies with in an elliptical orbit: where .

Consider a change in the viewing angle by , so that and . It is then clear that:

(5) |

Moreover, since:

and since is periodic, the two integrals over a single period:

(6) |

and:

(7) |

are equal so that:

(8) |

for arbitrary angle . This implies that is constant so that (3) is generally true with:

(9) |

for both circular and elliptical orbits.

The planet will be observed as transiting only for inclinations between where:

and is the star-planet distance. Therefore, the probability of a transit:

(10) |

which is:

Now, using the fact that the star-planet distance varies as a function of the true anomaly as described in (1), the equation above becomes:

so that in the case of an eccentric orbit, the probability of observing the transit of a planet for which the eccentricity is:

(11) |

whereas for a circular orbit, or an orbit for which the eccentricity is zero, , the transit probability is:

(12) |

One can now compute the probability that a distant observer would view a planet in a circular orbit as transiting its host star. Table 1 lists these probabilities for several planets in our solar system for which the Sun has a radius of AU. According to (12), these probabilities fall off as , indicating that the transit method of exoplanet detection is heavily biased toward the detection of closely-orbiting exoplanets.

Planet | Semi-Major Axis (AU) | Transit Probability |
---|---|---|

Mercury | 0.39 | 1.18% |

Venus | 0.72 | 0.64% |

Earth | 1.00 | 0.46% |

Mars | 1.53 | 0.30% |

Jupiter | 5.20 | 0.09% |

Saturn | 9.54 | 0.05% |

Solar Radius | 0.0046 AU |

The probability of an alien civilization being able to observe Earth transiting the Sun (within the Earth transit zone (ETZ)) is only 0.46%. Since the average stellar density in the galactic neighborhood surrounding the Sun is approximately 0.004 stars per cubic light year Gregersen (2010), there are about 17,000 stars within 100 LY of the Sun. Of these 17,000 stars, only about 78 of them are positioned so that Earth could be observed transiting the Sun. A recent study utilized the Hipparcus database to identify 37 K and 45 G dwarf stars in the ETZ within 1000 parsecs (3261 LY) Heller and Pudritz (2016).

The secondary eclipse is similar to the transit. It occurs when the planet passes behind the host star so that the planetary flux, which is composed of both reflected light and thermal emissions, is blocked. The magnitude of the secondary eclipse is dependent on the magnitude of the planetary flux, which, for a relatively cool planet, is bounded by the maximum amount of reflected light (14):

(13) |

where we are considering a circular orbit of radius , is the geometric albedo, is the radius of the planet, is the radius of the host star and is the stellar flux. Since typically , secondary eclipses are not always detectable in the recorded data. EXONEST models secondary eclipses using the methodology developed by Mandel and Agol Mandel and Agol (2002) applied to a uniform source. An example of a secondary eclipse can be seen in Figure 3 at the orbital phase of 1.1 days.

### 4.2 Reflected Light

Planets reflect light both from cloud tops and from solid or liquid surfaces, and the degree to which they reflect light is quantified by their albedo Jenkins and Doyle (2003); Seager (2010); Perryman (2011); Placek (2014); Placek et al. (2014). EXONEST employs the simplest model for reflected light in which the planet is assumed to be diffusely isotropically reflecting, so that the brightness does not depend on the angle of observation. In addition, it is assumed that the planet is situated sufficiently distant from the host star that the starlight impinges on the planet in parallel rays so that half of the planet is in daylight and the other half of the planet is in night.

Given these assumptions, the stellar flux reflected by the planet can be written as Sobolev (1975); Carter (ress):

(14) |

where is the radius of the planet, is the geometric albedo and is the time-varying orbital phase angle. This result appears to work well for many exoplanets.

#### New Efforts to Properly Model Reflected Light

However, it is known that a significant number of planets (mainly hot Jupiters) closely orbit giant stars so that the assumption that the starlight impinges on the planet in the form of parallel rays is false. Instead, the planet can be divided into three zones. The first zone is that of full daylight where the entire apparent disk of the star can be seen from the surface. The second zone is a penumbral zone where only a portion of the star is visible as it appears to be setting or rising. The third zone is the night zone where the star is not visible. There are a number of cases in which the penumbral zone extends well into what would normally be the night zone. For example, as much as 70% of Kepler-91b’s surface is lit by starlight. This situation is illustrated in Figure 4B. This has significant implications for the primary transit since a fraction of the backside of the planet may be partially illuminated in the penumbral zone, thus compensating somewhat for the reduced starlight due to the transit. As a result, unless the illumination geometry is properly modeled, primary transit depths will be less extensive resulting in an under-estimation of the planetary radius or an over-estimation of the nightside thermal flux followed by an over-estimation of nightside temperatures.

We are currently working to develop reflected light models for EXONEST that can accommodate the case in which planets are illuminated by parallel rays and the case in which planets are illuminated while closely-orbiting giant stars Ḱopal (1954, 1959); Carter (ress). The fractional surface areas in the fully-illuminated zone, partially-illuminated penumbral zone and night zone are given by Carter (ress):

(15) | ||||

(16) | ||||

(17) |

Since these surface areas depend on the star-planet separation distance , eccentric orbits will require the use of the more general reflectance model.

A quick survey of the exoplanet archive Caltech (2017) reveals 1306 exoplanets for which the radius of the star, , the radius of the planet, , and the semi-major axis, , are reported. Figure 5 illustrates histograms for these 1306 exoplanets revealing the percent of the surface area that is in full daylight, in a penumbral zone and in completely dark night. It is most common for planets to have a 5% penumbral zone, and more than one quarter of the planets have penumbral zones that are larger than 5%. There are nine planets for which the penumbral zone covers about 30% of the planet’s surface and two planets for which the penumbral zone covers 40% of the planet’s surface! Of these 1306 exoplanets, 39 of them have more than 60% of their surface illuminated (dayside zone plus penumbral zone) by their host star. Another survey of the 1073 confirmed Kepler exoplanets that have orbital periods shorter than 10 days indicates that 88 of those planets (about 8% of them) have more than 60% of their surface illuminated (dayside zone plus penumbral zone) by their host star. The significance of this effect in the context of the reflected light and the analysis of Kepler data will be explored in a future paper.

We believe that improperly modeling stellar illumination was a problem in our earlier study involving a proposed Trojan partner to the planet Kepler-91b. In that analysis, the Bayesian evidence weighed heavily in favor of a Trojan partner, but the modeled planetary temperatures were unphysically high, which led us to discount the possibility of a Trojan partner Placek et al. (2015). The star Kepler-91 is an asymptotic giant branch star with a radius of approximately (where is the solar radius), and the semi-major axis of the hot Jupiter Kepler-91b is AU Esteves et al. (2015). The planet Kepler-91b has a radius of , where represents Jupiter’s radius Esteves et al. (2015). By Equations (15)–(17), this implies that of Kepler-91b is in full daylight, is in partial penumbral light and is in full darkness, so that approximately of the planet is illuminated at any point in time (as is illustrated in Figure 4B). By not properly accounting for this illumination, the modeled thermal emissions would be forced to compensate with increased day-side and night-side temperatures possibly explaining the unphysically high temperatures obtained in the analysis Placek et al. (2015), which led to the Trojan hypothesis being discounted. This is supported by the fact that other studies of this system have led to unphysically high temperatures. For example, a study of KOI-2133b (Kepler-91b) by Esteves et al. Esteves et al. (2013) found that the nightside temperature was K, which was greater than the expected equilibrium temperature of K, leading them to hypothesize that Kepler-91b, which was later re-confirmed as a planet Sliski and Kipping (2014), was self-luminous.

### 4.3 Thermal Emissions

The modeling of thermal emissions is similar to that of reflected light. Current models assume that there is a day-side and a night-side and that these two sides are characterized by two effective temperatures and , respectively.

The planet is assumed to radiate as a thermal blackbody so that the detected thermal flux is given by Charbonneau et al. (2005):

(18) |

where is the effective temperature of the radiating region, represents the spectral radiance of a blackbody at temperature :

(19) |

is Planck’s constant, is Boltzmann’s constant, and is the Kepler response as a function of wavelength Van Cleve and Caldwell (2009).

The relative day-side and night-side fluxes are given by:

(20) | ||||

(21) |

with a total relative flux of:

(22) |

The similarity of the thermal emission phase curve to the reflected light phase curve means that for a circular orbit, the models are degenerate so that it is not possible, using single-bandpass photometry, to distinguish thermal emission from reflected light. Only for significantly eccentric orbits, with , in which the planet-star distance changes, thus affecting reflected light, can thermal flux be modeled independently of reflected light Placek et al. (2014), assuming that the temperatures on the planet are relatively constant. Observations in multiple bandpasses, such as combining data from both Kepler and TESS, will enable thermal emissions to be distinguished from reflected light Placek et al. (2016). Just as in the case of reflected light, a closely-orbiting planet will have three zones: a dayside zone, a penumbral zone and a nightside zone. While one could model each zone with a separate average temperature, previous studies have found that the contribution of the penumbral zone to thermal emissions is negligible Léger et al. (2011); Maurin et al. (2012).

### 4.4 New Efforts to Model Refracted and Forward Scattered Light

The situation is significantly more complicated for a planet or moon that possesses an atmosphere. As the planet approaches the primary transit, it becomes more and more backlit. In addition to reflected light from atmospheric clouds or the surface, one may potentially detect light that is both forward scattered by the atmosphere and refracted through the atmosphere Sidis (2010); Muñoz et al. (2012); Misra et al. (2014); Misra (2014); Bétrémieux and Kaltenegger (2015); Muñoz and Cabrera (2017). One of these effects is illustrated in Figure 6 where in Panel (A), we see an image taken by the Cassini probe of Saturn’s moons Rhea (upper left), Titan (center right) and Mimas (lower center). Both Rhea and Mimas are airless worlds, and one sees only reflected light from their surfaces. Titan, on the other hand, is a large moon with an extended nitrogen and methane atmosphere and, more critically, rich in strongly forward scattering haze. Here, one sees that the starlight entering the atmosphere is forward scattered by the haze particles, thereby forming an illuminated arc much brighter and more extensive than on the other two moons for which such a ring is not possible Muñoz et al. (2017). Forward scattered by hazy exoplanet atmospheres, if occurring, will be more important for close-in planets, whereas the refraction signal is likely more important for far-out transiting planets and immediately before/after transit, where refraction in the upper atmosphere may cause shoulders in the pre-/post-transit light curves Muñoz et al. (2012); Misra et al. (2014); Misra (2014); Dalba (2017).

The increase in photometric flux due to refraction is on the order of Sidis (2010):

(23) |

where is the scale height of the planetary atmosphere, which is given by:

(24) |

where is Boltzmann’s constant, is the average atmospheric temperature, is the mean molecular weight of the atmospheric constituents and is the gravitational acceleration at the surface of the planet Sidis (2010). Typically, , so that by (23), the refractive effects are much smaller than the primary transit depth:

(25) |

For example, in the case of the Earth ( K, m s, g mol), the characteristic atmospheric height is estimated to be m, which is reasonable considering that the troposphere has an average height of m. As a result, the photometric flux of light refracting through Earth’s atmosphere at the point of occultation is about times, or three orders of magnitude, smaller than the primary transit depth.

In the event that the photometric effects of refracted light can be detected, this would provide information about the nature of the exoplanetary atmosphere.

### 4.5 New Efforts to Model Atmospheric Effects

Atmospheric effects can play a significant role in photometry. For example, atmospheric clouds can shift a planet’s highly reflective bright spot away from the substellar point. Similarly, atmospheric superrotation can shift a planet’s thermally bright hotspot away from the substellar point. Such effects, which can be modeled as an angular shift in the reflectance or thermal phase curves, can reveal a great deal about the exoplanetary atmosphere Showman and Polvani (2011); Esteves et al. (2015); Faigler and Mazeh (2015); Placek et al. (2017). More exotic effects are possible. For example, it is thought that supersonic flow and shocks can develop on the dayside of irradiated hot Jupiters upstream of the substellar point and that these shocks will affect the position of the peak of the thermal phase curve Heng (2012). Heavily-irradiated hot Jupiters are expected to have an electrically-conductive atmosphere due to ionization. As a result, exoplanetary magnetic fields are expected to affect atmospheric flow Batygin et al. (2013) and in some cases, magnetically-dragged winds are expected to reduce the longitudinal offset of the thermal hotspot Menou (2012); Batygin et al. (2013); Rauscher and Menou (2013).

Other deviations in the reflectance and thermal phase curves may reveal more about planetary atmospheres, such as inhomogeneous cloud cover across the planet’s surface Demory et al. (2013); Heng and Demory (2013), atmospheric variability Armstrong et al. (2016), evolving weather patterns Pallé et al. (2008); Artigau et al. (2009); Gillon et al. (2013), diurnal atmospheric effects Ford et al. (2001) or seasonal atmospheric effects Van Eylen et al. (2013).

### 4.6 Doppler Boosting/Beaming

The gravitational influence of a planet causes both the star and the planet(s) to orbit a common center of mass. As a result, stars with planets tend to wobble. This is the basis of the radial velocity technique in which the radial velocity of the star is carefully measured using spectroscopy and knowledge about the Doppler shift to detect the presence of planets. While the radial velocity technique relies on spectroscopy, there is a photometric component that arises due to relativistic considerations related to stellar aberration in which an object that radiates light uniformly at rest will, in motion, tend to radiate more light in the direction of motion and less light opposite that direction Rybicki and Lightman (2008). This effect is called Doppler boosting or Doppler beaming.

The effect can be derived using special relativity resulting in the following formula for the boosted flux Rybicki and Lightman (2008); Placek (2014); Placek et al. (2014):

(26) |

where , and where is the speed of light and is the stellar flux in the reference frame of the star. Since the motion of the star will involve velocities no more than 100–1000 m/s (), it is perhaps surprising that this relativistic effect could be detectable. However, the fact that the star is approximately -times brighter than the reflected light from a planet means that this effect can be on the order of that of reflected light. Given these speeds, a non-relativistic approximation is reasonable:

(27) |

where represents the radial component of the stellar velocity as viewed from Earth scaled by the speed of light so that:

(28) |

where is the radial velocity of the planet from the perspective of Earth:

(29) |

in which is the true anomaly as a function of time, is the argument of the periastron, is the eccentricity (see Figure 1) and is the radial velocity semi-amplitude Placek (2014) given by:

(30) |

in which is Newton’s gravitational constant, is the period of the planet, is the mass of the planet, is the inclination of the orbital plane with respect to Earth and is the mass of the star. Note that when the inclination of the orbital plane is zero, , the orbit is face on (in the projected plane of the sky), and both the radial velocity semi-amplitude and the radial velocity are zero, resulting in no Doppler effect. Furthermore, as the planet orbits the star, the true anomaly , which is an angle, advances. There will be two times per orbital period at which the trigonometric factor goes to zero, resulting in zero radial velocity and therefore zero Doppler effect. This is most easily seen by considering a circular orbit for which the eccentricity is zero, , and the argument of the periastron is zero, . In the case of this circular orbit, the radial velocity is zero when the true anomaly is and , which occurs when the planet is moving in the projected plane of the sky. Since the Doppler boosting effect involves both the mass of the planet and the inclination angle of the orbital plane, this phenomenon facilitates the estimation of these two quantities.

### 4.7 Tidal Forces and Ellipsoidal Variations

Tidal interactions between a massive object and a star create distortions on the star (see Figure 7), which are potentially observable by the Kepler Space Telescope. Known as ellipsoidal variations, this photometric effect appears with two maxima per orbit, each occurring when the largest cross-sectional area of stellar surface is observed. Bayesian model testing provides a rigorous framework to compare the various representations of ellipsoidal variations. A preferred representation of ellipsoidal variations will become increasingly important for exoplanet classification with higher precision photometry.

#### 4.7.1 Trigonometric Models

The BEaming, Ellipsoidal and Reflection (BEER) model, which was developed by Faigler and Mazeh Faigler and Mazeh (2007), models the ellipsoidal variations as being proportional to the cosine of twice the phase angle:

(31) |

in which on the left-hand side, the flux due to the ellipsoidal variations is normalized by the stellar flux , is the planet’s mass, is the star’s mass, is the star’s radius, is the star-planet distance, is the inclination of the orbital plane and is the orbital phase [0, ]. The BEER model allows for positive and negative values for the amount of ellipsoidal variation observed. The effects of linear gravity-darkening and limb-darkening, , are modeled by Morris (1985):

(32) |

in which the linear limb-darkening and gravity-darkening coefficients are represented by and , respectively. Estimates of the coefficients and are provided by modeling metallicity and effective stellar temperature Heyrovský (2007); Sing (2010).

The Kane and Gelino model hypothesizes ellipsoidal variations to be proportional to the projected separation distance between the star and planet Kane and Gelino (2012), which, when substituted for the orbital phase, become:

(33) |

However, this model suffers from a discontinuity in the first derivative, which is not seen in simulations of ellipsoidal variations using gravitational isopotentials.

A modified form, labeled Kane and Gelino (Mod), was introduced by Placek et al. Placek et al. (2014), which removes the discontinuity by setting the effect proportional to the square of the projected separation distance:

(34) |

Unlike the BEER model, which produces negative values, the two Kane and Gelino models have a minimum value of zero.

#### 4.7.2 Direct Modeling: EVIL-MC

The Ellipsoidal Variations Induced by a Low-Mass Companion (EVIL-MC) model was developed by Jackson et al. Jackson et al. (2012) to directly model tidal effects on the stellar atmosphere. EVIL-MC was modified for use in EXONEST Gai (2016); Gai and Knuth (2017) in which a geodesic structure of 640 triangles is used to represent the stellar surface as illustrated in Figure 7. The deviations from sphericity are computed for each point on the surface model. The observed flux, , is computed by summing the blackbody flux within the Kepler bandpass (18) over the observed surface area. The ellipsoidal variation effect is then computed with:

(35) |

where is the flux from a spherical star and is the flux from the ellipsoidal star. While direct modeling should provide a more accurate representation, EVIL-MC is significantly (1000–10,000 times) more computationally expensive than the trigonometric models. The Kane and Gelino (Mod) model provides a computationally fast approximation to EVIL-MC.

#### 4.7.3 Model Testing

We have recently performed model testing on Kepler light curve data using EXONEST Gai (2016); Gai and Knuth (2017). EVIL-MC provided the best representation of ellipsoidal variations in the Kepler-13 system followed by the Kane and Gelino (Mod) and BEER models, which had similar evidence values. The values are reported in Table 2. The Kane and Gelino (Mod) model is closest to the EVIL-MC results and has slightly higher evidence over the others suggesting that, given its straight-forward analytic form, the Kane and Gelino (Mod) is the most effective model for this effect.

Model | Log (Evidence) | EVIL-MC’s Bayes’ Factor |
---|---|---|

EVIL-MC | 14,256.91 0.91 | 1 |

Kane and Gelino (Mod) | 14,256.71 0.73 | 2.20 |

BEER | 14,256.26 0.74 | 2.65 |

Kane and Gelino (2012) | 14,253.66 0.72 | 5.25 |

## 5 Bayesian Inference Engine

For a given exoplanetary model, the Bayesian inference engine allows us to compute the Bayesian evidence and the model parameter values along with their associated uncertainties. The current implementation allows for a choice of sampling methods: Metropolis–Hastings Markov chain Monte Carlo (MCMC) Metropolis et al. (1953); Hastings (1970) (which does not provide the Bayesian evidence), the original nested sampling algorithm Skilling (2004, 2006); Sivia and Skilling (2006) or the more recent variant called MultiNest Feroz et al. (2009, 2011, 2013). While the MultiNest algorithm has proven to shorten run times from several days to 10s of hours, it is usually not possible to operate with more than 75 samples, which severely reduces the precision to which the Bayesian evidence and the parameter values can be estimated.

Unfortunately, even by integrating a single period in cases with one exoplanet, the time required to evaluate the log likelihood is too great to allow for a rapid analysis. It may be of benefit to re-parameterize the eccentricity and argument of periapsis as suggested by Ford Ford (2006). In addition, since we have found almost no evidence of phase changes in our analyses to date, it has been suggested (P. Goggans, personal communication) that EXONEST be designed to utilize simulated annealing.

Model testing is a significant part of the EXONEST analysis pathway. By turning on-and-off various photometric models and comparing the Bayesian evidence in each case, one can build a case that a candidate exoplanet is a real planet and not a false positive Placek et al. (2014); Knuth et al. (2015).

### Priors and Likelihoods

EXONEST allows priors on the parameter values to be specified. However, since we are still learning a great deal about exoplanets, the assignment of an informative prior based on known exoplanets is potentially detrimental as it may bias against the discovery of new phenomena. For this reason, uniform or non-informative priors are most often used.

Likelihoods can also be specified for EXONEST. While Gaussian likelihoods are typically employed, there have been situations involving asymptotic giant branch stars in which likelihoods that accommodate correlated noise (red noise) are necessary Placek et al. (2015), prompting us to implement Sivia’s nearest-neighbor approach Sivia and Skilling (2006):

(36) |

where is the number of data points, is the noise variance, is the correlation strength, which is a parameter to be estimated, and depends on the sum of the square of the residuals by:

(37) |

The quantity is the sum of the first and the last squared residuals, and is the sum of the nearest neighbor squared residuals:

(38) |

As usual, the residuals are given by the difference between the modeled relative flux and the observed relative flux. In practice, one divides the observed flux by the average of the observed flux to obtain the observed relative flux. The residuals are typically given by:

(39) |

where is the observed or measured photometric flux at time and is the sum of the modeled relative photometric fluxes at time . Unfortunately, this practice is not ideal since it is not generally true that since the photometric effects (transits, etc.) do not all integrate to zero. In principle, one could subtract the flux during the secondary eclipse, since ellipsoidal variations will be minimal and there will be no planetary effects. However, it would be better to model the stellar flux and to marginalize over it.

It is a straightforward problem to incorporate additional data, such as radial velocity measurements, by adding terms to the log likelihood. This has been done in our study of the Kepler-91 system Placek et al. (2015).

## 6 Application: Kepler-76b

Kepler-76b is a transiting hot-Jupiter with a radius of approximately that orbits a star with an orbital period of days Faigler et al. (2013). This planet is known to exhibit reflection/thermal emission, Doppler boosting, as well as ellipsoidal variations (Faigler et al., 2013; Angerhausen et al., 2015; Esteves et al., 2015). This section outlines EXONEST’s model selection capabilities with fits of four different forward models that each account for different combinations of photometric effects. The simplest model includes only the reflection component (R), which in the case of a circular orbit subsumes thermal emission. The second model includes reflection and ellipsoidal variations (RE). The third model accounts for reflection, ellipsoidal variations and Doppler beaming (REB). Finally, the last model accounts for each of the three effects in the REB model in addition to an angular shift in the phase curve (reflection component), which is indicative of atmospheric superrotation (REB-SR). The data, which were obtained from the Mikulski Archive for Space Telescopes (MAST) managed by the Space Telescope Science Institute (STScI), consists of all available quarters of long-cadence Kepler Space Telescope data for Kepler-76b. The data were preprocessed to remove the primary transit so that only the overall phase curve and the secondary eclipse were fit by the EXONEST models. Note that even without considering the primary transit, excellent parameter estimates were obtained.

The results of the model testing are listed in Table 3. In this example, the most complex model (REB-SR), illustrated in Figure 8, was favored by both the Bayesian evidence () and the maximum log-likelihood () to explain the data. We can thus claim the detection of a shift in the phase curve by 8.53 0.12 degrees westward of the sub-stellar point, which is most likely caused by atmospheric superrotation. This result is consistent with previous studies of Kepler-76b, which yielded estimated westward shifts of 10.3 2.0 Faigler et al. (2013) and 8.28 1.22 (Esteves et al., 2015). The parameter estimates for the REB-SR model are listed in Table 4. Note that these estimates are comparable to those reported Faigler et al. (2013) when Kepler-76b was discovered.

Model | Z (10 | |
---|---|---|

R | 9.5987 | 9.5987 |

RE | 9.5822 | 9.5822 |

REB | 9.5685 | 9.5684 |

REB + SR | 9.5676 | 9.5676 |

Kepler-76b Parameter Estimates (REB + SR Model) | ||

Model Parameters | REB + SR Parameter Estimates | Literature Values |

1.35 0.05 | 1.25 0.08 | |

1.97 0.12 | 2.0 0.26 | |

77.9 0.03 | 78.0 0.2 | |

0.31 0.22 | 0.23 0.02, 0.22 0.02 |

## 7 Conclusions

EXONEST: The Exoplanetary Explorer is a software package for detecting and characterizing exoplanets from Kepler and CoRoT data. The package models a large number of photometric effects, each of which provides valuable information about the exoplanet and the exoplanetary system. New efforts include developing more precise models of reflected light, refracted light and atmospheric effects, as well as re-coding to transform EXONEST into a Python-based open source project with the capability to employ third-party plug-and-play models of exoplanet-related photometric effects. In addition, new planetary models are being developed that are relevant to rapidly rotating planets with shapes ranging from that of an oblate ellipsoid to pyriform shapes Abraham and Shaw (1983) and synestias Lock and Stewart (2017). These features are summarized in Table 5.

Modeled Effect | EXONEST | Under Study |
---|---|---|

Circular Orbits | X | |

Elliptical Orbits | X | |

Primary Transit | X | |

Secondary Eclipse | X | |

Limb Darkening | X | M |

Reflected Light | X | R |

Refracted Light | S | |

Thermal Emissions | X | R |

Doppler Beaming | X | |

Ellipsoidal Variations | X | M |

Atmospheric Effects | S | |

Multiple Planets | D | |

Planetary Shapes | S |

While it is true that by adding more parameters, one can over-fit the data by increasing the likelihood, for this reason, EXONEST employs Bayesian model testing, which enables one to evaluate the significance of each photometric effect. The model testing capabilities of EXONEST were demonstrated on Kepler Space Telescope of Kepler-76. Kepler-76b was demonstrated to exhibit several photometric effects: reflected light, Doppler boosting, ellipsoidal variations and atmospheric superrotation. Parameter estimates were reported that are in agreement with previous studies Faigler et al. (2013); Angerhausen et al. (2015); Esteves et al. (2015) using both radial velocity and photometric data. In addition, the model testing capabilities of EXONEST enable one to test photometric models against one another.

Future missions promise instruments of increasing sensitivity and detail with multiple spectral bands. Despite this, photometry will continue to be a valuable resource until detailed direct imaging is possible. In the meantime, we aim to include more relevant photometric effects into EXONEST and to develop the software package to work in conjunction with future multispectral studies.

## Acknowledgements

This paper was presented as an invited talk at the 37th International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering (MaxEnt 2017) in Jarinu SP, Brasil. Knuth would like to thank the organizers, Adriano Polpo and Julio Stern, for both the kind invitation and their generous hospitality. We also wish to thank Jon Jenkins, Doug Caldwell, John Skilling, Paul Goggans, Wesley Henderson, Ariel Caticha, Keith Earle, Oleg Lunin and Matthew Szydagis for interesting discussions and helpful advice, questions and comments. We would also like to thank the two anonymous reviewers whose comments served to improve the quality of this manuscript.

## Author Contributions

K.H.K. and B.P. conceived of EXONEST. K.H.K. advised the research and work performed by B.P., J.L.C., B.D., A.D.G., and B.C., and wrote the greater part of the manuscript. B.P. developed, coded, tested, and utilized EXONEST as a major component of his Ph.D. research. B.P. performed the analysis of Kepler-76b and wrote the corresponding section of the manuscript. D.A. contributed to the development of EXONEST and the inclusion of several photometric effects. As the focus of her Ph.D. research, J.L.C. has led the effort to more accurately model reflected light, and she contributed to the writing of the section on reflected light. As the focus of his M.S. thesis, B.D. implemented the three-body orbital code. As the focus of his M.S. thesis, A.D.G. coded and tested the efficacy of stellar ellipsoidal variation models, and wrote the corresponding section on ellipsoidal variations. B.C. assisted J.L.C. in the modeling of reflected light, and as the focus of his Ph.D. research, he has worked to characterize planetary shapes and their transit signatures.

## References

- Abraham and Shaw (1983) Abraham, R. H. and Shaw, C. D. (1983). Dynamics-The Geometry of Behavior: Vol.: 2: Chaotic Behavior. Aerial Press, Incorporated.
- Angerhausen et al. (2015) Angerhausen, D., DeLarme, E., and Morse, J. A. (2015). A Comprehensive Study of Kepler Phase Curves and Secondary Eclipses: Temperatures and Albedos of Confirmed Kepler Giant Planets. Publications of the Astronomical Society of the Pacific, 127, 1113–1130.
- Armstrong et al. (2016) Armstrong, D. J., de Mooij, E., Barstow, J., Osborn, H. P., Blake, J., and Saniee, N. F. (2016). Variability in the atmosphere of the hot giant planet HAT-P-7 b. Nature Astronomy, 1, 0004.
- Artigau et al. (2009) Artigau, É., Bouchard, S., Doyon, R., and Lafrenière, D. (2009). Photometric variability of the T2. 5 brown dwarf SIMP J013656. 5+ 093347: Evidence for evolving weather patterns. ApJ, 701(2), 1534.
- Barnes (2007) Barnes, J. W. (2007). Effects of orbital eccentricity on extrasolar planet transit detectability and light curves. Publications of the Astronomical Society of the Pacific, 119(859), 986.
- Batygin et al. (2013) Batygin, K., Stanley, S., and Stevenson, D. J. (2013). Magnetically controlled circulation on hot extrasolar planets. ApJ, 776, 53.
- Beichman et al. (2014) Beichman, C., Benneke, B., Knutson, H., Smith, R., Lagage, P.-O., Dressing, C., Latham, D., Lunine, J., Birkmann, S., Ferruit, P., et al. (2014). Observations of transiting exoplanets with the James Webb Space Telescope (JWST). PASP, 126(946), 1134.
- Bétrémieux and Kaltenegger (2015) Bétrémieux, Y. and Kaltenegger, L. (2015). Refraction in planetary atmospheres: improved analytical expressions and comparison with a new ray-tracing algorithm. MNRAS, 451(2), 1268–1283.
- Borucki and Summers (1984) Borucki, W. J. and Summers, A. L. (1984). The photometric method of detecting other planetary systems. Icarus, 58, 121–134.
- Borucki et al. (2009) Borucki, W. J., Koch, D., Jenkins, J., Sasselov, D., Gilliland, R., Batalha, N., Latham, D. W., Caldwell, D., Basri, G., Brown, T., Christensen-Dalsgaard, J., Cochran, W. D., DeVore, E., Dunham, E., Dupree, A. K., Gautier, T., Geary, J., Gould, A., Howell, S., Kjeldsen, H., Lissauer, J., Marcy, G., Meibom, S., Morrison, D., and Tarter, J. (2009). Kepler’s optical phase curve of the exoplanet HAT-P-7b. Science, 325, 709.
- Borucki et al. (2010) Borucki, W. J., Koch, D., Basri, G., Batalha, N., Brown, T., Caldwell, D., Caldwell, J., Christensen-Dalsgaard, J., Cochran, W. D., DeVore, E., Dunham, E. W., Dupree, A. K., Gautier, T. N., Geary, J. C., Gilliland, R., Gould, A., Howell, S. B., Jenkins, J. M., Kondo, Y., Latham, D. W., Marcy, G. W., Meibom, S., Kjeldsen, H., Lissauer, J. J., Monet, D. G., Morrison, D., Sasselov, D., Tarter, J., Boss, A., Brownlee, D., Owen, T., Buzasi, D., Charbonneau, D., Doyle, L., Fortney, J., Ford, E. B., Holman, M. J., Seager, S., Steffen, J. H., Welsh, W. F., Rowe, J., Anderson, H., Buchhave, L., Ciardi, D., Walkowicz, L., Sherry, W., Horch, E., Isaacson, H., Everett, M. E., Fischer, D., Torres, G., Johnson, J. A., Endl, M., MacQueen, P., Bryson, S. T., Dotson, J., Haas, M., Kolodziejczak, J., Van Cleve, J., Chandrasekaran, H., Twicken, J. D., Quintana, E. V., Clarke, B. D., Allen, C., Li, J., Wu, H., Tenenbaum, P., Verner, E., Bruhweiler, F., Barnes, J., and Prsa, A. (2010). Kepler planet-detection mission: Introduction and first results. Science, 327, 977–.
- Broeg et al. (2013) Broeg, C., Fortier, A., Ehrenreich, D., Alibert, Y., Baumjohann, W., Benz, W., Deleuil, M., Gillon, M., Ivanov, A., Liseau, R., Meyer, M., Oloffson, G., Pagano, I., Piotto, G., Pollacco, D., Queloz, D., Ragazzoni, R., Renotte, E., Steller, M., Thomas, N., and the CHEOPS Team (2013). CHEOPS a transit photometry mission for ESA’s small mission programme. EPJ Web of Conferences, 47.
- Brown (2009) Brown, R. A. (2009). Photometric orbits of extrasolar planets. The Astrophysical Journal, 702(2), 1237.
- Caltech (2017) Caltech (2017). Nasa exoplanet archive. https://exoplanetarchive.ipac.caltech.edu/. Accessed: August 1, 2017.
- Carter (ress) Carter, J. (In Progress). Estimation of Planetary Photometric Emissions for Extremely Close-in Exoplanets. Ph.D. thesis, University at Albany.
- Charbonneau et al. (2005) Charbonneau, D., Allen, L. E., Megeath, S. T., Torres, G., Alonso, R., Brown, T. M., Gilliland, R. L., Latham, D. W., Mandushev, G., O’Donovan, F. T., et al. (2005). Detection of thermal emission from an extrasolar planet. ApJ, 626(1), 523.
- CHEOPS Science Team - University of Bern (2017) CHEOPS Science Team - University of Bern (2017). CHEOPS mission homepage. http://cheops.unibe.ch/. Accessed: 2017-03-28.
- Dalba (2017) Dalba, P. A. (2017). Out-of-transit refracted light in the atmospheres of transiting and non-transiting exoplanets. ApJ, 848(2), 91.
- Demory et al. (2013) Demory, B.-O., De Wit, J., Lewis, N., Fortney, J., Zsom, A., Seager, S., Knutson, H., Heng, K., Madhusudhan, N., Gillon, M., et al. (2013). Inference of inhomogeneous clouds in an exoplanet atmosphere. ApJ Lett, 776(2), L25.
- ESA (2016a) ESA (2016a). Science & Technology: CHEOPS. http://sci.esa.int/cheops/. Accessed: 2016-09-30.
- ESA (2016b) ESA (2016b). Science & technology: PLATO. http://sci.esa.int/plato/. Accessed: 2016-09-30.
- ESA Study Team and PLATO Science Study Team (2010) ESA Study Team and PLATO Science Study Team (2010). PLATO science requirements document. ESA Public Documents.
- Esteves et al. (2015) Esteves, L., De Mooij, E., and Jayawardhana, R. (2015). Changing phases of alien worlds: Probing atmospheres of Kepler planets with high-precision photometry. ApJ, 804(2), 150.
- Esteves et al. (2013) Esteves, L. J., De Mooij, E. J. W., and Jayawardhana, R. (2013). Optical phase curves of kepler exoplanets. ApJ, 772(1), 51.
- Faigler and Mazeh (2007) Faigler, S. and Mazeh, T. (2007). Photometric detection of non-transiting short-period low-mass companions through the beaming, ellipsoidal, and reflection effects in Kepler and CoRoT lightcurves. Monthly Notices of the Royal Astronomical Society, 1-19.
- Faigler and Mazeh (2015) Faigler, S. and Mazeh, T. (2015). BEER analysis of Kepler and CoRoT lightcurves. II. evidence for superrotation in the phase curves of three Kepler hot jupiters. The Astrophysical Journal, 800(1).
- Faigler et al. (2013) Faigler, S., Tal-Or, L., Mazeh, T., Latham, D. W., and Buchhave, L. A. (2013). BEER analysis of Kepler and CoRoT light curves. i. discovery of Kepler-76b: A hot Jupiter with evidence for superrotation. ApJ, 771(1), 26.
- Feroz et al. (2009) Feroz, F., Hobson, M. P., and Bridges, M. (2009). MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics. MNRAS, 398(4), 1601–1614.
- Feroz et al. (2011) Feroz, F., Balan, S. T., and Hobson, M. P. (2011). Detecting extrasolar planets from stellar radial velocities using Bayesian evidence. MNRAS, 415(4), 3462–3472.
- Feroz et al. (2013) Feroz, F., Hobson, M. P., Cameron, E., and Pettitt, A. N. (2013). Importance nested sampling and the multinest algorithm. arXiv preprint arXiv:1306.2144.
- Ford (2006) Ford, E. B. (2006). Improving the efficiency of Markov Chain Monte Carlo for analyzing the orbits of extrasolar planets. ApJ, 642(1), 505.
- Ford et al. (2001) Ford, E. B., Seager, S., and Turner, E. L. (2001). Characterization of extrasolar terrestrial planets from diurnal photometric variability. arXiv preprint astro-ph/0109054.
- Gai (2016) Gai, A. (2016). Bayesian Model Testing of Models for Ellipsoidal Variation on Stars Due to Hot Jupiters. Master’s thesis, University at Albany.
- Gai and Knuth (2017) Gai, A. and Knuth, K. H. (2017). Bayesian model testing of ellipsoidal variations on stars due to hot Jupiters. In Prep.
- Gillon et al. (2013) Gillon, M., Triaud, A. H. M. J., Jehin, E., Delrez, L., Opitom, C., Magain, P., Lendl, M., and Queloz, D. (2013). Fast-evolving weather for the coolest of our two new substellar neighbours. Astronomy & Astrophysics, 555, L5.
- Gregersen (2010) Gregersen, E. (2010). The Milky Way and Beyond: Stars, Nebulae, and Other Galaxies. The Rosen Publishing Group.
- Hastings (1970) Hastings, W. K. (1970). Monte Carlo sampling methods using Markov chains and their applications. Biometrika, 57(1), 97109.
- Heller and Pudritz (2016) Heller, R. and Pudritz, R. E. (2016). The search for extraterrestrial intelligence in Earth’s solar transit zone. Astrobiology, 16(4), 259–270.
- Heng (2012) Heng, K. (2012). On the existence of shocks in irradiated exoplanetary atmospheres. ApJ Lett, 761(1), L1.
- Heng and Demory (2013) Heng, K. and Demory, B.-O. (2013). Understanding trends associated with clouds in irradiated exoplanets. ApJ, 777(2), 100.
- Heyrovský (2007) Heyrovský, D. (2007). Computing limb-darkening coefficients from stellar atmosphere models. The Astrophysical Journal, 656(1), 483–492.
- Hippke and Angerhausen (2015) Hippke, M. and Angerhausen, D. (2015). Photometry’s bright future: Detecting solar system analogs with future space telescopes. ApJ, 810(1), 29.
- Jackson et al. (2012) Jackson, B. K., Lewis, N., Barnes, J., Deming, L., Showman, A., and Fortney, J. (2012). The EVIL-MC model for ellipsoidal variations of planet-hosting stars and applications to the HAT-P-7 system. The Astrophysical Journal, 751:112 (13pp).
- Jaynes (2003) Jaynes, E. T. (2003). Probability Theory: The Logic of Science. Cambridge Univ. Press, Cambridge.
- Jenkins and Doyle (2003) Jenkins, J. M. and Doyle, L. R. (2003). Detecting reflected light from close-in extrasolar giant planets with the Kepler photometer. The Astrophysical Journal, 595(1), 429–445.
- Kane and Gelino (2012) Kane, S. R. and Gelino, D. M. (2012). Distinguishing between stellar and planetary companions with phase monitoring. MNRAS, 424(1), 779788.
- Knuth et al. (2015) Knuth, K. H., Habeck, M., Malakar, N. K., Mubeen, A. M., and Placek, B. (2015). Bayesian evidence and model selection. Digital Signal Processing, 47, 50–67. (arXiv:1411.3013).
- Ḱopal (1954) Ḱopal, Z. (1954). Photometric effects of reflection in close binary systems. MNRAS, 114(1), 101–117.
- Ḱopal (1959) Ḱopal, Z. (1959). Close Binary Systems. New York John Wiley & Sons Inc.
- Léger et al. (2011) Léger, A., Grasset, O., Fegley, B., Codron, F., Albarede, A. F., Barge, P., Barnes, R., Cance, P., Carpy, S., Catalano, F., et al. (2011). The extreme physical properties of the CoRoT-7b super-Earth. Icarus, 213(1), 1–11.
- Lock and Stewart (2017) Lock, S. J. and Stewart, S. T. (2017). The structure of terrestrial bodies: Impact heating, corotation limits, and synestias. Journal of Geophysical Research: Planets, 122(5), 950–982. (arXiv:1705.07858 [astro-ph.EP]).
- Mandel and Agol (2002) Mandel, K. and Agol, E. (2002). Analytic light curves for planetary transit searches. ApJ Lett, 580(2), L171.
- Maurin et al. (2012) Maurin, A. S., Selsis, F., Hersant, F., and Belu, A. (2012). Thermal phase curves of nontransiting terrestrial exoplanets-ii. characterizing airless planets. Astronomy & Astrophysics, 538, A95.
- Menou (2012) Menou, K. (2012). Magnetic scaling laws for the atmospheres of hot giant exoplanets. ApJ, 745(2), 138.
- Metropolis et al. (1953) Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., and Teller, E. (1953). Equations of state calculations by fast computing machines. Journal of Chemical Physics, 21(6), 1087–1092.
- Misra (2014) Misra, A. (2014). The Effects of Refraction and Forward Scattering on Exoplanet Transit Transmission Spectroscopy. Ph.D. thesis, University of Washington.
- Misra et al. (2014) Misra, A., Meadows, V., and Crisp, D. (2014). The effects of refraction on transit transmission spectroscopy: application to Earth-like exoplanets. The Astrophysical Journal, 792(1), 61.
- Morris (1985) Morris, S. L. (1985). The ellipsoidal variable stars. The Astrophysical Journal, 295, 143–152.
- Muñoz and Cabrera (2017) Muñoz, A. G. and Cabrera, J. (2017). Exoplanet phase curves at large phase angles. Diagnostics for extended hazy atmospheres. MNRAS. In Press.
- Muñoz et al. (2012) Muñoz, A. G., Osorio, M. R. Z., Barrena, R., Montañés-Rodríguez, P., Martín, E. L., and Pallé, E. (2012). Glancing views of the Earth: from a lunar eclipse to an exoplanetary transit. ApJ, 755(2), 103.
- Muñoz et al. (2017) Muñoz, A. G., Lavvas, P., and West, R. A. (2017). Titan brighter at twilight than in daylight. Nature Astronomy, 1, 0114. arXiv:1704.07460 [astro-ph.EP].
- NASA (2014) NASA (2014). NASA Strategic Plan 2014. NASA.
- NASA (2016) NASA (2016). Wide field infrared survey telescope. https://wfirst.gsfc.nasa.gov/index.html. Accessed: 2016-10-29.
- Pallé et al. (2008) Pallé, E., Ford, E. B., Seager, ., Montañés-Rodríguez, P., and Vazquez, M. (2008). Identifying the rotation rate and the presence of dynamic weather on extrasolar Earth-like planets from photometric observations. ApJ, 676(2), 1319.
- Perryman (2011) Perryman, M. (2011). The exoplanet handbook. Cambridge University Press.
- Placek (2014) Placek, B. (2014). Bayesian detection and characterization of extra-solar planets via photometric variations. Ph.D. thesis, University at Albany (SUNY).
- Placek et al. (2014) Placek, B., Knuth, K. H., and Angerhausen, D. (2014). EXONEST: Bayesian model selection applied to the detection and characterization of exoplanets via photometric variations. Astrophys. J., 795(2), 112. arXiv:1310.6764 [astro-ph.EP].
- Placek et al. (2015) Placek, B., Knuth, K., Angerhausen, D., and Jenkins, J. (2015). Characterization of Kepler-91b and the investigation of a potential trojan companion using EXONEST. The Astrophysical Journal, 814(2), 147.
- Placek et al. (2016) Placek, B., Knuth, K. H., and Angerhausen, D. (2016). Combining photometry from Kepler and TESS to improve short-period exoplanet characterization. PASP, 128(7), 074503.
- Placek et al. (2017) Placek, B., Angerhausen, D., and Knuth, K. H. (2017). Optimizing observing strategies for exoplanet secondary eclipses and phase curves. ApJ. In Press.
- Rauscher and Menou (2013) Rauscher, E. and Menou, K. (2013). Three-dimensional atmospheric circulation models of HD 189733b and HD 209458b with consistent magnetic drag and ohmic dissipation. ApJ, 764, 103.
- Ricker et al. (2010) Ricker, G. R., Latham, D. W., Vanderspek, R. K., Ennico, K. A., Bakos, G., Brown, T. M., Burgasser, A. J., Charbonneau, D., Clampin, M., Deming, L. D., et al. (2010). Transiting exoplanet survey satellite (tess). In Bulletin of the American Astronomical Society, volume 42, page 459.
- Rybicki and Lightman (2008) Rybicki, G. B. and Lightman, A. P. (2008). Radiative processes in astrophysics. John Wiley & Sons.
- Seager (2010) Seager, S. (2010). Exoplanet atmospheres: physical processes. Princeton University Press.
- Seager and Mallěn-Ornelas (2003) Seager, S. and Mallěn-Ornelas, G. (2003). A unique solution of planet and star parameters from an extrasolar planet transit light curve. ApJ, 585(2), 1038.
- Showman and Polvani (2011) Showman, A. P. and Polvani, L. M. (2011). Equatorial superrotation on tidally locked exoplanets. ApJ, 738(1), 71.
- Sidis (2010) Sidis, O. (2010). Transits of transparent planetsatmospheric lensing effects. The Astrophysical Journal, 720(1), 904.
- Sing (2010) Sing, D. (2010). Stellar limb-darkening coefficients for CoRoT and Kepler. Astronomy & Astrophysics, 510(1), A21.
- Sivia and Skilling (2006) Sivia, D. S. and Skilling, J. (2006). Data Analysis. A Bayesian Tutorial. Oxford University Press, Oxford, second edition.
- Skilling (2004) Skilling, J. (2004). Nested sampling. In R. Fischer, V. Dose, R. Preuss, and U. von Toussaint, editors, Bayesian Inference and Maximum Entropy Methods in Science and Engineering, Garching, Germany 2004, number 735 in AIP Conf. Proc., pages 395–405. AIP, New York.
- Skilling (2006) Skilling, J. (2006). Nested sampling for general Bayesian computation. Bayesian Analysis, 1(4), 833–859.
- Sliski and Kipping (2014) Sliski, D. H. and Kipping, D. M. (2014). A high false positive rate for Kepler planetary candidates of giant stars using asterodensity profiling. ApJ, 788(2), 148.
- Sobolev (1975) Sobolev, V. V. (1975). Light Scattering in Planetary Atmospheres. Oxford: Pergamon.
- Van Cleve and Caldwell (2009) Van Cleve, J. E. and Caldwell, D. A. (2009). Kepler Instrument Handbook. NASA/Ames Research Center, Moffett Field, California.
- Van Eylen et al. (2013) Van Eylen, V., Nielsen, M. L., Hinrup, B., Tingley, B., and Kjeldsen, H. (2013). Investigation of systematic effects in kepler data: Seasonal variations in the light curve of HAT-P-7b. ApJ Lett, 774(2), L19.