Showing 1–35 of 35 results for author: Stam, D M

A general polarimetric model for transiting and non-transiting ringed exoplanets

Authors: Allard K. Veenstra, Jorge I. Zuluaga, Jaime A. Alvarado-Montes, Mario Sucerquia, Daphne M. Stam

Abstract: We explore the potential of polarimetry as a tool for detecting and characterizing exorings. For that purpose, we have improved the publicly available photometric code Pryngles by adding the results of radiative transfer calculations that fully include polarization and scattering by irregularly shaped particles. With this improved code, we compute the total and polarized fluxes and the degree of p… ▽ More We explore the potential of polarimetry as a tool for detecting and characterizing exorings. For that purpose, we have improved the publicly available photometric code Pryngles by adding the results of radiative transfer calculations that fully include polarization and scattering by irregularly shaped particles. With this improved code, we compute the total and polarized fluxes and the degree of polarization of a ringed gas giant along its orbit. We vary key model parameters such as the orbit inclination, ring size and orientation, particle albedo and optical thickness, and demonstrate the versatility of our code by predicting the total and polarized fluxes of the "puffed-up" planet HIP41378f assuming this planet has an opaque dusty ring. We find that spatially unresolved dusty rings can significantly modify the flux and polarization signals of the light that is reflected. Rings are expected to have a low polarization signal and will generally decrease the degree of polarization as the ring casts a shadow on the planet and/or blocks part of the light the planet reflects. During ring-plane crossings, when the thin ring is illuminated edge-on, a ringed exoplanet's flux and degree of polarization are close to those of a ring-less planet and generally appear as sharp changes in the flux and polarization curves. Ringed planets in edge-on orbits tend to be difficult to distinguish from ring-less planets in reflected flux and degree of polarization. We show that if HIP41378f is surrounded by a ring, its reflected flux (compared to the star) will be of the order of $10^{-9}$, and the ring would decrease the degree of polarization in a detectable way. The improved version of the photometric code Pryngles that we present here shows that dusty rings may produce distinct polarimetric features in light curves across a wide range of orbital configurations, orientations and ring optical properties. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 19 pages, 15 figures

Chasing rainbows and ocean glints: Inner working angle constraints for the Habitable Worlds Observatory

Authors: Sophia R. Vaughan, Timothy D. Gebhard, Kimberly Bott, Sarah L. Casewell, Nicolas B. Cowan, David S. Doelman, Matthew Kenworthy, Johan Mazoyer, Maxwell A. Millar-Blanchaer, Victor J. H. Trees, Daphne M. Stam, Olivier Absil, Lisa Altinier, Pierre Baudoz, Ruslan Belikov, Alexis Bidot, Jayne L. Birkby, Markus J. Bonse, Bernhard Brandl, Alexis Carlotti, Elodie Choquet, Dirk van Dam, Niyati Desai, Kevin Fogarty, J. Fowler , et al. (19 additional authors not shown)

Abstract: NASA is engaged in planning for a Habitable Worlds Observatory (HabWorlds), a coronagraphic space mission to detect rocky planets in habitable zones and establish their habitability. Surface liquid water is central to the definition of planetary habitability. Photometric and polarimetric phase curves of starlight reflected by an exoplanet can reveal ocean glint, rainbows and other phenomena caused… ▽ More NASA is engaged in planning for a Habitable Worlds Observatory (HabWorlds), a coronagraphic space mission to detect rocky planets in habitable zones and establish their habitability. Surface liquid water is central to the definition of planetary habitability. Photometric and polarimetric phase curves of starlight reflected by an exoplanet can reveal ocean glint, rainbows and other phenomena caused by scattering by clouds or atmospheric gas. Direct imaging missions are optimised for planets near quadrature, but HabWorlds' coronagraph may obscure the phase angles where such optical features are strongest. The range of accessible phase angles for a given exoplanet will depend on the planet's orbital inclination and/or the coronagraph's inner working angle (IWA). We use a recently-created catalog relevant to HabWorlds of 164 stars to estimate the number of exo-Earths that could be searched for ocean glint, rainbows, and polarization effects due to Rayleigh scattering. We find that the polarimetric Rayleigh scattering peak is accessible in most of the exo-Earth planetary systems. The rainbow due to water clouds at phase angles of ${\sim}20-60^\circ$ would be accessible with HabWorlds for a planet with an Earth equivalent instellation in ${\sim}{46}$ systems, while the ocean glint signature at phase angles of ${\sim}130-170^\circ$ would be accessible in ${\sim}{16}$ systems, assuming an IWA${=}62$ mas ($3λ/D$). Improving the IWA${=}41$ mas ($2λ/D$) increases accessibility to rainbows and glints by factors of approximately 2 and 3, respectively. By observing these scattering features, HabWorlds could detect a surface ocean and water cycle, key indicators of habitability. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: MNRAS accepted, 9 pages, 8 figures, 3 tables

Modelling reflected polarised light from close-in giant exoplanet WASP-96b using PolHEx (Polarisation of Hot Exoplanets)

Authors: Katy L. Chubb, Daphne M. Stam, Christiane Helling, Dominic Samra, Ludmila Carone

Abstract: We present the Polarisation of Hot Exoplanets (PolHEx) code for modelling the total flux (F) and degree of linear polarisation (P) of light spectra reflected by close-in, tidally locked exoplanets. We use the output from a global climate model (GCM) combined with a kinetic cloud model of hot Jupiter WASP-96b as a base to investigate effects of atmospheric longitudinal-latitudinal inhomogeneities o… ▽ More We present the Polarisation of Hot Exoplanets (PolHEx) code for modelling the total flux (F) and degree of linear polarisation (P) of light spectra reflected by close-in, tidally locked exoplanets. We use the output from a global climate model (GCM) combined with a kinetic cloud model of hot Jupiter WASP-96b as a base to investigate effects of atmospheric longitudinal-latitudinal inhomogeneities on these spectra. We model F and P-spectra as functions of wavelength and planet orbital phase for various model atmospheres. We find different materials and sizes of cloud particles to impact the reflected flux F, and particularly the linear polarisation state P. A range of materials are used to form inhomogeneous mixed-material cloud particles (Al2O3, Fe2O3, Fe2SiO4, FeO, Fe, Mg2SiO4, MgO, MgSiO3, SiO2, SiO, TiO2), with Fe2O3, Fe, and FeO the most strongly absorbing species. The cloud particles near the relatively cool morning terminator are expected to have smaller average sizes and a narrower size distribution than those near the warmer evening terminator, which leads to different reflected spectra at the respective orbital phases .We also find differences in the spectra of F and P as functions of orbital phase for irregularly or spherically shaped cloud particles. This work highlights the importance of including polarisation in models and future observations of the reflection spectra of exoplanets. △ Less

Submitted 2 November, 2023; v1 submitted 13 July, 2023; originally announced July 2023.

Comments: Accepted for publication in MNRAS

From exo-Earths to exo-Venuses -- Flux and Polarization Signatures of Reflected Light

Authors: Gourav Mahapatra, Fouad Abiad, Loic Rossi, Daphne M. Stam

Abstract: Terrestrial exoplanets in habitable zones are ubiquitous. It is, however, unknown which have Earth-like or Venus-like climates. Distinguishing different planet-types is crucial for determining whether a planet could be habitable. We investigate the potential of polarimetry for distinguishing exo-Earths from exo-Venuses. We present computed fluxes and polarisation of starlight that is reflected by… ▽ More Terrestrial exoplanets in habitable zones are ubiquitous. It is, however, unknown which have Earth-like or Venus-like climates. Distinguishing different planet-types is crucial for determining whether a planet could be habitable. We investigate the potential of polarimetry for distinguishing exo-Earths from exo-Venuses. We present computed fluxes and polarisation of starlight that is reflected by exoplanets with atmospheres in evolutionary states from current Earth to current Venus, with cloud compositions ranging from pure water to 0.75 sulfuric acid solution, for wavelengths between 0.3 and 2.5 microns. The polarisation of the reflected light shows larger variations with the planetary phase angle than the total flux. Across the visible, the largest polarisation is reached for an Earth-like atmosphere with water clouds, due to Rayleigh scattering above the clouds and the rainbow near 40 deg phase angle. In the near-infrared, the planet with a Venus-like CO2 atmosphere and thin water clouds shows the most prominent polarisation features due to scattering by the small cloud droplets. A planet around Alpha Centauri A would leave temporal variations on the order of 10E-13 W/m3 in the reflected flux and 10E-11 in the degree of polarisation along the planet's orbit for a spatially unresolved star-planet system. Star-planet contrasts are on the order of 10E-10. Current polarimeters cannot distinguish between the possible evolutionary phases of spatially unresolved terrestrial exoplanets, as a sensitivity near 10E-10 is required to discern the planet signal on the background of unpolarised starlight. Telescopes capable of reaching planet-star contrasts lower than 10E-9 should be able to observe the variation of the planet's resolved degree of polarisation as a function of its phase angle and thus to discern an exo-Earth from an exo-Venus based on its clouds' unique polarisation signatures. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 671, A165 (2023)

