Dark Field DDM characterization of roto-translational diffusion

Dark Field Differential Dynamic Microscopy enables the accurate characterization of the roto-translational dynamics of bacteria and colloidal clusters

Roberto Cerbino, Davide Piotti, Marco Buscaglia, and Fabio Giavazzi 1 Dipartimento di Biotecnologie Mediche e Medicina Traslazionale, Universit� degli Studi di Milano, via F.lli Cervi 93, 20090 Segrate, Italy roberto.cerbino@unimi.it, fabio.giavazzi@unimi.it


Micro- and nanoscale objects with anisotropic shape are key components of a variety of biological systems and inert complex materials, and represent fundamental building blocks of novel self-assembly strategies. The time scale of their thermal motion is set by their translational and rotational diffusion coefficients, whose measurement may become difficult for relatively large particles with small optical contrast. Here we show that Dark Field Differential Dynamic Microscopy is the ideal tool for probing the roto-translational Brownian motion of shape anisotropic particles. We demonstrate our approach by successful application to aqueous dispersions of non-motile bacteria and of colloidal aggregates of spherical particles.

1 Introduction

Understanding and quantifying Brownian processes is relevant for soft condensed matter scientists as well as for a wider audience that ranges from biologists to economists [[1], [2]]. As far as colloidal particles are of interest, the erratic nature of their Brownian motion is reflected in the well-known fractal appearance of their trajectories as well as in the irregular change of their orientation in time [[3], [4]]. In the past, rotational Brownian motion has received considerably less attention than its translational counterpart, in part because characterizing the rotational Brownian motion is more challenging. Most of the characterization makes use of optical methods such as Video Particle Tracking (VPT) , Dynamic Light Scattering (DLS), and Fluorescence Correlation Spectroscopy (FCS) [[5], [6], [7], [8]]. In particular, Depolarized Dynamic Light Scattering (D-DLS) has been shown to be a powerful tool to assess the roto-translational dynamics of an ensemble of anisotropic (by shape and/or optically) particles by analyzing the fluctuations in the depolarized scattered light intensity [[6], [9], [10], [5], [11]]. Very recently, it was shown that D-DLS experiments, usually requiring a custom laser-based optical setup, can be performed successfully with an optical microscope [[12]]. This approach, termed polarized-Differential Dynamic Microscopy (p-DDM), builds on Differential Dynamic Microscopy (DDM) that extracts scattering information from the quantitative analysis of time-lapse microscope movies [[13]]. While DDM in its original implementation probes the translational diffusion coefficient of a colloidal suspension observed in bright field microscopy, p-DDM gives also access to the rotational diffusion of anisotropic colloidal particles, provided that the sample is observed between suitably oriented polarizers. Compared to DLS and D-DLS, DDM and p-DDM offer a simpler implementation, more robust performances, higher flexibility, and a better rejection of both stray and multiply-scattered light [[14], [15], [16]].

Both D-DLS and p-DDM work well with spherical particles made of a markedly birefringent material and with non-spherical (optically homogeneous) particles provided that both the shape anisotropy and the optical contrast with the dispersion medium are large enough. These requirements are not always met and, as a result, the rotational dynamics of important classes of particles cannot be probed with D-DLS and p-DDM. In particular, particles that are quasi-index matched with the dispersion liquid or sparse aggregates exhibit a depolarized component of the scattered light that is insufficiently large for computing reliably the temporal autocorrelation function of the scattered field (p-DDM) or intensity (D-DLS). A quick and reliable tool for the quantitative characterization of the roto-translational dynamics of such samples is not yet available but it would be useful not only for the great potential of anisotropic particles in the self assembly arena [[17]] but also because many systems of biological interest such as bacteria and prokaryotic cells cannot be easily studied with D-DLS [[18]].

An alternative to D-DLS for the characterization of the rotational dynamics of shape anisotropic particles is available, at least for particles that are not too small. In fact, for an anisotropic particle whose longest dimension is comparable or larger than the inverse scattering wave-vector ( , typically), the intensity of the light scattered at a large angle depends upon the orientation of the particle itself and from the study of the fluctuations in the scattered light intensity one can measure the so called dynamic form factor [[19]], which encodes information on the rotational dynamics of the particles, as well as on other internal degrees of freedom, if present [[6]].

Here we show that Dark Field Differential Dynamic Microscopy (d-DDM), a recently proposed method for the characterization of the translational diffusion of colloidal particles [[20]], enables also the effective quantification of the intensity fluctuation associated with the rotational dynamics. Our key observation is that, in a conventional dark field microscope, the image of a small particle is mainly formed by light scattered at relatively large angles, while forward scattered light is not collected. This implies that the intensity associated with the dark field image of an anisotropic particle depends on its orientation. Accordingly, if the particle rotates, for example because of thermal motion, the intensity of the collected light fluctuates, leading to an image that “blinks” over time. This makes d-DDM a valuable tool also for the quantitative characterization of the roto-translational dynamics of colloidal particles, especially in regimes where p-DDM does not work, namely the case of micron-sized particles closely index-matched with the solvent. To demonstrate this approach we first present measurements on a dispersion of non-motile rod-shaped bacteria. The determination of the rotational and translational diffusion coefficients allows to determine the relevant dimensions of the bacterial particles. In a second set of experiments, d-DDM is applied to a suspension of spherical colloidal particles. In this case d-DDM is demonstrated to be extremely sensitive to the presence of even a small fraction of aggregates of particles, whose anisotropic nature leads to a significant rotational signal.

