OSA's Digital Library

Optics Express

Optics Express

  • Editor: J. H. Eberly
  • Vol. 9, Iss. 11 — Nov. 19, 2001
  • pp: 567–578
« Show journal navigation

Morpho butterflies wings color modeled with lamellar grating theory

Boris Gralak, Gérard Tayeb, and Stefan Enoch  »View Author Affiliations


Optics Express, Vol. 9, Issue 11, pp. 567-578 (2001)
http://dx.doi.org/10.1364/OE.9.000567


View Full Text Article

Acrobat PDF (880 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We describe an approach for converting reflection coefficients of any structure into colors, taking into account human color perception. This procedure is applied to the study of the colors reflected by Morpho rhetenor butterflies wings. The scales of these wings have a tree-like periodic structure which is modeled with the help of a rigorous lamellar grating electromagnetic theory. In this way, we are able to determine the colors reflected by the wing under various illumination conditions.

© Optical Society of America

1. Introduction

The electromagnetic modeling of the reflection by a given structure generally provides reflection coefficients which depend on several parameters: wavelength, polarization, observation angles. The color of the object also depends on the light source, and must take into account human color perception. In the first part of this paper, we present a convenient way to convert these reflection coefficients into colors that can easily be displayed on a computer screen, then printed.

2. From reflection coefficient to color

In this section, our objective is to summarize the basic ideas that we have used for converting reflection coefficients into colormaps. Basically, we have picked out and put together several methods well known in color science and graphic imaging. We will therefore not give all the details, and this section should be considered as a collection of recipes. For a quick illustration of the main idea, we recommend the spectrum applet in Ref. [3

3. Educational color applets, The Spectrum Applet, http://www.cs.rit.edu/~ncs/color/

].

Let us suppose that the body under study is illuminated by an illuminant characterized by its energy distribution D(λ). In this paper, we use the CIE (Commission Internationale de l–Éclairage) normalized illuminant D65, which closely matches that of the sky daylight [4

4. Colourware, Colour FAQ, http://www.colourware.co.uk/cpfaq.htm

6]. For convenience, we have used a 5th order polynomial approximation of this function (Fig.1a). Assuming that the body has a reflectivity R(λ), we can compute the CIE XYZ tristimulus values [7

7. Charles A. Poynton, Color technology, http://www.inforamp.net/~poynton/ColorFAQ.html

] which are one possible way to characterize the reflected color, according to:

X=1kD(λ)R(λ)x¯(λ)dλ
Y=1kD(λ)R(λ)y¯(λ)dλ
Z=1kD(λ)R(λ)z¯(λ)dλ
(1)

where the weighting functions , ȳ, (Fig.1b) are the CIE 1964 spectral tristimulus values [5

5. CVRL Color & Vision database, http://www.cvrl.org/

,6,8] and k is a normalization factor defined in such a way that an object with a uniform reflectivity R(λ)=1 gives a luminance component Y equal to 1:

k=D(λ)y¯(λ)dλ
(2)

In the numerical treatment, all the integrals are replaced by discrete sums. Due to the human eye spectral bandwidth, the only values of the wavelength that we will retain lie in the range 0.38–0.78 µm. We use 11 equally spaced sampling values in this range, which are enough to get convergence in the color rendering.

Fig. 1. a: Relative energy distribution of the D65 illuminant (blue) and its 5th order polynomial approximation (red) versus the wavelength λ. b: Spectral tristimulus values.

We then transform the XYZ components into RGB components, which are more convenient for the imaging process. Since RGB values are device dependent, we use the BT.709 Recommendation, which is representative of contemporary monitors in studio video, computing and computer graphics. In this case, the conversion is given by the linear transformation [7

7. Charles A. Poynton, Color technology, http://www.inforamp.net/~poynton/ColorFAQ.html

]:

[RGB]=[3.2404791.5371500.4985350.9692561.8759920.0415560.0556480.2040431.057311][XYZ]
(3)

These RGB values are normalized and should stay between 0 and 1. In fact, they can in some cases be slightly less than 0 or greater than 1. It means that the color can not be accurately represented on such monitor devices by additive reproduction of these RGB values. In this case, the values are rounded to 0 or 1.

The last step is the color representation. For this purpose, the PPM (Portable PixMap) file format [6] is very simple to use: it is a text file in which one simply enter the RGB values for each pixel. Note that each RGB value should be an integer between 0 and 255, obtained by scaling the values obtained from (3). These PPM files can be easily displayed on a computer screen using classical image editors, and easily printed.

At this stage, we have replaced a complicated set of data (i.e. the spectrum R(λ)) by a color point in a colormap. As illustrated in the following sections, this representation can be very useful for obtaining a global insight on the reflection properties of a body. A two-dimensional colormap based on this representation provides two free parameters (for example the two angles of diffraction).

3. Modeling the wing structure

Fig. 2. a: Transmission electron microscope image showing the cross-section through a single Morpho rhetenor scale. Image reprinted from [1], with permission from P. Vukusic and the Royal Society. b: Modeled structure; the two red lines define a grating layer, the optical index of black regions is n, and that of white regions is 1.
Fig. 3. a: k i is the incident wave vector, with angles θ i and φ i . s i and p i are respectively normal and parallel to the incidence plane Π i . The polarization of the incident electric field E i is defined by the angle δi . b: Schematic representation of two conical subregions of the upper hemisphere.

Let us focus on the following problem: we suppose that the butterfly wing is illuminated by the skylight, and we want to determine the color reflected by the wing in the direction (θ rr ). Of course, it is a particular problem and, using the same technique, we could also model the color reflected by the wing if it was for instance illuminated by a focussed beam. Since the wing is opaque and supposed to lie in the horizontal plane, we only need to consider the incident light coming from the upper halfspace. We divide this halfspace in conical subregions centered around a direction (θ,φ), with solid angle ΔΩθ,φ (Fig.3b). In order to model the uniform luminance of the skylight, we consider that two plane waves impinge on the wing with a wave vector k i colinear to the axis of each cone. These two plane waves have orthogonal polarizations (δi =0 and 90°) since daylight is unpolarized. Depending on θ i , φ i and the grating period, they produce one or several reflected orders flowing inside one or several solid angles ΔΩθnr,φnr,. For each wavelength λ, and each conical subregion centered around a reflected direction (θ rr ), we add incoherently the contributions of all the incident waves giving a reflected order in this subregion, regardless of polarization, and obtain the power reflected in this subregion, which we denote by Rθr,φr(λ). Using the recipe described in section 2 where R(λ) is now replaced by Rθr,φr(λ), we obtain a color (including the brightness information) which can be denoted by a function RGB(θ rr ), and which can be easily plotted.

Note for instance that if the halfspace is divided in conical subregions with 10° aperture (i.e. 9 values of θ i and 36 values of φ i , with 2 values of δi and 11 values of λ), one has to compute 9×36×2×11=7128 grating problems). Our computer code, running on a Pentium IV 1.4 GHz computer, deals this problem in about 18 minutes.