Polarized Signatures of a Habitable World: Comparing Models of an Exoplanet Earth with Visible and Near-infrared Earthshine Spectra

Authors: Kenneth E. Gordon, Theodora Karalidi, Kimberly M. Bott, Paulo A. Miles-Páez, Willeke Mulder, Daphne M. Stam

Abstract: In the JWST, Extremely Large Telescopes, and LUVOIR era, we expect to characterize a number of potentially habitable Earth-like exoplanets. However, the characterization of these worlds depends crucially on the accuracy of theoretical models. Validating these models against observations of planets with known properties will be key for the future characterization of terrestrial exoplanets. Due to i… ▽ More In the JWST, Extremely Large Telescopes, and LUVOIR era, we expect to characterize a number of potentially habitable Earth-like exoplanets. However, the characterization of these worlds depends crucially on the accuracy of theoretical models. Validating these models against observations of planets with known properties will be key for the future characterization of terrestrial exoplanets. Due to its sensitivity to the micro- and macro-physical properties of an atmosphere, polarimetry will be an important tool that, in tandem with traditional flux-only observations, will enhance the capabilities of characterizing Earth-like planets. In this paper we benchmark two different polarization-enabled radiative-transfer codes against each other and against unique linear spectropolarimetric observations of the earthshine that cover wavelengths from $\sim$0.4 to $\sim$2.3 $μ$m. We find that while the results from the two codes generally agree with each other, there is a phase dependency between the compared models. Additionally, with our current assumptions, the models from both codes underestimate the level of polarization of the earthshine. We also report an interesting discrepancy between our models and the observed 1.27 $μ$m $O_2$ feature in the earthshine, and provide an analysis of potential methods for matching this feature. Our results suggest that only having access to the 1.27 $μ$m $O_2$ feature coupled with a lack of observations of the $O_2$ A and B bands could result in a mischaracterization of an Earth-like atmosphere. Providing these assessments is vital to aid the community in the search for life beyond the solar system. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 24 pages, 18 figures

Spectropolarimetry of life: airborne measurements from a hot air balloon

Authors: Willeke Mulder, C. H. Lucas Patty, Stefano Spadaccia, Antoine Pommerol, Brice-Olivier Demory, Christoph U. Keller, Jonas G. Kühn, Frans Snik, Daphne M. Stam

Abstract: Does life exist outside our Solar System? A first step towards searching for life outside our Solar System is detecting life on Earth by using remote sensing applications. One powerful and unambiguous biosignature is the circular polarization resulting from the homochirality of biotic molecules and systems. We aim to investigate the possibility of identifying and characterizing life on Earth by us… ▽ More Does life exist outside our Solar System? A first step towards searching for life outside our Solar System is detecting life on Earth by using remote sensing applications. One powerful and unambiguous biosignature is the circular polarization resulting from the homochirality of biotic molecules and systems. We aim to investigate the possibility of identifying and characterizing life on Earth by using airborne spectropolarimetric observations from a hot air balloon during our field campaign in Switzerland, May 2022. In this work we present the optical-setup and the data obtained from aerial circular spectropolarimetric measurements of farmland, forests, lakes and urban sites. We make use of the well-calibrated FlyPol instrument that measures the fractionally induced circular polarization ($V/I$) of (reflected) light with a sensitivity of $<10^{-4}$. The instrument operates in the visible spectrum, ranging from 400 to 900 nm. We demonstrate the possibility to distinguish biotic from abiotic features using circular polarization spectra and additional broadband linear polarization information. We review the performance of our optical-setup and discuss potential improvements. This sets the requirements on how to perform future airborne spectropolarimetric measurements of the Earth's surface features from several elevations. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: 13 pages, 10 figures, to be submitted in SPIE Proceedings 12214-4

Ocean signatures in the total flux and polarization spectra of Earth-like exoplanets

Authors: Victor J. H. Trees, Daphne M. Stam

Abstract: Numerical simulations of starlight that is reflected by Earth-like exoplanets predict habitability signatures that can be searched for with future telescopes. We explore signatures of water oceans in the flux and polarization spectra of this reflected light. With an adding-doubling algorithm, we compute the total flux F, polarized flux Q and degree of polarization P of starlight reflected by dry a… ▽ More Numerical simulations of starlight that is reflected by Earth-like exoplanets predict habitability signatures that can be searched for with future telescopes. We explore signatures of water oceans in the flux and polarization spectra of this reflected light. With an adding-doubling algorithm, we compute the total flux F, polarized flux Q and degree of polarization P of starlight reflected by dry and ocean model planets with Earth-like atmospheres and patchy clouds. The oceans consist of Fresnel reflecting surfaces with wind-ruffled waves, foam and wave shadows, above natural blue seawater. Our results are presented as functions of wavelength (from 300 to 2500 nm with 1 nm resolution) and as functions of the planetary phase angle from 90 to 170 degrees. The ocean glint increases F, |Q| and P with increasing phase angle at non-absorbing wavelengths, and causes the spectra of F and |Q| for the various phase angles to intersect. In the near-infrared, Q is negative, i.e. the direction of polarization is perpendicular to the plane through the star, planet, and observer. In the P-spectra, the glint leaves dips (instead of peaks) in gaseous absorption bands. All those signatures are missing in the spectra of dry planets. The dips in P, and the negative Q in the near-infrared, can be searched for at a phase angle of 90 degrees, where the planet-star separation is largest. Those ocean signatures in polarized light do not suffer from false positive glint signals that could be due to clouds or reflecting dry surfaces. For heavily cloudy planets, ocean detection is possible when the glint is (partially) cloud-free. When modelling signals of planets with oceans, using horizontally inhomogeneous cloud covers is thus crucial. Observations spread over time would increase the probability of catching a cloud-free glint and detecting an ocean. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 664, A172 (2022)

Exoplanet Cartography using Convolutional Neural Networks

Authors: K. Meinke, D. M. Stam, P. M. Visser

Abstract: In the near-future, dedicated telescopes observe Earth-like exoplanets in reflected light, allowing their characterization. Because of the huge distances, every exoplanet will be a single pixel, but temporal variations in its spectral flux hold information about the planet's surface and atmosphere. We test convolutional neural networks for retrieving a planet's rotation axis, surface and cloud map… ▽ More In the near-future, dedicated telescopes observe Earth-like exoplanets in reflected light, allowing their characterization. Because of the huge distances, every exoplanet will be a single pixel, but temporal variations in its spectral flux hold information about the planet's surface and atmosphere. We test convolutional neural networks for retrieving a planet's rotation axis, surface and cloud map from simulated single-pixel flux and polarization observations. We investigate the assumption that the planets reflect Lambertian in the retrieval while their actual reflection is bidirectional, and of including polarization in retrievals. We simulate observations along a planet's orbit using a radiative transfer algorithm that includes polarization and bidirectional reflection by vegetation, desert, oceans, water clouds, and Rayleigh scattering in 6 spectral bands from 400 to 800 nm, at various photon noise levels. The surface-types and cloud patterns of the facets covering a model planet are based on probability distributions. Our networks are trained with simulated observations of millions of planets before retrieving maps of test planets. The neural networks can constrain rotation axes with a mean squared error (MSE) as small as 0.0097, depending on the orbital inclination. On a bidirectionally reflecting planet, 92% of ocean and 85% of vegetation, desert, and cloud facets are correctly retrieved, in the absence of noise. With realistic noise, it should still be possible to retrieve the main map features with a dedicated telescope. Except for face-on orbits, a network trained with Lambertian reflecting planets, yields significant retrieval errors when given observations of bidirectionally reflecting planets, in particular, brightness artefacts around a planet's pole. Including polarization improves retrieving the rotation axis and the accuracy of the retrieval of ocean and cloud facets. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 38 pages, 25 figures, 5 tables, accepted for publication in Astron. Astrophys

Journal ref: A&A 664, A59 (2022)

Biosignatures of the Earth I. Airborne spectropolarimetric detection of photosynthetic life

Authors: C. H. Lucas Patty, Jonas G. Kühn, Petar H. Lambrev, Stefano Spadaccia, H. Jens Hoeijmakers, Christoph Keller, Willeke Mulder, Vidhya Pallichadath, Olivier Poch, Frans Snik, Daphne M. Stam, Antoine Pommerol, Brice O. Demory

