Dynamic Light Scattering

Dynamic light scattering (also known as Photon Correlation Spectroscopy or Quasi-Elastic Light Scattering) is a technique in physics, which can be used to determine the size distribution profile of small particles in solution


When light hits small particles the light scatters in all directions (Rayleigh scattering) so long as the particles are small compared to the wavelength (< 250 nm). If the light source is a laser, and thus is monochromatic and coherent, then one observes a time-dependent fluctuation in the scattering intensity. These fluctuations are due to the fact that the small molecules in solutions are 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.

There are several ways to derive dynamic information about particles' movement in solution by Brownian motion. One such method is dynamic light scattering, also known as quasi-elastic laser light scattering. The dynamic information of the particles is derived from an autocorrelation of the intensity trace recorded during the experiment. The second order autocorrelation curve is generated from the intensity trace as follows:


where g2(q;τ) is the autocorrelation function at a particular wave vector, q, and delay time, τ, and I is the intensity. 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 starts to exponentially decay to zero, meaning that after a long time 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 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:


where 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. Once the autocorrelation data have been generated, different mathematical approaches can be employed to determine from it the z-averaged translational diffusion coefficient. The simplest approach is to treat the first order autocorrelation function as the sum of single exponential decays with fractions G(Γ), where Γ is the decay rate.


Data analysis

Cumulant method

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


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




where  is the incident laser wavelength,  is the refractive index of the sample and  is angle at which the detector is located with respect to the sample cell.

Depending on the anisotropy and polydispersity of the system, a resulting plot of  vs.  may or may not show an angular dependence. Spherical particles will show no angular dependence, hence no anisotropy. A plot of  vs.  will result in a horizontal line. Particles with a shape other than a sphere will show anisotropy and thus an angular dependence when plotting of  vs. .[3] The intercept will be in any case the . The  can be calculated into the hydrodynamic radius for a sphere through the Stokes-Einstein equation.

One must note that the cumulant method is valid for small  and sufficiently narrow .[4] One should seldom use parameters beyond , because overfitting data with many parameters in a power-series expansion will render all the parameters including  and , less precise [5].

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.[6].[7] The CONTIN analysis is ideal for heterodisperse, polydisperse and multimodal systems which 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 : 1E-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 a number of iterative steps to minimize the deviation of the fitted data from the experimental data and subsequently reducing the  of the fitted data.


For more information for Dynamic light scattering see our Product pages for Nanosight

Ankersmid M&C bvba          info@ankersmid-lab.be,  (32) - 3/820 80 01