Open Access
Issue
A&A
Volume 663, July 2022
Article Number A55
Number of page(s) 12
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202142521
Published online 13 July 2022

© M. Lietzow and S. Wolf 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.

1 Introduction

Clouds are an essential part of atmospheres of most planets in the Solar System, such as sulfuric acid clouds on Venus (Young 1973), water clouds on Earth, or ammonia, ammonium hydro-sulfide, and water-ice clouds on Jupiter or Saturn (Atreya et al. 1999). Cloud particles can form and grow when supersaturation occurs. If the saturation vapor pressure curve of a cloud species crosses the planetary pressure-temperature profile, the particles condensate and form clouds (Sánchez-Lavega et al. 2004). Condensation of particles includes simple phase changes (e.g., water, ammonia) or thermochemical reactions, for example, in the case of ammonium hydrosulfide in the atmosphere of Jupiter (Lodders & Fegley 2002). In addition to Solar System objects, the existence of clouds is predicted (e.g., Burrows et al. 1997; Ackerman & Marley 2001; Helling et al. 2001; Lodders & Fegley 2002; Morley et al. 2012). Their signature has been observed in the atmosphere of brown dwarfs and extrasolar planets (e.g., Charbonneau et al. 2002; Pont et al. 2008; Sing et al. 2009; Demory et al. 2011; Gibson et al. 2013; Morley et al. 2013).

Cloud compositions do not only leave an imprint on the transmission spectra from the optical to infrared (Morley et al. 2014; Wakeford & Sing 2015), they also affect the planetary albedo and therefore the observed reflected flux as a function of the planetary phase angle (e.g., Seager et al. 2000). Furthermore, radiation scattered in the planetary atmosphere is usually polarized, whereas the net polarization of inactive solar-type stars is assumed to be negligible (Kemp et al. 1987; Cotton et al. 2017). Consequently, measuring the scattered polarized flux of a planet at various phase angles and wavelengths allows determining atmospheric properties. Hansen & Hovenier (1974) demonstrated this method by characterizing the size distribution and real part of the refractive index of the cloud particles in the atmosphere of Venus using polarimetric observations by Lyot (1929) and Coffeen & Gehrels (1969). Thus, polarime-try also represents a promising tool for the characterization of exoplanetary atmospheres.

With modern instruments, such as HIPPI (Bailey et al. 2015), HIPPI-2 (Bailey et al. 2020), POLISH2 (Wiktorowicz & Nofi 2015), POLLUX (Muslimov et al. 2018), or SPHERE (Beuzit et al. 2019), exoplanet polarimetry has come into reach. Various studies of scattered polarized radiation of exoplanets have been made. Selected examples are predictive studies of the reflected flux and polarization of close-in extrasolar giant planets (Seager et al. 2000), cloud-free Rayleigh scattering planetary atmospheres (Buenzli & Schmid 2009), Jupiter-like gas giants with homogeneous (Stam et al. 2004) and inhomogeneous cloud coverage (Karalidi et al. 2013), Earth-like exoplanets with homogeneous (Stam 2008; Karalidi et al. 2011) and inhomogeneous cloud coverage (Karalidi & Stam 2012), the retrieval of cloud coverage of Earth-like planets (Rossi & Stam 2017), and the polarization of self-luminous exoplanets (Stolker et al. 2017).

The first detection of polarized scattered radiation of an exo-planet was reported by Berdyugina et al. (2008). The authors were able to obtain polarimetric measurements of the HD 189733 system in B band. The system harbors a planet with an orbital semimajor axis of only about 0.03 AU and a period of about 2.2 days (Bouchy et al. 2005). Assuming Rayleigh scattering, the authors concluded that the observed modulation was due to scattered radiation from the atmosphere of HD 189733b. This first detection was also confirmed by Berdyugina et al. (2011). New observations of HD 189733 by Wiktorowicz & Nofi (2015) and Bott et al. (2016) also reported a linear polarization signal and confirmed the previously reported observations, but with a lower amplitude of polarization variations. Bailey et al. (2021), however, assumed that the highly variable polarization has its origin in the magnetic activity of the host star.

Another system for which linear polarization was reported by Bott et al. (2018) is WASP-18. This system harbors a planet with an orbital distance of about 0.02 AU and a period of 0.94 days (Hellier et al. 2009). Although the polarization is dominated by the additional influence of the interstellar medium, the authors set an upper limit of 40 ppm on the amplitude of reflected polarized radiation from WASP-18b.

Our study is organized as follows: in Sect. 2 we start with a brief overview of the applied methods. Our model of an exoplan-etary atmosphere is described in Sect. 3. Subsequently, the first goal of this study was to review characteristic properties of the polarization state of the radiation that is scattered once (hereafter: single-scattered) by cloud particles. For this purpose, in Sect. 4.1, we calculate optical properties of selected cloud particle compositions that are expected to condensate in various types of planetary atmospheres, ranging from ice giants to hot Jupiters. Our second goal was to investigate the scattered polarized flux of cloudy planetary atmospheres with various cloud compositions, which we report in Sect. 4.2. Therein, we compare the considered cloud compositions, discuss the polarization, and evaluate the feasibility of using characteristic features in the net polarization to distinguish between these cloud particles. To study the polarization at various phase angles and wavelengths ranging from 0.3 μm to 1 μm, the publicly available1 three-dimensional Monte Carlo radiative transfer code POLARIS (Reissl et al. 2016) was applied. It is well tested and has been used to investigate a broad range of astrophysical models such as molecular clouds (Reissl et al. 2017; Pellegrini et al. 2020; Seifried et al. 2020), Bok globules (Brauer et al. 2016; Zielinski et al. 2021), and protoplanetary disks (Heese et al. 2020; Brunngräber & Wolf 2020, 2021). For the purpose of our study, it has been optimized to handle the radiative transfer in planetary atmospheres (Lietzow et al. 2021). Finally, a discussion and a brief study of the observability of these polarization features is given in Sect. 5. Our conclusions are summarized in Sect. 6.

2 Methods

To calculate the scattered (polarized) flux of an exoplanet, we simulated the radiative transfer in its atmosphere with a Monte Carlo method. We used the three-dimensional Monte Carlo radiative transfer code POLARIS (Reissl et al. 2016), which has recently been optimized to calculate the scattered radiation of planetary atmospheres (Lietzow et al. 2021). Therein, the radiation field is represented by photon packages that are emitted by a spatially extended radiation source, propagate through the three-dimensional model space, are scattered by atmospheric particles, and are finally detected by an observer.

In the past, radiative transfer calculations applied to planetary atmospheres with the Monte Carlo approach have also been performed, for instance, by Buenzli & Schmid (2009), García Muñoz & Mills (2015), Emde et al. (2016), Stolker et al. (2017), or Bailey et al. (2018), while Stam et al. (2006) or Rossi et al. (2018), for example, performed calculations based on the adding-doubling method (de Haan et al. 1987).

To describe the state and degree of polarization of the photon packages, we used the Stokes formalism (e.g., Bohren & Huffman 1983). Here, each photon package carries a wavelength-dependent Stokes vector S = (I, Q, U, V)T. Its four components are the total flux I, the components Q and U, describing the linear polarization, and the circularly polarized flux V. At each interaction of the photon package with an atmospheric particle, the Stokes vector is multiplied by a scattering matrix F(θ, ϕ) to account for the change in polarization, SF(θ,ϕ)S',${\bf{S}} \propto {\bf{F}}\left( {\theta ,\phi } \right){\bf{S}}',$(1)

where S′ is the Stokes vector of the incoming radiation, and θ and ϕ are the scattering angles. When the photon package leaves the atmosphere, the resulting Stokes vector can be measured by an observer. The degree P1 and angle χ of linear polarization is P1=Q2+U2I,tan(2χ)=UQ.${P_1} = {{\sqrt {{Q^2} + {U^2}} } \over I},\quad \tan \left( {2\chi } \right) = {U \over Q}.$(2)

