A. Mikhailovskaya*a,
J. Crassousb,
A. Salonena and
D. Langevina
aLaboratoire de Physique des Solides, UMR 8502 Université Paris-Saclay, bât 510, 91405 Orsay cedex, France. E-mail: alesya.mikhailovskaia@u-psud.fr
bInstitut de Physique de Rennes, UMR 6251, Université de Rennes 1, Campus de Beaulieu Bâtiment 11A, 35042 Rennes Cedex, France
First published on 15th February 2016
Diffusing wave spectroscopy (DWS) was used to determine the size and volume fraction of nanoparticles (NP) within a foam taken as an example of a turbid media. The presence of two types of scatterers – dispersed NPs and liquid/gas interfaces – leads to two-fold dynamics in the system. Two characteristic decays of the temporal autocorrelation function are observed, their position and amplitude are dependent on particle concentration and foam age. Since only those NPs which are performing Brownian motion decorrelate the signal, one can follow the dynamics of particles' trapping into the foam matrix during the evolution of the system. This technique is a direct and noninvasive way to detect NPs in turbid media. The detection limit strongly depends on the NP nature and an algorithm for its estimation is provided.
An additional difficulty in the case of foam is that their study requires the application of non-invasive methods: treatments such as dilution, drying, freezing etc. result in sample damage. The non-invasive method of diffusing wave spectroscopy (DWS)3,4 has therefore be widely used for the characterization of foam evolution. With this method, the frequency spectrum of the transmitted or reflected light can be related to the multiple scattering events of photons. Information about structural rearrangements can be obtained.5,6 Instead of the light spectrum, and as in classical dynamic light scattering (DLS),7,8 the correlation function of the scattered light is generally measured, this function being the Fourier transform of the light power spectrum. Recently DWS was also applied to more complex dispersions containing colloidal particles dispersed in a foam,9,10 where the system evolution was determined both by Brownian motion of the freely diffusing colloids and by aging of the foam. Taking this foam mixed with nanoparticles as a model system, it was shown that the interpretation of the results can be considerably biased because of the traces of dispersed colloids.9 This can be particularly important in food science since DWS is a traditional tool for food microrheology11 and it is commonly used to follow the dynamics of gelling systems.12
The corresponding theoretical model was proposed in ref. 9, however the limits of application of the method remained open for discussion. We used a specially designed measuring cell which allowed to carry out the experiments for different thicknesses but in the same foam and at the same age of the sample. This method allows to prove the validity of the model and to determine the range of applicability of this DWS method. We studied a Gillette shaving foam with two different types of nanoparticles, silica and latex.
Another motivating question was the feasibility of the measurement of particle concentration in turbid media. This is not a trivial problem for the light scattering methods, although the use of DLS13 and diffused photon density wave spectroscopy (DPDWS)14 are reported in the literature. However DLS requires considerable dilution of the sample, which is complicated with delicate samples. The use of DPDWS is not very common although it has been used to measure micellar aggregation in a turbid solution. In the present work, we have studied in detail the issue of measuring the particle concentration within the evolving matrix.
The article is structured as follows: we start with an introduction to DWS, which finishes with a short summary. The experiments are then described, and we show how the concentration of nanoparticles can be measured in Gillette foam.
(1) |
A photon traveling through highly scattering media undergoes many scattering events and changes its initial direction of propagation. Once the photons have travelled over the distance l*, the transport mean free path, the direction is assumed to be totally randomized and light propagation can be modeled with a diffusion equation.3 Thus the photons arriving to the detector have taken random paths of various lengths. The paths are distinct from each other with a different number of scattering events and, consequently, they create different phase variations. We note Δφs(t, t + τ) the phase difference for a path of length s between times t and t + τ. Temporal fluctuations reflect normalized path distribution P(s)
The autocorrelation function can be expressed as
gE(τ) = ∫P(s)gE(s)(τ)ds | (2) |
If light is scattered by a dispersion of particles, their Brownian motion causes temporal fluctuations of intensity and the rate of signal decorrelation is related to particle size: the autocorrelation function decays faster in the case of smaller particles.3 The character of scattering in foams is different. Light is scattered by the liquid/gas interfaces and the decorrelation of the signal arises from the reorganization of the foam skeleton during the evolution of the system.16 The characteristic time between structural rearrangement events determines the position of the autocorrelation function decay.
In the case of a mixture of foam and NPs, the presence of two kinds of scatterers results in a twofold dynamics in the system. The solution of the problem has been demonstrated in detail elsewhere,9 and considers the impact of both dynamical processes in the decorrelation rate. The characteristic timescale of bubble reorganization is taken as τR = 1/Al*, where A is a constant depending on the bulk properties of the dispersion and l* is of the order of a few bubble diameters. For the Brownian motion of NPs in the continuous liquid phase the characteristic time is τB = 1/4k2D, where D is the Brownian diffusion coefficient of colloidal particles, k = 2πn/λ, with λ the wavelength inside the material and n the refractive index of the continuous media. The contribution of colloidal particles to the total scattering can be characterized by introducing the mean distance lC between two scattering events by NPs. The definition is related to the number of NPs in a sample, lC = 1/nCσC, where nC is the number of particles per unit volume and σC is the scattering cross-section of the particles. Considering the probabilities of indicated scattering events, it has been shown9 that the field autocorrelation function in the case of a parallel slab of thickness L is
(3) |
(4) |
To summarize briefly, the dynamical processes in the system cause temporal fluctuations of the signal – the intensity of scattered light. The rate of these fluctuations is characterized by the autocorrelation function. This function describes the extent of correlation (or simply similarity) between signals separated by a certain time segment. The decay of the function depends on the rate of the corresponding internal dynamics of the system, i.e. on the size and number of scatterers.
To prepare a sample with a specific volume fraction of NPs ΦC, defined as the volume of colloidal particles over the volume of the foam, we followed the procedure described in ref. 9. A controlled volume of NP suspension was added to a fixed mass of foam and the mixture was stirred for 1–2 minutes.
Afterwards the mixture was placed in the measuring cell, the glass walls of which are at an angle of about 2° (Fig. 1). This geometry makes it possible to carry out measurements for various thicknesses in the same foam at similar ages without preparing numerous samples for the same system. The length of the cell is 28 cm and the height of the walls is 10 cm. With this cell one can vary the sample thickness roughly between 5 and 10 mm while avoiding any leakage of scattered light and being sure that the measurements are made in the multiple scattering regime (L ≫ l*).
The scheme of the experimental setup is presented in Fig. 1. The measurements were carried out in the transmission configuration. Since the method is based on the measurement of the averaged intensity of the scattered light, the probing area was extended by expanding the laser beam (532 nm, 100 mV, Compass 315M from COHERENT) with a lens. The resulting diameter of the light spot was around 5 mm. Doing this one can be sure that even for aged foams, with bubbles of diameter up to about 200 μm,16 the signal is averaged over sufficient numbers of bubbles. The intensity of scattered light was measured with a single-mode optical fiber with collimation optics. The fiber output is sent to a photomultiplier, and the normalized intensity autocorrelation function gI(τ) is measured with a digital correlator (Flex03LQ-12 from http://Correlator.com) operating in multi-tau mode. The typical photon rate is 7–10 kHz and correlation functions are obtained after 1–5 min averages depending on signal statistics.
The lower plateau at long times corresponds to a complete loss of correlation. The decay rates τB and τR correspond to the two dynamical processes. At an early age of the system, τB and τR are quite close to each other and the decays are merged and indistinguishable. Upon foam aging, the bubble size increases, τR shifts towards longer times and is about 0.2 s eight hours after sample preparation. Then the two decays related to the different dynamical processes become well separated and an intermediate plateau appears. The same tendency is observed in Fig. 3(b) for the normalized gE(τ). In a log–lin scale: one can distinguish between the slopes of the “fast” decay due to the Brownian motion of particles and the “slow” decay related to bubble rearrangements, which is getting more pronounced.
In ref. 9 it is shown that the intercept of the exponential decay of gE(τ) (and consequently the height of the intermediate plateau of gI(τ) after the fast decay) depends directly on the volume fraction of colloidal particles. In the limit τ ≫ τB, g(B)E ≅ 0 and the electric field autocorrelation function (3) can be written as
(5) |
Fig. 4 Autocorrelation functions (a) gI(τ) and (b) gE(τ) (zoomed in the inset) for the mixture of Gillette foam with silica NP (ΦC = 1.3 × 10−4) at various cell thickness and at the same sample age. |
The dependencies of the correlation functions gE(τ) with the geometry of the cell L and with the properties of the material may be separated by rewriting (3) as
Fig. 5 Master curves for a mixture of Gillette foam with (a) silica NP (ΦC = 1.3 × 10−4) aged 24 hours and (b) latex beads (ΦC = 1.4 × 10−4) aged 27 hours. |
Since the duration of the experiment was a few minutes, we can be sure that the system didn't evolve significantly and that the scattering properties didn't change dramatically during the measurements. The cell we use allows us to vary the sample thickness in the range of 5–10 mm. Due to the cell filling method, air cavities within the sample are sometimes present. This makes it difficult to carry out measurements properly over the full range of cell thicknesses on the same sample. Nevertheless, in Fig. 5(a) one can see that the coincidence of the curves is good within experimental error. Furthermore the experiments were repeated using latex nanobeads as NPs (Fig. 5(b)). In this case the scaling of data is also observed.
Thus it is shown that the developed theory adequately describes dynamics in the system and can be used to monitor its evolution.
We assume that the values of the scattering cross-section σC of the spherical colloidal particles dispersed in the foam can be calculated with Mie theory.21 The refractive index of particles is 1.46 and 1.60 for silica and latex respectively, the refractive index of the medium surrounding the particles is 1.344. This assumption may be tested with a complementary experiment where we performed measurements on the same system (same cell thicknesses, same foam age, same NPs and their concentrations), with two different optical wavelengths. Data are plotted in Fig. 6. We see that the slow decays are the same, but that the values of g*E are different. The values of nC and L are independent of the wavelength, and we expect the values of l* to be very close at the two wavelengths. Indeed scattering in foam is mainly due to reflection and refraction which depends on liquid refractive indices, which do not appreciably depend on wavelength. This is supported by the fact that the slow time decays are the same for the two samples.
Fig. 6 Normalized autocorrelation functions for a Gillette foam with volume fraction of silica NP 1.47 × 10−3 measured at laser wavelengths 633 nm (red squares) and 532 nm (green squares). |
We then expect the ratio (lng*E)/σC to be independent of the optical wavelength. We found for our systems that the values of lng*E vary by a factor of two, but the values of (lng*E)/σC vary by less than 10% ((lng*E)/σC = 5.1 × 10−3 nm at 532 nm, (lng*E)/σC = 4.7 × 10−3 nm at 633 nm). This analysis supports the fact that Mie scattering may be used to estimate the scattering cross section of a NP in the liquid phase. At the laser wavelength of 532 nm, we calculated for silica NP σC = 112 nm2 and for latex nanobeads σC = 108 nm2.
The values of l* may be measured by numerous techniques. However as Gillette foam has been used in many DWS studies, we use here the empirical formulas for the values of l*2 = l0*2 + atw with tw the age of the foam, and l*0 and a two numerical constants.22 The added NP should not modify noticeably the value of l* of the material.
We can now calculate the volume fraction of NPs for our two samples. For the foam with an input concentration of silica NP ΦC = 1.3 × 10−4 the volume fraction became ΦC = 9.4 × 10−6 after 24 hours. In the case of the foam with latex NP (input concentration ΦC = 1.4 × 10−4) the volume fraction was ΦC = 4.8 × 10−5 at a foam age of 27 hours. Thus we observed a decrease of NP concentration within the foam during the system aging. Because the estimation of the amount of NPs in the system is made on the basis of the signal decorrelation due to their Brownian motion, the decrease of NP volume fraction means that some of them are immobilized and cannot move randomly anymore. One can also see that the silica NP concentration is decreased more strongly than that of latex NP. To understand the origin of this change we followed the kinetic behaviour of the system.
The first time decay is related to the foam dynamics and the second one to the NP dynamics. For the dilute NP systems the foam dynamics is slower (see Fig. 3) and then (1/τf) ≅ L2/4lCl*τB. Eventually one can obtain
(1 − gE)l* ≅ CnCτ |
Thus the movement of some NPs is blocked and they don't contribute in the fast decay of autocorrelation function.
Since titanium dioxide (TiO2) nanoparticles are used for photocatalytic application and as additives for the enhancement of the refractive index in health-care products (sunscreens, some medications, cosmetics), there is a challenge for the detection and characterization of TiO2 nanoparticles in turbid media.23 Therefore we will make the calculations by taking NPs of TiO2 with the diameter of 25 nm as an example.
For simplicity we discuss the experiment in a transmission cell, but the results may be easily extended to other scattering geometries. The basic idea of this method is that the dynamics due to NP is fast compared to the emulsion or foam intrinsic dynamics. The two associated time scales are τR(l*)2/L2 for the matrix dynamics, and τBlCl*/L2 for the Brownian dynamics. So we expect the two dynamics to be well separated if (τB/τR) ≤ 0.1(l*/lC). It should be noted that this relation is independent of the cell thickness L and more generally of the scattering geometry. For a typical numerical application, we consider the case of NP of 25 nm diameter dispersed in water. The Brownian time is then τB ≅ 0.1 ms. If the scattering matrix has an intrinsic reorganization time τR ≅ 1 s and a transport mean free path of l* ≅ 0.1 mm, we should have lC ≤ 0.1 m for the two dynamics to be well separated. The proper ratio of τB/τR can be attained by slowing down the dynamics of the matrix, which in the case of a foam simply means to wait until a certain system age.
The value of lC is related to the concentration of free particles and to their scattering cross section. It may be calculated using Mie scattering theory.21,24 We plotted in Fig. 8 the values of lC for spherical NP of radius r for a solid concentration of ΦC = 10−3. For TiO2 NP of refractive index 2.67 and of diameter 25 nm, we find that lC ≅ 15 mm. Since lC ∝ Φ−1C, the condition lC ≤ 0.1 m becomes ΦC ≥ 1.5 × 10−4. This concentration condition ensures that the two timescales associated with the two dynamics will be well separated.
Fig. 8 Values of the scattering lengths as a function of the NP radius. Solutions are at a concentration ΦC = 10−3. Solid curves are results of Mie scattering for NP: blue – TiO2, refractive index n = 2.67; red – fused quartz, refractive index n = 1.46. Optical wavelength is 532 nm, and the surrounding medium is water. Corresponding colored dotted lines are Rayleigh scattering limits, and black dotted line is the limit for large spheres (efficiency factors Q = 2).24 |
Coming back to the geometry of the experiment, the cell thickness must be chosen such that the correlation decay due to NP is measurable, but does not mask the decay due to the internal dynamics. The upper limit of NP volume fraction is determined only by the ratio between the sample thickness L and characteristic scattering length scales l* and lC and not by the two timescales. Thus if we consider that we may measure accurately 0.1 < g*E < 0.99, we obtain in a transmission geometry a thickness range of 0.07l*lC < L < 1.1l*lC (i.e. 0.2 mm < L < 3.5 mm for the example of TiO2 NP). The time scale of the decay due to Brownian motion of NP is the same than in the DLS technique and is in the range of usual digital correlators.
We have discussed the possibilities of using the method to study dopant particles, however it could have applications where the particles followed are taking an active part in the dynamics and ageing of the turbid media. Such an example would be particle stabilized foams or emulsions, where the concentrations of adsorbed and free particles could be followed, similarly to that done in ref. 10.
This journal is © The Royal Society of Chemistry 2016 |