Issue 
A&A
Volume 578, June 2015



Article Number  A59  
Number of page(s)  9  
Section  Numerical methods and codes  
DOI  https://doi.org/10.1051/00046361/201424202  
Published online  04 June 2015 
Predicting exoplanet observability in time, contrast, separation, and polarization, in scattered light
^{1} LESIA, Observatoire de Paris, CNRS/UMR 8109, UPMC, Université Paris Diderot, 5 place J. Janssen, 92195 Meudon, France
email: guillaume.schworer@obspm.fr
^{2} Sydney Institute for Astronomy (SIfA), School of Physics, The University of Sydney, NSW 2006, Australia
Received: 14 May 2014
Accepted: 14 March 2015
Context. Polarimetry is one of the keys to enhanced direct imaging of exoplanets. Not only does it deliver a differential observable providing extra contrast, but when coupled with spectroscopy, it also reveals valuable information on the exoplanetary atmospheric composition. Nevertheless, angular separation and contrast ratio to the hoststar make for extremely challenging observation. Producing detailed predictions for exactly how the expected signals should appear is of critical importance for the designs and observational strategies of tomorrow’s telescopes.
Aims. We aim at accurately determining the magnitudes and evolution of the main observational signatures for imaging an exoplanet: separation, contrast ratio to the hoststar and polarization as a function of the orbital geometry and the reflectance parameters of the exoplanet.
Methods. These parameters were used to construct a polarizedreflectance model based on the input of orbital parameters and two albedo values. The model is able to calculate a variety of observational predictions for exoplanets at any orbital time.
Results. The interdependency of the three main observational criteria – angular separation, contrast ratio, polarization – result in a complex timeevolution of the system. They greatly affect the viability of planet observation by direct imaging. We introduce a new generic display of the main observational criteria, which enables an observer to determine whether an exoplanet is within detection limits: the SeparationPOlarizationContrast (SPOC) diagrams.
Conclusions. We explore the complex effect of orbital and albedo parameters on the visibility of an exoplanet. The code we developed is available for public use and collaborative improvement on the python package index, together with its documentation. It is another step towards a full comprehensive simulation tool for predicting and interpreting the results of future observational exoplanetary discovery campaigns.
Key words: planets and satellites: atmospheres / planets and satellites: detection / methods: observational
© ESO, 2015
1. Introduction
Observational technologies for detecting the light reflected from an exoplanet are reaching a level of precision that makes direct imaging of exoplanets a realistic possibility in about the coming decade. Because of the very challenging contrast ratio between exoplanets and hoststars in the optical, there is great interest in differential methods, such as polarimetry, to deliver an extended reach to imaging instruments (Hough et al. 2006; Keller 2006; Schmid et al. 2006). Built upon previous detections of polarized signals from exoplanets Berdyugina et al. (2008), Wiktorowicz (2009), Lucas et al. (2009), Berdyugina et al. (2011), the polarimetric imaging modules that are currently being integrated on the 10m class telescopes (SPHEREZIMPOL on the VLT Thalmann et al. (2008), GPI on Gemini Wiktorowicz et al. (2012), and VAMPIRES on Subaru Telescope Norris et al. 2015) promise to achieve a differential polarized contrast down to ≈ 10^{6} in the visible or near infrared. It is therefore an opportune moment to produce detailed predictions for exactly how the expected signals should appear, which will be of critical importance for the designs and observational strategies of these instruments. We aim at accurately determining the magnitudes and evolution of the main observational signatures as a function of the basic parameters of the exoplanetary system: the starexoplanet orbital parameters and the optical properties of the planet. We incorporate the complexity arising from effects such as polarization based on Rayleigh scattering. Additionally and as a firstorder observability estimator, the absolute flux in photon per unit time and surface received from the exoplanet target is computed assuming a blackbodied star. We finally provide a querier and parser of the http://exoplanet.eu/ exoplanet database for searching and importing any starplanet target.
Several models were already developed to predict the light signature of an exoplanet as it should appear to an observer, for example see Cahoy et al. (2010), Buenzli & Schmid (2009), or for the case of Earthlike planets with varying atmospheric parameters, see Karalidi & Stam (2012) and Zugger et al. (2010). However, to the best of our knowledge, the interdependency of the main observational criteria is never taken carefully into account. The relative evolution of angular separation, polarization, and contrast ratio is of critical importance because their respective maxima do not occur at the same time in the general case: this makes the bestcase scenario for detection very unlikely. We correct for this by showing the relative evolution of the angular separation, contrast ratio, and polarization as a function of the orbital and reflectance parameters of the planet. The findings resulting from this integrated treatment highlight dependencies that are much more complex than previously reported. The three main observational criteria – angular separation, contrast ratio, and polarization – in general do not exhibit a strong positive correlation and must be analyzed separately to determine the direct visibility of a potential planet target.
2. Model
We developed a polarized reflectance model. Its algorithm relies on the orbital parameters of the exoplanet and the linear combination of two albedo values; it is briefly described below.
As our prime observables are in any event differential between the star and exoplanet (such as the contrast ratio) or inherent to the planet (such as the polarization fraction), absolute fluxes are not described in detail. A simple blackbody modeling of the host stars is used to give a firstorder idea on the absolute fluxes that are expected from a exoplanet. For the sake of conciseness, “polarization” is used instead of “polarization fraction”.
2.1. Orbital parameters
Each exoplanet orbit is defined by the usual parameters listed in Table 1. Note that the shape of the orbit as seen by the observer is solely defined by e, i, and ω_{p}.
Orbital parameters for an exoplanet.
The phase angle of the exoplanet and the distance between the star and the planet are obtained from these orbital parameters. The phase angle α is defined as the vector angle between the star, an exoplanet, and the observer, as seen from the exoplanet. It is 0° when the planet is at full phase (superior conjunction) and 180° when the exoplanet is at new phase (inferior conjunction, transiting).
2.2. Rayleigh scattering polarization
For the sake of conciseness, in the following discussions we use the quantity called phased albedo, which is defined as the product of the phase function and the geometric albedo of a planet, for a given wavelength: (1)for an incident flux πF, with I(α) the emerging flux from a body at the phase angle α. It represents the fraction of the incident irradiance that is reflected by the planet when it is seen at phase angle α by the observer, so that E_{Planet, out} = A_{α}(α)·E_{Planet, in}, with E being the spectral irradiance arriving at or leaving the planet.
The polarization induced by an exoplanet on a reflected beam of light is described by many different models, such as Mie, Rayleigh, or Raman scattering. It is assumed here that single Rayleigh scattering is the predominant source of polarization in planetary atmospheres. This is a good approximation as long as Mie scattering (and especially the primary rainbow polarizing effect of liquid droplets in clouds, see Bailey 2007) remains lower than Rayleigh scattering. This scenario is usually achieved for wavelengths shorter than approximately 1 μm, which also corresponds to the most favorable wavelengths for Rayleigh polarization measurements of scattered light, as it is function of λ^{4}. Integrated over the stellar disk, the flux from the star can be considered to be unpolarized (Kemp et al. 1987). Hence, only the light from the exoplanet carries polarization.
The results of Madhusudhan & Burrows (2012) for the Stokes parameters Q/I and U/I were adopted as input data to carry out more polarization and reflectance computations. Their model assumes an unresolved semiinfinite homogeneous atmosphere (hence cloud free), dominated by Rayleigh scattering. The atmosphere is assumed to be spherical, which is a satisfactory assumption when the planetary rotation remains slow. Given the strong depolarizing effect of multiple scattering, they used the singlescattering albedo as a unique reflectance parameter from which a geometric albedo A_{g,Rayleigh} arising from the Rayleigh singlescattering albedo was calculated; their result is reproduced in Fig. 1.
Fig. 1 Geometric albedo A_{g,Rayleigh} arising from singlescattering as a function of the Rayleigh singlescattering albedo w_{R}. 
StokesV is found to be zero for Rayleigh scattering in planet atmospheres. Hence, the degree of polarization Pz_{Rayleigh} is defined by (2)and it depends on the phase angle α of the planet and the singlescattering albedo w_{R}, which represents the amount of absorption versus scattering in a given atmosphere. The existing literature provides solutions to several phase functions. Madhusudhan & Burrows (2012) provided an analytical solution for scattering models for Lambertian, Rayleigh, isotropic, and asymmetric scattering. Zugger et al. (2010) addressed the liquid surface scattering model. Some phase functions are reproduced in Fig. 2. Given the symmetry of the phase angle α along the orbit, only α = [0,180] ° values are represented. The illuminated fraction values are also plotted for comparison.
Fig. 2 Phase functions φ as a function of phase angle α for several scattering models. The illuminated fraction is also displayed (thick line). The Rayleighscattering model corresponds to a vectorial Rayleigh phase matrix. 
Rayleigh and Lambertian surface models follow the same general Scurve shape between zero flux and full flux. A Lambertian or Rayleigh planet is faint at small α because the illuminated fraction f_{I} is small, hence fewer photons are reflected.
2.3. Polarized and unpolarized albedo
It is now apparent that clouds and condensates are very common in planet atmospheres. An atmosphere solely described by a singlescattering albedo could only take into account Rayleigh scattering of the atmosphere in the optically thick case, which is quite restrictive. Indeed, clouds, rocky surfaces or liquid surfaces do not follow Rayleigh scattering, but they often strongly contribute to the reflectivity of the exoplanet.
In this section we describe a novel approach to the reflectance properties of an exoplanet atmosphere to include effects of a greater atmospheric variety.
We assumed that any nonRayleigh singlescattering from the atmosphere follows Lambertian scattering. This assumption is discussed below. A planet is better defined by two combined albedo quantities:

a total geometric albedo A_{g,Total} that represents the total fraction ofincoming irradiance reflected by the exoplanet, taking intoaccount all scattering processes,

a Rayleigh singlescattering albedo w_{R} that represents the fraction of incoming irradiance that is scattered by the atmosphere according to Rayleigh singlescattering alone, from which arises a Rayleigh geometric albedo A_{g,Rayleigh}. Only this part of reflected light carries a defined degree of polarization in this model. This albedo can be approximated by the singlescattering albedo of the exoplanet multiplied by the fraction of Rayleighscattered emerging light.
These albedo values can be added linearly so that (3)Following this, we can express the outgoing irradiances using the phased albedo A_{α}: The total emerging spectral irradiance of the exoplanet I_{Total,out} including thermal emission is (6)Note that this thermal emission is much weaker than the reflected light from the hoststar unless the exoplanet is a hot Jupiter and the observing wavelength lies in the infrared. Below, we assume that E_{Thermal,out} is negligible for λ_{Obs} ≲ 1 μm.
The different Stokes parameters and polarization degree can be easily obtained with where Q/I_{Rayleigh,out} and U/I_{Rayleigh,out} are the Stokes ratios from the Rayleighscattering polarization section of this model; they are used to calculate Pz_{Raylgeigh} using Eq. (2).
According to previous sections, the following complete relation for the exoplanet reflected irradiance, given the hoststar surface irradiance, is (11)where R are radii. The observerexoplanet distance is here approximated by the distance observerstar d. The contrast ratio CR between a planet and its hoststar is then obtained from Eq. (11): (12)We highlight here that the polarization of the exoplanet Pz_{Total} and the contrast ratio CR do not depend on the incoming irradiance from the hoststar, they are intrinsic to the planet. However, they are strongly dependent on the observation wavelength λ_{Obs} through Rayleigh singlescattering and geometric albedo, even though for clarity of notation, λ was not explicitly written in the relations.
Thanks to the separation of the contributions from the different scattering processes on an exoplanet, an original model for planetary polarized reflectance has now been set up. It yields the phased albedo A_{α}, which indicates the fraction of reflected irradiance from an observed planet, as a function of geometric albedo A_{g}, phase angle α, and Rayleigh singlescattering albedo w_{R}. While α can be easily determined as described in the orbital computations part of the model, the albedo parameters A_{g} and w_{R} depend on factors too numerous to be modeled here: they therefore need prior computation. If A_{g} and w_{R} are obtained separately, the albedo and phase functions previously defined are successful in describing the emerging irradiance that is reflected by a planet. Its most accurate results are obtained for wavelengths shorter than approximately 1 μm unless thermal emission and Mie scattering are also added to the model.
We assumed a polarizedreflectance model that incorporates Rayleigh and Lambertian scattering as two linear scattering phenomena. The approximation made is that any scattering that does not relate to Rayleigh singlescattering is assumed to be Lambertian. As a consequence, a beam of light undergoing multipleRayleigh scattering or a combination of Lambertian and Rayleigh scattering is assumed to be multiple Lambertian scattered. This approximation mainly affects the angle at which a photon emerges from the exoplanet and in turn the phasefunction of integrated emerging light. In our model, this integrated phase function is a linear combination of Rayleigh and Lambertian phase functions. Adding these secondorder scattering phenomena to the computation of the integrated phase function adds terms to its computation: a multiple Rayleigh phase function term and a coupled RayleighLambertian phase function term. We assumed that these secondorder phase functions are similar to the Rayleigh and Lambertian phase functions, which differ only slightly (see Fig. 2). Furthermore, we assumed that the weight of these additional phase functions in the computation of the integrated phase function is smaller than Rayleigh singlescattering and Lambertian single or multiple scatterings. Note also that the lower the albedo values, the less likely multiplescattering becomes and the smaller the errors induced by this approximation.
The results were benchedmarked against several other models providing photopolarimetric curves as a function of one or several orbital parameters, such as models developed by Buenzli & Schmid (2009), Madhusudhan & Burrows (2012), Fluri & Berdyugina (2010), and Zugger et al. (2010). The locations of the polarization peaks and minima were reproduced with very good agreement; they mostly depend on the phase functions and the processing of orbital parameters. The intensity of the polarization peaks were reproduced with good accuracy; they mostly depended on the finetuning of the two albedo values. More specifically, the shift of the polarization peak to phase angles greater than 90° reported in the last benchmark reference was also observed.
2.4. Model completeness
The calculation code that implements this model takes up to 11 input parameters and N orbital positions for which the calculations are performed. It generates Nelement vectors for up to 11 output parameters, listed in Table 2. Table 3 shows the mapping between the input and output parameters, where the input parameters are listed as column headers and output parameters are row headers of Table 3.
Output parameters from the calculation code.
Mapping between the input (columns) and output (lines) parameters of the model.
In the way the code is written, the only mandatory input is the eccentricity, the inclination, and the argument at periapsis. When compared with exoplanetary detections reported in the literature, these three orbital parameters are rarely known with any precision. This is currently the case for 132 exoplanets out of the 1790 (≈7.4%, source: http://exoplanet.eu/). However, where more complete data exist, reasonably good estimates can be made.
3. Results
3.1. SPOC diagrams
An important tool introduced here is the SeparationPOlarizationContrast (SPOC) diagram which presents all the useful information to enable an observer to evaluate whether a planet target is observable at given instrumental limits and how this signal will evolve with orbit and time. Perhaps most importantly, in the event of a detection, SPOC diagrams provide a powerful mechanism to constrain exoplanet properties given observational imaging data.
Figure 3 shows an example of a SPOC diagram applied to Alpha Centaury Bb, an Earthmass planet discovered in 2012 (Dumusque et al. 2012) that was later debated (Hatzes 2013). Its probable orbital parameters are shown in Table 4. The exoplanet was discovered with the radial velocities method, hence its inclination and argument of ascending node are unknown. For i, values in [15, 35, 60, 85]° were explored, whereas Ω was neglected because it does not change the shape of the orbit: it only defines a rotation of the orbit locus as seen by an observer along the line of sight (i.e., it is only useful for projecting the exoplanet location around the star onto RADec axes). Its radius was calculated from its mass·sin(i) = 0.0036·M_{Jupiter} assuming the same density as Earth, which leads to M = [4.42,2.00,1.32,1.15]M_{Earth} and R = [1.64,1.26,1.10,1.05]R_{Earth} for the four previous i. In this example, the light from Alpha Centaury A received by the exoplanet is neglected because it is roughly 2 × 10^{5} fainter than that from Alpha Centaury B. Reflectance parameters were chosen to be similar to those of a Venuslike planet (A_{g} = 0.67, ω_{p} = 0.85) in Vband.
In Fig. 3, Cen_{α,0.34} is shown with i =15°. This figure displays the polarization encoded with color, while the two axes plot the separation against contrast ratio; finally, a set of different markers give temporal information. The first striking feature of the plot is that the exoplanet does not run a “back and forth” locus in the separationcontrast ratio phase diagram. This is because the apparent symmetry of the orbit is broken. When the orbital parameters i and e are nonzero and ω_{p} is different from [0,90,180,270]°, the geometry of the orbit as it appears to the observer has no symmetry. In other words, the star is not at a focal point of the apparent elliptical locus of the orbit. Note that the semimajor axis a does not play a role in the shape of the orbit; it acts as a simple scaling factor. The photon flux on the righthand side yaxis is given in photons per collectingarea per hour, integrated over Vband. It was calculated assuming that Alpha Centaury B is a blackbody, T_{eff} = 5214 and M_{v} = 1.33. The collecting area is 51.7 m^{2}, which corresponds to a diameter of 8.2 m with a central obstruction of 2% (same as VLT telescopes).
Fig. 3 SPOC diagram for Cen_{α,0.34} with i = 15°, see Table 4, assumed to be Venuslike regarding its reflectance parameters. It represents the polarization values as gray scale, plotted against the angular separation and the contrast ratio axes. The exoplanet’s orbit is therefore displayed in a separationcontrast phase diagram. The wavelength is the Vband. The “x” indicates the position of the exoplanet at an arbitrary observation time. The diamond indicates the position of the planet one sixth of the orbital period T later (giving the forward direction of time). The 10 plus signs indicate the evolution of the planet along the curve; they are linearly spaced in time. The empty and filled triangles are periapsis and apoapsis: r_{min} = 0.026 AU and r_{max} = 0.054 AU. The empty and filled circles are minimum and maximum phase angle: α_{min} = 75° and α_{max} = 105°. The bottom left panel shows the orbit as seen by the observer, where the shaded area corresponds to the portion of the orbit that is behind the plane of the paper. The dashed lined corresponds to the limit of the same instrument sensitivities as in the preceding example. The righthand side yaxis is given in photons per VLT collecting area per hour integrated over Vband. Two successive plus signs correspond to 7h46min. 
Two sets of probable orbital parameters for Alpha Centaury Bb, see Sect. 9 of the supplementary information of Dumusque et al. (2012).
Figure 3 shows that the lowest and highest contrast ratios are reached close to the maximum and minimum phase angle (black and white disks). This corresponds to the inferior and superior conjunction of the exoplanet with its star. We note that the periapsis and apoapsis are not reached at these phase angle extrema because of the nonnull value of ω_{p}, which reorients the orbit with respect to its host star. Polarization reaches two maxima near α = 90° (the limit between the shaded and the bright area in the bottom left subplot of the same figure). This is linked to the fact that polarization is observed from Rayleigh scattering, which has a peak polarization for 90°. However, in the case of nonentirely Rayleigh planets, the peak polarization occurs at slightly different phase angles.
An even more interesting diagram is shown in Fig. 4. It represents the apparent orbital motion in mas per hour plotted against the angular separation as xaxis and the polarized contrast ratio as yaxis. The apparent orbital motion is found to be an important factor to take into account here because the period of Alpha Centaury Bb is only 3.24 days: it reaches 1.8 to 3.5 mas per hour, which limits the observer to a few hours of exposures, depending on the plate scale of its detector. The photon flux is extremely high as a result of the brightness of the star (1.3 mag in V) and the large bandwidth of the V filter (550 ± 44 nm). The SPHEREZIMPOL performance curve for 1 h exposure time using a Vfilter with the doubledifference polarization calibration was added to the plot (see Fig. 4 in Roelfsema et al. 2014). It shows that this instrument could observe Alpha Centaury Bb during ≈35% of its orbit if i = 15°. Longer integrations could increase the sensitivity at the expense of larger apparent orbital motion on the detector, which will significantly limit the angular differential imaging (ADI) processing capability.
Fig. 4 Same caption as Fig. 3 except that it represents the apparent orbital motion in mas per hour as gray scale, plotted against the angular separation and the polarized contrast ratio axes. The red dotted line corresponds to the SPHEREZIMPOL performance for 1 h exposure (observability domain being above the line). 
Figure 5 shows SPOC diagrams for the two sets of orbital parameters Cen_{α,0.34} and Cen_{α,0}, and for all four i.
Fig. 5 Contrast ratio versus angular separation plots for both orbital parameters of Alpha Centaury Bb as described in Fig. 3 and for inclinations in [15, 35, 60, 85]°. Dotted lines correspond to e = 0 and ω_{p} = 0° and full lines to e = 0.34 and ω_{p} = 246°. A color version is available in the online journal. 
This figure shows that for higher orbital inclinations, the radius inferred from the orbital velocity profile decreases, which leads to a smaller reflecting surface and consequently, a lower contrast ratio. The effect of eccentricity is clearly seen: curves for Cen_{α,0} all run a “back and forth locus” because of their orbital symmetry. Eccentric orbits display an increasingly complex shape as the inclination increases: this is due to the projection of the orbit locus on the sky as an observer would see it.
We refer to Appendix A for more detailed studies of diverse SPOC diagrams that are based on known exoplanets.
3.2. Separations and contrast ratios for the solar system
An interesting illustration of the usefulness of the SPOC diagrams is provided by our solar system as viewed by an external observer.
Fig. 6 Contrastseparation diagram for the planets of the solar system as seen from 1pc distance in visible light. Small black curves correspond to an inclination of the ecliptic of 0° (observer reference frame); colored curves show the solar system at an inclination of 45°. Note that both axes have logarithmic scales. A color version is available online. 
Planetary and orbital parameters for the planets of the solar system, in the Vband.
Figure 6 gives the contrastseparation diagram for the planets of the solar system, whose parameters are listed in Table 5. In this figure, polarization is not shown (for clarity). It displays values for two different angles of inclination for the observer, 0° and 45°. The observational wavelength is visible light. Note that no rings were modeled for Saturn; they are anticipated to have a variable effect on both contrast ratio and polarization (see discussion in Sect. 3.3).
Note that the angular separation of Earth to the Sun is ≈1 arcsec when the ecliptic inclination is 0°, consistent with the definition of a parsec. Earth’s average angular separation decreases as the orbit appears more inclined to the observer.
Figure 6 highlights the strong dependency of the contrast ratio and angular separation on inclination, argument at periapsis, and eccentricity. Mercury especially, but also Mars and Saturn do not follow a “back and forth locus” because the symmetry of their orbit as it appears to the observer is broken.
Venus is seen to be the second brightest planet in the solar system (for a ringless Saturn). Significantly brighter than Earth, but fainter than Jupiter. A rule of thumb can explain this easily: Venus receives twice as much solar flux as Earth and has an albedo nearly twice as high. Because the diameters of the two planets are close to each other, they present almost the same reflective area (Venus has ≈90% of that of Earth): Venus is therefore on average nearly four times brighter than Earth to an external observer. The same rules can be applied to compare the relative brightness of Jupiter and Venus’: Jupiter is only twice as bright as Venus to the observer. The same rule of thumbs can explain the surprising faintness of Mars, Uranus, and Neptune. They show a relatively small reflective area because of their distance to the Sun, which makes them receive little flux; furthermore, Mars has a significantly low albedo. Note that these comparative values are averaged over the entire orbits.
Finally, this same figure highlights once more that focusing the discussion on averaged contrast ratios and separation is futile: depending on the respective configuration of the planets, Earth might easily become the brightest planet in the solar system at some epochs. As an example, let us consider the case when the observer sees the Solar System with an inclination of 45°. Figure 7 shows the contrast ratios of Jupiter, Venus, and Earth as a function of time during a whole Jovian orbit.
It can be calculated that, during the Jovian revolution (11.9 years):

CR(Earth) > CR(Venus) for 2.55 years (21.53% ofT_{Jup}),

CR(Earth) > CR(Jupiter) for 1.04 years (8.82%),

CR(Venus) > CR(Jupiter) for 3.16 years (26.61%),

Venus is the brightest of the 3 planets for 3.11 years (26.16%),

Earth is the brightest of the 3 planets for 4.83 months (3.39%),

Jupiter is the brightest of the 3 planets for 8.36 years (70.44%),
where CR stands for contrast ratio. These values become 28.28%, 18.89%, 32.14%, 30.30%, 4.83%, and 64.86% respectively when the observer sees the solar system at an inclination of 75°. Note also that for an inclination of 45°, Mercury is brighter than Jupiter during 1.75% of T_{Jup}. This duration increases to 11.62% for an inclination of 75°.
Despite the large variability in apparent brightness, it remains true that all solar system planets are extremely challenging: Earth’s relative magnitude to the Sun is about ≈ +23 mag, at best. Gas giant planets at the snow line such as Jupiter do not always offer dramatic gains in observability, as might naively be thought.
3.3. Effect of extensive planetary rings
Rings were not included in this model. To do so would at least double the number of parameters describing a planet (exoplanet obliquity, obliquity at periapsis, ring radii (inner and outer), reflectance parameters, etc.). The ring orientation to the observer is critical for determining reflected light from the hoststar. Beyond the model of our own solar system, there is little observational data to constrain the speculative range of ring properties. However, their presence around an exoplanet would significantly change the observable properties of planets. Rings can either act as reflectors with a potential polarized enhancement or obstruct of the exoplanet illumination or irradiance, depending on their apparent inclination to the star, the observer and their polarizedreflectance characteristics. Both reflection and obstruction effects are coupled for extensive systems of rings that project shadows onto the exoplanet atmosphere. We refer to Dyudina et al. (2005) or Arnold & Schneider (2006) for a more detailed discussion of the impacts of rings on light curves.
3.4. Distinguishing radius and albedo with polarization
Knowing both the contrast ratio (or absolute flux) of the planet and its polarization degree, preferably at several orbital positions, allows the observer to distinguish between several atmospheric models of the atmosphere and finally determine its radius and albedo. For a given exoplanet atmosphere and the same unpolarized flux measured, a planet with large radius and low albedo will exhibit a higher polarization degree than a small exoplanet with high albedo. We refer to the relation between peak polarization and singlescattering albedo in Madhusudhan & Burrows (2012). Such radiusalbedo measurements were carried out in Berdyugina et al. (2008, 2011).
4. Conclusions
4.1. Model and calculation code
We have constructed a model to fully describe the emerging radiation field from an exoplanet with given orbital geometry and reflectance parameters. This model aclculates the three main observables that are relevant for direct imaging of an exoplanet: polarization, contrast ratio, and angular separation as a function of date. A key strength of this model is the relatively restricted number of free input parameters despite the complex processes addressed. Only 3 of these input parameters – eccentricity, inclination, and argument of periapsis – are mandatory for performing a first assessment of the variability over time of the exoplanet observability.
Fig. 7 Evolution of the contrast ratio to the Sun for Venus, Jupiter and the Earth, during an entire Jovian revolution (11.9 years). 
The calculation code developed is available for public use and collaborative improvement on the python package index https://pypi.python.org/pypi under “exospoc”, together with its documentation. This code is implemented on the http://exoplanet.eu/ exoplanet database (see Schneider et al. 2011), where SPOC diagrams are accessible in the exoplanet information sheets at the link “Observability Predictor”.
4.2. SPOC diagrams
We have introduced a novel tool − the SPOC diagram. It highlights the interdependency of polarization, contrast ratio, and angular separation to the host star for an input exoplanet and gives the observer critical information for predicting expected exoplanetary signal from a minimum set of parameters. The complex shape of the SPOC diagram curves highlights the fact that in the general case (inclination higher than ≈10°), the critical observables strongly depend on the geometry of the orbit as it appears to the observer. As a consequence, we stress that restricting consideration to values averaged along the entire orbit for angular separation, contrast ratio, or polarization may be a misleading oversimplification in many applications. An illustration of this is that an external observer would report Earth to be the brightest of the eight planets for a significant fraction of random observations.
The variation of the main observational criteria is critical in timing a direct observation of an exoplanet. This variation mainly relies on the combined effect of the inclination, the eccentricity and the argument of periapsis, hence the shape of the orbit as it appears to the observer. This latter parameter has a surprisingly important role to play in the visibility prediction for a planet in acting to either amplify or cancel the effects of inclination and eccentricity over the contrast ratio and angular separation values over time. The albedo parameters only shifts the lowest and highest values of the contrast ratio; however they do not significantly change the span of its minimum and maximum. The planet radius, semimajor axis and, observerhoststar distance parameters simply scale all observables to higher or lower contrast ratios and angular separations.
The simultaneity of the maxima of three main observational criteria can also be studied with SPOC diagrams. The simultaneity of contrast ratio and polarization (or of contrast ratio and angular separation) maxima is unlikely (or very unlikely) in the t_{CR,max} ± (20% of T) temporal window. Therefore, the bestcase scenario for detection, maxima of contrast, separation, and polarization, is extremely unlikely.
Acknowledgments
We thank Vincent Coudé du Foresto and Sylvestre Lacour for their comments and advice. We acknowledge that this work was partly supported by the French National Agency for Research (ANR13JS05000501). We thank the exoplanet.eu team and especially Marco Mancini, who adapted this work to the online database.
References
 Arnold, L., & Schneider, J. 2006, in Direct Imaging of Exoplanets: Science and Techniques, eds. C. Aime, & F. Vakili, IAU Colloq. 200, 105 [Google Scholar]
 Bailey, J. 2007, Astrobiology, 7, 320 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Berdyugina, S. V., Berdyugin, A. V., Fluri, D. M., & Piirola, V. 2008, ApJ, 673, L83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Berdyugina, S. V., Berdyugin, A. V., Fluri, D. M., & Piirola, V. 2011, ApJ, 728, L6 [NASA ADS] [CrossRef] [Google Scholar]
 Buenzli, E., & Schmid, H. M. 2009, A&A, 504, 259 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cahoy, K. L., Marley, M. S., & Fortney, J. J. 2010, ApJ, 724, 189 [NASA ADS] [CrossRef] [Google Scholar]
 Dumusque, X., Pepe, F., Lovis, C., et al. 2012, Nature, 491, 207 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Dyudina, U. A., Sackett, P. D., Bayliss, D. D. R., et al. 2005, ApJ, 618, 973 [NASA ADS] [CrossRef] [Google Scholar]
 Fluri, D. M., & Berdyugina, S. V. 2010, A&A, 512, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hatzes, A. P. 2013, ApJ, 770, 133 [NASA ADS] [CrossRef] [Google Scholar]
 Hough, J. H., Lucas, P. W., Bailey, J. A., et al. 2006, PASP, 118, 1302 [NASA ADS] [CrossRef] [Google Scholar]
 Karalidi, T., & Stam, D. M. 2012, A&A, 546, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Keller, C. U. 2006, in SPIE Conf. Ser., 6269 [Google Scholar]
 Kemp, J. C., Henson, G. D., Steiner, C. T., & Powell, E. R. 1987, Nature, 326, 270 [NASA ADS] [CrossRef] [Google Scholar]
 Lucas, P. W., Hough, J. H., Bailey, J. A., et al. 2009, MNRAS, 393, 229 [NASA ADS] [CrossRef] [Google Scholar]
 Madhusudhan, N., & Burrows, A. 2012, ApJ, 747, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Norris, B., Schworer, G., Tuthill, P., et al. 2015, MNRAS, 447, 2894 [NASA ADS] [CrossRef] [Google Scholar]
 Roelfsema, R., Bazzon, A., Schmid, H. M., et al. 2014, in SPIE Conf. Ser., 9147, 3 [Google Scholar]
 Schmid, H. M., Beuzit, J.L., Feldt, M., et al. 2006, in Direct Imaging of Exoplanets: Science and Techniques, eds. C. Aime, & F. Vakili, IAU Colloq. 200, 165 [Google Scholar]
 Schneider, J., Dedieu, C., Le Sidaner, P., Savalle, R., & Zolotukhin, I. 2011, A&A, 532, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Thalmann, C., Schmid, H. M., Boccaletti, A., et al. 2008, in SPIE Conf. Ser., 7014 [Google Scholar]
 Wiktorowicz, S. J. 2009, ApJ, 696, 1116 [NASA ADS] [CrossRef] [Google Scholar]
 Wiktorowicz, S. J., MillarBlanchaer, M., Perrin, M. D., et al. 2012, in SPIE Conf. Ser., 8446 [Google Scholar]
 Zugger, M. E., Kasting, J. F., Williams, D. M., Kane, T. J., & Philbrick, C. R. 2010, ApJ, 723, 1168 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Observability predictor examples
Fig. A.1 Contrastseparation diagram for Vband reflectance parameters (no polarization information added for clarity reasons) showing the famous HR 8799 exoplanetary system with its four confirmed candidates. The planet e was assumed to be the same size as the three other planets (1.2R_{Jupiter}). All four planets but planet d follow a “back and forth” locus because of their null eccentricity. Despite their large reflecting area, the contrast ratios are extremely low because of their distance to the host star. 
Fig. A.2 SPOC diagram for HD 80606 b, for which all orbital parameters are known down to a very high accuracy. The very high eccentricity (0.93) of this transiting planet gives a peculiar phase diagram. The higest contrast ratio is not reached at periapsis because of the combination of i, e, and ω_{p} orbital parameters. Photon count is given per hour per VLT collecting area. The reflectance parameters were assumed to be Jupiterlike, in Vband. 
All Tables
Two sets of probable orbital parameters for Alpha Centaury Bb, see Sect. 9 of the supplementary information of Dumusque et al. (2012).
Planetary and orbital parameters for the planets of the solar system, in the Vband.
All Figures
Fig. 1 Geometric albedo A_{g,Rayleigh} arising from singlescattering as a function of the Rayleigh singlescattering albedo w_{R}. 

In the text 
Fig. 2 Phase functions φ as a function of phase angle α for several scattering models. The illuminated fraction is also displayed (thick line). The Rayleighscattering model corresponds to a vectorial Rayleigh phase matrix. 

In the text 
Fig. 3 SPOC diagram for Cen_{α,0.34} with i = 15°, see Table 4, assumed to be Venuslike regarding its reflectance parameters. It represents the polarization values as gray scale, plotted against the angular separation and the contrast ratio axes. The exoplanet’s orbit is therefore displayed in a separationcontrast phase diagram. The wavelength is the Vband. The “x” indicates the position of the exoplanet at an arbitrary observation time. The diamond indicates the position of the planet one sixth of the orbital period T later (giving the forward direction of time). The 10 plus signs indicate the evolution of the planet along the curve; they are linearly spaced in time. The empty and filled triangles are periapsis and apoapsis: r_{min} = 0.026 AU and r_{max} = 0.054 AU. The empty and filled circles are minimum and maximum phase angle: α_{min} = 75° and α_{max} = 105°. The bottom left panel shows the orbit as seen by the observer, where the shaded area corresponds to the portion of the orbit that is behind the plane of the paper. The dashed lined corresponds to the limit of the same instrument sensitivities as in the preceding example. The righthand side yaxis is given in photons per VLT collecting area per hour integrated over Vband. Two successive plus signs correspond to 7h46min. 

In the text 
Fig. 4 Same caption as Fig. 3 except that it represents the apparent orbital motion in mas per hour as gray scale, plotted against the angular separation and the polarized contrast ratio axes. The red dotted line corresponds to the SPHEREZIMPOL performance for 1 h exposure (observability domain being above the line). 

In the text 
Fig. 5 Contrast ratio versus angular separation plots for both orbital parameters of Alpha Centaury Bb as described in Fig. 3 and for inclinations in [15, 35, 60, 85]°. Dotted lines correspond to e = 0 and ω_{p} = 0° and full lines to e = 0.34 and ω_{p} = 246°. A color version is available in the online journal. 

In the text 
Fig. 6 Contrastseparation diagram for the planets of the solar system as seen from 1pc distance in visible light. Small black curves correspond to an inclination of the ecliptic of 0° (observer reference frame); colored curves show the solar system at an inclination of 45°. Note that both axes have logarithmic scales. A color version is available online. 

In the text 
Fig. 7 Evolution of the contrast ratio to the Sun for Venus, Jupiter and the Earth, during an entire Jovian revolution (11.9 years). 

In the text 
Fig. A.1 Contrastseparation diagram for Vband reflectance parameters (no polarization information added for clarity reasons) showing the famous HR 8799 exoplanetary system with its four confirmed candidates. The planet e was assumed to be the same size as the three other planets (1.2R_{Jupiter}). All four planets but planet d follow a “back and forth” locus because of their null eccentricity. Despite their large reflecting area, the contrast ratios are extremely low because of their distance to the host star. 

In the text 
Fig. A.2 SPOC diagram for HD 80606 b, for which all orbital parameters are known down to a very high accuracy. The very high eccentricity (0.93) of this transiting planet gives a peculiar phase diagram. The higest contrast ratio is not reached at periapsis because of the combination of i, e, and ω_{p} orbital parameters. Photon count is given per hour per VLT collecting area. The reflectance parameters were assumed to be Jupiterlike, in Vband. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.