In studies of Solar System objects, a common approach is to express the degree of polarization in a reference system. In particular, the degree of polarization is defined as Ps=P1cos(2ϑr),ϑr=χ=(ψ±π2),${P_{\rm{s}}} = {P_1}\cos \left( {2{\vartheta _{\rm{r}}}} \right),\quad {\vartheta _{\rm{r}}} = \chi = \left( {\psi \pm {\pi \over 2}} \right),$(3)

where ψ is the position angle of the scattering plane, that is, a plane containing star, planet (or comet), and observer. The term inside the brackets must satisfy the condition 0 ≤ (ψ ± π/2) ≤ π (Chernova et al. 1993). As a consequence, the degree of polarization has a sign, and we can distinguish between positive and negative polarization, which means that the radiation is polarized perpendicular and parallel to the scattering plane, respectively. Furthermore, if the planet is symmetric with respect to the scattering plane, it is U = 0, and with ψ = 0, we define a signed degree of linear polarization such as Ps=QI.${P_{\rm{s}}} = - {Q \over I}.$(4)

In this study, the total flux only contains the net scattered planetary flux and not the stellar contribution. The planet can therefore be spatially resolved from its star.

In the case of a homogeneous spherical particle, the scattering matrix given in Eq. (1) is reduced to F=(F11F1200F12F110000F33F3400F34F33).${\bf{F}} = \left( {\matrix{{{F_{11}}} & {{F_{12}}} & 0 & 0 \cr {{F_{12}}} & {{F_{11}}} & 0 & 0 \cr 0 & 0 & {{F_{33}}} & {{F_{34}}} \cr 0 & 0 & { - {F_{34}}} & {{F_{33}}} \cr } } \right).$5

The four essential elements depend on the complex refractive index nc = nr + ini; of the particle and the size parameter x = 2πr/λ, where r/λ is the ratio of particle radius and wavelength.

For a collection of particles with a size distribution n(r), we define an averaged matrix element as Fij =rminrmaxn(r)Fij(r)  drrminrmaxn(r)  dr.$\left\langle {{F_{ij}}} \right\rangle = {{\int_{{r_{\min }}}^{{r_{\max }}} {n\left( r \right)} {F_{ij}}\left( r \right)\,\,{\rm{d}}r} \over {\int_{{r_{\min }}}^{{r_{\max }}} {n\left( r \right)} \,\,{\rm{d}}r}}.$(6)

Considering incident unpolarized radiation, the degree of polarization of single-scattered radiation by a collection of particles following Eqs. (1) and (4) becomes Ps= F12 F11 .${P_{\rm{s}}} = - {{\left\langle {{F_{12}}} \right\rangle } \over {\left\langle {{F_{11}}} \right\rangle }}.$7

3 Model setup

We used a spherical model space that contains the planet and the illuminating source (star). The planet is located at the center of coordinate space, and the star is located at a distance d from the planet. The planetary atmosphere was divided into a logarithmically spaced radial pressure grid, ranging from a top pressure of 10−5bar to a pressure of 100 bar. The atmosphere consisted of molecular hydrogen. The total scattering optical depth of the gaseous layers therefore was approximately 18 at 0.5 μm.

Assuming an absorbing surface underneath, the lower boundary represents a thick cloud layer that absorbs all incoming radiation. However, the impact of the lower boundary of the atmosphere is negligible, especially at shorter wavelengths, since the information of the polarization state is lost at this optical depth (Buenzli & Schmid 2009). For the gaseous particles, we only considered scattering and ignored absorption (see Sect. 5 for further discussion). The scattering cross-section, the refractive index, and the scattering matrix were calculated using the formula given by Sneep & Ubachs (2005), Cox (2000), and Hansen & Travis (1974), respectively. For molecular hydrogen, a depolarization factor of 0.02 was used (Hansen & Travis 1974).

In addition to the gaseous layers, we also considered clouds in the atmosphere. The clouds consist of small spherical particles whose distribution is characterized by a power law with exponential decay, such as n(r)rp1er/p2,$n\left( r \right) \propto {r^{{p_1}}}{{\rm{e}}^{{{ - r} \mathord{\left/ {\vphantom {{ - r} {{p_2}}}} \right. \kern-\nulldelimiterspace} {{p_2}}}}},$(8)

where r is the particle radius. Because of its simplicity, it is commonly applied to describe various types of clouds (Hansen 1971). This size distribution is sometimes also referred to as the standard (gamma) distribution (Hansen & Travis 1974). The two parameters p1 and p2 depend on the effective radius reff and effective variance veff of the distribution p1=13veffveff,p2=reffveff.${p_1} = {{1 - 3{v_{{\rm{eff}}}}} \over {{v_{{\rm{eff}}}}}},\quad {p_2} = {r_{{\rm{eff}}}}{v_{{\rm{eff}}}}.$(9)

In contrast to the gaseous particles, we considered both scattering and absorption by the cloud particles. The absorption and scattering cross-section as well as the scattering matrix of the cloud particles were determined using the program miex (Wolf & Voshchinnikov 2004), which is based on the Mie scattering theory (Mie 1908). An underlying assumption of this approach is that the particles are chemically homogeneous spheres.

Although cloud condensates in planetary atmospheres, such as water droplets, water-ice particles or solid particles in general, can be expected not to be perfectly spherical, the assumption of spherical particles can be used to mimic optical properties, such as cross-sections, of nonspherical particles (Kitzmann et al. 2010). However, polarization is much more sensitive to the particle shape. Linear polarization of randomly oriented nonspherical particles strongly depends on the aspect ratio (Mishchenko & Travis 1994). In addition, for size parameters ≳5 and large scattering angles, calculations assuming spheres usually do not agree with the measurements of nonspherical particles, but rounded particles are closely approximated by the Mie scattering theory (Perry et al. 1978). Furthermore, while liquid cloud condensates can be approximated by spheres, such as water clouds on Earth (Goloub et al. 1994) or sulfuric acid clouds on Venus (Hansen & Hovenier 1974; García Muñoz et al. 2014), solid particles, such as ice crystals in cirrus clouds, show different polarization signals compared to liquid water clouds due to their large variability in shape, size, and density (Hess 1998; Goloub et al. 2000). These ice particles in cirrus clouds, however, which were observed to be columns, rosettes, and bullet rosettes in shape, are usually several hundreds of micrometers in size, while smaller micrometer-sized particles appear to be spherical (Lawson et al. 1998).

For this study, we therefore assumed small spherical particles with p1 = 7 and p2 = 0.1 μm, defining the particle size distribution to apply the Mie scattering theory. These parameters correspond to an effective particle radius of 1 μm and an effective variance of 0.1 (see Eq. (9)). Probing larger and irregularly shaped particles, such as water-ice particles (see Karalidi et al. 2012), was beyond the scope of this paper.

The cloud layer was located between 1 bar and 0.1 bar. The choice of the size parameter and pressure values was based on previously motivated assumptions of micrometer-sized particles at pressures of several hundred millibar in cloudy atmospheres, for example, in the atmosphere of Jupiter (West et al. 1986) and Venus (Hansen & Hovenier 1974; Knibbe et al. 1997). The parameter values for the standard size distribution are also in agreement with the cloud model by Stam et al. (2004) for a Jupiter-like planet. In addition, in the case of HD 189733b, three-dimensional radiative-hydrodynamics with kinetic microphysi-cal mineral cloud formation processes show that mean cloud particle sizes are typically submicron to micron at pressures lower than 1 bar (Lee et al. 2016).

Pérez-Hoyos et al. (2005) presented a study of the vertical structure of clouds in the atmosphere of Saturn and found a strong latitudinal dependence on its optical thickness in the visible wavelength region, ranging from 20 to 45 at the equator, but 5 at the pole. Since the main goal of this study was to investigate the influence of different cloud compositions, and to avoid introducing too many variables, our planetary model had a homogenous cloud cover. The vertical optical depth of the cloud layer was chosen to be equal to 10 at 0.5 μm, which falls in the range assumed for the clouds of Saturn (Pérez-Hoyos et al. 2005), and also represents the approximate visible mean optical depth of the clouds of Jupiter as proposed by Sato & Hansen (1979).

