## Abstract

We show spectroscopic Mueller-matrix data measured at multiple incidence angles of the scarab beetle *C. aurata*. A method of regression decomposition can decompose the Mueller matrix into a set of two matrices representing one polarizer and one dielectric reflector. We also report on a tentative decomposition of the beetle *C. argenteola* using the same method.

© 2016 Optical Society of America

## 1. INTRODUCTION

A Mueller matrix is a complete description of the linear optical response of a sample when light is reflected specularly off it or transmitted through it [1]. With Mueller-matrix spectroscopic ellipsometry (MMSE), the spectral and angle-of-incidence dependence of a Mueller matrix for a sample can be measured accurately [2] and can be used to extract optical and structural properties of samples composed of complex nanostructured materials. Mueller matrices are based on a Stokes vector description of light and thus have the benefit of including interaction not only with polarized light but also with unpolarized and partially polarized light. MMSE is frequently used in the analysis of bulk and thin film structures, and during the last decade it has also been employed to study natural photonic structures [3,4]. These results may be helpful in biomimicry as nature has long been a source of inspiration for developing and refining technologies in materials science. Structural coloring may serve as one example, and mimicking of structures of *morpho* butterflies could lead to environmentally friendly long-lasting colors [5]. In this report, we present studies on beetles in the Scarabaeidae family. The selected beetles show brilliant colors and interesting polarization features as well. Mueller matrices measured on such beetles are of large interest for exploring biomimetics and understanding the biological relevance of the observed polarization phenomena. Several species in the Scarabaeidae family have been studied by Hodgkinson *et al.* [6], Goldstein [3], and our group [4,7], to mention a few. Ellipticity, degree of polarization, and other derived parameters have been reported [4], and Arwin *et al.* [7] also did optical modeling to determine structural parameters of the scutellum part of the exoskeleton of *Cetonia aurata*.

Mueller matrices are very rich in information about the sample’s properties and can also be analyzed by addressing depolarization. Cloude [8] showed that a depolarizing Mueller matrix can be represented by a sum of up to four nondepolarizing Mueller matrices weighted by the eigenvalues of the covariance matrix of the Mueller matrix. These eigenvalues are all positive for a physically realizable Mueller matrix, and this so-called sum decomposition can be used to filter matrices and obtain a measure of experimental fidelity [9]. The result of the decomposition can also be used to describe a Mueller matrix as a set of basic optical elements having direct physical meaning such as polarizers and retarders. Pioneering work on decomposition of Mueller-matrix images, including studies of beetles, was performed by Ossikovski *et al.* [10]. We have also previously [11] demonstrated this with Cloude as well as regression decomposition of Mueller-matrix spectra and images measured at near-normal incidence on *C. aurata*. Using Cloude decomposition, we found that the experimentally determined Mueller matrix of *C. aurata* at near-normal incidence decomposes into a set of a mirror and a circular polarizer. Those results were then the basis for a more stable regression decomposition where the result was confirmed.

Our objective in this work is to further develop the parameterization of Mueller matrices of beetles in terms of basic optical elements using regression decomposition. We generalize the previous decomposition by including angle of incidence depending Mueller matrices, which leads to more complicated decompositions.

## 2. EXPERIMENTAL DETAILS

We have measured Mueller matrices of two species of scarab beetles: *Cetonia aurata* (Linnaeus, 1758) collected on the island of Öland in Sweden and *Chrysina argenteola* (Bates, 1888), which is on loan from the Natural Museum of History in Stockholm, Sweden, but originates from Colombia. MMSE measurements were performed on the scutellum, which is a small, triangular, and relatively flat area on the dorsal side of the beetle. The incidence angles were in the range of 15°–55°. The ellipsometer used is a dual rotating compensator ellipsometer (RC2) from J.A. Woollam Co., Inc. that is equipped with focusing optics to reduce the spot size to $<100\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{\mu m}$. The spectral range of the instrument is 245–1690 nm, but only data in the range of 245–1000 nm are presented here. All calculations regarding decomposition were made using the lsqcurvefit function in Matlab.

## 3. THEORY

When the polarization state of light is described with a Stokes vector, the incoming and outgoing waves are related by a $4\times 4$ Mueller matrix $\mathbf{M}$ as [1]

where ${\mathbf{S}}_{\mathbf{i}}$ and ${\mathbf{S}}_{\mathbf{o}}$ are Stokes vectors for the incoming and the outgoing beams of light, respectively. If we expand the equation, we get#### A. Regression Sum Decomposition

