Download Brownian Motion Analysis: Understanding Particle Motion and Light Scattering and more Study notes Physics in PDF only on Docsity!
Physics
A STUDY OF DIFFUSION BY LIGHT SCATTERING
I. INTRODUCTION .....................................................................................................................................................
II. THEORETICAL CONSIDERATIONS ................................................................................................................
III. MEASUREMENTS................................................................................................................................................
IV. DATA ANALYSIS................................................................................................................................................
I. Introduction
In 1827 an English botanist, Robert Brown, noticed that pollen grains in water perform peculiarly erratic movements. Some simple experiments convinced him that the phenomenon, now called Brownian motion, was physical, not biological, and could be observed for any small particle freely suspended in a fluid (gas or liquid). With the development of kinetic theory in the late nineteenth century it became evident that Brownian motion was due to the continual buffeting of the suspended particle by the atoms of the surrounding fluid. The first reasonably satisfactory theory was advanced by A. Einstein in 1905 and subsequently refined by many authors. It would obviously be very tedious and difficult to follow the paths of a large number of individual particles to study Brownian motion. Instead, one can use the time-varying intensity of light scattered from the suspension to gain some information about the average motions of the particles. The procedure is to create a mathematical model of the particle motion, calculate the light scattering that would result, and compare the calculations with experiment. The basic thermally-induced motions are now well enough understood that light scattering measurements are used to determine many properties of particles or large molecules in suspensions, including size and shape distributions. Interactions between particles will also modify the individual motions in characteristic ways, allowing one to determine some features of the inter-particle forces. It has even been possible to determine the number of living bacteria in a culture medium by detecting how many are swimming actively. The present experiment utilizes light scattering to examine the Brownian motion in a particularly simple system, a dilute suspension of identical polystyrene spheres. In the next section we present a simple model of our system and calculate the scattered intensity I ( t ). The intensity is expected to vary randomly with time, so we also develop a probabilistic description of such fluctuating signals. Succeeding sections describe the hardware and software needed for the experiment.
I ( t ) = E 0 '^
2 exp (^) [ i q! • (! r (^) j!! r (^) k )] j , k
" (4)
where the j,k sum runs over the N particles. Evidently I ( t ) will vary as the particles move, changing! r ( t ) randomly. In order to describe the motion, we define a probability density P (! r' , t ' ;! r, t ). This is the probability of finding a given
particle in a small volume at! r' at time t ', given that it was at! r at time t. Since we assume that the particles do not interact, P (! r' , t ' ;! r, t ) must depend only on the coordinates of the chosen
particle, and must be the same function for all the particles. Further, there is nothing about the system which picks out a special origin of coordinates or time, so the function must depend only on spatial and time differences
R =! r' !! r and! = t' - t. A non-rigorous derivation of P
( R ,^ !;0, 0) proceeds as follows. Consider a region of fluid in
which the number of particles per unit volume n ( ! r ) varies with! r , as suggested by Fig. 2. The motions of the particles are random, so the number of particles crossing any plane in the fluid is proportional to the number density at that point. In Fig. 2, then, the number crossing in the + x direction is proportional to n ( x ) while the number crossing in the - x direction is proportional to n ( x + dx ). Defining
J n as the net number of particles crossing the plane per unit area per unit time,
one finds that
! J n =! D
" n (5)
The parameter D which makes this an equality is called the diffusion coefficient. It must also be true that the number of particles in a fixed volume changes only because particles enter or leave the volume at the boundaries (conservation of particles) so
n(x)
n(x + dx)
x^ Fig. 2 Net particle flux across a plane due to a concentration gradient.
! n ! t
J n = 0 (6)
Putting these together leads to the basic equation for the variation in number density
! n ! t " D #^2 n = 0 (7)
Since the probability of finding a particle in a small volume is directly proportional to the number density, this leads to the equation for P
( R ,^ !;0, 0) :
! P
! t
" D #^2 P = 0 (8)
We want a solution to this equation with the initial condition that the particle is at! r = 0 at t = 0,
or P (0, 0;0, 0 ) !"(! r ). One can show by direct differentiation that
P
( R ,^ !;0, 0) =^ (^4 " Dt )#3/ 2^ exp^ [ #^ R^2 / 4 Dt ] (9)
is a solution of Eq. 8 and fits the required initial condition. The diffusion coefficient D can be related to other properties of the system by invoking the Einstein relation
D = k (^) BT μ (10)
The quantity μ is the mobility, the ratio of drift velocity to driving force when the particle is pushed through the fluid by an external force. (Note that this equation relates a thermal equilibrium quantity, D , to the non-equilibrium response to an external force. This sort of relationship is a fundamental property of many-particle systems.) The mobility can be calculated from classical hydrodynamics for spheres of radius a which are large compared to the surrounding molecules, with the result
vdrift = Fext / 6 !" a (11)
and therefore
S ( !) = I *
+"
(^ $) I (0) e
% i !$ (^) d $ (16)
This is the expression that must be averaged over particle positions. The formal average of S ( !) is written
S ( !) = P (
R ,^ ";0, 0)^ S (^!^ ) d^3
R (17)
where the integration is over the illuminated volume of the sample. It is advantageous to interchange the order of integration on d " and d^3
R to get
S ( !) = I *( ") I (0)
#$
+$ % e
i !" (^) d " (18)
Writing out the average I *^ (! ) I (0) gives
I *^ (! ) I (0) = E 0 '^
4 exp[" i q! • (! r l ( !) "! r m (! ))] exp[ i! q • (! r j (0) "! r k (0))] j , k
l , m
(19)
Using our assumption that the motions of different particles are independent allows us to factor most of the terms into a product of one or more single-particle averages exp(! i q! •! r ). Since
there is no preferred position, these single-particle averages are zero, eliminating most of the N^4 terms in Eq. 19. The only terms which survive are those containing two pairs of equal indices, which are:
l = m j = k N! N = N^2 terms l = k m = j m " l N! ( N -1) # N^2 terms (20) l = j m = k m " j N! ( N -1) # N^2 terms
The first set of terms is just <1! 1> = 1. The second group can be rewritten as
exp[! i q! •
R ( " )] exp[! 2 i q! •! r l (0)] 2 = 0 (21)
since
R (! ) is independent of! r. The third group is non-zero:
exp[! i! q • (! r l ( " )!! r l (0))] exp[ i! q • (! r m ( ")!! r m (0))] = exp[! i! q •
R ( " )]
2 (22)
The averaging integral can be done analytically, yielding
exp[! i! q •
R ( ")]
2 = exp(! 2 Dq^2 " ) (23)
Combining the terms we at last get
I *^ (! ) I (0) = N^2 E 0 '^
4
[^1 +^ exp("^2 Dq^2 !)] (24)
which can be substituted into Eq. 18 to give
S ( !) = N^2 E 0 '^
4 2 "#(! ) + 4 Dq
2 (2 Dq^2 ) 2 +! 2
This indicates that the power spectrum of I ( t ) consists of a delta function at! = 0, representing the average intensity, and a broad peak with a frequency width proportional to 2 Dq^2. There are several points of comparison between the measured and calculated power spectra. First, one can verify that < S ( !)> has the predicted form by fitting the expected function to the data and checking the "^2. The width, found in the fitting procedure, should vary as sin^2 ( #/2) as the scattering angle # is varied. Finally, the extracted value of D should be approximately equal to that estimated from Eq. 12. If one or more of these tests fail, it would tell us that the model must include other factors, such as particle interactions or multiple-scattering of the light. A numerical estimate of the expected width is useful at this point. The viscosity of water at 20oC is about 0.01 poise, which gives aD! 2.2 " 10 -13^ sec -1^. For laser light with $ = 632nm, an index of refraction n = 1.33 and # = 90o we find q! 1.8 " 105 cm-1^. A typical particle radius would be a = 0.1μm so 2 Dq^2 /2 %! 230 Hz (The 2 % is included because the apparatus is calibrated in Hz, not angular frequency.) Accordingly, we must plan to measure fluctuation rates up to several hundred hertz in order to adequately capture the input signal.
contain a small amount of detergent, 0.1 wt% sodium dodecyl sulfate, and fungicide, 0.1 wt% sodium azide) Each sample is marked with the nominal sphere diameter and type of solvent. You will need to know the actual viscosity of the solvent at the measuring temperature, so note the room temperature when you take the data. Some pertinent solvent properties are tabulated below. An external supply and standard divider chain are used for the photo multiplier tube. The output from the photo multiplier goes to a current to voltage converter which drives a DC meter to indicate the average anode current. Do not exceed the maximum ratings for the tube: 1250V or 50 μA current (approximately full scale on the meter), whichever occurs first. After the meter, the DC portion of the signal is blocked, while the AC portion goes through a lowpass filter and then to the A/D converter. The filter starts to cut off the signal at about 2. kHz, and passes less than 1% of the input above 5.0 kHz. The output of the filter should be connected to the CH1 input on the acquisition input panel.
B. Data acquisition The raw data at each angle consists of several series of samples of I ( t ) over a time range 0! T. Each series is Fourier transformed and the average of the transforms is stored in a file for later analysis. The sampling procedure is repeated as necessary for other angles. These data are acquired with a PC, using an A/D converter board and a LabVIEW virtual instrument (VI). The remainder of this section explains how to configure the software and acquire these data. Several steps are needed to set up the apparatus: a. Position the arm in the desired position, turn off the room lights and adjust the HV to obtain an average current below half-scale on the IDC meter. b. Start the VI. Double-click the LabVIEW icon and then open the file NoiseSpectAve.vi in My Documents\LabVIEW Data. This should give a LabVIEW front panel. Go to View > Tools Palette and select the 'pointing finger' icon to get the control cursor. c. Set configuration parameters to match the hardware connections. Use the pull-downs or click in a box with the 'hand' cursor, type a value, and click ENTER or press the Enter key (Return won't work.). Typical values (with explanations) are: Physical Channel: Dev1/ai1 (analog input channel 1 on device 1) Minimum Value -5 (-5V expected minimum input) Maximum Value +5 (+5V expected maximum input) d. Set desired FFT and averaging parameters. Samples per channel: 2048 (must be power of 2) Sample rate (Hz): 10000 (to get Nyquist frequency > 5 kHz) Number of sweeps: 100 (desired number for averaging) Window Hann (general-purpose window) File output: OFF (until you want to save data)
e. Set FFT display parameters Display unit Vrms^2 (power is voltage squared) Log/Linear Scale Linear (but the semilog plot shows more detail) f. With the VI set for a small number of sweeps, click on Run. Windowed voltage vs time data will be displayed in the upper window for each sweep. After the specified number of sweeps the averaged power spectrum will be displayed in the lower window. Once the apparatus is adjusted you can obtain data at the chosen scattering angles by turning on the file output and running the VI. When the averaging is complete you will be prompted to enter a file name to save the data. (Pick distinctive names, so other students do not erase your files. If you intend to analyze with MATLAB on this machine, it is convenient to use the folder My Documents\MATLAB.)
Viscosity (centipoise) as a function of Optical index of refraction as a function temperature for water and methyl alcohol. of temperature for water and methyl alcohol
T(K) methanol water 280 0.7201 1. 285 0.6669 12409 290 0.6189 1. 295 0.5755 0. 300 0.5360 0. 305 0.5000 0. 310 0.4670 0.
T(K) methanol water 293.15 1.3284 1. 298.15 1.3265 1.
