## Abstract

Particle-doped thin films that are translucent and diffusive have applications in cosmetics, coatings, and display technologies, but finding material combinations that produce these effects simultaneously is challenging: formulations tend to be either transparent or opaque. Using a combination of Mie scattering calculations and spectral transmission measurements on monodisperse colloidal suspensions, we demonstrate that the two characteristic optical properties of the films, total transmittance and haze, scale with the effective backscattering and forward scattering cross sections, both of which are properties of single particles. These scalings enable an efficient computational search for combinations of particle sizes, concentrations, and refractive indices that break the trade-off between translucency and diffusion. The optimum particle sizes and concentrations obey power-law dependences on the refractive index difference, a result of the interference condition for resonances in the scattering cross sections. The power laws serve as design equations for formulating particle-doped thin films.

© 2014 Optical Society of America

## 1. Introduction

Coatings and thin films of small particles suspended in a medium are important for many applications, including personal care products, paints, paper, display technologies, and polymer composites. Many of the most interesting applications require coatings that produce competing optical effects simultaneously. For instance, light diffuser plates in liquid crystal displays (LCD) [1] and ‘soft focus’ cosmetic products [2] require thin films that are both translucent and diffusive. That is, they must allow light to pass through, but that transmitted light should be scattered so as to appear diffuse. However, achieving translucency and diffusion simultaneously is challenging: as one increases the other tends to decrease. This trade-off, coupled with the wide variety of available ingredients, makes it difficult to formulate materials with the desired combination of translucency and diffusion by trial and error. Therefore, models and design rules able to predict bulk optical properties from single-particle attributes are needed.

A number of models describing the radiative transfer through particle-doped coatings have been developed over the past century [3]. However, all such models struggle to provide a connection between experimental parameters such as particle size, refractive index, and concentration, and their effects on the bulk optical properties of the composite film. On one hand, phenomenological models such as that of Kubelka and Munk [4, 5] can relate reflectance and transmittance of thin films to scattering and absorbing “constants,” but the connection between these constants and single-particle scattering parameters is unclear. Moreover, the constants are not directly related to physical properties of the material and depend on external factors such as whether the illumination geometry is diffuse or collimated [6]. On the other hand, Monte Carlo-based models [7] that use Mie scattering theory [8] to explicitly calculate the trajectories of incident photons can predict the full angular and spatial light intensity distributions from single-particle attributes, but they are computationally expensive and obscure simple physical insight into the most relevant variables.

Here we take an intermediate approach that combines spectral transmission measurements and Mie scattering calculations. With applications such as light diffuser plates and soft focus cosmetics in mind, we focus on coatings and thin films that maximize both total transmittance (translucency) and haze (diffusion). To systematically explore the relationship between transmittance and haze, we first measure the bulk optical properties of monodisperse colloidal suspensions, investigating the role of particle size, volume fraction, and refractive index. As we show, our measurements of transmittance and haze collapse when rescaled by two single-particle scattering quantities, the effective backscattering and forward scattering cross sections, showing that the film’s appearance is primarily due to single-scattering processes. We also show that particle-size-dependent resonances in these two cross sections can be exploited to break the trade-off between translucency and diffusion, yielding a set of design rules for formulating particle-doped thin films.

## 2. Methods

#### 2.1. Sample preparation

Monodisperse suspensions of polystyrene (PS) spheres with a variety of diameters are obtained from Life Technologies Corporation; poly(methyl methacrylate) (PMMA) particles are synthesized by surfactant-free emulsion polymerization (320-nanometer-diameter) or dispersion polymerization (2-micrometer-diameter)[9, 10]. All colloids are washed three times in Millipore water by centrifugation/resuspension and then resuspended in 10 mM NaCl/ Millipore water. Thin samples are then prepared by loading 20–30 microliters of suspension into a sample cell consisting of a precleaned glass slide and a No. 1 glass coverslip that are sealed on two sides with optical-grade, UV-curable epoxy (Norland Optical Adhesive 68, Norland Products, Inc.). The chamber thickness is controlled by Mylar® spacers (200 gauge, DuPont) epoxied between the slide and coverslip. Once fabricated, sample chambers are mounted onto the integrating sphere setup shown in Fig. 1.

#### 2.2. Spectral transmission measurements

The setup consists of a halogen lamp (DH-2000, Ocean Optics), an integrating sphere (RT-060-SF, Labsphere, Inc.), and a spectrometer (HR2000+, Ocean Optics). Unpolarized light from the halogen lamp is coupled to a multimode optical fiber (QP600-2-UV-Vis, Ocean Optics) and collimated by an achromatic fiber collimator (74-ACR, Ocean Optics) just upstream of the integrating sphere. After passing through the sample chamber, transmitted light is collected by the integrating sphere and sent to the spectrometer by an identical optical fiber.

Using the integrating sphere setup, we measure total transmittance *T* (the fraction of incident light that is transmitted through the sample) and diffuse transmittance *T*_{diffuse} (the fraction of incident light that is transmitted though the sample and scattered more than 4.75° from the direction of incident illumination) by changing the position of a light trap (Fig. 1). The recorded spectra are dark-count subtracted and normalized at each wavelength to a spectrum transmitted through a glass chamber filled with 10 mM NaCl/Millipore water. In an effort to minimize systematic errors due to variations in the average sphere reflectance, the sample, light trap, and reference chamber are all attached to the integrating sphere in each of three measurements taken per sample, only at different positions. For example, we place the light trap at the port hole opposite the sample chamber for measurements of diffuse transmittance, and at the port hole adjacent to the sample chamber for measurements of total transmittance. Removing the light trap from the integrating sphere entirely and replacing it with a port-hole plug for measurements of total transmittance would raise the average sphere reflectance, resulting in a systematic overestimation of *T* relative to *T*_{diffuse}. Using these measurements, we compute the haze, defined as the ratio of diffuse transmittance to total transmittance:

## 3. Results and discussion

#### 3.1. Relating bulk optical properties to single-particle scattering quantities

Our measurements of total transmittance and haze demonstrate a trade-off between translucency and diffusion. We measure *T* and *H* of aqueous suspensions of polystyrene or PMMA spheres of different diameters as a function of volume fraction *φ* (Figs. 2(a)–2(b)). For a given particle size, we find that both total transmittance and haze are monotonic functions of the particle concentration, with *T* decreasing and *H* increasing with increasing *φ*. This trend is consistent with our intuition: as the volume fraction (or number density) of particles increases, so does the likelihood of an incident photon being scattered by a particle.

Interestingly, however, we find that haze scales non-monotonically with particle size at fixed volume fraction: it first increases and then decreases with increasing particle diameter *D* (Fig. 2(b)). That is, there appear to be specific particle sizes for which haze is maximized (Fig. 2(b), inset). Total transmittance, in contrast, is only weakly dependent on particle size for the polystyrene samples studied here.

To quantify the relationship between particle size, volume fraction, and refractive index further, we construct a simple model that relates the observed bulk optical properties to single-particle scattering quantities. We first derive a expression for total transmittance, assuming that there are no structural correlations between particles and that effects due to multiple scattering are negligible. In the absence of absorption, attenuation of totally transmitted light is due to one of two processes: (*1*) scattering of light by more than 90° from the direction of incident illumination, or (*2*) total internal reflection of scattered light from the glass-air interface (Fig. 3).

That is, no light scattered at an angle larger than the critical angle for total internal
reflection *θ _{c}* =
sin

^{−1}(

*n*

_{air}/

*n*) is captured by the integrating sphere, where

_{m}*n*

_{air}is the refractive index of air and

*n*is the refractive index of the medium [15]. We assume that contributions due to specular reflection are negligible since we normalize the measured transmission to that of an identical glass sample cell. Therefore, light of intensity

_{m}*I*should be attenuated by an amount

*dx*, where

*π*/2 <

*θ*<

*π*) and total internal reflection (

*θ*<

_{c}*θ*<

*π*/2);

*ρ*is the particle number density. Like the Beer-Lambert law, Eq. (2) comes from considering the scattering from a thin slab of thickness Δ

*x*with area

*A*. In our case, we are concerned only with scattering at angles larger than the critical angle associated with total internal reflection. Therefore, the attenuation of intensity due to scattering is Rearranging and taking the limit as Δ

*x*goes to 0 yields

*dσ*

_{sca}/

*d*Ω is proportional to the amount of light scattered into a unit solid angle about a given direction and can be evaluated exactly for scattering of electromagnetic radiation by sphere, given by Mie's solution to Maxwell's equations [8]. Integrating Eq. (2) over the film thickness

*L*yields an expression for the total transmittance in which

*T*scales exponentially with a dimensionless length scale ${\sigma}_{b}^{\text{eff}}\rho L$. The expression that we find is analogous to the Beer-Lambert law [16], except that the total scattering cross section has been replaced by the effective backscattering cross section.

We derive a similar expression for haze, relying again on the same assumptions. Our original definition of haze (given by Eq. (1)) can be rewritten as *H*/100% = 1 − *T*_{direct}/*T*, where *T* = *T*_{diffuse} + *T*_{direct} and *T*_{direct} is the fraction of incident light transmitted at *θ* = 0. Determining *H*, then, reduces to calculating the fraction of transmitted light that is transmitted directly through the particle-doped film. In this case, in a slab of thickness *dx*, directly transmitted light of intensity *I* is attenuated by an amount

*σ*

_{sca}is the scattering cross section, which is proportional to the total amount of scattered light. Integrating Eq. (4), we find that the fraction of light transmitted directly through the film is

*T*

_{direct}/100% = exp(−

*σ*

_{sca}

*ρL*), recovering exactly the Beer-Lambert law and yielding the following expression for haze: Plugging in our solution for transmittance, we find that haze should scale exponentially with a different dimensionless length scale ${\sigma}_{f}^{\text{eff}}\rho L$ according to in which we define an effective forward scattering cross section

*σ*

_{sca}, the effective forward scattering cross section is roughly equivalent to the total scattering cross section, and the dimensionless length scale ${\sigma}_{f}^{\text{eff}}\rho L\approx {\sigma}_{\text{sca}}\rho L$ represents the average number of scattering events experienced by a photon as it traverses the film. The length scale 1/(

*σ*

_{sca}

*ρ*) is often called the scattering length or scattering mean free path and can again be calculated exactly for suspensions of spheres.

Remarkably, we find that all of our measurements of total transmittance and haze collapse onto two master curves when plotted against the two dimensionless length scales
${\sigma}_{b}^{\text{eff}}\rho L$ and
${\sigma}_{f}^{\text{eff}}\rho L$ (Figs. 2(c)–2(d)). The collapse demonstrates that *T* is a function of
${\sigma}_{b}^{\text{eff}}\rho L$ only, and *H* is a function of
${\sigma}_{f}^{\text{eff}}\rho L$ only. In other words, total transmittance is determined almost entirely by backscattering and total internal reflection, while haze is determined almost entirely by forward scattered light that is not totally internally reflected. Moreover, because our data are captured quantitatively by Eqs. (3) and (5) without any free parameters, we conclude that single scattering does indeed dominate to a first approximation, even for films with a white appearance. Only above
${\sigma}_{b}^{\text{eff}}\rho L\approx 1$ do our data of total transmittance begin to deviate from Eq. (3), likely because of multiple scattering. Indeed, calculations of the transport length, defined as (*σ*_{sca}*ρ*)^{−1}/(1 − *g*), where *g* is the asymmetry parameter [16], show that the transport length can be two to three times shorter than the sample thickness for the most concentrated suspensions. Models using the diffusion approximation have been developed to describe transmittance in the multiple scattering regime [17], but they are only applicable when the transport length is more than ten times shorter than the sample thickness, which is not the case in our experiments.

Ultimately, Eqs. (3) and (5) allow us to calculate the bulk optical properties of particle-doped films and evaluate their performance computationally. Because we hope to design films that maximize total transmittance and haze simultaneously, we seek suspensions that maximize *σ*_{sca}*ρL* while minimizing
${\sigma}_{b}^{\text{eff}}\rho L$. That is, we seek combinations of particles and media that scatter strongly (maximize *σ*_{sca}), but do so primarily in the forward direction (minimize
${\sigma}_{b}^{\text{eff}}$). Breaking the trade-off between total transmittance and haze, then, requires not only maximizing the amount of scattered light, but also engineering the angular distribution of that scattered light.

#### 3.2. Developing design rules

Following our experimental result that particle size is key to breaking the trade-off between translucency and diffusion, we use Eqs. (3) and (5), along with Mie scattering calculations to calculate the optical properties of thin films composed of different particle-medium combinations. Specifically, we seek to maximize the geometric mean of total transmittance and haze,
$\sqrt{H\times T}$, by varying the volume fraction and particle size for particle-medium combinations with different refractive index contrasts. We choose the geometric mean as the figure of merit because it is more sensitive to extreme differences between *T* and *H* than the arithmetic mean. For example, a thin film characterized by a total transmittance of 0% could have an arithmetic mean of *T* and *H* as high as 50% even though it would be practically useless for many applications. The geometric mean, however, would be 0%, which is much more characteristic of the appearance of the film.

Representative calculations, shown in Fig. 4, highlight the challenge of designing particle-doped films that are both translucent and diffusive:
$\sqrt{H\times T}$ is large only in a small region of parameter space. At fixed particle size,
$\sqrt{H\times T}$ first increases and then decreases with increasing volume fraction. As one might expect, the film is transparent if there are too few particles, and opaque if there are too many. This competition gives rise to an optimum volume fraction, which we call *φ*^{*}. At fixed volume fraction, the geometric mean
$\sqrt{H\times T}$ also displays a maximum with respect to particle size, reproducing the trend that we observe in experiment. This trend is less intuitive than the volume-fraction dependence, but still gives rise to an optimum particle diameter *D*^{*}. The (*D*^{*}, *φ*^{*}) coordinates are marked with an orange x in Fig. 4.

To explore the physical origin of the optimum particle size, we calculate
*D*^{*} for particle-medium combinations
spanning a range of refractive index differences Δ*n*
= *n _{p}* −

*n*from 0.03 to 0.3, where

_{m}*n*is the particle’s refractive index; examples for polystyrene in water (Δ

_{p}*n*= 0.26) and PMMA in water (Δ

*n*= 0.16) are shown in Fig. 4. Comparing Fig. 4(a) to Fig. 4(b), we see immediately that the region of parameter space corresponding to high $\sqrt{H\times T}$, and thus the value of the optimum diameter, depends on the specific particle-medium combination: the best particle size for one particle-medium combination might differ from another. However, we find that the computed optimum diameters follow a master curve, scaling as a power law with the refractive index difference Δ

*n*for all particle-medium combinations (Fig. 5(a)). The scaling law that we find is fitted by

Although the scaling demonstrates that the optimum particle size is dominated by the refractive index difference, it does not provide a physical explanation. To understand the basis of the inverse relation between *D*^{*} and Δ*n*, we return to our qualitative design criteria: to maximize total transmittance and haze simultaneously, we need suspensions that scatter strongly, particularly in the forward direction. The scattering efficiency *Q*_{sca}, defined as *σ*_{sca}/(*πa*^{2}), is enhanced for particles of particular sizes relative to the wavelength of light due to constructive interference between the incident and scattered light, where *a* is the particle radius. According to Mie scattering theory, these resonances occur when

*m*=

*n*is the refractive index contrast,

_{p}/n_{m}*p*is an integer, and

*x*= 2

*πan*is the size parameter, which represents the size of the particle relative to the wavelength of light

_{m}/λ*λ*in the medium [16]. Rearranging Eq. (7), we see immediately that, for a fixed wavelength, the particle sizes corresponding to the resonances scale as

*D*∝ Δ

*n*

^{−1}, confirming that the non-monotonic dependence of haze on particle size is due to resonances in the scattering cross section. The effective backscattering cross section does not contribute significantly to the scaling since it is 1–2 orders of magnitude smaller than

*σ*

_{sca}for particle sizes near

*D*

^{*}.

We also find that the optimum volume fraction *φ*^{*} follows a power law scaling with the refractive index difference and is fit well by

*σ*

_{sca}and ${\sigma}_{b}^{\text{eff}}$, depend only on Δ

*n*. The total scattering cross section scales as Δ

*n*

^{−2}, consistent with scattering from spheres with low refractive index contrast (

*m*< 2) [18]. We also find that the effective backscattering cross section at

*D*

^{*}is constant and independent of Δ

*n*. As a result, to maintain balance between the two dimensionless length scales, ${\sigma}_{b}^{\text{eff}}\rho L$ and

*σ*

_{sca}

*ρL*, we would expect

*φ*

^{*}to scale with an exponent somewhere between −1 and −3, which it does. The

*L*

^{−1}dependence of

*φ*

^{*}in Eq. (8) comes simply from the fact that the number of scatterers increases linearly with the sample thickness

*L*.

Taken together, Eqs. (6) and (8) serve as convenient design equations for formulating thin, particle-doped films with both translucency and diffusion. For a particle-medium combination of known Δ*n*, these equations can be used to select the optimum particle size and volume fraction needed for formulation. Instead, if it is desirable to minimize the particle loading, as is the case in many cosmetic applications, Eqs. (6) and (8) indicate that small particles with a higher refractive index difference are best.

## 4. Conclusions

We have shown that total transmittance and haze of colloidal suspensions depend only on two microscopic length scales under conditions where single scattering dominates. These length scales can be computed accurately using Mie scattering theory and controlled by varying the particle size, volume fraction, and refractive index, as well as the refractive index of the medium. Such quantitative relationships between the attributes of individual scatterers and the resulting bulk optical properties enable the rational design of particle-doped materials. For particular applications that simultaneously require high transmittance and high haze, seemingly contradictory requirements, we have shown that the optimum particle size and concentration follow power-law dependences on the refractive index difference between the particle and medium. These power-law scalings serve as simple design equations.

Although our analysis has considered scattering from spheres only, these concepts could also be extended to scattering from non-spherical particles [19, 20, 21]. Moreover, while we focused on designing suspensions that are both translucent and diffusive, a similar approach to deriving design criteria from Mie theory could also be used to choose particle specifications that maximize total transmittance while minimizing diffusion [22]. This might be desirable in composite materials where particles are added for their mechanical benefit, and light scattering is an unwanted side effect. More generally, having models that relate bulk optical properties to single-particle attributes should enable the rational design of other coatings whose appearance is dominated by light scattering.

## Acknowledgments

We thank Christian Mahé for helpful discussions, and Solène Fortun, Alice Lefebvre, and Rodrigo Guerra for providing materials. This work was supported by Chanel USA Inc. and by the Harvard MRSEC funded by the National Science Foundation through grant number DMR-0820484.

## References and links

**1. **P. Yeh and C. Gu, *Optics of Liquid Crystal Displays*, vol. 67 (John Wiley & Sons, 2010).

**2. **R. Emmert, “Quantification of the soft-focus effect,” Cosmetics and Toiletries **111**, 57–61 (1996).

**3. **S. Chandrasekhar, *Radiative Transfer* (Courier Dover Publications, 1960).

**4. **P. Kubelka and F. Munk, “Ein Beitrag zur Optik der Ein Beitrag zur Optik der Farbanstriche,” Z. Tech. Phys. (Leipzig) **12**, 593–601 (1931).

**5. **P. Kubelka, “New contributions to the optics of intensely light-scattering materials. Part I,” J. Opt. Soc. Am. **38**, 448 (1948). [CrossRef] [PubMed]

**6. **L. Yang and B. Kruse, “Revised Kubelka–Munk theory. I. Theory and application,” J. Opt. Soc. Am. A **21**, 1933–1941 (2004). [CrossRef]

**7. **M. Yamada, M. D. Butts, and K. K. Kalla, “Spatial and angular distribution of light incident on coatings using Mie-scattering Monte Carlo simulations,” J. Cosmet. Sci. **56**, 193–204 (2005). [PubMed]

**8. **G. Mie, “Beiträge zur optik trüber medien, speziell kolloidaler metallösungen,” Annalen der Physik (1908). [CrossRef]

**9. **M. Egen and R. Zentel, “Surfactant-free emulsion polymerization of various methacrylates: Towards monodisperse colloids for polymer opals,” Macromolecular Chemistry and Physics **205**, 1479–1488 (2004). [CrossRef]

**10. **S. Jiang, E. D. Sudol, V. L. Dimonie, and M. S. El-Aasser, “Seeding as a means of controlling particle size in dispersion polymerization,” Journal of Applied Polymer Science **108**, 4096–4107 (2008). [CrossRef]

**11. **S. N. Kasarova, N. G. Sultanova, C. D. Ivanov, and I. D. Nikolov, “Analysis of the dispersion of optical plastic materials,” Optical Materials **29**, 1481–1490 (2007). [CrossRef]

**12. **F. C. Cheong, K. Xiao, D. J. Pine, and D. G. Grier, “Holographic characterization of individual colloidal spheres’ porosities,” Soft Matter **7**, 6816–6819 (2011). [CrossRef]

**13. **J. Fung, K. E. Martin, R. W. Perry, D. M. Kaz, R. McGorty, and V. N. Manoharan, “Measuring translational, rotational, and vibrational dynamics in colloids with digital holographic microscopy,” Optics Express **19**, 8051–8065 (2011). [CrossRef] [PubMed]

**14. **R. W. Perry, G. Meng, T. G. Dimiduk, J. Fung, and V. N. Manoharan, “Real-space studies of the structure and dynamics of self-assembled colloidal clusters,” Faraday Discussions **159**, 211–234 (2012). [CrossRef]

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

**16. **C. F. Bohren and D. R. Huffman, *Absorption and Scattering of Light by Small Particles* (Wiley-VCH Verlag GmbH, 2004).

**17. **P. Kaplan, A. Dinsmore, A. Yodh, and D. Pine, “Diffuse-transmission spectroscopy: A structural probe of opaque colloidal mixtures,” Physical Review E **50**, 4827 (1994). [CrossRef]

**18. **H. C. Hulst and H. Van De Hulst, *Light Scattering by Small Particles* (Courier Dover Publications, 1957).

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

**20. **M. I. Mishchenko, J. W. Hovenier, and L. D. Travis, *Light Scattering by Nonspherical Particles: Theory, Measurements, and Applications* (Academic Press, 1999).

**21. **Z. Gimbutas and L. Greengard, “Fast multi-particle scattering: A hybrid solver for the maxwell equations in microstructured materials,” Journal of Computational Physics **232**, 22–32 (2013). [CrossRef]

**22. **A. Small, S. Hong, and D. Pine, “Scattering properties of core–shell particles in plastic matrices,” Journal of Polymer Science Part B: Polymer Physics **43**, 3534–3548 (2005). [CrossRef]