Holographic characterization of colloidal fractal aggregates
In-line holographic microscopy images of micrometer-scale fractal aggregates can be interpreted with an effective-sphere model to obtain each aggregate’s size and the population-averaged fractal dimension. We demonstrate this technique experimentally using model fractal clusters of polystyrene nanoparticles and fractal protein aggregates composed of bovine serum albumin and bovine pancreas insulin.
Holograms of micrometer-scale colloidal spheres obtained with in-line holographic video microscopy can be analyzed with the Lorenz-Mie theory of light scattering to obtain each sphere’s radius and refractive index, typically with part-per-thousand precision [1, 2]. Characterizing a single sphere requires a few milliseconds on a standard computer , which is fast enough that several thousand spheres can be characterized in under ten minutes. Analyzing holograms of aspherical particles and colloidal clusters is substantially more challenging , particularly if no information is available a priori about the particles’ geometry. We previously have demonstrated that the Lorenz-Mie analysis developed for homogeneous spheres also yields useful characterization data for porous spheres [5, 6], dimpled spheres , and protein aggregates . The last of these applications treats each irregularly-shaped protein aggregate as an effective sphere composed of the aggregate itself and the fluid medium filling its pores. Here, we develop an effective-sphere formalism for Lorenz-Mie microscopy of fractal aggregates and demonstrate its efficacy through measurements on model systems.
Detecting and characterizing micrometer-scale aggregates is useful both for fundamental research and also for solving real-world problems. Protein aggregation, for example, is a critical concern for the biopharmaceutical industry because it limits the efficacy of protein-based drugs and can induce harmful immunogenic responses in patients . Information on the concentration, size distribution and morphology of protein aggregates provides guidance for formulating stable products and for avoiding adverse clinical outcomes. Conventional light-scattering techniques do not work well for particles in the relevant size range , and cannot distinguish aggregates of interest from other contaminants commonly found in commercial formulations. Similar detection and characterization challenges arise in the precision slurries used by the semiconductor manufacturing for chemical-mechanical planarization . As a particle-resolved measurement technique, holographic characterization naturally differentiates micrometer-scale particles by size and composition . The effective-sphere model extends these capabilities to include assessment of particle morphology without sacrificing speed or ease of use.
II Lorenz-Mie characterization
Lorenz-Mie characterization, depicted schematically in Figure 1(a), is based on in-line holographic video microscopy , in which the sample is illuminated with a collimated laser beam. Light scattered by a particle interferes with the remainder of the illumination in the focal plane of a microscope. The intensity of the magnified interference pattern is recorded with a conventional video camera for analysis. A typical example is shown in Figure 1(b). Each holographic snapshot is corrected for the camera’s dark count, , and is normalized by the background intensity in the field of view, , to obtain [13, 1, 2]
The normalized hologram then is fit to the prediction 
where is the wavenumber of light in the medium, is the position of the particle’s center relative to the center of the microscope’s focal plane, and is the Lorenz-Mie function that describes scattering of the incident wave by a sphere of radius and refractive index [14, 15]. The form of Eq. (2) is appropriate for a beam propagating along that is linearly polarized along . Fitting Eq. (2), pixel by pixel, to the normalized hologram of a sphere yields the sphere’s three-dimensional position, its radius and its refractive index . Figure 1(c) presents the result of fitting to the experimental hologram from Figure 1(b).
The custom-built holographic microscope used for this study illuminates the sample with the collimated beam from a solid state laser (Coherent Cube) operating at a vacuum wavelength of . The sample flows through the observation volume in microfluidic channel created by bonding the edges of a number 1.5 glass microscope cover slip to the face of a standard glass microscope slide. This sample cell is mounted on a translation stage (Prior, Proscan II) in the focal plane of an oil-immersion objective lens (Nikon, Plan Apo, , numerical aperture 1.45). Light collected by the objective lens is relayed by an achromatic tube lens to a video camera (NEC, TI-324AII), which records its intensity 30 times per second with an effective magnification of . Each video frame is a measurement of with a resolution of . Features associated with dispersed particles are identified in digitized holographic images  and are analyzed with Eqs. (1) and (2) using methods that have been described in detail elsewhere [3, 2]. The example in Figure 1(b) is a region of interest cropped from the normalized hologram, , obtained from according to Eq. (1).
All principal results were reproduced using a second instrument based on a air objective (Nikon Plan Fluor, numerical aperture 0.75) operating at a vacuum wavelength of (Thorlabs, CPS532 ) with an effective magnification of on an Allied Vision Mako U-130B camera. This camera yields images with . Samples flow through this instrument in prefabricated microfluidic channels with path length (Ibidi, Slide VI, uncoated). Flow is driven by a syringe pump (New Era Systems, NE 100). This reduces the possibility that instrumental artifacts might have influenced the scaling relationships reported here.
A colloidal sample is characterized by placing a aliquot in the reservoir at one end of the microfluidic channel and drawing it through with a small pressure gradient. The resulting Poiseuille flow has a peak speed along its axis of , which is small enough to avoid artifacts due to motion blurring [17, 18] and allows each particle to be recorded several times during its transit. Given a concentration on the order of , a few thousand particles will pass through the observation volume in .
A single snapshot of an individual colloidal particle can be analyzed in several milliseconds using standard computer hardware [17, 3]. Characterization data therefore can be acquired in real time as particles flow down the microfluidic channel . The images in Figure 2(a) show eight stages of the transit of a typical polystyrene aggregate at intervals. The resulting time series of position and characterization data can be linked into a trajectory using maximum likelihood methods [19, 17]. The scatter plot in Figure 2(b) shows the estimated values for the radius and refractive index obtained at each stage of the aggregate’s trajectory, recorded at intervals. These results can be combined into a trajectory-averaged estimate for the associated particle’s characteristics.
When this measurement technique is applied to spherical particles, the standard deviation of the trajectory-averaged characteristics is comparable to the single-measurement precision [17, 2]. Results for the irregularly-shaped aggregate in Figure 2 vary more substantially. The standard deviation of the radius, is a factor of 10 larger than the single-measurement precision. The standard deviation of the refractive index, by contrast, is comparable to the single particle precision, . It is possible that the variation in apparent size occurs because the aggregate is irregularly shaped, tumbles as it travels, and so presents different-sized projections to the instrument. We develop this interpretation in the next section.
III Effective sphere model
Rather than attempting to generalize Eq. (2) to account for the detailed structure of random colloidal aggregates, we instead analyze their holograms with Eq. (2) itself and interpret the results with effective medium theory [20, 21]. The radius, , and refractive index, , obtained from such a fit then characterize an effective sphere enclosing both the fractal aggregate and also the intercalated fluid medium. The goal of the present work is to establish a relationship between the measured properties of the effective sphere and the underlying properties of the actual aggregate. This relationship then enables Lorenz-Mie microscopy to probe the morphology of fractal aggregates without incurring the computational burden of detailed modeling.
For simplicity, we assume that a fractal aggregate of fractal dimension is composed of identical spherical monomers, each of radius and refractive index . The aggregate is immersed in a medium of refractive index that fills the pores. Provided that both the monomers and the pores are substantially smaller than the wavelength of light, this composite structure may be modeled as a continuous medium whose refractive index, , is given by the Maxwell Garnett relation ,
|where is the volume fraction of monomers in the sphere and|
|is the Lorentz-Lorenz factor.|
The number of monomers within radius of a fractal aggregate’s center is
and the associated volume fraction is
For a particle of radius , the overall volume fraction is
From this and the Maxwell Garnett relation, we obtain an expression for the cluster’s effective refractive index,
This result also may be expressed as a scaling relationship between the radius of a fractal aggregate and its effective refractive index,
In a population of aggregates grown under comparable conditions, Eq. (8) can be used to estimate the population-averaged fractal dimension, .
Equation (7) treats a cluster as a homogeneous medium. In fact, the density of monomers decreases with scale, and therefore with radius within the cluster. As discussed in the Appendix, the associated radial gradient in the refractive index has little influence on practical measurements of clusters’ effective characteristics. We, therefore, ignore the clusters’ spatial inhomogeneity in the discussion that follows.
IV Experimental studies of model fractal aggregates
IV.1 Fractal aggregates of polystyrene nanospheres
The data in Figure 3 were obtained for colloidal fractal aggregates grown under conditions conducive to diffusion-limited cluster aggregation (DLCA) [22, 23, 24]. The primary particles are monodisperse polystyrene spheres with a mean radius of (Thermo Scientific, catalog number 5016A, 10% w/w).
Aggregation was initiated by dispersing these particles in NaCl solution at a concentration of by weight. After one hour, the dispersion was diluted by a factor of with deionized water to stop further aggregation. The resulting sample was then analyzed immediately, before the aggregates might have time to restructure [25, 24].
Each data point in Figure 3(a) represents the characteristics of a single colloidal particle obtained from a trajectory such as the example in Figure 2. The plot symbols’ size is comparable to the numerical uncertainty in the fit parameters. Each point is colored according to the relative probability density, , of measurements according to the color bar inset into Figure 3(a).
The effective sphere model works well for aggregates with radii smaller than . Despite their irregular shapes, these objects are small enough that their holograms display the radial symmetry typical of spheres. This can be seen in Figure 3(b). The reduced statistics for fits such as Figure 3(c) typically fall within ten percent of unity, suggesting that the fit parameters reliably reflect the aggregates properties.
Larger aggregates are reliably detected and counted by the feature identification algorithm , but are poorly characterized by the effective-sphere model . Their holograms are more substantially asymmetric, and the reduced statistic for these fits typically exceeds . It is not surprising, therefore, that the estimated characteristics for aggregates with do not follow the trend expected for fractal aggregates, their refractive indexes falling below the scaling prediction.
The solid curve in Figure 3(a) shows the prediction from Eq. (8) for , with no adjustable parameters. In addition to the monomers’ radius, the effective-sphere model is parametrized by the monomers’ refractive index and the clusters’ fractal dimension. The former, , was obtained from holographic characterization studies of emulsion-polymerized polystyrene spheres [1, 17, 26, 5]. The latter, , was obtained from independent light-scattering studies on aggregates grown under comparable conditions [25, 23, 24, 27, 28], and is consistent with expectations for DLCA [29, 27].
Figure 3(b) shows the same data replotted to emphasize the scaling prediction from Eq. (8). Agreement with the effective sphere model is quite good for particles with apparent radii smaller than . Because fractals’ pore size increases with scale, larger aggregates presumably do not satisfy the requirements of effective medium theory and so are not so well described by the effective-sphere model.
The scanning electron microscope image inset into Figure 3(b) shows a typical vacuum-dried aggregate. This image resolves the individual spheres, whose arrangement is consistent with the irregular branched structure inferred from holographic characterization of similar samples. Although details of the structure undoubtedly were altered during drying, the presence of voids at multiple scales within the cluster is consistent with a fractal dimension smaller than .
The model’s success for smaller clusters supports the contention that Eq. (8) can be useful for measuring the population-averaged fractal dimension of micrometer-scale fractal clusters. We next apply this approach to characterize aggregates of two model proteins whose cluster morphologies have been independently established. This application not only serves to verify the effective-sphere model for fractal clusters, but also illustrates the utility of Lorenz-Mie microscopy for measuring the size distribution and morphology of protein aggregates, a subject of considerable interest in biology  and of substantial practical importance in pharmaceutical manufacturing .
IV.2 Protein aggregates
Like the colloidal nanoparticles considered in the previous section, proteins in solution also have a tendency to aggregate . Some of the resulting macromolecular structures perform important biological functions. Others cause diseases. Protein aggregation is a principal failure mechanism for biopharmaceutical formulations not only because clustered proteins are less effective as therapeutic agents, but also because they can elicit dangerous immune responses . We previously have demonstrated that inline holographic video microscopy can detect protein aggregates in solution, and can distinguish them from such common contaminants as silicone oil droplets . Here, we apply the effective-sphere model to study the aggregates’ morphology by estimating their fractal dimension.
IV.2.1 Bovine insulin
The data in Figure 4(a) were obtained for aggregates of bovine pancreas insulin (Mw: , Sigma-Aldrich, CAS number: 11070-73-8) that were prepared according to previously published methods [32, 33]. Insulin was dissolved at a concentration of in Tris buffer (Life Technologies, CAS number 77-86-1). The pH of the buffer was adjusted to 7.4 with hydrochloric acid (Sigma Aldrich, CAS number: 7647-01-0). The solution then was centrifuged at for to induce aggregation, at which time the sample still appeared transparent to visual inspection.
Running of this sample through the holographic characterization instrument reveals a concentration of , including the examples whose properties are plotted in Figure 4(a). Because the effective monomer radius is not known a priori, radii in Figure 4 are scaled by an arbitrary factor, . This choice does not affect the estimate for .
The main distribution of single-particle characteristics follows the scaling prediction quite well, and has a slope consistent with a fractal dimension of . This is denoted in Figure 4(a) by a solid (red) line superimposed on the data.
Dashed lines in this plot show equivalent results for fractal dimensions and . The proposal that bovine insulin forms branched fractal aggregates is consistent with independent measurements of such aggregates’ morphology using atomic-force microscopy .
In addition to the main distribution of points, Figure 4(a) features an outlying cluster of large-size aggregates comparable to those in Figure 3. It also includes a cluster of small particles with low refractive indexes. These latter features appear to correspond to globular aggregates that are distinct from the fractal clusters of interest here.
IV.2.2 Bovine serum albumin
Solutions of bovine serum albumin (BSA) (Mw: , Sigma Aldrich, CAS number: 9048-46-8) were aggregated by complexation with poly(allylamine hydrochloride) (PAH) (Mw: , CAS number: 71550-12-4, average degree of polymerization: ) [35, 36]. BSA and PAH were dissolved in Tris-HCl buffer (pH 7.4) (Life Technologies, CAS number: 77-86-1) to achieve concentrations of and , respectively. The reagents were mixed by vortexing to ensure dissolution, and aggregates formed after the sample was allowed to equilibrate for one hour.
Holographic characterization data for a sample of BSA-PAH complexes, shown in Figure 4(b), reveal in the range of radii running from to , and a peak radius of . Although our implementation of holographic characterization is capable of detecting aggregates with radii up to , no aggregates were observed with radii exceeding . The results for BSA aggregates with radii smaller than agree well with the scaling prediction from Eq. (8), this time with an apparent fractal dimension of , as indicated by the solid (red) line superimposed on the data in Figure 4(b). The nearly linear structure suggested by this low fractal dimension is consistent with atomic force microscopy images of BSA aggregates .
The experimental results presented in Sec. IV demonstrate that Lorenz-Mie microscopy can provide useful insights into the properties of micrometer-scale fractal aggregates. Holograms of micrometer-scale colloidal fractal aggregates can be interpreted with the effective-sphere model presented in Sec. III to estimate an aggregate’s size and effective refractive index. These particle-resolved data, in turn, can be pooled to estimate the population-averaged fractal dimension. The effective-sphere model therefore extends the particle-characterization capabilities of Lorenz-Mie microscopy to irregularly branched objects.
The success of the effective-sphere model in characterizing ramified protein aggregates lends additional support to the earlier proposal  that Lorenz-Mie characterization meaningfully assesses such aggregates’ sizes. It therefore establishes Lorenz-Mie microscopy as a method for sizing protein aggregates, characterizing their morphology, and differentiating them from other types of colloidal particles.
This work also provides a baseline against which more detailed approaches to holographic characterization of fractal structures may be compared. Future extensions based on machine-learning techniques  or direct modeling of the spatial distribution of dielectric material  thus can be tested directly using the methods described here.
This work was funded primarily by the National Institutes of Health under Grant Number 1R43TR001590. Additional support was provided by the MRSEC program of the National Science Foundation under Award Number DMR-1420073. The holographic characterization instrument used for this study was constructed with support of the MRI program of the NSF under Award Number DMR-0922680. The scanning electron microscope was purchased with financial support from the MRI program of the NSF under Award DMR-0923251. We are grateful to Prof. Andrew Hollingsworth for assistance with the SEM.
Appendix: Effective-sphere model for gradient-index spheres
The effective-sphere model treats an aggregate as if the monomers were distributed uniformly within it. In fact, the marginal volume fraction decreases with distance from the center of the aggregate as
This corresponds to a radial variation of the effective refractive index described by
Unless the entire aggregate is smaller than the wavelength of light, this radial structure might be expected to influence results obtained by applying effective medium theory to holograms of fractal aggregates.
To assess this influence, we use the effective-sphere model to analyze synthetic holograms of gradient-index particles with refractive-index profiles described by Eq. (10). These holograms are computed by replacing in Eq. (2) with the corresponding generalized Lorenz-Mie result for a stratified sphere [38, 39, 40] whose layers have refractive indexes given by Eq. (10). The number of layers is chosen to converge the computed intensities to within at each pixel. This typically occurs with layer thicknesses comparable to . A similar approach has proved successful for Luneburg spheres and other particles with continuous radial refractive index profiles . The hologram is then fit to the Lorenz-Mie model from Eq. (2) for the effective sphere’s radius, , and refractive index, . These parameters then can be compared with the true radius of the stratified sphere, , and the sphere-averaged refractive index, , obtained from Eq. (8).
Figure 5 shows the performance of the effective-sphere model for particles with , and in a medium with . These parameters are chosen to model fractal polystyrene aggregates dispersed in water [23, 24, 28]. The effective sphere’s radius, Figure 5(a), and refractive index, Figure 5(b), both track the true values, albeit with systematic offsets. They quite closely satisfy the anticipated scaling form predicted by Eq. (8) with a slope consistent with the input fractal dimension, as can be seen in Figure 5(c).
Comparable results are obtained with different values of the fractal dimension. The gradient-index structure of fractal aggregates therefore does not substantially diminish the ability of the effective-sphere model to estimate such particles’ fractal dimension. Although the gradient-index model does not address the influence of real aggregates’ branched structure, it lends additional confidence to the proposal that Lorenz-Mie characterization usefully assesses the properties of such objects.
-  S.-H. Lee, Y. Roichman, G.-R. Yi, S.-H. Kim, S.-M. Yang, A. van Blaaderen, P. van Oostrum and D. G. Grier. “Characterizing and tracking single colloidal particles with video holographic microscopy.” Opt. Express 15, 18275–18282 (2007).
-  B. J. Krishnatreya, A. Colen-Landy, P. Hasebe, B. A. Bell, J. R. Jones, A. Sunda-Meya and D. G. Grier. “Measuring Boltzmann’s constant through holographic video microscopy of a single sphere.” Am. J. Phys. 82, 23–31 (2014).
-  A. Yevick, M. Hannel and D. G. Grier. ‘‘Machine-learning approach to holographic particle characterization.” Opt. Express 22, 26884–26890 (2014).
-  J. Fung, R. W. Perry, T. G. Dimiduk and V. N. Manoharan. “Imaging multiple colloidal particles by fitting electromagnetic scattering solutions to digital holograms.” J. Quant. Spectr. Rad. Trans. 113, 2482–2489 (2012).
-  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).
-  C. Wang, H. W. Moyses and D. G. Grier. “Stimulus-responsive colloidal sensors with fast holographic readout.” Appl. Phys. Lett. 107, 051903 (2015).
-  M. Hannel, C. Middleton and D. G. Grier. “Holographic characterization of imperfect colloidal spheres.” Appl. Phys. Lett. 107, 141905 (2015).
-  C. Wang, X. Zhong, D. B. Ruffner, A. Stutt, L. A. Philips, M. D. Ward and D. G. Grier. “Holographic characterization of protein aggregates.” J. Pharm. Sci. 105, 1074–1085 (2016).
-  W. Wang. “Protein aggregation and its inhibition in biopharmaceutics.” Int. J. Pharma. 289, 1–30 (2005).
-  J. Panchal, J. Kotarek, E. Marszal and E. M. Topp. “Analyzing subvisible particles in protein drug products: a comparison of Dynamic Light Scattering (DLS) and Resonant Mass Measurement (RMM).” AAPS J. 16, 440–451 (2014).
-  G. B. Basim and B. M. Moudgil. “Effect of soft agglomerates on CMP slurry performance.” J. Colloid Interface Sci. 256, 137–142 (2002).
-  J. Sheng, E. Malkiel and J. Katz. “Digital holographic microscope for measuring three-dimensional particle distributions and motions.” Appl. Opt. 45, 3893–3901 (2006).
-  S.-H. Lee and D. G. Grier. “Holographic microscopy of holographically trapped three-dimensional structures.” Opt. Express 15, 1505–1512 (2007).
-  C. F. Bohren and D. R. Huffman. Absorption and Scattering of Light by Small Particles (Wiley Interscience, New York, 1983).
-  M. I. Mishchenko, L. D. Travis and A. A. Lacis. Scattering, Absorption and Emission of Light by Small Particles (Cambridge University Press, Cambridge, 2001).
-  B. J. Krishnatreya and D. G. Grier. ‘‘Fast feature identification for holographic tracking: The orientation alignment transform.” Opt. Express 22, 12773–12778 (2014).
-  F. C. Cheong, B. Sun, R. Dreyfus, J. Amato-Grill, K. Xiao, L. Dixon and D. G. Grier. “Flow visualization and flow cytometry with holographic video microscopy.” Opt. Express 17, 13071–13079 (2009).
-  L. Dixon, F. C. Cheong and D. G. Grier. “Holographic particle-streak velocimetry.” Opt. Express 19, 4393–4398 (2011).
-  J. C. Crocker and D. G. Grier. “Methods of digital video microscopy for colloidal studies.” J. Colloid Interface Sci. 179, 298–310 (1996).
-  D. E. Aspnes. “Local-field effects and effective-medium theory: A microscopic perspective.” Am. J. Phys. 50, 704–709 (1982).
-  T. C. Choy. Effective Medium Theory: Principles and Applications (Oxford University Press, New York, 1999).
-  P. Meakin. “Fractal aggregates.” Adv. Colloid Interface Sci. 28, 249–331 (1988).
-  D. Majolino, F. Mallamace, P. Migliardo, N. Micali and C. Vasi. “Elastic and quasielastic light-scattering studies of the aggregation phenomena in watersolutions of polystyrene particles.” Phys. Rev. A 40, 4665–4674 (1989).
-  Z. Zhou and B. Chu. “Fractal study of polystyrene latex and silica particle aggregates.” Physica A 177, 93–100 (1991).
-  C. Aubert and D. S. Cannell. “Restructuring of colloidal silica aggregates.” Phys. Rev. Lett. 56, 738–741 (1986).
-  K. Xiao and D. G. Grier. “Sorting colloidal particles into multiple channels with optical forces: Prismatic optical fractionation.” Phys. Rev. E 82, 051407 (2010).
-  C. M. Sorensen. “Light scattering by fractal aggregates: A review.” Aerosol Sci. Tech. 35, 648–687 (2001).
-  H. Wu, M. Lattuada and M. Morbidelli. “Dependence of fractal dimension of DLCA clusters on size of primary particles.” Adv. Colloid Interface Sci. 195-196, 41–49 (2013).
-  D. A. Weitz, J. S. Huang, M. Y. Lin and J. Sung. “Limits of the fractal dimension for irreversible kinetic aggregation of gold colloids.” Phys. Rev. Lett. 54, 1416–1419 (1985).
-  A. M. Morris, M. A. Watzky and R. G. Finke. “Protein aggregation kinetics, mechanism and curve fitting: A review of the literature.” Biochimica et Biophysica Acta 1794, 375–397 (2009).
-  W. Wang. “Instability, stabilization, and formulation of liquid protein pharmaceuticals.” Int. J. Pharm. 185, 129–188 (1999).
-  V. Sluzky, J. A. Tamada, A. M. Klibanov and R. Langer. “Kinetics of insulin aggregation in aqueous solutions upon agitation in the presence of hydrophobic surfaces.” Proc. Natl. Acad. Sci. U.S.A. 88, 9377–9381 (1991).
-  H. R. Costantino, R. Langer and A. M. Klibanov. “Moisture-induced aggregation of lyophilized insulin.” Pharm. Res. 11, 21–29 (1994).
-  K. Siposova, M. Kubovcikova, Z. Bednarikova, M. Koneracka, V. Savisova, A. Antosova, P. Kopcansky, Z. Daxnerova and Z. Gazova. “Depolymerization of insulin amyloid fibrils by albumin-modified magnetic field.” Nanotech. 23, 055101 (2012).
-  V. Ball, M. Winterhalter, P. Schwinte, P. Lavalle, J.-C. Voegel and P. Schaaf. “Complexation mechanism of bovine serum albumin and poly(allylamine hydrochloride).” J. Phys. Chem. B 106, 2357–2364 (2002).
-  T. Hagiwara, H. Kumagai and K. Nakamura. “Fractal analysis of aggregates formed by heating dilute BSA solutions using light scattering methods.” Biosci. Biotech. Biochem. 60, 1757–1763 (1996).
-  M. Omichi, A. Asano, S. Tsukuda, K. Takano, M. Sugimoto, A. Saeki, D. Sakamaki, A. Onoda, T. Hayashi and S. Seki. “Fabrication of enzyme-degradable and size-controlled protein nanowires using single particle nano-fabrication technique.” Nature Commun. 5, 3717 (2014).
-  W. Yang. “Improved recursive algorithm for light scattering by a multilayered sphere.” Appl. Opt. 42, 1710–1720 (2003).
-  O. Pena and U. Pal. “Scattering of electromagnetic radiation by a layered sphere.” Comp. Phys. Commun. 180, 2348–2354 (2009).
-  G. Gouesbet and G. Gréhan. Generalized Lorenz-Mie Theories (Springer-Verlag, Berlin, 2011).
-  M. Selmke. “Mie and potential scattering by a refractive inhomogeneity: Electromagnetic scattering by an infinite Coulomb-like scatterer.” J. Quant. Spectr. Rad. Trans. 162, 175–183 (2015).