4. Numerical results

We first need to obtain convenient values for the various geometrical dimensions involved in the Morpho rhetenor scales. For this purpose, we have used the image of Fig.2a, and determined the scaling factor in order to match some measured values reported in Ref. [1

1. P. Vukusic, J.R. Sambles, C.R. Lawrence, and R.J. Wootton, “Quantified interference and diffraction in single Morpho butterfly scales,” Proceedings: Biological Sciences, The Royal Society of London 266, 1403–1411 (1999). [CrossRef]

]. The grating period is determined in order to have the same diffraction directions as in Fig.4a of [1

1. P. Vukusic, J.R. Sambles, C.R. Lawrence, and R.J. Wootton, “Quantified interference and diffraction in single Morpho butterfly scales,” Proceedings: Biological Sciences, The Royal Society of London 266, 1403–1411 (1999). [CrossRef]

], the scales being surrounded by air as well as by isopropyl alcohol (IPA), and illuminated in normal incidence with λ=488 nm. Since the scale structure differs from a real periodic grating, some diffusion is added and enlarges the diffraction peaks, which makes the determination of the grating period inaccurate. A more accurate value is obtained by matching the spectral behavior of the scale reflectivity (Fig.5a of [1

1. P. Vukusic, J.R. Sambles, C.R. Lawrence, and R.J. Wootton, “Quantified interference and diffraction in single Morpho butterfly scales,” Proceedings: Biological Sciences, The Royal Society of London 266, 1403–1411 (1999). [CrossRef]

]). Combining these two fits, we obtain a grating period of 746 nm, which will be used throughout this section. The other characteristic dimensions are shown on Fig.4.

