Backscattering properties of small layered plates: a model for iridosomes

doi:10.1364/OE.20.025111

« Show journal navigation

Backscattering properties of small layered plates: a model for iridosomes

Meng Gao, Yu You, Ping Yang, and George W. Kattawar »View Author Affiliations

Optics Express, Vol. 20, Issue 22, pp. 25111-25120 (2012)
http://dx.doi.org/10.1364/OE.20.025111

View Full Text Article

Acrobat PDF (1958 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Numerical model & theoretical background
3. Numerical results and discussion
4. Conclusion
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (26)
Cited By
Figures (9)
Metrics

Abstract

The backscattering properties of small layered plates are studied for various size parameter values with respect to the plate thickness, plate aspect ratio, number of layers, incident direction, and polarization states of the incident light. The results are compared with the analytical results for semi-infinite plates. The phase functions and the corresponding backscattering efficiencies of the small plates are computed with the Discrete Dipole Approximation method. The angular width of the reflection peak is found to depend on both the size parameter and the aspect ratio. The criteria for using the reflectance of semi-infinite plates to approximate the backscattering efficiencies of finite size plates are quantified with respect to the number of layers, incident angle, and polarization state.

1. Introduction

Multiple-layered plates with finite sizes are abundant in animal skin and play an important role in producing color and polarized reflection [1

M. F. Land, “The physics and biology of animal reflectors,” Prog. Biophys. Mol. Biol. 24, 75–106 (1972). [CrossRef] [PubMed]

]. A type of small layered structure, called an iridosome, in the skin of cephalopods has been found to play a key role in both camouflage and communication [2

L. M. Mäthger, E. J. Denton, N. J. Marshall, and R. T. Hanlon, “Mechanisms and behavioural functions of structural coloration in cephalopods,” J. R. Soc. Interface 6(2), S149–S163 (2009). [PubMed]

]. Meanwhile, the analytical reflectance results of semi-infinite multiple-layered plates provide an elegant, accurate, and efficient way to study reflection properties [3

M. Born and E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffraction of light (Cambridge University Press, 1999).

–6

B. E. A. Saleh and M. C. Teich, Fundamentals of photonics (John Wiley & Sons, 2007).

]. This semi-infinite model has been widely used to approximate the reflectance of finite multiple-layered plates [1

M. F. Land, “The physics and biology of animal reflectors,” Prog. Biophys. Mol. Biol. 24, 75–106 (1972). [CrossRef] [PubMed]

, 7

A. L. Holt, A. M. Sweeney, S. Johnsen, and D. E. Morse, “A highly distributed Bragg stack with unique geometry provides effective camouflage for Loliginid squid eyes,” J. R. Soc. Interface 8(63), 1386–1399 (2011). [CrossRef] [PubMed]

–9

L. M. Mäthger and E. J. Denton, “Reflective properties of iridophores and fluorescent ‘eyespots’ in the loliginid squid Alloteuthis subulata and Loligo vulgaris,” J. Exp. Biol. 204(Pt 12), 2103–2118 (2001). [PubMed]

], but neither the backscattering properties of these plates nor the criteria for the validity of the semi-infinite approximation have been studied. We will numerically study the backscattering properties of small layered plates, based on the structure of a single iridosome, and compare the results with their semi-infinite plate counterparts. From these comparisons, we will provide the criteria for using the semi-infinite approximations.

Analysis based on semi-infinite plates can offer physical insights into the understanding of finite plates. This analysis approximately predicts the angle of specular reflection and the size parameter that leads to maximum constructive and destructive interference, and establishes the Brewster angle (i.e., the incident angle at which the reflected light is completely polarized) for a finite size structure. However, iridosomes of different species are substantially different in size [10

J. Arnold, R. Young, and M. King, “Ultrastructure of a cephalopod photophore. II. Iridophores as reflectors and transmitters,” Biol. Bull. 147(3), 522–534 (1974). [CrossRef]

–16

S. Mirow, “Skin color in the squids Loligo pealii and Loligo opalescens. II. Iridophores,” Z. Zellforsch. Mikrosk. Anat. 125(2), 176–190 (1972). [CrossRef] [PubMed]

], including both the thickness and the aspect ratio. Therefore, we need to evaluate the influence of the size parameter, aspect ratio, and other parameters, such as number of layers and polarization state, on the backscattering properties of a single iridosome.

To accurately account for all the scattering effects beyond the semi-infinite approximation, we use the Discrete Dipole Approximation (DDA) method [17

E. M. Purcell and C. R. Pennypacker, “Scattering and absorption of light by nonspherical dielectric grains,” Astrophys. J. 186, 705–714 (1973). [CrossRef]

, 18

M. A. Yurkin and A. G. Hoekstra, “The discrete-dipole-approximation code ADDA: Capabilities and known limitations,” J. Quant. Spectrosc. Radiat. Transf. 112(13), 2234–2247 (2011). [CrossRef]

] in this study. The results can be generalized to include other systems with similar structures, such as hexagonal ice crystal plates [19

P. Yang and K. N. Liou, “Light scattering by hexagonal ice crystals: solutions by a ray-by-ray integration algorithm,” J. Opt. Soc. Am. A 14(9), 2278–2289 (1997). [CrossRef]

]. Furthermore, the edge effect [20

L. Bi, P. Yang, and G. W. Kattawar, “Edge-effect contribution to the extinction of light by dielectric disks and cylindrical particles,” Appl. Opt. 49(24), 4641–4646 (2010). [CrossRef] [PubMed]

] by the small scatterers could be one reason the reflection properties of small plates deviate from those of semi-infinite plates.

2. Numerical model & theoretical background

An iridosome is modeled as a cylinder with a layered inner structure. Two examples are shown in Fig. 1 for both single-layer and 5-layer plates. Plate diameter, thickness, and the total cylinder height are respectively denoted by

D

d

, and

L

. The direction of the incident light is at an angle β relative to the symmetry axis of the iridosome and is shown in Fig. 1(a).

Fig. 1 (a) A single layer plate (b) A 5-layer plate.

The iridosome plate material is found to be a kind of protein, called reflectin, and its refractive index is measured to be

1.591 \pm 0.002

with a negligible dispersion in the visual spectrum [21

T. H. Chiou, L. M. Mäthger, R. T. Hanlon, and T. W. Cronin, “Spectral and spatial properties of polarized light reflections from the arms of squid (Loligo pealeii) and cuttlefish (Sepia officinalis L.),” J. Exp. Biol. 210(20), 3624–3635 (2007). [CrossRef] [PubMed]

, 22

R. M. Kramer, W. J. Crookes-Goodson, and R. R. Naik, “The self-organizing properties of squid reflectin protein,” Nat. Mater. 6(7), 533–538 (2007). [CrossRef] [PubMed]

]. The material of the ambient tissue and that between plates are mainly cytoplasm, whose refractive index is close to that of water (

n_{w} = 1.33

). Thus, the relative refractive index of the plate (

n_{r}

) is 1.2 referenced to the cytoplasm. In our model, we choose the same relative refractive index for iridosomes. To further simplify the modeled structure, we consider ideal multiple-layered plates in which the optical length of the plates and the spacing between plates are the same, and therefore the spacing thickness is

n_{r} d

The size parameter is defined as

x = 2 π d / λ

, where

λ

is the wavelength in the medium. The plate thickness d of iridosomes can range from 50nm to 200nm [10

J. Arnold, R. Young, and M. King, “Ultrastructure of a cephalopod photophore. II. Iridophores as reflectors and transmitters,” Biol. Bull. 147(3), 522–534 (1974). [CrossRef]

–16

S. Mirow, “Skin color in the squids Loligo pealii and Loligo opalescens. II. Iridophores,” Z. Zellforsch. Mikrosk. Anat. 125(2), 176–190 (1972). [CrossRef] [PubMed]

], and, in the visual spectrum, the incident wavelength

λ_{0}

in air is from 0.4 to 0.75 μm. Therefore, the size parameter of the plate thickness in the medium is

x = 2 π n_{w} d / λ_{0}

, and can vary approximately from 0.5 to 5. This size parameter range is used in our simulations in steps of 0.1.

The aspect ratio is defined as

a = D / d

for one layer of the plate. For example, both single-layer and 5-layer plates in Figs. 1(a) and (b) have the same aspect ratio of 10. A range of aspect ratios from 1 to 20 is chosen for our study. Rather than discussing the aspect ratio variations among certain cephalopod species, we focus on the aspect ratio’s influence on the optical properties, from which a particular type of iridisome can be evaluated later on.

The phase function

P (Ω)

describes the angular distribution of scattered light for an unpolarized light source, where

Ω = (θ_{s}, ϕ)

is the solid angle for the scattering direction.

ϕ

is the azimuthal angle, and

θ_{s}

is the scattering angle between the scattering direction and the incident direction. As shown in Fig. 1(a), the upper solid angle hemisphere

Ω^{+}

is defined as the backscattering region. Sine the scattered light has an angular distribution described by

P (Ω)

, and as we will discuss later, even for iridosomes with large aspect ratios, there is still a finite angular width associated with the reflection peak in the phase function. In order to compare the reflectance of irradiance for the semi-infinite plates, we will integrate the phase function over the whole backscattering region

Ω^{+}

. Therefore, we can define the backscattering efficiency for unpolarized incident light as

Q_{b} = \frac{1}{k^{2} C_{g}} \iint_{Ω^{+}} P (Ω) d Ω,

(1)

where

k = 2 π / λ

is the wavenumber, and

C_{g} = π D^{2} \cos β / 4 + D L \sin β

is the projected geometric area of the cylinder. Note that the projected geometric area is perpendicular to the incident direction. Two sets of examples of the projected areas are shown in Fig. 2 for both single-layer and 5-layer plates. Q_b indicates the irradiance reflectance for a finite iridosome. When the aspect ratio approaches infinity, Q_b becomes equal to the reflectance of the semi-infinite plate. Therefore, we use the backscattering efficiency as an indicator to show the reflection deviation of the small plates from that of the semi-infinite plates.

Fig. 2 The projected geometric areas (

C_{g}

) which are perpendicular to the incident direction are shown for (a) a single-layer plate and (b) a 5-layer plate. Incident angle

β

is from 0° to 90°.

Besides representing the reflectance of a single iridosome, Q_b can also be used to calculate the diffuse irradiance reflectance

R_{d i f f}

for a plane parallel system with a collection of iridosomes. We assume that

n_{d}

is the number density of the iridosomes and

Z

is the thickness of the system. For an optically thin system, the multiple scattering between iridosomes is ignored and only the scattering by each single iridosome dominates. Considering an area of A and thickness Z in this system, the total number of iridosomes would be

N = n_{d} A Z

. The total reflectance of the whole system would be the ratio of the total backscattering cross section and the total area A. Therefore, we have

R_{d i f f} = N C_{b} / A = n_{d} C_{b} Z

, where

C_{b} = Q_{b} C_{g}

is the backscattering cross section of a single iridosome. Iridosomes can have arbitrary orientations and C_b needs to be averaged according to the orientation distributions of these iridosomes. For an optically thick system with random oriented iridosomes,

R_{d i f f}

can still be estimated using the randomly averaged C_b. For such a system, there is significant multiple scattering which tends to drive the reflected radiance toward an isotropic distribution, and we can use a one-dimensional model to calculate the irradiance reflectance as proposed by Kattawar et al. [23

G. W. Kattawar and G. N. Plass, “Interior radiances in optically deep absorbing media. 1.Exact solutions for one-dimensional model,” J. Quant. Spectrosc. Radiat. Transf. 13(11), 1065–1080 (1973). [CrossRef]

], where the diffuse reflectance can be obtained as

R_{d i f f} = n_{d} C_{b} Z / (1 + n_{d} C_{b} Z)

. Therefore, the backscattering efficiency of one scatterer together with its geometric cross section can be related to the diffuse reflection of a bulk system. The diffuse reflectance is important in determining the optical appearance of the system, such as the color and brightness of the cephalopod skin. However, if the angular distribution of the diffuse reflection is considered, the radiative transfer equation must be solved using the full phase function [24

C. Bohren and E. Clothiaux, Fundamentals of atmospheric radiation (Wiley-VCH, 2006).

To calculate the reflectance of a single-layer semi-infinite plate (R₁) (Fig. 3 ), Airy’s formula is used for different size parameters, polarization states, and both normal and oblique incidences [3

M. Born and E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffraction of light (Cambridge University Press, 1999).

R_{1 s, p} = \frac{4 r_{s, p}^{2} \sin^{2} α}{{(1 - r_{s, p}^{2})}^{2} + 4 r_{s, p}^{2} \sin^{2} α},

(2)

r_{s} = \frac{\cos β - n_{r} \cos θ}{\cos β + n_{r} \cos θ},

(3)

r_{p} = \frac{n_{r} \cos β - \cos θ}{n_{r} \cos β + \cos θ},

(4)

where

α = 2 π n_{r} d \cos β / λ

is a phase factor and related to the size parameter by

α = n_{r} x \cos β

. As shown in Eqs. (3), (4),

r_{s}

and

r_{p}

are the Fresnel reflectance formulas for the perpendicular and parallel polarized components relative to the principal plane defined by the direction of the incident and reflected light, where

β

is the incident angle and

θ

is the refracted angle. The reflectance is a periodic function of the size parameter

x

. Our analytical results for an arbitrary number of layers are derived by using the transfer matrix method [4

E. Hecht, Optics (Addison-Wesley, 2002).

, 6

B. E. A. Saleh and M. C. Teich, Fundamentals of photonics (John Wiley & Sons, 2007).

], and the resulting reflectance remains a periodic function of the size parameter.

Fig. 3 Illustration of the interference of the light reflected from the top and bottom surfaces of a semi-infinite plate with thickness d and relative refractive index

n_{r}

in the principal plane.

Additionally, the maximum constructive and destructive interferences determine the maximum and minimum reflection of a semi-infinite plate. As shown in Fig. 3, the difference in the optical lengths of the two reflection paths is

Δ l = 2 n_{r} s_{1} - s_{0} = 2 n_{r} d \cos θ

, where

θ

is the refracted angle. The phase factor α in Eq. (2) can be related to the path difference by

α = π Δ l

under normal incidence. Due to the

λ / 2

phase shift at the interface, the condition for maximum constructive interference is

Δ l = (m + 1 / 2) λ

, where

m

is an integer. Thus, to have maximum constructive interference under normal incidence, we have

α = π / 2

and

3 π / 2

, which corresponds to

x = 1.3

and 3.9; to have maximum destructive interference, we have

α = π

, which corresponds to

x = 2.6

. However, for a finite size system, the reflection maxima and minima deviate from the prediction of these size parameters, and to compare the difference we indicate

x = 1.3

, 2.6, and 3.9 in the reflection spectrum graphs.

3. Numerical results and discussion

We used the ADDA code developed by Yurkin and Hoekstra [18

M. A. Yurkin and A. G. Hoekstra, “The discrete-dipole-approximation code ADDA: Capabilities and known limitations,” J. Quant. Spectrosc. Radiat. Transf. 112(13), 2234–2247 (2011). [CrossRef]

] for the DDA simulation to calculate the scattering properties of the iridosome. In this simulation, 10 dipoles are used along the thickness of every single plate, and the required relative residual norm is specified as 10⁻⁵ for the involved quasi minimal residual solver. An angular resolution of 0.25°

\times

2° for the zenith and azimuthal angle is used to calculate the backscattering efficiency. The phase functions are calculated with various size parameters, aspect ratios, incident directions, and polarization states. We discuss the angular distribution of the reflection peak in the phase functions and compare the backscattering efficiencies of the iridosome with the reflectance of the semi-infinite plates.

For a size parameter x = 1.3, the phase functions are shown in Fig. 4 with different aspect ratios under normal incidence for a single-layer plate. For convenience in comparing backscattering properties, each phase function is normalized by

k^{2} C_{g}

as in Eq. (1). Because the light is at normal incidence on the plate, the phase function is azimuthally symmetric. For example, as shown in Fig. 4, the phase function is symmetric around

θ_{s}

= 180°, which corresponds to the direct backscattering direction. The full width at half maximum (FWHM) of the reflection peaks is denoted by

Δ θ

and provides a way to quantify the angular spreading of the reflection. When the aspect ratio a increases, the reflection peak around

θ_{s}

= 180° becomes narrower and

Δ θ

decreases. The FWHM of the reflection peak is shown by the horizontal red solid bar in Fig. 4. For the aspect ratio from 4 to 40,

Δ θ

is found to decrease from 73.5° to 7.5° with an angular resolution of 0.5°. For an aspect ratio of 1 or 2, no angle corresponds to the phase function with a value half of that at 180°, and, therefore, no value for

Δ θ

exists; the phase function has less variations compared with the one with large aspect ratio particles.

Fig. 4 Phase function

P (θ_{s})

with aspect ratio a ranging from 1 to 40 for a single-layer plate under normal incidence and size parameter x = 1.3. The FWHM of the reflection peaks (

Δ θ

) is indicated as a solid horizontal red bar. The aspect ratio, a, is shown close to its corresponding curve.

For size parameters in the range from 0.5 to 5, the values of

Δ θ

are summarized in Fig. 5 with the same 0.5° angular resolution for the single-layer plate as just discussed.

Δ θ

decreases as the aspect ratio increases throughout all the considered size parameters. Note that

Δ θ

is sensitive to the size parameter around x = 2.6, which corresponds to maximum destructive interference for the semi-infinite plate. A jump of the value can be observed for a = 4 around x = 2.6. We further studied

Δ θ

with a finer size parameter resolution, and observed that a continuous peak existed for the value of

Δ θ

around x = 2.6. For aspect ratios a = 1 and 2 and size parameters approximately less than 2.6, no value exists for

Δ θ

. For size parameters larger than 2.6 and a = 1,

Δ θ

can be found between 44° to 96°, and for a = 2,

Δ θ

decreases to between 36° and 46°. Furthermore, when a = 4,

Δ θ

is between 130° and 18.5° for size parameters from 0.5 to 5. When a = 40,

Δ θ

further decreases to between 18.5° and 2° for the same range of size parameters. For the plates with a thickness of 5 and a = 40, the reflection peak is at least 2° in width. The spreading becomes less prominent when the aspect ratio is large. However, for the semi-infinite plates, the reflection is confined to a single direction determined by Snell’s law. The angular broadening of the reflection can be explained by the ray-spreading effect [25

P. Yang and K. N. Liou, “Geometric-optics-integral-equation method for light scattering by nonspherical ice crystals,” Appl. Opt. 35(33), 6568–6584 (1996). [CrossRef] [PubMed]

]. When the angular distribution is considered either for a phase function or the diffuse reflection of a random system, the influence of the spreading of the reflection peak needs to be evaluated.

Fig. 5 The FWHM of the reflection peak (

Δ θ

) versus the plate thickness size parameter for a single-layer plate. The aspect ratio is a. Vertical dashed lines are for size parameters x = 1.3, 2.6, and 3.9.

The backscattering efficiency Q_b (Eq. (1)) is sensitive to the aspect ratio and the number of layers of an iridosome. As the aspect ratio increases, Q_b will undergo a metamorphosis from a small particle to a semi-infinite layered plate. In Fig. 6 , Q_b is shown for various numbers of layers and aspect ratios at x = 1.3. When the radius increases, Q_b approaches the reflection of the semi-infinite plates (red asterisks in the graph) for each number of layers. However, the dependency of Q_b on the aspect ratio is not monotonic. Q_b first increases for the aspect ratio from a = 1 to 6, then decreases from a = 6 to 20. This can be understood by looking at a particular size parameter, x = 1.3, in the spectrum of Q_b shown in Fig. 7(b) ; with the decreasing of the aspect ratio, the size parameter for the maximum Q_b will shift toward larger size parameters. When a = 1 or 2 and the number of layers is larger than 1, Q_b is much smaller than the reflection of semi-infinite plates. When a = 4, Q_b will decrease as the number of layers increase when the number of layers is larger than 7. In general for a>4, as shown in Fig. 6, adding more layers of plates into the scattering will increase the backscattering efficiencies as predicted by the semi-infinite plates. The size effect plays a significant role for small scatterers in determining the backscattering efficiency. Meanwhile, Q_b can be larger than 1 due to the interference effects of the scattered waves similar to the extinction efficiency [26

H. C. van de Hulst, Light scattering by small particles (Dover Publications, 1981).

Fig. 6 Backscattering efficiency Q_b versus number of layers for a size parameter of the plate thickness

x = 1.3

with different aspect ratios and under normal incidence. The reflection for semi-infinite plates is indicated by asterisks.

Fig. 7 Backscattering efficiency Q_b versus size parameter

x = 2 π d / λ

of the plate thickness at different aspect ratios under normal incidence. The reflection for semi-infinite plates is indicated by solid red lines. (a) Results for a single-layer plate. Each curve is successively moved upward by 0.05. Maximum reflection of the semi-infinite plates is 0.033 at x = 1.3 and 3.9. (b) Results for 5-layer plate. Each curve is successively moved upward by 1.0. Maximum reflection of the semi-infinite plates is 0.52 at x = 1.3 and 3.9. The step for the size parameter in the calculation is 0.1. Vertical dashed lines are for the interference maxima with size parameter x = 1.3, 2.6, and 3.9

For normal incidence, the backscattering efficiency Q_b versus the size parameter for a single-layer and a 5-layer plate is compared with the semi-infinite results in Fig. 7. With large aspect ratios, such as a = 20, the backscattering efficiency can be approximated with great accuracy by using the results from the semi-infinite plates. When the aspect ratio decreases, two maximum reflection peaks remain observed, but the positions shift to shorter wavelengths and the values generally differ more from the semi-infinite results. For example, when a = 1, as shown in Fig. 7(a) and (b), one of the Q_b maxima around x = 2.6 becomes the minimum for the semi-infinite plate. Since the radius of the plates is very small, unlike the semi-infinite plates, the interference effects between the top and bottom interface of the plates become unimportant in determining Q_b.

To quantify the error of using the semi-infinite results to approximate the backscattering efficiency (Q_b), we defined a quantity Q_error by comparing the value of the maxima between the finite and semi-infinite plates: Q_error = |(Q_b,max-R_max)/ (Q_b,max + R_max)|, where Q_b,max is the maximum backscattering efficiency for the finite plates and R_max is the maximum reflectance for the semi-infinite plates. A smaller Q_error means less difference between the maxima of these two kinds of plates, and, as observed from Fig. 7, when the maxima are close, the wavelengths for the locations of those maxima are also close to each other. We choose a Q_error of less than 10% as an acceptable criterion to use the semi-infinite approximation. Therefore, from Fig. 7, for single-layer plates, the aspect ratio must be

a \geq 4

, and for 5-layer plates, the aspect ratio must be

a \geq 8

. A larger layer number requires a larger aspect ratio.

With an aspect ratio a = 10, for a single-layer plate, R_max = 0.033 at both x = 1.3 and 3.9, Q_b,max = 0.035 at x = 1.4, and Q_error = 4%; for a 5-layer plate, R_max = 0.52 at both x = 1.3 and 3.9, Q_b,max = 0.60 at x = 4.0, and Q_error = 7%. At this aspect ratio, both the single-layer plate and the 5-layer plate have Q_error <10% and, therefore, to further study the dependency of backscattering efficiency on the incident angles we use a = 10 for both the single and 5-layer plates.

The angular dependence of the backscattering efficiency (Q_b) can be used to study the response of different incident polarizations on the small plates backscattering. The results are shown in Figs. 8 and 9 for angle intervals of 10° and an aspect ratio of 10 for both perpendicular and parallel polarization states. To have Q_error ≤ 10% for the single-layer plate, we need an incident angle β≤ 30° for a parallel polarized component, and β≤ 50° for a perpendicular component. However, for the 5-layer plate to have the same Q_error≤10%, we need β≤10° for both polarizations due to the large fraction of edge compared to the total projected area. Therefore, the perpendicular component provides a better approximation compared with the parallel polarized component for both the single-layer and 5-layer results. The reflectance of semi-infinite plates approaches unity for all size parameters as the incident angle approaches 90°. The flat line in the top of Figs. 8 and 9 indicates a 100% reflection for β = 90° calculated by using Eq. (2). Since the light is incident directly on the edges of the plates, the backscattering efficiency is completely different from the reflection. The large value of Q_b is mainly due to the small projected area as shown in Fig. 2. We cannot use the reflection of semi-infinite plates to approximate Q_b for large incident angles.

Fig. 8 Same as Fig. 7 but for a single-layer plate with an aspect ratio 10 and varying incident angles

β

. Incident light is polarized in (a) parallel direction and (b) perpendicular direction relative to the principal plane. Each curve is successively moved upward by 0.1 except for

β

= 70°, 80°, and 90° which are plotted in the upper panels. The solid red lines indicate the reflection for semi-infinite plates for

β

= 10° to 90° from bottom to top.

Fig. 9 Same as Fig. 8 but for a five-layer plate. Each curve is successively moved upward by 1.0.

The semi-infinite plates will respond differently to different incident polarization states; the reflectance of the parallel-polarized incident light will vanish at the Brewster angle, leaving the reflected light totally polarized in the direction perpendicular to the principal plane. As shown in Figs. 8 and 9, we found similar behaviors for both the single and 5-layer finite size plates. We compared the maximum reflectance for each incident angle, and found the minimum value for single-layer plates to be 0.013 and for 5-layer plates to be 0.06, and both to occur at β = 40°. The angle is smaller than the Brewster angle of 50.2° for n_r = 1.2. Moreover, these two values are 37% and 10% of their corresponding maximum Q_b at normal incidence, while for the semi-infinite plates the value of the minimum reflectance of the parallel, polarized light will be zero at the Brewster angle. This constitutes another prominent difference between finite and semi-infinite plates.

4. Conclusion

We studied the backscattering properties of small layered plates as a model for iridosomes, including both the backscattering efficiency and the shape of the reflection peak in the phase function. The width of the reflection peak is shown to depend on both the aspect ratio and the size parameter. The criteria for using the reflection of the semi-infinite plate to approximate the backscattering efficiency of small plates are discussed with respect to the aspect ratio, number of layers, and incident angle. For normal incidence and with an accuracy Q_error ≤ 10%, the semi-infinite plate assumption can work for a single-layer plate when the aspect ratio satisfies

a \geq 4

and for 5-layer plates when

a \geq 8

. For plates with small aspect ratios and large incident angles, the scattering properties must be calculated beyond the semi-infinite approximation. Furthermore, the backscattering efficiency can be used to estimate the diffuse reflectance of a random system, and therefore, determine its optical appearance.

Acknowledgments

This research was partially supported by the ONR MURI program N00014-09-1-1054. Drs. George Kattawar and Ping Yang also acknowledge support from the National Science Foundation (OCE-1130906) and (ATMO-0803779) respectively. We also thank the two anonymous reviewers for their contributions on improving the readability of the paper.

References and links

1.	M. F. Land, “The physics and biology of animal reflectors,” Prog. Biophys. Mol. Biol. 24, 75–106 (1972). [CrossRef] [PubMed]
2.	L. M. Mäthger, E. J. Denton, N. J. Marshall, and R. T. Hanlon, “Mechanisms and behavioural functions of structural coloration in cephalopods,” J. R. Soc. Interface 6(2), S149–S163 (2009). [PubMed]
3.	M. Born and E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffraction of light (Cambridge University Press, 1999).
4.	E. Hecht, Optics (Addison-Wesley, 2002).
5.	A. Huxley, “A theoretical treatment of the reflexion of light by multilayer structures,” J. Exp. Biol. 48, 227–245 (1968).
6.	B. E. A. Saleh and M. C. Teich, Fundamentals of photonics (John Wiley & Sons, 2007).
7.	A. L. Holt, A. M. Sweeney, S. Johnsen, and D. E. Morse, “A highly distributed Bragg stack with unique geometry provides effective camouflage for Loliginid squid eyes,” J. R. Soc. Interface 8(63), 1386–1399 (2011). [CrossRef] [PubMed]
8.	R. L. Sutherland, L. M. Mäthger, R. T. Hanlon, A. M. Urbas, and M. O. Stone, “Cephalopod coloration model. I. Squid chromatophores and iridophores,” J. Opt. Soc. Am. A 25(3), 588–599 (2008). [CrossRef] [PubMed]
9.	L. M. Mäthger and E. J. Denton, “Reflective properties of iridophores and fluorescent ‘eyespots’ in the loliginid squid Alloteuthis subulata and Loligo vulgaris,” J. Exp. Biol. 204(Pt 12), 2103–2118 (2001). [PubMed]
10.	J. Arnold, R. Young, and M. King, “Ultrastructure of a cephalopod photophore. II. Iridophores as reflectors and transmitters,” Biol. Bull. 147(3), 522–534 (1974). [CrossRef]
11.	S. L. Brocco and R. A. Cloney, “Reflector cells in the skin of Octopus dofleini,” Cell Tissue Res. 205(2), 167–186 (1980). [CrossRef] [PubMed]
12.	R. Cloney and S. Brocco, “Chromatophore organs, reflector cells, iridocytes and leucophores in cephalopods,” Integr. Comp. Biol. 23(3), 581–592 (1983). [CrossRef]
13.	K. M. Cooper, R. T. Hanlon, and B. U. Budelmann, “Physiological color change in squid iridophores. II. Ultrastructural mechanisms in Lolliguncula brevis,” Cell Tissue Res. 259(1), 15–24 (1990). [CrossRef] [PubMed]
14.	D. Froesch and J. B. Messenger, “On leucophores and the chromatic unit of Octopus vulgaris,” J. Zool. 186(2), 163–173 (1978). [CrossRef]
15.	L. M. Mäthger, T. F. Collins, and P. A. Lima, “The role of muscarinic receptors and intracellular Ca2+ in the spectral reflectivity changes of squid iridophores,” J. Exp. Biol. 207(11), 1759–1769 (2004). [CrossRef] [PubMed]
16.	S. Mirow, “Skin color in the squids Loligo pealii and Loligo opalescens. II. Iridophores,” Z. Zellforsch. Mikrosk. Anat. 125(2), 176–190 (1972). [CrossRef] [PubMed]
17.	E. M. Purcell and C. R. Pennypacker, “Scattering and absorption of light by nonspherical dielectric grains,” Astrophys. J. 186, 705–714 (1973). [CrossRef]
18.	M. A. Yurkin and A. G. Hoekstra, “The discrete-dipole-approximation code ADDA: Capabilities and known limitations,” J. Quant. Spectrosc. Radiat. Transf. 112(13), 2234–2247 (2011). [CrossRef]
19.	P. Yang and K. N. Liou, “Light scattering by hexagonal ice crystals: solutions by a ray-by-ray integration algorithm,” J. Opt. Soc. Am. A 14(9), 2278–2289 (1997). [CrossRef]
20.	L. Bi, P. Yang, and G. W. Kattawar, “Edge-effect contribution to the extinction of light by dielectric disks and cylindrical particles,” Appl. Opt. 49(24), 4641–4646 (2010). [CrossRef] [PubMed]
21.	T. H. Chiou, L. M. Mäthger, R. T. Hanlon, and T. W. Cronin, “Spectral and spatial properties of polarized light reflections from the arms of squid (Loligo pealeii) and cuttlefish (Sepia officinalis L.),” J. Exp. Biol. 210(20), 3624–3635 (2007). [CrossRef] [PubMed]
22.	R. M. Kramer, W. J. Crookes-Goodson, and R. R. Naik, “The self-organizing properties of squid reflectin protein,” Nat. Mater. 6(7), 533–538 (2007). [CrossRef] [PubMed]
23.	G. W. Kattawar and G. N. Plass, “Interior radiances in optically deep absorbing media. 1.Exact solutions for one-dimensional model,” J. Quant. Spectrosc. Radiat. Transf. 13(11), 1065–1080 (1973). [CrossRef]
24.	C. Bohren and E. Clothiaux, Fundamentals of atmospheric radiation (Wiley-VCH, 2006).
25.	P. Yang and K. N. Liou, “Geometric-optics-integral-equation method for light scattering by nonspherical ice crystals,” Appl. Opt. 35(33), 6568–6584 (1996). [CrossRef] [PubMed]
26.	H. C. van de Hulst, Light scattering by small particles (Dover Publications, 1981).

OCIS Codes
(290.1350) Scattering : Backscattering
(310.6860) Thin films : Thin films, optical properties
(160.1435) Materials : Biomaterials
(290.5855) Scattering : Scattering, polarization

ToC Category:
Scattering

History
Original Manuscript: August 24, 2012
Revised Manuscript: October 10, 2012
Manuscript Accepted: October 10, 2012
Published: October 18, 2012

Citation
Meng Gao, Yu You, Ping Yang, and George W. Kattawar, "Backscattering properties of small layered plates: a model for iridosomes," Opt. Express 20, 25111-25120 (2012)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-22-25111

Sort: Author | Year | Journal | Reset

References

M. F. Land, “The physics and biology of animal reflectors,” Prog. Biophys. Mol. Biol.24, 75–106 (1972). [CrossRef] [PubMed]
L. M. Mäthger, E. J. Denton, N. J. Marshall, and R. T. Hanlon, “Mechanisms and behavioural functions of structural coloration in cephalopods,” J. R. Soc. Interface6(2), S149–S163 (2009). [PubMed]
M. Born and E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffraction of light (Cambridge University Press, 1999).
E. Hecht, Optics (Addison-Wesley, 2002).
A. Huxley, “A theoretical treatment of the reflexion of light by multilayer structures,” J. Exp. Biol.48, 227–245 (1968).
B. E. A. Saleh and M. C. Teich, Fundamentals of photonics (John Wiley & Sons, 2007).
A. L. Holt, A. M. Sweeney, S. Johnsen, and D. E. Morse, “A highly distributed Bragg stack with unique geometry provides effective camouflage for Loliginid squid eyes,” J. R. Soc. Interface8(63), 1386–1399 (2011). [CrossRef] [PubMed]
R. L. Sutherland, L. M. Mäthger, R. T. Hanlon, A. M. Urbas, and M. O. Stone, “Cephalopod coloration model. I. Squid chromatophores and iridophores,” J. Opt. Soc. Am. A25(3), 588–599 (2008). [CrossRef] [PubMed]
L. M. Mäthger and E. J. Denton, “Reflective properties of iridophores and fluorescent ‘eyespots’ in the loliginid squid Alloteuthis subulata and Loligo vulgaris,” J. Exp. Biol.204(Pt 12), 2103–2118 (2001). [PubMed]
J. Arnold, R. Young, and M. King, “Ultrastructure of a cephalopod photophore. II. Iridophores as reflectors and transmitters,” Biol. Bull.147(3), 522–534 (1974). [CrossRef]
S. L. Brocco and R. A. Cloney, “Reflector cells in the skin of Octopus dofleini,” Cell Tissue Res.205(2), 167–186 (1980). [CrossRef] [PubMed]
R. Cloney and S. Brocco, “Chromatophore organs, reflector cells, iridocytes and leucophores in cephalopods,” Integr. Comp. Biol.23(3), 581–592 (1983). [CrossRef]
K. M. Cooper, R. T. Hanlon, and B. U. Budelmann, “Physiological color change in squid iridophores. II. Ultrastructural mechanisms in Lolliguncula brevis,” Cell Tissue Res.259(1), 15–24 (1990). [CrossRef] [PubMed]
D. Froesch and J. B. Messenger, “On leucophores and the chromatic unit of Octopus vulgaris,” J. Zool.186(2), 163–173 (1978). [CrossRef]
L. M. Mäthger, T. F. Collins, and P. A. Lima, “The role of muscarinic receptors and intracellular Ca2+ in the spectral reflectivity changes of squid iridophores,” J. Exp. Biol.207(11), 1759–1769 (2004). [CrossRef] [PubMed]
S. Mirow, “Skin color in the squids Loligo pealii and Loligo opalescens. II. Iridophores,” Z. Zellforsch. Mikrosk. Anat.125(2), 176–190 (1972). [CrossRef] [PubMed]
E. M. Purcell and C. R. Pennypacker, “Scattering and absorption of light by nonspherical dielectric grains,” Astrophys. J.186, 705–714 (1973). [CrossRef]
M. A. Yurkin and A. G. Hoekstra, “The discrete-dipole-approximation code ADDA: Capabilities and known limitations,” J. Quant. Spectrosc. Radiat. Transf.112(13), 2234–2247 (2011). [CrossRef]
P. Yang and K. N. Liou, “Light scattering by hexagonal ice crystals: solutions by a ray-by-ray integration algorithm,” J. Opt. Soc. Am. A14(9), 2278–2289 (1997). [CrossRef]
L. Bi, P. Yang, and G. W. Kattawar, “Edge-effect contribution to the extinction of light by dielectric disks and cylindrical particles,” Appl. Opt.49(24), 4641–4646 (2010). [CrossRef] [PubMed]
T. H. Chiou, L. M. Mäthger, R. T. Hanlon, and T. W. Cronin, “Spectral and spatial properties of polarized light reflections from the arms of squid (Loligo pealeii) and cuttlefish (Sepia officinalis L.),” J. Exp. Biol.210(20), 3624–3635 (2007). [CrossRef] [PubMed]
R. M. Kramer, W. J. Crookes-Goodson, and R. R. Naik, “The self-organizing properties of squid reflectin protein,” Nat. Mater.6(7), 533–538 (2007). [CrossRef] [PubMed]
G. W. Kattawar and G. N. Plass, “Interior radiances in optically deep absorbing media. 1.Exact solutions for one-dimensional model,” J. Quant. Spectrosc. Radiat. Transf.13(11), 1065–1080 (1973). [CrossRef]
C. Bohren and E. Clothiaux, Fundamentals of atmospheric radiation (Wiley-VCH, 2006).
P. Yang and K. N. Liou, “Geometric-optics-integral-equation method for light scattering by nonspherical ice crystals,” Appl. Opt.35(33), 6568–6584 (1996). [CrossRef] [PubMed]
H. C. van de Hulst, Light scattering by small particles (Dover Publications, 1981).

Appl. Opt.

Astrophys. J.

E. M. Purcell and C. R. Pennypacker, “Scattering and absorption of light by nonspherical dielectric grains,” Astrophys. J.186, 705–714 (1973). [CrossRef]

Biol. Bull.

J. Arnold, R. Young, and M. King, “Ultrastructure of a cephalopod photophore. II. Iridophores as reflectors and transmitters,” Biol. Bull.147(3), 522–534 (1974). [CrossRef]

Cell Tissue Res.

S. L. Brocco and R. A. Cloney, “Reflector cells in the skin of Octopus dofleini,” Cell Tissue Res.205(2), 167–186 (1980). [CrossRef] [PubMed]
K. M. Cooper, R. T. Hanlon, and B. U. Budelmann, “Physiological color change in squid iridophores. II. Ultrastructural mechanisms in Lolliguncula brevis,” Cell Tissue Res.259(1), 15–24 (1990). [CrossRef] [PubMed]

Integr. Comp. Biol.

R. Cloney and S. Brocco, “Chromatophore organs, reflector cells, iridocytes and leucophores in cephalopods,” Integr. Comp. Biol.23(3), 581–592 (1983). [CrossRef]

J. Exp. Biol.

L. M. Mäthger, T. F. Collins, and P. A. Lima, “The role of muscarinic receptors and intracellular Ca2+ in the spectral reflectivity changes of squid iridophores,” J. Exp. Biol.207(11), 1759–1769 (2004). [CrossRef] [PubMed]
A. Huxley, “A theoretical treatment of the reflexion of light by multilayer structures,” J. Exp. Biol.48, 227–245 (1968).
L. M. Mäthger and E. J. Denton, “Reflective properties of iridophores and fluorescent ‘eyespots’ in the loliginid squid Alloteuthis subulata and Loligo vulgaris,” J. Exp. Biol.204(Pt 12), 2103–2118 (2001). [PubMed]
T. H. Chiou, L. M. Mäthger, R. T. Hanlon, and T. W. Cronin, “Spectral and spatial properties of polarized light reflections from the arms of squid (Loligo pealeii) and cuttlefish (Sepia officinalis L.),” J. Exp. Biol.210(20), 3624–3635 (2007). [CrossRef] [PubMed]

J. Opt. Soc. Am. A

J. Quant. Spectrosc. Radiat. Transf.

M. A. Yurkin and A. G. Hoekstra, “The discrete-dipole-approximation code ADDA: Capabilities and known limitations,” J. Quant. Spectrosc. Radiat. Transf.112(13), 2234–2247 (2011). [CrossRef]
G. W. Kattawar and G. N. Plass, “Interior radiances in optically deep absorbing media. 1.Exact solutions for one-dimensional model,” J. Quant. Spectrosc. Radiat. Transf.13(11), 1065–1080 (1973). [CrossRef]

J. R. Soc. Interface

A. L. Holt, A. M. Sweeney, S. Johnsen, and D. E. Morse, “A highly distributed Bragg stack with unique geometry provides effective camouflage for Loliginid squid eyes,” J. R. Soc. Interface8(63), 1386–1399 (2011). [CrossRef] [PubMed]
L. M. Mäthger, E. J. Denton, N. J. Marshall, and R. T. Hanlon, “Mechanisms and behavioural functions of structural coloration in cephalopods,” J. R. Soc. Interface6(2), S149–S163 (2009). [PubMed]

J. Zool.

D. Froesch and J. B. Messenger, “On leucophores and the chromatic unit of Octopus vulgaris,” J. Zool.186(2), 163–173 (1978). [CrossRef]

Nat. Mater.

R. M. Kramer, W. J. Crookes-Goodson, and R. R. Naik, “The self-organizing properties of squid reflectin protein,” Nat. Mater.6(7), 533–538 (2007). [CrossRef] [PubMed]

Prog. Biophys. Mol. Biol.

M. F. Land, “The physics and biology of animal reflectors,” Prog. Biophys. Mol. Biol.24, 75–106 (1972). [CrossRef] [PubMed]

Z. Zellforsch. Mikrosk. Anat.

S. Mirow, “Skin color in the squids Loligo pealii and Loligo opalescens. II. Iridophores,” Z. Zellforsch. Mikrosk. Anat.125(2), 176–190 (1972). [CrossRef] [PubMed]

Other

M. Born and E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffraction of light (Cambridge University Press, 1999).
E. Hecht, Optics (Addison-Wesley, 2002).
B. E. A. Saleh and M. C. Teich, Fundamentals of photonics (John Wiley & Sons, 2007).
C. Bohren and E. Clothiaux, Fundamentals of atmospheric radiation (Wiley-VCH, 2006).
H. C. van de Hulst, Light scattering by small particles (Dover Publications, 1981).

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.

Click here to see a list of articles that cite this paper