Improving Inferences about Exoplanet Habitability Improving Inferences about Exoplanet Habitability

: Assessing the habitability of exoplanets (planets orbiting other stars) is of great importance in deciding which planets warrant further careful study. Planets in the habitable zones of stars like our Sun are sufﬁciently far away from the star so that the light rays from the star can be assumed to be parallel, leading to straightforward analytic models for stellar illumination of the planet’s surface. However, for planets in the close-in habitable zones of dim red dwarf stars, such as the potentially habitable planet orbiting our nearest stellar neighbor, Proxima Centauri, the analytic illumination models based on the parallel ray approximation do not hold, resulting in severe biases in the estimates of the planetary characteristics, thus impacting efforts to understand the planet’s atmosphere and climate. In this paper, we present our efforts to improve the instellation (stellar illumination) models for close-in orbiting planets and the signiﬁcance of the implementation of these improved models into EXONEST, which is a Bayesian machine learning application for exoplanet characterization. The ultimate goal is to use these improved models and parameter estimates to model the climates of close-in orbiting exoplanets using planetary General Circulation Models (GCM). More speciﬁcally, we are working to apply our instellation corrections to the NASA ROCKE-3D GCM to study the climates of the potentially habitable planets in the Trappist-1 system.


Introduction
Typically, the illumination of exoplanets is determined under the assumption that the incident stellar radiation reaches the exoplanet as plane parallel rays.In this situation, half of the planet is fully illuminated and the other half is dark.For extremely close-in exoplanets for which the star-planet separation is on the order of only a few stellar radii, it is particularly important to take into account the finite angular size of the host star.Doing so reveals that two additional zones are necessary to describe the illumination of the exoplanet.Previous work has shown that the planet will have four distinct illumination zones: the fully-illuminated zone, the first penumbral zone, the second penumbral zone, and the un-illuminated (or night) zone [1] (Figure 1).
The fully-illuminated zone (FZ) is the region on the planet in which the entire apparent disk of the star is visible in the sky.However, due to the close-in geometry, the surface area of the visible star is less than one-half of the total area.The first penumbral zone (PZ 1 ) is the region on the planet in which a portion of the star (less than one-half of the potentially visible area) is occluded by the horizon; that is, the star is on the verge of setting.The second penumbral zone (PZ 2 ) is the region on the planet in which more than one-half of the potentially visible area of the star is occluded by the horizon and the un-illuminated zone (UZ) is the region on the planet from which the star cannot be seen.
The result is that the rather abrupt terminator between the day-side and night-side of a distant planet is, in the case of a close-in planet, spread out into the two penumbral zones, which impinge on both the fully illuminated and the un-illuminated zones.So while a smaller surface area of the close-in planet is fully-illuminated, more of the planet's surface receives light because the un-illuminated zone is also smaller.For example, in the case of the close-in orbiting hot Jupiter, Kepler-91b, about 30% of the surface is in full daylight, 40% of the surface is in partial penumbral light, and 30% of the surface is in full darkness, so that approximately 70% of the planet is receiving some light at any point in time [2].These differences in illumination are expected to impact the planet's climate.The main goal of this paper is to construct analytic expressions for the flux incident on an exoplanet to be implemented in EXONEST, a Bayesian machine learning application developed for exoplanet characterization [2,3].

Methods
In this section we build analytic expressions for the light received by a point on the exoplanet taking into account the finite angular size of the host star.

Instellation Setup
Suppose that the star emits P photons per unit time.Integrating the photons over the blackbody spectrum will give the amount of emitted energy, and therefore this P is related to the power emitted by the star.If each area element of the star, dσ S , emits dp photons per unit time, the total number of photons emitted by the star in unit time, P, is given by: If p is the number of photons emitted by an area element in a unit time, then and Let the radius of the star be R S .Assuming the star to be a sphere, we can write dσ S in terms of the surface area element as Integrating over the whole star yields which allows us to write Some of these photons will hit the planet.
Next, let us consider photons emitted from the area element on the stellar surface at point A, and determine how many of these photons will be incident at a point B on the planet.Let d f A→B denote the number of photons emitted by A and arriving at B in a unit time.Then, from Figure 2: Substituting ρ = ρ ρ and dp using Equations ( 2)-( 6) gives:

Coordinate Systems
When we have a star-planet dynamical system, it is convenient to introduce three coordinate systems: stellar, planetocentric, and overall.For the scope of this paper, it suffices to use just two, stellar and overall, by choosing the origin of the overall coordinate system at the center of the exoplanet such that it would coincide with the planetocentric coordinate system (Figure 3).Furthermore, the axial symmetry of stellar illumination about the axis joining the centers of the planet and the star can be exploited by orienting the z-axis of the overall coordinate system in this direction.With this coordinate system, the problem simplifies to having to determine the illumination along one semicircular arc.
The stellar and overall Cartesian coordinate systems are related to one another as follows: Now that we have chosen appropriate coordinate systems, next let us derive the vectors required for the illumination in (8).O O , respectively.The stellar coordinate system, which employs spherical coordinates (θ, ω), is used to define points on the star, while the overall coordinate system employs spherical coordinates (η, φ).

The Vector ρ
The vector ρ is given by where in the overall frame The sum then gives which has a magnitude of We can now verify that if the point B is located at (η, φ) = (0, π/2), then as expected.

The Vector ρ
Let us consider the triangle O S AB in Figure 2.

− →
The vector − − → O S A in the stellar coordinate system is given by Constructing the vector − − → O S A in the overall coordinate system requires some care.We can use the fact that where by the relation in Equation (10), and such that which simplifies to The squared-magnitude of ρ is The distance ρ is then found by taking the square root of the result above.

The Unit Vectors
Now let us obtain the relevant unit vectors, m and n, in the overall coordinate system.Note that m is the unit vector in the − − → With an understanding of how to construct the relevant vectors and express them in the overall coordinate system, we can now write out the necessary integrals to obtain the total instellation.

Total Instellation
To compute the total instellation, we must integrate over the visible surface of the star, which is described by the spherical coordinates θ and ω.The total instellation at a point B on the planet is given by integrating over all possible points, A, on the star that is visible to point B on the planet: where θ 1 , θ 2 , ω 1 and ω 2 are the limits of integration that depend on the zone in which point B is located.Upon substitution of ρ , its magnitude and the unit vectors m and n, into the total instellation in (31), it simplifies to 5 distinct double integrals with respect to θ and ω, which must be integrated over the limits of each zone to obtain an analytic model.

Discussion
In this section, we discuss the applications and implications of improved instellation models to illustrate the benefits of having accurate models for stellar illumination in order to determine the habitability of close-in exoplanets.

Bayesian Exoplanet Characterization and Model Testing
The EXONEST software suite performs Bayesian model-based exoplanet characterization using photometric and radial velocity data [2][3][4][5].EXONEST relies on predictive photometric models that generate a predicted stellar light curve based on a set of hypothesized model parameters.Bayesian probability theory [6][7][8][9] uses this predicted stellar light curve to compute the posterior probability that the set of hypothetical model parameters could have been responsible for the recorded light curve.This allows detailed astrophysical models to be evaluated by the data analysis algorithm, as well as statistical testing among competing models [3,8,10,11].EXONEST allows the user to employ one of several variants of the Nested Sampling Algorithm, such as Nested Sampling [8,12], Multi-Nest [13,14], and the state-of-the-art Nested Sampling with Galilean Monte Carlo [15][16][17].These Nested Sampling algorithms provide estimates of the model parameter values as well as the Bayesian evidence, which allows for statistical model testing.Using evidencebased Bayesian techniques [10], exoplanetary models can be ranked based on the Bayesian evidence providing the capability for scientific hypothesis testing.One can test various photometric effects and compute the evidence to determine the probability that these effects are present in any given system.Relying on data from the NASA Exoplanet Archive, we can utilize EXONEST and its ability to perform Bayesian model testing to evaluate the efficacy of our improved instellation models by computing the Bayesian evidence of the parallel-ray approximation models versus our proposed models.
In any problem involving model parameter estimation, there are legitimate concerns about potential degeneracies among the model parameters, especially as the number of parameters increases.Degeneracies can be expected when more than one parameter affects the predicted light curve in a similar way.In reality, this is a warning that additional, different data are needed to effectively estimate the model parameters.Examining the cross-correlations between the parameter estimates is a very effective way to identify degeneracies.Model testing, as described above, can be used to detect degeneracies as well, since a degenerate set of model parameters will not significantly increase the likelihood function, but it will significantly decrease the prior probability, resulting in lower evidence.