Fig. 4. Characteristic dimensions of the modeled scale (in nanometers). The white regions are filled with air (or with a liquid solvent), and the gray regions are filled with a material with optical index n. In this figure, and using the tree image, we will say that each tree has four branches. The lower layer models the bottom membrane of the scale. Apart from the bottom air layer, all the dimensions are scaled.

Using a structure with the same dimensions as those given in Fig.4, but with 16 branches for each tree, we compute the diffraction efficiencies (which are the fractions of the incident energy diffracted in the various grating orders) when the scale is illuminated in normal incidence by a plane wave with λ=488 nm. We take an average value of the efficiencies obtained for the two fundamental polarizations in order to simulate an unpolarized incident light. The computations are made when the scale is in air and when it is immersed in IPA (with optical index 1.38) and are plotted on Fig.5. When the scale is in air, the main part of the diffracted light is reflected in the ±1 orders, whereas when the scale is in IPA, the index contrast between the scale and the external medium is weaker, and almost all the light is transmitted in the 0 order. Note that the missing part of the energy is absorbed by the structure.

Fig. 5. Diffracted efficiencies by a scale in air and in IPA. The reflected angles are between 0 and 180° and the transmitted angles are between 180 and 360°. The normals to the scale correspond to 90° and 270°.

Figure 6 shows the total reflection and transmission of the same structure (16 branches, lying in air) illuminated by a plane wave in normal incidence. These total reflection and transmission coefficients are obtained by adding all the reflected or transmitted efficiencies, for the two fundamental polarizations. We call E‖ (resp. H‖) the case where the electric (resp. magnetic) field is parallel to the scale ridges (x-axis in Fig.2b). Note that these cases are respectively denoted by TM and TE in Ref.[1

1. P. Vukusic, J.R. Sambles, C.R. Lawrence, and R.J. Wootton, “Quantified interference and diffraction in single Morpho butterfly scales,” Proceedings: Biological Sciences, The Royal Society of London 266, 1403–1411 (1999). [CrossRef]

]. The reflection curves are very similar to the experimental data reported in the Fig.5 of Ref.[1

1. P. Vukusic, J.R. Sambles, C.R. Lawrence, and R.J. Wootton, “Quantified interference and diffraction in single Morpho butterfly scales,” Proceedings: Biological Sciences, The Royal Society of London 266, 1403–1411 (1999). [CrossRef]

], and show a pronounced peak of reflectivity for the smaller wavelengths, giving to the wing its blue color. Note that the transmission curves predict that a sizeable amount of energy is transmitted for larger wavelengths, which does not agree with the experience. This phenomenon is quite general in our modeling: we generally find a good agreement where the amount of energy is significant, and a poor agreement where the amount of energy is low. This discrepancy can be due to several reasons, and we give below those which appear to be the most relevant.

First, we can observe that when the amount of energy is significant, the phenomenon is well described by a periodic structure and its interference process, whereas when the amount of energy is low, the diffusion by the random characteristics of the structure (which is not taken into account here) seems to give a significant contribution.

Secondly, the wing structure in the region between the ridges and the bottom membrane is actually rather complex and disordered, with probable inclusions of pigments. It will be shown below (Fig. 13) that for shorter wavelengths the field does not go deeply inside the structure, contrary to longer wavelengths. It means that the longer wavelengths are much more sensitive to the bottom region of the scale. If the absorption phenomenon in the actual wing is mainly due to localized effects inside the pigments (or other absorbing defaults) located in this bottom region, it explains the discrepancy between our computations and the experimental data for longer wavelengths. Note that in our model, the absorption is only due to the imaginary part of the ridges’ permittivity.