2 Materials and methods

2.1 Sample preparation

In order to assess the potential of d-DDM in characterizing the roto-translational dynamics of anisotropic colloidal particles, we have employed two different kind of samples that are difficult to study with other tools such as p-DDM and/or D-DLS. In particular, we studied aqueous dispersions of non-motile bacteria and of aggregates of spherical nanoparticles, both presenting shape anisotropy and a moderate refractive index difference with the dispersion medium (water).

Non-motile bacteria

One of the most frequently used E. coli strain for routine biological cloning applications is the DH5 strain, which are non flagellated bacteria. They are almost non-motile, less fragile to handle and easy to grow, which makes them the perfect anisotropic particles for our purposes. We grew single colonies from frozen stocks on Luria broth (LB) agar plates at °C overnight. A single colony was transferred from a plate to ml of liquid LB and incubated overnight ( h) at while shaken (for aeration) at 200 rpm. The next step was to transfer cells from LB to a minimal medium with no exogenous nutrients in order not to have bacterial reproduction and to minimize the growth rate. To this aim, the medium underwent a washing process consisting of a 2 min centrifugation (6000 rpm and 2000 g), the expulsion of supernatant and the resuspension of bacterial pellet in 1 ml of PBS (phosphate-buffered saline, a water-based salt solution containing sodium hydrogen phosphate and sodium chloride). This washing procedure was repeated 3 times (2 centrifugations). After this, the system was sucked into a mm capillary and the capillary was fixed to a glass slide with vaseline petroleum jelly, a non-toxic glue. This operation was carefully performed in order not to let air get inside the capillary that could lead to unwanted sample drifts. It was also important to carefully avoid liquid residuals between the capillary and the glass slide, which would have caused a dynamic drying front expanding during the microscope acquisition and affecting the DDM experiments.

To avoid bacteria sedimentation, which impact on the dynamics, the density mismatch between bacteria and the PBS physiological medium was compensated by adding Percoll to the dispersion. Percoll is routinely used for density gradient centrifugation of cells, viruses, and sub-cellular particles. It consists of colloidal silica particles of nm diameter prepared at % w/w in water. The silica particles are coated with polyvinylpyrrolidone (PVP), which makes them completely non-toxic and ideal for use with biological materials. PVP is randomly bound to the silica particles as a monomolecular layer. The size of these particles is so small that the intensity of the scattered light is negligible when compared to the intensity of the light scattered by bacteria, as we checked experimentally in preliminary DDM experiments. We found that a sample made of % of PBS medium with bacteria and 80% of Percoll was stable upon centrifugation for 2 hours at 2000 g, which is a time longer than the 30-40 minutes needed to acquire microscope videos in bright field and dark field. We note that for large concentrations of bacteria, bacterial aggregation was found to occur in the presence of Percoll , an effect presumably attributable to depletion interactions. For this reason we worked at a low bacterial concentration (about bacteria/ml, corresponding to a volume fraction of about ), below the threshold needed to trigger appreciable aggregation during our experiments. The viscosity of the PBS-Percoll solution has been measured by a capillary viscosimeter for different temperatures in the range °C. For °C, the temperature at which the experiments with bacteria have been performed, we found .

Polystyrene particles

We used spherical polystyrene particles (Spherotech (TM) SPHERO (TM) Biotin Polystyrene Particles), with a certified mean diameter equal to (intensity-weighted Nicomp distribution rescaled to number density). The samples were prepared by gently vortexing the bottle in order to resuspend the colloidal particles. Serial dilutions in PBS buffer lead to a final concentration of about particles/ml, corresponding to a volume fraction of about . The solution was then sonicated for min and confined into a rectangular glass capillary (VITROCOM, internal size: mm) for the microscopy observations. The capillaries were sealed on both sides with UV glue (UV30-20, Loxeal s.r.l., Cesano Maderno, Italy), cured for min under UV light (VL-6.M, Vilbert Lourmat, Marne la Vall�e, France) and successively loaded on the microscope. The particles concentration was low enough to ensure that single particles could be identified and tracked when observed under the microscope, which allowed a direct-space based characterization of the particles trajectories, in addition to DDM. The viscosity of the PBS solution at the temperature °C at which the experiments have been performed has been estimated by using literature value for pure water at the same temperature ( ).

2.2 Differential Dynamic Microscopy