Depending on the temperature profile in the atmosphere and its composition, clouds can condensate in various pressure regions (e.g., Lodders & Fegley 2002; Visscher et al. 2006, 2010), with condensates grown up to several hundreds of micrometers in size (Helling et al. 2008; Lee et al. 2015). However, investigating the latitudinal and longitudinal dependence on the height and optical depth of clouds as well as the dependence on the particle size and its influence on the reflected polarized flux were beyond the scope of the current study. Therefore, the planetary model discussed above served as a reference model for all cloud types, regardless of the composition.

The compositions considered here (see Table 1) are expected to condensate in different types of atmospheres of extrasolar planets and brown dwarfs (Kitzmann & Heng 2018). Methane (CH4), ammonia (NH3), or water (H2O) condensate at temperatures of about 80K to ~350 K (Carlson et al. 1988; Lodders & Fegley 2002). Sulfide or chloride clouds, such as zinc sulfide (ZnS), sylvite (KCl), or sodium sulfide (Na2S), condensate in the range of about 700 K to ~1100 K and are therefore expected in atmospheres of Y and T dwarfs (Morley et al. 2012). Iron, corundum (Al2O3), or silicate clouds, such as enstatite (MgSiO3) or forsterite (Mg2SiO4), usually condensate at high temperatures (≳1300K) in atmospheres of L dwarfs and deep in the atmospheres of T dwarfs (Lodders 2003; Visscher et al. 2010).

Clouds can condensate in either liquid or solid phase. In the case of iron and silicates, liquid clouds are usually found at much higher temperatures and pressures. For example, liquid iron forms at ≥1809K, liquid enstatite at ≥1851K, and liquid forsterite at ≥2163K (Lodders 1999). The liquid condensates may also dissolve other elements and have complex compositions (Lodders 1999). In addition, optical properties of these iron or silicates compositions exist mainly for the solid phase in the literature.

When we compare the refractive index of the solid and liquid states of CH4, H2O, or NH3, for instance, measurements show a difference in the real part of the refractive index of about 0.05 at 0.5 μm for CH4 (Martonchik & Orton 1994, interpolated), 0.03 at 0.5 μm for H2O (Hale & Querry 1973; Warren & Brandt 2008), or 0.06 for the sodium D line for NH3 (Marcoux 1969). Although these differences in the real part of the refractive index are minimal, they have an impact on the resulting optical properties, such as cross-sections and the scattering matrix. For H2O, the reflected radiation of Earth-like exoplanets with liquid and solid water clouds was studied, for example, by Karalidi et al. (2012). In this study, we assumed Uranus- and Neptune-like or Jupiter-like cloud compositions, that is, methane-ice clouds (Lindal et al. 1987, 1990; Lindal 1992) or ammonia-ice and water-ice clouds (Atreya et al. 1999) are expected, respectively. Furthermore, Morley et al. (2014) modeled atmospheres of Y dwarfs and exoplanets and found that water condenses in the upper atmosphere to form ice clouds. Thus, we investigated the case of solid cloud condensates in the atmosphere.

Table 1

Different cloud condensates, and the reference for the complex refractive index.

4 Results

Our goal was to determine the feasibility of distinguishing various cloud condensates expected in exoplanetary atmospheres using polarimetric observations. Therefore, we first reviewed the characteristic properties of the polarization resulting from single scattering by individual cloud particles. Subsequently, we investigated the scattered radiation of a three-dimensional planetary atmosphere consisting of gas and cloud layers as described in Sect. 3, taking multiple scattering events as well as absorption of the cloud particles and the ground layer into account.

4.1 Polarization characteristics of individual particles

As a preparatory study for the investigation of the polarization due to scattering in exoplanetary clouds, we analyzed and reviewed characteristic features that can be found in the polarization due to single scattering by potential cloud particles (Hansen & Travis 1974; Bailey 2007; Karalidi et al. 2011). For this purpose, we calculated the corresponding scattering matrix, based on the wavelength-dependent complex refractive index of the various condensates (see Table 1). Figure 1 shows the degree of polarization given by Eq. (7) for single-scattered radiation as a function of wavelength and scattering angle for various compositions. In general, materials with similar refractive indices also show a similar polarization pattern. Furthermore, materials with a high imaginary part of the refractive index (≳10−1) mainly show a positive polarization, while materials with a low imaginary part of the refractive index (≲10−3) also show a negative polarization.

For H2O and CH4 (see Fig. 1) with a real part of the refractive index of about 1.31 and 1.32, respectively, at a wavelength of 0.5 μm, there is a broad plateau of positive polarization at a scattering angle of approximately 150°. This is the result of rays that are reflected inside the particle once, also called the primary rainbow. This has been discussed, for example, by Liou & Hansen (1971), Hansen & Travis (1974), or Bailey (2007). The plateau of positive polarization ranges over the entire considered wavelength range. For scattering angles below 90°, the degree of polarization mainly has a negative sign, which is a consequence of rays that are refracted twice at the particle (Hansen & Travis 1974). At scattering angles of about 20°, there is a second plateau of positive polarization. This is a result of rays that are externally reflected once at the particle (Hansen & Travis 1974). This plateau ranges up to a wavelength of about 0.5 μm.

For increasing wavelengths, the rainbow is shifted toward shorter wavelengths and larger scattering angles (e.g., Liou & Hansen 1971; Hansen & Travis 1974; Bailey 2007). This is the case for NH3 with a real part of the refractive index of about 1.44 at a wavelength of 0.5 μm, or for KCl with a real part of the refractive index of about 1.5 at 0.5 μm. For increasing wavelengths, a plateau of negative polarization forms at large scattering angles of about 170° due to rays that are incident on the edge of the particle and internally reflected once as well (Hansen & Travis 1974). In addition, at small scattering angles of about 20°, the plateau of positive polarization ranges up to a wavelength of about 0.6 μm for NH3 and 0.7 μm for KCl. Similar features were discussed by Hansen & Hovenier (1974) in the context of the analysis of scattered light polarimetry of Venus. The real part of the refractive index of sulfuric acid clouds in the atmosphere of Venus is about 1.44 and is therefore similar to the real refractive index of ammonia.

This shift of the rainbow with increasing real part of the refractive index is valid up to a refractive index of about 1.547 (Nussenzveig 1969). Hence, for compositions such as MgAl2O4, MgO, or Fe2SiO4, the positive polarization shifts back to longer wavelengths, as also discussed by Hansen & Travis (1974). In addition, for a real part of the refractive index of ≳ 1.7, the plateau of positive polarization at small scattering angles ranges over the entire considered wavelength region. However, a peak of positive polarization at the end of this plateau, which is described by Hansen & Travis (1974) and is a result of anomalous diffraction (van de Hulst 1957), is not visible here.

For a real part of the refractive index >2, such as CaTiO3, TiO2, or ZnS, two plateaus of positive polarization form at large scattering angles. One plateau is located at a scattering angle of about 150° and ranges from a wavelength of 0.4 to ~0.7 μm. The second plateau is located at scattering angles of about 180° and ranges over the entire considered wavelength region.

For an increasing imaginary part of the refractive index, the amount of radiation traveling through the particle decreases. As a consequence, features such as negative polarization due to double refraction, or positive polarization due to internal reflections are lost, and radiation reflected from outside the particle with a positive degree of polarization dominates the polarization pattern (Hansen & Travis 1974). Thereby, the maximum degree of polarization is shifted toward 180°−2θB, where θB is the Brew-ster angle (Brewster 1815). This is the case for C, Cr, Fe, FeO, FeS, or Fe2O3. Furthermore, for compositions with a very large imaginary part of the refractive index (≳3), such as Fe or Cr, the degree of linear polarization of single-scattered radiation decreases compared to compositions such as FeO or FeS with an imaginary part of the refractive index of ≲1.4 at 0.5 μm. This applies to compositions with a high absolute value of the complex refractive index in general.