Nevertheless, since the color is generated by the stronger reflected components, we can be confident in the relevance of our approach as regards the color rendering.

Fig. 6. Total reflection and transmission for both polarizations.

In order to present the color appearance of the scale, we now implement the approach exposed in section 2. We suppose that the wing is illuminated by a wide source modeling the sky, with the spectral energy distribution of the D65 illuminant, and assume that the incident light is coming from the whole upper halfspace. We divide this halfspace in conical subregions centered around directions (θ,φ) separated with angular variations of 10° degrees, then we collect and add the energies of the diffracted orders reflected in each conical subregion centered around a reflected direction (θ rr ). Finally, we transform this information in a RGB(θ rr ) map.

Figure 7 shows this map for the structure described previously (16 branches, lying in air). On these maps, the horizontal axis gives the θ r value, from 0 to 80°, and the vertical axis gives the φ r value, from 0 to 180°. Note that the definition of θ r and φ r follows that of Fig.3a. It means that the left side of the maps gives the reflected color in the direction normal to the scale (θ r =0) and the right side gives the reflected color for near grazing angles.

In figure 7, the left map color scale intensity goes from black (nothing reflected) to almost white (in fact the D65 illuminant color) if all wavelengths were totally reflected. This map appears dark since the reflection is generally much lower than unity. Since the human eye can distinguish between much more different levels of illumination than a computer monitor or a printer can render, this dark map do not really agree with the visual perception. In the right map, the global brightness is enhanced in such a way that the greater R, G or B value among all the map reaches the maximum value of 255. This map retains the relative brightness information between square subregions and agrees better with the visual perception.

Fig. 7. Reflected colors by a scale with 16 branches in air. The horizontal axis gives the θ r value, from 0 to 80°, and the vertical axis gives the φ r value, from 0 to 180°. The angular shift between two neighboring squares is equal to 10°.

When the scale is immersed in IPA, the index contrast between the scale and the external medium changes, thus producing drastic modifications in the reflected colors. Figure 8 shows that in this case, the global appearance of the scale becomes green, as reported in [16

16. OLE (Opto & Laser Europe), feature article, June 1998 issue.

,17

17. S. Berthier, Les couleurs des papillons ou l’impérative beauté, Springer-Verlag France (2000).

]. We can also notice that the total amount of reflected energy is much lower than when the scale is in air.

Fig. 8. Same as Fig.7, but the scale is immersed in IPA.

Fig. 9. Characteristic dimensions of the second modeled scale (in nanometers). The nine lower lamellae have the same dimensions as those described in Fig.4, and the upper lamellae have decreasing widths. Apart from the bottom air layer, all the dimensions are scaled.

The colors reflected by the structure of Fig.9 are represented in Figs.10 and 11. Comparison with Figs. 7 and 8 shows that the amount of reflected energy is lower than when all the lamellae have the same size. This feature appears clearly on the maps when the scale is in air. When the scale is in IPA, both left maps are so dark that the difference is not obvious, but the difference can be seen on the data from which the maps are obtained. Anyway, the overall colors are almost the same for both structures.

Fig. 10. Reflected colors by the structure of Fig.9 lying in air.
Fig. 11. Reflected colors by the structure of Fig.9 lying in IPA.

Another question that arises is linked to the left-right dissymmetry of the structure. Since both structures (Fig.4 or Fig.9) are not symmetric, it is interesting to know if the reflected colors are different when the incident light comes from the left or from the right. To answer this question, we have modeled a collimated beam with mean incidence (see Fig.3a) corresponding to θ i =45°, and φ i =0 (incidence from the “left”) or φ i =180° (incidence from the “right”). The incident light comes in a solid angle which spreads 10° around this mean incidence, and the reflected light is collected in a solid angle which spreads 10° around the normal (θ r =0). We have used the structure of Fig.4, with 16 branches for each tree, and lying in air. In this case, there is no significant difference between the reflected colors for φ i =0 and φ i =180° : they both give the same blue shade.