Microscopy measurements were performed with a Nikon Eclipse Ti-E commercial microscope equipped with a Hamamatsu Orca Flash 4.0 v2 camera (pixel size m). Dark-field images are collected with a standard microscopy objective ), while the sample is illuminated with a condenser stage () coupled with a PH3 phase-contrast ring mask. Bright field images are also collected with the same objective lens and bright field illumination, by using the same condenser stage with a standard diaphragm. For bacteria we also used a phase-contrast objective (), the sample being illuminated through the same condenser stage used in the previous cases, using a proper ring mask in order to achieve the phase-contrast condition. Each acquisition typically corresponds to a sequence of images acquired with a frame rate equal to fps in the case of the bacterial suspension and to fps for the latex particles. Movies acquired in dark field microscopy exhibit a characteristic blinking due to the rotation of anisotropic particles, which is not present in bright field or phase contrast movies (see Supplementary Movies SM00 and SM01, respectively). To extract quantitative information from these movies we analyzed them by using the standard DDM algorithm [[21], [13], [15]], which is based on calculating the difference between two images acquired at times and . Once this quantity is obtained, its spatial Fourier power spectrum is computed by using a Fast Fourier Transform (FFT) routine and, in the presence of stationary or quasi-stationary statistical processes, an average over power spectra with the same but different reference time is obtained, which increases the statistical accuracy of the data. This leads to the so called image structure function (ISF)


that captures the dynamics of the sample as a function of the two-dimensional scattering wave-vector and of the delay time . The ISF is connected to the (normalized) intermediate scattering function [[6]] by the relation


where is a term that accounts for the camera noise, is an amplitude term that contains information about the static scattering from the sample and details about the imaging system [[15]]. The form of the intermediate scattering function for a sample observed under dark-field imaging, undergoing both translational and rotational Brownian motion will be discussed below in Section 3.1.

In general, the two-dimensional nature of the ISF provides a powerful means to probe the sample dynamics along different directions in the plane that may be of particular interest for the problem under study [[22]]. However, whenever the ISF bears a circular symmetry, as for all the experiment described here, azimuthal averaging of is often used to obtain the one-dimensional function , of the radial wave-vector [[13], [23], [24], [25], [26], [27], [28]].

In order to obtain a reliable determination of the sample dynamics, which is encoded in the intermediate scattering function, a robust procedure for estimating the amplitude and the noise contribution in Eq. 2 is required. Dark-field microscopy experiments are very sensitive to non-idealities such as like dust particles and scratches on the optical elements and/or on the sample cell, which can give a significant background signal on top of which the signal from the particle is superimposed. This background signal appears as an additive, positive term to the image intensity distribution. The intensity distribution associated with each image can be thus written as the sum of three independent terms:


where is a background image (i.e. the intensity distribution that would be observed in absence of the sample), is the contribution to the image from the particles and is the camera noise. In our case, the main contribution to the camera noise is from the shot noise and we can safely assume that it has zero average: and that it is delta-correlated in both space and time: . If the amplitude of the noise is small compared with the background intensity and if the density of the scatterers is low enough to ensure that only a finite fraction of the image is covered by the particles, an accurate reconstruction of the background image can be obtained as . In fact, by picking up for each pixel the lowest intensity value registered during a suitably large time window - larger than the diffusion time of a single particle over its image - allows to minimize the additive contribution from the particles themselves. Once an estimate for has been obtained, an estimate of the total amplitude of the fluctuating parts can be extracted from , where we have made use of Eq. 3. Eq. 2 can be thus rewritten in terms of as


can be estimated with high accuracy as the intercept for of . This could be obtained in practice by fitting over a small interval to a polynomial function and by taking the -th order coefficient. This procedure sets the value of the amplitude of the first term, allowing thus a reliable estimate of the relaxation times in even if they exceed the width of the acquisition window.

3 Theory

In this Section we will first provide a brief summary of the scattering theory from optically anisotropic particles [[6]]. In the second part we will describe the features of a dark-field imaging system and the imaging process of anisotropic particles.

3.1 Scattering by anisotropic particles

Small anisotropic particles: Rayleigh scattering

The description of the scattering of light by a particle much smaller than the wavelength is usually based on the so-called Rayleigh approximation [[29]]. For simplicity, we will restrict our discussion to the case of uniaxial particles, whose polarizability tensor admits a diagonal form with diagonal elements , , , where . For this kind of particles, the anisotropy parameter is defined as and the average (excess) polarizability as . Here is the particle volume and is the refractive index of the dispersion medium. If a plane wave electric field of wave-number impinges on such particle, the latter emits a scattered field that, in addition to the component that is parallel to , also bears a perpendicular component . One has [[29]]




where is the imaginary unit, is the distance from the particle, and and are dimensionless amplitudes that depend on the scattering angle measured with respect to the direction of the incident radiation, and on the orientation of the particle. For small scattering angles, the scattering amplitudes and are given by the following expression [[6]]


Here is the spherical harmonic function of order , .

Given the shape and the refractive index distribution within a particle, the calculation of its polarizability tensor is not in general a trivial task. Even in the Rayleigh regime, closed-form expressions can be obtained only in presence of particularly simple geometries. In order to better elucidate how the scattering properties of a particle depend on its shape and refractive index, it can be worth considering a specific model system for which a simple analytical description is available, namely a homogeneous ellipsoid of revolution (spheroid) of semi-axes . In this case, the diagonal elements of the (excess) polarizability tensor are given by , where is the ratio between the refractive index of the particle and the refractive index of the solvent and , , with . is a parameter ranging from 0 (flat elliptical disk) to (infinitely long ellipsoid). The value corresponds to the isotropic case. We can also define an anisotropy factor in such a way that corresponds to a sphere, while one has for an infinitely long ellipsoid. If we also define , we obtain the following simple expression for the amplitude depolarization ratio [[6]]