Abstract: Context. Homochirality is a generic and unique property of life on Earth and is considered a universal and agnostic biosignature. Homochirality induces fractional circular polarization in the incident light that it reflects. Because this circularly polarized light can be sensed remotely, it can be one of the most compelling candidate biosignatures in life detection missions. While there are also o… ▽ More Context. Homochirality is a generic and unique property of life on Earth and is considered a universal and agnostic biosignature. Homochirality induces fractional circular polarization in the incident light that it reflects. Because this circularly polarized light can be sensed remotely, it can be one of the most compelling candidate biosignatures in life detection missions. While there are also other sources of circular polarization, these result in spectrally flat signals with lower magnitude. Additionally, circular polarization can be a valuable tool in Earth remote sensing because the circular polarization signal directly relates to vegetation physiology. Aims. While high-quality circular polarization measurements can be obtained in the laboratory and under semi-static conditions in the field, there has been a significant gap to more realistic remote sensing conditions. Methods. In this study, we present sensitive circular spectropolarimetric measurements of various landscape elements taken from a fast-moving helicopter. Results. We demonstrate that during flight, within mere seconds of measurements, we can differentiate (S/N>5) between grass fields, forests, and abiotic urban areas. Importantly, we show that with only nonzero circular polarization as a discriminant, photosynthetic organisms can even be measured in lakes. Conclusions. Circular spectropolarimetry can be a powerful technique to detect life beyond Earth, and we emphasize the potential of utilizing circular spectropolarimetry as a remote sensing tool to characterize and monitor in detail the vegetation physiology and terrain features of Earth itself. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: 7 pages, 6 figures

Journal ref: A&A 651, A68 (2021)

Design of the Life Signature Detection Polarimeter LSDpol

Authors: Christoph U. Keller, Frans Snik, C. H. Lucas Patty, Dora Klindžić, Mariya Krasteva, David S. Doelman, Thomas Wijnen, Vidhya Pallichadath, Daphne M. Stam, Brice-Olivier Demory, Jonas G. Kühn, H. Jens Hoeijmakers, Antoine Pommerol, Olivier Poch

Abstract: Many biologically produced chiral molecules such as amino acids and sugars show a preference for left or right handedness (homochirality). Light reflected by biological materials such as algae and leaves therefore exhibits a small amount of circular polarization that strongly depends on wavelength. Our Life Signature Detection polarimeter (LSDpol) is optimized to measure these signatures of life.… ▽ More Many biologically produced chiral molecules such as amino acids and sugars show a preference for left or right handedness (homochirality). Light reflected by biological materials such as algae and leaves therefore exhibits a small amount of circular polarization that strongly depends on wavelength. Our Life Signature Detection polarimeter (LSDpol) is optimized to measure these signatures of life. LSDpol is a compact spectropolarimeter concept with no moving parts that instantaneously measures linear and circular polarization averaged over the field of view with a sensitivity of better than 1e-4. We expect to launch the instrument into orbit after validating its performance on the ground and from aircraft. LSDpol is based on a spatially varying quarter-wave retarder that is implemented with a patterned liquid-crystal. It is the first optical element to maximize the polarimetric sensitivity. Since this pattern as well as the entrance slit of the spectrograph have to be imaged onto the detector, the slit serves as the aperture, and an internal field stop limits the field of view. The retarder's fast axis angle varies linearly along one spatial dimension. A fixed quarter-wave retarder combined with a polarization grating act as the disperser and the polarizing beam-splitter. Circular and linear polarization are thereby encoded at incompatible modulation frequencies across the spectrum, which minimizes the potential cross-talk from linear into circular polarization. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 10 pages, 10 figures, SPIE Proceedings 11443-167

Journal ref: Proc. SPIE 11443, Space Telescopes and Instrumentation 2020: Optical, Infrared, and Millimeter Wave, 114433R (15 December 2020)

LOUPE: Observing Earth from the Moon to prepare for detecting life on Earth-like exoplanets

Authors: Dora Klindžić, Daphne M. Stam, Frans Snik, H. Jens Hoeijmakers, Michelle Willebrands, Theodora Karalidi, Vidhya Pallichadath, Chris N. van Dijk, Marco Esposito

Abstract: LOUPE, the Lunar Observatory for Unresolved Polarimetry of the Earth, is a small, robust spectro-polarimeter with a mission to observe the Earth as an exoplanet. Detecting Earth-like planets in stellar habitable zones is one of the key challenges of modern exoplanetary science. Characterising such planets and searching for traces of life requires the direct detection of their signals. LOUPE provid… ▽ More LOUPE, the Lunar Observatory for Unresolved Polarimetry of the Earth, is a small, robust spectro-polarimeter with a mission to observe the Earth as an exoplanet. Detecting Earth-like planets in stellar habitable zones is one of the key challenges of modern exoplanetary science. Characterising such planets and searching for traces of life requires the direct detection of their signals. LOUPE provides unique spectral flux and polarisation data of sunlight reflected by the Earth, the only planet known to harbor life. This data will be used to test numerical codes to predict signals of Earth-like exoplanets, to test algorithms that retrieve planet properties, and to fine-tune the design and observational strategies of future space observatories. From the Moon, LOUPE will continuously see the entire Earth, enabling it to monitor the signal changes due to the planet's daily rotation, weather patterns, and seasons, across all phase angles. Here, we present both the science case and the technology behind LOUPE's instrumental and mission design. △ Less

Submitted 3 August, 2020; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 13 pages, 5 figures. Accepted for publication in Royal Society Philosophical Transactions A. Corrected typos in v2

Colors of an Earth-like exoplanet -- Temporal flux and polarization signals of the Earth

Authors: A. Groot, L. Rossi, V. J. H. Trees, J. C. Y. Cheung, D. M. Stam

Abstract: Understanding the total flux and polarization signals of Earth-like planets and their spectral and temporal variability is essential for the future characterization of such exoplanets. We provide computed total (F) and linearly (Q and U) and circularly (V) polarized fluxes, and the degree of polarization P of sunlight that is reflected by a model Earth, to be used for instrument designs, optimizin… ▽ More Understanding the total flux and polarization signals of Earth-like planets and their spectral and temporal variability is essential for the future characterization of such exoplanets. We provide computed total (F) and linearly (Q and U) and circularly (V) polarized fluxes, and the degree of polarization P of sunlight that is reflected by a model Earth, to be used for instrument designs, optimizing observational strategies, and/or developing retrieval algorithms. We modeled a realistic Earth-like planet using one year of daily Earth-observation data: cloud parameters (distribution, optical thickness, top pressure, and particle effective radius), and surface parameters (distribution, surface type, and albedo). The Stokes vector of the disk-averaged reflected sunlight was computed for phase angles alpha from 0 to 180 degrees, and for wavelengths lambda from 350 to 865 nm. The total flux F is one order of magnitude higher than the polarized flux Q, and Q is two and four orders of magnitude higher than U and V, respectively. Without clouds, the peak-to-peak daily variations due to the planetary rotation increase with increasing lambda for F, Q, and P, while they decrease for U and V. Clouds modify but do not completely suppress the variations that are due to rotating surface features. With clouds, the variation in F increases with increasing lambda, while in Q, it decreases with increasing lambda, except at the largest phase angles. In earlier work, it was shown that with oceans, Q changes color from blue through white to red. The alpha where the color changes increases with increasing cloud coverage. Here, we show that this unique color change in Q also occurs when the oceans are partly replaced by continents, with or without clouds. The degree of polarization P shows a similar color change. Our computed fluxes and degree of polarization will be made publicly available. △ Less

Submitted 30 July, 2020; originally announced July 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

The polarimetric imaging mode of VLT/SPHERE/IRDIS II: Characterization and correction of instrumental polarization effects

Authors: R. G. van Holstein, J. H. Girard, J. de Boer, F. Snik, J. Milli, D. M. Stam, C. Ginski, D. Mouillet, Z. Wahhaj, H. M. Schmid, C. U. Keller, M. Langlois, K. Dohlen, A. Vigan, A. Pohl, M. Carbillet, D. Fantinel, D. Maurel, A. Origné, C. Petit, J. Ramos, F. Rigal, A. Sevin, A. Boccaletti, H. Le Coroller , et al. (9 additional authors not shown)

Abstract: Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and corr… ▽ More Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and correct the instrumental polarization effects of the complete optical system, i.e. the telescope and SPHERE/IRDIS. Methods. We create a detailed Mueller matrix model in the broadband filters Y-, J-, H- and Ks, and calibrate it using measurements with SPHERE's internal light source and observations of two unpolarized stars. We develop a data-reduction method that uses the model to correct for the instrumental polarization effects, and apply it to observations of the circumstellar disk of T Cha. Results. The instrumental polarization is almost exclusively produced by the telescope and SPHERE's first mirror and varies with telescope altitude angle. The crosstalk primarily originates from the image derotator (K-mirror). At some orientations, the derotator causes severe loss of signal (>90% loss in H- and Ks-band) and strongly offsets the angle of linear polarization. With our correction method we reach in all filters a total polarimetric accuracy of <0.1% in the degree of linear polarization and an accuracy of a few degrees in angle of linear polarization. Conclusions. The correction method enables us to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks, and is a vital tool for detecting unresolved (inner) disks and measuring the polarization of substellar companions. We have incorporated the correction method in a highly-automatic end-to-end data-reduction pipeline called IRDAP which is publicly available at https://irdap.readthedocs.io. △ Less

