Version of November 14, 2018 A Preprint typeset using LTEX style emulateapj v. 11/10/09 DERIVING TELESCOPE MUELLER MATRICES USING DAYTIME SKY POLARIZATION OBSERVATIONS David M. Harrington1, J. R. Kuhn2, and Shannon Hall3 1Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI, 96822 2Institute for Astronomy Maui, University of Hawaii, 34 Ohia Ku St., Pukalani, HI, 96768 and 3Department of Astronomy, Whitman College, 345 Boyer Ave., Walla Walla, WA, 99362 Version of November 14, 2018 ABSTRACT Telescopes often modify the input polarization of a source so that the measured circular or linear output state of the optical signal can be signficantly different from the input. This mixing, or polar- ization “cross-talk”, is defined by the optical system Mueller matrix. We describe here an efficient method for recovering the input polarization state of the light and the full 4 × 4 Mueller matrix of the telescope with an accuracy of a few percent without external masks or telescope hardware mod- ification. Observations of the bright, highly polarized daytime sky using the Haleakala 3.7m AEOS telescope and a coud´espectropolarimeter demonstrate the technique. Subject headings: Astronomical Techniques, Astronomical Instrumentation, Data Analysis and Tech- niques 1. INTRODUCTION illumination of the telescope optics during target obser- Spectropolarimetry is a powerful tool which may be vations. enhanced by efficient techniques for telescope-detector The new technique presented here is designed to yield polarization calibration. Unfortunately, it is common the full Mueller matrix calibration of many telescope- for modern altitude-azimuth telescopes with multiple polarimeter systems with an accuracy of a few percent. mirror reflections to scramble the input polarization It requires only two or more observations of the linearly of the light before it reaches the analyzer of a coud´e, polarized sky due to solar Rayleigh scattering. The day- Nasmyth, or Gregorian polarimeter. Solar astronomers time sky is normally easily observable, highly polarized, in particular have developed a repretoire of tech- and is a relatively well characterized calibration source. niques for calibrating telescopes (Elmore et al. 1992; It is observable without using precious dark observing Kuhn et al. 1993; Giro et al. 2003; Socas-Navarro time and requires no hardware calibration modifications 2005a,b; Socas-Navarro et al. 2005; Selbing 2005; of the telescope. As long as multiple linear polarization Socas-Navarroet al. 2006; Ichimoto et al. 2008; input states are detected at each telescope pointing, the Keller & Snik 2009; Elmore et al. 2010. Our ex- full telescope Mueller matrix can be recovered with an ac- perience with the Advanced Electro-Optical System curacy of a few percent. This is realized by observing the (AEOS) 3.67m telescope is that linearly polarized solar illuminated sky with the Sun at different zenith and light can even, in some circumstances, be transformed azimuth positions. This paper illustrates the technique to nearly 100% circularly polarized light. There are using the Haleakala Air Force AEOS 3.7m telescope on 7 highly oblique reflections before polarization is Maui. analyzed at our High Resolution Visible and Infra- red Spectropolarimeter (HiVIS) (cf. Hodapp et al. 1.1. Polarization 2000, Thornton et al. 2003, Harrington et al. The following discussion of polarization formalism 2006, 2009, 2010, Harrington & Kuhn 2007, 2008; closely follows Collet 1992 and Clarke 2009. In the Harringtion & Kuhn 2009a,b). Stokes formalism, the polarization state of light is de- arXiv:1105.2615v1 [astro-ph.IM] 13 May 2011 Effective calibration of the telescope and instrument noted as a 4-vector: minimizes these systematic errors and enhances the po- larimetric precision of an observation. Moving optical T Si = [I,Q,U,V ] (1) elements, detector effects such as non-uniform sensitiv- ity or nonlinearity, atmospheric seeing, and transparency In this formalism, I represents the total intensity, Q variations can all induce polarization errors that are of- and U the linearly polarized intensity along polarization ten larger than the photon noise. Strategies to minimize, position angles 0◦ and 45◦ in the plane perpendicular to stabilize and calibrate the instrument and telescope of- the light beam, and V is the right-handed circularly po- ten involve use of non-moving optics and emphasize rapid larized intensity. Note that according to this definition, chopping and common-path differential techniques. Cali- linear polarization along angles 90◦ and 135◦ will be de- bration techniques include observing unpolarized and po- noted as −Q and −U, respectively. The typical conven- larized standard stars, and use of polarized optical cal- tion for astronomical polarimetry is for the +Q electric ibration units that are designed to mimic the nominal field vibration direction to be aligned North-South while +U has the electric field vibration direction aligned to [email protected] North-East and South-West. [email protected] The normalized Stokes parameters are denoted with [email protected] lower case and are defined as: 2 Harrington et al. is also available using industry standard atmospheric ra- [1,q,u,v]T = [I,Q,U,V ]T /I (2) diative transfer software packages such as MODTRAN (cf. Fetrow et al. 2002, Berk et al. 2006). The degree of polarization can be defined as a ratio of polarized light to the total intensity of the beam: Q2 + U 2 + V 2 P = = q2 + u2 + v2 (3) p I For details on polarization ofp light and stellar spec- tropolarimetry, see Collet 1992 and Clarke 2009. To describe how polarized light propagates through any optical system, the Mueller matrix is constructed which specifies how the incident polarization state is transferred to the output polarization state. The Mueller matrix is a 4 × 4 set of transfer coefficients which when S multiplied by the input Stokes vector ( iinput ) gives the S output Stokes vector ( ioutput ): S M S ioutput = ij iinput (4) If the Mueller matrix for a system is known, then one inverts the matrix and deprojects a set of measurements to recover the inputs. One can represent the individ- ual Mueller matrix terms as describing how one incident state transfers to another. In this work we will use the notation: Fig. 1.— The celestial triangle representing the geometry for the sky polarization at any telescope pointing. γ is the angu- II QI UI VI lar distance between the telescope pointing and the sun. θs is the solar angular distance from the zenith. θ is the angular M IQ QQ UQ VQ ij = IU QU UU VU (5) distance between the telescope pointing and the zenith. φ is IV QV UV VV the angle between the zenith direction and the solar direction at the telescope pointing. ψ is the angle between the tele- 1.2. Deriving Telescope Mueller Matrix Elements scope pointing and the solar direction at the zenith. The law of Cosines is used to solve for any angles needed to compute Typical calibration schemes on smaller telescopes of- degree of polarization and position angle of polarization. ten use fixed polarizing filters placed over the telescope aperture to provide known input states that are de- tected and yield terms of the Mueller matrix. This ap- There are many atmospheric and geometric con- proach has been used for solar telescopes (Socas-Navarro siderations that determine the skylight polarization 2005a,b, Socas-Navarro et al. 2005, 2006). Another ap- at a particular observatory site. The linear po- proach, also used in the solar observations has been larization amplitude angle can depend on the so- to image known sources and use spectropolarimet- lar elevation, atmospheric aerosol content, aerosol ric data and Zeeman effect symmetry properties to vertical distribution, aerosol scattering phase func- isolate and determine terms in the Mueller matrix tion, wavelength of the observation and secondary (Kuhn et al. 1993, Elmore et al. 2010). Night-time sources of illumination (cf. Horv´ath et al. 2002a,b, observations can use unpolarized and polarized stan- Lee 1998, Liu & Voss 1997, Suhai & Horv´ath 2004, dard stars to measure polarization properties of tele- G´al et al. 2001, Vermeulen et al. 2000, Pomozi et al. scopes (cf. Hsu & Breger 1982, Schmidt et al. 1992, 2001, Cronin et al. 2005, 2006, Heged¨us et al. 2007). Gil-Hutton & Benavidez 2003, Fossati et al. 2007). Anisotropic scattered sunlight from reflections off land Many studies have either measured and calibrated tele- or water can be highly polarized and temporally scopes, measured mirror properties or attempted to variable (Litvinov et al. 2010, Peltoniemi et al. 2009, design instruments with minimal polarimetric defects He et al. 2010, Salinas & Liew 2007, Ota et al. 2010, (cf. S´anchez Almeida et al. 1991, Giro et al. 2003, Kisselev & Bulgarelli 2004). Aerosol particle op- Patat & Romaniello 2006,Tinbergen 2007, Joos et al. tical properties and vertical distributions also vary 2008, van Harten et al. 2009, Roelfsema et al. 2010). (cf. Wu & Jin 1997, Shukurov & Shukurov 2006, Vermeulen et al. 2000, Ougolnikov & Maslov 2002, 1.3. The Polarized Sky as a Calibration 2005a,b, 2009a,b, Ugolnikov et al. 2004). The polar- The observed polarization of the daytime sky is a ization can change across atmospheric absorption bands useful calibration source. Scattered sunlight is bright, or can be influenced by other scattering mechanisms highly (linearly) polarized and typically illuminates (cf. Boesche et al. 2006, Zeng et al. 2008 Aben et al. the telescope optics more realistically than calibration 1999, 2001). screens. A single-scattering atomic Rayleigh calcula- Deviations from a single scattering Rayleigh model tion (cf. Coulson 1980, 1998) is often adequate to de- grow as the aerosol, cloud, ground or sea-surface scat- scribe the skylight polarization. More realistic model- tering sources affect the telescope line-of-sight. Clear, ing, involving multiple scattering and aerosol scattering, cloudless, low-aerosol conditions should yield high lin- Mueller Matrices Using Sky Polarization 3 ear polarization amplitudes and small deviations in the polarization direction from a Rayleigh model.