Technical note: Characterizing individual milk fat globules with holographic video microscopy
Abstract.
We use inline holography to create images of individual milk fat globules in diluted samples of milk. Analyzing these images with the exact LorenzMie light scattering theory then yields the droplets' radii with nanometer resolution and their refractive indexes to within one part in a thousand. This procedure rapidly and directly characterizes both the quantity and quality of fat in milk.
Optical characterization methods are widely used to monitor and optimize the properties of milk during processing and packaging ((5)). Although rapid and noninvasive, such methods generally provide only an indirect view of such properties as the concentration, size distribution and composition of the milk fat globules in a sample. Holographic video microscopy ((14); (7); (8)) offers a complementary and highly quantitative view of individual globules' optical properties, yet is rapid enough to provide a statistical overview of a bulk sample's properties. In particular, this method yields the radius of individual globules with nanometer resolution and the refractive index to within one part in a thousand. A statistical sampling of the globule radii provides a direct measurement of a sample's fat concentration and size distribution. The distribution of refractive indexes is useful for characterizing the globules' composition, and thus the nature and quality of the milk.
Our holographic microscope, depicted schematically in Fig. 1, is based on a commercial inverted light microscope (Zeiss Axiovert TV 100 S). The conventional incandescent illuminator is replaced with the 3 mm diameter collimated beam from a 10 mW HeNe laser (Uniphase) operating at a vacuum wavelength of . Such weak laser irradiation will not appreciably raise the temperature of an aqueous sample.
An individual milk fat globule of radius at position scatters a small portion of the planewave illumination. The scattered light then propagates to the focal plane of the microscope, where it interferes with the unscattered portion of the laser beam. The resulting interference pattern is magnified by the microscope's objective lens (S Plan Apo, , NA 1.4, oil immersion) and projected by a video eyepiece () onto a CCD camera (NEC TI324AII) before being recorded as uncompressed digital video with a digital video recorder (Panasonic DVRH110).
Each holographic image in the video stream is a timeresolved snapshot of the threedimensional distribution of scatterers in the microscope's field of view. We then use results of the LorenzMie theory of light scattering by small particles ((2); (1); (11)) to measure each particle's radius and index of refraction ((8)).
We used this method to analyze milk fat droplets from a range of commercial milk products including several grades of homogenized pasteurized cows' milk and goats' milk. In each case, the sample was diluted by 1000:1 with deionized water before being sealed in the 50 gap between a microscope slide and a glass cover slip and mounted on the stage of the microscope. Given the imaging system's calibrated magnification of 101 nm/pixel, a typical image, , captures roughly 10 resolvable globules.
Unprocessed holograms suffer from large intensity variations due to speckle, interference effects in the microscope's optics and scattering by dust and other imperfections. We correct for these by normalizing with a background hologram obtained with no sample in the field of view. The normalized hologram then can be fit to the prediction of LorenzMie theory ((2); (1); (11)),
(1) 
where is the wavenumber of light in a medium of refractive index , and where is the LorenzMie scattering function describing how light of polarization is scattered by a sphere located at . In practice, the illuminating beam is not perfectly uniform, and the factor can be used to account for positiondependent variations. Assuming the scatterer to be a uniform and homogeneous dielectric sphere illuminated by light that is linearly polarized in the direction,
(2) 
where , and where and are vector spherical harmonics ((2); (1)),
(3) 
and
(4) 
Here, is the associated Legendre polynomial of the first kind, and is the spherical Bessel function of the first kind of order . The expansion coefficients are given by by ((2))
(5) 
and
(6) 
where is the particle's refractive index relative to the medium, is the spherical Bessel function of the first type of order , is the spherical Hankel function of the first type of order , and where primes denote derivatives with respect to the argument. The sum in Eq. (2) converges after a number of terms, , which depends on the particle's size ((1); (11)). The only challenge in this computation is to calculate the Bessel functions and their ratios both accurately and efficiently. We use the accurate, but computationally intensive continued fraction algorithm due to Lentz to compute the and coefficients ((9)), and the more efficient recursive algorithm due to Wiscombe for the spherical Bessel functions ((15)). This tradeoff ensures that we can obtain accurate results using our apparatus for spheres ranging in diameter from 10 to more than 10 and refractive indexes exceeding .
To characterize a sphere, we fit its normalized hologram to Eqs. (1) through (6) for the sphere's radius , its refractive index , its threedimensional position , and also , using a standard LevenbergMarquardt leastsquares algorithm ((13); (12); (4)). Despite the fairly large number of free parameters, these fits converge rapidly and robustly, and typically yield the particle's size and position with nanometerscale resolution, and its refractive index to within one part in a thousand ((8)). Such extreme subpixel spatial resolution is possible because information from many pixels contributes to each fit.
We applied this technique to 5 samples of commercially processed milk obtained from a local supermarket. These include pasteurized homogenized cows' milk from Elmherst Dairy (Jamaica, NY) with designated fat contents ranging from fatfree to whole milk, as well as pasteurized homogenized goats' milk from Meyenberg Goat Milk (Turlock, CA). Samples were allowed to equilibrate on the microscope stage to ambient temperature, . Up to 100 randomly selected fat droplets were analyzed for each sample to obtain estimates for the size and refractive index distributions for the fat droplets in each sample. The results are summarized in Fig. 2, Fig. 3 and Table 1. Figure 2 presents histograms of the milk fat globule radii and refractive indexes for each of the samples. The mean values presented in Table 1 are computed as averages of the singledroplet values. Reported errors are the quadrature sum of the standard deviation of the values and the measurement error estimated from the fits.
Sample  Radius []  Refractive Index 
Cow Whole  
Cow 2%  
Cow 1%  
Cow Fat Free  
Goat Whole  
Goat 1 month  
Type B Immersion Oil 
As might be expected, the average radius of the fat globules in cows' milk is smaller in samples with lower fat content. This is consistent with the creaming process by which fat is removed from whole milk. The mean refractive index of the cows' milk fat globules varies far less from sample to sample. The mean refractive index of , averaged over all of the droplets from all of the cows' milk samples, agrees with bulk values obtained by Abbé refractometry on melted fat samples ((6)) and with sampleaveraged singledroplet refractive indexes obtained by light scattering ((10)). This value substantially exceeds the range of 1.3444 to 1.3525 obtained for the overall refractive index of bulk samples ((5)), which is dominated by the optical properties of water, whose refractive index at room temperature is .
Estimating the milk fat globules' refractive index from its contribution to light scattering or conventional refractometry measurements requires accurate information regarding the droplets' concentration and size distribution. Holographic characterization, by contrast, unambiguously quantifies individual droplets' optical properties without recourse to models or ancillary measurements. Because a droplet's refractive index depends on its composition, moreover, these measurements also can be used to characterize the composition and quality of milk samples. Indeed, the data in Table 1 and Fig. 3 demonstrate that it is possible to distinguish goats' milk from cows' milk on this basis, the goats' milk having a resolvably smaller mean refractive index, the difference being confirmed at the 95 percent confidence level by a Welch's ttest.
The high resolution afforded by holographic characterization allows us to track subtle changes in a sample as the milk ages. Table 1 shows that the mean refractive index of goat milk droplets increases over the course of a month, rising from 1.43 to 1.45. This coincides with a decrease in the mean droplet radius from 0.58 to 0.44 . Both changes are resolved at the 95 confidence level according to Welch's ttest. This result demonstrates that holographic characterization can be useful not only for distinguishing types of milk, but also may be used to assess the age of a sample.
Having access to particleresolved data also reveals an interesting correlation between the size of the fat globules and the range of their estimated refractive indexes. Larger particles' apparent refractive indexes are consistent with each other to within the resolution of the measurement technique ((8)). Smaller droplets, by contrast, display substantially larger variations in refractive index from droplet to droplet.
This is not an inherent limitation of the measurement technique, as the data in Fig. 4 demonstrate. The circular data points in Fig. 4 were obtained for spherical droplets of Cargille Type B microscope immersion oil that were transiently dispersed in deionized water by vigorous shearing and then diluted with deionized water to obtain droplet concentrations and size distributions qualitatively similar to those of the milk samples. Although these surfactantfree dispersions are thermodynamically unstable, their recombination kinetics are sufficiently slow that we are able to obtain clear holograms of the individual droplets before they coalesce. The immersion oil has a nominal bulk refractive index of 1.515 for red light at a temperature of 25. The measured values differ from this by no more than 2 percent over the entire range of droplet radii that were sampled.
Results more reminiscent of those for milk fat droplets are obtained when the oil droplets are stabilized with surfactant. The square points in Fig. 4 were obtained for Type B oil with the addition of 0.1% (v/v) Tergitol NP9, a nonionic surfactant whose bulk refractive index is 1.491. The addition of this lowerindex surfactant systematically reduces the singledroplet refractive index for larger particles. It also increases the variability of measured refractive indexes to more than 10 percent for particles with radii smaller than 0.6 .
Because the excess of surfactant guarantees complete coverage of the oil droplets' surfaces, these observations demonstrate that holographic particle characterization is sensitive to surface coverage by surfactants. In the case of milk droplets, this suggests that holographic microscopy is influenced by the optical properties of the milk fat globule membrane (MFGM) ((3)). This sensitivity is noteworthy because, at just 10 to 20 thickness, the MFGM is much thinner than the wavelength of light and constitutes a very small proportion of the droplets' volume.
The LorenzMie light scattering formulas in Eq. (2) through (6) are appropriate for a homogeneous isotropic sphere with an abrupt interface. Using this result to interpret holograms of coated spheres consequently can lead to inconsistencies in the extracted parameters. This effect should be more pronounced for smaller spheres whose surfaceareatovolume ratio is higher. Applying a more sophisticated form of the scattering function that accounts for coreshell structure ((2); (1)) should reduce this variability at the expense of additional computational complexity. The benefit would be the ability to characterize the MFGM.
The basic implementation of holographic droplet characterization presented here provides a dropletbydroplet analysis of milk fat droplets' sizes and also their refractive indexes at a particular wavelength of light. This information, in turn, offers insights into the droplets' composition, for example enabling us to distinguish goats' milk from cows' milk. Holographic characterization can be extended to multiple wavelengths, which would provide true spectroscopic information on milk fat globules' composition. This technique requires little specialized equipment, and so can be easily adapted for process control and quality assurance applications.
We are grateful to Osman Akcakir for bringing (9) to our attention. This work was supported by the National Science Foundation under Grant Number DMR0606415.
References

(1)
Barber, P. W. and S. C. Hill, 1990. Light Scattering by Particles: Computational Methods, volume 2 of Advanced Series in Applied Physics. World Scientific, New Jersey.

(2)
Bohren, C. F. and D. R. Huffman, 1983. Absorption and Scattering of Light by Small Particles. Wiley Interscience, New York.

(3)
Dewettinck, K., R. Rombaut, N. Thienpont, T. T. Le, K. Messens, and J. Van Camp, 2007. Nutritional and technological aspects of milk fat globule membrane. Int. Dairy J. 18:436–457.

(4)
Gill, P. E. and W. Muray, 1978. Algorithms for the solution of the nonlinear leastsquares problem. SIAM J. Numer. Anal. 15:977–992.

(5)
Jääskeläinen, A. J., K.E. Peiponen, and J. A. Räty, 2001. On reflectometric measurement of a refractive index of milk. J. Dairy Sci. 84:38–43.

(6)
Krukovsky, V. N., 1961. Vitamin A, carotenoid, iodine, and thiocyanogen values, and the refractive index of milk fat as influenced by feed, and by individual and breed differences. J. Ag. Food Chem. 9:326–330.

(7)
Lee, S.H. and D. G. Grier, 2007. Holographic microscopy of holographically trapped threedimensional structures. Opt. Express 15:1505–1512.

(8)
Lee, S.H., Y. Roichman, G.R. Yi, S.H. Kim, S.M. Yang, A. van Blaaderen, P. van Oostrum, and D. G. Grier, 2007. Characterizing and tracking single colloidal particles with video holographic microscopy. Opt. Express 15:18275–18282.

(9)
Lentz, W. J., 1976. Generating Bessel functions in Mie scattering calculations using continued fractions. Appl. Opt. 15:668–671.

(10)
Michalski, M.C., V. Briard, and F. Michel, 2001. Optical parameters of milk fat globules for laser light scattering measurements. Lait 81:787–796.

(11)
Mishchenko, M. I., L. D. Travis, and A. A. Lacis, 2001. Scattering, Absorption and Emission of Light by Small Particles. Cambridge University Press, Cambridge.

(12)
Moré, J. J., 1977. The LevenbergMarquardt algorithm: Implementation and theory. In G. A. Watson, editor, Numerical Analysis, Lecture Notes in Mathematics 630. SpringerVerlag, Berlin.

(13)
Moré, J. J., B. S. Garbow, and K. E. Hillstrom, 1980. User guide for MINPACK1. Technical Report ANL8074, Argonne National Laboratory, Argonne, IL.

(14)
Sheng, J., E. Malkiel, and J. Katz, 2006. Digital holographic microscope for measuring threedimensional particle distributions and motions. Appl. Opt. 45:3893–3901.

(15)
Wiscombe, W. J., 1980. Improved Mie scattering algorithms. Appl. Opt. 19:1505–1509.