Implications for Potentially Habitable M-Dwarf Worlds
The exoplanet TRAPPIST-1d (R star = 0.00056 AU = 0.121R SUN , a = 0.022 AU = 39R star ), which closely orbits its host M-dwarf star in the habitable zone (HZ), has a penumbral zone that covers about 2.5% of its surface or 8.25 million km 2 .While this is a small percentage of the surface of the planet, this is greater than the area of Central and South America combined, which contains several important biogeographical regions on Earth.(e.g., [18][19][20][21]).Our estimates reveal that planets orbiting M-dwarf stars on the inner edge of the HZ will have penumbral zones ranging in area from around 1 to 3% of the surface area of the planet.While this may seem like a small value, it may have its greatest climatological impact at the poles of low obliquity worlds where the ice-albedo feedback manifests itself consistently and tends to drive mean surface temperatures lower, which can affect the fractional habitability metrics like that proposed in [22].A detailed understanding of the incident stellar radiation in these penumbral regions is expected to have significant relevance to understanding habitability on these worlds.

Climate Modeling of Close-In Orbiting Planets
As mentioned before, corrections are required to the illumination models used in planetary General Circulation Models (GCM) in order to properly model the climate of a close-in exoplanet.The analytical expressions for the incident stellar flux can be applied to the illumination model in the ROCKE-3D GCM [23].Not only have such corrections never been applied in a GCM, but their effects have never been quantified.For example, differences in illumination at the poles of zero obliquity worlds due to these corrected illumination models could have large effects on ice coverage at the poles.Recent work [24] has shown that polar supercontinents and ocean circulation can affect global mean surface temperatures by several degrees because of the extent of polar ice and the resulting icealbedo feedback.We expect corrected illumination models to have similar effects.Using the modified ROCKE-3D code for close-in planets, we could study the climates of potentially habitable planets.Critically, we can make one-to-one comparisons of GCM results between penumbral and parallel ray (the current default) illuminations and derive statistical measures of difference (e.g., systematic latitudinal-and global-surface temperature differences, ice-cover, precipitation patterns).We could look at a variety of topographies, land/sea masks, and surface water content used in previous ROCKE-3D publications such as Proxima Centauri b [25], Future Earth [24], and ancient Venus [26].

Conclusions
Developing accurate analytical stellar instellation models for close-in orbiting exoplanets will increase both the accuracy and precision of inferred model parameter estimates obtained from data, thereby improving exoplanet characterization and potentially the detection of habitable exoplanets.This is a multidisciplinary project, founded on astrophysics with machine learning.Studying the implications of these improved illumination models in the context of exoplanetary climate modeling will provide a better understanding of the climates and atmospheric physics of close-in orbiting exoplanets.

Figure 1 .
Figure 1.Extremely close-in exoplanets exhibit four zones of illumination: a fully illuminated zone (FZ) defined by the inner tangents, two penumbral zones (PZ 1 , PZ 2 ) defined by whether more than half or less than half of the star is visible, and a dark un-illuminated zone (UZ) defined by the outer tangents.

Figure 2 .
Figure 2. Star-planet geometry involving light emission at point A on the star and incident on the planet at point B: The star with radius R S is centered at O S and the exoplanet with radius R P is centered at O O .The star-planet separation is denoted by the vector r.

Figure 3 .
Figure 3. Orientations of the stellar and overall coordinate systems with origins located at O S andO O , respectively.The stellar coordinate system, which employs spherical coordinates (θ, ω), is used to define points on the star, while the overall coordinate system employs spherical coordinates (η, φ).
and − − → O P B is the position vector of the point B in the overall frame: − − → O P B = R P sin η cos φ xo + R P sin η sin φ ŷo + R P cos η ẑo .