which in fact is the ratio between the amplitudes of orientation-dependent component and the orientation-independent component of the scattered field in Eqs. 6, 8. It is clear from Eq. 9 that the contrast in the fluctuation of the depolarized scattering component depends both on the anisotropy of the particle and on the refractive index mismatch with respect to the solvent. Even a strongly anisotropic particle () cannot produce a significant depolarized signal if it is quasi index-matched with the solvent (). We note that, as far as the Rayleigh approximation is satisfied, the result expressed by Eq. 9 holds independently of the particle size.

Weakly scattering anisotropic particles of larger size: Rayleigh-Gans-Debye description

The case of a homogeneous, quasi-index-matched particle can be adequately described within the Rayleigh-Gans-Debye (RGD) approximation if the overall phase delay associated with the the particle is small: [[29]]. In this approximation, the amplitude of the wave scattered by the particle can be calculated as the sum of independent contributions from each portion of particle. In the RGD approximation the amplitude of the depolarized scattering is negligible and the scattering amplitude entails only a polarized component




is the so-called form factor amplitude and is the transferred momentum, i.e. the difference between the scattering wave-vector and the wave-vector of the incident light. For a particle of a given shape, has in general an implicit dependence on the particle orientation. Within the RGD approximation, the scattering in the forward direction () does not depend on the details of the particle shape or orientation, as . A non-trivial dependence of the scattered amplitude on the particle orientation can be observed only for and if the size of the particle is not too small compared with the wavelength of light [[6]].

For a spheroidal particle, such as the one considered in the previous paragraph, the form factor amplitude can be calculated explicitly:


where . Here is the unit vector oriented along the particle axis and the angle between and the transferred momentum . In contrast with the forward depolarized scattering in the Rayleigh regime, here the relative amplitude of the fluctuation in the scattered intensity due to a rotation of the particle does not depend at all on the optical contrast, while is strongly influenced by the overall particle size and by the collection angle, which in turn determines the transferred momentum .

This property, namely that fact that fluctuation in the intensity of the light scattered at a large angle by an ensemble of anisotropic colloidal particles reflects also its rotational dynamics, has been exploited in DLS for accessing its roto-translational dynamics, as a complementary approach with respect to low angle depolarized scattering measurements [[19], [30], [31]]. Eq. 12 takes a particularly simple form if the particle is not too large compared to the inverse transferred momentum. In fact, for a relatively small uniaxial anisotropic particle (for which , where is its longer dimension), the expression in Eq. ,11, can be expanded in , leading to:


or, to the same order in , to


where and are the diagonal elements of the particle tensor of inertia [[6]], and is the angle between the axis of the particle and the transferred momentum . The simple harmonic dependence of the scattered intensity on described by Eq. 14 will be used in the following paragraphs to link the correlations of the scattered intensity to the statistical properties of the rotational motion of the particle.

3.2 Dark-field microscopy

Dark field imaging

In the context of optical microscopy, the term dark field (DF) indicates a family of microscopy configurations characterized by the fact that the transmitted illumination beam is not collected by the imaging optics and thus only the light scattered from the sample contributes to the image. In practice, this can be obtained in a number of different ways, for example by using dedicated illumination stages.

Figure 1: a) Schematic representation of an incoherent dark field microscope (IDFM) with Koehler illumination. b) Elementary scattering process contributing to the image formation in an IDFM. denotes the incident wave-vector, is the scattering wave-vector, is the transferred momentum, and the wave-vector oriented along the optical axis. The transferred momentum can be also written as , where and . c) Schematic representation of a coherent dark field microscope (CDFM) in a custom bench-top homodyne near-field scattering set-up (see text for additional details). d) Scattering process in CDFM. The incident wave-vector is parallel to optical axis. The transferred momentum is given by , where is the scattering wave vector.

A common implementation, which is compatible with most commercial inverted microscopes is the one reported in Fig. 1 a), where a schematic representation of a microscope with Koehler illumination is shown. In this configuration, a circular aperture or radius carved into an opaque mask is placed in the back focal plane of the condenser lens, of focal length . In this condition the object plane is illuminated only by rays forming an angle , where . The dark field condition is achieved by using an objective lens with numerical aperture , so that the transmitted beam is not collected by the objective. In practice, this is commonly achieved by coupling a high numerical aperture phase contrast ring with a low power objective. In this configuration, the illumination beam can be thought of as the incoherent superposition of many plane waves propagating at angles with respect to optical axis i.e. the sample is illuminated by a collection of uncorrelated coherent patches of transverse size . Since , the size of these coherent patches is smaller than the transverse size of the objective PSF. For this reason, the imaging process is to all effects incoherent, as interference effects between different scattering centers are negligible. For this reason, we refer to this configuration as Incoherent Dark-Field Microscopy (IDFM).

A different implementation of a dark field imaging system is the one shown in 1 c), which is typically adopted in homodyne imaging and/or near-field scattering set-ups [[32], [14]]. In this case the sample is illuminated by a highly coherent beam and the transmitted light is blocked by a small opaque patch placed in the back focal plane of the objective lens. We call this configuration Coherent Dark Field Microscopy (CDFM), since the sample is ideally illuminated by a single plane wave propagating along the optical axis and the light scattered by different points of the sample always bears a well defined relative phase. The superposition of the scattering patterns from different points of the sample takes thus place on a coherent basis, since interference effects cannot be neglected and the imaging process is linear in the complex amplitude rather than in the intensity, as it is for IDFM [[33]].

According to the nomenclature introduced in [[23]], IDFM is thus a linear imaging system, whereas CDFM is not, in that the image intensity recorded in presence of two particles is not given by just the sum of the intensities associated with the two particles imaged separately. For this reason, IDFM is to be preferred for Digital Fourier Microscopy experiments, where a linear space-invariant imaging system turns out to beneficial [[15]]. We note that, as pointed out in Ref. [[20]], dark field imaging can in some cases lead to a non-homogeneous illumination pattern, which violates the space-invariant assumption. In our optical set-up this effect was negligible, in that the illumination was found to be uniform across the entire field of view.

The incoherent dark field microscope as a fixed, non-zero angle scattering set-up

In a dark field microscope the main contribution to the particle image intensity comes from light scattered at a non-zero angle, roughly corresponding to the illumination numerical aperture. As a consequence, if the particle is anisotropic, its reorientation produces a modulation in the image intensity. This property opens to the possibility of exploiting IDFM to study the rotational dynamics of anisotropic particles. The main elements of typical dark field microscopy setup are shown in the simplified scheme in Fig. 1. The elementary scattering process contributing to the image formation is shown in Fig. 1 b), where denotes the incident wave-vector, is the scattering wave-vector and is the transferred momentum. By introducing as the wave-vector oriented along the optical axis, having the same amplitude of and we can express the transferred wave-vector as , where and . In the image plane, the total intensity associated with the particle’s image can be calculated as the (incoherent) sum of all such processes


where the function , representative of the angular distribution of the illumination beam, weights the contribution of each incoming plane wave, the particle form factor accounts for the scattering properties of the sample, and the incoherent transfer function quantifies the collection efficiency of the objective lens.

Our choice of a thin phase contrast ring of numerical aperture makes the evaluation of the integral on the right hand side of Eq. 15 particularly simple, as , where . In the general, is an azimuthally symmetric function of width centered around . If we take the (rather crude) approximation: , corresponding to considering only the light scattered along the optical axis, we obtain the following simple result, where the total intensity of the particle’s image is written as an azimuthal average of the particle form factor performed for : , where and are the projections of the transferred momentum in the direction perpendicular and parallel to the optical axis, respectively . In IDFM, the main contribution to the particle image intensity comes from light scattered at a non-zero angle, roughly corresponding to the illumination numerical aperture. If we assume for the form factor the expression given in Eq. 14, can be easily integrated leading to:


or, equivalently


where is the spherical harmonic of order 2,0, is the angle between the axis of the particle and the optical axis, is a constant amplitude and is a constant. Eq. 17 shows that the intensity signal due to an anisotropic particle depends on the particle orientation with respect to the optical axis, which explains why a particle that undergoes a rotational Brownian motion appears as a randomly blinking object. We note that considering a finite value for the objective numerical aperture does not affect this result. In fact, in that case is still given by the integral of over an annular region of radius but with a finite width , which does not affect the angular dependence of in Eq. 17 but only the prefactor . We note that Eq. 17 holds under the same hypotheses under which Eq. 14 is valid, namely that the , where is longer dimension of the particle. For very large particles a more complex expression is expected, involving higher order spherical harmonics.

Dark field DDM probes the roto-translational dynamics of anisotropic particles

The fact that the image intensity corresponding to one particle depends on the particle orientation can be in principle used to assess the rotational dynamics by studying the characteristic time of blinking in movies obtained by IDFM (see Supplementary Movie SM01 and SM02). For the case of interest in this work, i.e. when the particle undergoes a rotational Brownian motion, its rotational diffusion coefficient can be obtained by calculating the intensity temporal auto-correlation function


whose characteristic time would give immediate access to the rotational diffusion coefficient [[6]]. However, measuring is usually quite difficult: for instance, particles can disappear from the image when they exit the focal region or two particles might superimpose along the optical axis and give rise to spurious effects. These difficulties can be bypassed by working in the wave-vector space, as done in Dark Field Differential Dynamic Microscopy (d-DDM)[[20]]. The image intensity distribution of a collection of identical particles observed in IDFM is given by


where is a function that describes the distribution of the individual scatterers in real space. If is the spatial Fourier transform of , an alternative quantification of the particle dynamics is provided by the intermediate scattering function


where is the translational diffusion coefficient of the particles. Eq. 20 can be rewritten by introducing the normalized intermediate scattering function [[6]]