The reason why this structure behaves the same way when illuminated from the left or from the right is illustrated on Fig.12. Here, the structure is illuminated by a plane wave with an incidence angle θ i =45°, and we study the two cases where φ i =0 or φ i =180°. In these conditions, one reflected order escapes near the normal to the structure, and that order carries the most part of the reflected energy. The figure shows the averaged efficiencies between the two fundamental polarizations (δi =0 and 90°), in order to model an unpolarized incident light. The two curves are almost identical, with a pronounced peak for the shorter wavelengths, thus generating the same blue color.

Fig. 12. Reflected energy near the normal when the structure with 16 branches in air is illuminated from the left (φ i =0) and from the right (φ i =180°).

Finally, Fig.13 shows field maps when the structure of Fig.9 is illuminated in normal incidence by a plane wave. For λ=0.42 µm, it shows that the field does not go deeply inside the structure, and is reflected by the scale. For λ=0.60 µm, we see that the field propagates inside the structure following the dielectric regions, thus explaining a larger absorption due to the strong field value inside the lossy medium.

Fig. 13. Field maps of the squared modulus of the total electric field for the structure of Fig.9 illuminated by a E//polarized plane wave in normal incidence.

5. Conclusion

The understanding of the iridescent colors reflected by butterflies wing is an attractive subject which has aroused several studies in the recent years. We think that the approach exposed in this paper brings two new interesting points of view.

First, the way to convert reflection coefficients into colors provides a smart technique able to bring together many numerical data into easily interpretable colormaps. Since all the content of a reflection spectrum R(λ) can be replaced by a color point in a colormap, such a colormap offers easily two free parameters in a 2D plot (and even more when using more sophisticated plottings). Such representations can be very useful to get a global insight on the reflection properties of a body.

Using these techniques, we have recovered several properties observed experimentally, such as the color reflected by the scale lying in air or in IPA, or the spatial and spectral distributions of the energy diffracted by the wing.

At the present time, we have not considered the fact (clearly apparent in Fig. 2) that the ridges are tilted (toward the left on the Fig. 2). This asymmetry could be taken into account using our modeling, using a staircase approximation of the tilted ridges. It would lead to much larger computation times, but such a study should be affordable.

Finally, our approach could be used for various biologic species possessing an almost periodic micrometric structure.

References and links

1.

P. Vukusic, J.R. Sambles, C.R. Lawrence, and R.J. Wootton, “Quantified interference and diffraction in single Morpho butterfly scales,” Proceedings: Biological Sciences, The Royal Society of London 266, 1403–1411 (1999). [CrossRef]

2.

H. Tada, S. Mann, I. Miaoulis, and P. Wong, “Effects of a butterfly scale microstructure on the iridescent color observed at different angles,” Opt. Express 5, 87 (1999), http://epubs.osa.org/oearchive/source/11782.htm. [CrossRef] [PubMed]

3.

Educational color applets, The Spectrum Applet, http://www.cs.rit.edu/~ncs/color/

4.

Colourware, Colour FAQ, http://www.colourware.co.uk/cpfaq.htm

5.

CVRL Color & Vision database, http://www.cvrl.org/

6.

Color resources, http://institut.fresnel.free.fr

7.

Charles A. Poynton, Color technology, http://www.inforamp.net/~poynton/ColorFAQ.html

8.

Color Science, http://www.physics.sfasu.edu/astro/color.html

9.

H. Ghiradella, “Light and color on the wing: structural colors in butterflies and moths,” Applied Optics 30, 3492–3500 (1991). [CrossRef] [PubMed]

10.

P. Vukusic, J.R. Sambles, and H. Ghiradella, “Optical classification of microstructure in butterfly wing-scales,” Photonics Science News 6, 61–66 (2000).

11.

P. Vukusic, J. R. Sambles, C.R. Lawrence, and R.J. Wootton, “Structural colour: now you see it - now you don’t,” Nature 410, 36 (2001). [CrossRef] [PubMed]

12.

L.C. Botten, M.S. Craig, R.C. McPhedran, J.L. Adams, and J.R. Andrewartha, “The finitely conducting lamellar diffraction grating,” Optica Acta 28, 1087–1102 (1981). [CrossRef]