For compositions such as Al2O3, MgFeSiO4, MnS, or SiO, there is a high imaginary part of the refractive index at short wavelengths and a small imaginary part of the refractive index at long wavelengths. At short wavelengths, the scattered radiation is therefore dominated by a high degree of positive polarization, while with increasing wavelength, the polarization decreases and becomes negative at a characteristic wavelength.

thumbnail Fig. 1

Degree of signed polarization Ps of single-scattered radiation (see Eq. (7)) for various cloud types as a function of wavelength and scattering angle (θ = 0°: forward scattering). In addition, the approximate value of the complex refractive index at a wavelength of 0.5 μm is shown in the top left corner of each figure. See Sect. 4.1 for details.

4.2 Polarization of an exoplanetary atmosphere

Although it is possible to characterize the refractive index of the material by its single-scattered polarized radiation, as discussed in Sect. 4.1, multiple scattering and additional Rayleigh scattering by gaseous particles influence the polarized radiation scattered in a cloudy exoplanetary atmosphere (e.g., Karalidi et al. 2011). It therefore remained to be investigated whether the considered cloud compositions with their characteristic polarization features identified in Sect. 4.1 due to, for instance, a wavelength-dependent imaginary part of the refractive index, are distinguishable in this significantly more complex environment. Based on the model described in Sect. 3, we calculated the scattered polarized radiation of a cloudy planetary atmosphere. In Sects. 4.2.1 to 4.2.4, we investigated and compared different cloud compositions depending on the properties of their refractive index.

At shorter wavelengths (~0.3 μm), the degree of polarization for all atmospheric models is in general positive, with a maximum value around a phase angle of 90° due to Rayleigh scattering of the gaseous layers above the clouds, as discussed, for example, by Hansen & Hovenier (1974), Stam (2008), or Karalidi et al. (2011). However, with increasing wavelength, the impact of the gaseous particles weakens because of the strong wavelength dependence on the Rayleigh scattering cross-section (∝λ−4). Therefore, the degree of polarization of the cloud particles starts to dominate the net scattered radiation (Karalidi et al. 2011). As a result of multiple scattering in the atmosphere, the degree of polarization is lower than that of single scattering (Hansen & Travis 1974; Karalidi et al. 2011). However, although multiple scattering decreases the degree of polarization, it does not change the sign of polarization because the reflected radiation is dominated by single-scattered radiation from the upper layers of the clouds (Karalidi et al. 2011). Thus, the features discussed in Sect. 4.1 are still imprinted on the scattered polarized flux. In the following and in contrast to Fig. 1, we describe the scattered polarized flux as a function of the planetary phase angle α. For single-scattered radiation, it is α = 180° − θ.

4.2.1 Clouds with various real parts of the refractive index

In Fig. 2, the signed degree of linear polarization (Ps) of various model atmospheres is shown as a function of wavelength and phase angle. The model atmospheres have clouds consisting of H2O, NH3, MgSiO3, MgAl2O4, Fe2SiO4, or CaTiO3. These different cloud compositions have a small imaginary part of the refractive index (≲10−3 at 0.5 μm) in common and cover a broad range of the real part of the refractive index (~1.34 to ~2.36 at 0.5 μm). In this case, our simulations confirm that the characteristic polarization features due to the real part of the refractive index can be used to distinguish these different materials and to characterize the refractive index. Especially the locations of the positive degree of polarization at small phase angles (or at large scattering angles) have been discussed previously by various authors (e.g., Hansen & Hovenier 1974; Bailey 2007; Karalidi et al. 2011).

For clouds with a real part of the refractive index of about 1.55 to 1.75 such as MgSiO3 and MgAl2O4, the plateau of negative polarization at small phase angles (≲10°) is visible for the entire considered wavelength region. In addition, for cloud condensates with a real part of the refractive index of ≳1.75 such as MgAl2O4, Fe2SiO4, and CaTiO3, the plateau of the positive degree of polarization at large phase angles ranges over the entire considered wavelength region as well. For CaTiO3 clouds, there is a small spot of zero polarization at a phase angle of about 15° and at a wavelength of about 0.38 μm, surrounded by a positive degree of polarization. Furthermore, the degree of polarization increases at small wavelengths (~0.3 μm) for CaTiO3 because the imaginary part of the refractive index increases slightly at these wavelengths.

Although Karalidi et al. (2011) considered liquid-water clouds, the results of our simulations for water clouds are in agreement with those obtained by Karalidi et al. (2011) because the authors used similar parameter values for the size distribution and the difference in the refractive index comparing liquid and solid water is small (see Sect. 3). In contrast, in the study by Karalidi et al. (2012), the plateau of positive polarization at phase angles about 30° solely arises from internal reflections of liquid-water clouds and not from water-ice particles because the authors assumed nonspherical ice crystals.

The polarization curves of close-in extrasolar giant planets by Seager et al. (2000) or the polarization curves applied to HD 189733b by Kopparla et al. (2016) and Bailey et al. (2018) assumed cloud condensates, such as iron or silicates. Seager et al. (2000) considered a homogeneous MgSiO3-Al2O3-Fe cloud, and Kopparla et al. (2016) assumed a silicate cloud with a refractive index of 1.68 + 10−4i. The authors of these two studies observed a peak of polarization at small phase angles, which agrees with our peak of negative polarization at small phase angles. Bailey et al. (2018) compared clouds composed of MgSiO3, Mg2SiO4, Al2O3, or Fe. However, the authors considered particle sizes of reff ≤ 0.1 μm, which results in size parameters of ≲1.4. The scattering is therefore almost Rayleigh scattering. Thus, the polarization shows a maximum centered around a phase angle of 90°.

4.2.2 Clouds with a large imaginary part of the refractive index

In Fig. 3, the signed degree of linear polarization (Ps) is shown for model atmospheres with clouds consisting of FeS, Fe, FeO, Fe2O3, Cr, or C. These different cloud compositions cover a broad range of the real part of the refractive index (~1.36 to ~2.99 at 0.5 μm) and have a large imaginary part of the refractive index (≳0.7 at 0.5 μm) in common. Due to the high imaginary part of the refractive index, the single-scattering albedo of the cloud condensates decreases, resulting in a lower impact of multiple scattering and, thus, an increase in the degree of polarization. However, for Fe and Cr clouds, the polarization decreases because the degree of polarization after single scattering decreases (see Sect. 4.1). The maximum value of the degree of polarization is ~0.59 and ~0.47 for Fe and Cr clouds, respectively. For atmospheres composed of FeS, FeO and Fe2O3 clouds, the maximum value of the degree of polarization is ~0.75. In contrast to purely Rayleigh scattering, the maximum degree of polarization is shifted to larger phase angles for these compositions (see Sect. 4.1). At long wavelengths, the maximum degree of polarization is shifted to phase angles of about 125° for FeS, Fe, and Cr clouds, while for Fe2O3 and C clouds, the maximum is shifted to phase angles of about 115°, and to about 110° for FeO clouds. Since the imaginary part of the refractive index decreases at longer wavelengths for Fe2O3, there is a plateau of negative polarization for atmospheres composed of Fe2O3 at phase angles of about 30°.

thumbnail Fig. 2

Degree of signed polarization Ps (see Eq. (4)) as a function of wavelength and phase angle for various cloud compositions. As indicated by the spherical symbols at the top, a phase angle of 0° corresponds to the planet in full phase. See Sect. 4.2.1 for details.

thumbnail Fig. 3

Same as Fig. 2, but for different cloud compositions. See Sect. 4.2.2 for details.

thumbnail Fig. 4

Same as Fig. 2, but for different cloud compositions. See Sect. 4.2.3 for details.

4.2.3 Clouds with a variable imaginary part of the refractive index

