## Backscattering of light from disk-like particles with aperiodic angular fine structure

Optics Express, Vol. 15, Issue 25, pp. 16424-16430 (2007)

http://dx.doi.org/10.1364/OE.15.016424

Acrobat PDF (401 KB)

### Abstract

Recent computations of the backscattering cross section (σ_{
b
}) of randomly-oriented disk-like particles (refractive index, 1.20) with small-scale periodic angular internal structure, have been repeated for similarly sized particles, but with the periodic structure replaced by an *aperiodic* structure. The latter is formed by randomly perturbing a periodic structure. Although σ_{
b
} for individual realizations of an aperiodic disk can differ significantly from that of its periodic counterpart, averaging over several realizations brings the two into confluence, unless the aperiodicity is too large. These computations suggest that using disks with perfectly periodic (as opposed to quasi-periodic) fine structure for modeling the backscattering of detached coccoliths from *E. huxleyi* is justified.

© 2007 Optical Society of America

## 1. Introduction

1. H.R. Gordon and A.Y. Morel, *Remote Assessment of Ocean Color for Interpretation of Satellite Visible Imagery: A Review* (Springer-Verlag, 1983). [CrossRef]

2. D. Stramski, E. Boss, D. Bogucki, and K.J. Voss, “The role of seawater constituents in light backscattering in the ocean,” Prog. Oceanogr. **61**, 27–56 (2004). [CrossRef]

4. H.R. Gordon and Tao Du, “Light scattering by nonspherical particles: application to coccoliths detached from *Emiliania huxleyi*,” Limnol. Oceanogr. **46**, 1438–1454 (2001). [CrossRef]

*E. huxleyi*, which has a well-defined shape (resembling a disk or two roughly parallel disks) and a known composition (Calcite, refractive index relative to water ~1.20). (See Ref. 5 for scanning electron micrographs of

*E. huxleyi*coccoliths.) However,