13.

S.E. Sandström, G. Tayeb, and R. Petit, “Lossy multistep lamellar gratings in conical diffraction mountings: an exact eigenfunction solution,” Journal of Electromagnetic Waves and Applications 7, No.5, 631–649 (1993). [CrossRef]

14.

L. Li, “A modal analysis of lamellar diffraction gratings in conical mountings,” J. Modern Optics 40, 553–573, 1993. [CrossRef]

15.

B. Gralak, “Étude théorique et numérique des propriétés des structures à bandes interdites photoniques,” PhD thesis, Université Aix-Marseille (2001), http://institut.fresnel.free.fr.

16.

OLE (Opto & Laser Europe), feature article, June 1998 issue.

17.

S. Berthier, Les couleurs des papillons ou l’impérative beauté, Springer-Verlag France (2000).

OCIS Codes
(050.1950) Diffraction and gratings : Diffraction gratings
(260.2110) Physical optics : Electromagnetic optics
(330.1690) Vision, color, and visual optics : Color

ToC Category:
Research Papers

History
Original Manuscript: October 12, 2001
Published: November 19, 2001

Citation
Boris Gralak, Gerard Tayeb, and Stefan Enoch, "Morpho butterflies wings color modeled with lamellar grating theory," Opt. Express 9, 567-578 (2001)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-9-11-567


Sort:  Journal  |  Reset  

References

  1. P. Vukusic, J.R. Sambles, C.R. Lawrence and R.J. Wootton, "Quantified interference and diffraction in single Morpho butterfly scales," Proceedings: Biological Sciences, The Royal Society of London 266, 1403 -1411 (1999). [CrossRef]
  2. H. Tada, S. Mann, I. Miaoulis, and P. Wong, "Effects of a butterfly scale microstructure on the iridescent color observed at different angles," Opt. Express 5, 87 (1999), http://www.opticsexpress.org/oearchive/source/11782.htm. [CrossRef] [PubMed]
  3. Educational color applets, The Spectrum Applet, http://www.cs.rit.edu/~ncs/color/
  4. Colourware, Colour FAQ, http://www.colourware.co.uk/cpfaq.htm
  5. CVRL Color & Vision database, http://www.cvrl.org/
  6. Color resources, http://institut.fresnel.free.fr
  7. Charles A. Poynton, Color technology, http://www.inforamp.net/~poynton/ColorFAQ.html
  8. Color Science, http://www.physics.sfasu.edu/astro/color.html
  9. H. Ghiradella, "Light and color on the wing: structural colors in butterflies and moths," Appl. Opt. 30, 3492-3500 (1991). [CrossRef] [PubMed]
  10. P. Vukusic, J.R.Sambles, and H. Ghiradella, "Optical classification of microstructure in butterfly wing-scales," Photonics Science News 6, 61-66 (2000).
  11. P. Vukusic, J. R. Sambles, C.R. Lawrence, R.J. Wootton, "Structural colour: now you see it - now you don't," Nature 410, 36 (2001). [CrossRef] [PubMed]
  12. L.C. Botten, M.S. Craig, R.C. McPhedran, J.L. Adams, and J.R. Andrewartha, "The finitely conducting lamellar diffraction grating," Optica Acta 28, 1087-1102 (1981). [CrossRef]
  13. S. E. Sandstr�m, G. Tayeb, R. Petit, "Lossy multistep lamellar gratings in conical diffraction mountings: an exact eigenfunction solution," J. Electromagnetic Waves and Applications 7, 631-649 (1993). [CrossRef]
  14. L. Li, "A modal analysis of lamellar diffraction gratings in conical mountings," J. Modern Opt. 40, 553-573, (1993). [CrossRef]
  15. B. Gralak, "�tude th�orique et num�rique des propri�t�s des structures � bandes interdites photoniques," PhD thesis, Universit� Aix-Marseille (2001), http://institut.fresnel.free.fr.
  16. OLE (Opto & Laser Europe), feature article, June 1998 issue.
  17. S. Berthier, Les couleurs des papillons ou l'imp�rative beaut�, Springer-Verlag France (2000).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited