Feasibility of spectro-polarimetric characterization of exoplanetary atmospheres with direct observing instruments
1 Nishi-Harima Astronomical Observatory, Center for Astronomy, University of Hyogo, 407-2, Nishigaichi, Sayo, 679-5313 Hyogo, Japan
2 Department of Earth and Space Science, Graduate School of Science, Osaka University 1–1, Machikaneyama-Cho, Toyonaka, 560-0043 Osaka, Japan
Received: 28 January 2016
Accepted: 9 November 2016
Context. Spectro-polarimetry of reflected light from exoplanets is anticipated to be a powerful method for probing atmospheric composition and structure.
Aims. We aim to evaluate the feasibility of the search for a spectro-polarimetric feature of water vapor using a high-contrast polarimetric instrument on a 30–40 m-class ground-based telescope.
Methods. Three types of errors are considered: (a) errors from the difference between efficiencies for two orthogonally polarized states; (b) errors caused by speckle noises; and (c) errors caused by photon noise from scattered starlight. Using the analytically derived error formulas, we estimate the number of planets for which feasible spectro-polarimetric detection of water vapor is possible.
Results. Our calculations show that effective spectro-polarimetric searches for water vapor are possible for 5 to 14 known planets. Spectro-polarimetric characterization of exoplanetary atmospheres is feasible with an extremely large telescope and a direct observing spectro-polarimeter.
Key words: planetary systems / techniques: polarimetric / instrumentation: polarimeters
© ESO, 2017
Spectro-polarimetric observations of visible or near-infrared reflected light from exoplanets can be a powerful tool for probing their atmospheres. At molecular absorption wavelengths, the polarization degree of planetary reflected light increases owing to the decrease in intensity of the multiply scattered component as compared with that of the continuum wavelengths. Stam et al. (2004) and Stam (2008) calculated polarization degree spectra for Jupiter-like and Earth-like planets, respectively, and derived enhanced features due to absorption by atmospheric molecules such as CH4, H2O, and O2. In fact, Sterzik et al. (2012) and Miles-Páez et al. (2014) clearly detected H2O (at wavelengths 0.72, 0.82, 0.93, and 1.12 μm) and O2 (0.69, 0.76, and 1.25 μm) features on the polarization spectra of lunar Earthshine. Hence, spectro-polarimetry for exoplanets can be used for the detection of these atmospheric molecules.
Spectro-polarimetry facilitates not only the detection of atmospheric molecules, but also further characterization of the atmosphere. Stam et al. (2004) and Stam (2008) showed that the polarization of planets is sensitive to the structure of the atmosphere. Stam et al. (2004) simulated three different model atmospheres of Jupiter-like planets: (1) an atmosphere that contains only gaseous molecules (“clear atmosphere”); (2) an atmosphere with a tropospheric cloud layer (“low cloud atmosphere”); and (3) an atmosphere with a tropospheric cloud layer and a stratospheric haze layer (“high haze atmosphere”). The resulting differences between the polarization spectra from the models were significant in overall wavelength-dependence and strength of the molecular features, whereas the differences in the intensity spectra were less clear. Similarly, Stam (2008) pointed out that the degree of polarization at the continuum wavelengths near a molecular feature is more sensitive to the altitude of the cloud’s top compared with intensity. Therefore, spectro-polarimetry provides information that cannot be obtained solely by intensity spectroscopy.
In addition to these scientific benefits, polarimetry has some technical advantages. Because polarimetry is a relative measurement, the degree (and position angle) of polarization is not affected by isotropic transmittance between celestial objects and observers. This advantage will be especially effective for ground-based detection of H2O, O2, and other components common with the Earth’s atmosphere, if the telluric transmittance is isotropic.
The degree of polarization P of a planet is derived by the equation: (1)where and represent planetary intensities of light oscillating perpendicular and parallel to the scattering plane (the plane which includes the central star, the planet, and the observer), respectively, and λ is the wavelength. If the instrumental reference plane is aligned to the scattering plane (using a half-wave plate), we have Stokes parameters U = 0 and Q = P, because scattering or reflected light is polarized perpendicular or parallel to the scattering plane. We take and as the values that are not altered by any absorption. When observations are conducted at the ground, intensities are affected by telluric transmittance.
Schematic illustration of the hypothetical observational system considered in this paper. It consists of a) a telescope; b) a half-wave plate; c) extreme adaptive optics; d) a coronagraph or a nulling interferometer; e) an integral field spectrograph; f) a polarizing beam splitter; and g) a detector.
|Open with DEXTER|
Previously, it was believed that telluric transmittance was isotropic, and thus did not affect the polarization state of light from an astronomical object. However, Bailey et al. (2008) found that dust in the telluric atmosphere can produce polarization up to ~ 5 × 10-5. In addition, the polarization state is affected by interstellar transmittance, instrumental transmittance, and any offset error in the intensities. Therefore, we should express measured intensity values (I′p) as where rχ is the total efficiency from the object to the detector for ordinary (χ = 0°) or extraordinary (χ = 90°) components. It is the product of interstellar transmittance (), telluric transmittance (), and instrumental throughput including detector efficiency (kχ). ϵχ is any offset error in the measurement of intensity values. In this paper, we examine how rχ and ϵχ affect the measured degree of polarization.
In the 2020s, 30–40 m-class ground-based telescopes such as the Thirty-Meter Telescope (TMT), the European Extremely Large Telescope (E-ELT), and the Giant Magellan Telescope (GMT) will begin operation. Direct spectro-polarimetry or narrow-band imaging polarimetry for exoplanets is anticipated to be realized by a combination of high-contrast and polarimetric instruments on these extremely large telescopes (ELTs). Polarimetric differential imaging (PDI) is one of the highly effective techniques used to enhance the contrast range of an instrument and facilitate the detection of reflected light from planets around stars. PDI allows not only the detection of planets but also the measurement of polarization. PDI is widely utilized for currently operating instruments, such as HiCIAO (e.g., Tamura et al. 2006; Hashimoto et al. 2011; Tanii et al. 2012) on the Subaru telescope, and the SPHERE-ZIMPOL (e.g., Thalmann et al. 2008, 2015) and NACO (e.g., Lenzen et al. 2003; Rousset et al. 2003; Garufi et al. 2013) on the Very Large Telescopes (VLTs). A high-contrast direct-imaging polarimeter for exoplanets, EPICS-EPOL, is designed based on the ZIMPOL and is proposed for the E-ELT (Kasper et al. 2010; Keller et al. 2010). Furthermore, a concept of the spectro-polarimetric option for EPICS-EPOL has also been studied (Rodenhuis et al. 2012). A spectro-polarimeter based on their concept will enable an effective suppression of the residual speckles by combining PDI with simultaneous spectral differential imaging (SSDI) or chromatic speckle suppression (Marois et al. 2000; Sparks & Ford 2002). Therefore, at this point, it is important to assess the feasibility of polarimetry using a high-contrast spectro-polarimeter with SSDI and PDI techniques.
In this paper, we evaluate the feasibility of spectro-polarimetric characterization of exoplanetary atmospheres by taking into account the errors expected in high-contrast instruments. As shown Fig. 1, we consider a hypothetical spectro-polarimeter based on the concept proposed by Rodenhuis et al. (2012). The instrument shall be equipped with an integral field spectrograph (IFS) and a polarizing beam splitter (PBS) to obtain a polarization degree spectrum of a planet. Although the EPIS-EPOL is a visible light instrument, we assume that our instrument works in the near-infrared because we focus on the detection of the strongest spectro-polarimetric feature of H2O at 1.12 μm on a visible-to-near-infrared spectrum (Miles-Páez et al. 2014). We assume that our instrument has a medium wavelength resolution R of a few thousand. Although our target spectral feature is relatively broad (Δλ ≅ 0.05 μm, Miles-Páez et al. 2014), we follow the setting adopted by Barman et al. (2011) and Konopacky et al. (2013), who detected CH4 and H2O in exoplanet atmospheres by using the SSDI method for R ~ 4000 data.
We focus on H2O vapor, which is a fundamental component in planetary atmospheres. H2O plays a major role in meteorology and habitability of Earth and other Solar System planets. Detection of H2O vapor has already been reported for several hot-Jupiters based on transit spectroscopy (e.g., Barman 2007; Tinetti et al. 2007; McCullough et al. 2014). However, transit spectroscopy is only possible for close-in planets. Direct observations can investigate planets in different orbital areas such as those near the habitable zone and the snow line. In addition, successful spectro-polarimetric searches for a H2O feature will provide information on the structure of planetary atmospheres, such as cloud/haze height. Note that the framework we construct for H2O detectability evaluation can easily be applied to other atmospheric species.
We discuss the detectability of the H2O spectro-polarimetric feature with ΔP ≡ P(λabs)−P(λcnt) = 10% and P(λcnt) = 10%, where λabs and λcnt represent absorption and its nearby continuum wavelengths, respectively (λabs = 1.12μm). These values are based on lunar Earthshine observations by Miles-Páez et al. (2014) with an approximate correction of depolarization at the lunar surface (Earth pol. ≅ Earthshine pol. × 3, Dollfus 1957). Although our target exoplanets include gas-giants, the setting above is justified for the following reasons; (1) H2O vapor can be abundant if a planet is located within the snow line; and (2) as long as scattering of the atmosphere has a significant contribution to the total planetary intensity (i.e., clouds or hazes are not too high in altitude), a polarization degree spectrum of the planet is expected to have enhancement features at molecular absorption wavelengths (Stam et al. 2004; Stam 2008).
This paper is organized as follows. In Sect. 2, we formulate different types of errors in polarization degrees. Based on the derived equations, requirements for detection of a spectro-polarimetric feature are formulated and evaluated in Sect. 3. Then, we examine the feasibility of spectro-polarimetric searches for H2O vapor in Sect. 4. After discussions on unconsidered error factors and remaining future works in Sect. 5, we summarize and conclude this paper in Sect. 6.
We consider three types of errors: (a) errors arising from the difference between the efficiencies (throughputs) of two orthogonally polarized states, ; (b) errors caused by speckle noises, ; and (c) errors arising from photon noise, . (a) and (b) are systematic biases, whereas (c) is a random error. Each type of error is formulated below by expanding Eq. (1). Note that, in this section, we do not write out wavelength dependence in functions so as to avoid complicating the equations.
We analyze the effect of the difference between r0 and r90. By substituting I′p in Eq. (2) into Ip in Eq. (1), the total error in the polarization degree is written as (4)where (5)The derivation of the equations is given in Appendix A, where and | r90/r0−1 | ≪ 1 are assumed in the approximations. is the bias caused by different efficiencies as it is caused by r90/r0. The barred values (e.g., ) satisfy and δX (e.g., , ) is defined by δX = (X0−X90)/2. The terms , , and are the polarization degrees through interstellar, telluric, and instrumental (including telescope) transmission, respectively. The difference of total efficiencies, r90/r0−1, can be written as the sum of the three individual polarization terms.
Keller (1996) analyzed the effect of nonlinearity in detector responses. Following his formulation, this effect is implemented into as (6)where is the nonlinearity term and has a typical value of ~ 0.01 (Keller 1996).
The second term in Eq. (4) is an error in P caused by offset errors in intensity measurement. It can be either systematic or random: it is systematic if ϵχ is systematic (e.g., due to speckle noises) and random if ϵχ is random (e.g., due to photon noise). Each type of error is analyzed below.
In high-contrast instruments, stellar intensity is reduced with a coronagraph. Even with an extreme AO system, a certain roughness of wavefront cannot be corrected, and hence, the speckle noises are generated on an image plane. The speckle noises disturb the detection and measurement of planetary photons. In this section, we examine the systematic bias in the planetary polarization degree due to the speckle noises.
We examine the second term of Eq. (4). In this paper, we consider a case where the PDI technique is applied after the SSDI procedure. The typical strength of ϵχ (speckle noises) after the SSDI but before PDI is , where Cdet is the detection contrast of the instrument achieved with the SSDI but without the PDI. It describes the ratio of the strength of speckles to the non-coronagraphic stellar intensity, (see Fig. 2 for a schematic explanation and Appendix C for a more rigid description). Thus, we express where σϵ0 + ϵ90 and σϵ0−ϵ90 represent the standard deviations of (ϵ0 + ϵ90) and (ϵ0−ϵ90), respectively. Hence, (ϵ0 + ϵ90) and (ϵ0−ϵ90) are typically within ± σϵ0 + ϵ90 and ± σϵ0−ϵ90, respectively.
Schematic illustration of an image profile before (upper line) and after (lower line) the speckle subtraction procedures. In the situation assumed for this paper, the speckle subtraction procedures include the SSDI but not the PDI at this stage. The PDI is conducted further on.
|Open with DEXTER|
If both of the ordinary and extraordinary components of speckle noises have the same pattern and intensity, then we have ϵ0−ϵ90 = 0. However, it is not true due to, for example, different throughputs and wavefront errors for orthogonal polarized states; and polarization aberrations (deviation of the point-spread functions; Breckinridge & Oppenheimer 2004; Breckinridge et al. 2015). Besides, it is not known yet how well the PDI after SSDI works because a spectro-polarimeter based on the combination of an IFS and a PBS has not been realized. Thus, we introduce a parameter, g, which expresses how well the PDI suppresses the speckle noise. The detection contrast improves from Cdet to gCdet by applying the PDI. We consider the best and worst cases with regard to g. In the best case, g is determined only by the difference in the throughputs, that is, . In the worst case, PDI does not work at all, that is, g = 1.0.
A typical value of the error caused by the speckle noise () is derived from the second term of Eq. (4) using Eqs. (7)–(8). If the spatial patterns of (ϵ0−ϵ90) and (ϵ0 + ϵ90) are similar, we have (9)where is the intensity contrast of a planet to its host star. Contrarily, if (ϵ0−ϵ90) and (ϵ0 + ϵ90) are not correlated, we should express (10)Because the highly correlated components between ϵ0 and ϵ90 will largely disappear in (ϵ0−ϵ90), we expect that Eq. (10) is the case.
is a function of Cdet/Cp/s. It is understandable that is smaller when Cdet is better (smaller) and Cp/s (the planetary brightness relative to the star) is greater. When , we have . This means that we have a smaller if we have a better (smaller) g, as expected. On the other hand, when , we have . This indicates that, even if the PDI were perfect (g = 0), P would suffer from an error caused by the speckles, because we need to use not only a differential (I′0−I′90) image but also a sum (I′0 + I′90) image to derive P.
It may be puzzling that we have a larger error when the planet is more polarized. This is a result of the fact that the change of a P value due to virtually unpolarized speckle noises is greater when the planet is more polarized. Hence, most high-contrast polarimetry data is expressed in terms of polarized intensity, not degree of polarization. However, for ground-based spectro-polarimetric detection of H2O (and other components common to Earth’s atmosphere, e.g., O2), avoidance of deriving the degree of polarization (i.e., dividing by the sum data) may pose a problem because one must cancel out telluric features using a standard star spectrum, which will be observed at a different time and sky position from the target. Therefore, we continue to study the errors in degree of polarization, rather than polarized intensity.
As we assume the scaling factor at the head of Eq. (10) ranges from the order of 10-1 to unity, Cdet/Cp/s should be less than the order of 10-1 to 10-2 so as to achieve a of a few percent.
If ϵ0 and ϵ90 in Eq. (2) are caused by photon noise and are thus random, the second term of Eq. (4) represents a random error () in P. We consider scattered starlight remaining after a coronagraph as the dominating source of the photon noise. As the intensity of scattered starlight is , we have (11)Hence, (12)where is expressed in photon numbers. Naturally, the error from photon noise is smaller when the host star is brighter ( is greater), the planet is brighter (Cp/s is greater), and the instrument raw contrast is better (Craw is less).
In summary, the following types of errors in P are formulated:
In this section, requirements for detecting a spectro-polarimetric feature are formulated based on the derived equations for , , and (Eqs. (13)–(15)). Then, we calculate these errors to obtain the requirements of the instruments.
A spectrum has a spurious feature at the absorption wavelength λabs, as described later. Therefore, the absolute strength of the spurious feature must be significantly smaller than that of the real feature on a P spectrum for its detection. In other expressions, the following requirement should be satisfied: (16)where ΔP( ≡ P(λabs)−P(λcnt)) is the polarization degree of the feature and w is a coefficient determined by the detection definition. w is set to 1/3 throughout the paper.
Input values to the calculations.
We verify that the spectrum has a feature at λabs by examining the three factors in Eq. (13). The following properties are found:
The first factor, , is presumably featureless. The reasons are givenin the following paragraph.
The second factor, (1−P2), has a feature if the P spectrum has a real feature. If P does not have a feature, neither does (1−P2), which is quite evident.
The third factor, , has a feature, even if Ip (and thus P) does not have a feature. This is because the detector nonlinearity term, is proportional to and telluric transmittance, , always has an H2O absorption feature at λabs.
At the end, the spectrum of will have a spurious feature at λabs. Therefore, the requirement Eq. (16) must be satisfied.
We describe the reasons why is presumed to be featureless by looking at each term separately. Although may have a weak feature at λabs, interstellar polarization does not have a feature at λabs because the polarization is caused by aligned dust particles and is not due to H2O molecules (see P spectra of polarized standard stars, which are polarized through interstellar transmission, e.g., Fossati et al. 2007). has a strong feature at λabs. However, similar to the case of , telluric polarization should be featureless at λabs because the polarization is produced by dust particles (Bailey et al. 2008). Since k is determined by optical elements, k and do not have a feature at λabs. Thus, we regard the spectrum of as featureless.
Here we address the , , and values, which are used for the calculation. All of them are wavelength-dependent. However, in order to simplify the calculation, we treat , , and as constants within a minimum wavelength range necessary to detect a feature, that is, the range between λabs ± Δλ (Δλ is the FWHM of the feature, ~0.05 μm). This assumption is allowed because they are expected to be featureless, as described above.
The interstellar polarization () is set to be 1% as a safe assumption. Note that most strongly polarized standard stars, which are polarized through a relatively long (a few hundred parsec) interstellar path, have a polarization degree of one to a few percent (Whittet et al. 1992). As direct observation of exoplanets is limited to relatively nearby systems (within ~30 pc), we expect that the interstellar polarization for many of the targets is likely to be lower than the assumed value. As for polarization through telluric transmission (), we take a value 5.0 × 10-5, measured by Bailey et al. (2008). We assume instrumental (including telescope) polarization () to be 5%, based on modeling for the E-ELT by de Juan Ovelar et al. (2012). We do not know the signatures (position angles of polarization) of , , and , but assume all of them to be positive so as to maximize the value of their sum.
For the purpose of understanding how is determined, it is written as (17)where Δ [X] represents X(λabs)−X(λcnt). It is seen that the contribution from is negligible. We also find that the ratio of contributions from the and P is approximately 1:5.
Bias in P with respect to Cdet/Cp/s(λabs). The speckle error for g = 0.05 (solid line) and g = 1.0 (dashed line) are plotted. The efficiency error (dotted line) is also shown. The shaded areas show where detection is difficult because the bias is larger than wΔP. The arrows point to the maximum limit of Cdet/Cp/s(λabs) for the feature detection. w is set to one third.
|Open with DEXTER|
with respect to for different stellar magnitudes: m∗ = 0.0 (dotted line), 5.0 (solid line), 10.0 (dashed line). Magnitude m∗ determines in Eq. (15). The shaded area shows where detection is difficult because is larger than wΔP. The arrows point to the maximum limit of for detection. w is set to one third.
|Open with DEXTER|
Speckles move outward from the host star on an image plane as the wavelength increases (see Fig. 23 in Sparks & Ford 2002), whereas the position of a planet image is independent of the wavelength. The SSDI technique removes the speckles using this property; however, some residuals will remain. Because the speckles “cross” a planet on a space of the separation (θ) and the wavelength (λ; see Fig. 2 in Rodenhuis et al. 2012), they will appear as spurious features on the I′p and P spectra. It is difficult to distinguish the spurious speckle features from the real molecular feature based on their spectral widths (see Appendix B). Therefore, it must be guaranteed that the strength of the speckle features on the P spectrum is sufficiently small compared to that of the real feature, ΔP.
Unlike the efficiency error, , we only know the standard deviation of the speckle error, . Thus, we need to define the detection requirement differently: both and must be significantly smaller than ΔP. Owing to a greater P and a smaller Cp/s at λabs than those at λcnt, it is expected that , assuming Eq. (16) is satisfied (see Eq. (14)). Hence, the requirement regarding is expressed as (18)Using Eq. (14), the values of are calculated as a function of Cdet/Cp/s(λabs). As described, we consider two cases with regard to g: g = 0.05 as the best case and g = 1.0 as the worst. All the other input values are given in Table 1.
Figure 3 displays together with . As shown in Eq. (14), is linearly correlated to Cdet/Cp/s(λabs). For the best case (g = 0.05), Cdet should be better than ~ 1.3 × 10-1Cp/s(λabs), which means that a contrast performance (before PDI) to achieve ~8-sigma detection of the planet signal is required. For the worst case (g = 1.0), Cdet should be better than ~ 3.2 × 10-2Cp/s(λabs), meaning a requirement of ~30-sigma detection performance. Naturally, the requirement for Cdet is less strict for the case of a better g.
The detectability is limited by the speckle error rather than the efficiency error because is small enough compared to ΔP. Thus, hereafter, the requirement pertaining to the efficiency error (Eq. (16)) is neglected in the evaluation of the detectability, though in Eq. (14) is still taken into account.
The error due to photon noise should be smaller than ΔP both at λcnt and λabs. Because of the smaller and Cp/s at λabs than at λcnt, we have . Thus, the requirement regarding the photon noise is written as (19)Using Eq. (15), the values of are calculated as a function of for stellar magnitudes m∗ = 0.0, 5.0, and 10.0 (Fig. 4). Telescope diameter D = 39 m (E-ELT) and total integration time Δt = 15 h (5 h × 3 nights) are taken to derive the stellar photon counts . The other settings are summarized in Table 1.
As seen in Fig. 4, for a host star of magnitude m∗ = 5.0, should be smaller than ~ 1 × 105. In other words, Craw should be smaller than ~ 1 × 1010(Cp/s(λabs))2. As might be expected, the requirement for Craw is less strict if the host star is brighter.
We discuss the feasibility of spectro-polarimetry through estimating the number of planets for which a H2O feature with P(λcnt) = 10% and ΔP = 10% is detectable. The detectability is estimated according to Eqs. (18) and (19) by calculating and for the 2041 discovered planets listed in the Exoplanet.eu (the Extrasolar Planets Encyclopaedia, Schneider et al. 2011) database as of January 7, 2016. Assuming we use the E-ELT at Cerro Armazones, Chile, we exclude objects with declination >+ 35°.
In Eq. (15) Craw is taken from the values delivered by the AO system (called XAO) for EPICS, a planned direct observing instrument for exoplanets (Fig. 3 in Kasper et al. 2010). We choose EPICS-IFS as the reference of Cdet in Eq. (14) (Fig. 4 in Kasper et al. 2010), because: (1) it is actually being studied for the E-ELT; (2) it covers the λabs (=1.12 μm); and (3) its Cdet can be applied directly to our formulae because the values without the PDI procedure are provided in Kasper et al. (2010). Note that EPICS-IFS itself will not have a spectro-polarimetry capability. We consider a hypothetical spectro-polarimeter with a contrast performance (before PDI) comparable to that of EPICS-IFS. The Cdet values of EPICS-IFS are converted from 5-sigma to 1-sigma values to fit our formulation because we use an external parameter w to define the criteria of the detection.
In Eqs. (14) and (15), Cp/s for each planet is calculated by where p(λ) is the geometric albedo of a planet at wavelength λ, Φα is the phase law at phase angle α for a Lambert sphere, Rp is the planet radius, ap is the distance from the host star to the planet, and d is the distance to the star from an observer (Traub & Oppenheimer 2010). p(λ) and α are fixed for all the planets (Table 1). Rp, ap (assumed to be close to the semi-major axis), and d is obtained primarily from the Exoplanet.eu database. If d is not available, we transform it from the parallax data in the SIMBAD database (Wenger et al. 2000). If Rp is not found in the Exoplanet.eu database, we predict it based on the empirical mass-radius relations studied by Weiss et al. (2013) and Weiss & Marcy (2014), where the planet mass (Mp) is acquired from the Exoplanet.eu. Different relations are applicable, depending on the different ranges of the planet mass (Mp): (1) Mp> 150 MEarth; (2) 10 MEarth<Mp ≤ 150 MEarth (Weiss et al. 2013); (3) 4 MEarth<Mp ≤ 10 MEarth, and (4) Mp< 4 MEarth (Weiss & Marcy 2014). Weiss et al. (2013) show that the radius is weakly (exponents −0.03 or 0.094) correlated with the incident stellar flux on the planet for the Mp ranges of (1) and (2). The time-averaged incident flux is calculated using the stellar effective temperature, the stellar radius, the semi-major axis, and the orbital eccentricity, which are retrieved primarily from the Exoplanet.eu database. If the eccentricity is not provided in the database, it is assumed to be zero. If either the stellar temperature or radius is not in the database, we use the mean flux in Weiss et al. (2013; 8.6 × 108 erg s-1 cm-2).
Mass and orbital distances of the observable planets. The open circles are planets observable with a setting of g = 1.0 and Δt = 15 h. The filled larger circles are those with a setting of g = 0.05 and Δt = 15 h. The filled smaller circles are those with a setting of g = 0.05 and Δt = 50 h.
|Open with DEXTER|
List of planets for which spectro-polarimetric detection of H2O vapor is possible with g = 0.05.
The SIMBAD database (Wenger et al. 2000) is used for the reference of the J-band magnitudes (m ∗) of the host stars, which determine in Eq. (15). All the other fixed input values to the calculations are shown in Table 1.
Table 2 lists the planets for which we predict the spectro-polarimetric detection of the H2O feature is possible with the setting g = 0.05. We have extracted 14 planets, of which 11 are estimated to be located within the snow line. The snow line distance is calculated according to Eq. (17) in Ida & Lin (2005). Among the 14 observable planets with g = 0.05, five planets are observable even with g = 1.0. Figure 5 shows the mass and orbital distances of the observable planets. The mass of the planets range from 7.5 MJup (HD 60532 c) to 0.027 MJup ≅ 0.5 MNep (GJ 682 c).
In order to illustrate our results graphically, we transform the expressions for the detection requirements (Eqs. (18) and (19)) into suitable forms: The right-hand side of the inequations above is interpreted as the contrast of the spectro-polarimetric feature signals to the host star signals. The left hand side of Eqs. (22) and (23) are contrast performances for the feature detection, constrained by the speckle noises and the photon noise, respectively.
In Fig. 6, we plot the (i) speckle-limited contrast (the left-hand side of Eq. (22)); (ii) photon-limited contrast (the left-hand side of Eq. (23)); and (iii) spectro-polarimetric feature contrast (ΔPCp/s) against the separation from the star. When both lines for (i) and (ii) are below (iii) in the figure, the feature is detectable. Four different cases are shown: (a) a “typical” planet with averaged properties of the observable planets with g = 1.0; (b) same as (a) except for g = 0.05. (c) Gliese 876 b with g = 1.0; and (d) Gliese 876 b with g = 0.05. Except for (a), detectability is determined by the photon-limited contrast for a certain θ range. Thus, the detectability can be improved by obtaining more photons (e.g., with a longer Δt or a better ). In fact, if we extend Δt to 50 h (with g = 0.05), the number of observable planets increases to 20 (Fig. 5).
Note that we have extracted observable planets from the already known planets. Further, observable planets are not necessarily limited to those in Table 2, given that transit and radial velocity detection of exoplanets is biased toward edge-on systems, whereas direct observations can also find those in face-on systems.
Among the listed observable planets, Aldebaran b orbits the brightest star (m∗(J) = −2.1 mag). The existence of a planetary companion with a mass of 6.47 ± 0.53MJup is considered to be the most plausible explanation for the observed radial velocity (RV) variations with a 629-day period (Hatzes et al. 2015). Owing to the extreme brightness of the host star, for the planet will be suppressed to a very low level (~0.2%). Aldebaran b, if it truly exists, may be one of the first targets for the spectro-polarimetry.
Gliese 876 (=GJ 876) is currently known to host four planets. Discoveries of planets “b” and “e” were made based on RV observations by Delfosse et al. (1998) and Marcy et al. (1998); and Rivera et al. (2010); respectively. It is expected that both of the planets can be observable even with g = 1.0. The contrast plot for Gliese 876 b is shown in Figs. 6c and d for g = 1.0 and 0.05, respectively. The orbital distance of the snow line in the system is evaluated to be 0.30 AU based on the mass of the host star (0.334 M⊙). Fortunately, the orbital distance of planet “b” (ap = 0.21 AU) is estimated to be within the snow line, whereas that of planet “e” (0.33 AU) is predicted to be beyond the line. Therefore, this system may be an optimal target to investigate the existence of H2O vapor on planets inside and outside the snow line.
Comparison of contrasts. Contrast performance limited by the speckle noises (bold solid lines) and by the photon noise (dashed lines) is compared with ΔPCp/s (filled circles), that is, contrast of the polarized feature signal. a) a “typical” planet with averaged properties of the observable planets with g = 1.0. b) same as a) except for g = 0.05. c) Gliese 876 b with g = 1.0. d) Gliese 876 b with g = 0.05. The thin solid line in each panel is ΔPCp/s(θ) for a planet with the same radius. The polarization feature is detectable in the ranges of θ designated by the arrows. 3-sigma Cdet (dashed-dotted lines) and Craw (dotted lines) are also plotted for the reference.
|Open with DEXTER|
In this paper, we explore the feasibility of spectro-polarimetric detection of H2O vapor on exoplanets by evaluating the P errors, which arise mainly from the efficiency differences between polarization states, speckle noises, and photon noise. In this section, some unconsidered factors are discussed.
Influences from the circumstellar disks were not discussed. Min et al. (2012) pointed out that the central resolution element on a disk image can be polarized because of scattering from the innermost regions of the disk. This makes the central star appear to have an intrinsic polarization. However, the polarization degree was calculated to be less than 0.1% for a spatial element obtained through 30–40 m telescopes (Min et al. 2012). As we have assumed much greater interstellar and instrumental polarization (several percent, Table 1), consideration of this effect will not change the result of our feasibility estimates.
As another influence from the disks, we consider exo-zodical light around a planet image. For low-spatial-resolution images, exo-zodical light may obscure a planet image. Simulations by Cash et al. (2008) showed that intensity of an Earth twin in a solar an clases/aa/full4>Fiws where deExoplanet.eu database. If the eccentricity is not provided in the database, it is assumed to be zero and a s>1), conrIo3a>), conrI4-16. Fouef="/arti206-16.ha 4.in en, s notihus, the det al. 201>)). Ht8206- nois1)ent obtaineand <2017/0InR7">s,sub>Jup makessharMin et al. 201,ase. If the ecand ouo3a>), conrI4-16pe_html/2017/0ost regions of the table>
Inf(a), dec bq33.png"> Currefu usage tichics show cumulated icountor oAre sup Views (full-uexteere sup views includicieHTML views, PDFcenceePub downloads, according dond (2avaia6ble data) cenceAhae nots Views ckiVipape4P28200aint lin. Data asrr spondeuleusage oand (2aint lineafcma m15. The2currefu usage tichics is avaia6ble 48-96nhourseafcma onr te publPICaeck and is updbted daiay2ck week days.
100tella/aaafull_ 03/aa28ro and a stricity istml/2017/ws less sng from luence from the disfull_somef="/artiing the
) left hanmbtmany notbe impr. abs). For the best case ( ss="simple-math">~ 3.2 × 10-2C ticlBecausca, ifn7/0 claf trast performance (bs.org/apple/articlrI4-16lutioe mainly f,sub>O vapor oEqs. (P errors, which aThs b for ats insidedine. > to defrti206. Htyppan clpeckleh"span>).spaue requiremd through 30–40 mith
Inf(a), dec bq33.png">for the planet will be suppressed to a very low level (16pecklehl/201dottep/s
Currefu usage tichics show cumulated icountor oAre sup Views (full-uexteere sup views includicieHTML views, PDFcenceePub downloads, according dond (2avaia6ble data) cenceAhae nots Views ckiVipape4P28200aint lin."
Data asrr spondeuleusage oand (2aint lineafcma m15. The2currefu usage tichics is avaia6ble 48-96nhourseafcma onr te publPICaeck and is updbted daiay2ck week days."