EDP Sciences
Free Access
Issue
A&A
Volume 599, March 2017
Article Number A56
Number of page(s) 12
Section Astronomical instrumentation
DOI https://doi.org/10.1051/0004-6361/201628206
Published online 28 February 2017

© ESO, 2017

1. Introduction

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.

thumbnail Fig. 1

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 (Ip) 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 ΔPP(λ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.

2. Formulation of errors

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.

2.1. Different efficiencies

We analyze the effect of the difference between r0 and r90. By substituting Ip 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 = (X0X90)/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.

2.2. Speckle noises

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.

thumbnail Fig. 2

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 (I0I90) image but also a sum (I0 + I90) 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.

2.3. Photon noise

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

2.4. Summary of error formulations

In summary, the following types of errors in P are formulated:

  • Error from different efficiencies(13)

  • Error from speckle noises (14)

  • Error from photon noise (15)

3. Requirements for detecting a spectro-polarimetric feature

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.

3.1. Different efficiencies

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.

Table 1

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.

The value of is calculated using Eq. (13) and the values listed in Table 1. As a result, we obtain , which satisfies the requirement of Eq. (16).

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.

thumbnail Fig. 3

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

thumbnail Fig. 4

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

3.2. Speckle noises

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 Ip 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.

3.3. photon noise

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.

4. Evaluation of feasibility

4.1. Expected number of observable planets

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

thumbnail Fig. 5

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

Table 2

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.

4.2. Notes on some of the observable planets

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.

thumbnail Fig. 6

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

5. Discussions

5.1. Unconsidered factors

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> 100tella/aaafull_ 03/aa28ro and a stricity istml/2017/ws less sng from luence from the disfull_somef="/artiing the ) against s), that is, contrast of the polarized featurarticles/aa/fnet image.5). Owing to the extreme brightness of the host star, 1C. -hould be smallh">~MEarth -hould be smallclass="me="I/aa28206-16/ale> <)). As we have assumed much greater interstellar and instrumental psfulKe and 8206articles/06-16-eub>seale> <) < tomeande">14) (sired. Fohe lew g = ref="/artimath">λcnt and P&mates. <_html/2017/differences between polarization states, h">~ <)). As we have assumed much greater interstellar and instrumental psfulKe and 8206arti;">Cash et al.; Oppenheimer 2010). p(λ)<; OBagnuloan>-band ma9 or a beath"sn claf R15">sub>ticlare nR15">sub>a href">~cnt 6.html#FD13)06-16/aahalf-4;Jute/aa282fpeckelng >) lse id rimpbandath"> = 1.0r influesorknet image.Obv Appb>pticlhusskyg star = 50 h.

) left hanmbtmany notbe impr. ) /span>017. t = 50 h.

) imclanhresprequiremtiing the veraged p206-16.ha anmbt6-16ain. imclPn. Because tanet image.Despite considmor a befeinst sorkcle-math"tudy) a. MO vapor on exoplanets by evaluating the ) cencesrulentedic polaa28206-16/a16">18
) a/i> = 0.05abs). For the best case (ss="simple-math">~ 3.2 × 10-2C~ 3.2 × 1gon states, h"16-eq182.png"> for the 2041 discovered planets listed in the Exoplanet.eu (the Extrasolto 50 h (with 4.1. Expecte (2008)6. S>1ticlvern Eqrbitidmorkcrea2mage.Motived-dotted lin wittr 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 for the planet will be suppressed to a very low level (16pecklehlutio the obserthis sect sub>ticlare nR15">sub>rI4-16p> for the planet will be suppressed to a very low level (16pecklehl/201dotted.

Inf(a), dec bq33.png"> for the planet will be suppressed to a very low level (16pecklehl/201dottep/s)l/2017cla to the th">18) as.orlaa28206-16/ahe photon_html/ten>O feaanhrespe observabP errors, wle-math">λcnt, we have . Tspan> can be apple">ΔP = 10% is detectable. The det%14) imclhat is, contrast of the polarized featurarticles/aa/fneles/aa/mo; planet3/aa28206-16/aa28206-16.html# is predicte16) is les/aa/O vapor oEqs.ple-matle-matFS. The to definef,sub> =/aa28206-16.ht, sarti left4.1. Expecack(2008)Acknowl ) rea2mage.W"me="Ideep4.1. Expec45e(2008)Rlet.php? rea2ma="S153/aaid="he pn>)"a name="S15d se="S12"plet.php?s"e"> [NASA ADS]ppen"> 86/605654">[CrossRle]ppen"> [NASA ADS]ppen"> /j.1365-2966. ma8.13088.x">[CrossRle]ppen"> [NASA ADS]ppen"> 86/518736">[CrossRle]ppen">Weioth aC Chile, ApJ, 733, 65"> [NASA ADS]ppen"> 88

7X/733/1/65">[CrossRle]ppen"> class="sim, B. Rnd ma4, ApJ, 600, 1091"> [NASA ADS]ppen"> 86/379872">[CrossRle]ppen">Chipnes, RndAimg sr, PASP, 127, 445"> [NASA ADS]ppen"> 86/681280">[CrossRle]ppen"> 7/12.7897e (2[CrossRle]ppen">Wenagost star iAa28onomy V,nProc. nPIE, 8449, 844912">Weyo1,5M. 8208, A#R45"A, 33e, 581"> [NASA ADS]ppen"> [NASA ADS]ppen">Lenci Degl’Inno, cois E.d ma7, ven (2leinst of P/i> [NASA ADS]ppen"> [CrossRle]ppen"> [EDP S with D]ppen"> [NASA ADS]ppen"> 88/ m41-8205/729/2/Le (2[CrossRle]ppen"> [NASA ADS]ppen"> [CrossRle]ppen"> [EDP S with D]ppen"> [NASA ADS]ppen"> 86/429953(2[CrossRle]ppen"> [NASA ADS]ppen"> 07/BF0014663 (2[CrossRle]ppen">Weioth aC Chil3, S with , 339, 1398"> [NASA ADS]ppen"> 126/n with .123 m03(2[CrossRle]ppen">),-/sub>A. F. M. MoorwoodclSPIEnoisf CSer., 4841, 944"> tuassuuficietrici’spatmo clasaluceas th pape) f near-"16-efar-iafrircl cles/peck,tTech. rep. (NASA)"> [NASA ADS]ppen"> 86/311623(2[CrossRle]ppen"> [NASA ADS]ppen"> 86/316492">[CrossRle]ppen">Wedhusudhes, N Chil4, ApJ, 79e, 55"> [NASA ADS]ppen"> 88

7X/791/1/55">[CrossRle]ppen">
Zapai>ro Osorio,5M. R Chil4, A#R45"A, 56a, L5"> [NASA ADS]ppen"> [CrossRle]ppen"> [EDP S with D]ppen">Ke and aC CU. vera, A#R45"A, 537, A75"> [NASA ADS]ppen"> [CrossRle]ppen"> [EDP S with D]ppen">Koon0., S. E.d ma.30ApJ, 629, 1175"> [NASA ADS]ppen"> 86/431420">[CrossRle]ppen"> [NASA ADS]ppen"> 88

7X/719/1/890">[CrossRle]ppen">
Ke and aC CU. vera, in G/a>),-/sub><16-eAirborneeIaa28206-16pecklar iAa28onomy IV,nProc. nPIE, 8446, 84469">D. BonacciniclSPIEnoisf CSer., 4839, 140"> [NASA ADS]ppen"> [CrossRle]ppen"> [EDP S with D]ppen">Ford, H. C Chi02, ApJ, 578, 543"> [NASA ADS]ppen"> 86/342401">[CrossRle]ppen"> [NASA ADS]ppen"> [CrossRle]ppen"> [EDP S with D]ppen">Wai>rs, L. B. F. M. ma4, A#R45"A, 428, 663"> [NASA ADS]ppen"> [CrossRle]ppen"> [EDP S with D]ppen">Pacal,eE. vera, Nserva, 483, 64"> [NASA ADS]ppen"> 38/nserva10778">[CrossRle]ppen"> [NASA ADS]ppen">-band msr, ApJ, 8a8, L41"> [NASA ADS]ppen"> 88/ m41-8205/8a8/2/L41">[CrossRle]ppen"> [NASA ADS]ppen"> 38/nserva06002">[CrossRle]ppen"> class="sim, B. Rnd ms0, Dir hr>Imagicie820e mainly f (Tucssk AZ: Unervasvapor oArizona P2820),tid.aS.eSeagim, 111">TraubleW. A.d ma7, Nserva, 446, 771"> [NASA ADS]ppen"> 38/nserva05729 O[CrossRle]ppen"> [NASA ADS]ppen"> 86/503322">[CrossRle]ppen">Wei82, Gine.Chil4, ApJ, 783, L6"> [NASA ADS]ppen"> 88/ m41-8205/783/1/L6">[CrossRle]ppen"> [NASA ADS]ppen"> 88/

7X/768/1/14">[CrossRle]ppen">
[NASA ADS]ppen"> [CrossRle]ppen"> [EDP S with D]ppen"> [NASA ADS]ppen"> 86/17 39 O[CrossRle]ppen"> ticl">By suhaeituficiehat is, contrimg-inr te(20img src="R2012)) /fulehat is, contrimg-inr te(20img src="R1012)),nd (2p="/arzapeck degrle with errors is writtln asl.1. ExpecFD20206-16hat is, contrimg-equapeck(20img src="R(A.1)01/aa/f01/aa/fwhereehat is, contrimg-inr te(20img src="R) cereeused:l.1. Expec206-16hat is, contrimg-equapeck(20img src="R(A.8)01/aa/f01/aa/fwhereehat is, contrsimple-maph(20i>ρ ≪ 101/aa/f is assumed in d (2approximapeck . Iand (2 cotuceas linepeck,t.1. Expec206-16hat is, contrimg-equapeck(20img src="R(A.11)01/aa/f01/aa/fwhereeh1. Expec206-16hat is, contrimg-equapeck(20img src="R(A.12)01/aa/f01/aa/f (2underr tedepart is d (2error in hat is, contrsimple-maph(20i>P01/aa/f,nd at is, hat is, contrsimple-maph(2ℰp01/uh201/aa/f. Sis) ihat is, contrimg-inr te(20img src="R(A.13)01/aa/f01/aa/f n clwe c isalso writ ihat is, contrimg-inr te(20img src="R(A.14)01/aa/f01/aa/f

">
(2lWHM 820a2s(I′p01/aa/f ancehat is, contrsimple-maph(20i>P01/aa/f2s(λ(B.1)01/aa/f01/aa/fwhereehat is, contrsimple-maph(2Δθ01/aa/f2is d (2lWHM 820a2po/fu sour) iimage ancehat is, contrsimple-maph(20i>θ0/uh2p01/uh201/aa/f is d (2separapeckl820d (2pinly from d (2s"> .eEquapeck (.1.06" taRn from d (2s"> .eIand (2case 820d (2E-ELT (apererva dia clasehat is, contrsimple-maph(20i>D = 39l.s/aa/fm), we havel.at is, contrsimple-maph(2Δθ ≅ 0i>λD ≅ 0.006′′<1/aa/f. When hat is, contrsimple-maph(20i>θ0/uh2p01/uh201/aa/f is dakln asld (2innaseworkicieangl(2(IWA) 820d (2EPICS-IFS, hat is, contrsimple-maph(20.03′′<1/aa/f (.1.06" taR-band ms0012)) ar iexample, we havel.at is, contrsimple-maph(2Δλs( ≅ 0.22 0i>μ01/aa/fmclwhich is wid(2en-16. ulebeedis"iciuishedefrom d (2real2feserva with a2lWHM (.at is, contrsimple-maph(2Δλ<1/aa/f) 820.at is, contrsimple-maph(2~ 0.05 0i>μ01/aa/fm (.1.06" taR-band ms4012)). Howevim, when hat is, contrsimple-maph(20i>θ0/uh2p01/uh201/aa/f is hat is, contrsimple-maph(20.1′′<1/aa/f, hat is, contrsimple-maph(2Δλs(01/aa/f is hat is, contrsimple-maph(2~ 0.07 0i>μ01/aa/fmclwhich is assuarabla with .at is, contrsimple-maph(2Δλ<1/aa/f. As locieas0d (2s(P01/aa/f2s(P01/aa/f).

">
Iand is aclasdix,0d (2raw asn haste(.at is, contrsimple-maph(2C0/uh2raw01/uh201/aa/f)cenced (2dicmcteck asn haste(.at is, contrsimple-maph(2C0/uh2dic01/uh201/aa/f)c820a2high-asn hasteiaa28206-1 ereeexpin tede16-eformulatedeusiciewavefro-1 rippl(s 820d (2elea hic field.

".p> (2/fuensvapo(equival6-1 dond e numbma 820photck c-ma uni1 ereac-ma uni1 trti)c820s( hat is, contrsimple-maph(20i>θ01/aa/f anceposvaeck angl(2hat is, contrsimple-maph(20i>ψ<1/aa/f) encedrti2hat is, contrsimple-maph(20i>t<1/aa/f. S(φ01/aa/f c isb(2separaped indonresdome(.at is, contrsimple-maph(2φ0/uh2rand01/uh201/aa/f)cencequasi-s16pece(.at is, contrsimple-maph(2φ0/uh2qsc01/uh201/aa/f)cassuon6-1s.2hat is, contrsimple-maph(20i>φ0/uh2rand01/uh201/aa/f is mn tly2caused by0atmo clasalucurbulith cencev/ares resdomly2ck shor1 trti-scal6sc820a2few donden millrseasnds.2hat is, contrsimple-maph(20i>φ0/uh2qsc01/uh201/aa/f comas from ocEPICS impmaf hreck c820d (2iaa28206-1 enceevolvas oandrti-scal6sc820minut(s donhours.2hat is, contrsimple-maph(20i>φ0/uh2qsc01/uh201/aa/f c isb(2divid(d indontwocassuon6-1s: oae from anderation dona asronagraph (includicied (2csronagraph;2hat is, contrsimple-maph(20i>φ0/uh2qsc101/uh201/aa/f)cenced (2od ( iafcma d (2csronagraphe(.at is, contrsimple-maph(2φ0/uh2qsc201/uh201/aa/f;2see Fig. .1.06" taR1012)).

".p>Weiexam te d (2s"el"/a2elea hic field immees/pely2afcma d (2csronagraph, hat is, contrsimple-maph(2E0/uh2c01/uh201/aa/f. Heva, we dakl0a2shearicie/fueaf rome hic nu and as0a csronagraphear ia2simplecenalyEPICS discu pape, but d (2discu papesb(low c isb(2aclliedndonod ( ityp6sc820csronagraph . Iand (2nu and, incid(-1 ligh is separaped indontwocrays hat is, contrsimple-maph(2E0/uh2101/uh2 = 0i>Aeiφ0/uh2101/uh201/up>01/aa/f ancehat is, contrsimple-maph(20i>E0/uh2201/uh2 = 0i>Aei(φ0/uh2201/uh2 + 0i>π)01/up>01/aa/f, whereehat is, contrsimple-maph(20i>A01/aa/f2is a csns">n , ancehat is, contrsimple-maph(20i>φ0/uh2101/uh201/aa/f ancehat is, contrsimple-maph(20i>φ0/uh2201/uh201/aa/f are d (2phase perervbapeck clwhich are assumed ulebeesmall2(.at is, contrsimple-maph(2| 0i>φ0/uh2101/uh2 | , | 0i>φ0/uh2201/uh2 | ≪ 101/aa/f).2hat is, contrsimple-maph(20i>E0/uh2c01/uh201/aa/f c isb(2writtln asl.1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.1)01/aa/f01/aa/fwhereehat is, contrsimple-maph(2Δφ ≡ 0i>φ0/uh2201/uh2−0i>φ0/uh2101/uh201/aa/f. As for hat is, contrsimple-maph(20i>φ0/uh2c01/uh201/aa/f, we finceh1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.2)01/aa/f01/aa/fwhereewe sllit hat is, contrsimple-maph(2Δφ01/aa/f indonresdome(.at is, contrsimple-maph(2Δφ0/uh2rand01/uh201/aa/f)cencequasi-s16pece(.at is, contrsimple-maph(2Δφ0/uh2asc101/uh201/aa/f)cassuon6-1s.2When d (2s">peceperervbapeck2hat is, contrsimple-maph(20i>φ0/uh2qsc201/uh201/aa/f is added ulehat is, contrsimple-maph(20i>E0/uh2c01/uh201/aa/f, it becomas h1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.3)01/aa/f01/aa/f n clhat is, contrsimple-maph(20i>φ0/uh2qsc201/uh201/aa/f is negligible. The2elea hic field at d (2followicieimage pinly is .1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.4)01/aa/f01/aa/fwhereeFT()erepresn>) cd (2louri( iteas lin.

".p>Fr iak incid(-1 field hat is, contrsimple-maph(20i>Aeiφ01/up>01/aa/f, d (2rapecc820s(C01/aa/f, is given by0.1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.5)01/aa/f01/aa/f

".p>whereed (2peakscencev/ any c820d (2wavefro-1 rippl( havelv/ uas (length) hat is, contrsimple-maph(2± 0i>h01/aa/f (.1.06" taR class="simd ms0012),eEq. (121)). The2phase perervbapeck is datiribede1s .1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.6)01/aa/f01/aa/fwhereehat is, contrsimple-maph(20i>x01/aa/f2is a posvaeck at d (2pupil pinly ancehat is, contrsimple-maph(20i>x0/uh2001/uh201/aa/f is aitypPICS /aatiCS period 820d (2rippl(s.2hat is, contrsimple-maph(20i>a′<1/aa/f, hat is, contrsimple-maph(2α′<1/aa/f, hat is, contrsimple-maph(2b′<1/aa/f, ancehat is, contrsimple-maph(20i>β′<1/aa/f are csns">n s with uni1 c820cles/ns.2hat is, contrsimple-maph(20i>a′<1/aa/f is relatedeulehat is, contrsimple-maph(20i>h01/aa/f by0.1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.7)01/aa/f01/aa/fSee eqs. (108), (109)cenced (2relatedeuexteia Traub #R45"> class="simd(.1.06" taR".p>Iand (2case be0.g studied,0d (2raw asn hasteis .1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.8)01/aa/f01/aa/fwhereehat is, contrsimple-maph(20i>h0/uh2rand01/uh201/aa/f ancehat is, contrsimple-maph(20i>h0/uh2qsc101/uh201/aa/f areehat is, contrsimple-maph(20i>h01/aa/f for hat is, contrsimple-maph(2Δφ0/uh2rand01/uh201/aa/f ancehat is, contrsimple-maph(2Δφ0/uh2qsc101/uh201/aa/f, vas(
h0/uh2qsc101/uh201/aa/f c isb(2suppressed grearly2by0aede linebla mirror in d (2AO, it is saf donassumeehat is, contrsimple-maph(20i>h0/uh2rand01/uh2 ≫ 0i>h0/uh2qsc101/uh201/aa/f in d (iritypPICS v/ uas. Thn clwe finceh1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.9)01/aa/f01/aa/fC0/uh2raw01/uh201/aa/f is airapecc820ins">n nlyous2/fuensvares or photckirapes.2When w(2/fuegrape0photck cck a2dicmctor ar ia2cereain durapeck2hat is, contrsimple-maph(2Δt01/aa/f, d e effh0/uh2rand01/uh201/aa/f (.at is, contrsimple-maph(2h0/uh2rand,e01/uh201/aa/f)cwillsb(2h1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.10)01/aa/f01/aa/fwhereehat is, contrsimple-maph(20i>t0/uh2001/uh201/aa/f is aitypPICS drti-scal6 820.at is, contrsimple-maph(2Δφ0/uh2rand01/uh201/aa/fclwhich is a2few donden millrseasnds.2With a2suffi witt hat is, contrsimple-maph(2Δt01/aa/f, .at is, contrsimple-maph(2h0/uh2rand,e01/uh201/aa/f is much smallma d ak2hat is, contrsimple-maph(20i>h0/uh2qsc101/uh201/aa/f (i20.at is, contrsimple-maph(20i>h0/uh2rand01/uh2 ~ s00/up>-401/up>0i>λ<1/aa/f ancehat is, contrsimple-maph(20i>h0/uh2qsc101/uh2 ~ s00/up>-601/up>0i>λ<1/aa/f, d ek2hat is, contrsimple-maph(20i>h0/uh2rand,e01/uh2<0i>h0/uh2qsc101/uh201/aa/f for hat is, contrsimple-maph(2Δt> s0001/aa/f s). The2dicmcteck asn haste.at is, contrsimple-maph(2C0/uh2dic01/uh201/aa/f is daf tede1s d (2csn hasteacheeved with a2suffi witt hat is, contrsimple-maph(2Δt01/aa/f, which is h1. Expec206-16hat is, contrimg-equapeck(20img src="R(C.11)01/aa/f01/aa/fAccording donexperimn>) cby0.1.06" taRTraub ( ma7)012),eit is possible donreach .at is, contrsimple-maph(2C0/uh2dic01/uh2 = s00/up>-1001/up>01/aa/f in /aace. Iand is paper, .at is, contrsimple-maph(2C0/uh2dic01/uh201/aa/f is 6" erredefrom d (2v/ uas ar ia2g/a>),-/sub>-ban012) (.1.06" taR"ticl">< Iaput2v/ uas dond e calculateck .

".divs, contrin-txt(201.06" taR h/div>h/div> hdivs, contrinset(20divs, contrligte(2 .1.06" taR
List 820ainly f ar iwhich s( h/div>h/div> hh2s, contrsec O".1. Expecfigs(2012)AllsFigrvas uSchem>peceilluae npeckl820d (2hypod (EPICS observ6peckCS /netsa) anderation ;0b) anhalf-wave0aint ;0c) ex heti2adapsed iocEPIs;0d) ancsronagrapher ia2nu aicie/fueaf rome er;0e) an2/fuegral field s(f) anp="/arziciebeam sllit er;0ancehb>g) andicmctor.

"./tdO"./trO".tr>.tce,ol/aa/="2(201.06" target://dex 6r.Open with DEXTER 51/

61/ m16l8206s.">
Open with DEXTER6-16h/tdO./trO".tr>.tce,ol/aa/="2(s, contrin-txt(201.06" taR hdivs, contrinset(20teblaO".trO".tcev/ igt="middle(201.06" taR Schem>peceilluae npeckl820an2/mage profile bear a (up-ma r te)cenceafcma (lower r te)cd (2s("./tdO"./trO".tr>.tce,ol/aa/="2(201.06" target://dex 6r.Open with DEXTER 51/

61/ m16l8206s.">
Open with DEXTER6-16h/tdO./trO".tr>.tce,ol/aa/="2(s, contrin-txt(201.06" taR hdivs, contrinset(20teblaO".trO".tcev/ igt="middle(201.06" taR Bias in hat is, contrsimple-maph(20i>P01/aa/f with vas(C0/uh2dic01/uh2/0i>C0/uh2p/a01/uh2(0i>λ P01/aa/f. The2arrows2po/fu dond e maximum limit 820.at is, contrsimple-maph(20i>C0/uh2dic01/uh2/0i>C0/uh2p/a01/uh2(0i>λw01/aa/f2is setndonote d ird.

"./tdO"./trO".tr>.tce,ol/aa/="2(201.06" target://dex 6r.Open with DEXTER 51/

61/ m16l8206s.">
Open with DEXTER6-16h/tdO./trO".tr>.tce,ol/aa/="2(s, contrin-txt(201.06" taR hdivs, contrinset(20teblaO".trO".tcev/ igt="middle(201.06" taR hat is, contrimg-inr te(20img src="Rm∗01/up> =01/aa/f 0.0 (dot eder te), 5.0 (solid r te), 10.0 (dasheder te). Maglitudel.at is, contrsimple-maph(20i>m∗01/up>01/aa/f dicmrm tesehat is, contrimg-inr te(20img src="Rw01/aa/f2is setndonote d ird.

"./tdO"./trO".tr>.tce,ol/aa/="2(201.06" target://dex 6r.Open with DEXTER 51/

61/ m16l8206s.">
Open with DEXTER6-16h/tdO./trO".tr>.tce,ol/aa/="2(s, contrin-txt(201.06" taR hdivs, contrinset(20teblaO".trO".tcev/ igt="middle(201.06" taR MconcenceorbitCS distath Dl820d (2observ6ble ainly f. The2open cirsuppcere0ainly f observ6ble with a2s6-16cie820hat is, contrsimple-maph(20i>g = 1.001/aa/f ancehat is, contrsimple-maph(2Δt = 1501/aa/f h. The2filleder> g = 0.0501/aa/f ancehat is, contrsimple-maph(2Δt = 1501/aa/f h. The2filledesmallma cirsuppcere0those with a2s6-16cie820hat is, contrsimple-maph(20i>g = 0.0501/aa/f ancehat is, contrsimple-maph(2Δt = 5001/aa/f h.

"./tdO"./trO".tr>.tce,ol/aa/="2(201.06" target://dex 6r.Open with DEXTER 51/

61/ m16l8206s.">
Open with DEXTER6-16h/tdO./trO".tr>.tce,ol/aa/="2(s, contrin-txt(201.06" taR hdivs, contrinset(20teblaO".trO".tcev/ igt="middle(201.06" taR Cssuarisckl820csn hastf. Csn hastepmaflineth climited by0d (2s(P0i>C0/uh2p/a01/uh201/aa/f (filledecirsupp),nd at is, csn haste820d (2p="/arzedefeserva sigtbanda) an“typPICS”2pinly with averagedeproperei Dl820d (2observ6ble ainly f with .at is, contrsimple-maph(2g = 1.001/aa/f.0b) sati2asda) except ar i.at is, contrsimple-maph(2g = 0.0501/aa/f.0c) Gli De 876 b with .at is, contrsimple-maph(2g = 1.001/aa/f.0d) Gli De 876 b with .at is, contrsimple-maph(2g = 0.0501/aa/f.0The2d in solid r te in each aa/el is hat is, contrsimple-maph(2ΔP0i>C0/uh2p/a01/uh2(0i>θ)01/aa/f for a2pinly with d (2sati2clesuf. The2p="/arzapeck feserva is dacmct6ble in d (2rang Dl820hat is, contrsimple-maph(20i>θ01/aa/f desigtbted by0d (2arrows. 3-sigma .at is, contrsimple-maph(2C0/uh2dic01/uh201/aa/f (dashed-dot eder tes)cencehat is, contrsimple-maph(20i>C0/uh2raw01/uh201/aa/f (dot eder tes)cere also aiot ed ar id (26" erith .

"./tdO"./trO".tr>.tce,ol/aa/="2(201.06" target://dex 6r.Open with DEXTER 51/

61/ m16l8206s.">
Open with DEXTER6-16h/tdO./trO".tr>.tce,ol/aa/="2(s, contrin-txt(201.06" taR h/div> h/div> a hrCurrefu usage tichicsAbou1 ere sup tichicsReervn donare sup hdivs, contraa/els.data-ar ="tichics-siq"> h/div> hdivs, contraa/els.data-ar ="tichics-alm">h/div> hdivs, contraa/els.data-ar ="info">"

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.

" " .at is, contrico ico-info">01/aa/f" .div> hp>IaitiCS downloadl820d (2tichics may dakl0a2while.

" h/div> h/div> h/div> h/div>
h/div> ="articl
ticlticles/a> ticlticl ticlticles/ah/div> h/div> /ahfoocma> a ticlticles ticlticl a h/div> lehat is, contrteas lateck s.data- Expecseparapor_16-s.data-v/ uata 6-s>01/aa/f"hat is, contrteas lateck s.data- Expecshow_shor1_summarys.data-v/ uataShow shor1 summarys>01/aa/f"hat is, contrteas lateck s.data- Expechide_shor1_summarys.data-v/ uataHid(2shor1 summarys>01/aa/f"hat is, contrteas lateck s.data- Expeccookie_texts.data-v/ uataByeusicied is website, you agrle d at EDP S with D may storeeweb audieth cmeasuremn>) cookiencenc,2ck soti2pages, csokiencfrom sociCS ly works. .1.06" t'Rassuon6-1/cookie_p="icy/'>Moreeinfoinepeck and seerp<116s>01/aa/f"hat is, contrteas lateck s.data- Expeccookie_closes.data-v/ uataC"ickeulecloseed is notifPICaecks>01/aa/f"hat is, contrteas lateck s.data- Expecloadicis.data-v/ uataLoadici author...s>01/aa/f"hat is, contrteas lateck s.data- Expecview_fulltireens.data-v/ uataC"ickeuleview fulltireens>01/aa/f"hat is, contrteas lateck s.data- Expechide_fulltireens.data-v/ uataC"ickeanywhereedo hideid (2fulltireen 8verrays>01/aa/f"hat is, contrteas lateck s.data- Expecainy_movies.data-v/ uataC"ickeuleainy0d (2tovies>01/aa/f"hat is, contrteas lateck s.data- Expeccookie_is_mobiles.data-v/ uatas>01/aa/f"hat is, contrteas lateck s.data- Expecadd_to_ceres.data-v/ uataAddeuleceres>01/aa/f"hat is, contrteas lateck s.data- Expecremod _from_ceres.data-v/ uataRemod ifrom ceres>01/aa/f"hat is, contrteas lateck s.data- Expecloadici_ceres.data-v/ uataLoadici...s>01/aa/f"hat is, contrteas lateck s.data- Expecematy_ceres.data-v/ uataYour cart is ematys>01/aa/f""hat is, contrteas lateck s.data- Expecmakl6" _urls.data-v/ uataRassuon6-1/makl6" /?task=show#R45"typ6tntra#R45"doi=">01/aa/f"hat is, contrteas lateck s.data- Expecset_cere_ere sup_s">pus_urls.data-v/ uataRassuon6-1/servih D/?task=set_cere_ere sup_s">pus">01/aa/f"hat is, contrteas lateck s.data- Expec01/aa/f"hat is, contrteas lateck s.data- Expeccupar_shoppici_cere_urls.data-v/ uataRassuon6-1/servih D/?task=cupar_shoppici_ceres>01/aa/f"hat is, contrteas lateck s.data- Expecauthors_urls.data-v/ uataRassuon6-1/ere sup/?task=01/aa/f""hdivsid="waits.styla28disainy:note(20/div> hdivsstyla28disainy:note(20at isid="wait-phogr820s>01/aa/fh/div> /ahat isid="googup_enalyEPI s.data-hos nacds( data-domn tnacds.a16-0EXTE" data-t1="UA-39520934-1" data-t2="UA-5187131-2(s, contrhiddens>01/aa/f"haiript src="Rtemaint s/sour) /js/cookie.js(201tiript> haiript src="Rtemaint s/sour) /js/common.js(201tiript>