Physics:Dynamic light scattering

From HandWiki
Short description: Technique for determining size distribution of particles
Hypothetical dynamic light scattering of two samples: Larger particles on the top and smaller particles on the bottom

Dynamic light scattering (DLS) is a technique in physics that can be used to determine the size distribution profile of small particles in suspension or polymers in solution.[1] In the scope of DLS, temporal fluctuations are usually analyzed using the intensity or photon autocorrelation function (also known as photon correlation spectroscopyPCS or quasi-elastic light scatteringQELS). In the time domain analysis, the autocorrelation function (ACF) usually decays starting from zero delay time, and faster dynamics due to smaller particles lead to faster decorrelation of scattered intensity trace. It has been shown that the intensity ACF is the Fourier transform of the power spectrum, and therefore the DLS measurements can be equally well performed in the spectral domain.[2][3] DLS can also be used to probe the behavior of complex fluids such as concentrated polymer solutions.

Setup

A monochromatic light source, usually a laser, is shot through a polarizer and into a sample. The scattered light then goes through a second polarizer where it is collected by a photomultiplier and the resulting image is projected onto a screen. This is known as a speckle pattern (Figure 1).[4]

Figure 1. Typical speckle pattern.

All of the molecules in the solution are being hit with the light and all of the molecules diffract the light in all directions. The diffracted light from all of the molecules can either interfere constructively (light regions) or destructively (dark regions). This process is repeated at short time intervals and the resulting set of speckle patterns is analyzed by an autocorrelator that compares the intensity of light at each spot over time. The polarizers can be set up in two geometrical configurations. One is a vertical/vertical (VV) geometry, where the second polarizer allows light through that is in the same direction as the primary polarizer. In vertical/horizontal (VH) geometry the second polarizer allows light that is not in the same direction as the incident light.

Description

When light hits small particles, the light scatters in all directions (Rayleigh scattering) as long as the particles are small compared to the wavelength (below 250 nm). Even if the light source is a laser, and thus is monochromatic and coherent, the scattering intensity fluctuates over time. This fluctuation is due to small particles in suspension undergoing Brownian motion, and so the distance between the scatterers in the solution is constantly changing with time. This scattered light then undergoes either constructive or destructive interference by the surrounding particles, and within this intensity fluctuation, information is contained about the time scale of movement of the scatterers. Sample preparation either by filtration or centrifugation is critical to remove dust and artifacts from the solution.

The dynamic information of the particles is derived from the autocorrelation of the intensity trace recorded during the experiment. The second order autocorrelation curve is generated from the intensity trace as follows:

[math]\displaystyle{ g^2(q;\tau) = \frac{\langle I(t)I(t+\tau)\rangle}{\langle I(t)\rangle^2} }[/math]

where g2(q;τ) is the autocorrelation function at a particular wave vector, q, and delay time, τ, and I is the intensity. The angular brackets [math]\displaystyle{ \langle\cdot\rangle }[/math] denote the expected value operator, which in some texts is denoted by a capital E.

At short time delays, the correlation is high because the particles do not have a chance to move to a great extent from the initial state that they were in. The two signals are thus essentially unchanged when compared after only a very short time interval. As the time delays become longer, the correlation decays exponentially, meaning that, after a long period has elapsed, there is no correlation between the scattered intensity of the initial and final states. This exponential decay is related to the motion of the particles, specifically to the diffusion coefficient. To fit the decay (i.e., the autocorrelation function), numerical methods are used, based on calculations of assumed distributions. If the sample is monodisperse (uniform) then the decay is simply a single exponential. The Siegert equation relates the second-order autocorrelation function with the first-order autocorrelation function g1(q;τ) as follows:

[math]\displaystyle{ g^2(q;\tau)= 1+\beta\left[g^1(q;\tau)\right]^2 }[/math]

where the first term of the sum is related to the baseline value (≈1) and the parameter β is a correction factor that depends on the geometry and alignment of the laser beam in the light scattering setup. It is roughly equal to the inverse of the number of speckle (see Speckle pattern) from which light is collected. A smaller focus of the laser beam yields a coarser speckle pattern, a lower number of speckle on the detector, and thus a larger second-order autocorrelation. The most important use of the autocorrelation function is its use for size determination.

Multiple scattering

Dynamic light scattering provides insight into the dynamic properties of soft materials by measuring single scattering events, meaning that each detected photon has been scattered by the sample exactly once. In principle, the DLS measurements can be performed with the detector positioned at any angle. The choice of the best angle depends on the sample properties, such as turbidity and particle size.[5] Back scattering detection (e.g., 173° or 175°) is particularly interesting for turbid and highly concentrated samples, which contain large particles. Side scattering detection (90°) is recommended for weakly scattering samples, including small particles and transparent samples. Finally, forward scattering detection (e.g., 13° or 15°) is suitable for detection of samples containing small particles with few large particles. Some DLS instruments in the market also allow automatic angle selection based on a continuous transmittance measurement.