In Fig. 4, the signed degree of linear polarization (Ps) is shown for model atmospheres with clouds consisting of Al2O3, MgFeSiO4, SiO, MnS, Na2S, or ZnS. These different cloud compositions cover a broad range of the real part of the refractive index (~1.58 to ~2.95 at 0.5 μm) and have a varying imaginary part of the refractive index in the considered wavelength region. For example, for MnS, the imaginary part of the refractive index decreases from about 1.3 at 0.3 μm to about 2.7 × 10−5 at 1 μm. In this case, the reflected radiation has solely a high degree of positive polarization at short wavelengths, but characteristic positive as well as negative polarization features at long wavelengths. Thus, at wavelengths ≳0.6 μm, the plateaus of positive and negative polarization determine the real part of the refractive index of the composition. At short wavelengths, the behavior of the positive polarization is linked to the imaginary part of the refractive index. In contrast to Al2O3 and MgFeSiO4 clouds, the positive polarization for an atmosphere composed of SiO, MnS, Na2S, and ZnS clouds decreases rapidly with increasing wavelength. At a characteristic wavelength of about 0.5 to 0.6 μm, depending on the material, the sign of the polarization changes. This is the result of the strong decrease of the imaginary part of the refractive index at these wavelengths. For atmospheres composed of Al2O3 and MgFeSiO4 clouds, no such characteristic wavelength can be defined. Instead, the sign of polarization changes at different wavelengths for different planetary phase angles. In addition, for Al2O3, MgFeSiO4, and Na2S clouds, the degree of negative polarization at long wavelengths is larger than for cloud compositions such as SiO, MnS, or ZnS. This is due to the larger imaginary part of the refractive index of Al2O3, MgFeSiO4, and Na2S at these wavelengths, resulting in a smaller single-scattering albedo and, thus, a lower impact of multiple scattering.

4.2.4 Clouds with an equal real part of the refractive index

In Fig. 5, the signed degree of linear polarization (Ps) is shown for model atmospheres with clouds consisting of KCl, MgSiO3, NaCl, SiO2, Al2O3, and Mg2SiO4, which have real part of the refractive index of about 1.50 to 1.62 at 0.5 μm. Except for Al2O3, the compositions have a small imaginary part of the refractive index (≲10−4 at 0.5 μm) as well. The plateau of positive polarization at large phase angles and the plateau of negative polarization at small phase angles show similar characteristics, regardless of the material. The most significant difference is found for an atmosphere consisting of Al2O3 clouds. The plateau of positive polarization at large phase angles and at long wavelengths reveals a real part of the refractive index, similar to Mg2SiO4 clouds. However, the higher degree of polarization for Al2O3 clouds compared to the other compositions indicates a smaller single-scattering albedo and, thus, a larger imaginary part of the refractive index.

Compositions with a similar behavior of real as well as imaginary part of the refractive index in the considered wavelength region, such as MgSiO3, NaCl, or SiO2, cannot be distinguished from each other. However, NaCl and MgSiO3 have different condensation temperatures of ~800 (Burrows & Sharp 1999) and ~1600K (Visscher et al. 2010), respectively. Thus, NaCl potentially condensates in atmospheres of T dwarfs, for example, while MgSiO3 will more likely be present in atmospheres of M dwarfs or deep in atmospheres of L dwarfs. In addition, as discussed by Visscher et al. (2010), if MgSiO3 condensates in the atmosphere, it removes silicon from the gas phase, preventing SiO2 (~1550K, Visscher et al. 2010) from condensating.

thumbnail Fig. 5

Same as Fig. 2, but for different cloud compositions. See Sect. 4.2.4 for details.

5 Discussion

In this section, we address selected parameters and corresponding questions that also have an impact on the net polarization, but were not covered in our simulations. The optical depth of the cloud layer in this study had a fixed value of 10 at 0.5 μm. An increasing optical depth of the cloud layer affects the degree of net polarization because the resulting increasing multiple scattering decreases the net degree of polarization (Karalidi et al. 2011; Bailey et al. 2018). On the other hand, a decreasing cloud optical depth would cause the Rayleigh scattering to dominate the net degree of polarization. This is similar to the case of a cloud layer that is located in higher pressure regions because the optical depth of gaseous particles increases above the cloud layer. Information about the polarization of deep cloud layers is lost at optical depths ≳2 of the upper layers (Buenzli & Schmid 2009). For a detailed discussion of cloud top pressure and cloud optical depth, see Karalidi et al. (2011).

The size distribution of the cloud particles is another crucial parameter determining the scattering properties of the particles. While the effective radius causes the degree of polarization to shift along the wavelength axis, increasing the effective variance causes various polarization features to smooth out (Hansen & Hovenier 1974; Hansen & Travis 1974). For smaller particles for a given observing wavelength, that is, smaller size parameters, the scattering is in the Rayleigh regime, resulting in a loss of the characteristic features in the polarization. For a detailed study of various particle sizes, see, for instance, Hansen & Hovenier (1974), Hansen & Travis (1974), Seager et al. (2000), Karalidi et al. (2011, 2012).

In our study, the planetary atmosphere had a homogeneous cloud coverage that covered the entire planet. However, the cloud coverage fraction and an asymmetric coverage of hemispheres have an impact on the scattered flux as well. The influence of the cloud coverage was studied in detail by Karalidi & Stam (2012), and Karalidi et al. (2013), and the retrieval of the cloud coverage of Earth-like exoplanets was studied, for example, by Rossi & Stam (2017). In addition, horizontally inhomogeneous cloud coverage also causes a nonzero circular polarization degree, which was studied by Rossi & Stam (2018).

As already mentioned in Sect. 3, we ignored absorption by gaseous particles. In general, absorption will decrease the effect of multiple scattering, thus, the degree of polarization increases. However, the net degree of polarization can decrease as well if there are, for instance, reflecting clouds in the lower layers. For a detailed discussion and radiative transfer computation with absorption by gaseous particles, see, for example, Stam et al. (2004) or Buenzli & Schmid (2009).

For the degree of polarization shown in Figs. 25, we assumed that the exoplanet and its host star are spatially separated observationally, thus, the high fraction of unpolarized stellar radiation is not included in the net polarization. Considering an exoplanet-star system where both components cannot be separated observationally, the ratio of the polarized planetary flux and the total stellar flux, that is, the polarization contrast, is given by (Hunziker et al. 2020) Cpol=P(α,λ)I(α,λ)r2d2.${C_{{\rm{pol}}}} = P\left( {\alpha ,\lambda } \right)I\left( {\alpha ,\lambda } \right){{{r^2}} \over {{d^2}}}.$10

Here, P and I are the phase angle and wavelength-dependent degree of polarization and reflectivity of the planet, respectively. At a phase angle of 0°, the reflectivity is equivalent to the planetary geometric albedo. Furthermore, r and d are the planetary radius and distance between planet and star, respectively. To calculate the planetary temperature profile, we followed Guillot (2010) and took parameter values based on HD 209458b (κth = 10−2 cm2 g−1, κv = 4 × 10−3 cm2 g−1, Tint = 500 K, Guillot 2010; r = 1.45 Rj, T = 6092 K, Boyajian et al. 2015; m = 0.64 Mj, Snellen et al. 2010). However, we set d = 1 AU with an observational distance of 5 pc to be suitable for SPHERE/ZIMPOL (Hunziker et al. 2020). Thus, we obtained temperatures of about 500 to ~1000K in the pressure region of about 10mbar to ~1 bar, which means that compositions such as KCl, NaCl, Na2S, or ZnS potentially form cloud condensates.

Figure 6 shows the polarization contrast (Eq. (10)) as a function of wavelength and planetary phase angle for various compositions that are expected to condensate in the atmosphere, and a cloud-free atmosphere for comparison. For the cloud-free atmosphere, the polarization contrast is highest with maximum values of about 3.7 × 10−8 at phase angles of about 65° across the considered wavelength region. For the cloudy atmospheres, the polarization contrast is highest at short wavelengths (~0.35 μm) with values of ~2 × 10−8 due to Rayleigh scattering. In addition, for NaCl clouds, the polarization contrast reaches values of ≳2 × 10−8 at small phase angles and long wavelengths, making the plateau of negative polarization detectable, for instance, with SPHERE/ZIMPOL (Hunziker et al. 2020). For Na2S and ZnS clouds, the high degree of polarization at short wavelength with a rapidly decreasing degree of polarization for increasing wavelength, as discussed in Sect. 4.2.3, produces a polarization contrast of ≲108 because the total reflected flux decreases with decreasing single-scattering albedo as well.