where the contribution of the two relaxation modes with decay rates and can be appreciated. This form for the intermediate scattering function is pretty common, as is found also in D-DLS and p-DDM experiments [[5], [12]].

4 Results and discussion

In this Section, we demonstrate d-DDM as an effective powerful tool for the simultaneous determination of the translational and rotational dynamics by presenting results obtained with two representative samples: a suspension of quasi-index matched anisotropic colloidal particles (non-motile rod-shaped bacteria) and a suspension of quasi-monodisperse spherical colloids both in a non-aggregated state and forming small clusters.

4.1 Bacteria

A first set of measurements was performed on a suspension of density-matched non-motile coliform bacteria, prepared as described in the Material and Methods Section. Movies of the sample are acquired with dark field, bright field and phase contrast microscopy. A representative dark field image of the sample is shown in Fig. 2 (see also Supplementary Movie SM01): the bacterial particles appear as bright spots diffusing on a dark background and whose intensity fluctuates with a characteristic time of about . We will show that these two processes, namely, the concentration fluctuation due to the translational motion of the particles and the intensity fluctuations caused by the rotational dynamics of single particles, are well captured by the d-DDM analysis.

Figure 2: Main panel: representative dark field image of the bacterial suspension, where the optical background has been subtracted as described in the Material and Methods Section. Inset: representative phase-contrast image of the same sample.

In Fig. 3 we show for some representative values of the wave-vector , logarithmically spaced in the in range , the intermediate scattering functions obtained with d-DDM analysis. As expected from Eq. 21, the observed relaxation exhibits two distinct decays. This is particularly evident in the low- regime, where the time-scale separation between the two decays is more pronounced (see Fig. 3 b). The intermediate scattering functions are well described by a sum of two simple-exponential decays (continuous lines in Fig. 3). Fitting the obtained curves to Eq. 21 provides thus an estimate for the two -dependent relaxation rates and . The so-obtained and are shown in Fig. 4 together with the best fitting curves: , and . From this last fit we obtain an estimate for the translational and the rotational diffusion coefficient of the bacteria: and .

Figure 3: Symbols: intermediate scattering function obtained for different wave vectors , ,, , , , , , in d-DDM experiments. Continuous lines: best fitting curves with the model function given in Eq. 21.

The simultaneous measurement of both and allows estimating the size of the bacterial particles. We consider two simple models, for which analytical expressions for the diffusion coefficients are available, namely a spheroid of semiaxis [[34]] and a cylinder of length and radius [[35], [36]]. Explicit expression for and in the two cases are reported in Appendix A and can be inverted numerically to determine the best estimates for the geometrical parameters in the two cases: , for the spheroid and , for the cylinder. We note that both models provide meaningful results, that are in good relative agreement and are fully compatible with both literature values [[37], [38]] and direct, high magnification microscopy observations (see Fig. 2).

As a consistency check, we also analyzed bright field and phase-contrast movies of the same sample. In both cases, the obtained intermediate scattering functions show a single decay that is very well fitted by a simple exponential function (data not shown). The corresponding -dependent relaxation rates are reported in Fig. 4 and exhibit a very clean scaling and the estimated translational diffusion coefficient is in excellent agreement with the results of d-DD. Of note, the wave-vector range over which a reliable estimate of the dynamics can be obtained is not the same for all the imaging conditions. In particular, we observe that, although dark field and bright field measurements are performed with the same objective, frame rate and duration of the acquisition, dark field DDM is much more effective in probing the dynamics in the the low- regime. This can explained by inspecting the static scattering amplitudes obtained from the DDM analysis in the two datasets show in the nset of Fig. 4, which outline the effect of the different transfer function of the two methods. The bright field amplitude shows the characteristic depression at low , reflecting the fact that bacteria, being quasi-index matched with the solvent, behave as phase objects [[23]]. On the contrary, the dark field amplitude is a monotonically decreasing function of . The phase contrast data, being obtained with a different microscope objective, with larger magnification and numerical aperture (40 vs 10 and 0.6 vs 0.15), cover a range shifted by approximately half a decade to larger wave-vectors.

Figure 4: Decorrelation rates as a function of the wave-vector obtained from the fit of the intermediate scattering functions obtained in dark field (open blue circles and squares), in bright field (orange crosses) and in phase-contrast (purple circles) microscopy experiments on a dilute bacterial dispersion. The dash-dotted and the continuous lines are the best fits to the curves and , respectively, which lead to the estimates and . Bright field and phase-contrast data provided the value , obtained by fitting the experimental data to the curve (fitting curve not shown). Inset: static amplitude obtained with the same microscope objective from DDM analysis of bright field (red squares) and dark field (blue circles) movies of the same bacterial dispersion. While the dark field amplitude is a decreasing function , the optical transfer function for the bright field experiments is characterized by a visible depression at small scattering wave-vectors.

4.2 Spherical colloids