At the lower end of the turbidity range, the Cavity Amplified Scattering Spectroscopy method[6] makes use of an integrating cavity to elongate photon paths through quasi non-scattering samples. As opposed to conventional DLS instruments, this method is angle independent as it probes samples isotropically from all directions.

Even though the DLS measurement using a single-angle detection has been the most diffuse technique, the application to many systems of scientific and industrial relevance has been limited due to often-encountered multiple scattering, wherein photons are scattered multiple times by the sample before being detected. Accurate interpretation becomes exceedingly difficult for systems with non-negligible contributions from multiple scattering. Especially for larger particles and those with high refractive index contrast, this limits the technique to very low particle concentrations, and a large variety of systems are, therefore, excluded from investigations with dynamic light scattering. However, as shown by Schaetzel,Cite error: Closing </ref> missing for <ref> tag[7][8] Therefore, previous refractive index data from the scattering medium should be evaluated with dedicated instruments, known as refractometers. Alternatively, DLS instruments containing a refractive index measurement module allow a good estimative for this important parameter within ±0.5%, which is the accuracy defined by ISO 22412:2017[9][10] for refractive index values required for DLS. Besides the refractive index of the medium, the refractive index of the particles is only necessary when analyzing larger particle size (usually above 100 nm) and volume- or number-weighted size distributions are needed. In these cases, prior knowledge of the refractive index and absorbance of the material is required in order to apply the Mie scattering.[11][12]

Depending on the anisotropy and polydispersity of the system, a resulting plot of (Γ/q2) vs. q2 may or may not show an angular dependence. Small spherical particles will show no angular dependence, hence no anisotropy. A plot of (Γ/q2) vs. q2 will result in a horizontal line. Particles with a shape other than a sphere will show anisotropy and thus an angular dependence when plotting (Γ/q2) vs. q2.[13] The intercept will be in any case the Dt. Thus, there is an optimum angle of detection θ for each particle size. A high-quality analysis should always be performed at several scattering angles (multiangle DLS). This becomes even more important in a polydisperse sample with an unknown particle size distribution. At certain angles the scattering intensity of some particles will completely overwhelm the weak scattering signal of other particles, thus making them invisible to the data analysis at this angle. DLS instruments which only work at a fixed angle can only deliver good results for some particles. Thus, the indicated precision of a DLS instrument with only one detection angle is only ever true for certain particles.

Dt is often used to calculate the hydrodynamic radius of a sphere through the Stokes–Einstein equation. It is important to note that the size determined by dynamic light scattering is the size of a sphere that moves in the same manner as the scatterer. So, for example, if the scatterer is a random coil polymer, the determined size is not the same as the radius of gyration determined by static light scattering. It is also useful to point out that the obtained size will include any other molecules or solvent molecules that move with the particle. So, for example, colloidal gold with a layer of surfactant will appear larger by dynamic light scattering (which includes the surfactant layer) than by transmission electron microscopy (which does not "see" the layer due to poor contrast).

In most cases, samples are polydisperse. Thus, the autocorrelation function is a sum of the exponential decays corresponding to each of the species in the population.

[math]\displaystyle{ g^1(q;\tau)= \sum_{i=1}^n G_i(\Gamma_i)\exp(-\Gamma_i\tau) = \int G(\Gamma)\exp(-\Gamma\tau)\,d\Gamma. }[/math]

It is tempting to obtain data for g1(q;τ) and attempt to invert the above to extract G(Γ). Since G(Γ) is proportional to the relative scattering from each species, it contains information on the distribution of sizes. However, this is known as an ill-posed problem. The methods described below (and others) have been developed to extract as much useful information as possible from an autocorrelation function.

Cumulant method

One of the most common methods is the cumulant method,[14][15] from which in addition to the sum of the exponentials above, more information can be derived about the variance of the system as follows:

[math]\displaystyle{ \ g^1(q;\tau) = \exp\left(-\bar{\Gamma}\left(\tau - \frac{\mu_2}{2!}\tau^2 + \frac{\mu_3}{3!}\tau^3 + \cdots\right)\right) }[/math]

where Γ is the average decay rate and μ2/Γ2 is the second order polydispersity index (or an indication of the variance). A third-order polydispersity index may also be derived but this is necessary only if the particles of the system are highly polydisperse. The z-averaged translational diffusion coefficient Dz may be derived at a single angle or at a range of angles depending on the wave vector q.

