We investigate the physical mechanism for color production by isotropic nanostructures with short-range order in bird feather barbs. While the primary peak in optical scattering spectra results from constructive interference of singly-scattered light, many species exhibit secondary peaks with distinct characteristic. Our experimental and numerical studies show that these secondary peaks result from double scattering of light by the correlated structures. Without an analog in periodic or random structures, such a phenomenon is unique for short-range ordered structures, and has been widely used by nature for non-iridescent structural coloration.
© 2010 Optical Society of America
Structural colors have attracted much attention in a wide variety of disciplines . They originate from physical interactions of light with nanostructures. Recent studies have focused on ordered structures in the natural world which produce iridescent colors [2, 3, 4]. However, nature has also used extensively quasi-ordered structures to create vivid colors that are weakly iridescent, e.g., the spongy nanostructures in the feather barbs of numerous birds . In the 19th century, lack of pronounced iridescence led to the hypothesis that the colors were produced by wavelength selective scattering from single particles. Raman was the first to cast doubt on this hypothesis , which was later falsified by Dyck . Prum et al. proved the existence of structural order by Fourier analysis of electron micrographs of medullary keratin in feather barbs . In the last ten years, the previously unappreciated class of quasi-ordered nanostructures have been identified also in bird skins, mammal skins, dragonfly cuticles, and butterfly scales [1, 5].
Although it is now recognized that light scattering is organized to contribute effectively to coloration of quasi-ordered structures, there are significant spectral features whose origin remains elusive [5, 9, 10, 11]. In previous studies, optical reflectance spectra are predicted from the structural Fourier spectra with the assumption of single scattering . The dominant length scale of structural correlation gives one peak in the predicted reflectance spectrum. However, many species we have measured feature two reflection peaks, in agreement with prior observations [5, 9, 12]. What is the origin of the second peak? Does multiple scattering play a significant role in coloration? How does the coloration of quasi-ordered structures differ from that of ordered structures? Answering these questions is essential to identify the physical mechanism for coloration of the quasi-ordered structures that are widely distributed in nature and clearly function in animal communication. Moreover, understanding the effects of short-range structural order on multiple scattering of light is fundamental to the physics of light interaction with complex media. In practical terms, the wavelength dependence of light scattering from quasi-ordered structures has not been taken advantage of technologically, and may have useful advantages for photonic coatings across many areas such as the production of color for paint, cosmetics and textile industries and photon management in solar cells.
In this paper, we investigate the physical mechanism for coloration of nanostructures with short-range order in bird feather barbs by performing angle- and polarization-resolved scattering spectrometry. The primary peak in optical scattering spectra maintains the polarization of incident light and its spectral shift with observation angle agrees with single scattering predictions. The short-range structural order leads to phase correlation of scattered waves and their constructive interference produces a narrow spectral peak. For many avian species, the scattering spectra have a secondary peak which is depolarized and has different angular dispersion from the primary peak. We show that the secondary peak is caused by interference of light scattered twice by the nanostructures. This result is unexpected from the general belief that multiple scattering of light simply broadens the single scattering peak. There is no analogous phenomenon in periodic structures where multiple diffraction or high-order diffraction does not create additional peaks.
The color-producing quasi-ordered nanostructures within the medullary cells of bird feather barb rami belong to two morphological classes[5, 13]. Channel-type nanostructures consist of β-keratin bars and air channels in elongate and tortuous forms. Sphere-type nanostructures consist of spherical air cavities in a β-keratin matrix. Our recent studies suggest these nanostructures are self-assembled during phase separation of β-keratin protein from the cytoplasm of the cell . In this work we have performed small angle X-ray scattering (SAXS) and optical measurements on both types of structures from many bird species. Since the results are similar, we present here the data from turquoise wing feathers of Coracias benghalensis (C. benghalensis), the same bird species that Raman studied over seventy years ago [Fig. 1(a)].
The barb nanostructures of C. benghalensis are a 3D network of keratin rods of ~100 nm diameter and air channels of similar width [Fig. 1(b)]. The SAXS pattern exhibits a ring [inset of Fig. 1(b)], implying the nanostructure is isotropic. Exploiting the rotational symmetry, we average the scattering intensity azimuthally to obtain the intensity I as a function of the magnitude of spatial vector q [Fig. 2(a)]. The strong peak of I(q) at q 0 = 0.028 nm-1 reveals the existence of a dominant length scale s for structural correlation. s = 2π/q 0 = 225 nm is approximately half of the wavelength of blue light. The full width at half maximum (FWHM) Δq of the peak in I(q) reflects the range ξ of spatial order, ξ = 2π/Δq = 1.16 mm. Since ξ is only 5 times of the spatial period s, the order is short-ranged. Such structures can be called amorphous photonic materials, in analog to amorphous electronic materials with only short-range order.
We have measured the scattered light from the feather barbs as function of wavelength, polarization, sample orientation, incident angle and observation angle. The cylindrical-shaped feather barbs are mounted at the rotation center of a goniometer, with their axes perpendicular to the rotation axis of the goniometer. Collimated white light is incident on the sample at an angle ϕ from the surface normal. Scattered light is collected by a lens and focused to an optical fiber that is connected to a spectrometer. The measured spectra of scattered light are normalized by the spectrum of incident light. The angular resolution of detection is about 5°. Linear polarizers select the scattered light whose polarization is either parallel (co-polarization) or perpendicular (cross-polarization) to that of incident light. Figure 2(b) shows the spectra of co- and cross-polarized light scattered from the feather barbs. The scattering spectrum of co-polarization has two peaks: the primary peak at 568 nm and the secondary peak at 396 nm. In the scattering spectrum of cross-polarization the primary peak disappears, while the secondary peak remains with reduced intensity. Since it is linearly polarized with the polarization direction parallel to that of the incident light, the primary peak is produced by single scattering process. The secondary peak is depolarized, suggesting it is from high-order scattering process. However, the memory of incident light polarization is not completely lost, as the cross-polarization component is still weaker than the co-polarization one.
To further confirm the origins of the two peaks, we have measured their angular dependences in a series of experiments. In the first experiment, the direction of incident light and the detector position are fixed. As the sample rotates, the frequencies of the primary peak and secondary peak in the scattering spectra do not change. Their independence of the sample orientation is consistent with the isotropic nature of the nanostructures. Next we fix the sample orientation and the incident angle of white light, and move the detector to measure the spectra of light scattered into different directions. Both the primary peak and secondary peak shift in wavelength with the detection angle, but in the opposite directions. We also change the incident angle and repeat the measurement. Analysis of all the data reveals that the frequencies of primary peak and secondary peak do not depend separately on the incident angle or the observation angle, but on their difference, that is, the angle θ between the directions of illumination and observation [inset of Fig. 3(a)]. As shown in Fig. 3(a), the primary peak shifts to shorter wavelength and the secondary peak to longer wavelength when θ increases. We have also measured the reflection spectra of white light for different incident angle ϕ. A comparison to the scattering spectra reveals that the reflection peaks at an incident angle ϕ have the same frequencies as the scattering peaks at angle θ = 2ϕ [Fig. 3(a)]. Since the angle between the incident beam and the reflected beam is θ = 2ϕ, their spectral coincidence implies the scattering peaks and reflection peaks have the same physical origin. Therefore, specular reflection can be regarded as a special case of light scattering from the quasi-ordered structures.
For single scattering of light from the incident wavevector k i to a wavevector k o, k o-k i = q, where q is a spatial vector of the structure. θm is the angle between -k i and k o [inset of Fig. 3 (b)]. The scattering is elastic, ∣k i∣ = ∣k o∣ ≡ k. According to the inset diagram of Fig. 3 (b), the magnitude of q is
where k = 2πne/λ, ne is the effective index of refraction, and λ is the vacuum wavelength. We compute the q values for the primary peak and secondary peak in the measured scattering spectra. Light refraction at sample surface is taken into account by converting the experimental angle θ outside the sample to θm inside the sample . As shown in Fig. 3 (b), the q value for the primary peak of light scattering is the same for different θ. It is equal to q 0 of the SAXS data [Fig. 2(a)] when ne = 1.25. This value corresponds to a volume fraction of air of 57% in the keratin matrix, which matches the estimate from electron micrographs. Therefore, the primary peak in the optical scattering spectra is produced by single scattering process involving the dominant spatial vectors of the nanostructure.
Unlike the primary peak, the q value for the secondary peak changes with θ. We also compare the optical scattering spectra to X-ray scattering results by plotting the nearly backscattered light intensity (θ = 10°) against q in Fig. 2(a). While the primary peak of light scattering matches that of SAXS, the secondary peak finds no correspondence in the SAXS data. This discrepancy is attributed to stronger scattering of light than X-ray by the nanostructures, which results in multiple scattering of light. The lowest order multiple scattering is double scattering, i.e., the incident light is scattered successively twice before exiting the sample. The two single scattering events occur at different locations and are considered independent of each other. As shown schematically in Fig. 4 (a), light is scattered from k i to k m, and then from k m to k o. The two successive single scattering events are considered to be uncorrelated. k m - k i = q 1, and k o - k m = q 2, where q 1 and q 2 are spatial vectors of the structure. Light scattering is elastic, ∣k i∣ = ∣k m∣ = ∣k o∣ = k. When the directions of k i and k o are fixed, k m can be in multiple directions, as illustrated in Figs. 4(a). α 1 and α 2 denote the angles for two successive single scattering processes. As α 1 and α 2 increase, k must increase in order to keep q 1 and q 2 close to q 0.
We calculate the double scattering spectra of light from the Fourier power spectrum I(q) of the structure. I(q), which is measured directly by SAXS, tells the spatial vectors existing in the structure and their strength. Since the structure is isotropic, I(q) = I(q). For simplicity, we ignore the polarization of light. The probability of incident light with k i being scattered to k m is proportional to I(q 1 = k m - k i), and the probability of successive scattering from k m to k o is proportional to I(q 2 = k o - k m). The summation of I(q 1)I(q 2) over all possible directions of k m gives the probability of incident light being scattered to direction θm. In the calculation we sample 1000 directions of k m which are uniformly distributed over 4π solid angle . For comparison with the experimental data, we take into account light refraction at sample surface and sum the scattering spectra over a 5° collection angle.
Figure 4(b) shows the calculated double scattering spectra when the incident light is normal to the sample surface. As θ increases, the double scattering peak is broadened and split into two subpeaks. The subpeak at longer wavelength is stronger, and it red-shifts with increasing θ. The broadening of double scattering spectrum can be understood qualitatively by approximating I(q) as a delta function δ(q-q 0). When θm = 0, α 1 and α 2 can only be 90°. Only one wavelength is allowed and its value can be obtained by replacing θm in Eq. 1 by α 1 or α 2. Once θm > 0, α 1 and α 2 can take multiple values. The wavelength is no longer unique. The larger the θm, the wider the range of possible values for α 1 and α 2, and the broader the double scattering spectrum.
Experimentally the efficiency of our detector decreases rapidly when the wavelength of light is shorter than 340 nm. Thus we could not resolve the subpeak at shorter wavelength of the double scattering spectra. In Fig. 3 (b) we compare the q value for the subpeak at longer wavelength to that of the secondary peak in the measured scattering spectrum. They exhibit similar dependence on θ. Therefore, the secondary peak in the measured light scattering spectra is mostly from the subpeak at longer wavelength of the double scattering spectra in Fig. 4(b). As shown in Fig. 2(b), the measured secondary peak in the cross-polarized scattering spectrum is broadened and red-shifted compared to that in the co-polarized one. This may be due to triple scattering of light that can cause stronger depolarization than double scattering. A qualitative analysis of three successive single scattering events reveals that the triple scattering spectrum is broader than that of double scattering and its center is in between those of single scattering and double scattering.
Note that the secondary peak cannot be generated by the second-order scattering, which can be considered intuitively as two scattering events occur simultaneously at the same location and gives k o = k i+q 1+q 2. This process does not have such a constraint ∣k i+q 1∣ = k, which applies to the double scattering due to k conservation in each of two independent single scattering events. Consequently the second-order scattering peak is much broader in frequency. For example, in the backward direction k o = -k i, with the approximation I(q) ≈ δ (q-q 0), the value of k for the second-order scattering ranges from 0 and q 0. Hence, the second-order scattering spectrum spreads from λ = 280 nm to infinity, completely different from the measured secondary peak in Fig. 2(b).
Finally we compare our results to the previous ones. There have been a series of studies on butterfly scales with single-crystalline, poly-crystalline and amorphous photonic structures by measuring the spectrally and angularly resolved reflected and scattered light [17, 18, 19]. The optical reflectance spectra of the quasi-ordered structures can be explained by single scattering of light using the structural Fourier spectra . In this study, we demonstrate for the first time that double scattering in a short-range ordered nanostructures produces additional peaks in the light scattering spectrum that contribute significantly to structural colors. This phenomenon has no analog in either periodic structures or random structures. In a perfectly ordered structure, multiple diffraction or high-order diffraction does not create additional peaks, because a set of primary Bragg peaks would diffract into the same set of primary peaks. Vukusic et al. discovered new colors being generated by double reflection from the concaved multilayers of butterfly wings . The normally incident light, reflected from on one 45° inclined surface is directed across the concavity to the opposite orthogonal surface from where it returns parallel to the incident direction. This process strongly depends on the macroscopic geometry, and the direction of incident light. In contrast, the double scattering phenomenon we report here originates from the intrinsic nanostructure with short-range order. It is insensitive to the surface profile, and exists for any direction of incident light. Therefore, the creation of additional spectral feature by double scattering is a unique property of short-range ordered structures, that has been widely utilized by nature for structural coloration.
This work was supported with seed funding from the Yale NSF-MRSEC (DMR-0520495) and NSF grants to HC (EECS-0823345), ERD (CBET-0547294), and SGJM (DMR), and the W. R. Coe Fund, Yale University (to ROP). Feather specimens were provided by the Yale Peabody Museum of Natural History and the University of Kansas Natural History Museum and Biodiversity Research Center. SAXS data were collected at beam line 8-ID-I at the Advanced Photon Source at Argonne National Labs with the help of Drs. Alec Sandy and Suresh Narayanan, and supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357.
References and links
1. S. Kinoshita, S. Yoshioka, and J. Miyazaki, “Physics of structural colors,” Rep. Prog. Phys. , 71, 076401 (2008) [CrossRef]
2. H. Ghiradella, “Hairs, Bristles and Scales” in Microscopic Anatomy of Invertebrates, edited by M. Locke (Wiley-Liss, New York, 1998), Vol. 11A, 257–287
5. R. O. Prum, “Anatomy, Physics, and Evolution of Structural Colors” in Bird coloration, (Harvard University Press, Boston, 2006), Vol 1, 295–353
6. C. V. Raman, “The Origin of the Colours in the Plumage of Birds,” Proc. Ind. Acad. Sci. (A) 1, 1–7 (1934)
8. R. O. Prum, R. H. Torres, S. Williamson, and J. Dyck, “Coherent light scattering by blue feather barbs,” Nature 396, 28–29 (1998) [CrossRef]
9. D. Osorio and A. D. Ham, “Spectral reflectance and directional properties of structural coloration in bird plumage,” J. Exp. Biol. 205, 2017–2027 (2002) [PubMed]
10. M. D. Shawkey, V. Saranathan, H. Pálsdóttir, J. Crum, M. H. Ellisman, M. Auer, and R. O. Prum, “Electron tomography, three-dimensional Fourier analysis and colour prediction of a three-dimensional amorphous biophotonic nanostructure,” J. R. Soc. Interface 6, S213–S220 (2009) [PubMed]
11. E. R. Dufresne, H. Noh, V. Saranathan, S. G. J. Mochrie, H. Cao, and R. O. Prum,“Self-assembly of amorphous biophotonic nanostructures by phase separation,” Soft Matter 5, 1792 (2009) [CrossRef]
12. A. S. Barreira, G. Garcia, D. A. Lijtmaer, S. C. Lougheed, and P. L. Tubaro, “Blue males and green females:Sexual dichromatism in the Blue dacnis (Dacnis cayania) and the Swallow tanager (Tersina viridis),” Ornitologia Neotropical 19, 441–450 (2008)
13. J. Dyck, “Structural Colors,” Proc. Int. Orthnithol. Congr. 16, 426–437 (1976)
14. H. Noh, “Optical Properties of Nanostructures: from Random to Periodic,” Ph.D thesis (Northwestern Univ., Evanston, 2009)
15. Source code can be downloaded at http://local.wasp.uwa.edu.au/pbourke/geometry
16. P. Vukusic and R. J. Sambles, “Shedding light on butterfly wings,” Proc. SPIE 4438, 85–95 (2001) [CrossRef]
17. K. Kertész, Z. Bálint, Z. Vértesy, G. I. Márk, V. Lousse, J. P. Vigneron, M. Rassart, and L. P. Biró, “Gleaming and dull surface textures from photonic-crystal-type nanostructures in the butterfly Cyanophrys remus,” Phys. Rev. E 74, 021922 (2006) [CrossRef]
19. G. I. Márk, Z. Vértesy, K. Kertész, Z. Bálint, and L. P. Biró, “Order-disorder effects in structure and color relation of photonic-crystal-type nanostructures in butterfly wing scales,” Phys. Rev. E 80, 051903 (2009) [CrossRef]