A second set of d-DDM measurements was performed on a suspension of monodisperse spherical latex particles. A first sample (Sample 1) was prepared with a 5 minutes sonication stage before the measurement, as described in the Materials and Methods Section. Surprisingly, we found that d-DDM analysis showed two distinct decays and for the ISF (5), as previously found for the bacteria. This finding contrasted our expectation to observe a single relaxation mode due to translational diffusion of the particles. In order to better understand the reason of such unexpected behavior we thus performed a bright field DDM experiment on the same sample: as expected for a reasonably monodisperse sample, the experimental ISFs were well fitted to a single exponential decay and the so-obtained relaxation rate was well fitted to the function (Fig. 5). The estimate µm/s for the translational diffusion coefficient of the particles is fully compatible with the value found in dark field experiments, which gave µm/s and for the translational and rotational diffusion coefficients, respectively. However, none of them was found to be compatible with the value , certified by the producer for the hydrodynamic radius of the particles. Both values and , obtained from and , respectively, via the Stokes-Einstein relation are about larger than the nominal value. We interpreted all these results as consequences of the presence of anisotropic particles in suspension, originated from aggregation of the spherical particles. Careful inspection of dark field movies (Supplementary Movie SM02) further supported this hypothesis in that it pointed to the presence of a small number of blinking particles.

To obtain final confirmation, we performed measurements on a second sample (Sample 2) that was carefully prepared from the same batch of particles by using a longer sonication stage ( minutes instead of ). With Sample 2, the results of DDM analysis of bright field and dark field movies (see Supplementary Movie SM03) confirm the absence of aggregates. In both cases, the ISFs exhibit a single exponential decay, with a relaxation rate displaying a clean quadratic scaling with , from which we obtain the estimate µm/s for the translational diffusion coefficient (Fig. 5). The extracted hydrodynamic radius is now fully compatible with the nominal size of the particle, as certified by the producer. We could therefore safely conclude that Sample 1 contained aggregated particles and that d-DDM provides a very powerful means to spot the presence of small aggregates in colloidal dispersions of spherical particles. In particular, when sizing unknown samples, it would be recommended to complement bright field DDM experiments with dark field ones, to check whether a rotational diffusion decay mode due to aggregates is present or not.

We then turned to assessing how far can d-DDM be brought in obtaining information on the aggregates. To this aim we evaluated the values expected for the translational and rotational diffusion coefficients of small clusters of particles (dimers, trimers, tetramers, etc.) [[39]]. In particular, the value µm/s, corresponding to the hydrodynamic radius certified by the producer, gives the following values for small clusters: µm/s (dimers), µm/s (trimers), and µm/s (tetramers). We note that these values are very close to each other, which explains why translational diffusion does not discriminate very effectively between these contributions. This is confirmed by our d-DDM results, as we found that the translational diffusion coefficient extracted with d-DDM exhibited an intermediate value between the expectation for monomers and for the small clusters. By contrast, the rotational diffusion coefficient is not sensitive to the presence of monomers. Consistently, the value obtained with d-DDM ( ) was found to be within the range of the rotational diffusion coefficient expected for small clusters: (dimers), (trimers), and (tetramers).

Figure 5: Decorrelation rates obtained from the fit of the intermediate scattering functions obtained in dark field (open blue circles and squares) and in bright field (orange dots) microscopy experiments on a partially aggregated suspension of spherical colloidal particles (Sample 1). The dash-dotted and the dashed lines are the best fits to the curves and , respectively, which lead to the estimates and . Bright field data provided the value , obtained by fitting the experimental data to the curve (fitting curve not shown). Upward purple triangles correspond to the single decorrelation rate measured in a dark field microscopy experiment on a non-aggregated suspension of the same particles (Sample 2), The continous line is the best fit to the curve , leading to the estimate .

In order to check the consistency of this picture, we performed a detailed Video Particle Tracking (VPT) analysis on the same movies, which is very delicate and time consuming but provides a more detailed information about the sub-populations of which the sample is composed. Particles trajectories are obtained with the Particle Tracker Plugin, included in the Mosaic Suite for ImageJ. With a custom software written in Matlab, we extracted from each trajectory the particles mean squared displacement and, by fitting the resulting curve as a function of the time delay we obtain and estimate for the translational diffusion coefficient of each single particle. By rejecting trajectories shorter than time steps, we obtained for about particles for each sample. In Fig. 6, we report the histograms representing the distributions of the values of obtained for each of the two samples. For Sample 2, the distribution of is fairly symmetric and it is well described by a Gaussian function with mean value and standard deviation . This result is compatible with a moderately dispersed distribution peaked around a value that is in excellent agreement with estimated µm/s obtained for the translational diffusion coefficient from the d-DDM analysis. For Sample 1, the same SPT-based analysis provides a completely different result. The distribution of is broader and, beside a peak centered about a value compatible with , a secondary peak for is also clearly visible. As shown in Fig. 6, the range covered by this secondary peak is compatible with the translational diffusion coefficients expected for small clusters of monomers, assuming as the diffusion coefficient of a monomer. These results strongly corroborate our interpretation and confirm that d-DDM can be used to measure with high sensitivity the translational dynamics of spherical colloids and to spot aggregation of spherical particles in a simple and effective way.