[math]\displaystyle{ \ \bar{\Gamma}=q^2D_z\, }[/math]

One must note that the cumulant method is valid for small τ and sufficiently narrow G(Γ).[16] One should seldom use parameters beyond µ3, because overfitting data with many parameters in a power-series expansion will render all the parameters, including [math]\displaystyle{ \scriptstyle \bar{\Gamma} }[/math] and µ2, less precise.[17] The cumulant method is far less affected by experimental noise than the methods below.

Size-distribution function

The particle size distribution can also be obtained using the autocorrelation function. However, polydisperse samples are not well resolved by the cumulant fit analysis. Thus, the combination of non-negative least squares (NNLS) algorithms with regularization methods, such as the Tikhonov regularization, can be used to resolve multimodal samples.[11] An important feature of the NNLS optimization is the regularization term used to identify specific solutions and minimize the deviation between the measure data and the fit. There is no ideal regularization term that is suitable for all samples. The shape of this term can determine if the solution will represent a general broad distribution with small number of peaks or if narrow and discrete populations will be fit. Alternatively, the calculation of the particle size distribution is performed using the CONTIN algorithm.

CONTIN algorithm

An alternative method for analyzing the autocorrelation function can be achieved through an inverse Laplace transform known as CONTIN developed by Steven Provencher.[18][19] The CONTIN analysis is ideal for heterodisperse, polydisperse, and multimodal systems that cannot be resolved with the cumulant method. The resolution for separating two different particle populations is approximately a factor of five or higher and the difference in relative intensities between two different populations should be less than 1:10−5.

Maximum entropy method

The Maximum entropy method is an analysis method that has great developmental potential. The method is also used for the quantification of sedimentation velocity data from analytical ultracentrifugation. The maximum entropy method involves several iterative steps to minimize the deviation of the fitted data from the experimental data and subsequently reduce the χ2 of the fitted data.

Scattering of non-spherical particles

If the particle in question is not spherical, the rotational motion must be considered as well because the scattering of the light will be different depending on orientation. According to Pecora, rotational Brownian motion will affect the scattering when a particle fulfills two conditions; they must be both optically and geometrically anisotropic.[20] Rod-shaped molecules fulfill these requirements, so a rotational diffusion coefficient must be considered in addition to a translational diffusion coefficient. In its most succinct form, the equation appears as

[math]\displaystyle{ \frac{A}{B} = \frac{5}{4} \frac{4\Mu_p+2\Nu\Mu_l\Mu_p+\Mu_l}{\Mu_p-\Nu+\Mu_l} }[/math]

Where A/B is the ratio of the two relaxation modes (translational and rotational), Mp contains information about the axis perpendicular to the central axis of the particle, and Ml contains information about the axis parallel to the central axis.

In 2007, Peter R. Lang and his team decided to use dynamic light scattering to determine the particle length and aspect ratio of short gold nanorods.[21] They chose this method since it does not destroy the sample and it has a relatively easy setup. Both relaxation states were observed in VV geometry and the diffusion coefficients of both motions were used to calculate the aspect ratios of the gold nanoparticles.

Applications

DLS is used to characterize the size of various particles including proteins,[22] polymers, micelles,[23] Protein cages and virus-like particles,[24][25] vesicles,[26] carbohydrates, nanoparticles, biological cells,[27] and gels.[28] If the system is not disperse in size, the mean effective diameter of the particles can be determined. This measurement depends on the size of the particle core, the size of surface structures, particle concentration, and the type of ions in the medium.

Since DLS essentially measures fluctuations in scattered light intensity due to diffusing particles, the diffusion coefficient of the particles can be determined. DLS software of commercial instruments typically displays the particle population at different diameters. If the system is monodisperse, there should only be one population, whereas a polydisperse system would show multiple particle populations. If there is more than one size population present in a sample then either the CONTIN analysis should be applied for photon correlation spectroscopy instruments, or the power spectrum method should be applied for Doppler shift instruments.

Stability studies can be done conveniently using DLS. Periodical DLS measurements of a sample can show whether the particles aggregate over time by seeing whether the hydrodynamic radius of the particle increases. If particles aggregate, there will be a larger population of particles with a larger radius. In some DLS machines, stability depending on temperature can be analyzed by controlling the temperature in situ.

See also


References

  1. Berne, B.J.; Pecora, R. Dynamic Light Scattering. Courier Dover Publications (2000) ISBN:0-486-41155-9
  2. Chu, B. (1 January 1970). "Laser Light Scattering". Annual Review of Physical Chemistry 21 (1): 145–174. doi:10.1146/annurev.pc.21.100170.001045. Bibcode1970ARPC...21..145C. 
  3. Pecora., R. (1964). "Doppler Shifts in Light Scattering from Pure Liquids and Polymer Solutions". The Journal of Chemical Physics 40 (6): 1604. doi:10.1063/1.1725368. Bibcode1964JChPh..40.1604P. 
  4. Goodman, J (1976). "Some fundamental properties of speckle". J. Opt. Soc. Am. 66 (11): 1145–1150. doi:10.1364/josa.66.001145. Bibcode1976JOSA...66.1145G. 
  5. "Multiple detection angles in dynamic light scattering analysis :: Anton Paar Wiki" (in de). https://wiki.anton-paar.com/at-de/dynamische-lichtstreuung-mit-mehreren-messwinkeln/. 
  6. Graciani, Guillaume; King, John T.; Amblard, François (2022-08-30). "Cavity-Amplified Scattering Spectroscopy Reveals the Dynamics of Proteins and Nanoparticles in Quasi-transparent and Miniature Samples" (in en). ACS Nano 16 (10): 16796–16805. doi:10.1021/acsnano.2c06471. ISSN 1936-0851. PMID 36039927. https://pubs.acs.org/doi/10.1021/acsnano.2c06471. 
  7. Strutt, J.W. (1871-04-01). "XXXVI. On the light from the sky, its polarization and colour". The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 41 (273): 274–279. doi:10.1080/14786447108640479. ISSN 1941-5982. https://doi.org/10.1080/14786447108640479. 
  8. Stetefeld, Jörg; McKenna, Sean A.; Patel, Trushar R. (2016-12-01). "Dynamic light scattering: a practical guide and applications in biomedical sciences" (in en). Biophysical Reviews 8 (4): 409–427. doi:10.1007/s12551-016-0218-6. ISSN 1867-2469. PMID 28510011. PMC 5425802. https://doi.org/10.1007/s12551-016-0218-6. 
  9. "ISO 22412:2017. Particle size analysis — Dynamic light scattering (DLS)" (in en). https://www.iso.org/cms/render/live/en/sites/isoorg/contents/data/standard/06/54/65410.html. 
  10. Tareq, Syed Mohammed; Boutchuen, Armel; Roy, Shuvashish; Zimmerman, Dell; Jur, Gitapun; Bathi, Jejal Reddy; Palchoudhury, Soubantika (August 2021). "A Dynamic Light Scattering Approach for Detection of Nanomaterials in Tennessee River" (in en). Water Resources Research 57 (8). doi:10.1029/2020WR028687. ISSN 0043-1397. Bibcode2021WRR....5728687T. https://onlinelibrary.wiley.com/doi/10.1029/2020WR028687. 
  11. 11.0 11.1 "White Paper: Understanding the Concept of Dynamic Light Scattering :: Anton-Paar.com" (in en). https://www.anton-paar.com/corp-en/services-support/document-finder/application-reports/white-paper-understanding-the-concept-of-dynamic-light-scattering/. 
  12. Alfano, Brigida; Barretta, Luigi; Del Giudice, Antonio; De Vito, Saverio; Di Francia, Girolamo; Esposito, Elena; Formisano, Fabrizio; Massera, Ettore et al. (2020-11-29). "A Review of Low-Cost Particulate Matter Sensors from the Developers' Perspectives" (in en). Sensors 20 (23): 6819. doi:10.3390/s20236819. ISSN 1424-8220. PMID 33260320. Bibcode2020Senso..20.6819A. 
  13. Gohy, Jean-François; Varshney, Sunil K.; Jérôme, Robert (2001). "Water-Soluble Complexes Formed by Poly(2-vinylpyridinium)-block-poly(ethylene oxide) and Poly(sodium methacrylate)-block-poly(ethylene oxide) Copolymers". Macromolecules 34 (10): 3361. doi:10.1021/ma0020483. Bibcode2001MaMol..34.3361G. 
  14. Koppel, Dennis E. (1972). "Analysis of Macromolecular Polydispersity in Intensity Correlation Spectroscopy: The Method of Cumulants". The Journal of Chemical Physics 57 (11): 4814–4820. doi:10.1063/1.1678153. Bibcode1972JChPh..57.4814K. 
  15. Frisken, Barbara J. (2001). "Revisiting the Method of Cumulants for the Analysis of Dynamic Light-Scattering Data". Applied Optics 40 (24): 4087–91. doi:10.1364/AO.40.004087. PMID 18360445. Bibcode2001ApOpt..40.4087F. https://www.sfu.ca/biophysics/publications/2001/ApplOpt2001cumulants.pdf. 
  16. Hassan, Pa; Kulshreshtha, Sk (August 2006). "Modification to the cumulant analysis of polydispersity in quasielastic light scattering data". Journal of Colloid and Interface Science 300 (2): 744–8. doi:10.1016/j.jcis.2006.04.013. ISSN 0021-9797. PMID 16790246. Bibcode2006JCIS..300..744H. 
  17. Chu, B (1992). Laser Light scattering: Basic Principles and Practice. Academic Press. ISBN 978-0-12-174551-6. 
  18. Provencher, S (1982). "CONTIN: A general purpose constrained regularization program for inverting noisy linear algebraic and integral equations". Computer Physics Communications 27 (3): 229–242. doi:10.1016/0010-4655(82)90174-6. Bibcode1982CoPhC..27..229P. http://s-provencher.com/pub/contin/cpc2.pdf. 
  19. Provencher, S. W. (1982). "A constrained regularization method for inverting data represented by linear algebraic or integral equations". Comput. Phys. Commun. 27 (3): 213–227. doi:10.1016/0010-4655(82)90173-4. Bibcode1982CoPhC..27..213P. http://s-provencher.com/pub/contin/cpc1.pdf. 
  20. Aragón, S. R.; Pecora, R. (1976). "Theory of dynamic light scattering from polydisperse systems". The Journal of Chemical Physics 64 (6): 2395. doi:10.1063/1.432528. Bibcode1976JChPh..64.2395A. 
  21. Rodríguez-Fernández, J.; Pérez−Juste, J.; Liz−Marzán, L. M.; Lang, P. R. (2007). "Dynamic Light Scattering of Short Au Rods with Low Aspect Ratios". The Journal of Physical Chemistry 111 (13): 5020–5025. doi:10.1021/jp067049x. http://juser.fz-juelich.de/record/56159/files/DLS_aurods_revised.pdf. 
  22. Dolinska, Monika B.; Young, Kenneth L.; Kassouf, Claudia; Dimitriadis, Emilios K.; Wingfield, Paul T.; Sergeev, Yuri V. (2020-01-03). "Protein Stability and Functional Characterization of Intra-Melanosomal Domain of Human Recombinant Tyrosinase-Related Protein 1" (in en). International Journal of Molecular Sciences 21 (1): 331. doi:10.3390/ijms21010331. ISSN 1422-0067. PMID 31947795. 
  23. Bhut, Parth Rajeshkumar; Pal, Nilanjan; Mandal, Ajay (2019-12-03). "Characterization of Hydrophobically Modified Polyacrylamide in Mixed Polymer-Gemini Surfactant Systems for Enhanced Oil Recovery Application" (in en). ACS Omega 4 (23): 20164–20177. doi:10.1021/acsomega.9b02279. ISSN 2470-1343. PMID 31815217. 
  24. "Cytochrome C with peroxidase-like activity encapsulated inside the small DPS protein nanocage". Journal of Materials Chemistry B 9 (14): 3168–3179. March 2021. doi:10.1039/d1tb00234a. PMID 33885621. 
  25. "Virus-Like Particles (VLPs) as a Platform for HierarchicalCompartmentalization". Biomacromolecules 21 (6): 2060–2072. April 2020. doi:10.1021/acs.biomac.0c00030. PMID 32319761. 
  26. Velu, Sabareesh K. P.; Yan, Minhao; Tseng, Kuo-Pi; Wong, Ken-Tsung; Bassani, Dario M.; Terech, Pierre (6 February 2013). "Spontaneous Formation of Artificial Vesicles in Organic Media through Hydrogen-Bonding Interactions". Macromolecules 46 (4): 1591–1598. doi:10.1021/ma302595g. Bibcode2013MaMol..46.1591V. 
  27. Jena, Sidhartha S.; Joshi, Hiren M.; Sabareesh, K.P.V.; Tata, B.V.R.; Rao, T.S. (2006). "Dynamics of Deinococcus radiodurans under Controlled Growth Conditions". Biophysical Journal 91 (7): 2699–2707. doi:10.1529/biophysj.106.086520. PMID 16829564. Bibcode2006BpJ....91.2699J. 
  28. Sabareesh, K. P. V.; Jena, Sidhartha S.; Tata, B. V. R. (5 May 2006). "Dynamic Light Scattering Studies on Photo Polymerized and Chemically Cross‐linked Polyacrylamide Hydrogels". AIP Conference Proceedings 832 (1): 307–310. doi:10.1063/1.2204513. ISSN 0094-243X. Bibcode2006AIPC..832..307S. 

External links