Cloude has shown that a depolarizing Mueller matrix can be represented as a sum of up to four nondepolarizing matrices [8]. As an alternative to Cloude’s method, regression decomposition can be used, as recently shown for spectral data and Mueller-matrix images on *C. aurata* [11]. An experimentally determined normalized Mueller matrix $\mathbf{M}$ is then decomposed into a sum of Mueller matrices ${\mathbf{M}}_{i}^{\text{reg}}$ representing optical devices such as polarizers, retarders, and rotators

#### B. Mueller Matrices for Rotated Polarizing Components

If an optical device is rotated an angle $\alpha $ counterclockwise with respect to the $x$ axis in a reference Cartesian coordinate system when looking into the source, the Mueller matrix ${\mathbf{M}}^{\prime}$ of the rotated device will be

where $\mathbf{M}$ is the Mueller matrix of the unrotated device, and $\mathbf{R}(\alpha )$ is given by [14]#### C. Polarizers

A polarizer will attenuate two orthogonal components of the electric field in a plane wave by different amounts. An ideal linear polarizer with its transmission axis in the $x$ direction is given by [1]

*et al.*[18]) and we obtain

#### D. Isotropic Reflector

An isotropic reflector used at oblique incidence introduces retardation $\delta $ and diattenuation $\mathrm{cos}\text{\hspace{0.17em}}2\psi $ between the electric field components in two orthogonal directions parallel and perpendicular to the plane of incidence [2]

## 4. RESULTS AND DISCUSSION

#### A. Regression Sum Decomposition of Mueller Matrices of C. Aurata

Mueller matrices measured on *C. aurata* at incidence angles in the range of 15° to 55° with a step size of 1° in the spectral range from 245 to 1000 nm are summarized in a contour plot in Fig. 1. Below 500 nm and above 600 nm, all off-diagonal block elements as well as ${m}_{34}$ and ${m}_{43}$ are almost constant and close to zero, and the rest of the elements have monotonously changing values with increasing incidence angle. Between 500 and 600 nm, however, many elements show dramatic changes. In this region, the Mueller matrix is depolarizing and can therefore be sum decomposed according to Cloude [8], as shown in our previous report [11]. We call this region the circularly polarizing regime and we will refer to it as the ${P}_{C}$-regime. Within the ${P}_{C}$-regime, the cuticle will reflect incident unpolarized light as near-circular polarized at small angles of incidence $\theta $ with an ellipticity decreasing with $\theta $. In our previous work based on near-normal incidence data, we used an ansatz with a decomposition into a circular polarizer and a dielectric mirror. A circular polarizer is a special case of an elliptical polarizer, and a dielectric mirror can be represented by a retarder with $\delta =180\xb0$ and with the conventions used. Here, we expand the circular polarizer to an elliptical polarizer and the mirror to an isotropic reflector, according to Eqs. (10) and (13), respectively.

We apply a regression analysis according to Eqs. (3) and (4) with the ansatz

### 1. Outside the ${\mathsf{P}}_{\mathsf{C}}$-regime

Outside the ${P}_{C}$-regime, the contribution from the elliptical polarizer to the optical properties is approximately zero, and the values of $\epsilon $ and $\alpha $ in Figs. 2(b) and 2(c), respectively, are therefore of little relevance.

The values of $\psi $ and $\delta $ of the isotropic reflector can be seen in Figs. 2(e) and 2(f), respectively. The value of $\psi $ is a function of wavelength and varies from 43° at near-normal incidence to 0° at larger incidence angles, and $\delta $ is constant at 180°.

When the values at near-normal incidence are put into Eq. (13), the isotropic reflector becomes a mirror as in Eq. (14), and for large incidence angles the same isotropic reflector is a linear polarizer described in Eq. (15).

### 2. Inside the ${\mathsf{P}}_{\mathsf{C}}$-regime

The optical properties inside the ${P}_{C}$-regime can be described by a sum of the elliptical polarizer and the retarder. Here, the values of $\psi $ and $\delta $ of the retarder show some small variations, as seen in Figs. 2(e) and 2(f). Compared to the value outside the ${P}_{C}$-regime, there is a small increase in $\psi $ of approximately 10° for incidence angles above 35° in a very narrow band in the center of the ${P}_{C}$-regime. $\delta $ varies from 180° down to 160° for incidence angles below approximately 40°, and it drops to 110° in the same narrow band in the center of the ${P}_{C}$-regime close to 55°.

The angle $\epsilon $ and the azimuthal rotation $\alpha $ of the polarizer both have a strong dependence on wavelength and oscillate rapidly if viewed at a fixed angle of incidence. $\epsilon $ [Fig. 2(b)] varies from approximately 130°–150° with a tendency to lower values at smaller incidence angles. $\alpha $ [Fig. 2(c)] shows small variations between 30° and 60° throughout the ${P}_{C}$-regime. When these values are put into Eq. (12) and multiplied by the weight factor, coefficient $a$, the matrix (Fig. 3) will be similar to that of a left-handed circular polarizer. However, all elements will have non-zero values inside the ${P}_{C}$-regime, particularly in rows 1 and 4 and in columns 1 and 4. This deviation from a pure circular polarizer is likely due to in-plane anisotropy in the sample under inspection.

The values from Figs. 2(e) and 2(f) are put into Eq. (13) and multiplied by its weight factor $b$, and the resulting matrix can be seen in Fig. 4. Here, we see very small values inside the ${P}_{C}$-regime, and outside it has the characteristics of a reflection from a dielectric material, i.e., Eq. (13) with $\delta =180\xb0$.

When the matrices in Figs. 3 and 4 are added according to Eq. (16), the result is very close to the experimental matrix. The differences between the experimental data and ${\mathbf{M}}^{\text{reg}}$ after decomposition can be seen in Fig. 5. Given the complexity of the biological reflector, the agreement is generally very good, and the differences between the experimental data and the result of the regression is within $\pm 0.08$ for all $\theta $ and $\lambda $. From Figs. 2(a) and 2(d), the conclusion can be drawn that outside the ${P}_{C}$-regime the polarizing properties for all $\lambda $ and $\theta $ of *C. aurata* are completely described by an isotropic reflector with $\psi $, according to Fig. 2(e) and $\delta =180\xb0$. Inside the ${P}_{C}$-regime, the situation is slightly more complex, and the polarizing properties cannot be completely described by a single, easily identifiable Mueller matrix. Here, the properties are described by a sum of an elliptical polarizer and an isotropic reflector where the contribution from the reflector is only in the diagonal of the Mueller matrix and by a small amount. The main contribution comes from the polarizer, where the values are such that it is close to a circular polarizer. However, the elements ${m}_{12}$, ${m}_{13}$, ${m}_{21}$, ${m}_{24}$, ${m}_{31}$, ${m}_{34}$, ${m}_{42}$, and ${m}_{43}$ are non-zero, which indicates that some anisotropy in the surface-plane is involved. The sensitivity to such anisotropy depends on the spot size, which will increase along the plane of incidence with larger incidence angles. This has been discussed in previous work [11], and the effect is seen in Fig. 5 as an increasing value in the matrix elements with increasing incidence angle.

#### B. Regression Sum Decomposition of Mueller Matrices of C. Argenteola

Mueller matrices of *C. argenteola* were recorded at incidence angles in the range of 15°–55° with a step size of 1° in the spectral range of 245–1000 nm and are summarized in a contour plot in Fig. 6. The broadband reflector *C. argenteola* exhibits left-handed polarization features for small $\theta $, as seen as the large blue regions in Mueller matrix elements ${m}_{14}$ and ${m}_{41}$, and right-handed polarization features for large $\theta $ are seen as the small red regions close to 55° in the same elements. A notable feature in Fig. 6 is also the interference oscillations that are seen in most elements as rapid variations in wavelength dependence. We here present only rudimentary trials on regression decomposition of *C. argenteola* using the same ansatz as was used for *C. aurata* [i.e., Eq. (16)]. The resulting parameters are presented in Fig. 7. However, the fit is not satisfactory as the Frobenius norm using Eq. (4) has values up to 0.3 in some regions for several of the elements and up to 0.5 in element ${m}_{44}$. We also find that a Cloude decomposition results in three (or more depending on the incidence angle; see our previous work [11]) non-zero eigenvalues for *C. argenteola* compared to only two for *C. aurata*, as seen in Fig. 8. This means that more than two matrices should be included in a regression analysis and not only two, as was used to derive the parameters in Fig. 7. One reason for the larger complexity in the analysis of *C. argenteola* data is that the pronounced interference patterns seen both in primary data in Fig. 6 and in eigenvalues in Fig. 8 show that its cuticle is more transparent compared to that of *C. aurata*. The probe depth is therefore larger, and structures deeper inside the cuticle are probed. In a broadband reflector like *C. argenteola*, the Mueller matrix then becomes more complex as structures have different pitches at different depths. Similar effects have been observed and analyzed in the beetle *Cotinis mutabilis* [19]. Work is in progress to determine additional matrices suitable for the regression decomposition of *C. argenteola* data.

#### C. Final Remarks

We have presented a regression-based decomposition analysis of Mueller matrices measured as a function of wavelength and incidence angle. A sum of an elliptical polarizer and an isotropic reflector has been fitted to these experimental Mueller matrices. In the case of *C. aurata*, this works well. In the case of *C. argenteola*, however, the ansatz presented in Eq. (16) does not work, which is explained by the fact that a Cloude decomposition reveals that additional matrices are needed in the ansatz, and further investigations will be performed.

## Funding

Swedish Government Strategic Research Area in Materials Science on Functional Materials at Linköping University (Faculty Grant SFO Mat LiU 2009 00971); Vetenskapsrådet (VR) (621-2011-4283); Knut och Alice Wallenbergs Stiftelse (2004.0233); Carl Tryggers Stiftelse för Vetenskaplig Forskning (CTS12:31).

## Acknowledgment

The authors thank Jan Landin for providing specimens of *C. aurata* and the Swedish Museum of Natural History for lending a specimen of *C. argenteola*.

## REFERENCES

**1. **E. Collett, *Polarized Light: Fundamentals and Applications* (Marcel Dekker, 1993).

**2. **H. Fujiwara, *Spectroscopic Ellipsometry: Principles and Applications* (Wiley, 2007).

**3. **D. H. Goldstein, “Polarization properties of scarabaeidae,” Appl. Opt. **45**, 7944–7950 (2006). [CrossRef]

**4. **H. Arwin, R. Magnusson, J. Landin, and K. Järrendahl, “Chirality-induced polarization effects in the cuticle of scarab beetles: 100 years after Michelson,” Philos. Mag. **92**(12), 1583–1599 (2012). [CrossRef]

**5. **A. Saito, “Material design and structural color inspired by biomimetic approach,” Sci. Technol. Adv. Mater. **12**, 064709 (2011). [CrossRef]

**6. **I. Hodgkinson, S. Lowrey, L. Bourke, A. Parker, and M. W. McCall, “Mueller-matrix characterization of beetle cuticle: polarized and unpolarized reflections from representative architectures,” Appl. Opt. **49**, 4558–4567 (2010). [CrossRef]

**7. **H. Arwin, T. Berlind, B. Johs, and K. Järrendahl, “Cuticle structure of the scarab beetle Cetonia aurata analyzed by regression analysis of Mueller-matrix ellipsometric data,” Opt. Express **21**, 22645–22656 (2013). [CrossRef]

**8. **S. R. Cloude, “Conditions for the physical realisability of matrix operators in polarimetry,” Proc. SPIE **1166**, 177–187 (1990). [CrossRef]

**9. **R. Ossikovski, M. Anastasiadou, S. Ben Hatit, E. Garcia-Caurel, and A. De Martino, “Depolarizing Mueller matrices: how to decompose them?” Phys. Status Solidi A **205**, 720–727 (2008). [CrossRef]

**10. **R. Ossikovski, M. Foldyna, C. Fallet, and A. De Martino, “Experimental evidence for naturally occurring nondiagonal depolarizers,” Opt. Lett. **34**, 2426–2428 (2009). [CrossRef]

**11. **H. Arwin, R. Magnusson, E. Garcia-Caurel, C. Fallet, K. Järrendahl, M. Foldyna, A. De Martino, and R. Ossikovski, “Sum decomposition of Mueller-matrix images and spectra of beetle cuticles,” Opt. Express **23**, 1951–1966 (2015). [CrossRef]

**12. **J. J. Gil, “Polarimetric characterization of light and media,” Eur. Phys. J. **40**, 1–47 (2007). [CrossRef]

**13. **G. H. Golub and C. F. V. Loan, *Matrix Computations*, 3rd ed. (JHU, 1996).

**14. **R. M. A. Azzam and N. M. Bashara, *Ellipsometry and Polarized Light* (Elsevier, 1999).

**15. **D. S. Kliger, J. W. Lewis, and C. E. Randall, *Polarized Light in Optics and Spectroscopy*, 2nd ed. (Academic, 1990).

**16. **H. C. van de Hulst, *Light Scattering by Small Particles* (Dover, 1981), reprinted.

**17. **L. Li, “Symmetries of cross-polarization diffraction coefficients of gratings,” J. Opt. Soc. Am. A **17**, 881–887 (2000). [CrossRef]

**18. **E. Garcia-Caurel, R. Ossikovski, M. Foldyna, A. Pierangelo, B. Drevillon, and A. De Martino, “Advanced Mueller ellipsometry instrumentation and data analysis,” in *Ellipsometry at the Nanoscale*, M. Losurdo and K. Hingerl, eds. (Springer, 2013).

**19. **A. Mendoza-Galván, E. Muñoz-Pineda, K. Järrendahl, and H. Arwin, “Evidence for a dispersion relation of optical modes in the cuticle of the scarab beetle *Cotinis mutabilis*,” Opt. Mater. Express **4**, 2484–2496 (2014). [CrossRef]