Submitted 28 September, 2019; originally announced September 2019.

Comments: 28 pages, 25 figures

Journal ref: A&A 633, A64 (2020)

The polarimetric imaging mode of VLT/SPHERE/IRDIS I: Description, data reduction and observing strategy

Authors: J. de Boer, M. Langlois, R. G. van Holstein, J. H. Girard, D. Mouillet, A. Vigan, K. Dohlen, F. Snik, C. U. Keller, C. Ginski, D. M. Stam, J. Milli, Z. Wahhaj, M. Kasper, H. M. Schmid, P. Rabou, L. Gluck, E. Hugot, D. Perret, P. Martinez, L. Weber, J. Pragt, J. -F. Sauvage, A. Boccaletti, H. Le Coroller , et al. (9 additional authors not shown)

Abstract: Context. Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of protoplanetary disks, and has the potential to be instrumental in characterizing exoplanets. VLT/SPHERE contains the InfraRed Dual-band Imager and Spectrograph (IRDIS) with a dual-beam polarimetric imaging (DPI) mode, which offers the capability to obtain linear polarization imag… ▽ More Context. Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of protoplanetary disks, and has the potential to be instrumental in characterizing exoplanets. VLT/SPHERE contains the InfraRed Dual-band Imager and Spectrograph (IRDIS) with a dual-beam polarimetric imaging (DPI) mode, which offers the capability to obtain linear polarization images at high contrast and resolution. Aims. We aim to provide an overview of IRDIS/DPI and study its optical design to improve observing strategies and data reduction. Methods. For H-band observations of TW Hya, we compare two data reduction methods that correct for instrumental polarization effects in different ways: a minimization of the noise image, and a polarimetric-model-based correction method that we present in Paper II of this study. Results. We use observations of TW Hya to illustrate the data reduction. In the images of the protoplanetary disk around this star we detect variability in the polarized intensity and angle of linear polarization with pointing-dependent instrument configuration. We explain these variations as instrumental polarization effects and correct for these effects using our model-based correction method. Conclusions. IRDIS/DPI has proven to be a very successful and productive high-contrast polarimetric imaging system. However, the instrument performance depends on the specific instrument configuration. We suggest adjustments to future observing strategies to optimize polarimetric efficiency in field tracking mode by avoiding unfavourable derotator angles. We recommend reducing on-sky data with the pipeline called IRDAP that includes the model-based correction method (described in Paper II) to optimally account for the remaining telescope and instrumental polarization effects and to retrieve the true polarization state of the incident light. △ Less

Submitted 28 September, 2019; originally announced September 2019.

Comments: 17 pages, 10 figures

Journal ref: A&A 633, A63 (2020)

A snapshot full-Stokes spectropolarimeter for detecting life on Earth

Authors: Frans Snik, Christoph U. Keller, David S. Doelman, Jonas Kühn, C. H. Lucas Patty, H. Jens Hoeijmakers, Vidhya Pallichadath, Daphne M. Stam, Antoine Pommerol, Olivier Poch, Brice-Olivier Demory

Abstract: We present the design of a point-and-shoot non-imaging full-Stokes spectropolarimeter dedicated to detecting life on Earth from an orbiting platform like the ISS. We specifically aim to map circular polarization in the spectral features of chlorophyll and other biopigments for our planet as a whole. These non-zero circular polarization signatures are caused by homochirality of the molecular and su… ▽ More We present the design of a point-and-shoot non-imaging full-Stokes spectropolarimeter dedicated to detecting life on Earth from an orbiting platform like the ISS. We specifically aim to map circular polarization in the spectral features of chlorophyll and other biopigments for our planet as a whole. These non-zero circular polarization signatures are caused by homochirality of the molecular and supramolecular configurations of organic matter, and are considered the most unambiguous biomarker. To achieve a fully solid-state snapshot design, we implement a novel spatial modulation that completely separates the circular and linear polarization channels. The polarization modulator consists of a patterned liquid-crystal quarter-wave plate inside the spectrograph slit, which also constitutes the first optical element of the instrument. This configuration eliminates cross-talk between linear and circular polarization, which is crucial because linear polarization signals are generally much stronger than the circular polarization signals. This leads to a quite unorthodox optical concept for the spectrograph, in which the object and the pupil are switched. We discuss the general design requirements and trade-offs of LSDpol (Life Signature Detection polarimeter), a prototype instrument that is currently under development. △ Less

Submitted 5 September, 2019; originally announced September 2019.

Journal ref: Proc. SPIE 11132-10 (2019)

Blue, white, and red ocean planets - Simulations of orbital variations in flux and polarization colors

Authors: Victor J. H. Trees, Daphne M. Stam

Abstract: An exoplanet's habitability will depend strongly on the presence of liquid water. Flux and/or polarization measurements of starlight that is reflected by exoplanets could help to identify exo-oceans. We investigate which broadband spectral features in flux and polarization phase functions of reflected starlight uniquely identify exo-oceans. We compute total fluxes F and polarized fluxes Q of starl… ▽ More An exoplanet's habitability will depend strongly on the presence of liquid water. Flux and/or polarization measurements of starlight that is reflected by exoplanets could help to identify exo-oceans. We investigate which broadband spectral features in flux and polarization phase functions of reflected starlight uniquely identify exo-oceans. We compute total fluxes F and polarized fluxes Q of starlight reflected by cloud-free and (partly) cloudy exoplanets, for wavelengths from 350 to 865 nm. The ocean surface has waves composed of Fresnel reflecting wave facets and whitecaps, and scattering within the water body is included. Total flux F, polarized flux Q, and degree of polarization P of ocean planets change color from blue, through white, to red at phase angles alpha ranging from 134-108 deg for F, and from 123-157 deg for Q, with cloud coverage fraction fc increasing from 0.0 to 1.0 for F, and to 0.98 for Q. The color change in P only occurs for fc ranging from 0.03-0.98, with the color crossing angle alpha ranging from 88-161 deg. The total flux F of a cloudy, zero surface albedo planet can also change color, and for fc=0.0, an ocean planet's F will not change color for surface pressures ps > 8 bars. Polarized flux Q of a zero surface albedo planet does not change color for any fc. The color change of P of starlight reflected by an exoplanet, from blue, through white, to red with increasing alpha above 88 deg, appears to identify a (partly) cloudy exo-ocean. The color change of polarized flux Q with increasing alpha above 123 deg appears to uniquely identify an exo-ocean, independent of surface pressure or cloud fraction. At the color changing phase angle, the angular distance between a star and its planet is much larger than at the phase angle where the glint appears in reflected light. The color change in polarization thus offers better prospects for detecting exo-oceans. △ Less

Submitted 26 April, 2019; v1 submitted 18 April, 2019; originally announced April 2019.

Comments: Accepted for publication in Astron. Astrophys; multicolumn version

Journal ref: A&A 626, A129 (2019)

Spectral and Temporal Variability of Earth Observed in Polarization

Authors: Michael F. Sterzik, Stefano Bagnulo, Daphne M. Stam, Claudia Emde, Mihail Manev

Abstract: We present a comprehensive set of spectropolarimetric observations of Earthshine as obtained by FORS2 at the VLT for phase angles from 50degree to 135degree (Sun-Earth-Moon angle), covering a spectral range from 430nm to 920nm. The degree of polarization in BVRI passbands, the differential polarization vegetation index, and the equivalent width of the O2A polarization band around 760nm are determi… ▽ More We present a comprehensive set of spectropolarimetric observations of Earthshine as obtained by FORS2 at the VLT for phase angles from 50degree to 135degree (Sun-Earth-Moon angle), covering a spectral range from 430nm to 920nm. The degree of polarization in BVRI passbands, the differential polarization vegetation index, and the equivalent width of the O2A polarization band around 760nm are determined with absolute errors around 0.1 percent in the degree of polarization. Earthshine polarization spectra are corrected for the effect of depolarization introduced by backscattering on the lunar surface, introducing systematic errors of the order of 1 percent in the degree of polarization. Distinct viewing sceneries such as observing the Atlantic or Pacific side in Earthshine yield statistically different phase curves. The equivalent width defined for the O2A band polarization is found to vary from -5nm to +2nm. A differential polarized vegetation index is introduced and reveals a larger vegetation signal for those viewing sceneries that contain larger fractions of vegetated surface areas. We corroborate the observed correlations with theoretical models from the literature, and conclude that the Vegetation Red Edge(VRE) is a robust and sensitive signature in polarization spectra of planet Earth. The overall behaviour of polarization of planet Earth in the continuum and in the O2A band can be explained by existing models. Biosignatures such as the O2A band and the VRE are detectable in Earthshine polarization with a high degree of significance and sensitivity. An in-depth understanding of Earthshines temporal and spectral variability requires improved models of Earths biosphere, as a prerequisite to interpret possible detections of polarised biosignatures in earthlike exoplanets in the future. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: 19 pages, 14 figures, 3 tables

Journal ref: A&A 622, A41 (2019)

Traces of exomoons in computed flux and polarization phase curves of starlight reflected by exoplanets

Authors: Javier Berzosa Molina, Loic Rossi, Daphne M. Stam

Abstract: Context: Detecting moons around exoplanets is a major goal of current and future observatories. Moons are suspected to influence rocky exoplanet habitability, and gaseous exoplanets in stellar habitable zones could harbour abundant and diverse moons to target in the search for extraterrestrial habitats. Exomoons will contribute to exoplanetary signals but are virtually undetectable with current me… ▽ More Context: Detecting moons around exoplanets is a major goal of current and future observatories. Moons are suspected to influence rocky exoplanet habitability, and gaseous exoplanets in stellar habitable zones could harbour abundant and diverse moons to target in the search for extraterrestrial habitats. Exomoons will contribute to exoplanetary signals but are virtually undetectable with current methods. Aims: We identify and analyse traces of exomoons in the temporal variation of total and polarised fluxes of starlight reflected by an Earth-like exoplanet and its spatially unresolved moon across all phase angles, with both orbits viewed in an edge-on geometry. Methods: We compute the total and linearly polarised fluxes, and the degree of linear polarization P of starlight that is reflected by the exoplanet with its moon along their orbits, accounting for the temporal variation of the visibility of the planetary and lunar disks, and including effects of mutual transits and mutual eclipses. Our computations pertain to a wavelength of 450 nm. Results: Total flux F shows regular dips due to planetary and lunar transits and eclipses. Polarization P shows regular peaks due to planetary transits and lunar eclipses, and P can increase and/or slightly decrease during lunar transits and planetary eclipses. Changes in F and P will depend on the radii of the planet and moon, on their reflective properties, and their orbits, and are about one magnitude smaller than the smooth background signals. The typical duration of a transit or an eclipse is a few hours. Conclusions: Traces of an exomoon due to planetary and lunar transits and eclipses show up in F and P of sunlight reflected by planet-moon systems and could be searched for in exoplanet flux and/or polarisation phase functions. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 618, A162 (2018)

Circular polarization signals of cloudy (exo)planets

Authors: Loïc Rossi, Daphne M. Stam

Abstract: The circular polarization of light that planets reflect is often neglected because it is very small compared to the linear polarization. It could, however, provide information on a planet's atmosphere and surface, and on the presence of life, because homochiral molecules that are the building blocks of life on Earth are known to reflect circularly polarized light. We compute $P_c$, the degree of… ▽ More The circular polarization of light that planets reflect is often neglected because it is very small compared to the linear polarization. It could, however, provide information on a planet's atmosphere and surface, and on the presence of life, because homochiral molecules that are the building blocks of life on Earth are known to reflect circularly polarized light. We compute $P_c$, the degree of circular polarization, for light that is reflected by rocky (exo)planets with liquid water or sulfuric acid solution clouds, both spatially resolved across the planetary disk and, for planets with patchy clouds, integrated across the planetary disk, for various planetary phase angles $α$. The optical thickness and vertical distribution of the atmospheric gas and clouds, the size parameter and refractive index of the cloud particles, and $α$ all influence $P_c$. Spatially resolved, $P_c$ varies between $\pm 0.20\%$ (the sign indicates the polarization direction). Only for small gas optical thicknesses above the clouds do significant sign changes (related to cloud particle properties) across the planets' hemispheres occur. For patchy clouds, the disk--integrated $P_c$ is typically smaller than $\pm 0.025\%$, with maximums for $α$ between $40^\circ$ and $70^\circ$, and $120^\circ$ to $140^\circ$. As expected, the disk--integrated $P_c$ is virtually zero at $α=0^\circ$ and 180$^\circ$. The disk--integrated $P_c$ is also very small at $α\approx 100^\circ$. Measuring circular polarization signals appears to be challenging with current technology. The small atmospheric circular polarization signal could, however, allow the detection of circular polarization due to homochiral molecules. Confirmation of the detectability of such signals requires better knowledge of the strength of circular polarization signals of biological sources. △ Less

Submitted 22 May, 2018; originally announced May 2018.

Comments: 15 pages, 11 figures, Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 616, A117 (2018)

PyMieDAP: a Python--Fortran tool to compute fluxes and polarization signals of (exo)planets

Authors: Loïc Rossi, Javier Berzosa-Molina, Daphne M. Stam

Abstract: PyMieDAP (the Python Mie Doubling-Adding Programme) is a Python--based tool for computing the total, linearly, and circularly polarized fluxes of incident unpolarized sun- or starlight that is reflected by, respectively, Solar System planets or moons, or exoplanets at a range of wavelengths. The radiative transfer computations are based on an adding--doubling Fortran algorithm and fully include po… ▽ More PyMieDAP (the Python Mie Doubling-Adding Programme) is a Python--based tool for computing the total, linearly, and circularly polarized fluxes of incident unpolarized sun- or starlight that is reflected by, respectively, Solar System planets or moons, or exoplanets at a range of wavelengths. The radiative transfer computations are based on an adding--doubling Fortran algorithm and fully include polarization for all orders of scattering. The model (exo)planets are described by a model atmosphere composed of a stack of homogeneous layers containing gas and/or aerosol and/or cloud particles bounded below by an isotropically, depolarizing surface (that is optionally black). The reflected light can be computed spatially--resolved and/or disk--integrated. Spatially--resolved signals are mostly representative for observations of Solar System planets (or moons), while disk--integrated signals are mostly representative for exoplanet observations. PyMieDAP is modular and flexible, and allows users to adapt and optimize the code according to their needs. PyMieDAP keeps options open for connections with external programs and for future additions and extensions. In this paper, we describe the radiative transfer algorithm that PyMieDAP is based on and the code's principal functionalities. And we provide benchmark results of PyMieDAP that can be used for testing its installation and for comparison with other codes. PyMieDAP is available online under the GNU GPL license at http://gitlab.com/loic.cg.rossi/pymiedap △ Less

Submitted 23 April, 2018; originally announced April 2018.

Comments: 15 pages, 7 figures, 4 tables. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 616, A147 (2018)

Combining angular differential imaging and accurate polarimetry with SPHERE/IRDIS to characterize young giant exoplanets

Authors: Rob G. van Holstein, Frans Snik, Julien H. Girard, Jozua de Boer, Christian Ginski, Christoph U. Keller, Daphne M. Stam, Jean-Luc Beuzit, David Mouillet, Markus Kasper, Maud Langlois, Alice Zurlo, Remco J. de Kok, Arthur Vigan

Abstract: Young giant exoplanets emit infrared radiation that can be linearly polarized up to several percent. This linear polarization can trace: 1) the presence of atmospheric cloud and haze layers, 2) spatial structure, e.g. cloud bands and rotational flattening, 3) the spin axis orientation and 4) particle sizes and cloud top pressure. We introduce a novel high-contrast imaging scheme that combines angu… ▽ More Young giant exoplanets emit infrared radiation that can be linearly polarized up to several percent. This linear polarization can trace: 1) the presence of atmospheric cloud and haze layers, 2) spatial structure, e.g. cloud bands and rotational flattening, 3) the spin axis orientation and 4) particle sizes and cloud top pressure. We introduce a novel high-contrast imaging scheme that combines angular differential imaging (ADI) and accurate near-infrared polarimetry to characterize self-luminous giant exoplanets. We implemented this technique at VLT/SPHERE-IRDIS and developed the corresponding observing strategies, the polarization calibration and the data-reduction approaches. By combining ADI and polarimetry we can characterize planets that can be directly imaged with a very high signal-to-noise ratio. We use the IRDIS pupil-tracking mode and combine ADI and principal component analysis to reduce speckle noise. We take advantage of IRDIS' dual-beam polarimetric mode to eliminate differential effects that severely limit the polarimetric sensitivity (flat-fielding errors, differential aberrations and seeing), and thus further suppress speckle noise. To correct for instrumental polarization effects, we apply a detailed Mueller matrix model that describes the telescope and instrument and that has an absolute polarimetric accuracy $\leq0.1\%$. Using this technique we have observed the planets of HR 8799 and the (sub-stellar) companion PZ Tel B. Unfortunately, we do not detect a polarization signal in a first analysis. We estimate preliminary $1σ$ upper limits on the degree of linear polarization of $\sim1\%$ and $\sim0.1\%$ for the planets of HR 8799 and PZ Tel B, respectively. The achieved sub-percent sensitivity and accuracy show that our technique has great promise for characterizing exoplanets through direct-imaging polarimetry. △ Less

Submitted 16 October, 2017; v1 submitted 21 September, 2017; originally announced September 2017.

Comments: 16 pages, 8 figures, 2 tables; v2: added acknowledgement and corrected two typos

Journal ref: Techniques and Instrumentation for Detection of Exoplanets VIII, Proc. of SPIE Vol. 10400, 1040015 (2017)

Polarized scattered light from self-luminous exoplanets. Three-dimensional scattering radiative transfer with ARTES

Authors: Tomas Stolker, Michiel Min, Daphne M. Stam, Paul Mollière, Carsten Dominik, Rens Waters

Abstract: Direct imaging has paved the way for atmospheric characterization of young and self-luminous gas giants. Scattering in a horizontally-inhomogeneous atmosphere causes the disk-integrated polarization of the thermal radiation to be linearly polarized, possibly detectable with the newest generation of high-contrast imaging instruments. We aim to investigate the effect of latitudinal and longitudinal… ▽ More Direct imaging has paved the way for atmospheric characterization of young and self-luminous gas giants. Scattering in a horizontally-inhomogeneous atmosphere causes the disk-integrated polarization of the thermal radiation to be linearly polarized, possibly detectable with the newest generation of high-contrast imaging instruments. We aim to investigate the effect of latitudinal and longitudinal cloud variations, circumplanetary disks, atmospheric oblateness, and cloud particle properties on the integrated degree and direction of polarization in the near-infrared. We have developed a three-dimensional Monte Carlo radiative transfer code (ARTES) for scattered light simulations in (exo)planetary atmospheres. The code is applicable to calculations of reflected light and thermal radiation in a spherical grid with a parameterized distribution of gas, clouds, hazes, and circumplanetary material. The disk-integrated degree of polarization of a horizontally-inhomogeneous atmosphere is maximal when the planet is flattened, the optical thickness of the equatorial clouds is large compared to the polar clouds, and the clouds are located at high altitude. For a flattened planet, the integrated polarization can both increase or decrease with respect to a spherical planet which depends on the horizontal distribution and optical thickness of the clouds. The direction of polarization can be either parallel or perpendicular to the projected direction of the rotation axis when clouds are zonally distributed. Rayleigh scattering by submicron-sized cloud particles will maximize the polarimetric signal whereas the integrated degree of polarization is significantly reduced with micron-sized cloud particles as a result of forward scattering. The presence of a cold or hot circumplanetary disk may also produce a detectable degree of polarization ($\lesssim$1%) even with a uniform cloud layer in the atmosphere. △ Less

Submitted 8 November, 2018; v1 submitted 28 June, 2017; originally announced June 2017.

Comments: 18 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 607, A42 (2017)

The O2 A-band in fluxes and polarization of starlight reflected by Earth-like exoplanets

Authors: Thomas Fauchez, Loic Rossi, Daphne M. Stam

Abstract: Earth-like, potentially habitable exoplanets are prime targets in the search for extraterrestrial life. Information about their atmosphere and surface can be derived by analyzing light of the parent star reflected by the planet. We investigate the influence of the surface albedo $A_{\rm s}$, the optical thickness $b_{\rm cloud}$ and altitude of water clouds, and the mixing ratio $η$ of biosignatur… ▽ More Earth-like, potentially habitable exoplanets are prime targets in the search for extraterrestrial life. Information about their atmosphere and surface can be derived by analyzing light of the parent star reflected by the planet. We investigate the influence of the surface albedo $A_{\rm s}$, the optical thickness $b_{\rm cloud}$ and altitude of water clouds, and the mixing ratio $η$ of biosignature O$_2$ on the strength of the O$_2$ A-band (around 760 nm) in flux and polarization spectra of starlight reflected by Earth-like exoplanets. Our computations for horizontally homogeneous planets show that small mixing ratios ($η$ < 0.4) will yield moderately deep bands in flux and moderate to small band strengths in polarization, and that clouds will usually decrease the band depth in flux and the band strength in polarization. However, cloud influence will be strongly dependent on their properties such as optical thickness, top altitude, particle phase, coverage fraction, horizontal distribution. Depending on the surface albedo, and cloud properties, different O$_2$ mixing ratios $η$ can give similar absorption band depths in flux and band strengths in polarization, in particular if the clouds have moderate to high optical thicknesses. Measuring both the flux and the polarization is essential to reduce the degeneracies, although it will not solve them, in particular not for horizontally inhomogeneous planets. Observations at a wide range of phase angles and with a high temporal resolution could help to derive cloud properties and, once those are known, the mixing ratio of O$_2$ or any other absorbing gas. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: 21 pages, 20 figures, accepted for publication in ApJ

Flux and polarization signals of spatially inhomogeneous gaseous exoplanets

Authors: T. Karalidi, D. M. Stam, D. Guirado

Abstract: We present numerically calculated, disk--integrated, spectropolarimetric signals of starlight that is reflected by vertically and horizontally inhomogeneous gaseous exoplanets. We include various spatial features that are present on Solar System's gaseous planets: belts and zones, cyclonic spots, and polar hazes, to test whether such features leave traces in the disk--integrated flux and polarizat… ▽ More We present numerically calculated, disk--integrated, spectropolarimetric signals of starlight that is reflected by vertically and horizontally inhomogeneous gaseous exoplanets. We include various spatial features that are present on Solar System's gaseous planets: belts and zones, cyclonic spots, and polar hazes, to test whether such features leave traces in the disk--integrated flux and polarization signals. Broadband flux and polarization signals of starlight that is reflected by gaseous exoplanets are calculated using an efficient, adding--doubling radiative transfer code, that fully includes single and multiple scattering and polarization. The planetary model atmospheres are vertically inhomogeneous and can be horizontally inhomogeneous, and contain gas molecules and/or cloud and/or aerosol particles. △ Less

Submitted 3 July, 2013; originally announced July 2013.

Comments: 12 pages, 11 figures

Looking for the rainbow on exoplanets covered by liquid and icy water clouds

Authors: T. Karalidi, D. M. Stam, J. W. Hovenier

Abstract: Looking for the primary rainbow in starlight that is reflected by exoplanets appears to be a promising method to search for liquid water clouds in exoplanetary atmospheres. Ice water clouds, that consist of water crystals instead of water droplets, could potentially mask the rainbow feature in the planetary signal by covering liquid water clouds. Here, we investigate the strength of the rainbow fe… ▽ More Looking for the primary rainbow in starlight that is reflected by exoplanets appears to be a promising method to search for liquid water clouds in exoplanetary atmospheres. Ice water clouds, that consist of water crystals instead of water droplets, could potentially mask the rainbow feature in the planetary signal by covering liquid water clouds. Here, we investigate the strength of the rainbow feature for exoplanets that have liquid and icy water clouds in their atmosphere, and calculate the rainbow feature for a realistic cloud coverage of Earth. We calculate flux and polarization signals of starlight that is reflected by horizontally and vertically inhomogeneous Earth--like exoplanets, covered by patchy clouds consisting of liquid water droplets or water ice crystals. The planetary surfaces are black. On a planet with a significant coverage of liquid water clouds only, the total flux signal shows a weak rainbow feature. Any coverage of the liquid water clouds by ice clouds, however, dampens the rainbow feature in the total flux, and thus the discovery of liquid water in the atmosphere. On the other hand, detecting the primary rainbow in the polarization signal of exoplanets appears to be a powerful tool for detecting liquid water in exoplanetary atmospheres, even when these clouds are partially covered by ice clouds. In particular, liquid water clouds covering as little as 10%-20% of the planetary surface, with more than half of these covered by ice clouds, still create a polarized rainbow feature in the planetary signal. Indeed, calculations of flux and polarization signals of an exoplanet with a realistic Earth--like cloud coverage, show a strong polarized rainbow feature. △ Less

Submitted 6 November, 2012; originally announced November 2012.

Comments: accepted for publication in Astronomy & Astrophysics

Modeled flux and polarisation signals of horizontally inhomogeneous exoplanets, applied to Earth--like planets

Authors: T. Karalidi, D. M. Stam

Abstract: We present modeled flux and linear polarisation signals of starlight that is reflected by spatially unresolved, horizontally inhomogeneous planets and discuss the effects of including horizontal inhomogeneities on the flux and polarisation signals of Earth-like exoplanets. Our code is based on an efficient adding--doubling algorithm, which fully includes multiple scattering by gases and aerosol/cl… ▽ More We present modeled flux and linear polarisation signals of starlight that is reflected by spatially unresolved, horizontally inhomogeneous planets and discuss the effects of including horizontal inhomogeneities on the flux and polarisation signals of Earth-like exoplanets. Our code is based on an efficient adding--doubling algorithm, which fully includes multiple scattering by gases and aerosol/cloud particles. We divide a model planet into pixels that are small enough for the local properties of the atmosphere and surface (if present) to be horizontally homogeneous. Given a planetary phase angle, we sum up the reflected total and linearly polarised fluxes across the illuminated and visible part of the planetary disk, taking care to properly rotate the polarized flux vectors towards the same reference plane. We compared flux and polarisation signals of simple horizontally inhomogeneous model planets against results of the weighted sum approximation, in which signals of horizontally homogeneous planets are combined. Apart from cases in which the planet has only a minor inhomogeneity, the signals differ significantly. In particular, the shape of the polarisation phase function appears to be sensitive to the horizontal inhomogeneities. The same holds true for Earth-like model planets with patchy clouds above an ocean and a sandy continent. Our simulations clearly show that horizontal inhomogeneities leave different traces in flux and polarisation signals. Combining flux with polarisation measurements would help retrieving the atmospheric and surface patterns on a planet. △ Less

Submitted 11 October, 2012; originally announced October 2012.

Journal ref: published in Astronomy & Astrophysics, volume 546, A56, October, 2012

The influence of forward-scattered light in transmission measurements of (exo)planetary atmospheres

Authors: R. J. de Kok, D. M. Stam

Abstract: [Abridged] The transmission of light through a planetary atmosphere can be studied as a function of altitude and wavelength using stellar or solar occultations, giving often unique constraints on the atmospheric composition. For exoplanets, a transit yields a limb-integrated, wavelength-dependent transmission spectrum of an atmosphere. When scattering haze and/or cloud particles are present in the… ▽ More [Abridged] The transmission of light through a planetary atmosphere can be studied as a function of altitude and wavelength using stellar or solar occultations, giving often unique constraints on the atmospheric composition. For exoplanets, a transit yields a limb-integrated, wavelength-dependent transmission spectrum of an atmosphere. When scattering haze and/or cloud particles are present in the planetary atmosphere, the amount of transmitted flux not only depends on the total optical thickness of the slant light path that is probed, but also on the amount of forward-scattering by the scattering particles. Here, we present results of calculations with a three-dimensional Monte Carlo code that simulates the transmitted flux during occultations or transits. For isotropically scattering particles, like gas molecules, the transmitted flux appears to be well-described by the total atmospheric optical thickness. Strongly forward-scattering particles, however, such as commonly found in atmospheres of Solar System planets, can increase the transmitted flux significantly. For exoplanets, such added flux can decrease the apparent radius of the planet by several scale heights, which is comparable to predicted and measured features in exoplanet transit spectra. We performed detailed calculations for Titan's atmosphere between 2.0 and 2.8 micron and show that haze and gas abundances will be underestimated by about 8% if forward-scattering is ignored in the retrievals. At shorter wavelengths, errors in the gas and haze abundances and in the spectral slope of the haze particles can be several tens of percent, also for other Solar System planetary atmospheres. We also find that the contribution of forward-scattering can be fairly well described by modelling the atmosphere as a plane-parallel slab. △ Less

Submitted 14 August, 2012; originally announced August 2012.

Comments: Icarus, accepted for publication

Atmospheric characterization of cold exoplanets using a 1.5-m coronagraphic space telescope

Authors: A. -L. Maire, R. Galicher, A. Boccaletti, P. Baudoz, J. Schneider, K. L. Cahoy, D. M. Stam, W. A. Traub

Abstract: Context. High-contrast imaging is currently the only available technique for the study of the thermodynamical and compositional properties of exoplanets in long-period orbits. The SPICES project is a coronagraphic space telescope dedicated to the spectro-polarimetric analysis of gaseous and icy giant planets as well as super-Earths at visible wavelengths. So far, studies for high-contrast imaging… ▽ More Context. High-contrast imaging is currently the only available technique for the study of the thermodynamical and compositional properties of exoplanets in long-period orbits. The SPICES project is a coronagraphic space telescope dedicated to the spectro-polarimetric analysis of gaseous and icy giant planets as well as super-Earths at visible wavelengths. So far, studies for high-contrast imaging instruments have mainly focused on technical feasibility because of the challenging planet/star flux ratio of 10-8-10-10 required at short separations (200 mas or so) to image cold exoplanets. However, the analysis of planet atmospheric/surface properties has remained largely unexplored. Aims. The aim of this paper is to determine which planetary properties SPICES or an equivalent direct imaging mission can measure, considering realistic reflected planet spectra and instrument limitation. Methods. We use numerical simulations of the SPICES instrument concept and theoretical planet spectra to carry out this performance study. Results. We find that the characterization of the main planetary properties (identification of molecules, effect of metallicity, presence of clouds and type of surfaces) would require a median signal-to-noise ratio of at least 30. In the case of a solar-type star \leq 10 pc, SPICES will be able to study Jupiters and Neptunes up to ~5 and ~2 AU respectively. It would also analyze cloud and surface coverage of super-Earths of radius 2.5 RE at 1 AU. Finally, we determine the potential targets in terms of planet separation, radius and distance for several stellar types. For a Sun analog, we show that SPICES could characterize Jupiters (M \geq 30 ME) as small as 0.5 Jupiter radii at ~2 AU up to 10 pc, and super-Earths at 1-2 AU for the handful of stars that exist within 4-5 pc. Potentially, SPICES could perform analysis of a hypothetical Earth-size planet around alpha Cen A and B. △ Less

Submitted 13 March, 2012; originally announced March 2012.

Comments: 18 pages, 20 figures, 9 tables, accepted for publication in A&A

Observing the Earth as an exoplanet with LOUPE, the Lunar Observatory for Unresolved Polarimetry of Earth

Authors: T. Karalidi, D. M. Stam, F. Snik, S. Bagnulo, W. B. Sparks, C. U. Keller

Abstract: The detections of small, rocky exoplanets have surged in recent years and will likely continue to do so. To know whether a rocky exoplanet is habitable, we have to characterise its atmosphere and surface. A promising characterisation method for rocky exoplanets is direct detection using spectropolarimetry. This method will be based on single pixel signals, because spatially resolving exoplanets is… ▽ More The detections of small, rocky exoplanets have surged in recent years and will likely continue to do so. To know whether a rocky exoplanet is habitable, we have to characterise its atmosphere and surface. A promising characterisation method for rocky exoplanets is direct detection using spectropolarimetry. This method will be based on single pixel signals, because spatially resolving exoplanets is impossible with current and near-future instruments. Well-tested retrieval algorithms are essential to interpret these single pixel signals in terms of atmospheric composition, cloud and surface coverage. Observations of Earth itself provide the obvious benchmark data for testing such algorithms. The observations should provide signals that are integrated over the Earth's disk, that capture day and night variations, and all phase angles. The Moon is a unique platform from where the Earth can be observed as an exoplanet, undisturbed, all of the time. Here, we present LOUPE, the Lunar Observatory for Unresolved Polarimetry of Earth, a small and robust spectropolarimeter to observe our Earth as an exoplanet. △ Less

Submitted 1 March, 2012; originally announced March 2012.

Comments: 14 pages, 3 figures, submitted in special Issue of Planetary and Space Science on Scientific Preparations for Lunar Exploration

Characterizing exoplanetary atmospheres through infrared polarimetry

Authors: R. J. de Kok, D. M. Stam, T. Karalidi

Abstract: Planets can emit polarized thermal radiation, just like brown dwarfs. We present calculated thermal polarization signals from hot exoplanets, using an advanced radiative transfer code that fully includes all orders of scattering by gaseous molecules and cloud particles. The code spatially resolves the disk of the planet, allowing simulations for horizontally inhomogeneous planets. Our results show… ▽ More Planets can emit polarized thermal radiation, just like brown dwarfs. We present calculated thermal polarization signals from hot exoplanets, using an advanced radiative transfer code that fully includes all orders of scattering by gaseous molecules and cloud particles. The code spatially resolves the disk of the planet, allowing simulations for horizontally inhomogeneous planets. Our results show that the degree of linear polarization, P, of an exoplanet's thermal radiation is expected to be highest near the planet's limb and that this P depends on the temperature and its gradient, the scattering properties and the distribution of the cloud particles. Integrated over the disk of a spherically symmetric planet, P of the thermal radiation equals zero. However, for planets that appear spherically asymmetric, e.g. due to flattening, cloud bands or spots in their atmosphere, differences in their day and night sides, and/or obscuring rings, P is often larger than 0.1 %, in favorable cases even reaching several percent at near-infrared wavelengths. Detection of thermal polarization signals can give access to planetary parameters that are otherwise hard to obtain: it immediately confirms the presence of clouds, and P can then constrain atmospheric inhomogeneities and the flattening due to the planet's rotation rate. For zonally symmetric planets, the angle of polarization will yield the components of the planet's spin axis normal to the line-of-sight. Finally, our simulations show that P is generally more sensitive to variability in a cloudy planet's atmosphere than the thermal flux is, and could hence better reveal certain dynamical processes. △ Less

Submitted 5 August, 2011; originally announced August 2011.

Comments: 9 pages, 5 figures, accepted for publication in ApJ

The influence of non-isotropic scattering of thermal radiation on spectra of brown dwarfs and hot exoplanets

Authors: R. J. de Kok, Ch. Helling, D. M. Stam, P. Woitke, S. Witte

Abstract: (abridged) We calculate near-infrared thermal emission spectra using a doubling-adding radiative transfer code, which includes scattering by clouds and haze. Initial temperature profiles and cloud optical depths are taken from the drift-phoenix brown dwarf model. As is well known, cloud particles change the spectrum compared to when clouds are ignored. The clouds reduce fluxes in the near-infrared… ▽ More (abridged) We calculate near-infrared thermal emission spectra using a doubling-adding radiative transfer code, which includes scattering by clouds and haze. Initial temperature profiles and cloud optical depths are taken from the drift-phoenix brown dwarf model. As is well known, cloud particles change the spectrum compared to when clouds are ignored. The clouds reduce fluxes in the near-infrared spectrum and make it redder than for the clear sky case. We also confirm that not including scattering in the spectral calculations can result in errors on the spectra of many tens of percent, both in magnitude and in variations with wavelength. This is especially apparent for particles that are larger than the wavelength and only have little iron in them. Scattering particles will show deeper absorption features than absorbing (e.g. iron) particles and particle size will also affect the calculated infrared colours. Large particles also tend to be strongly forward-scattering, and we show that assuming isotropic scattering in this case also leads to very large errors in the spectrum. Thus, care must be taken in the choice of radiative transfer method for heat balance or spectral calculations when clouds are present in the atmosphere. Besides the choice of radiative transfer method, the type of particles that are predicted by models will change conclusions about e.g. infrared colours and trace gas abundances. As a result, knowledge of the scattering properties of the clouds is essential when deriving temperature profiles or gas abundances from direct infrared observations of exoplanets or brown dwarfs and from secondary eclipse measurements of transiting exoplanets, since scattering clouds will change the depth of gas absorption features, among other things. Thus, ignoring the presence of clouds can yield retrieved properties that differ significantly from the real atmospheric properties. △ Less

Submitted 3 June, 2011; v1 submitted 16 May, 2011; originally announced May 2011.

Comments: Accepted for publication in Astronomy and Astrophysics. The abstract and a part of the introduction have been re-worded compared to the accepted version to avoid misinterpretation of the paper as much as possible

Journal ref: Astronomy and Astrophysics, 531, A67, 2011

Behavior of the reflection function of a plane-parallel medium for directions of incidence and reflection tending to horizontal directions

Authors: J. W. Hovenier D. M. Stam

Abstract: The atmospheres of (exo) planets and moons, as well as reflection nebulae, contain in general independently scattering particles in random orientation and are often supposed to be plane-parallel. Relations are presented for the (bidirectional) reflection function and several related functions of such a medium in case the directions of incidence and reflection both tend to horizontal directions.… ▽ More The atmospheres of (exo) planets and moons, as well as reflection nebulae, contain in general independently scattering particles in random orientation and are often supposed to be plane-parallel. Relations are presented for the (bidirectional) reflection function and several related functions of such a medium in case the directions of incidence and reflection both tend to horizontal directions. The results are quite general. The medium may be semi-infinite or finite, with or without a reflecting surface underneath, and vertically homogeneous or inhomogeneous. Some approximative formulae for the reflection function of a plane-parallel medium with independently scattering particles in random orientation, including Lambert's law, may be very inaccurate if the directions of incidence and reflection are both nearly horizontal. △ Less

Submitted 22 September, 2008; originally announced September 2008.

Comments: Accepted for publication in Astron. Astrophys

Scattering matrices and expansion coefficients of Martian analogue palagonite particles

Authors: E. C. Laan, H. Volten, D. M. Stam, O. Munoz, J. W. Hovenier, T. L. Roush

Abstract: We present measurements of ratios of elements of the scattering matrix of Martian analogue palagonite particles for scattering angles ranging from 3 to 174 degrees and a wavelength of 632.8 nm. To facilitate the use of these measurements in radiative transfer calculations we have devised a method that enables us to obtain, from these measurements, a normalized synthetic scattering matrix coverin… ▽ More We present measurements of ratios of elements of the scattering matrix of Martian analogue palagonite particles for scattering angles ranging from 3 to 174 degrees and a wavelength of 632.8 nm. To facilitate the use of these measurements in radiative transfer calculations we have devised a method that enables us to obtain, from these measurements, a normalized synthetic scattering matrix covering the complete scattering angle range from 0 to 180 degrees. Our method is based on employing the coefficients of the expansions of scattering matrix elements into generalized spherical functions. The synthetic scattering matrix elements and/or the expansion coefficients obtained in this way, can be used to include multiple scattering by these irregularly shaped particles in (polarized) radiative transfer calculations, such as calculations of sunlight that is scattered in the dusty Martian atmosphere. △ Less

Submitted 16 September, 2008; originally announced September 2008.

Comments: 34 pages 7 figures 1 table

Using SPICA Space Telescope to characterize Exoplanets

Authors: J. R. Goicoechea, B. Swinyard, G. Tinetti, T. Nakagawa, K. Enya, M. Tamura, M. Ferlet, K. G. Isaak, M. Wyatt, A. D. Aylward, M. Barlow, J. P. Beaulieu, A. Boccaletti, J. Cernicharo, J. Cho, R. Claudi, H. Jones, H. Lammer, A. Leger, J. Martín-Pintado, S. Miller, F. Najarro, D. Pinfield, J. Schneider, F. Selsis , et al. (4 additional authors not shown)

Abstract: We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to… ▽ More We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to the next stage of ESA's Cosmic Vision 2015-2025 selection process. In this White Paper we present the main specifications of the three instruments currently baselined for SPICA: a mid-infrared (MIR) coronagraph (~3.5 to ~27 um) with photometric and spectral capabilities (R~200), a MIR wide-field camera and high resolution spectrometer (R~30,000), and a far-infrared (FIR ~30 to ~210 um) imaging spectrometer - SAFARI - led by a European consortium. We discuss their capabilities in the context of MIR direct observations of exo-planets (EPs) and multiband photometry/high resolution spectroscopy observations of transiting exo-planets. We conclude that SPICA will be able to characterize the atmospheres of transiting exo-planets down to the super-Earth size previously detected by ground- or space-based observatories. It will also directly detect and characterize Jupiter/Neptune-size planets orbiting at larger separation from their parent star (>5-10 AU), by performing quantitative atmospheric spectroscopy and studying proto-planetary and debris disks. In addition, SPICA will be a scientific and technological precursor for future, more ambitious, IR space missions for exo-planet direct detection as it will, for example, quantify the prevalence exo-zodiacal clouds in planetary systems and test coronographic techniques, cryogenic systems and lightweight, high quality telescopes. (abridged) △ Less

Submitted 15 September, 2008; v1 submitted 1 September, 2008; originally announced September 2008.

Comments: A White Paper for ESA's Exo-Planet Roadmap Advisory Team, submitted on 2008 July 29

Spectropolarimetric signatures of Earth-like extrasolar planets

Authors: D. M. Stam

Abstract: We present results of numerical simulations of the flux (irradiance), F, and the degree of polarization (i.e. the ratio of polarized to total flux), P, of light that is reflected by Earth-like extrasolar planets orbiting solar-type stars, as functions of the wavelength (from 0.3 to 1.0 micron, with 0.001 micron spectral resolution) and as functions of the planetary phase angle. We use different… ▽ More We present results of numerical simulations of the flux (irradiance), F, and the degree of polarization (i.e. the ratio of polarized to total flux), P, of light that is reflected by Earth-like extrasolar planets orbiting solar-type stars, as functions of the wavelength (from 0.3 to 1.0 micron, with 0.001 micron spectral resolution) and as functions of the planetary phase angle. We use different surface coverages for our model planets, including vegetation and a Fresnel reflecting ocean, and clear and cloudy atmospheres. Our adding-doubling radiative transfer algorithm, which fully includes multiple scattering and polarization, handles horizontally homogeneous planets only; we simulate fluxes and polarization of horizontally inhomogeneous planets by weighting results for homogeneous planets. Like the flux, F, the degree of polarization, P, of the reflected starlight is shown to depend strongly on the phase angle, on the composition and structure of the planetary atmosphere, on the reflective properties of the underlying surface, and on the wavelength, in particular in wavelength regions with gaseous absorption bands. The sensitivity of P to a planet's physical properties appears to be different than that of F. Combining flux with polarization observations thus makes for a strong tool for characterizing extrasolar planets. The calculated total and polarized fluxes will be made available through the CDS. △ Less

Submitted 26 July, 2007; originally announced July 2007.

Comments: 31 pages text, 17 figures, 1 table Submitted to A&A