thumbnail Fig. 6

Polarization contrast Cpol (see Eq. (10)) as a function of wavelength and phase angle for various cloud compositions and a cloud-free atmosphere (top left). As indicated by the spherical symbols at the top, a phase angle of 0° corresponds to the planet in full phase. See Sect. 5 for details.

6 Conclusions

In the first part of this study (Sect. 4.1), we investigated the single-scattered polarized radiation of various cloud condensates, which are expected in planetary atmospheres. We reviewed characteristic properties of the polarization state depending on the scattering angle, wavelength, and refractive index of the cloud particle. As demonstrated in selected earlier studies, we confirmed the significant impact of the chemical composition and, thus, the complex refractive index of the cloud particles on the state of polarization of the scattered radiation.

If the imaginary part of the refractive index is small, refracted rays dominate the polarization pattern (Hansen & Travis 1974). Characteristic features, such as the rainbow, depend on the real part of the refractive index (e.g., Liou & Hansen 1971; Hansen & Travis 1974; Bailey 2007). However, most of the considered cloud condensates have a nonnegligible imaginary part. For an increasing imaginary part of the refractive index, polarization features resulting from rays that are refracted or reflected inside the particle are lost, and radiation that is reflected externally dominates (Hansen & Travis 1974). For compositions with a large imaginary part of the refractive index over the entire considered wavelength region, the polarization is positive, with the maximum value shifted to smaller scattering angles (see Hansen & Travis 1974). Compositions such as Al2O3, MgFeSiO4, SiO, MnS, Na2S, and SiO have a strongly varying imaginary part of the refractive index. At short wavelengths, the polarization is therefore mainly positive and becomes negative at longer wavelengths. In between, there is a change of sign of the polarization due to the decrease of the imaginary part. The wavelength at which the sign changes depends on the chemical composition and, thus, allows for the identification of the chemical composition of the respective condensate.

The goal of the second part of this study (Sect. 4.2) was to evaluate the feasibility to distinguish between the different considered cloud materials in a planetary atmosphere. For this purpose, we simulated the reflected polarized flux of cloudy planetary atmospheres with selected cloud compositions investigated in Sect. 4.1. These simulations were performed with the three-dimensional Monte Carlo radiative transfer code POLARIS, which was optimized to handle the radiative transfer in planetary atmospheres. The degree of linear polarization was studied at planetary phase angles ranging from 0° to 180° at wavelengths from 0.3 μm to 1 μm.

Building on the results of the single-scattering polarization obtained in the first part of our study, we first presented results for compositions with a small imaginary, but various real parts of the refractive index in Sect. 4.2.1. In agreement with earlier studies (Hansen & Hovenier 1974), our simulations confirmed that the polarization features such as the rainbow can be used to determine the refractive index of the cloud particles.

For high values of the imaginary part of the refractive index, considered in Sect. 4.2.2, the polarization is mainly positive and the overall degree of polarization increases. However, each of the considered cloud composition shows a characteristic behavior of the polarization, depending on the individual imaginary part of the refractive index.

  • Cloud compositions with an imaginary part of ≳3, such as Fe and Cr, show a decrease in the degree of polarization compared to compositions such as FeS and FeO.

  • For FeS, Fe, and Cr clouds, the maximum degree of polarization at long wavelengths is shifted to phase angles of about 125°, while for FeO clouds, the maximum degree of polarization is shifted to phase angles of about 110°.

  • For Fe2O3, there is a plateau of negative polarization at wavelengths of ≳0.9 μm because the imaginary part of the refractive index is decreasing rapidly with increasing wavelength.

If the imaginary part of the refractive index strongly varies in the considered wavelength region (Sect. 4.2.3), the behavior of the degree of polarization and a change of the sign of the polarization degree can be linked to the behavior of the imaginary part of the refractive index of the condensate, allowing one to characterize the cloud composition.

  • Clouds composed of Al2O3 and MgFeSiO4 show a smooth decrease of positive polarization with increasing wavelength.

  • Clouds composed of SiO, MnS, Na2S, and ZnS show a rapidly decreasing degree of polarization with increasing wavelength. In addition, at a characteristic wavelength of about 0.5 to 0.6 μm, depending on the imaginary part of the refractive index, the sign of the polarization changes.

Finally, for condensates with similar values of the real part of the refractive index (Sect. 4.2.4), a distinction of various compositions is almost impossible.

  • The most significant difference in the polarization was identified in the case of Al2O3 clouds. Here, the underlying imaginary part of the refractive index results in a high degree of positive and negative polarization at short and long wavelengths, respectively.

  • However, chlorides or sulfides with similar refractive indices as silicates have a different condensation temperature. Thus, compositions with a similar net polarization across the considered wavelength and phase angle region, but different condensation temperatures, can be distinguished indirectly by comparison of the temperatures of the planetary atmosphere and, thus, its potential to host certain condensates/clouds.

In summary, most of the considered compositions that are expected to condensate in exoplanetary atmospheres, such as chlorides, sulfides, or silicates, can be distinguished via polarization measurements due to their unique wavelength-dependent real as well as imaginary part of the refractive index and/or different condensation temperatures.

Acknowledgements