*E. huxleyi*in fact has a rather complex fine structure that might influence backscattering and should be addressed. I tried to examine this in an earlier paper [5

5. H.R. Gordon, “Backscattering of light from disk-like particles: is fine-scale structure or gross morphology more important?,” Appl. Opt. **45**, 7166–7173
(2006). [CrossRef] [PubMed]

6. B.T. Draine, “The discrete-dipole approximation and its application to interstellar graphite grains,” Astrophys. J. **333:**848–872 (1988). [CrossRef]

7. B.T. Draine and P. Flatau, “Discrete-dipole approximation for scattering calculations,” J. Opt. Soc. Am. A **ll**, 1491–1499 (1994). [CrossRef]

*α*=2

*π*/2

^{n}, and removing the dipoles from alternate sectors. In that study the diameter of the disks ranged from 1.50 to 2.75 µm and the thickness from 0.05 to 0.15 µm. The values used for

*n*were 4, 5, 6, and 7, providing pinwheel-looking objects (Fig. 1) with 8, 16, 32, and 64 vanes, respectively. The principal result of the study was that when the scale of the periodicity, s (defined to be the length of an open or closed sector measured along the circumference of the disk), was <

*λ*/4, where

*λ*is the wavelength of the light

*in*the medium (water) the backscattering was found to be nearly identical to that of a homogeneous disk possessing a reduced refractive index. In contrast, significant increases in backscattering were observed when the scale of the periodicity was greater than

*λ*/4, reaching a maximum when the scale becomes ~

*λ*/2 [8

8. H.R. Gordon, “Rayleigh-Gans scattering approximation: surprisingly useful for understanding backscattering from disk-like particles,” Opt. Express **15**, 5572–5588 (2007). [CrossRef] [PubMed]

5. H.R. Gordon, “Backscattering of light from disk-like particles: is fine-scale structure or gross morphology more important?,” Appl. Opt. **45**, 7166–7173
(2006). [CrossRef] [PubMed]

*α*of individual sectors are introduced.

## 2. Model of an aperiodic pinwheel

*angles*

^{n}*αP*. The individual boundary angles are then perturbed to

*α*according to

_{I}*ε*≤1 is a constant and -1/2≤

*ρ*≤1/2 is a random number with a uniform probability density. Then, the material of the disk is removed from ever other sector, yielding a pinwheel with a quasi-periodic structure. Four realizations (each based on a difference sequence of pseudorandom numbers) of such pinwheels for

*n*=5 are provided in Fig. 1 for

*ε*=0.5 and 1.0. Defining ∑

_{I}to be the standard deviation in the angle

*α*, we find ∑

_{I}*=*

_{I}^{12}-½εΔ

*α*≈0.3

*ε*Δ

*α*, where Δ

*α*=2

*π*/2

*. Likewise defining ∑*

^{n}*to be the standard deviation of the removed (or occupied) sector angles,*

_{Δα}*∑*=6-½

_{Δα}*ε*Δα≈0.4

*ε*Δα. Thus, for

*ε*=0.5 and 1.0, the relative standard deviation in angle of the removed (or occupied) sectors is 20 and 40%, respectively. I examined two aperiodic pinwheels. The first has a diameter (

*D*) of 1.50 µm, a thickness (

*t*) of 0.15 µm, and

*n*=5 (Fig. 1). The second has

*D*=2.75 µm,

*t*=0.05 µm, and

*n*=6. The larger disk is similar in size to the distal shield of individual

*E. huxleyi*coccoliths [5

5. H.R. Gordon, “Backscattering of light from disk-like particles: is fine-scale structure or gross morphology more important?,” Appl. Opt. **45**, 7166–7173
(2006). [CrossRef] [PubMed]

*ρ*begins. Thus, for realization 0000 the sampling begins with the first number, for realization 1000 it begins with the 1001

^{th}number, etc. The variation of the volume for a given

*ε*can be as much as 25% for the smaller (

*D*=1.5 µm,

*n*=5) and 7% for the larger (

*D*=2.75 µm,

*n*=6) pinwheels. The reduction in dispersion from the smaller to the larger is due to the increase in

*n*, which doubles the number of sectors, increasing the probability that the individual realizations have a volume closer to the mean.

## 3. Operation of the discrete-dipole scattering code

*d*, the more dipoles are required to fill the volume of the particle. A convenient measure of the spacing in regard to the wavelength is |

*m*|

*kd*, where

*m*is the refractive index and

*k*=2

*π/λ*. The orientation of a disk-like object is specified by three angles:

*θ*the angle the axis of the disk makes with the incident beam,

*ϕ*, the azimuth of the axis relative to a laboratory-fixed plane containing the incident beam, and

*β*, the angle of rotation around the axis required to place the disk in a specified orientation given

*θ*and

*ϕ*. In general, for an object with no rotational symmetry, e.g., an aperiodic pinwheel, 0°≤

*θ*≤180°, 0°≤

*ϕ*≤360° and 0°≤

*β*≤360° ; however, the high symmetry of a uniform disk reduces these to 0°≤

*θ*≤90°, 0°

*≤*ϕ≤180°, and requires only one value of

*β*, e.g.,

*β*=0. For periodic pinwheels (

*ε*=0) the angle

*β*is required, however, its range need only be enough to completely cover one open sector and one adjacent occupied sector. The DDA code performs orientational averaging by computing the scattering at discrete angles equally spaced in

*ϕ*and

*β*, while the angle

*θ*is divided in uniform increments of cos

*θ*. An important consideration in the averaging is that the computation time is roughly proportional to the number of angles (

*N*×

_{θ}*N*×

_{ϕ}*N*) used in the averaging.

_{β}8. H.R. Gordon, “Rayleigh-Gans scattering approximation: surprisingly useful for understanding backscattering from disk-like particles,” Opt. Express **15**, 5572–5588 (2007). [CrossRef] [PubMed]

*D*=2.7 µm using ~5000 orientations (

*N*=51,

_{θ}*N*=99,

_{ϕ}*N*=1) for the orientational average, the error in the backscattering cross section (

_{β}*σb*) was of the order of 5% for |

*m*|

*kd*<0.5, and decreased rapidly for smaller values of |

*m*|

*kd*. In the present work, I have always used enough dipoles to keep |

*m*|

*kd*<0.5 (and often <0.4) and a number of orientations that would provide approximately the same averaging accuracy as the uniform disk with

*D*=2.7 µm. For the periodic pinwheels I used

*N*=51,

_{θ}*N*=99 and

_{ϕ}*N*=4. The four

_{β}*β*increments were equally spaced over one open sector and an adjacent occupied sector. Because of the symmetry, this would be equivalent to choosing

*N*=128 and 0°≤

_{β}*β*≤360°. For the aperiodic pinwheels (with

*ε*=1), I used

*N*=101,

_{θ}*N*=99 and

_{ϕ}*N*=9. These latter orientations were chosen by extensive testing in which I examined the differential scattering cross section for scattering angles Θ>90°. The region 120°≤Θ≤180°, the phase functions of periodic and aperiodic pinwheels (

_{β}*D*=2.75 µm) are highly oscillatory in Θ, and for the above choice of the orientational averaging, the observed oscillations in the periodic and aperiodic cases were similar, facilitating a fair comparison between the two. Actually,

*σ*appears to become stable with smaller values of

_{b}*N*and

_{θ}, N_{φ}*N*, where the orientationally-averaged phase function becomes more strongly dependent on the number of orientations, providing confidence that a sufficient number of orientations have been used. In the final analysis, the choice of the number of orientations must be balanced against computational time required.

_{β}## 4. Results of the computations

*D*and

*t*). Note that the homogeneous disk has twice the volume (mass) of the periodic pinwheel (

*ε*=0) and approximately twice the volume of the aperiodic pinwheels (

*ε*>0).

*s*=

*λ*/4 occurs when

*t/λ*=0.255 for the smaller and 0.093 for the larger pinwheel. For

*t/λ*larger than these values, the backscattering increases rapidly with decreasing

*λ*and then undergoes a series of maxima and minima with progressively increasing backscattering at each maximum. The backscattering at the maxima for the smaller pinwheel is, in magnitude, approximately that at the maxima for the homogeneous disk (twice the volume or mass of the pinwheel). For the larger pinwheel, it is approximately 75% of the maxima for a uniform disk. These pinwheel maxima are the result of interference of the fields scattered by the individual vanes of the pinwheel as they occur in the Rayleigh-Gans approximation as well (although only the first maximum occurs at the same position [8

8. H.R. Gordon, “Rayleigh-Gans scattering approximation: surprisingly useful for understanding backscattering from disk-like particles,” Opt. Express **15**, 5572–5588 (2007). [CrossRef] [PubMed]

*ε*=0.5 its backscattering closely follows the periodic pinwheel, with the dispersion of backscattering reaching 20% at the smallest wavelength, while when

*ε*=1.0 there is more significant deviation from the periodic case, and the dispersion is somewhat larger. Recalling that

*∑Δα*≈0.4

*ε*Δα=0. 4×2

*πε*/2

*, we note that the smaller pinwheel with*

^{n}*ε*=0.5 (

*n*=5) has the same value of

*∑*as the larger pinwheel with

_{Δα}*ε*=1.0 (

*n*=6). Interestingly, Figs. 2 (Left Panel) and 3 show that the behavior of

*σ*with decreasing

_{b}*λ*up to its first minimum are similar in the two cases (

*ε*=0.5,

*n*=5 and

*ε*=1.0,

*n*=6): there are only minor deviations from the periodic pinwheels; and there is small dispersion between the various realizations of the aperiodic pinwheels. In contrast, when

*ε*is increased from 0.5 to 1.0 for the smaller pinwheel, the dispersion increases, and

*σ*near its first maximum (

_{b}*t/λ*≈0.4) shows a significant decrease from the

*ε*=0 and 0.5 cases. This behavior would be expected under the hypothesis that the maxima in the periodic case results from constructive interference of light interacting with the individual vanes of the pinwheel — when the spacing and angular size of the vanes becomes random the constructive interference is reduced.

*σ*for the several realizations is ordered with increasing volume, i.e., the realization with the smallest

_{b}*σ*has the smallest volume, etc. However, with increasing

_{b}*λ*(decreasing

*t/λ*) the order can reverse, particularly near the first maximum near

*t/λ*≅0.4. Thus, there is no way to try to reduce the dispersion throughout the whole wavelength range by normalizing for volume differences. In the Rayleigh-Gans domain,

*t/λ*<0.2 [8

**15**, 5572–5588 (2007). [CrossRef] [PubMed]

*σ*is proportional to the square of the volume. This would explain the envelope of the variation near

_{b}*t/λ*=0.2 in Fig. 2.

*E huxleyi*coccoliths) would display similar variations in volume. In fact, if pinwheels were to represent real biological particles, samples would be expected to consist of a number of realizations of their aperiodicity. In this regard, the average

*σ*(denoted by 〈

_{b}*σ*〉) is more important than that for any given realization. Figure 4 compares the 〈

_{b}*σ*〉 for the four realizations of the aperiodicity examined here with the associated periodic pinwheel. It clearly shows that the main difference between 〈

_{b}*σ*〉 for the small periodic and aperiodic pinwheels (left panel) occurs near the maxima in the backscattering, and that near the first maximum (but not the second) the difference increases as deviation from periodicity increases (i.e., as

_{b}*ε*increases). Considering the large departures from periodicity for the

*ε*=1 realizations in Fig. 1, it is remarkable that, when averaged over realizations, their 〈

*σ*〉 is so close to that of periodic pinwheels (

_{b}*ε*=0). As suggested in Fig. 3, Fig. 4 (right panel) shows that the 〈

*σ*〉 for the large aperiodic pinwheel is very close to its periodic counterpart.

_{b}## 5. Discussion

*σ*changes much less than compared to the transition from uniform disk to periodic pinwheel. Second, the dispersion in

_{b}*σ*among realizations of the aperiodic pinwheels is associated with the dispersion in

_{b}*Δα*(or

*s*), which increases with increasing

*ε*. Third, the aperiodic 〈

*σ*〉 will usually be somewhat smaller than the periodic

_{b}*σ*, at least near the position of the first (long-wave) maximum, and this decrease increases with increasing aperiodicity (

_{b}*ε*).

**15**, 5572–5588 (2007). [CrossRef] [PubMed]

*∑Δα/Δα*~0.27, and there were 40 open angular sectors. This coccolith shield is similar in size and shape to the larger (2.75 µm) pinwheel (

*n*=6, 32 open sectors) examined here. The computations for the larger pinwheel show that, for the purpose of computing backscattering, the periodic pinwheel is a good approximation to aperiodic pinwheel as long as

*∑*≤0.4 (

_{Δα}/Δα*ε*≤1). This suggests that replacing the aperiodic fine structure of the distal shield of

*E. huxleyi*coccoliths with a strictly periodic fine structure will not degrade the modeling of their backscattering, especially for natural samples containing large numbers of coccoliths.

## Acknowledgments

## References

1. | H.R. Gordon and A.Y. Morel, |

2. | D. Stramski, E. Boss, D. Bogucki, and K.J. Voss, “The role of seawater constituents in light backscattering in the ocean,” Prog. Oceanogr. |

3. | M.I. Mishchenko, L.D. Travis, and A.A. Lacis, |

4. | H.R. Gordon and Tao Du, “Light scattering by nonspherical particles: application to coccoliths detached from |

5. | H.R. Gordon, “Backscattering of light from disk-like particles: is fine-scale structure or gross morphology more important?,” Appl. Opt. |

6. | B.T. Draine, “The discrete-dipole approximation and its application to interstellar graphite grains,” Astrophys. J. |

7. | B.T. Draine and P. Flatau, “Discrete-dipole approximation for scattering calculations,” J. Opt. Soc. Am. A |

8. | H.R. Gordon, “Rayleigh-Gans scattering approximation: surprisingly useful for understanding backscattering from disk-like particles,” Opt. Express |

**OCIS Codes**

(010.4450) Atmospheric and oceanic optics : Oceanic optics

(290.1350) Scattering : Backscattering

**ToC Category:**

Scattering

**History**

Original Manuscript: October 31, 2007

Revised Manuscript: November 20, 2007

Manuscript Accepted: November 20, 2007

Published: November 27, 2007

**Citation**

Howard R. Gordon, "Backscattering of light from disk-like particles with aperiodic angular fine structure," Opt. Express **15**, 16424-16430 (2007)

http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-25-16424

Sort: Year | Journal | Reset

### References

- H. R. Gordon and A. Y. Morel, Remote Assessment of Ocean Color for Interpretation of Satellite Visible Imagery: A Review (Springer-Verlag, 1983). [CrossRef]
- D. Stramski, E. Boss, D. Bogucki, and K. J. Voss, "The role of seawater constituents in light backscattering in the ocean," Prog. Oceanogr. 61, 27-56 (2004). [CrossRef]
- M. I. Mishchenko, L. D. Travis, and A. A. Lacis, Scattering, Absorption, and Emission of Light by Small Particles (Cambridge, 2002).
- H. R. Gordon and T. Du, "Light scattering by nonspherical particles: application to coccoliths detached from Emiliania huxleyi," Limnol. Oceanogr. 46, 1438-1454 (2001). [CrossRef]
- H. R. Gordon, "Backscattering of light from disk-like particles: is fine-scale structure or gross morphology more important?," Appl. Opt. 45, 7166-7173 (2006). [CrossRef] [PubMed]
- B. T. Draine, "The discrete-dipole approximation and its application to interstellar graphite grains," Astrophys. J. 333, 848-872 (1988). [CrossRef]
- B. T. Draine and P. Flatau, "Discrete-dipole approximation for scattering calculations," J. Opt. Soc. Am. A ll, 1491-1499 (1994). [CrossRef]
- H. R. Gordon, "Rayleigh-Gans scattering approximation: surprisingly useful for understanding backscattering from disk-like particles," Opt. Express 15, 5572-5588 (2007). [CrossRef] [PubMed]

## 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.