Figure 6: Histograms of the distribution of translational diffusion coefficient for Sample 1 (large blue bars) and for Sample 2 (thin orange bars) as determined from VPT. The continous vertical line correspond to , obtained by fitting the histogram obtained from Sample 2 to a Gaussian function (continous curve). The histogram for Sample 1 is well described as the sum of two Gaussian functions (heavy dashed-dotted curve) centered in and , respectively. This second value falls within the range of the translational diffusion coefficients expected for small clusters of monomers, assuming as the diffusion coefficient of a monomer. The vertical dashed-dotted, dashed and dotted lines correspond to the expected values for dimers, trimers and tetramers, respectively.

5 Conclusions

In this work, we have shown that the recently introduced Dark Field Differential Dynamic Microscopy [[20]] is a simple and powerful tool for the simultaneous determination of the roto-translational dynamics of anisotropic microparticles in suspension. Our experiments with bacterial suspensions showed that d-DDM is somehow complementary to the recently proposed polarized-DDM [[12]], in that the latter may fail with particles characterized by moderate shape anisotropy and good refractive index-matching with the dispersion medium.

We also showed that the peculiar nature of the dark field signal associated with the rotational dynamics, makes d-DDM very effective in spotting and quantifying the presence of anisotropic aggregates of isotropic particles, even though whenever the particles are large enough, particle tracking may provide a better tool for a detailed analysis of aggregated samples in which clusters are present. On the other hand, dDDM analysis is statistically more robust and does not require the intervention of an experienced user for the fine tuning of the parameters involved in the image processing procedure.

These promising results open to the possibility of studying with d-DDM more complex and challenging systems. For example, it is know that, for motile bacteria, the rotational dynamics is very different from the purely Brownian one [[40]] and that a non-trivial interplay exists between rotational and transational degrees of freedom [[41]]. d-DDM could allow the simple high-throughput, characterization of this complex dynamics, possibly combination with other quantitative microscopy methods, like for example standard phase contrast DDM [[42]] or the so-called dark field flicker microscopy [[43]], that has been used for monitoring the rapid periodic fluctuation associated with the beating of flagella.

Another field of potential application is the optical characterization of the mechanical properties of soft materials, the realm of microrheology [[44], [45]]. Passive microrheology, in particular, exploits the thermally excited positional fluctuations of immersed tracer particles to probe the viscoelastic moduli of the hosting fluid [[46]]. Very recently, DDM has been demonstrated to be a reliable route to microrheology, enabling the accurate, tracking-free determination of the mean squared displacement of probe particles in a variety of imaging conditions [[47], [48]]. In view of the results presented in this work, we expect that d-DDM, in combination with calibrated anisotropic tracers, could provide the ideal ground to extend these ideas also to the rotational degree of freedom, enabling the simultaneous execution of translational and rotational microrheology experiments [[49], [50]]. While this combination could appear redundant in the case of a perfectly homogeneous fluid - where the two approaches are expected to give equivalent results - it could quite valuable in the presence of sources of non-ideality in the system (e.g. inhomogeneity of the matrix or specific tracer-fluid interactions altering the boundary conditions at the surface of the particles). Since these effects, that can seriously compromise the reliability microrheology results, are expected to have a different impact on rotational and translational degrees of freedom, the execution of a combined experiment could allow to spot them effectively.

Moreover, since d-DDM is very sensitive in detecting the presence of anisotropic particles, it can be exploited as a quality control step during the preparation or the execution of experiments involving allegedly spherical colloidal particles or for real-time monitoring of aggregation processes and self-assembly. Another intersting application of the method could be in the on-line moniting of water quality [[51]]. In this case, the ability of d-DDM to spot the presence of particles in solution and, by providing an estimate of their dimensions and optical contrast, to discriminate between them, could make it an useful screening tool for the automatic identification of particularly important classes of contaminants, in primis bacteria.

We thank Silvia Biffi, Filippo Saglimbeni and Giovanni Tagliabue for help in the preparation of the samples and for early contributions to the project. We also thank Roberto di Leonardo and Claudio Maggi for insightful comments and discussions. RC, DP and FG acknowledge funding by the Italian Ministry of Education and Research, Futuro in Ricerca Project ANISOFT (RBFR125H0M) and by Fondazione CARIPLO-Regione Lombardia Project Light for Life (2016-0998). MB acknowledges funding from the European Union’s Seventh Framework Programme (FP7) for Research, Technological Development and Demonstration through the NAPES project (grant agreement no. 604241).

Appendix A Roto-translational diffusion of spheroids and cylinders 21

We report here analytic expression available in the literature for the rotational and translational diffusion coefficient for a spheroidal particle [[34]] and for a cylinder [[35], [36]]:


The rotational and translational diffusion coefficients are given by, respectively

Here , , and are the minor and the major semi-axis of the particle, respectively, is the Boltzmann constant, is the absolute temperature and is the solvent viscosity.


The rotational and translational diffusion coefficients are given by, respectively

where and are the translational diffusion coefficients along the direction perpendicular and parallel to the axis of the particle, respectively:

In the above expression , . and are, respectively, the diameter and the length of the cylinder. All other quantities are defined as in case of the spheroid.


Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description