We wish to thank the anonymous referee, for suggestions improving the presentation of the results of this study. This research made use of Astropy (https://www.astropy.org), a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018), Matplotlib (https://matplotlib.org) (Hunter 2007), Numpy (https://numpy.org) (Harris et al. 2020), NASA’s Astrophysics Data System (https://ui.adsabs.harvard.edu), and a modified A&A bibliography style file with clickable link in the bibliography2.

References

  1. Ackerman, A. S., & Marley, M. S. 2001, ApJ, 556, 872 [Google Scholar]
  2. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  4. Atreya, S. K., Wong, M. H., Owen, T. C., et al. 1999, Planet. Space Sci., 47, 1243 [CrossRef] [Google Scholar]
  5. Bailey, J. 2007, Astrobiology, 7, 320 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bailey, J., Kedziora-Chudczer, L., Cotton, D. V., et al. 2015, MNRAS, 449, 3064 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bailey, J., Kedziora-Chudczer, L., & Bott, K. 2018, MNRAS, 480, 1613 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bailey, J., Cotton, D. V., Kedziora-Chudczer, L., De Horta, A., & Maybour, D. 2020, PASA, 37, e004 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bailey, J., Bott, K., Cotton, D. V., et al. 2021, MNRAS, 502, 2331 [NASA ADS] [CrossRef] [Google Scholar]
  10. Berdyugina, S. V., Berdyugin, A. V., Fluri, D. M., & Piirola, V. 2008, ApJ, 673, L83 [NASA ADS] [CrossRef] [Google Scholar]
  11. Berdyugina, S. V., Berdyugin, A. V., Fluri, D. M., & Piirola, V. 2011, ApJ, 728, L6 [NASA ADS] [CrossRef] [Google Scholar]
  12. Beuzit, J. L., Vigan, A., Mouillet, D., et al. 2019, A&A, 631, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Bohren, C. F., & Huffman, D. R. 1983, Absorption and scattering of light by small particles (Hoboken: John Wiley & Sons) [Google Scholar]
  14. Bott, K., Bailey, J., Kedziora-Chudczer, L., et al. 2016, MNRAS, 459, L109 [NASA ADS] [CrossRef] [Google Scholar]
  15. Bott, K., Bailey, J., Cotton, D. V., et al. 2018, AJ, 156, 293 [NASA ADS] [CrossRef] [Google Scholar]
  16. Bouchy, F., Udry, S., Mayor, M., et al. 2005, A&A, 444, L15 [EDP Sciences] [Google Scholar]
  17. Boyajian, T., von Braun, K., Feiden, G. A., et al. 2015, MNRAS, 447, 846 [Google Scholar]
  18. Brauer, R., Wolf, S., & Reissl, S. 2016, A&A, 588, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Brewster, D. 1815, Phil. Trans. R. Soc., 105, 125 [CrossRef] [Google Scholar]
  20. Brunngräber, R., & Wolf, S. 2020, A&A, 640, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Brunngräber, R., & Wolf, S. 2021, A&A, 648, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Buenzli, E., & Schmid, H. M. 2009, A&A, 504, 259 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Burrows, A., & Sharp, C. M. 1999, ApJ, 512, 843 [NASA ADS] [CrossRef] [Google Scholar]
  24. Burrows, A., Marley, M., Hubbard, W. B., et al. 1997, ApJ, 491, 856 [Google Scholar]
  25. Carlson, B. E., Rossow, W. B., & Orton, G. S. 1988, J. Atmos. Sci. 45, 2066 [NASA ADS] [CrossRef] [Google Scholar]
  26. Charbonneau, D., Brown, T. M., Noyes, R. W., & Gilliland, R. L. 2002, ApJ, 568, 377 [Google Scholar]
  27. Chernova, G. P., Kiselev, N. N., & Jockers, K. 1993, Icarus, 103, 144 [Google Scholar]
  28. Coffeen, D. L., & Gehrels, T. 1969, AJ, 74, 433 [NASA ADS] [CrossRef] [Google Scholar]
  29. Cotton, D. V., Marshall, J. P., Bailey, J., et al. 2017, MNRAS, 467, 873 [NASA ADS] [Google Scholar]
  30. Cox, A. N. 2000, Allen's Astrophysical Quantities (Berlin: Springer) [Google Scholar]
  31. de Haan, J. F., Bosma, P. B., & Hovenier, J. W. 1987, A&A, 183, 371 [Google Scholar]
  32. Demory, B., Seager, S., Madhusudhan, N., et al. 2011, ApJ, 735, L12 [NASA ADS] [CrossRef] [Google Scholar]
  33. Dorschner, J., Begemann, B., Henning, Th., Jäger, C., & Mutschke, H. 1995, A&A, 300, 503 [Google Scholar]
  34. Draine, B. T. 2003, ApJ, 598, 1017 [NASA ADS] [CrossRef] [Google Scholar]
  35. Emde, C., Buras-Schnell, R., Kylling, A., et al. 2016, Geophys. Mod. Dev., 9, 1647 [CrossRef] [Google Scholar]
  36. Fabian, D., Henning, Th., Jäger, C., et al. 2001, A&A, 378, 228 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. García Muñoz, A., & Mills, F. P. 2015, A&A, 573, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. García Muñoz, A., Pérez-Hoyos, S., & Sánchez-Lavega, A. 2014, A&A, 566, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Gibson, N. P., Aigrain, S., Barstow, J. K., et al. 2013, MNRAS, 436, 2974 [CrossRef] [Google Scholar]
  40. Goloub, P., Deuze, J. L., Herman, M., & Fouquart, Y. 1994, IEEE Trans. Geosci. Remote Sens., 32, 78 [NASA ADS] [CrossRef] [Google Scholar]
  41. Goloub, P., Herman, M., Chepfer, H., et al. 2000, J. Geophys. Res., 105, 14747 [NASA ADS] [CrossRef] [Google Scholar]
  42. Guillot, T. 2010, A&A, 520, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Hale, G. M., & Querry, M. R. 1973, Appl. Opt., 12, 555 [NASA ADS] [CrossRef] [Google Scholar]
  44. Hansen, J. E. 1971, J. Atmos. Sci., 28, 1400 [NASA ADS] [CrossRef] [Google Scholar]
  45. Hansen, J. E., & Hovenier, J. W. 1974, J. Atmos. Sci., 31, 1137 [NASA ADS] [CrossRef] [Google Scholar]
  46. Hansen, J. E., & Travis, L. D. 1974, Space Sci. Rev., 16, 527 [Google Scholar]
  47. Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357 [NASA ADS] [CrossRef] [Google Scholar]
  48. Heese, S., Wolf, S., & Brauer, R. 2020, A&A, 634, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Hellier, C., Anderson, D. R., Collier Cameron, A., et al. 2009, Nature, 460, 1098 [NASA ADS] [CrossRef] [Google Scholar]
  50. Helling, Ch., Oevermann, M., Lüttke, M. J. H., Klein, R., & Sedlmayr, E. 2001, A&A, 376, 194 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Helling, Ch., Woitke, P., & Thi, W. 2008, A&A, 485, 547 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Henning, Th., Begemann, B., Mutschke, H., & Dorschner, J. 1995, A&AS, 112, 143 [NASA ADS] [Google Scholar]
  53. Hess, M. 1998, J. Quant. Spectr. Rad. Transf., 60, 301 [NASA ADS] [CrossRef] [Google Scholar]
  54. Huffman, D. R., & Wild, R. L. 1967, Phys. Rev., 156, 989 [NASA ADS] [CrossRef] [Google Scholar]
  55. Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [NASA ADS] [CrossRef] [Google Scholar]
  56. Hunziker, S., Schmid, H. M., Mouillet, D., et al. 2020, A&A, 634, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Jäger, C., Dorschner, J., Mutschke, H., Posch, T., & Henning, Th. 2003, A&A, 408, 193 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Karalidi, T., & Stam, D. M. 2012, A&A, 546, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Karalidi, T., Stam, D. M., & Hovenier, J. W. 2011, A&A, 530, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Karalidi, T., Stam, D. M., & Hovenier, J. W. 2012, A&A, 548, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Karalidi, T., Stam, D. M., & Guirado, D. 2013, A&A, 555, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Kemp, J. C., Henson, G. D., Steiner, C. T., & Powell, E. R. 1987, Nature, 326, 270 [NASA ADS] [CrossRef] [Google Scholar]
  63. Khachai, H., Khenata, R., Bouhemadou, A., et al. 2009, J. Phys. Condens. Matter, 21, 095404 [NASA ADS] [CrossRef] [Google Scholar]
  64. Kitzmann, D., & Heng, K. 2018, MNRAS, 475, 94 [NASA ADS] [CrossRef] [Google Scholar]
  65. Kitzmann, D., Patzer, A. B. C., von Paris, P., et al. 2010, A&A, 511, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Knibbe, W. J. J., de Haan, J. F., Hovenier, J. W., & Travis, L. D. 1997, J. Geophys. Res., 102, 10945 [NASA ADS] [CrossRef] [Google Scholar]
  67. Koike, C., Kaito, C., Yamamoto, T., et al. 1995, Icarus, 114, 203 [NASA ADS] [CrossRef] [Google Scholar]
  68. Kopparla, P., Natraj, V., Zhang, Xi., et al. 2016, ApJ, 817, 32 [NASA ADS] [CrossRef] [Google Scholar]
  69. Lawson, R. P., Heymsfield, A. J., Aulenbach, S. M., & Jensen, T. L. 1998, Geophys. Res. Lett., 25, 1331 [NASA ADS] [CrossRef] [Google Scholar]
  70. Lee, E., Helling, Ch., Dobbs-Dixon, I., & Juncher, D. 2015, A&A, 580, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Lee, E., Dobbs-Dixon, I., Helling, Ch., Bognar, K., & Woitke, P. 2016, A&A, 594, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Lietzow, M., Wolf, S., & Brunngräber, R. 2021, A&A, 645, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Lindal, G. F. 1992, AJ, 103, 967 [Google Scholar]
  74. Lindal, G. F., Lyons, J. R., Sweetnam, D. N., et al. 1987, J. Geophys. Res., 92, 14987 [NASA ADS] [CrossRef] [Google Scholar]
  75. Lindal, G. F., Lyons, J. R., Sweetnam, D. N., et al. 1990, Geophys. Res. Lett., 17, 1733 [NASA ADS] [CrossRef] [Google Scholar]
  76. Liou, K., & Hansen, J. E. 1971, J. Atmos. Sci., 28, 995 [NASA ADS] [CrossRef] [Google Scholar]
  77. Lodders, K. 1999, ApJ, 519, 793 [NASA ADS] [CrossRef] [Google Scholar]
  78. Lodders, K. 2003, ApJ, 591, 1220 [Google Scholar]
  79. Lodders, K., & Fegley, B. 2002, Icarus, 155, 393 [Google Scholar]
  80. Lyot, B. 1929, Ann. Obs. Meudon, 8, 1 [Google Scholar]
  81. Marcoux, J. E. 1969, J. Opt. Soc. Am., 59, 998 [NASA ADS] [CrossRef] [Google Scholar]
  82. Martonchik, J. V., & Orton, G. S. 1994, Appl. Opt., 33, 8306 [NASA ADS] [CrossRef] [Google Scholar]
  83. Martonchik, J. V., Orton, G. S., & Appleby, J. F. 1984, Appl. Opt., 23, 541 [NASA ADS] [CrossRef] [Google Scholar]
  84. Mie, G. 1908, Ann. Phys., 330, 377 [Google Scholar]
  85. Mishchenko, M., & Travis, L. D. 1994, Appl. Opt., 33, 7206 [NASA ADS] [CrossRef] [Google Scholar]
  86. Morley, C. V., Fortney, J. J., Marley, M. S., et al. 2012, ApJ, 756, 172 [NASA ADS] [CrossRef] [Google Scholar]
  87. Morley, C. V., Fortney, J. J., Kempton, E. M. R., et al. 2013, ApJ, 775, 33 [CrossRef] [Google Scholar]
  88. Morley, C. V., Marley, M. S., Fortney, J. J., et al. 2014, ApJ, 787, 78 [Google Scholar]
  89. Muslimov, E., Bouret, J., Neiner, C., et al. 2018, Proc. SPIE, 10699, 1069906 [NASA ADS] [Google Scholar]
  90. Nussenzveig, H. M. 1969, J. Math. Phys., 10, 125 [NASA ADS] [CrossRef] [Google Scholar]
  91. Palik, E. D. 1985, Handbook of Optical Constants of Solids (Cambridge: Academic Press) [Google Scholar]
  92. Palik, E. D. 1991, Handbook of Optical Constants of Solids II (Cambridge: Academic Press) [Google Scholar]
  93. Pellegrini, E. W., Reissl, S., Rahner, D., et al. 2020, MNRAS, 498, 3193 [NASA ADS] [CrossRef] [Google Scholar]
  94. Pérez-Hoyos, S., Sánchez-Lavega, A., French, R. G., & Rojas, J. F. 2005, Icarus, 176, 155 [CrossRef] [Google Scholar]
  95. Perry, R. J., Huffman, D. R., & Hunt, A. J. 1978, Appl. Opt., 17, 2700 [NASA ADS] [CrossRef] [Google Scholar]
  96. Pollack, J. B., Hollenbach, D., Beckwith, S., et al. 1994, ApJ, 421, 615 [Google Scholar]
  97. Pont, F., Knutson, H., Gilliland, R. L., Moutou, C., & Charbonneau, D. 2008, MNRAS, 385, 109 [Google Scholar]
  98. Reissl, S., Wolf, S., & Brauer, R. 2016, A&A, 593, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Reissl, S., Seifried, D., Wolf, S., Banerjee, R., & Klessen, R. S. 2017, A&A, 603, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Rossi, L., & Stam, D. M. 2017, A&A, 607, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  101. Rossi, L., & Stam, D. M. 2018, A&A, 616, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  102. Rossi, L., Berzosa-Molina, J., & Stam, D. M. 2018, A&A, 616, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Sánchez-Lavega, A., Pérez-Hoyos, S., & Hueso, R. 2004, Am. J. Phys., 72, 767 [Google Scholar]
  104. Sato, M., & Hansen, J. E. 1979, J. Atmos. Sci., 36, 1133 [NASA ADS] [CrossRef] [Google Scholar]
  105. Seager, S., Whitney, B. A., & Sasselov, D. D. 2000, ApJ, 540, 504 [NASA ADS] [CrossRef] [Google Scholar]
  106. Seifried, D., Walch, S., Weis, M., et al. 2020, MNRAS, 497, 4196 [Google Scholar]
  107. Siefke, T., Kroker, S., Pfeiffer, K., et al. 2016, Adv. Opt. Mater., 4, 1780 [CrossRef] [Google Scholar]
  108. Sing, D. K., Désert, J. M., Lecavelier Des Etangs, A., et al. 2009, A&A, 505, 891 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  109. Sneep, M., & Ubachs, W. 2005, J. Quant. Spectr. Rad. Transf., 92, 293 [Google Scholar]
  110. Snellen, I. A. G., de Kok, R. J., de Mooij, E. J. W., & Albrecht, S. 2010, Nature, 465, 1049 [Google Scholar]
  111. Stam, D. M. 2008, A&A, 482, 989 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  112. Stam, D. M., Hovenier, J. W., & Waters, L. B. F. M. 2004, A&A, 428, 663 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  113. Stam, D. M., de Rooij, W. A., Cornet, G., & Hovenier, J. W. 2006, A&A, 452, 669 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  114. Stolker, T., Min, M., Stam, D. M., et al. 2017, A&A, 607, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  115. Ueda, K., Yanagi, H., Noshiro, R., Hosono, H., & Kawazoe, H. 1998, J. Phys. Condens. Matter, 10, 3669 [NASA ADS] [CrossRef] [Google Scholar]
  116. van de Hulst, H. C. 1957, Light Scattering by Small Particles (Hoboken: John Wiley & Sons) [Google Scholar]
  117. Visscher, C., Lodders, K., & Fegley, B. 2006, ApJ, 648, 1181 [NASA ADS] [CrossRef] [Google Scholar]
  118. Visscher, C., Lodders, K., & Fegley, B. 2010, ApJ, 716, 1060 [NASA ADS] [CrossRef] [Google Scholar]
  119. Wakeford, H. R., & Sing, D. K. 2015, A&A, 573, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  120. Warren, S. G., & Brandt, R. E. 2008, J. Geophys. Res., 113, D14220 [NASA ADS] [CrossRef] [Google Scholar]
  121. West, R. A., Strobel, D. F., & Tomasko, M. G. 1986, Icarus, 65, 161 [NASA ADS] [CrossRef] [Google Scholar]
  122. Wiktorowicz, S. J., & Nofi, L. A. 2015, ApJ, 800, L1 [NASA ADS] [CrossRef] [Google Scholar]
  123. Wolf, S., & Voshchinnikov, N. V. 2004, Comput. Phys. Commun., 162, 113 [Google Scholar]
  124. Young, A. T. 1973, Icarus, 18, 564 [NASA ADS] [CrossRef] [Google Scholar]
  125. Zeidler, S., Posch, T., Mutschke, H., Richter, H., & Wehrhan, O. 2011, A&A, 526, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  126. Zielinski, N., Wolf, S., & Brunngräber, R. 2021, A&A, 645, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Different cloud condensates, and the reference for the complex refractive index.

All Figures

thumbnail Fig. 1

Degree of signed polarization Ps of single-scattered radiation (see Eq. (7)) for various cloud types as a function of wavelength and scattering angle (θ = 0°: forward scattering). In addition, the approximate value of the complex refractive index at a wavelength of 0.5 μm is shown in the top left corner of each figure. See Sect. 4.1 for details.

In the text
thumbnail Fig. 2

Degree of signed polarization Ps (see Eq. (4)) as a function of wavelength and phase angle for various cloud compositions. As indicated by the spherical symbols at the top, a phase angle of 0° corresponds to the planet in full phase. See Sect. 4.2.1 for details.

In the text
thumbnail Fig. 3

Same as Fig. 2, but for different cloud compositions. See Sect. 4.2.2 for details.

In the text
thumbnail Fig. 4

Same as Fig. 2, but for different cloud compositions. See Sect. 4.2.3 for details.

In the text
thumbnail Fig. 5

Same as Fig. 2, but for different cloud compositions. See Sect. 4.2.4 for details.

In the text
thumbnail Fig. 6

Polarization contrast Cpol (see Eq. (10)) as a function of wavelength and phase angle for various cloud compositions and a cloud-free atmosphere (top left). As indicated by the spherical symbols at the top, a phase angle of 0° corresponds to the planet in full phase. See Sect. 5 for details.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.